[gdal-dev] Python interface to wrapdateline

[Even Rouault]

Le mardi 10 mai 2016 16:31:44, DeTracey, Brendan a écrit :
> Hi,
> 
> Is there any way to access the wrap dateline capabilities of OGR from the
> Python interface? For example how to split the following polygon:
> 
> from osgeo import ogr, osr
> srs = osr.SpatialReference()
> srs.ImportFromEPSG(4326)
> p = ogr.CreateGeometryFromWkt('POLYGON ((160 1, -160 1, -160 -1, 160 -1,
> 160 1))')
> 
> 
> Looking at the OGRGeometryFactory  source I see :
> static void CutGeometryOnDateLineAndAddToMulti(OGRGeometryCollection*
> poMulti, const OGRGeometry* poGeom, double dfDateLineOffset)
> 
> but I can't see how to access it from Python. Would it be via an option to
> the ogr.CreateGeometryFromWkt?

No, this is unavailable through Python bindings.

With GDAL 2.1, you can use gdal.VectorTranslate() (librarified version of 
ogr2ogr), although it might not be very practical if you just need to do that 
on a single geometry.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com