<div dir="ltr">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.<br>
<br>Thanks,<br>
<br>
Kerry Halligan<br><br><br>The gdalinfo result for this file includes:<br><br>Coordinate System is:<br>PROJCS["NAD83 / New York Long Island",<br> GEOGCS["NAD83",<br> DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.257222101,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6269"]],<br> PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br> UNIT["degree",0.01745329251994328,<br> AUTHORITY["EPSG","9122"]],<br> AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Lambert_Conformal_Conic_2SP"],<br> PARAMETER["standard_parallel_1",41.03333333333333],<br> PARAMETER["standard_parallel_2",40.66666666666666],<br> PARAMETER["latitude_of_origin",40.16666666666666],<br>
PARAMETER["central_meridian",-74],<br> PARAMETER["false_easting",984250],<br> PARAMETER["false_northing",0],<br> UNIT["us_survey_feet",0.3048006096012192]]<br>Origin = (983922.937212173830000,197166.639292352130000)<br>
Pixel Size = (10.000000000000000,-10.000000000000000)<br>Corner Coordinates:<br>Upper Left ( 983922.937, 197166.639) ( 74d 0'4.25"W, 40d42'28.27"N)<br>Lower Left ( 983922.937, 185976.639) ( 74d 0'4.24"W, 40d40'37.70"N)<br>
Upper Right ( 992382.937, 197166.639) ( 73d58'14.40"W, 40d42'28.26"N)<br>Lower Right ( 992382.937, 185976.639) ( 73d58'14.44"W, 40d40'37.69"N)<br>Center ( 988152.937, 191571.639) ( 73d59'9.33"W, 40d41'32.98"N)<br>
<br>If I use dataset.GetProjection() on this file, it returns the projection string as:<br><br>PROJCS["NAD83 / New York Long Island",<br> GEOGCS["NAD83",<br> DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 80",6378137,298.257222101,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6269"]],<br> PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br> UNIT["degree",0.01745329251994328,<br> AUTHORITY["EPSG","9122"]],<br> AUTHORITY["EPSG","4269"]]<br>
PROJECTION["Lambert_Conformal_Conic_2SP"]<br> PARAMETER["standard_parallel_1",41.03333333333333]<br> PARAMETER["standard_parallel_2",40.66666666666666]<br> PARAMETER["latitude_of_origin",40.16666666666666]<br>
PARAMETER["central_meridian",-74],<br> PARAMETER["false_easting",984250]<br> PARAMETER["false_northing",0],<br> UNIT["us_survey_feet",0.3048006096012192]]<br><br>What I expected to see for this dataset is:<br>
<pre>PROJCS["NAD83 / New York Long Island (ftUS)",<br> GEOGCS["NAD83",<br> DATUM["North_American_Datum_1983",<br> SPHEROID["GRS 1980",6378137,298.257222101,<br>
AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6269"]],<br> PRIMEM["Greenwich",0,<br> AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.01745329251994328,<br> AUTHORITY["EPSG","9122"]],<br> AUTHORITY["EPSG","4269"]],<br> UNIT["US survey foot",0.3048006096012192,<br>
AUTHORITY["EPSG","9003"]],<br> PROJECTION["Lambert_Conformal_Conic_2SP"],<br> PARAMETER["standard_parallel_1",41.03333333333333],<br> PARAMETER["standard_parallel_2",40.66666666666666],<br>
PARAMETER["latitude_of_origin",40.16666666666666],<br> PARAMETER["central_meridian",-74],<br> PARAMETER["false_easting",984250.0000000002],<br> PARAMETER["false_northing",0],<br>
AUTHORITY["EPSG","2263"],<br> AXIS["X",EAST],<br> AXIS["Y",NORTH]]</pre><br>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.<br>
<br>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. <br>
<br><br><br></div>