dimanche 15 juin 2008
PgRouting 1.02: Test de la fonctionnalité shortest_path_shooting_star()
Par david techer, dimanche 15 juin 2008 à 14:56 :: PostGIS et PostgreSQL
Le descriptif de cette fonctionnalité est fournie à http://pgrouting.postlbs.org/wiki/ShootingStar. Cette fonctionnalité permet notamment de prendre en compte les sens de parcours. On peut si on veut dans la colonne _to_cost préciser si le coût de parcours du tronçon est fonction de la distance, du prix de péage etc...
Ici pour la compilation, on ne se prend pas la tête!
tar xzf pgRouting-1.02.tgz ; cd pgrouting; cmake -DWITH_TSP=ON -DWITH_DD=ON \ -DPOSTGRESQL_INCLUDE_DIR=/usr/local/pgsql/include/server/ -DPOSTGRESQL_LIBRARIES=/usr/local/pgsql/lib/ \ -DCGAL_INCLUDE_DIR=/usr/local/cgal/include/ -DCGAL_LIBRARIES=/usr/local/cgal/lib/ make make installPour la création de la base avec les fonctionnalités adéquates
dropdb routing; createdb -U postgres -E UNICODE routing; createlang -U postgres plpgsql routing; psql -U postgres -f /usr/local/pgsql/share/lwpostgis.sql routing; psql -U postgres -f /usr/local/pgsql/share/spatial_ref_sys.sql routing; psql -U postgres -f /usr/share/postlbs/routing_core.sql routing; psql -U postgres -f /usr/share/postlbs/routing_core_wrappers.sql routing; psql -U postgres -f /usr/share/postlbs/routing_tsp.sql routing; psql -U postgres -f /usr/share/postlbs/routing_tsp_wrappers.sql routing; psql -U postgres -f /usr/share/postlbs/routing_dd.sql routing; psql -U postgres -f /usr/share/postlbs/routing_dd_wrappers.sql routing;On télécharge le fichier le fichier suivant http://www.davidgis.fr/download/reseau.sql
wget http://www.davidgis.fr/download/reseau.sql psql routing -f reseau.sqlLe contenu sera le suivant
routing=# select id,sens,source,target,geometrytype(the_geom),x1,y1,x2,y2,cost,reverse_cost,rule,to_cost from reseau order by id limit 10; id | sens | source | target | geometrytype | x1 | y1 | x2 | y2 | cost | reverse_cost | rule | to_cost ----+-------------+--------+--------+-----------------+----+-----+----+-----+------------------+--------------+------+--------- 1 | double sens | 1 | 2 | MULTILINESTRING | 1 | 0 | 5 | 0 | 4 | 4 | | 2 | double sens | 2 | 3 | MULTILINESTRING | 5 | 0 | 5 | 6 | 6 | 6 | | 3 | double sens | 4 | 5 | MULTILINESTRING | 0 | 7.5 | 3 | 7.5 | 3 | 3 | | 4 | sens direct | 5 | 3 | MULTILINESTRING | 3 | 7.5 | 5 | 6 | 2.91421356237309 | -1 | | 5 | sens direct | 3 | 6 | MULTILINESTRING | 5 | 6 | 7 | 7.5 | 2.91421356237309 | -1 | | 6 | sens direct | 6 | 7 | MULTILINESTRING | 7 | 7.5 | 5 | 9 | 2.91421356237309 | -1 | | 7 | sens direct | 7 | 5 | MULTILINESTRING | 5 | 9 | 3 | 7.5 | 2.91421356237309 | -1 | | 8 | double sens | 6 | 8 | MULTILINESTRING | 7 | 7.5 | 11 | 7.5 | 4 | 4 | | 9 | double sens | 8 | 9 | MULTILINESTRING | 11 | 7.5 | 11 | 11 | 3.5 | 3.5 | | 10 | double sens | 8 | 10 | MULTILINESTRING | 11 | 7.5 | 14 | 7.5 | 3 | 3 | | (10 rows)Quand
sens = 'sens direct', celà signifie qu'on traverse un tronçon relatif à un rond-point
J'aurais donc l'image suivante
Fig 1. Le réseau
shooting.sh dont le contenu est
#!/bin/sh
source=$1
target=$2
echo "begin transaction;
drop table if exists aller,retour;
create table aller(gid int4) with oids;
SELECT AddGeometryColumn( 'aller', 'the_geom', -1, 'MULTILINESTRING', 2 );
insert into aller(the_geom) select the_geom from reseau where id in
(SELECT edge_id FROM shortest_path_shooting_star('SELECT id, source, target, cost,reverse_cost, x1, y1, x2, y2, rule, to_cost FROM reseau',
(select id from reseau where source=$source or target=$source limit 1),(select id from reseau where target=$target or source=$target limit 1), true, true));
create table retour(gid int4) with oids;
SELECT AddGeometryColumn( 'retour', 'the_geom', -1, 'MULTILINESTRING', 2 );
insert into retour(the_geom) select the_geom from reseau where id in
(SELECT edge_id FROM shortest_path_shooting_star('SELECT id, source, target, cost,reverse_cost, x1, y1, x2, y2, rule, to_cost FROM reseau',
(select id from reseau where target=$target or source=$target limit 1),(select id from reseau where source=$source or target=$source limit 1), true, true));
end transaction;"|psql routing
On le rend exécutable en accordant les droits nécessaires
chmod u+x shooting.sh
remarque: j'ai fait ce script à l'arrache pour tester rapidement mais on peut l'améliorer
Par exemple pour aller de 37 Ã 45./shooting 37 45Dans Qgis, j'ai alors l'image suivante
Fig 2 - Sens aller de 37 Ã 45
Fig 3 - Sens retour de 45 Ã 37
Fig 4 - Sens aller/retour pour (36,26)
sur la route, j'étais en nage. Cette semaine, j'ai pu faire une mise au point pour la suite sur Antibes et lundi et mardi je profite de ce congé exceptionel (merci boss 
) pour deux rendez-vous 
. Celà fait du monde.... Mais bon rançon de la gloire oblige 
