Strange looking gdalinfo output
Stephen Davies
sdc at SDC.COM.AU
Fri Dec 16 18:26:17 PST 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