[postgis-tickets] [PostGIS] #3447: PG::InternalError: ERROR: Error computing Voronoi diagram.

[PostGIS]

#3447: PG::InternalError: ERROR:  Error computing Voronoi diagram.
----------------------+---------------------------
  Reporter:  loarie   |      Owner:  dbaston
      Type:  defect   |     Status:  new
  Priority:  medium   |  Milestone:  PostGIS 2.3.0
 Component:  postgis  |    Version:  trunk
Resolution:           |   Keywords:  ST_Voronoi
----------------------+---------------------------

Comment (by dbaston):

 Thanks for the example loarie.  Here are my observations so far:

 - The problem occurs within GEOS.
 - JTS has no problem producing a Voronoi diagram for this input.
 - ST_DelaunayTriangles has no problem with this input; I'd have expected
 both or none to work

 I'll keep looking around to see if I can turn up anything obvious.

 GEOS error is {{{IllegalArgumentException: Points of LinearRing do not
 form a closed linestring}}}, from

 {{{
 #0  0x00007f93f2c498b0 in __cxa_throw () from /usr/lib/x86_64-linux-
 gnu/libstdc++.so.6
 #1  0x00007f93f2fae617 in geos::geom::LinearRing::validateConstruction
 (this=this at entry=0x7f94068cbc00) at LinearRing.cpp:70
 #2  0x00007f93f2fae71b in geos::geom::LinearRing::LinearRing
 (this=0x7f94068cbc00, newCoords=0x7f94068cbbe0, newFactory=<optimized
 out>, __in_chrg=<optimized out>, __vtt_parm=<optimized out>)
     at LinearRing.cpp:46
 #3  0x00007f93f2fac115 in geos::geom::GeometryFactory::createLinearRing
 (this=this at entry=0x7f94067fa4a0, newCoords=newCoords at entry=0x7f94068cbbe0)
 at GeometryFactory.cpp:469
 #4  0x00007f93f303531b in
 geos::triangulate::quadedge::QuadEdgeSubdivision::getVoronoiCellPolygon
 (this=this at entry=0x7f940689efa0, qe=<optimized out>, geomFact=...) at
 QuadEdgeSubdivision.cpp:582
 #5  0x00007f93f3036c0e in
 geos::triangulate::quadedge::QuadEdgeSubdivision::getVoronoiCellPolygons
 (this=0x7f940689efa0, geomFact=...) at QuadEdgeSubdivision.cpp:528
 #6  0x00007f93f3036cc9 in
 geos::triangulate::quadedge::QuadEdgeSubdivision::getVoronoiDiagram
 (this=<optimized out>, geomFact=...) at QuadEdgeSubdivision.cpp:505
 #7  0x00007f93f303265b in
 geos::triangulate::VoronoiDiagramBuilder::getDiagram
 (this=this at entry=0x7ffd18db06d0, geomFact=...) at
 VoronoiDiagramBuilder.cpp:106
 #8  0x00007f93f414042e in GEOSVoronoiDiagram_r (extHandle=0x7f94067fa050,
 g1=0x7f94068bb7e0, env=0x0, tolerance=0, onlyEdges=0) at
 geos_ts_c.cpp:6519
 #9  0x00007f93f43c7fab in lwgeom_voronoi_diagram
 (g=g at entry=0x7f940683a670, env=env at entry=0x0, tolerance=tolerance at entry=0,
 output_edges=0) at lwgeom_geos.c:2013
 #10 0x00007f93f4374977 in ST_Voronoi (fcinfo=0x7f94068159e0) at
 lwgeom_geos.c:3611
 #11 0x00007f94042882f3 in ?? ()
 #12 0x00007f940428debc in ExecProject ()
 #13 0x00007f9404293f90 in ExecAgg ()
 #14 0x00007f9404287068 in ExecProcNode ()
 #15 0x00007f94042845e2 in standard_ExecutorRun ()
 #16 0x00007f9404367edf in ?? ()
 #17 0x00007f94043694aa in PortalRun ()
 #18 0x00007f9404366e65 in PostgresMain ()
 #19 0x00007f940415a280 in ?? ()
 #20 0x00007f9404321ab1 in PostmasterMain ()
 #21 0x00007f940415af03 in main ()
 }}}

--
Ticket URL: <https://trac.osgeo.org/postgis/ticket/3447#comment:4>
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.