[gdal-dev] Re: gdalinfo: corner coords in other projection

[Hermann Peifer]

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)