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 <s...@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 >
