[postgis-devel] [PostGIS] #2055: Intersection with a result from Intersection fails

[PostGIS]

#2055: Intersection with a result from Intersection fails
--------------------------+-------------------------------------------------
 Reporter:  bscholten     |       Owner:  pramsey      
     Type:  defect        |      Status:  new          
 Priority:  medium        |   Milestone:  PostGIS 1.5.6
Component:  postgis       |     Version:  1.5.X        
 Keywords:  Intersection  |  
--------------------------+-------------------------------------------------
 Basic problem:

 Following query gives an error:
 ----
 SELECT ST_Intersection(ST_Intersection(first_column, second_column),
 third_column) FROM test_bug_situation
 ----
 The error returned is:
 ----
 NOTICE:  TopologyException: no outgoing dirEdge found at
 184311.56259997687 319320.99999998196


 ERROR:  GEOS Intersection() threw an error!

 ********** Fout **********

 ERROR: GEOS Intersection() threw an error!
 SQL status:XX000
 ----
 ST_IsValid and ST_IsSimple both return true for all columns and
 ST_Intersection(first_column, second_column).

 Switching the columns doesn't throw an error:
 ST_Intersection(ST_Intersection(first_column, third_column),
 second_column)
 or
 ST_Intersection(ST_Intersection(third_column, second_column),
 first_column)

 Extra information:
 ----
 Postgres version:
 "PostgreSQL 8.4.12, compiled by Visual C++ build 1400, 32-bit"
 Postgis full version:
 "POSTGIS="1.5.5" GEOS="3.3.5-CAPI-1.7.5" PROJ="Rel. 4.6.1, 21 August 2008"
 LIBXML="2.7.8" USE_STATS"
 Tested on windows platform.

 See accompanying bug_report_XYZ_column_base64 for each geometry encoded in
 base64 Encoded WKB.

 Same test with POSTGIS 1.4.2 on both Windows and Linux has same results.

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