[postgis-devel] [PostGIS] #1735: ST_MakeValid: exception on geometry

[PostGIS]

#1735: ST_MakeValid: exception on geometry
-----------------------+----------------------------------------------------
 Reporter:  aperi2007  |       Owner:  strk         
     Type:  defect     |      Status:  assigned     
 Priority:  medium     |   Milestone:  PostGIS 2.0.1
Component:  postgis    |     Version:  trunk        
 Keywords:             |  
-----------------------+----------------------------------------------------

Comment(by strk):

 The exception is caused by an invalid input passed to GEOSSymDifference,
 but the invalid input was the output of a previous GEOSSymDifference call
 !

 Check this out:

 {{{
 with inp as ( select
 '010300000001000000080000000000000087523A41000000402C62524100000000B0523A41000000C01E6252410000000084523A41000000C023625241858C4F2488523A4151F8EEAC2562524111A0F52288523A414F50F3AC25625241D59FF52288523A414F50F3AC25625241000000007F523A4100000080256252410000000087523A41000000402C625241'::geometry
 as a ,
 '01030000000100000004000000A871502388523A4158FEF2AC2562524100A0F52288523A414F50F3AC2562524111A0F52288523A414F50F3AC25625241A871502388523A4158FEF2AC25625241'::geometry
 as b ) select ST_IsValid(a) as a_valid, ST_IsValid(b) as b_valid,
 ST_IsValid(ST_SymDifference(a,b)) as sd_valid from inp;

 NOTICE:  Self-intersection at or near point 1725064.1365604463
 4819094.70235069
  a_valid | b_valid | sd_valid
 ---------+---------+----------
  t       | t       | f
 (1 row)
 }}}

 Good material for another ticket (on GEOS).
 Will file and report back!

-- 
Ticket URL: <http://trac.osgeo.org/postgis/ticket/1735#comment:10>
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.