8 | | {{{ |
9 | | #!sql |
10 | | select * from european_wise2008.rbd_f1v3 where name_engl like 'Meuse' |
11 | | |
12 | | |
13 | | -- selection of the meuse basin |
14 | | select * from european_wise2008.rbd_f1v3 where name_engl like 'Meuse' |
15 | | -- 5 lines |
16 | | |
17 | | |
18 | | select * from europe.wso limit 10;-- id, wso_id, area |
19 | | -- I just need the two second columns, below a request spatial joining riversegments and the wise layer |
20 | | select distinct on (wso_id) wso_id, name_engl as area from ccm21.riversegments r |
21 | | join (SELECT the_geom, name_engl FROM european_wise2008.rbd_f1v3 As f where name_engl like 'Meuse' ) as sub |
22 | | ON ST_Intersects(sub.the_geom,r.the_geom) limit 10; |
23 | | |
24 | | |
25 | | |
26 | | select * from ccm21.riversegments where wso1_id in |
27 | | (291110, |
28 | | 291112, |
29 | | 291115, |
30 | | 291130, |
31 | | 291133, |
32 | | 292843, |
33 | | 292848, |
34 | | 298657, |
35 | | 324473, |
36 | | 324479) |
37 | | --324473 |
38 | | --324479 |
39 | | --291130 --Meuse |
40 | | select * from ccm21.riversegments where wso_id in |
41 | | (291110,--rhin |
42 | | 291112,-- Rhone |
43 | | 291115,-- Seine |
44 | | 291130,-- Meuse |
45 | | 291133,--Scheldt |
46 | | 292843,-- Meuse du bas |
47 | | 292848,-- micro pipi près de l'esuaire |
48 | | 298657)-- petit pipi près de l'esuaire |
49 | | |
50 | | |
51 | | -- i'll insert this in europe.wso |
52 | | |
53 | | insert into europe.wso(wso_id,area) VALUES (291130,'Meuse');-- Meuse |
54 | | insert into europe.wso(wso_id,area) VALUES (292848,'Meuse');-- micro pipi près de l'esuaire |
55 | | insert into europe.wso(wso_id,area) VALUES (298657,'Meuse');-- petit pipi près de l'esuaire |
56 | | insert into europe.wso(wso_id,area) VALUES (292843,'Meuse');--- Meuse du bas |
57 | | |
58 | | |
59 | | drop table if exists clc.meusebelge; |
60 | | create table clc.meusebelge as ( |
61 | | select distinct on (gid) cl.* from clc.clc00_v2_europe cl |
62 | | join |
63 | | (select the_geom from ccm21.catchments c join europe.wso e on e.wso_id=c.wso_id where e.area='Meuse' |
64 | | except |
65 | | select the_geom from ccm21.catchments c join europe.wso1 e on e.wso1_id=c.wso1_id where e.area='France') as ccm |
66 | | on st_intersects(cl.the_geom,ccm.the_geom) |
67 | | |
68 | | ); |
69 | | alter table clc.meusebelge ADD CONSTRAINT meusebelge_pkey PRIMARY KEY (gid); |
70 | | alter table clc.meusebelge ADD CONSTRAINT enforce_dims_the_geom CHECK (ndims(the_geom) = 2); |
71 | | alter table clc.meusebelge ADD CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'MULTIPOLYGON'::text OR the_geom IS NULL); |
72 | | alter table clc.meusebelge ADD CONSTRAINT enforce_srid_the_geom CHECK (srid(the_geom) = 3035); |
73 | | |
74 | | INSERT INTO geometry_columns(f_table_catalog, f_table_schema, f_table_name, f_geometry_column, coord_dimension, srid, "type") |
75 | | SELECT '', 'clc', 'meusebelge', 'the_geom', ST_CoordDim(the_geom), ST_SRID(the_geom), GeometryType(the_geom) |
76 | | FROM clc.meusebelge LIMIT 1; |
77 | | }}} |