[gdal-dev] Tests if the geometry is within the other geometry

Ian gisplaybox at gmail.com
Fri Jun 22 23:01:43 PDT 2018


inlay.SetSpatialFilter(refpoy) did not meet my requirement of  "if
geomi.Within(refpoly)".

if I use both  inlay.SetSpatialFilter(refpoy) and  geomi.Within(refpoly),
speed is same as before.
.SetSpatialFilter tested intersect of the envelopes, which i do not need.
any idea is appreciated

On Sat, Jun 23, 2018 at 2:01 PM, Ian <gisplaybox at gmail.com> wrote:

> after spatial filter, should i use  if geomi.Within(refpoly)?
>
> On Sat, Jun 23, 2018 at 4:49 AM, mikeucfl <mikeucfl at gmail.com> wrote:
>
>> 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
>> _______________________________________________
>> gdal-dev mailing list
>> gdal-dev at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20180623/846697f5/attachment.html>


More information about the gdal-dev mailing list