<div dir="ltr">Discussion here: <a href="https://github.com/libgeos/geos/issues/968#issuecomment-2041134101">https://github.com/libgeos/geos/issues/968#issuecomment-2041134101</a><div> </div><div>This could also be raised as a separate issue.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 25, 2024 at 2:07 PM Sandro Santilli <<a href="mailto:strk@kbt.io">strk@kbt.io</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I'm struggling with a robustness bug in PostGIS topology [1] which<br>
is making me question the relate computer of GEOS. Look at this:<br>
<br>
  E=010200000002000000B02EA9C2AE5430400E3994256830514037F136E82055304039FAE86851305140<br>
  I=01020000000200000064AB380D714E3040A38E9F0E1F305140E798A2D4E85430405F5969945C305140<br>
<br>
  geosop -a $E -b $I relate<br>
  # F01FF0102<br>
<br>
  # Interior of $E intersects the boundary if $I<br>
  #    I B E<br>
  #  I F 0 1<br>
  #  B F F 0<br>
  #  E 1 0 2  <br>
<br>
  geosop -a $E -b $(geosop -f wkb -a $I boundary) relate<br>
  # FF1FF00F2<br>
<br>
  # Boundary of $I has no intersection with interior of $E<br>
  #    I B E<br>
  #  I F F 1<br>
  #  B F F 0<br>
  #  E 0 F 2  <br>
<br>
<br>
So, where is that "boundary" intersection with 'i' ?<br>
<br>
[1] <a href="https://trac.osgeo.org/postgis/ticket/5699" rel="noreferrer" target="_blank">https://trac.osgeo.org/postgis/ticket/5699</a><br>
<br>
--strk; <br>
<br>
  Libre GIS consultant/developer<br>
  <a href="https://strk.kbt.io/services.html" rel="noreferrer" target="_blank">https://strk.kbt.io/services.html</a><br>
</blockquote></div>