[postgis-tickets] [PostGIS] #3886: ST_Subdivide loses boxes on some inputs

[PostGIS]

#3886: ST_Subdivide loses boxes on some inputs
---------------------+---------------------------
 Reporter:  komzpa   |      Owner:  pramsey
     Type:  defect   |     Status:  new
 Priority:  medium   |  Milestone:  PostGIS 2.4.1
Component:  postgis  |    Version:  2.4.x
 Keywords:           |
---------------------+---------------------------
 There is landmass geometry that ST_Subdivide fails to correctly split. It
 loses a whole big box, producing a warning:


 {{{
 [local] gis at gis=# select ST_IsValid(geom) from bad_geom;
 ┌────────────┐
 │ st_isvalid │
 ├────────────┤
 │ t          │
 └────────────┘
 (1 row)

 Time: 686,142 ms
 [local] gis at gis=# select ST_IsSimple(geom) from bad_geom;
 ┌─────────────┐
 │ st_issimple │
 ├─────────────┤
 │ t           │
 └─────────────┘
 (1 row)

 Time: 37,460 ms
 [local] gis at gis=# select ST_Subdivide(geom) from bad_geom limit 1;
 NOTICE:  00000: Error performing rectangular clipping:
 IllegalArgumentException: Invalid number of points in LinearRing found 3 -
 must be 0 or >= 4
 LOCATION:  pg_notice, lwgeom_pg.c:184
 ┌────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────┐
 │
 st_subdivide
 │
 ├────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────┤
 │
 0103000020110F00000100000005000000DB916D3B01A260C100000060C3CE5141000000F00C9F60C100000060C3CE5141000000F00C9F60C1E7B8076F2BC85141F6285CC78CA160C148E17A94BFCD5141DB916D3B01A260C100000060C3CE5141
 │
 └────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────┘
 (1 row)

 Time: 2070,714 ms (00:02,071)

 }}}

--
Ticket URL: <https://trac.osgeo.org/postgis/ticket/3886>
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.