[geos-devel] Correct use of setPrecision

[Nyall Dawson]

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.

Nyall