| 1210 | |
| 1211 | {{{#!sql |
| 1212 | |
| 1213 | -- On refais la table observation_place_500 du coup avec la nouvelle version BD_MAP ( et IDEM pour la table station_500 avec le même script) |
| 1214 | |
| 1215 | --select * from dbeel.observation_places where op_gis_layername='BDMAP' -- 13059 |
| 1216 | --select * from dbeel.observation_places where op_gis_layername='BDMAP' and the_geom is NULL -- 113 (c'est toujours les mêmes répétées par ans. |
| 1217 | select st_srid(the_geom) from dbeel.observation_places group by st_srid(the_geom) -- y a 3035 et lambert 93 |
| 1218 | update dbeel.observation_places set the_geom=st_transform(the_geom,3035) |
| 1219 | |
| 1220 | DROP TABLE IF EXISTS dbeel.observation_places_ccm_500; |
| 1221 | CREATE TABLE dbeel.observation_places_ccm_500 as ( |
| 1222 | SELECT distinct on (op_id) op_id, gid, wso1_id, min(distance) as distance,op_gis_layername, the_geom FROM ( |
| 1223 | SELECT op_id, gid , wso1_id,op_gis_layername, CAST(distance(r.the_geom, s.the_geom) as decimal(15,1)) as distance,s.the_geom |
| 1224 | FROM dbeel.observation_places As s |
| 1225 | INNER JOIN ccm21.riversegments r ON ST_DWithin(r.the_geom, s.the_geom,500) --69931 observations sur 87204 lignes |
| 1226 | WHERE s.the_geom IS NOT NULL |
| 1227 | ) AS sub |
| 1228 | GROUP BY op_id,distance, gid, wso1_id, the_geom,op_gis_layername |
| 1229 | ); |
| 1230 | |
| 1231 | alter table dbeel.observation_places_ccm_500 add column id serial; |
| 1232 | |
| 1233 | -- mise à jour de la table geometry_columns |
| 1234 | -- INSERT INTO geometry_columns(f_table_catalog, f_table_schema, f_table_name, f_geometry_column, coord_dimension, srid, "type") |
| 1235 | -- SELECT '', 'dbeel', 'observation_places_ccm_500', 'the_geom', ST_CoordDim(the_geom), ST_SRID(the_geom), GeometryType(the_geom) |
| 1236 | -- FROM dbeel.observation_places_ccm_500 LIMIT 1; |
| 1237 | -- creation d'index, clé primaire, et constraintes qui vont bien |
| 1238 | alter table dbeel.observation_places_ccm_500 add CONSTRAINT enforce_dims_the_geom CHECK (ndims(the_geom) = 2); |
| 1239 | alter table dbeel.observation_places_ccm_500 add CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'POINT'::text OR the_geom IS NULL); |
| 1240 | alter table dbeel.observation_places_ccm_500 add CONSTRAINT enforce_srid_the_geom CHECK (srid(the_geom) = 3035); |
| 1241 | alter table dbeel.observation_places_ccm_500 ADD CONSTRAINT pk_id PRIMARY KEY(id); |
| 1242 | CREATE INDEX indexstation_ccm_500 ON dbeel.observation_places_ccm_500 |
| 1243 | USING GIST ( the_geom GIST_GEOMETRY_OPS ); |
| 1244 | |
| 1245 | }}} |