Advertisement
Guest User

Untitled

a guest
Aug 19th, 2017
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.74 KB | None | 0 0
  1. create function update_reduced_geom(text, integer)
  2. returns boolean as
  3. 'DECLARE tbl ALIAS FOR $1;
  4. create temporary table tempparts as
  5. select gid, ST_Area((ST_Dump(the_geom)).geom)::bigint as part_area, (ST_Dump(the_geom)).geom
  6. from tbl;
  7. alter table tempparts add column id serial;
  8. create temporary table tempreduced as
  9. select gid, ST_Collect(geom) as reduced_geom
  10. from tempparts t1
  11. where id in
  12. (select id from county_parts t2
  13. where t1.gid = t2.gid order by part_area desc limit $2)
  14. group by gid;
  15. SELECT AddGeometryColumn (''public'', tbl, ''reduced_geom'', 40000, ''MULTIPOLYGON'',2);
  16. update tbl set reduced_geom =
  17. (select reduced_geom from county_reduced where tempreduced.gid = tbl.gid);'
  18. language 'sql';
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement