<div dir="ltr">Thanks Even,<div><br></div><div>From the <a href="https://www.gdal.org/classOGRSpatialReference.html#a4e0f1e23207c157f72388717c35be88c">SetStatePlane docs</a> it says</div><div>"""

<table class="gmail-params" style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:14px;line-height:22px;font-family:Roboto,sans-serif;margin-left:0px;padding-left:0px;color:rgb(0,0,0);background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><tbody><tr><td class="gmail-paramname" style="color:rgb(96,32,32);white-space:nowrap;font-weight:bold;vertical-align:top">nZone</td><td>State plane zone number, in the USGS numbering scheme (as distinct from the Arc/Info and Erdas numbering scheme.</td></tr></tbody></table>"""</div><div>which would seem to indicate that using the FIPS code in place of the USGS numbering scheme will lead to erroneous results.  It appears that my best path forward with FIPS within GDAL with python would be to read the Lib/site-packages/osgeo/data/gdal/esri_StatePlane_extra.wkt file directly and build a small converter utility.</div><div><br></div><div>Glen</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 9, 2018 at 9:02 AM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On jeudi 9 août 2018 08:42:57 CEST Glen Rice - NOAA Federal wrote:<br>
> Hello,<br>
> <br>
> I am looking for any hints for how to access ImportFromESRIStatePlaneWKT in<br>
> the OGRSpatialReference class through the python bindings.  I am working<br>
> with xyz text files where the header contains a FIPS zone reference and I<br>
> would like to convert it to another system.<br>
<br>
</span>You can't from Python. This method is only available from C++.<br>
<br>
But perhaps the osr.SetStatePlane(int zone, int is_nad83 = 1, char const <br>
*unitsname = "", double units = 0.0) that corresponds to<br>
OGRSpatialReference::<wbr>SetStatePlane() will do the job ?<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font color="#999999">Glen Rice</font></div><div dir="ltr"><font color="#999999"><span style="font-size:12.8000001907349px">Hydrographic Systems and Technology Branch</span><br></font></div><div dir="ltr"><font color="#999999">Office of Coast Survey / Coast Survey Development Laboratory</font><div><font color="#999999">603-862-1397</font></div></div></div></div></div></div></div></div>
</div>