<!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> </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 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> Please consider the environment before printing this
e-mail.</FONT></SPAN></P>
<DIV> </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> </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"
<amanda.m.henneke@boeing.com>@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"
<gdal-dev@lists.osgeo.org></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> <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> <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> <BR><BR>Do you see any issues with that? 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> 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> 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. What
Datum is the data in the "Geographic Projection" in? 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. 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>> -----Original Message-----<BR>>
From: Henneke, Amanda M [mailto:amanda.m.henneke@boeing.com]<BR>> Sent:
Wednesday, September 23, 2009 9:28 AM<BR>> To: Craig Miller<BR>> Subject:
RE: [gdal-dev] Translating between UTM and Geographic<BR>> <BR>> Hi
Craig-<BR>> <BR>> I'm translating a shapefile (from UTM to Geographic) and
then <BR>> comparing it<BR>to<BR>> <BR>> #1 A shapefile that was
translated from UTM to Geographic using the<BR>software,<BR>>
GlobalMapper.<BR>> #2 A raster file that was originally projected in
Geographic, so that <BR>> for example I can see the points don't match up
with the features they <BR>> should represent on the image.<BR>> <BR>>
The shapefile that GlobalMapper translates matches up perfectly with <BR>>
my<BR>raster<BR>> file.<BR>> <BR>> So my thought process is, that if
Global Mapper has some way to <BR>> convert<BR>the<BR>> different
projections without the shifting, how come GDAL can't do that.<BR>I'm<BR>>
thinking it can, but maybe I need to give it some additional data
maybe?<BR>Or<BR>> am I off?<BR>> <BR>> Thanks for your help!<BR>>
<BR>> Amanda M. Henneke<BR>> IG/DB Tools - TSS Capability Center<BR>>
<BR>> <BR>> When One tugs at a single thing in Nature, he finds it
attached to the<BR>rest of<BR>> the world. -John Muir<BR>> 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:21 AM<BR>> To: Henneke, Amanda M<BR>> Subject: RE: [gdal-dev]
Translating between UTM and Geographic<BR>> <BR>> If you are transforming
the datum from NAD27 to WGS84 then it is <BR>> expected<BR>that<BR>> your
points would be shifted exactly as you describe. How are
you<BR>determining<BR>> that they are shifted?<BR>> <BR>> a.
Compared to a map? Is the map reprojected? Are you comparing
<BR>> them to landform features or to the lat/lon grid on the map? For
<BR>> example, if you<BR>have<BR>> a USGS Quad map with lat/lon lines on
it based on the NAD27 datum then<BR>those<BR>> lat/lon lines ONLY show the
lat/lon using that datum. If you plot <BR>> points<BR>that<BR>> have
been shifted to the new datum, they will not agree with those<BR>>
NAD27<BR>lines<BR>> since the new points are based on lat/lon NAD83.<BR>>
<BR>> b. Comparing lat/lon directly? The points have been
transformed.<BR>Lat/lon<BR>> are relative to a datum and lan/lon NAD27 coords
and lat/lon NAD83 <BR>> coords<BR>will<BR>> have different values.<BR>>
<BR>> Craig<BR>> <BR>> <BR>> > -----Original Message-----<BR>>
> From: gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev- <BR>> >
bounces@lists.osgeo.org] On Behalf Of Henneke, Amanda M<BR>> > Sent:
Wednesday, September 23, 2009 8:17 AM<BR>> > To:
gdal-dev@lists.osgeo.org<BR>> > Subject: [gdal-dev] Translating between
UTM and Geographic<BR>> ><BR>> > Hello all-<BR>> ><BR>>
> I'm trying to figure out how to reproject some points from a UTM<BR>>
projection to<BR>> > Geographic (WGS84) projection. This works fine
when I'm working in <BR>> > UTM<BR>> Zone<BR>> > 12 with NAD83.
When I try the same code for Zone 11 NAD27, I get a <BR>> > slight
shift in the points (to the east) and an extremely slight <BR>> > shift up
(to the north).<BR>> ><BR>> > I've included the following code as it
stands now, I've tried all <BR>> > sorts of things thinking that I need to
somehow include more <BR>> > information for the<BR>> new<BR>> >
spatial reference (the destination SR), and so far I'm coming
up<BR>empty...<BR>> ><BR>> > Anyone have any ideas what I'm
missing?<BR>> ><BR>> > Where sSpRefParams is a string that contains
the WKT projection<BR>> information<BR>> > from a valid
layer...<BR>> ><BR>> > srOld = new
OSGeo.OSR.SpatialReference(sSpRefParams);<BR>> ><BR>> >
srNew.SetGeogCS("Geographic Coordinate System",<BR>>
OSGeo.OSR.Osr.SRS_DN_WGS84,<BR>> > OSGeo.OSR.Osr.SRS_DN_WGS84,
OSGeo.OSR.Osr.SRS_WGS84_SEMIMAJOR,<BR>> >
OSGeo.OSR.Osr.SRS_WGS84_INVFLATTENING, "Greenwich", 0.0, "degree", <BR>> >
0.0174532925199433);<BR>> ><BR>> > ctTransform = new
OSGeo.OSR.CoordinateTransformation(srOld, srNew);<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 <BR>> > the<BR>>
rest of<BR>> > the world. -John Muir<BR>> > Please consider
the environment before printing this e-mail.<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><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>