<div dir="ltr">I noticed that qgis-2.0 also has a problem with this file, but places it around 54 degrees west longitude (if using on-the-fly reprojection to WGS84). <div><br></div><div>The CRS is different from that identified by gdal (the lon_0, units and datum are different)<div>
<br></div><div>qgis: </div><div>+proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs</div><div><p style="margin:0px">
<br></p><p style="margin:0px">gdal:</p></div><div>PROJ.4 : '+proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=1968500 +y_0=0 +datum=NAD83 +units=us-ft +no_defs '<br>
</div><div><br></div><div>OGC WKT :</div><div>PROJCS["NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet",</div><div> GEOGCS["NAD83",</div><div> DATUM["North_American_Datum_1983",</div>
<div> SPHEROID["GRS 1980",6378137,298.2572221010002,</div><div> AUTHORITY["EPSG","7019"]],</div><div> AUTHORITY["EPSG","6269"]],</div>
<div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div><div> AUTHORITY["EPSG","4269"]],</div><div> PROJECTION["Lambert_Conformal_Conic_2SP"],</div>
<div> PARAMETER["standard_parallel_1",40.96666666666667],</div><div> PARAMETER["standard_parallel_2",39.93333333333333],</div><div> PARAMETER["latitude_of_origin",39.33333333333334],</div>
<div> PARAMETER["central_meridian",-77.75],</div><div> PARAMETER["false_easting",6458320.416666666],</div><div> PARAMETER["false_northing",0],</div><div> UNIT["us_survey_feet",0.3048006096012192],</div>
<div> AUTHORITY["EPSG","32129"]]</div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Jan 17, 2014 at 8:37 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org" target="_blank">even.rouault@mines-paris.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Le vendredi 17 janvier 2014 21:28:35, Reynolds, Scott a écrit :<br>
<div><div class="h5">> Hi,<br>
><br>
> I've downloaded several DEMs (e.g 27002570PAS_dem.zip) from<br>
> <a href="ftp://pamap.pasda.psu.edu/pamap_lidar/cycle1/DEM/South/2008/20000000/" target="_blank">ftp://pamap.pasda.psu.edu/pamap_lidar/cycle1/DEM/South/2008/20000000/</a>.<br>
> Gdalinfo is reporting the corner coordinates to be in the vicinity of 91<br>
> degrees 31 minutes west longitude which should be more like 75 degrees 36<br>
> minutes west longitude. Can someone explain what is happening here?<br>
<br>
</div></div>GeoTIFF support for non-meter linear units has been a major source of<br>
headaches in GDAL / libgeotiff. People will appreciate looking at<br>
frmts/gtiff/gt_wkt_srs.cpp<br>
<br>
I'm not sure if that GeoTIFF file is correct or not regarding the values of the<br>
false easting/northing with respect to the GeoTIFF spec (many in the wild are<br>
broken), but when looking at what happens, I noticed that a transformation<br>
from meters to feet happened twice. The attached patch fixes that (and the<br>
corner coordinates are around 75°), but I'm not sure if it is completely<br>
correct in all cases.<br>
<br>
You should likely create a GDAL Trac ticket with your report and the patch,<br>
but I'm not sure who will dare hurting his head against this wall...<br>
<br>
Even<br>
<br>
><br>
> Thanks,<br>
> Scott<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Geospatial professional services<br>
<a href="http://even.rouault.free.fr/services.html" target="_blank">http://even.rouault.free.fr/services.html</a><br>
</font></span><br>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br></div>