On 24/05/2011 18:15, Lefman, Jonathan ERDC-TEC-VA wrote:
Hi all,
Is there an executable utility that takes the corner coordinates from a gtiff
and translates them into another projection system? For example, I have a
gtiff in UTM and I want to translate the coordinates of the bounding box (or
corners) into WGS84. Or is there a way to do this non-programmatically using
a combination of executable utilities?
Couldn't you just simply use gdalinfo's output, see below?
Hermann
$ gdalinfo highway_number_3035.tif
Driver: GTiff/GeoTIFF
Files: highway_number_3035.tif
Size is 10384, 7031
Coordinate System is:
PROJCS["ETRS89 / ETRS-LAEA",
GEOGCS["ETRS89",
DATUM["European_Terrestrial_Reference_System_1989",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4258"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
PARAMETER["latitude_of_center",52],
PARAMETER["longitude_of_center",10],
PARAMETER["false_easting",4321000],
PARAMETER["false_northing",3210000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3035"]]
Origin = (270937.427360458299518,5929526.137714201584458)
Pixel Size = (701.323508343892058,-701.323508343892058)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 270937.427, 5929526.138) ( 67d52'39.91"W, 53d13'23.14"N)
Lower Left ( 270937.427, 998520.551) ( 29d49'49.65"W, 22d27'34.90"N)
Upper Right ( 7553480.738, 5929526.138) ( 81d 0' 6.31"E, 59d42'39.88"N)
Lower Right ( 7553480.738, 998520.551) ( 42d23'31.15"E, 25d49'31.87"N)
Center ( 3912209.083, 3464023.344) ( 3d44'23.41"E, 54d 7'25.38"N)
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
