Relate computer inconsistencies ?

[Martin Davis]

Discussion here:
https://github.com/libgeos/geos/issues/968#issuecomment-2041134101

This could also be raised as a separate issue.

On Mon, Mar 25, 2024 at 2:07 PM Sandro Santilli <strk at kbt.io> wrote:

> I'm struggling with a robustness bug in PostGIS topology [1] which
> is making me question the relate computer of GEOS. Look at this:
>
>
> E=010200000002000000B02EA9C2AE5430400E3994256830514037F136E82055304039FAE86851305140
>
> I=01020000000200000064AB380D714E3040A38E9F0E1F305140E798A2D4E85430405F5969945C305140
>
>   geosop -a $E -b $I relate
>   # F01FF0102
>
>   # Interior of $E intersects the boundary if $I
>   #    I B E
>   #  I F 0 1
>   #  B F F 0
>   #  E 1 0 2
>
>   geosop -a $E -b $(geosop -f wkb -a $I boundary) relate
>   # FF1FF00F2
>
>   # Boundary of $I has no intersection with interior of $E
>   #    I B E
>   #  I F F 1
>   #  B F F 0
>   #  E 0 F 2
>
>
> So, where is that "boundary" intersection with 'i' ?
>
> [1] https://trac.osgeo.org/postgis/ticket/5699
>
> --strk;
>
>   Libre GIS consultant/developer
>   https://strk.kbt.io/services.html
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/geos-devel/attachments/20240406/3d5d5e34/attachment.htm>