<br><font size=2 face="sans-serif">As one more data point, I used some of our internal software to convert UTM, zone 16, NAD27 to Geographic, WGS84 and received the same results as OGR. &nbsp;The software is based on GCTP and our own datum layer.</font>
<br><font size=2 face="sans-serif"><br>
--------------------------------------------------------------------------<br>
Tim Beckmann &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; tbeckman@usgs.gov<br>
Software Project Lead<br>
SAIC<br>
EROS Data Center, Sioux Falls, SD 57198<br>
605-594-2521 &nbsp; &nbsp;Phone<br>
605-594-6940 &nbsp; &nbsp;Fax<br>
</font>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td>
<td><font size=1 face="sans-serif"><b>Frank Warmerdam &lt;warmerdam@pobox.com&gt;</b></font>
<br><font size=1 face="sans-serif">Sent by: gdal-dev-admin@remotesensing.org</font>
<p><font size=1 face="sans-serif">12/12/2003 08:18 AM</font>
<br><font size=1 face="sans-serif">Please respond to gdal-dev</font>
<br>
<td><font size=1 face="Arial">&nbsp; &nbsp; &nbsp; &nbsp; </font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; To: &nbsp; &nbsp; &nbsp; &nbsp;gdal-dev@remotesensing.org</font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; cc: &nbsp; &nbsp; &nbsp; &nbsp;</font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; Subject: &nbsp; &nbsp; &nbsp; &nbsp;Re: [Gdal-dev] USGS vs. ESRI vs. OGR</font></table>
<br>
<br>
<br><font size=2 face="Courier New">Grabowski, Hank wrote:<br>
&gt; I followed your advise to compare to the USGS tool. &nbsp;The conversion<br>
&gt; results are listed below. &nbsp;As we can see, the ESRI and USGS tools match<br>
&gt; identically (to the precision of ESRI which was lower than the USGS<br>
&gt; solution). &nbsp;There is a convergence parameter from the USGS calculation<br>
&gt; which is significantly larger than the error between OGR and ESRI/USGS.<br>
&gt; Being a novice at this reprojection business, I don't know if that is<br>
&gt; the bounds of the accuracy of the solution. &nbsp;I would tend to believe not<br>
&gt; however, since a 40 minute minimum accuracy would be a multi-kilometer<br>
&gt; error in position. &nbsp;On a hunch I looked at the two definitions used for<br>
&gt; WGS84 by both OGR and ESRI, and their coefficients match to reported<br>
&gt; accuracy. &nbsp;I'll continue investigating this issue if no one has any<br>
&gt; ideas what could be causing it.<br>
&gt; <br>
&gt; UTM X &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 597196.62 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Meters &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>
&gt; UTM Y &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 4116325.5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Meters &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>
&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <br>
&gt; Product &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Lat &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Lon &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Convergence<br>
&gt; Scale<br>
&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; DD MM SS.sssss &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;DDD MM SS.sssss &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;DD MM SS.ss<br>
&gt; Unitless<br>
&gt; USGS &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;37 11 24.57282 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;085 54 17.71610 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0 39 43.15<br>
&gt; 9.99716E-001<br>
&gt; ESRI &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;37 11 24.57 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;085 54 17.72 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>
&gt; OGR &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;37 11 24.76 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 085 54 17.62 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>
<br>
Hank,<br>
<br>
I am a bit confused. &nbsp;When I do UTM16 NAD27 to geographic NAD27 I get:<br>
<br>
cs2cs +proj=utm +zone=16 +datum=NAD27 +to +proj=latlong +datum=NAD27<br>
597196.62 4116325.5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (input)<br>
85d54'17.716&quot;W &nbsp;37d11'24.573&quot;N 0.000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(output)<br>
<br>
It would *appear* that your locations for the USGS and ESRI are the NAD27<br>
lat/long location, not the NAD83 location unless the datum shift at the<br>
location is less than 1/100 of a second.<br>
<br>
I broke down the conversion into steps from UTM 16 NAD27 to Latlong NAD27,<br>
and then used &quot;nad2nad&quot; explicitly with the conus file just to be sure<br>
that my code for cs2cs wasn't mixing things up, and I get the same final<br>
OGR output coordinates you did. &nbsp;While the ESRI software uses a somewhat<br>
different grid shift file than ESRI, I believe that the PROJ/OGR software<br>
is doing the right thing given the values in the conus file.<br>
<br>
In summary, I question whether your USGS/ESRI numbers really represent<br>
applying a datum shift at all.<br>
<br>
Best regards,<br>
-- <br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up &nbsp; | Frank Warmerdam, warmerdam@pobox.com<br>
light and sound - activate the windows | http://pobox.com/~warmerdam<br>
and watch the world go round - Rush &nbsp; &nbsp;| Geospatial Programmer for Rent<br>
<br>
<br>
_______________________________________________<br>
Gdal-dev mailing list<br>
Gdal-dev@remotesensing.org<br>
http://remotesensing.org/mailman/listinfo/gdal-dev<br>
</font>
<br>
<br>