<br><font size=2 face="sans-serif">What is the wkt you are using for zone
12 nad83?</font>
<p><font size=3>&nbsp; </font>
<p>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>&quot;Henneke, Amanda M&quot;
&lt;amanda.m.henneke@boeing.com&gt;@lists.osgeo.org</b> </font>
<br><font size=1 face="sans-serif">Envoyé par : gdal-dev-bounces@lists.osgeo.org</font>
<p><font size=1 face="sans-serif">23/09/2009 01:39 PM</font>
<td width=59%>
<table width=100%>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">A</font></div>
<td><font size=1 face="sans-serif">&quot;gdal-dev@lists.osgeo.org&quot;
&lt;gdal-dev@lists.osgeo.org&gt;</font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">Objet</font></div>
<td><font size=1 face="sans-serif">RE: [gdal-dev] Translating between UTM
and Geographic</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br>
<br>
<br><tt><font size=2>Here is some more information:<br>
<br>
So, here is the WKT I'm getting from my code for the &quot;problem&quot;
source file: <br>
 &nbsp;<br>
<br>
PROJCS[&quot;NAD_1927_UTM_Zone_11N&quot;, <br>
<br>
GEOGCS[&quot;GCS_North_American_1927&quot;, <br>
<br>
DATUM[&quot;North_American_Datum_1927&quot;, <br>
<br>
SPHEROID[&quot;Clarke_1866&quot;,6378206.4,294.9786982]], <br>
<br>
PRIMEM[&quot;Greenwich&quot;,0], <br>
<br>
UNIT[&quot;Degree&quot;,0.017453292519943295]], <br>
<br>
PROJECTION[&quot;Transverse_Mercator&quot;], <br>
<br>
PARAMETER[&quot;latitude_of_origin&quot;,0], <br>
<br>
PARAMETER[&quot;central_meridian&quot;,-117], <br>
<br>
PARAMETER[&quot;scale_factor&quot;,0.9996], <br>
<br>
PARAMETER[&quot;false_easting&quot;,500000], <br>
<br>
PARAMETER[&quot;false_northing&quot;,0], <br>
<br>
UNIT[&quot;Meter&quot;,1]] <br>
<br>
 &nbsp;<br>
<br>
Here's the WKT I'm creating for the destination file: <br>
<br>
GEOGCS[&quot;Geographic Coordinate System&quot;, <br>
<br>
DATUM[&quot;WGS_1984&quot;, <br>
<br>
SPHEROID[&quot;WGS_1984&quot;,6378137,298.257223563]], <br>
<br>
PRIMEM[&quot;Greenwich&quot;,0], <br>
<br>
UNIT[&quot;degree&quot;,0.0174532925199433]] <br>
<br>
 &nbsp;<br>
<br>
Do you see any issues with that? &nbsp;It looks like the source file is
picking up that it's NAD27... <br>
 <br>
<br>
<br>
Amanda M. Henneke<br>
IG/DB Tools - TSS Capability Center<br>
<br>
<br>
 <br>
When One tugs at a single thing in Nature, he finds it attached to the
rest of the world. -John Muir<br>
 &nbsp;Please consider the environment before printing this e-mail.<br>
<br>
-----Original Message-----<br>
From: gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev-bounces@lists.osgeo.org]
On Behalf Of Henneke, Amanda M<br>
Sent: Wednesday, September 23, 2009 11:58 AM<br>
To: Craig Miller<br>
Cc: gdal-dev@lists.osgeo.org<br>
Subject: RE: [gdal-dev] Translating between UTM and Geographic<br>
<br>
Craig-<br>
<br>
I'm trying to have it in Geograhphic/WGS84, at least that's what I'm specifying
in the code...I agree I think it is a datum error, but I don't know what
to do to make it translate the datum from NAD27 to WGS84 correctly. <br>
<br>
I don't know if I can send the raster, but I can send the two shapefiles...<br>
<br>
Amanda M. Henneke<br>
IG/DB Tools - TSS Capability Center<br>
<br>
<br>
 <br>
When One tugs at a single thing in Nature, he finds it attached to the
rest of the world. -John Muir<br>
 &nbsp;Please consider the environment before printing this e-mail.<br>
<br>
-----Original Message-----<br>
From: Craig Miller [mailto:craig.miller@spatialminds.com]<br>
Sent: Wednesday, September 23, 2009 11:51 AM<br>
To: Henneke, Amanda M<br>
Subject: RE: [gdal-dev] Translating between UTM and Geographic<br>
<br>
Geographic is a projection, not a datum. &nbsp;What Datum is the data in
the &quot;Geographic Projection&quot; in? &nbsp;E.g. you can have a map
in either<br>
Geographic/NAD27 or Geographic/NAD83 (or another datum for that matter..<br>
such as WGS84).<br>
<br>
What software are you using to display the data?<br>
<br>
If you want me to take a look, zip up the raster, and the good shapefile,
and the bad shapefile and send them over. &nbsp;I don't have a lot of time
as I'm getting ready to leave on a trip, but I'll try and squeeze it in
before I leave.<br>
<br>
Craig<br>
<br>
<br>
&gt; -----Original Message-----<br>
&gt; From: Henneke, Amanda M [mailto:amanda.m.henneke@boeing.com]<br>
&gt; Sent: Wednesday, September 23, 2009 9:28 AM<br>
&gt; To: Craig Miller<br>
&gt; Subject: RE: [gdal-dev] Translating between UTM and Geographic<br>
&gt; <br>
&gt; Hi Craig-<br>
&gt; <br>
&gt; I'm translating a shapefile (from UTM to Geographic) and then <br>
&gt; comparing it<br>
to<br>
&gt; <br>
&gt; #1 A shapefile that was translated from UTM to Geographic using the<br>
software,<br>
&gt; GlobalMapper.<br>
&gt; #2 A raster file that was originally projected in Geographic, so that
<br>
&gt; for example I can see the points don't match up with the features
they <br>
&gt; should represent on the image.<br>
&gt; <br>
&gt; The shapefile that GlobalMapper translates matches up perfectly with
<br>
&gt; my<br>
raster<br>
&gt; file.<br>
&gt; <br>
&gt; So my thought process is, that if Global Mapper has some way to <br>
&gt; convert<br>
the<br>
&gt; different projections without the shifting, how come GDAL can't do
that.<br>
I'm<br>
&gt; thinking it can, but maybe I need to give it some additional data
maybe?<br>
Or<br>
&gt; am I off?<br>
&gt; <br>
&gt; Thanks for your help!<br>
&gt; <br>
&gt; Amanda M. Henneke<br>
&gt; IG/DB Tools - TSS Capability Center<br>
&gt; <br>
&gt; <br>
&gt; When One tugs at a single thing in Nature, he finds it attached to
the<br>
rest of<br>
&gt; the world. -John Muir<br>
&gt; &nbsp; Please consider the environment before printing this e-mail.<br>
&gt; <br>
&gt; -----Original Message-----<br>
&gt; From: Craig Miller [mailto:craig.miller@spatialminds.com]<br>
&gt; Sent: Wednesday, September 23, 2009 11:21 AM<br>
&gt; To: Henneke, Amanda M<br>
&gt; Subject: RE: [gdal-dev] Translating between UTM and Geographic<br>
&gt; <br>
&gt; If you are transforming the datum from NAD27 to WGS84 then it is <br>
&gt; expected<br>
that<br>
&gt; your points would be shifted exactly as you describe. &nbsp;How are
you<br>
determining<br>
&gt; that they are shifted?<br>
&gt; <br>
&gt; a. &nbsp;Compared to a map? &nbsp;Is the map reprojected? &nbsp;Are
you comparing <br>
&gt; them to landform features or to the lat/lon grid on the map? &nbsp;For
<br>
&gt; example, if you<br>
have<br>
&gt; a USGS Quad map with lat/lon lines on it based on the NAD27 datum
then<br>
those<br>
&gt; lat/lon lines ONLY show the lat/lon using that datum. &nbsp;If you
plot <br>
&gt; points<br>
that<br>
&gt; have been shifted to the new datum, they will not agree with those<br>
&gt; NAD27<br>
lines<br>
&gt; since the new points are based on lat/lon NAD83.<br>
&gt; <br>
&gt; b. &nbsp;Comparing lat/lon directly? &nbsp;The points have been transformed.<br>
Lat/lon<br>
&gt; are relative to a datum and lan/lon NAD27 coords and lat/lon NAD83
<br>
&gt; coords<br>
will<br>
&gt; have different values.<br>
&gt; <br>
&gt; Craig<br>
&gt; <br>
&gt; <br>
&gt; &gt; -----Original Message-----<br>
&gt; &gt; From: gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev- <br>
&gt; &gt; bounces@lists.osgeo.org] On Behalf Of Henneke, Amanda M<br>
&gt; &gt; Sent: Wednesday, September 23, 2009 8:17 AM<br>
&gt; &gt; To: gdal-dev@lists.osgeo.org<br>
&gt; &gt; Subject: [gdal-dev] Translating between UTM and Geographic<br>
&gt; &gt;<br>
&gt; &gt; Hello all-<br>
&gt; &gt;<br>
&gt; &gt; I'm trying to figure out how to reproject some points from a
UTM<br>
&gt; projection to<br>
&gt; &gt; Geographic (WGS84) projection. &nbsp;This works fine when I'm
working in <br>
&gt; &gt; UTM<br>
&gt; Zone<br>
&gt; &gt; 12 with NAD83. &nbsp;When I try the same code for Zone 11 NAD27,
I get a <br>
&gt; &gt; slight shift in the points (to the east) and an extremely slight
<br>
&gt; &gt; shift up (to the north).<br>
&gt; &gt;<br>
&gt; &gt; I've included the following code as it stands now, I've tried
all <br>
&gt; &gt; sorts of things thinking that I need to somehow include more
<br>
&gt; &gt; information for the<br>
&gt; new<br>
&gt; &gt; spatial reference (the destination SR), and so far I'm coming
up<br>
empty...<br>
&gt; &gt;<br>
&gt; &gt; Anyone have any ideas what I'm missing?<br>
&gt; &gt;<br>
&gt; &gt; Where sSpRefParams is a string that contains the WKT projection<br>
&gt; information<br>
&gt; &gt; from a valid layer...<br>
&gt; &gt;<br>
&gt; &gt; srOld = new OSGeo.OSR.SpatialReference(sSpRefParams);<br>
&gt; &gt;<br>
&gt; &gt; srNew.SetGeogCS(&quot;Geographic Coordinate System&quot;,<br>
&gt; OSGeo.OSR.Osr.SRS_DN_WGS84,<br>
&gt; &gt; OSGeo.OSR.Osr.SRS_DN_WGS84, OSGeo.OSR.Osr.SRS_WGS84_SEMIMAJOR,<br>
&gt; &gt; OSGeo.OSR.Osr.SRS_WGS84_INVFLATTENING, &quot;Greenwich&quot;,
0.0, &quot;degree&quot;, <br>
&gt; &gt; 0.0174532925199433);<br>
&gt; &gt;<br>
&gt; &gt; ctTransform = new OSGeo.OSR.CoordinateTransformation(srOld, srNew);<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; Amanda M. Henneke<br>
&gt; &gt; IG/DB Tools - TSS Capability Center<br>
&gt; <br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; When One tugs at a single thing in Nature, he finds it attached
to <br>
&gt; &gt; the<br>
&gt; rest of<br>
&gt; &gt; the world. -John Muir<br>
&gt; &gt; &nbsp; Please consider the environment before printing this e-mail.<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; gdal-dev mailing list<br>
&gt; &gt; gdal-dev@lists.osgeo.org<br>
&gt; &gt; http://lists.osgeo.org/mailman/listinfo/gdal-dev<br>
<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
gdal-dev@lists.osgeo.org<br>
http://lists.osgeo.org/mailman/listinfo/gdal-dev<br>
_______________________________________________<br>
gdal-dev mailing list<br>
gdal-dev@lists.osgeo.org<br>
http://lists.osgeo.org/mailman/listinfo/gdal-dev<br>
</font></tt>
<br>
<br>
<br>