[postgis-tickets] [PostGIS] #4274: ST_PointOnSurface generated point produces ST_Intersects = false

PostGIS trac at osgeo.org
Thu Dec 13 06:09:13 PST 2018


#4274: ST_PointOnSurface generated point produces ST_Intersects = false
---------------------+---------------------------
 Reporter:  andehhh  |      Owner:  pramsey
     Type:  defect   |     Status:  new
 Priority:  medium   |  Milestone:  PostGIS 2.5.2
Component:  postgis  |    Version:  2.4.x
 Keywords:           |
---------------------+---------------------------
 Given is this query which returns 'false':

 {{{
 SELECT ST_Intersects(ST_PointOnSurface(poly), poly) FROM (SELECT
 'SRID=4326; POLYGON ((2.0875049999999997 49.038467, 2.087497642793396
 49.0384427212182, 2.087475 49.038368, 2.0875049999999997
 49.038467))'::geometry AS poly) AS foo;
 }}}


 Looking at the documentation ("Returns a POINT **guaranteed** to intersect
 a surface.") this case should be handled internally to avoid confusion.

 Can it be explained with rounding issues as the polygon in question is a
 very narrow sliver polygon?

 When transforming to a metric system this issue disappears.

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/4274>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.


More information about the postgis-tickets mailing list