<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.3603" name=GENERATOR></HEAD>
<BODY>
<DIV dir=ltr align=left><FONT size=2>
<P>PROJCS["NAD_1983_UTM_Zone_12N",</P>
<P>GEOGCS["GCS_North_American_1983",</P>
<P>DATUM["North_American_Datum_1983",</P>
<P>SPHEROID["GRS_1980",6378137,298.257222101]],</P>
<P>PRIMEM["Greenwich",0],</P>
<P>UNIT["Degree",0.017453292519943295]],</P>
<P>PROJECTION["Transverse_Mercator"],</P>
<P>PARAMETER["latitude_of_origin",0],</P>
<P>PARAMETER["central_meridian",-111],</P>
<P>PARAMETER["scale_factor",0.9996],</P>
<P>PARAMETER["false_easting",500000],</P>
<P>PARAMETER["false_northing",0],</P>
<P>UNIT["Meter",1]]</P></FONT></DIV>
<DIV>&nbsp;</DIV><!-- Converted from text/rtf format -->
<P><B><SPAN lang=en-us><FONT face="Kristen ITC" color=#006291 size=2>Amanda M. 
Henneke</FONT></SPAN></B><SPAN lang=en-us></SPAN> <BR><SPAN lang=en-us><FONT 
face="Comic Sans MS" color=#808080 size=1>IG/DB Tools - TSS&nbsp;Capability 
Center</FONT></SPAN> <BR><SPAN lang=en-us><FONT face="Comic Sans MS" 
color=#808080 size=1>Boeing - St. Louis</FONT></SPAN> <BR><SPAN lang=en-us><FONT 
face="Comic Sans MS" color=#808080 size=1>314-777-6856</FONT></SPAN> </P>
<P><SPAN lang=en-us><FONT face=Arial></FONT></SPAN> 
<P align=center><SPAN lang=en-us><FONT face="Comic Sans MS" color=#008000 
size=1>When One tugs at a single thing in Nature, he finds it attached to the 
rest of the world. -John Muir</FONT></SPAN></P>
<P align=center><SPAN lang=en-us><FONT face="Comic Sans MS" color=#008000 
size=1>&nbsp; Please consider the environment before printing this 
e-mail.</FONT></SPAN></P>
<DIV>&nbsp;</DIV><BR>
<DIV class=OutlookMessageHeader lang=en-us dir=ltr align=left>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>From:</B> Steve.Toutant@inspq.qc.ca 
[mailto:Steve.Toutant@inspq.qc.ca] <BR><B>Sent:</B> Wednesday, September 23, 
2009 12:41 PM<BR><B>To:</B> Henneke, Amanda M<BR><B>Cc:</B> 
gdal-dev@lists.osgeo.org; gdal-dev-bounces@lists.osgeo.org<BR><B>Subject:</B> 
RE: [gdal-dev] Translating between UTM and Geographic<BR></FONT><BR></DIV>
<DIV></DIV><BR><FONT face=sans-serif size=2>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%">
  <TBODY>
  <TR vAlign=top>
    <TD width="40%"><FONT face=sans-serif size=1><B>"Henneke, Amanda M" 
      &lt;amanda.m.henneke@boeing.com&gt;@lists.osgeo.org</B> </FONT><BR><FONT 
      face=sans-serif size=1>Envoyé par : 
      gdal-dev-bounces@lists.osgeo.org</FONT> 
      <P><FONT face=sans-serif size=1>23/09/2009 01:39 PM</FONT> </P>
    <TD width="59%">
      <TABLE width="100%">
        <TBODY>
        <TR vAlign=top>
          <TD>
            <DIV align=right><FONT face=sans-serif size=1>A</FONT></DIV>
          <TD><FONT face=sans-serif size=1>"gdal-dev@lists.osgeo.org" 
            &lt;gdal-dev@lists.osgeo.org&gt;</FONT> 
        <TR vAlign=top>
          <TD>
            <DIV align=right><FONT face=sans-serif size=1>cc</FONT></DIV>
          <TD>
        <TR vAlign=top>
          <TD>
            <DIV align=right><FONT face=sans-serif size=1>Objet</FONT></DIV>
          <TD><FONT face=sans-serif size=1>RE: [gdal-dev] Translating between 
            UTM and Geographic</FONT></TR></TBODY></TABLE><BR>
      <TABLE>
        <TBODY>
        <TR vAlign=top>
          <TD>
          <TD></TR></TBODY></TABLE><BR></TR></TBODY></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 "problem" source file: 
<BR>&nbsp;<BR><BR>PROJCS["NAD_1927_UTM_Zone_11N", 
<BR><BR>GEOGCS["GCS_North_American_1927", 
<BR><BR>DATUM["North_American_Datum_1927", 
<BR><BR>SPHEROID["Clarke_1866",6378206.4,294.9786982]], 
<BR><BR>PRIMEM["Greenwich",0], <BR><BR>UNIT["Degree",0.017453292519943295]], 
<BR><BR>PROJECTION["Transverse_Mercator"], 
<BR><BR>PARAMETER["latitude_of_origin",0], 
<BR><BR>PARAMETER["central_meridian",-117], 
<BR><BR>PARAMETER["scale_factor",0.9996], 
<BR><BR>PARAMETER["false_easting",500000], 
<BR><BR>PARAMETER["false_northing",0], <BR><BR>UNIT["Meter",1]] 
<BR><BR>&nbsp;<BR><BR>Here's the WKT I'm creating for the destination file: 
<BR><BR>GEOGCS["Geographic Coordinate System", <BR><BR>DATUM["WGS_1984", 
<BR><BR>SPHEROID["WGS_1984",6378137,298.257223563]], 
<BR><BR>PRIMEM["Greenwich",0], <BR><BR>UNIT["degree",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 "Geographic Projection" 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("Geographic Coordinate System",<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, "Greenwich", 0.0, "degree", <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></P></BODY></HTML>