diff --git a/scripts/create_tables.sql b/scripts/create_tables.sql index a532b44..9d05330 100644 --- a/scripts/create_tables.sql +++ b/scripts/create_tables.sql @@ -240,8 +240,8 @@ SELECT topology.CreateTopology('osm_roads_topo', 2154); DROP TABLE IF EXISTS phase_2_trace_pgr ; CREATE TABLE phase_2_trace_pgr ( + secteur_id integer, -- info de routage - id bigint, path_seq bigint, node bigint, cost double precision, @@ -255,7 +255,7 @@ CREATE TABLE phase_2_trace_pgr name_fr text, name_br text, the_geom geometry, - CONSTRAINT phase_2_trace_pkey PRIMARY KEY (id), + CONSTRAINT phase_2_trace_pkey PRIMARY KEY (secteur_id, path_seq), 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) ); diff --git a/scripts/traitements_phase_2.sh b/scripts/traitements_phase_2.sh index a00b5f8..2943493 100755 --- a/scripts/traitements_phase_2.sh +++ b/scripts/traitements_phase_2.sh @@ -1,4 +1,13 @@ -#!/bin/sh +#!/bin/bash + +set -e +set -u + +PSQL=/usr/bin/psql +DB_HOST=localhost +DB_NAME=redadeg +DB_USER=redadeg + cd /data/www/vhosts/ar-redadeg_openstreetmap_bzh/htdocs/scripts/ @@ -20,13 +29,93 @@ curl -sS http://umap.openstreetmap.fr/fr/datalayer/817222/ > data/phase_2_umap_ # note : les coordonnées sont en 3857 mais la déclaration de la table = 4326 -psql -U redadeg -d redadeg -c "DROP TABLE phase_2_pk_secteur_3857 CASCADE;" +$PSQL -U $DB_USER -d $DB_NAME -c "DROP TABLE phase_2_pk_secteur_3857 CASCADE;" ogr2ogr -f "PostgreSQL" PG:"host=localhost user=redadeg password=redadeg dbname=redadeg" data/phase_2_umap_pk_secteur.geojson -nln phase_2_pk_secteur_3857 -lco GEOMETRY_NAME=the_geom -explodecollections -overwrite - -# on crée les tables en 3948 et on fait bien d'autres choses : -psql -U redadeg -d redadeg < traitements_phase_2.1.sql +# on crée les tables en 3948 +# et bien d'autres choses : +# - recalage des PK secteurs sur un nœud du réseau routable +$PSQL -U $DB_USER -d $DB_NAME < traitements_phase_2.1.sql + +# ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +# ici on va calculer un itinéraire pour chaque secteur +# en utilisant les PK de début (ou fin) de chaque secteur + +# https://www.manniwood.com/postgresql_and_bash_stuff/index.html + + +# on commence par vider la table qui contiendra les calculs d'itinéraires +$PSQL -h $DB_HOST -U $DB_USER -c "TRUNCATE TABLE phase_2_trace_pgr ;" + + +# on fait la requête qui va donner une liste de PK de secteurs +# et on calcule un itinéraire entre le PK de début et le PK suivant + +$PSQL -X -h $DB_HOST -U $DB_USER $DB_NAME \ + -c "SELECT s.id, replace(s.nom_fr,' ','') AS nom_fr, replace(s.nom_br,' ','') AS nom_br, pk.pgr_node_id +FROM phase_2_pk_secteur pk JOIN secteur s ON pk.id = s.id +ORDER BY pk.id ;" \ + --single-transaction \ + --set AUTOCOMMIT=off \ + --set ON_ERROR_STOP=on \ + --no-align \ + -t \ + --field-separator ' ' \ + --quiet | while read -a Record ; do + + # ici commence la boucle sur les PK de secteurs + echo "----------------------------" + + #IFS="|" pour forcer un délimiteur mais ne fonctionne pas : les espaces sont compris comme des séparateurs + # alors la requête supprime les espaces. TODO + + # le premier PK = PK de début + secteur_id=${Record[0]} + secteur_nom_fr="${Record[1]}" + secteur_nom_br="${Record[2]}" + pk_id_start=${Record[3]} + + # maintenant il faut une 2e requête pour aller trouver le PK de fin + # ce PK = le PK de début du secteur suivant + read pk_id_end <<< $($PSQL -h $DB_HOST -U $DB_USER --no-align -t --quiet \ + -c "SELECT pgr_node_id FROM phase_2_pk_secteur WHERE id = $secteur_id + 1 ;") + + # on teste si on récupère qqch sinon ça veurt dire qu'on a pas de nœud de fin donc impossible de calculer un itinéraire + if [[ -n "$pk_id_end" ]]; + then + echo "calcul d'un itinéraire pour le secteur $secteur_nom_fr ($pk_id_start --> $pk_id_end)" + + $PSQL -h $DB_HOST -U $DB_USER -c \ + "INSERT INTO phase_2_trace_pgr + SELECT + $secteur_id AS secteur_id, + -- info de routage + a.path_seq, + a.node, + a.cost, + a.agg_cost, + -- infos OSM + b.osm_id, + b.highway, + b.\"type\", + b.oneway, + b.ref, + b.name_fr, + b.name_br, + b.the_geom + FROM pgr_dijkstra( + 'SELECT id, source, target, cost, reverse_cost FROM osm_roads_pgr', $pk_id_start, $pk_id_end) as a + JOIN osm_roads_pgr b ON a.edge = b.id ;" + + echo "fait" + fi + + # fin de la boucle + echo "" + +done +