create_osm_roads

scripts_v2/create_osm_roads.sh

@@ -0,0 +1,140 @@
+#! /bin/bash
+
+# exit dès que qqch se passe mal
+set -e
+# sortir si "unbound variable"
+#set -u
+
+
+if [ -z "$1" ]
+  then
+    echo "Pas de millésime en argument --> stop"
+    exit 1
+fi
+
+# lecture du fichier de configuration
+. config.sh
+
+
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo "  Création de la couche osm_roads"
+echo ""
+echo ""
+
+
+echo ""
+echo "import phase_1_trace dans la base OSM"
+echo ""
+
+# 1. export du tracé phase 1 depuis la base redadeg
+pg_dump --dbname=postgresql://$DB_USER:$DB_PASSWD@$HOST_DB_redadeg/$DB_NAME \
+    --format=p --no-owner --section=pre-data --section=data --no-privileges --no-tablespaces --no-unlogged-table-data --no-comments \
+    --table phase_1_trace $DB_REDADEG \
+    --file $rep_data/redadeg_trace.sql
+
+
+# 2. import dans la base OSM
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName -c "DROP TABLE IF EXISTS phase_1_trace_$millesime ;"
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName < $rep_data/redadeg_trace.sql
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName -c "ALTER TABLE phase_1_trace RENAME TO phase_1_trace_$millesime ;"
+
+echo ""
+echo "fait"
+echo ""
+
+
+# 3. calcul de la couche osm_roads = intersection buffer trace et routes OSM
+
+echo ">> calcul de la couche osm_roads"
+echo ""
+
+# on supprime puis on recrée la table
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName -c "DROP TABLE IF EXISTS osm_roads_$millesime ;"
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName -c "
+CREATE TABLE osm_roads_$millesime
+(
+  uid bigint NOT NULL,
+  osm_id bigint,
+  highway text,
+  type text,
+  oneway text,
+  ref text,
+  name_fr text,
+  name_br text,
+  the_geom geometry,
+  CONSTRAINT osm_roads_pkey_$millesime PRIMARY KEY (uid),
+  CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'LINESTRING'::text OR geometrytype(the_geom) = 'MULTILINESTRING'::text),
+  CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 2154)
+);"
+
+echo ""
+echo "  table osm_roads_$millesime créée"
+echo ""
+echo "  chargement des données"
+echo ""
+
+PGPASSWORD=$osmDBPassword $PSQL -h $osmDBHost -p $osmDBPort -U $osmDBUser -d $osmDBName -c "WITH trace_buffer AS (
+  SELECT
+    secteur_id,
+    ST_Union(ST_Buffer(the_geom, 25, 'quad_segs=2')) AS the_geom
+  FROM phase_1_trace_$millesime
+  GROUP BY secteur_id
+  ORDER BY secteur_id
+)
+INSERT INTO osm_roads_$millesime
+(
+  SELECT
+    row_number() over() as id,
+    osm_id,
+    highway,
+    CASE 
+        WHEN highway IN ('motorway', 'trunk') THEN 'motorway' 
+        WHEN highway IN ('primary', 'secondary') THEN 'mainroad' 
+        WHEN highway IN ('motorway_link', 'trunk_link', 'primary_link', 'secondary_link', 'tertiary', 'tertiary_link', 'residential', 'unclassified', 'road', 'living_street') THEN 'minorroad' 
+        WHEN highway IN ('service', 'track') THEN 'service' 
+        WHEN highway IN ('path', 'cycleway', 'footway', 'pedestrian', 'steps', 'bridleway') THEN 'noauto' 
+        ELSE 'other' 
+    END AS type,
+    oneway,
+    ref,
+    name AS name_fr,
+    COALESCE(tags -> 'name:br'::text) as name_br,
+    ST_Intersection(ST_Transform(o.way,2154), t.the_geom) AS the_geom
+  FROM planet_osm_line o, trace_buffer t
+  WHERE highway IS NOT NULL AND ST_INTERSECTS(t.the_geom, ST_Transform(o.way,2154))
+) ;"
+
+echo ""
+echo "fait"
+echo ""
+
+
+
+# 4. export de osm_roads depuis la base OSM
+
+echo "transfert de osm_roads_$millesime depuis la base OSM vers la base redadeg"
+echo ""
+
+pg_dump --dbname=postgresql://$osmDBUser:$osmDBPassword@$osmDBHost/$osmDBName \
+    --format=p --no-owner --section=pre-data --section=data --no-privileges --no-tablespaces --no-unlogged-table-data --no-comments \
+    --table osm_roads_$millesime $DB_OSM \
+    --file $rep_data/osm_roads.sql
+
+
+# 5. import dans la base redadeg
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME -c "DROP TABLE IF EXISTS osm_roads;"
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME < $rep_data/osm_roads.sql
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME -c "ALTER TABLE osm_roads_$millesime RENAME TO osm_roads ;"
+
+echo ""
+echo "fait"
+echo ""
+
+echo "fini"
+
+echo ""
+echo ""
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo "  F I N "
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo ""

scripts_v2/sql/create_osm_roads.sql

@@ -0,0 +1,101 @@
+/*
+==========================================================================
+
+    OpenStreetMap : création de la couche du réseau routier OSM
+    pour le calcul d'itinéraires
+
+==========================================================================
+*/
+
+--  A faire dans une base OSM contenant un import osm2pgsql
+
+--  Le but est de créer une table osm_roads qui contiendra
+--  le réseau routier dans une zone-tampon de 500 autour du tracé.
+
+-- on fait ça dans une base à part à cause de la volumétrie des données OSM
+
+
+
+-- 1. import du tracé dans la base OSM
+
+-- export du tracé depuis la base redadeg
+-- pg_dump --file redadeg_trace.sql --host localhost --username redadeg --no-password --format=p --no-owner --section=pre-data --section=data --no-privileges --no-tablespaces --no-unlogged-table-data --no-comments --table public.phase_1_trace redadeg
+
+-- import dans la base osm
+-- psql -U osmbr -d osm -c "DROP TABLE public.phase_1_trace;"
+-- psql -U osmbr -d osm < redadeg_trace.sql
+
+
+
+
+-- 2. création de la table qui va accueillir les tronçons de routes
+
+-- la table qui contient le graphe routier de OSM
+DROP TABLE IF EXISTS osm_roads ;
+CREATE TABLE osm_roads
+(
+  uid bigint,
+  osm_id bigint,
+  highway text,
+  type text,
+  oneway text,
+  ref text,
+  name_fr text,
+  name_br text,
+  the_geom geometry,
+  CONSTRAINT osm_roads_pkey PRIMARY KEY (uid),
+  CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'LINESTRING'::text OR geometrytype(the_geom) = 'MULTILINESTRING'::text),
+  CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 2154)
+);
+
+
+
+-- 3. remplissage de la table à l'aide d'un requête
+-- qui va faire un buffer de 25 m autour du tracé
+-- environ 5 min de traitement
+
+TRUNCATE TABLE osm_roads ;
+WITH trace_buffer AS (
+  SELECT
+    secteur_id,
+    ST_Union(ST_Buffer(the_geom, 25, 'quad_segs=2')) AS the_geom
+  FROM phase_1_trace
+  GROUP BY secteur_id
+  ORDER BY secteur_id
+)
+INSERT INTO osm_roads
+(
+  SELECT
+    row_number() over() as id,
+    osm_id,
+    highway,
+    CASE 
+        WHEN highway IN ('motorway', 'trunk') THEN 'motorway' 
+        WHEN highway IN ('primary', 'secondary') THEN 'mainroad' 
+        WHEN highway IN ('motorway_link', 'trunk_link', 'primary_link', 'secondary_link', 'tertiary', 'tertiary_link', 'residential', 'unclassified', 'road', 'living_street') THEN 'minorroad' 
+        WHEN highway IN ('service', 'track') THEN 'service' 
+        WHEN highway IN ('path', 'cycleway', 'footway', 'pedestrian', 'steps', 'bridleway') THEN 'noauto' 
+        ELSE 'other' 
+    END AS type,
+    oneway,
+    ref,
+    name AS name_fr,
+    COALESCE(tags -> 'name:br'::text) as name_br,
+    ST_Intersection(ST_Transform(o.way,2154), t.the_geom) AS the_geom
+  FROM planet_osm_line o, trace_buffer t
+  WHERE highway IS NOT NULL AND ST_INTERSECTS(t.the_geom, ST_Transform(o.way,2154))
+) ;
+
+
+-- 4. on l'exporte pour la recharger dans la base redadeg
+
+-- commande d'export
+-- pg_dump --file osm_roads.sql --host localhost --username osmbr --no-password --format=p --no-owner --section=pre-data --section=data --no-privileges --no-tablespaces --no-unlogged-table-data --no-comments --table public.osm_roads osm
+
+-- commande d'import
+-- psql -U redadeg -d redadeg -c "DROP TABLE public.osm_roads;"
+-- psql -U redadeg -d redadeg < osm_roads.sql
+
+
+
+