Use a spatial filter using your refpoly before iterating the features: http://gdal.org/python/osgeo.ogr.Layer-class.html#SetSpatialFilter -- Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html