<div dir="ltr"><div><div>Hi,</div><div>I'm having an issue testing the intersection between a point and two curved polygons.</div><div>Here is the situation visually, from QGIS:</div><div><br></div><div><img src="cid:ii_m4wn1atd0" alt="image.png" width="527" height="194" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"><br></div><div><br></div><div>and zooming in (a lot):</div><div><br></div><div><img src="cid:ii_m4wn1p1m1" alt="image.png" width="402" height="411" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"><br></div><div><br></div><div>In particular, the point is (barely) inside the polygon at the top, whose identifier is 12875.</div><div>Now if I run an intersection test, I get both polygons as a result:</div><div><br></div><div>> SELECT ogc_fid FROM testdata WHERE ST_Intersects(geom, ST_GeomFromText('POINT (25492818 6677399.98)', 3879));</div> ogc_fid<br>---------<br>    1258<br>   12875<br>(2 rows)<br><div><br></div><div>As you can see... both polygons tested positive for intersection. Visually that does not make sense, the point is in a hole inside polygon 1258, but maybe there is tolerance/linearization in play, or it's just a visualization fluke (how the curved polygons are rendered)... but I've also visualized the same data in GeoServer, and the visual result is the same.</div><div><br></div><div>I have verified the behavior in the following two version:</div><div><ul><li style="margin-left:15px">postgis 3.2 and postgresql 14.15</li><li style="margin-left:15px">postgis 3.5 and postgresql 17.2 (latest postgis docker image)</li></ul><div>Test data can be found in <a href="https://gist.github.com/aaime/dcf55b635c5bd09c2b9ebb31a2ab083a" target="_blank">this Github GIST</a>.</div></div><div><br></div><div>Can anyone shed some light on the behavior I'm seeing?</div><div><br></div><div>Best regards</div><font color="#888888"><div>Andrea</div></font></div></div>