[geos-devel] Correct use of setPrecision

Nyall Dawson nyall.dawson at gmail.com
Mon Sep 4 02:13:29 PDT 2017

On 4 September 2017 at 17:29, Sandro Santilli <strk at kbt.io> wrote:
> On Mon, Sep 04, 2017 at 04:28:32PM +1000, Nyall Dawson wrote:
>> However, after looking into this in more depth I think the original
>> python code is incorrectly used. It's being used by QGIS in order to
>> set a "tolerance" on a spatial join operation - i.e. to consider that
>> two geometries which are 100m distant actually do intersect.
> [...]
>> If this is correct, is there any other way in GEOS to handle this situation?
> How about using GEOSDistance ?

Thanks - that's a good option for intersects tests. I'm not sure
offhand if it would be translatable to the other geometry predicates,
but then I'd say trying to do a 'touches' or 'contains' test with a
tolerance would be quite an edge case anyway.


More information about the geos-devel mailing list