<div dir="ltr"><div><div>Fantastic! Adding the the "+wktext" argument makes ogr2ogr generate the values I was expecting. I will rebuild my copy of gdal with the fixes in the ticked mentioned.<br><br></div>Thanks again,<br>
</div> -Jay <br><div><div> </div></div></div><div class="gmail_extra"><br clear="all"><div>- Jay Cable - <a href="mailto:jay@alaska.edu">jay@alaska.edu</a> - (907)474-1856 -<br>- Geographic Information Network of Alaska - University of Alaska -</div>
<br><br><div class="gmail_quote">On Fri, Jun 6, 2014 at 11:18 AM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org" target="_blank">even.rouault@mines-paris.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Jay,<br>
<br>
I've reproduced and fixed your issue. See<br>
<a href="http://trac.osgeo.org/gdal/ticket/5511" target="_blank">http://trac.osgeo.org/gdal/ticket/5511</a> for details<br>
<br>
An easy workaround is to add " +wktext" at the end of your proj.4 string. This<br>
will instruct the GDAL SpatialReference engine to keep literaly your proj.4<br>
string to use it later. Otherwise the original proj.4 is converted into WKT<br>
and not preserved, so when needed it afterwards, GDAL rebuilds it from WKT.<br>
<br>
Best regards,<br>
<br>
Even<br>
<div class="HOEnZb"><div class="h5"><br>
> Hello,<br>
><br>
> Please forgive me if this is the wrong place for this question, but I am<br>
> trying to use Alaska State Plane Zone 1 with ogr2ogr, and can't seem to get<br>
> it to project the data correctly. I am probably doing something wrong, but<br>
> can't seem to figure it out. I am pretty sure that my projection string<br>
> is correct, as the results from proj match what the NGS's online calculator<br>
> gives.<br>
><br>
> When I use ogr2ogr like this:<br>
> ogr2ogr -t_srs "+proj=omerc +lat_0=57 +lonc=-133.6666666666667<br>
> +alpha=-36.86989764583333 +k=0.9999 +x_0=5000000 +y_0=-5000000 +no_uoff<br>
> +datum=NAD83 +units=m +no_defs" akzone1.zhp point.84.shp<br>
><br>
> On a shapefile with a single point - POINT (-131.811033901222117<br>
> 55.409950489560302)<br>
><br>
> I get:<br>
> POINT (4988757.639998742379248 -5210867.590697132050991)<br>
><br>
> What I get from a shape file reprojected with ESRI products is:<br>
> POINT (936203.399999997578561 399658.200000001117587)<br>
><br>
> I believe the ESRI result is correct, as it matches the the NGA online<br>
> converter.<br>
><br>
> Proj also generates the correct value:<br>
> Jays-MacBook-Air-2:zone1 Jay$ proj +proj=omerc +lat_0=57<br>
> +lonc=-133.6666666666667 +alpha=-36.86989764583333 +k=0.9999 +x_0=5000000<br>
> +y_0=-5000000 +no_uoff +datum=NAD83 +units=m +no_defs<br>
> -131.811033901222117 55.409950489560302<br>
> 936203.40 399658.20<br>
><br>
><br>
> So, my question is what am I missing to make the results from ogr2ogr match<br>
> up with proj?<br>
><br>
> My thanks,<br>
> -Jay Cable<br>
><br>
> Jays-MacBook-Air-2:zone1 Jay$ proj<br>
> Rel. 4.8.0, 6 March 2012<br>
> usage: proj [ -beEfiIlormsStTvVwW [args] ] [ +opts[=arg] ] [ files ]<br>
><br>
> Jays-MacBook-Air-2:zone1 Jay$ ogr2ogr --version<br>
> GDAL 1.11.0, released 2014/04/16<br>
><br>
><br>
> Oh, the online NGS converter I mentioned:<br>
> <a href="http://www.ngs.noaa.gov/cgi-bin/spc_getpc.prl" target="_blank">http://www.ngs.noaa.gov/cgi-bin/spc_getpc.prl</a><br>
><br>
> =======================================================<br>
> Latitude Longitude Datum Zone<br>
> INPUT = N552435.8217 W1314839.7220 NAD83 5001<br>
> =======================================================<br>
><br>
> NORTH(Y) EAST(X) AREA CONVERGENCE SCALE<br>
> METERS METERS DD MM SS.ss<br>
> ----------- ----------- ---- ----------- ----------<br>
> 399658.198 936203.401 AK 1 1 32 38.22 0.99990159<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> - Jay Cable - <a href="mailto:jay@alaska.edu">jay@alaska.edu</a> - <a href="tel:%28907%29474-1856" value="+19074741856">(907)474-1856</a> -<br>
> - Geographic Information Network of Alaska - University of Alaska -<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Geospatial professional services<br>
<a href="http://even.rouault.free.fr/services.html" target="_blank">http://even.rouault.free.fr/services.html</a><br>
</font></span></blockquote></div><br></div>