Skip to content

Instantly share code, notes, and snippets.

@ian29
Last active December 23, 2015 21:09
Show Gist options
  • Save ian29/6694803 to your computer and use it in GitHub Desktop.
Save ian29/6694803 to your computer and use it in GitHub Desktop.
create table test2 (
gid int,
qs_adm0 text
);
SELECT AddGeometryColumn('public', 'test2', 'geom', 4326, 'MULTIPOLYGON', 2);
insert into test2 (
SELECT
gid,
qs_adm0,
ST_Multi(
St_Union(ering_geom)
) AS geom
FROM (
SELECT
gid,
qs_adm0,
ST_MakePolygon(
ST_ExteriorRing(dump_geom)
) AS ering_geom
FROM (
SELECT
qs.gid,
qs_adm0,
(ST_Dump(
st_union(qs.geom,eez.geom)
)).geom AS dump_geom
FROM qs LEFT JOIN eez
ON qs_adm0_a3 = iso_3digit
WHERE st_x(st_centroid(qs.geom)) < 0 -- adjust for hemisphers
-- LIMIT 20
) AS foo
GROUP BY dump_geom,gid,qs_adm0
) as bar
GROUP BY gid,qs_adm0
);
@ian29
Copy link
Author

ian29 commented Sep 25, 2013

ow my brain

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment