I have been attempting to intersect on set of geometries (bio.evc) with another geometry (bio.bnd_buf). All the geometries are valid according to isvalid and the query below successfully returns a table.<br><br><div style="margin-left: 40px;">
create table bio.bnd_evc as<br>select a.buf_dist, b.bio_no, b.bio_code, b.evc, b.constcde, b.consstat, b.density, intersection(a.the_geom, b.the_geom) as the_geom<br>from bio.bnd_buf as a, bio.evc as b<br>where a.the_geom
 && b.the_geom<br>and a.buf_dist = 1<br>and overlaps(a.the_geom, b.the_geom)<br><br></div>When I view the resulting table (bnd_evc) in QGIS there are several missing polygons (<a href="http://www.bendigo-orienteers.com.au/study/evc1_over.png">
image 1a</a>). One of the polygons lies entirely within bnd_bnf while the other has a small portion which overlaps the boundary.<br><br>When I try intersecting the same geometries but this time selecting buf_dist = 1000 (the 1000m buffer), additional evc polygons (which are wholly within buf_dist=1 and therefore also within buf_dist=1000) are missing from the output table and the whole thing is a 'mess' nad is in no way a representation of what the intersection shoul dlook like (
<a href="http://www.bendigo-orienteers.com.au/study/evc1000_over.png">image 2a</a>).<br><br>The other thing that occurs in both cases is that some polygon geometries are created where there are no evc polygons (<a href="http://www.bendigo-orienteers.com.au/study/evc1_under.png">
image 1b</a> & <a href="http://www.bendigo-orienteers.com.au/study/evc1000_under.png">2b</a>).<br><br>Both tables have a GIST index on the geometry column (the_geom).<br><br>Are the issues associated with having multiple polygon geometries in the one table (unlikely in my opinion)?
<br>Is it something to do with the query?<br>Are these issues associated with PostGIS or GEOS?<br>Is there a work around to solve the issue?<br><br>Thanks in advance,<br><br>Craig<br><br><br>