[gdal-dev] Correct gdalwarp usage from UTM to Equidistant Cylindrical Projection

[Christopher Hunt]

On 16/07/2008, at 12:27 PM, Frank Warmerdam wrote:

> Everything looks fine, but I had actually been hoping to see the
> gdalinfo report for the *output* file, not the input.  The output
> should have a coordinate system attached, so we should actually see
> that lat/long corner coordinates in addition to the projected  
> coordinates.
Ah - sorry - here you are:

Driver: GTiff/GeoTIFF
Files: out.tif
Size is 3780, 2360
Coordinate System is:
PROJCS["unnamed",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.2572235630016,
                 AUTHORITY["EPSG","7030"]],
             AUTHORITY["EPSG","6326"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4326"]],
     PROJECTION["Equirectangular"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",0],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (-160532.836602925264742,5798698.057208529673517)
Pixel Size = (59.180187349429140,-59.180187349429140)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( -160532.837, 5798698.057) (  1d26'31.53"W, 52d 5'26.13"N)
Lower Left  ( -160532.837, 5659032.815) (  1d26'31.53"W, 50d50'9.44"N)
Upper Right (   63168.272, 5798698.057) (  0d34'2.82"E, 52d 5'26.13"N)
Lower Right (   63168.272, 5659032.815) (  0d34'2.82"E, 50d50'9.44"N)
Center      (  -48682.283, 5728865.436) (  0d26'14.35"W, 51d27'47.79"N)
Band 1 Block=3780x1 Type=Byte, ColorInterp=Red
   Mask Flags: PER_DATASET ALPHA
Band 2 Block=3780x1 Type=Byte, ColorInterp=Green
   Mask Flags: PER_DATASET ALPHA
Band 3 Block=3780x1 Type=Byte, ColorInterp=Blue
   Mask Flags: PER_DATASET ALPHA
Band 4 Block=3780x1 Type=Byte, ColorInterp=Alpha

Thanks again for the verification.

Cheers,
-C