<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">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>