r
> >>> from shapely.wkb import loads
> >>> source = ogr.Open("/tmp/world_borders.shp")
> >>> borders = source.GetLayerByName("world_borders")
> >>> feature = borders.GetNextFeature()
> >>> loads(feature.GetGeometryR
easier to test intersect if there are only handful of vertices as
oppose to hundreds of vertices in my case).
Right now I am trying with layer of hundreds of thousands of polygons
vs. thousands of polygon, and it seems to be taking 10s of hours to
finish in my environment.
Thank you,
Yosuke Kimura