[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