<div dir="ltr"><div>Fortunately I didn't give too many details about this to bother you too much...</div><div><br></div><div>This is quite embarrassing as in the end I was missing the postprocessation of GEOMETRYCOLLECTION. Everything works very much as expected, very sorry for the false positive. The missing linestring was there. Thanks to st_collectionextract and friends this is easier than ever. I have warned my students at the university countless times against leaving out unprocessed GEOMETRYCOLLECTION and see what happened!</div><div><br></div><div>Thank you for answering! For topology clean-up I use to rely on GRASS most of the time. I piped the features through it too.</div><div><br></div><div>Very sorry and thank you all,<br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">







<p>---<br></p><p><span>Juan Pedro Pérez Alcántara</span></p>
<p><span><a href="mailto:jp.perez.alcantara@gmail.com" target="_blank">jp.perez.alcantara@gmail.com</a></span></p></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 3 Nov 2022 at 18:18, Martin Davis <<a href="mailto:mtnclimb@gmail.com">mtnclimb@gmail.com</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"><div dir="ltr">That certainly seems odd.  We'll need some data to investigate this further.<div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 2, 2022 at 4:37 AM Juan Pedro Pérez Alcántara <<a href="mailto:jp.perez.alcantara@gmail.com" target="_blank">jp.perez.alcantara@gmail.com</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"><div dir="ltr"><div><br>Hi, I may have found a bug in polygon intersection. I've carefully checked everything I could think of and updated all libraries to the latest version, to rule out the possibility of a GEOS error or something, but the error still persists.</div><div><br></div><div>postgis_full_version() is:</div><div><br></div><div>POSTGIS="3.3.1 3786b21" [EXTENSION] PGSQL="150" GEOS="3.11.0-CAPI-1.17.0" PROJ="9.1.0" LIBXML="2.9.13" LIBJSON="0.15" LIBPROTOBUF="1.3.3" WAGYU="0.5.0 (Internal)"</div><div><br></div><div>The system is Linux.</div><div><br></div><div>Here is a brief summary of the problem:</div><div><br></div><div><img src="cid:ii_l9zk1lxd0" alt="Screenshot from 2022-11-02 12-19-38.png" width="578" height="455"></div><div><br></div><div>As you can see in the screenshot, I have several polygons with a clean topology. A st_intersection of them, isolating only LINESTRING and MULTILINESTRING solutions, returns the green lines. As you can see, the clean topology produces clean results: all frontiers are correctly identified but the segment in black. I've thoroughly analyzed this segment, studying the as_ewkt() representation of both polygons and identifying the common primitive, they are identical coordinate by coordinate, but I'm at a loss at why st_intersection is unable to compute this exactly common primitive.</div><div><br></div><div>I've all assets available to reproduce this issue if developers think it is of interest, but I'm not going to dump them all here for brevity. Please ask me for further details if you think this is of any interest.</div><div><br></div><div>Thanks for your attention and your time,</div><div dir="ltr"><p>---<br></p><p><span>Juan Pedro Pérez Alcántara</span></p>
<p><span><a href="mailto:jp.perez.alcantara@gmail.com" target="_blank">jp.perez.alcantara@gmail.com</a></span></p></div></div>
_______________________________________________<br>
postgis-devel mailing list<br>
<a href="mailto:postgis-devel@lists.osgeo.org" target="_blank">postgis-devel@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-devel" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/postgis-devel</a><br>
</blockquote></div>
_______________________________________________<br>
postgis-devel mailing list<br>
<a href="mailto:postgis-devel@lists.osgeo.org" target="_blank">postgis-devel@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-devel" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/postgis-devel</a><br>
</blockquote></div>