Changes between Version 13 and Version 14 of CookBook EpA

Timestamp:

Mar 30, 2010 2:24:02 PM (15 years ago)

Author:

cedric

Comment:

--

CookBook EpA

@@ -47 +47 @@
 = Projecting the files on the ccm =
 
+ * first create two index to speed up queries
+
+
+{{{
+CREATE INDEX indexiabret ON ia.iabret
+  USING GIST ( the_geom GIST_GEOMETRY_OPS );
+CREATE INDEX indexiaiav ON ia.iaiav
+  USING GIST ( the_geom GIST_GEOMETRY_OPS );
+}}}
+Joining the riversegment and stationsp2, the_geom is that of stationsp2
+
+-- creation de la table correspondance, je rajoute un the_geom a la requete pour pouvoir voir la table dans Qgis  
+DROP TABLE IF EXISTS bd_map.correspondance;
+CREATE TABLE bd_map.correspondance as (
+        SELECT distinct on (st_codecsp) st_codecsp, gid, min(distance) as distance, the_geom FROM (
+               SELECT st_codecsp, gid ,CAST(distance(r.the_geom, s.the_geom) as  decimal(15,1)) as distance,s.the_geom 
+               FROM bd_map.stationsp2 As s
+               INNER JOIN  riversegments r ON ST_DWithin(r.the_geom, s.the_geom,300)
+               WHERE s.the_geom IS NOT NULL
+               ORDER BY st_codecsp) AS sub 
+        GROUP BY st_codecsp, gid, the_geom
+);
+alter table bd_map.correspondance add column id serial;
+-- mise à jour de la table geometry_columns
+INSERT INTO geometry_columns(f_table_catalog, f_table_schema, f_table_name, f_geometry_column, coord_dimension, srid, "type")
+SELECT '', 'bd_map', 'correspondance', 'the_geom', ST_CoordDim(the_geom), ST_SRID(the_geom), GeometryType(the_geom)
+FROM bd_map.correspondance LIMIT 1;
+
+-- creation d'index, clé primaire, et constraintes qui vont bien
+alter table bd_map.correspondance add CONSTRAINT enforce_dims_the_geom CHECK (ndims(the_geom) = 2);
+alter table bd_map.correspondance add  CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'POINT'::text OR the_geom IS NULL);
+alter table bd_map.correspondance add  CONSTRAINT enforce_srid_the_geom CHECK (srid(the_geom) = 3035);
+alter table bd_map.correspondance ADD CONSTRAINT pk_id PRIMARY KEY(id);
+CREATE INDEX indexcorrespondance ON bd_map.correspondance
+  USING GIST ( the_geom GIST_GEOMETRY_OPS );
+