Thanks Even, >From the SetStatePlane docs <https://www.gdal.org/classOGRSpatialReference.html#a4e0f1e23207c157f72388717c35be88c> it says """ nZone State plane zone number, in the USGS numbering scheme (as distinct from the Arc/Info and Erdas numbering scheme.""" 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.
Glen On Thu, Aug 9, 2018 at 9:02 AM, Even Rouault <[email protected]> wrote: > On jeudi 9 août 2018 08:42:57 CEST Glen Rice - NOAA Federal wrote: > > Hello, > > > > I am looking for any hints for how to access ImportFromESRIStatePlaneWKT > in > > the OGRSpatialReference class through the python bindings. I am working > > with xyz text files where the header contains a FIPS zone reference and I > > would like to convert it to another system. > > You can't from Python. This method is only available from C++. > > But perhaps the osr.SetStatePlane(int zone, int is_nad83 = 1, char const > *unitsname = "", double units = 0.0) that corresponds to > OGRSpatialReference::SetStatePlane() will do the job ? > > Even > > -- > Spatialys - Geospatial professional services > http://www.spatialys.com > -- Glen Rice Hydrographic Systems and Technology Branch Office of Coast Survey / Coast Survey Development Laboratory 603-862-1397
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
