Changes between Version 7 and Version 8 of UsersWikiSimplifyPreserveTopology


Ignore:
Timestamp:
Apr 8, 2012, 7:26:16 AM (12 years ago)
Author:
nicolasribotosgeo
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • UsersWikiSimplifyPreserveTopology

    v7 v8  
    148148        select gid, code_dept, (st_dump(geom)).*
    149149        from departement
    150 ) select distinct on (d.code_dept) d.code_dept, qux.geom
     150) select distinct on (d.code_dept) d.code_dept, baz.geom
    151151 from (
    152         select (st_dump(st_polygonize(distinct baz.geom))).geom as geom
     152        select (st_dump(st_polygonize(distinct geom))).geom as geom
    153153        from (
    154                 select (st_dump(geom)).geom as geom
     154                select (st_dump(st_simplifyPreserveTopology(st_linemerge(st_union(geom)), 10000))).geom as geom
    155155                from (
    156                         select st_simplifyPreserveTopology(st_linemerge(st_union(geom)), 10000) as geom
    157                         from (
    158                                 select st_exteriorRing((st_dumpRings(geom)).geom) as geom
    159                                 from poly
    160                         ) as foo
    161                 ) as bar
    162         ) as baz
    163 ) as qux,
     156                        select st_exteriorRing((st_dumpRings(geom)).geom) as geom
     157                        from poly
     158                ) as foo
     159        ) as bar
     160) as baz,
    164161poly d
    165 where st_intersects(d.geom, qux.geom)
    166 order by d.code_dept, st_area(st_intersection(qux.geom, d.geom))/st_area(qux.geom) desc;
     162where st_intersects(d.geom, baz.geom)
     163order by d.code_dept, st_area(st_intersection(baz.geom, d.geom))/st_area(baz.geom) desc;
    167164}}}
    168165