phase 2 : calcul routage + patch

data/2022/ar-redadeg-2022.xml

@@ -54,6 +54,17 @@
                 queryable="true"
             >
             </layer>
+            <layer
+                id="osm_roads_patch"
+                name="Zone de patch du filaire de voies"
+                visible="false"
+                url="https://ar-redadeg.openstreetmap.bzh/geoserver/redadeg_2022/wms"
+                infoformat="application/vnd.ogc.gml"
+                tiled="false"
+                searchable="false"
+                queryable="true"
+            >
+            </layer>
             <layer
                 id="phase_2_trace_pgr"
                 name="Tracé recalé sur le réseau viaire OSM"

scripts_v2/phase_2_routing.sh

@@ -0,0 +1,177 @@
+#! /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
+
+
+
+# ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+# 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
+
+
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo "  Calcul des itinéraires (pgrouting)"
+echo ""
+
+
+echo "patch de la couche osm_roads_pgr pour les cas particuliers"
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME < sql/patch_osm_roads_pgr.sql
+echo "  fait"
+echo ""
+
+# on commence par vider la table qui contiendra les calculs d'itinéraires
+echo "vidage de la couche de routage"
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME -c "TRUNCATE TABLE phase_2_trace_pgr ;"
+echo "  fait"
+echo ""
+
+# ensuite : on supprime les tronçons ciblés par la couche de points de nettoyage
+# AVANT de calculer les itinéraires
+echo "nettoyage de la couche de routage par les points ciblés"
+PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME -c "UPDATE osm_roads_pgr SET cost = 1000000, reverse_cost = 1000000 WHERE id IN (SELECT r.id FROM osm_roads_pgr r JOIN phase_2_point_nettoyage p ON r.id = p.edge_id);"
+echo "  fait"
+echo ""
+
+
+# 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
+
+# on va utiliser un compteur pour pouvoir sauter un sous-secteur à un autre
+counter=1
+# autre variables de contrôle
+longueur_totale=0
+longueur_inseree=0
+
+
+PGPASSWORD=$DB_PASSWD $PSQL -X -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME \
+    -c "SELECT pk.id, s.id AS secteur_id, pk.pgr_node_id, replace(s.nom_fr,' ','') AS nom_fr, replace(s.nom_br,' ','') AS nom_br, replace(pk.name,' ','_') AS name 
+FROM phase_2_pk_secteur pk JOIN secteur s ON pk.secteur_id = s.id
+ORDER BY s.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
+    pk_id=${Record[0]}
+    secteur_id=${Record[1]}
+    pk_id_start=${Record[2]}
+    secteur_nom_fr="${Record[3]}"
+    secteur_nom_br="${Record[4]}"
+    troncon_name=${Record[5]}
+
+    echo "  $secteur_id | $secteur_nom_fr / $secteur_nom_br"
+    echo "  tronçon : $troncon_name"
+    echo "  PK ID = $pk_id"
+    echo "  start node = $pk_id_start"
+
+    # on fait une requête pour récupérer l'id du nœud de routage de fin
+    # ce nœud = le PK de début du secteur suivant
+    read pk_id_end <<< $(PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME --no-align -t --quiet \
+      -c "SELECT pgr_node_id FROM phase_2_pk_secteur ORDER BY id OFFSET $counter LIMIT 1 ;" )
+
+    echo "  end node = $pk_id_end"
+
+
+    # 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 de l'itinéraire"
+
+        PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME -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,
+          CASE
+          WHEN b.name_fr IS NULL AND b.ref IS NOT NULL THEN b.ref
+        ELSE b.name_fr
+          END AS name_fr,
+          CASE
+          WHEN b.name_br IS NULL AND b.name_fr IS NULL AND b.ref IS NOT NULL THEN b.ref
+        WHEN b.name_br IS NULL AND b.name_fr IS NOT NULL THEN '# da dreiñ e brezhoneg #'
+        ELSE b.name_br
+          END AS 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 ;" >> /dev/null
+
+        # on fait une requête pour voir la longueur insérée
+        # en fait : la longueur totale - la longueur totale lors du précédent calcul
+        read longueur_base <<< $(PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME --no-align -t --quiet \
+          -c "SELECT trunc(SUM(ST_Length(the_geom))/1000) as longueur_totale FROM phase_2_trace_pgr ;" )
+        longueur_inseree=$(($longueur_base-$longueur_totale))
+        longueur_totale=$longueur_base
+        
+        # une alerte si 0 km insérés
+        if [ $longueur_inseree -eq 0 ] ;
+        then
+          echo ""
+          echo "    E R R E U R   !!!!!!!!"
+          echo ""
+        else
+          echo "  fait : $longueur_inseree km (total = $longueur_totale km)"
+        fi
+
+    else
+        echo ""
+        echo "  E R R E U R   !!!!!!!!"
+        echo "  impossible de calculer un itinéraire pour ce secteur"
+        echo ""
+    fi
+
+
+    # fin de la boucle
+    # on incrémente le compteur
+    ((counter++))
+    echo ""
+    
+done
+
+echo "  Calcul des itinéraires terminé"
+echo ""
+
+
+
+
+echo ""
+echo ""
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo "  F I N traitements phase 2 : calcul des itinéraires"
+echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
+echo ""

scripts_v2/sql/patch_osm_roads_pgr.sql

@@ -0,0 +1,40 @@
+
+
+
+-- suppression des objets couche osm_roads_pgr qui intersectent avec les zones de boucles
+DELETE FROM osm_roads_pgr WHERE id IN
+(
+  SELECT a.id 
+  FROM osm_roads_pgr a, osm_roads_pgr_patch_mask m
+  WHERE ST_INTERSECTS(a.the_geom, m.the_geom)
+);
+
+-- collage des objets de la couche osm_roads_pgr_patch à la place des objets supprimés
+-- on crée un identifiant négatif
+INSERT INTO osm_roads_pgr
+  SELECT
+    0-a.id AS id,
+    a.osm_id, a.highway, a.type, a.oneway, a.ref, a.name_fr, a.name_br,
+    NULL, NULL, NULL, NULL,
+    a.the_geom
+  FROM osm_roads_pgr_patch a, osm_roads_pgr_patch_mask m
+  WHERE ST_INTERSECTS(a.the_geom, m.the_geom);
+
+
+-- calcul des 2 attributs de coût (= longueur)
+UPDATE osm_roads_pgr 
+SET cost = st_length(the_geom), reverse_cost = st_length(the_geom)
+WHERE id < 0 ;
+
+
+-- recrée des nœuds uniquement sur les zones de patch
+SELECT pgr_nodeNetwork('osm_roads_pgr', 0.001, rows_where := 'id < 0');
+
+-- recalcul la topologie pgRouting uniquement sur les zones de patch
+SELECT pgr_createTopology('osm_roads_pgr', 0.001, rows_where := 'id < 0', clean := false);
+
+
+-- recalcul toute la topologie pgRouting en renumérotant de zéro
+--SELECT pgr_createTopology('osm_roads_pgr', 0.001, clean := true);
+
+--SELECT pgr_analyzegraph('osm_roads_pgr', 0.001);

scripts_v2/update_osm_roads_pgr.sh

@@ -71,9 +71,6 @@ echo ""
 echo ">> maj de la couche osm_roads_pgr qui sert au routage depuis la topologie"
 PGPASSWORD=$DB_PASSWD $PSQL -h $DB_HOST -p $DB_PORT -U $DB_USER -d $DB_NAME < sql/update_osm_roads_pgr.sql
 
-#echo ">> patch de la couche osm_roads_pgr pour les cas particuliers"
-#$PSQL -h $DB_HOST -U $DB_USER -d $DB_NAME < patch_osm_roads_pgr.sql
-
 
 echo ""
 echo ""
@@ -81,7 +78,5 @@ echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
 echo "  F I N "
 echo "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
 echo ""
-echo "/!\ patcher le filaire de voie si nécessaire"
-echo ""
 echo ""