I need to check "points in polygon" over various thousands of points and regional scale polygons. A for loop iteration + intersect() over the points geometries is absolutely inefficient. Is there a faster method? I'm on the way to test the Shpely iterators, but I would like to remain with OGR... thanks again, giovanni