[gdal-dev] problem with ERDAS Imagine files in StatePlane projection

[Kerry Halligan]

I am using GDAL with C# bindings from gdal_csharp.dll and reading the
projection information from an ERDAS Imagine file in Stateplane Projection
(EPSG 2263).  GDAL doesn't seem to be able to parse the stateplane
projection information in the Imagine file correctly when writing out as a
GeoTIFF.  My assumption is that ERDAS is using a non standard syntax and
this is preventing a match of the PROJCS tag?  Has anyone encountered this
before?  Any clues or ideas at all would be much appreciated.

Thanks,

Kerry Halligan


The gdalinfo result for this file includes:

Coordinate System is:
PROJCS["NAD83 / New York Long Island",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",41.03333333333333],
    PARAMETER["standard_parallel_2",40.66666666666666],
    PARAMETER["latitude_of_origin",40.16666666666666],
    PARAMETER["central_meridian",-74],
    PARAMETER["false_easting",984250],
    PARAMETER["false_northing",0],
    UNIT["us_survey_feet",0.3048006096012192]]
Origin = (983922.937212173830000,197166.639292352130000)
Pixel Size = (10.000000000000000,-10.000000000000000)
Corner Coordinates:
Upper Left  (  983922.937,  197166.639) ( 74d 0'4.25"W, 40d42'28.27"N)
Lower Left  (  983922.937,  185976.639) ( 74d 0'4.24"W, 40d40'37.70"N)
Upper Right (  992382.937,  197166.639) ( 73d58'14.40"W, 40d42'28.26"N)
Lower Right (  992382.937,  185976.639) ( 73d58'14.44"W, 40d40'37.69"N)
Center      (  988152.937,  191571.639) ( 73d59'9.33"W, 40d41'32.98"N)

If I use dataset.GetProjection() on this file, it returns the projection
string as:

PROJCS["NAD83 / New York Long Island",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 80",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]]
    PROJECTION["Lambert_Conformal_Conic_2SP"]
    PARAMETER["standard_parallel_1",41.03333333333333]
    PARAMETER["standard_parallel_2",40.66666666666666]
    PARAMETER["latitude_of_origin",40.16666666666666]
    PARAMETER["central_meridian",-74],
    PARAMETER["false_easting",984250]
    PARAMETER["false_northing",0],
    UNIT["us_survey_feet",0.3048006096012192]]

What I expected to see for this dataset is:

PROJCS["NAD83 / New York Long Island (ftUS)",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]],
    UNIT["US survey foot",0.3048006096012192,
        AUTHORITY["EPSG","9003"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",41.03333333333333],
    PARAMETER["standard_parallel_2",40.66666666666666],
    PARAMETER["latitude_of_origin",40.16666666666666],
    PARAMETER["central_meridian",-74],
    PARAMETER["false_easting",984250.0000000002],
    PARAMETER["false_northing",0],
    AUTHORITY["EPSG","2263"],
    AXIS["X",EAST],
    AXIS["Y",NORTH]]


If I use the string returned by GetProjection() as the projection
information for an output GeoTIFF from GDAL, the image displays with a large
horizontal offset despite having the same origin coordinates.  It appears
that the PROJCS information is not being handled correctly, so it is just
coming up as Lambert Conformal Conic, but not stateplane New York Long
Island.

If I export the data as a GeoTIFF (not as an Imagine file) from ERDAS, then
the projection information is read in / written out just fine by GDAL.
Based on this, I am assuming that GDAL isn't parsing out the PROJCS portion
of the projection string correctly, likely because ERDAS isn't writing it
correctly.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20080909/379aaa79/attachment.html