[geos-devel] [GEOS] #374: side location conflict GEOS touches() threw an error

[GEOS]

#374: side location conflict GEOS touches() threw an error
------------------------+---------------------------------------------------
 Reporter:  davideps    |       Owner:  geos-devel@…              
     Type:  defect      |      Status:  new                       
 Priority:  minor       |   Milestone:                            
Component:  Default     |     Version:  3.0.3                     
 Severity:  Unassigned  |    Keywords:                            
------------------------+---------------------------------------------------
 Hello,

 I reported this on the postGIS list, but realize now that the error is
 with this underlying component. I thought I might get some advice here. I
 have a postGIS polygon layer with about 320,000 rows. I am trying to
 return a table with ogc_fid and the average area of the surrounding
 (touching) polygons. I am getting this error:


 NOTICE:  TopologyException: side location conflict 1.34593e+07 289017
 ERROR:  GEOS touches() threw an error!

 ********** Error **********

 ERROR: GEOS touches() threw an error!
 SQL state: XX000


 I have tried to simplify the geometry many different times. Nothing
 seems to solve the problem. I can process the first 20,000 rows
 (object_number < 20000) but something is causing problems after that. I
 have pasted the script below, any suggestions would be welcomed! I *may*
 be able to post some of the data, but I will need to get permission from
 the owner.

 Thank you,
 -david


 SELECT p1.ogc_fid, avg(p2.shape_area)
 FROM parcels09_d3_v2_simplify as p1, parcels09_d3_v2_simplify as p2
 WHERE st_touches(p1.wkb_geometry, p2.wkb_geometry)
 GROUP BY p1.ogc_fid
 ORDER BY p1.ogc_fid

-- 
Ticket URL: <http://trac.osgeo.org/geos/ticket/374>
GEOS <http://geos.refractions.net/>
GEOS (Geometry Engine - Open Source) is a C++ port of the Java Topology Suite (JTS).