[Gdal-dev] Problem reading projection information consistently
Luke Roth
roth.luke at gmail.com
Fri Jun 22 10:48:56 EDT 2007
I'm having an issue reading the projection information from an Ikonos
image. I know it to be UTM zone 17, and when I open it in OpenEV
(from FWTools 1.3.1) I get the WKT
PROJCS["WGS 84 / UTM zone 17N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-81],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1]]
This is what I would expect, but when I use gdalinfo or the API
command GetProjectionRef() (using GDAL from the same version of
FWTools) I get:
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["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-81],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32617"]]
Why is there a difference? I'm trying to write a script to convert a
large number of images and shapefiles to a common projection, but this
is causing a problem. Any ideas? If I should be sending this to the
FWTools mailing list instead, let me know.
Thanks,
Luke Roth
More information about the Gdal-dev
mailing list