<div dir="ltr">Hi,<div><br></div><div>I am running the current setup:</div><div> "POSTGIS="2.0.2 r10789" GEOS="3.4.0dev-CAPI-1.8.0" PROJ="Rel. 4.7.1, 23 September 2009" GDAL="GDAL 1.9.0, released 2011/12/29" LIBXML="2.7.8" (core procs from "2.0.1 r9979" need upgrade) RASTER (raster procs from "2.0.1 r9979" need upgrade)"</div>
<div><br></div><div>I have two tables with geometries. Both include polygons representing administrative units. One at the first level and one at the second level.</div><div>I want a table that includes all the second level polygons where available, and where not available I want the first level polygons.</div>
<div><br></div><div>To find the level 1 polygons for filling the gaps where level 2 is N/A, I thought about finding where the centroid of level one polygons does not intersects with the polygons for level 2.</div><div><br>
</div><div>Hence, I tried:</div><div><div><br></div><div>SELECT <a href="http://g1.name">g1.name</a>, g1.geom INTO gadm2_level1_union_selected FROM gadm2_level1_union g1, gadm2_level2_union g2 WHERE ST_Intersects(g1.centroid, g2.geom)=FALSE;</div>
<div><br></div><div>I also tried the ST_Disjoint(g1.centroid, g2.geom).</div><div><br></div><div>My problem is that this is running for ages without finishing. This should be such a simple query, but I have used half of the day running it. If I take the two geometries into PostGIS and run a spatial query to find the points from level 1 not intersecting with polygons of level 2, it takes 15 seconds.</div>
<div><br></div><div>What could be the issue with my query?</div><div><br></div><div>Thanks!</div><div>Andreas</div><div><br></div></div></div>