<br><font size=2 face="sans-serif">Maybe the problem is in your data.</font>
<br><font size=2 face="sans-serif">Di you try to reproject zone 11 nad27
with the wkt nad83 below?</font>
<p><font size=3> </font>
<p>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>"Henneke, Amanda M"
<amanda.m.henneke@boeing.com></b> </font>
<p><font size=1 face="sans-serif">23/09/2009 01:44 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">"Steve.Toutant@inspq.qc.ca"
<Steve.Toutant@inspq.qc.ca></font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td><font size=1 face="sans-serif">"gdal-dev@lists.osgeo.org"
<gdal-dev@lists.osgeo.org>, "gdal-dev-bounces@lists.osgeo.org"
<gdal-dev-bounces@lists.osgeo.org></font>
<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><font size=2>PROJCS["NAD_1983_UTM_Zone_12N",</font>
<p><font size=2>GEOGCS["GCS_North_American_1983",</font>
<p><font size=2>DATUM["North_American_Datum_1983",</font>
<p><font size=2>SPHEROID["GRS_1980",6378137,298.257222101]],</font>
<p><font size=2>PRIMEM["Greenwich",0],</font>
<p><font size=2>UNIT["Degree",0.017453292519943295]],</font>
<p><font size=2>PROJECTION["Transverse_Mercator"],</font>
<p><font size=2>PARAMETER["latitude_of_origin",0],</font>
<p><font size=2>PARAMETER["central_meridian",-111],</font>
<p><font size=2>PARAMETER["scale_factor",0.9996],</font>
<p><font size=2>PARAMETER["false_easting",500000],</font>
<p><font size=2>PARAMETER["false_northing",0],</font>
<p><font size=2>UNIT["Meter",1]]</font>
<p><font size=3> </font>
<p><font size=2 color=#006291 face="Kristen ITC"><b>Amanda M. Henneke</b></font><font size=3>
</font><font size=1 color=#808080 face="Comic Sans MS"><br>
IG/DB Tools - TSS Capability Center</font><font size=3> </font><font size=1 color=#808080 face="Comic Sans MS"><br>
Boeing - St. Louis</font><font size=3> </font><font size=1 color=#808080 face="Comic Sans MS"><br>
314-777-6856</font><font size=3> </font>
<div align=center>
<p><font size=1 color=#008000 face="Comic Sans MS">When One tugs at a single
thing in Nature, he finds it attached to the rest of the world. -John Muir</font>
<p><font size=1 color=#008000 face="Comic Sans MS"> Please consider
the environment before printing this e-mail.</font></div>
<p><font size=3> </font>
<br>
<br>
<hr><font size=2 face="Tahoma"><b>From:</b> Steve.Toutant@inspq.qc.ca [mailto:Steve.Toutant@inspq.qc.ca]
<b><br>
Sent:</b> Wednesday, September 23, 2009 12:41 PM<b><br>
To:</b> Henneke, Amanda M<b><br>
Cc:</b> gdal-dev@lists.osgeo.org; gdal-dev-bounces@lists.osgeo.org<b><br>
Subject:</b> RE: [gdal-dev] Translating between UTM and Geographic</font><font size=3><br>
</font>
<br><font size=2 face="sans-serif"><br>
What is the wkt you are using for zone 12 nad83?</font><font size=3> </font>
<p><font size=3> </font>
<p><font size=3><br>
<br>
</font>
<table width=100%>
<tr valign=top>
<td width=58%><font size=1 face="sans-serif"><b>"Henneke, Amanda M"
<amanda.m.henneke@boeing.com>@lists.osgeo.org</b> <br>
Envoyé par : gdal-dev-bounces@lists.osgeo.org</font><font size=3> </font>
<p><font size=1 face="sans-serif">23/09/2009 01:39 PM</font><font size=3>
</font>
<td width=41%>
<br>
<table width=100%>
<tr valign=top>
<td width=9%>
<div align=right><font size=1 face="sans-serif">A</font></div>
<td width=90%><font size=1 face="sans-serif">"gdal-dev@lists.osgeo.org"
<gdal-dev@lists.osgeo.org></font><font size=3> </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>
<br>
<table width=100%>
<tr valign=top>
<td width=50%>
<td width=50%></table>
<br></table>
<br><font size=3><br>
<br>
<br>
<br>
</font><tt><font size=2><br>
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</font></tt><font size=3><br>
<br>
<br>
</font>
<br>
<br>
<br>