Hi Yosuke,
If you use MultiPolygon geometries (with your polygons as the parts), geos
(via ogr or shapely) should be able to do the job in one go, which should
speed things up considerably...
Cheers,
Simon
On Sat, Nov 6, 2010 at 5:04 AM, yosuke kimura wrote:
> All,
>
> I think I found
All,
I think I found what I should have been using. pick a feature from
one layer, use its geometry as an argument of SetSpatialFilter method
of the other layer. Seems to be yielding same answer quicker...
Francis,
Thank you for suggestion. I looked at the Shapely documentation
http://gispyth