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