On jeudi 9 août 2018 09:27:43 CEST Glen Rice - NOAA Federal wrote: > Thanks Even, > > From the SetStatePlane docs > <https://www.gdal.org/classOGRSpatialReference.html#a4e0f1e23207c157f7238871 > 7c35be88c> 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.
The mapping from USGS (left column) to ESRI (right column) zone number is at: https://github.com/OSGeo/gdal/blob/master/gdal/ogr/ogr_srs_esri.cpp#L134 and ultimately this reads the data/stateplane.csv with the USGS zone number as a key. > 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. Sounds reasonable > > 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 -- Spatialys - Geospatial professional services http://www.spatialys.com _______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
