Hi Saranya, 2017-02-20 5:57 GMT+01:00 Saranya Kuniyil <saranyak...@gmail.com>: > >> Hi all, >> I am new to postgis and I need to find intersection of three polygons. >> In my table named "mytable" contains large no of polygon coordinates and >> grouped into three categories by adding a field called "filename".I >> labelled the category names as part1, part2, part3 in filename . Now I need >> to find all polygon which spacially intersect in three categories. >> I write my code as: >> SELECT ST_Intersection(a.geom, b.geom) as intersect_ab >> FROM mytable a INNER JOIN mytable b ON ST_Intersects(a.geom,b.geom) INNER >> JOIN mytable c ON ST_Intersects(a.geom,c.geom) >> WHERE ST_Overlaps(a.geom, b.geom) AND ST_Overlaps(b.geom, c.geom) AND >> ST_Overlaps(a.geom, c.geom) AND ST_isvalid(a.geom)='t' AND >> ST_isvalid(b.geom)='t' AND ST_isvalid(c.geom)='t' AND a.filename = 'part1' >> AND b.filename = 'part2' AND c.filename = 'part3') ; >> Is it correct????? >> > In principle, your query is fine, but are you interested to retrieve just the three intersecting geometries or the resulting intersection of these (currently, you are selecting the intersection between part1 and part2 geometries)? Moreover, are you interested on overlapping or intersection (I mean, is it fine if a part1 geometry contains a part2 or part3 ones)? I'm asking this because the intersects computation is faster than the overlaps one.
Regards, Giuseppe. -- Giuseppe Broccolo - 2ndQuadrant Italy PostgreSQL & PostGIS Training, Services and Support giuseppe.brocc...@2ndquadrant.it | www.2ndQuadrant.it
_______________________________________________ postgis-users mailing list postgis-users@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/postgis-users
