[postgis-users] Adjecent Polygons to lines, typology related

Peter webwiz at pl.net
Thu Mar 22 16:58:56 PDT 2012

> Its a pity that st_intersection cant work on sets then it could be used to
>  return linestrings where adjacent polygons intersect.

Well upon reflection, it can, kind of (to answer my own question):

CREATE TABLE sharedlines AS
   ST_intersection(a.the_geom, b.the_geom) AS the_geom
FROM polys a, polys b
where ST_touches(a.the_geom, b.the_geom) and a.gid!=b.gid;

> Probably slow though.

Surpisingly not. Much, much, quicker than st_union. st_union has given me quite a bit of grief lately, eg see

Ill add the above to it later.


> Am i wasting time working on this or should i wait for postgis 2.0?
> If the latter is their a topology primer somewhere?
> Peter
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users

More information about the postgis-users mailing list