Hello List.
At first i wanna thanks for the great work you are doing with the gdal
library. It is a source of great tools and opportunities for spatial issues
which whom i involed in.
Among others i am using the commandline tool gdal_translate for clipping
raster files via the following example command:
gdal_translate -projwin 13.6666664938 53.2000002428
13.8333334784 53.1000017686 source.tif target.tif
The source tif is already georeferenced (EPSG:4326). I would now assume that
the boundingbox of the following GTiff is somethink like this:
Lower Left: 13.6666664938 53.1000017686
Upper Right: 13.8333334784 53.2000002428
Instead it is something like this (the full gdalinfo dialog you can see in
the attach to this mail):
Lower Left: 13.6666610463 53.1000244672
Upper Right: 13.8333191458 53.2000151873
So there is a small imprecision in the clipping result. My question is now
if this imprecision is normal (maybe a result of rounding operations within
the clipping process) or if i do something wrong (do i miss something in my
command).
I am using libgdal1 version 1.9.0-3.1ubunut1. Looking forward to hearing
from you.
Kind regards,
JM
Driver: GTiff/GeoTIFF
Files: clipped.tif
Size is 8052, 4831
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (13.666661046340586,53.200015187338956)
Pixel Size = (0.000020697727207,-0.000020697727207)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 13.6666610, 53.2000152) ( 13d39'59.98"E, 53d12' 0.05"N)
Lower Left ( 13.6666610, 53.1000245) ( 13d39'59.98"E, 53d 6' 0.09"N)
Upper Right ( 13.8333191, 53.2000152) ( 13d49'59.95"E, 53d12' 0.05"N)
Lower Right ( 13.8333191, 53.1000245) ( 13d49'59.95"E, 53d 6' 0.09"N)
Center ( 13.7499901, 53.1500198) ( 13d44'59.96"E, 53d 9' 0.07"N)
Band 1 Block=8052x1 Type=Byte, ColorInterp=Red
Band 2 Block=8052x1 Type=Byte, ColorInterp=Green
Band 3 Block=8052x1 Type=Byte, ColorInterp=Blue
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
