Strange looking gdalinfo output

Stephen Davies sdc at SDC.COM.AU
Fri Dec 16 21:26:17 EST 2005


I recently purchased a collection of images from the local State 
Government. The images were provided in both ECW and TIF formats.

If I do a gdalinfo on the TIFs, I get reasonable looking output:

Driver: GTiff/GeoTIFF
Size is 4270, 4079
Coordinate System is:
PROJCS["GDA94 / MGA zone 53",
    GEOGCS["GDA94",
        DATUM["Geocentric_Datum_of_Australia_1994",
            SPHEROID["GRS 1980",6378137,298.2572221010042,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6283"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4283"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",135],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","28353"]]
Origin = (701028.967763,6051913.962025)
Pixel Size = (2.00000000,-2.00000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=ER Mapper 7.0
  TIFFTAG_IMAGEDESCRIPTION=ER Mapper GeoTiff raster translator V1.0: 
Band 1 = Red, Band 2 = Green, Band 3 = Blue
  TIFFTAG_SOFTWARE=ER Mapper 7.0
  TIFFTAG_XRESOLUTION=0
  TIFFTAG_YRESOLUTION=0
Corner Coordinates:
Upper Left  (  701028.968, 6051913.962) (137d13'14.56"E, 35d39'22.23"S)
Lower Left  (  701028.968, 6043755.962) (137d13'21.91"E, 35d43'46.83"S)
Upper Right (  709568.968, 6051913.962) (137d18'53.94"E, 35d39'15.84"S)
Lower Right (  709568.968, 6043755.962) (137d19'1.60"E, 35d43'40.42"S)
Center      (  705298.968, 6047834.962) (137d16'8.00"E, 35d41'31.37"S)
Band 1 Block=4270x20 Type=Byte, ColorInterp=Red
Band 2 Block=4270x20 Type=Byte, ColorInterp=Green
Band 3 Block=4270x20 Type=Byte, ColorInterp=Blue


If I do it on ECW version of the same image, I get decidedly strange 
results:

Driver: ECW/ERMapper Compressed Wavelets
Size is 4270, 4079
Coordinate System is:
PROJCS["MGA53",
    GEOGCS["GEOCENTRIC DATUM of AUSTRALIA",
        DATUM["GDA94",
            SPHEROID["GRS80",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",135],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",10000000],
    PARAMETER["false_northing",500000],
    UNIT["Meter",1]]
Origin = (701028.967763,6051913.962025)
Pixel Size = (2.00000000,-2.00000000)
Corner Coordinates:
Upper Left  (  701028.968, 6051913.962) ( 75d26'18.41"E,515d 9'42.93"N)
Lower Left  (  701028.968, 6043755.962) ( 60d36'10.10"E,508d43'14.95"N)
Upper Right (  709568.968, 6051913.962) ( 67d59'32.56"E,511d19'30.87"N)
Lower Right (  709568.968, 6043755.962) ( 53d15'25.56"E,504d56'0.67"N)
Center      (  705298.968, 6047834.962) ( 64d17'19.78"E,510d 1'12.50"N)
Band 1 Block=4270x1 Type=Byte, ColorInterp=Red
  Overviews: arbitrary
Band 2 Block=4270x1 Type=Byte, ColorInterp=Green
  Overviews: arbitrary
Band 3 Block=4270x1 Type=Byte, ColorInterp=Blue
  Overviews: arbitrary


In particular, the false easting/northing values are swapped and the 
extent lat/long values are quite different.

Could this be a problem with gdalinfo or is the data wrong?

TIA,
Stephen Davies
-- 
 =======================================================================
This email is for the person(s) identified above, and is confidential to
the sender and the person(s).  No one else is authorised to use or
disseminate this email or its contents.

Stephen Davies Consulting                            Voice: 08-8177 1595
Adelaide, South Australia.                             Fax: 08-8177 0133
Computing & Network solutions.                       Mobile:0403 0405 83



More information about the mapserver-users mailing list