[gdal-dev] Re: Intersection of polygon layers

2010-11-05 Thread yosuke kimura
r > >>> from shapely.wkb import loads > >>> source = ogr.Open("/tmp/world_borders.shp") > >>> borders = source.GetLayerByName("world_borders") > >>> feature = borders.GetNextFeature() > >>> loads(feature.GetGeometryR

[gdal-dev] Intersection of polygon layers

2010-11-05 Thread yosuke kimura
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