<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hi,<br><br>I'm trying to determine all the possible intersections in sets of polygons.<br><br>The approach I'm using is to load each set of polygon perimeters into a Postgis topology & extract all the faces as geometries. This results in some new polygons lying between the input ones. Which is the nature of this technique- that is not the problem.<br><br>I then determine the number of original polygons that the new ones overlap, which should go from 0 up. My problem arises in that I'm getting new polygons, derived from the boundaries of the old ones, which are actually between old polygons, but have a non-zero (but very small - like 1e-15) area of overlap - but all they have in common is the boundary line which should have a zero area.<br><br>So some of my new polygons supposedly overlap an old polygon, but in fact they do (or should) not. <br><br>My
 solution at the moment is to reject (treat as a zero area) any overlaps < 1e-15 in area, so the number of overlaps represents the "real" data. It works for my purposes, but is not ideal.<br><br>Is this a bug? An inevitable result with near zero values & finite precision?<br>Is there a better way of dealing with this issue?  <br><br>Thanks,<br><br>    Brent Wood<br></td></tr></table>