[gdal-dev] ImportFromESRIStatePlaneWKT access from python

[Even Rouault]

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 <even.rouault at spatialys.com>
> 
> 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