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
signature.asc
Description: PGP signature