Le blog de Jean David TECHER, un Réunionnais à Saint-Priest/Lyon

Aller au contenu | Aller au menu | Aller à la recherche


1 2 3 4 5 6 7 8 9 10 >

samedi 31 mars 2007

PgRouting - autres tests avec plus de rond-points - fonction A* - 2/3

Par david techer, samedi 31 mars 2007 à 19:56 :: PostGIS et PostgreSQL

Bon je me suis amusé à  reprendre le réseau que j'avais pris dans mon précédent billet - billet concernant l'installation de PgRouting sous Ubuntu - pour lui ajouter de nouveaux tronçons en utilisant OpenJump. Pour se faire, j'ai importé ma table postgis de mon graphe précédent. En utilisant depuis le menu de OpenJump - dans les menu - "Ajout de nouvelles", j'ai complété le dessin. J'ai ensuite tout exporté en shapefile que j'ai nommé troncon_route.shp - toujours grà¢ce à  OpenJump -. Il est cool cet outil .

Pour importer le shapefile, comme d'habitude rien de compiquer: 
shp2pgsql -dDI troncon_route.shp troncon_route| psql -h 192.168.0.5 -U postgres geocoding
N.B: Le shapefile en question peut être téléchargé à  http://www.davidgis.fr/download/troncon_route.zip.> Il ne contient rien d'extraordinaire pour le moment:
SELECT gid,sens,astext(the_geom) FROM  troncon_route ORDER BY  gid "
 gid |    sens     |                astext
-----+-------------+--------------------------------------
   1 | double sens | MULTILINESTRING((1 0,5 0))
   2 | double sens | MULTILINESTRING((5 0,5 6))
   3 | double sens | MULTILINESTRING((0 7.5,3 7.5))
   4 | sens direct | MULTILINESTRING((3 7.5,3 7,4 6,5 6))
   5 | sens direct | MULTILINESTRING((5 6,6 6,7 7,7 7.5))
   6 | sens direct | MULTILINESTRING((7 7.5,7 8,6 9,5 9))
   7 | sens direct | MULTILINESTRING((5 9,4 9,3 8,3 7.5))
   8 | double sens | MULTILINESTRING((7 7.5,11 7.5))
   9 | double sens | MULTILINESTRING((11 7.5,11 11))
  10 | double sens | MULTILINESTRING((11 7.5,14 7.5))
  11 | double sens | MULTILINESTRING((14 7.5,21 7.5))
  ... | double sens | MULTILINESTRING((..................))
Bon je suis pas un grand spéclialiste du dessin - ça on s'en doute - mais au final celà  ressemble à  ça:

Fig 1. Mon graphe simulant un réseau routier avec rond-points.
Pour mon réseau, j'ai les sources et les targets suivants

Fig 2. Noeuds de mon graphe qui serviront de sources et de targets .
Ma table troncon_route_edges, je l'ai créé en faisant, les requêtes SQL suivantes. Par rapport à  mon précédent billet, j'ai quand même réussi à  améliorer les requêtes. 
BEGIN TRANSACTION;
--SELECT dropgeometrytable('troncon_route');
SELECT drop_graph_tables('troncon_route');

/*
   Ajouter les colonnes adéquates
*/
ALTER TABLE troncon_route ADD column source_id int4;
ALTER TABLE troncon_route ADD column target_id int4;
ALTER TABLE troncon_route ADD column edge_id int4;
/*
   Mettre à  jour le srid=1 sinon pgdijkstra gueule 8-(
*/
SELECT UPDATEgeometrysrid('troncon_route','the_geom',-1);

SELECT assign_vertex_id('troncon_route',0.00001);
/*
   Ok...Je crée mon graphe
*/
SELECT create_graph_tables('troncon_route', 'int4');

--SELECT UPDATE_cost_FROM_distance('troncon_route');
ALTER TABLE troncon_route_edges ADD column sens text;
ALTER TABLE troncon_route_edges ADD column x1 double precision;
ALTER TABLE troncon_route_edges ADD column y1 double precision;
ALTER TABLE troncon_route_edges ADD column x2 double precision;
ALTER TABLE troncon_route_edges ADD column y2 double precision;
ALTER TABLE troncon_route_edges ADD column edge_id int4;
/*
   Mise à  jour des colonnes x1,y1,x2,y2 originaux par rapport aux données géométriques de la table troncon_route
   et mise à  jour des colonnes sens et edge_id
*/
UPDATE troncon_route_edges SET cost=(select length(the_geom) FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
UPDATE troncon_route_edges SET x1=(select x(startpoint(the_geom)) FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
UPDATE troncon_route_edges SET y1=(select y(startpoint(the_geom)) FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
UPDATE troncon_route_edges SET x2=(select x(endpoint(the_geom)) FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
UPDATE troncon_route_edges SET y2=(select y(endpoint(the_geom)) FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
UPDATE troncon_route_edges SET edge_id=(select edge_id FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.edge_id);
UPDATE troncon_route_edges SET sens=(select sens::text FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.sens);
SELECT AddGeometryColumn( 'troncon_route_edges', 'the_geom', -1, 'MULTILINESTRING', 2 );
UPDATE troncon_route_edges SET the_geom=(select the_geom FROM troncon_route g WHERE g.edge_id=id GROUP BY id,g.the_geom);
/*
    Tout ce qui est à  double sens je le garde
*/
UPDATE troncon_route_edges SET reverse_cost=cost;
/*
    Paramétrer le coà»t des tronçons à  sens unique
*/
UPDATE troncon_route_edges SET reverse_cost=-1 WHERE sens='sens direct';

END TRANSACTION;
VACUUM FULL ANALYZE ;
Maintenant, ma table troncon_route_edges est parfaitement complète car la requête 
SELECT id,sens,astext(the_geom),x1,y1,x2,y2,source,target,edge_id,cost,reverse_cost FROM  troncon_route_edges ORDER BY  id LIMIT 10
me renvoit 
 id |    sens     |                astext                | x1 | y1  | x2 | y2  | source | target | edge_id |       cost       | reverse_cost
----+-------------+--------------------------------------+----+-----+----+-----+--------+--------+---------+------------------+--------------
  1 | double sens | MULTILINESTRING((1 0,5 0))           |  1 |   0 |  5 |   0 |      1 |      2 |       1 |                4 |            4
  2 | double sens | MULTILINESTRING((5 0,5 6))           |  5 |   0 |  5 |   6 |      2 |      3 |       2 |                6 |            6
  3 | double sens | MULTILINESTRING((0 7.5,3 7.5))       |  0 | 7.5 |  3 | 7.5 |      4 |      5 |       3 |                3 |            3
  4 | sens direct | MULTILINESTRING((3 7.5,3 7,4 6,5 6)) |  3 | 7.5 |  5 |   6 |      5 |      3 |       4 | 2.91421356237309 |           -1
  5 | sens direct | MULTILINESTRING((5 6,6 6,7 7,7 7.5)) |  5 |   6 |  7 | 7.5 |      3 |      6 |       5 | 2.91421356237309 |           -1
  6 | sens direct | MULTILINESTRING((7 7.5,7 8,6 9,5 9)) |  7 | 7.5 |  5 |   9 |      6 |      7 |       6 | 2.91421356237309 |           -1
  7 | sens direct | MULTILINESTRING((5 9,4 9,3 8,3 7.5)) |  5 |   9 |  3 | 7.5 |      7 |      5 |       7 | 2.91421356237309 |           -1
  8 | double sens | MULTILINESTRING((7 7.5,11 7.5))      |  7 | 7.5 | 11 | 7.5 |      6 |      8 |       8 |                4 |            4
  9 | double sens | MULTILINESTRING((11 7.5,11 11))      | 11 | 7.5 | 11 |  11 |      8 |      9 |       9 |              3.5 |          3.5
 10 | double sens | MULTILINESTRING((11 7.5,14 7.5))     | 11 | 7.5 | 14 | 7.5 |      8 |     10 |      10 |                3 |            3
(10 lignes)

J'obtiens des résultats très intéressants:


Fig 3. Parcours à  l'aller.

Fig 4. Parcours au retour.
Pour réaliser le test précédent, je ne cache pas que je ne maà®trise pas encore l'ordre pour la fonction shortest_path_asta(SQL,source,target,true/false,true/false) pour avoir aller et le retour. Selon l'ordre celà  me fait donc 8 possibilités: 4 ordres possibles pour l'aller et 4 pour le retour. Pour le retour, il suffit donc d'intervertir la place du source et du target. J'ai créé 8 tables associées chacune à  un des ordre possibles. Par exemple ma première table sera associée à  l'ordre source,target,,false,false, la seconde table sera associée à  l'ordre source,target,false,true et ainsi de suite.... Je me suis fait un petit script shell - nommé test.sh - dont le contenu est le suivant 
echo "begin transaction;
drop table if exists test1,test2,test3,test4,test5,test6,test7,test8;
create table test1(gid int4) with oids;SELECT AddGeometryColumn( 'test1', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test2(gid int4) with oids;SELECT AddGeometryColumn( 'test2', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test3(gid int4) with oids;SELECT AddGeometryColumn( 'test3', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test4(gid int4) with oids;SELECT AddGeometryColumn( 'test4', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test5(gid int4) with oids;SELECT AddGeometryColumn( 'test5', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test6(gid int4) with oids;SELECT AddGeometryColumn( 'test6', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test7(gid int4) with oids;SELECT AddGeometryColumn( 'test7', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table test8(gid int4) with oids;SELECT AddGeometryColumn( 'test8', 'the_geom', -1, 'MULTILINESTRING', 2 );


insert into test1(the_geom) (select the_geom from troncon_route_edges  where edge_id in
 (select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$1,$2,false,false)));

insert into test2(the_geom) (select the_geom from troncon_route_edges  where edge_id in
 (select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$1,$2,false,true)));

insert into test3(the_geom) (select the_geom from troncon_route_edges  where edge_id in 
(select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$1,$2,true,true)));

insert into test4(the_geom) (select the_geom from troncon_route_edges  where edge_id in 
(select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$1,$2,true,false)));



insert into test5(the_geom) (select the_geom from troncon_route_edges  where edge_id in 
(select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$2,$1,false,false)));

insert into test6(the_geom) (select the_geom from troncon_route_edges  where edge_id in 
(select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$2,$1,false,true)));

insert into test7(the_geom) (select the_geom from troncon_route_edges  where edge_id in
 (select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$2,$1,true,true)));

insert into test8(the_geom) (select the_geom from troncon_route_edges  where edge_id in 
(select edge_id from shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',$2,$1,true,false)));

end transaction;"|psql -h 192.168.0.5 -U postgres geocoding
Pour l'exécuter, il me suffit de lui passer les paramètre suivants o๠source et target sont les noeuds de la figure 2. 
./test.sh source target
Pour l'exemple ici entre le noeud 38 et 48, j'ai simplement fait 
./test4.sh 38 48
L'affichage a ensuite eu lieu tout simplement sous QGIS .

un commentaire :: aucun trackback

Installation de mon graveur DVD-RW Philips SDP2411BD - comme quoi fallait lire la notice :)

Par david techer, samedi 31 mars 2007 à 02:01 :: Ubuntu

Purée tout fier d'avoir récemment fait l'acquisition de mon graveur - pour 45 € -, j'ai décidé ce soir de l'installer sur mon ordinateur de bureau. Tout fier de moi, je sors mon tournevis, un chiffon, ma bombe à  gaz dépoussiérant - pour enlever la poussière -, je débranche mon ancien graveur de cd qui a rendu l'à¢me depuis bien longtemps. Je branche le cà¢ble IDE - le fameux cable plat tout grisé -, le cà¢ble d'alimentation. Je rebranche tout. Allé zou on redémarre. .........Il est bien reconnu sous Dapper mais quand j'insère un CD ou un DVD, rien ne se produit !.... Bon pas de panique, direction le chan IRC de Ubuntu au cas o๠il y aurait un topic pour ce genre de graveur ou d'essayer de bidouiller quoique ce soit. Et là  je tombe sur The_Marauder - ze master of ze master - qui veut bien me donner un coup de main! A 00h30, faut le faire quand même! Il est cool ce mec !

Ben on a pas cherché longtemps! L'idiot que je suis aurait dà» lire la notice avant de tout monter comme un bourrin! ==> Les lecteur/graveur sont montés et vendus généralement en étant en mode master alors que sur mon PC, il fallait juste le mettre en slave! La honte! Ben oui étant donné que je remplaçais mon ancien graveur de cd, il faut mettre le nouveau dans la même configuration! CQFD. Comme quoi...no comment

un commentaire :: aucun trackback

jeudi 29 mars 2007

Avoir des noms de fichiers en lettres minuscules

Par david techer, jeudi 29 mars 2007 à 04:56 :: Linux

Supposons qu'on ait le répertoire contenant les fichiers suivants

 ls /home/postgres/RESEAU_ROUTIER/*.TAB
/home/postgres/RESEAU_ROUTIER/CARREFOUR_COMPLEXE.TAB  /home/postgres/RESEAU_ROUTIER/NON_COMMUNICATION.TAB
/home/postgres/RESEAU_ROUTIER/NOEUD_ROUTIER.TAB       /home/postgres/RESEAU_ROUTIER/TRONCON_ROUTE.TAB

Et qu'on veuille ensuite pouvoir obtenir

carrefour_complexe
noeud_routier
non_communication
troncon_route

Pour se faire, il suffit d'appliquer la boucle suivante

for file in $(ls /home/postgres/RESEAU_ROUTIER/*.TAB);do table=`basename ${file} | cut -d '.' -f 1|tr [:upper:] [:lower:]`;echo ${table};done

Je me garde cette astuce au chaud car j'en ai souvent besoin pour importer à  chaud des shapefiles. En effet, les noms des fichiers majuscules convertis en minuscules me servent par la suite de nom de table.

un commentaire :: aucun trackback

mardi 27 mars 2007

Installation et tests de pgRouting 0.99 sous Ubuntu Dapper/Edgy: Recherche du plus cours chemin sur un réseau routier/graphe

Par david techer, mardi 27 mars 2007 à 16:39 :: PostGIS et PostgreSQL

Bon je viens de me lancer dans l'installation de pgRouting sous Dapper pour voir. Pas trop de soucis en particulier, cette version semble s'installer sans problème

PRE-REQUIS

Pour mes tests d'installation, j'ai

OUTILVERSION
POSTGRESQL8.2.3
POSTGIS1.2.1
GEOS3.0.0rc4
PROJ4.5.0

Ensembles des outils déjà  installés sur ma machine.

Rappel: pour installer ces outils, il y a une petite doc pour Ubuntu Edgy à  http://www.davidgis.fr/documentation/win32/html/ch03.html

INSTALLATION

Pour suivre l'installation, rien ne vaut la traduction française faite par Gérald sur le site de postgis.fr. Pour les recommendations originales, on pourra se fier à  http://www.postlbs.org/postlbs-cms/en/node/21

GAUL

Il faut se rendre à  http://prdownloads.sourceforge.net/gaul/gaul-devel-0.1849-0.tar.gz?download pour télécharger gaul-devel-0.1849-0.tar.bz2. Ensuite 

tar xvjf gaul-devel-0.1849-0.tar.bz2
cd gaul-devel-0.1849-0
./configure --enable-slang=no
make
make install
BGL -

Rien de spécial à  faire

apt-get install libboost-graph-dev libboost-graph1.33.1
CGAL
wget ftp://ftp.mpi-sb.mpg.de/pub/outgoing/CGAL/CGAL-3.2.1.tar.gz
tar xvzf CGAL-3.2.1.tar.gz
cd CGAL-3.2.1
./install_cgal --prefix=/usr/local/cgal --with-boost=n --without-autofind -ni /usr/bin/g++
pgRouting
wget http://www.postlbs.org/postlbs-cms/files/downloads/pgRouting-0.9.9.tgz
tar xvzf pgRouting-0.9.9.tgz
cd routing/
export CGAL_MAKEFILE=/usr/local/cgal/make/makefile_i686_Linux-2.6_g++-4.0.3
./configure --with-cgal=/usr/local/cgal --with-gaul=/usr/local/
make
make install

TESTS DE CHARGEMENT DU FICHIER SQL ET SIMULATION SUR UN GRAPHE MINIMALISTE AVEC DEUX ROND-POINTS

Pour charger les fichiers requis par pgRouting - routing.sql et routing_postgis.sql -, il est nécessaire avant tout de faire

echo /usr/local/cgal/lib/i686_Linux-2.6_g++-4.0.3/ >> /etc/ld.so.conf
ldconfig
Celà  est nécessaire afin que routing.so puisse trouver libCGAL.so à  laquelle elle est liée. Pour la suite, on devient super-utilisateur de PostgreSQL, et on charge le tout: 
su postgres
createdb georoute
createlang plpgsql georoute
psql -d georoute -f /usr/local/pgsql/share/lwpostgis.sql 
psql -d georoute -f /usr/local/pgsql/share/spatial_ref_sys.sql
psql -d georoute -f /usr/local/pgsql/share/contrib/routing.sql
psql -d georoute -f /usr/local/pgsql/share/contrib/routing_postgis.sql
SIMULATION

Je vais me créer une table troncon_route qui va contenir aussi bien des routes à  double-sens mais aussi deux rond-points. Mon petit réseau ressemblera donc à  ceci


Mon petit réseau routier avec deux rond-points.
Pour information, le contenu de ma table est 
BEGIN;
-- LA LIGNE SUIVANTE EST A DECOMMENTER SI ON DOIT RECHARGER LE TOUT POUR EFFACER PROPREMENT LA TABLE AVANT RECHARGEMENT
-- SELECT dropgeometrytable('troncon_route');
CREATE TABLE "troncon_route" (gid serial,"orientation" varchar) WITH OIDS;
SELECT AddGeometryColumn('','troncon_route','the_geom','-1','MULTILINESTRING',2);
--INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING(())',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((5 0,5 6))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((0 7.5,3 7.5))',-1));
-- Premier Rond-Point
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((3 7.5,3 7,4 6,5 6))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((5 6,6 6,7 7,7 7.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((7 7.5,7 8,6 9,5 9))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((5 9,4 9,3 8,3 7.5))',-1));
-- Fin pour Premier Rond-Point
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((7 7.5,11 7.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((11 7.5,11 11))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((11 7.5,14 7.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((14 7.5,21 7.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((14 7.5,14 11))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((14 11,18 11))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((18 11,18  13))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((18 11,20 11))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((20 11,20 9))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((14 11,14 13))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((5 9,5 14.5,8 14.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((8 14.5,8 18))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((8 14.5,12 14.5))',-1));
-- Second Rond-Point
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((12 14.5,12 14,13 13,14 13))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((14 13,15 13,16 14,16 14.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((16 14.5,16 15,15.5 15.5))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((15.5 15.5,15 16,14 16))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('sens direct',GeometryFromText('MULTILINESTRING((14 16,13 16,12 15,12 14.5))',-1));
-- Fin pour Second Rond-Point
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((14 16,14 20))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((16 14.5,23.5 14.5,23.5 19))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((15.5 15.5,19 19,23.5 19))',-1));
INSERT INTO troncon_route(orientation,the_geom) values ('double sens',GeometryFromText('MULTILINESTRING((23.5 19,30 19))',-1));

END;
Pour la suite, on va appliquer les fonctionnalités nécessaires 
BEGIN TRANSACTION;
/*
   Ajouter les colonnes adéquates
*/
alter table troncon_route add column source_id int4;
alter table troncon_route add column target_id int4;
alter table troncon_route add column edge_id int4;
/*
 Ajout des colonnes x1,y1,x2,y2
*/
ALTER TABLE troncon_route ADD column x1 double precision;
ALTER TABLE troncon_route ADD column x2 double precision;
ALTER TABLE troncon_route ADD column y1 double precision;
ALTER TABLE troncon_route ADD column y2 double precision;
UPDATE troncon_route SET x1 = x(startpoint(the_geom)),y1=y(startpoint(the_geom)),x2=x(endpoint(the_geom)),y2=y(endpoint(the_geom));
/*
   Mettre à  jour le srid=1 sinon pgrouting gueule 8-(
*/
select updategeometrysrid('troncon_route','the_geom',-1);
SELECT assign_vertex_id('troncon_route',0.1);
/*
   Ok...Je crée mon graphe
*/
SELECT create_graph_tables('troncon_route', 'int4');

SELECT update_cost_from_distance('troncon_route');
/*
    Tout ce qui est à  double sens je le garde
*/
UPDATE troncon_route_edges SET reverse_cost=cost;
ALTER TABLE troncon_route_edges ADD column x1 double precision;
ALTER TABLE troncon_route_edges ADD column y1 double precision;
ALTER TABLE troncon_route_edges ADD column x2 double precision;
ALTER TABLE troncon_route_edges ADD column y2 double precision;
/*
   Copie des x1,y1,x2,y2 originaux dans la table troncon_route_edges
*/
UPDATE troncon_route_edges SET x1=tr.x1,y1=tr.y1,x2=tr.x2,y2=tr.y2 FROM (SELECT * FROM troncon_route)as tr WHERE source=source_id AND target=target_id;

/*
   Je mets maintenant à  jour les chemins à  sens direct
*/
UPDATEtroncon_route_edges SET reverse_cost=reverse_cost+1000000 
WHERE source IN SELECT id FROM troncon_route_vertices 
           WHERE geom_id IN (SELECTsource_id FROM troncon_route WHERE orientation='sens direct')) 
AND target IN (SELECT id FROM troncon_route_vertices 
             WHERE geom_id IN (SELECT target_id FROM troncon_route WHERE orientation='sens direct'));
END TRANSACTION;
Maintenant, je vérifie que tout est prêt pour la suite, un petit coup d'oeil rapide pour mon la suite 
georoute=# select * from troncon_route_edges limit 10;
 id | source | target | cost | reverse_cost | x1 | y1  | x2 | y2
----+--------+--------+------+--------------+----+-----+----+-----
  1 |      1 |      2 |    6 |            6 |  5 |   0 |  5 |   6
  2 |      3 |      4 |    3 |            3 |  0 | 7.5 |  3 | 7.5
  7 |      5 |      7 |    4 |            4 |  7 | 7.5 | 11 | 7.5
  8 |      7 |      8 |  3.5 |          3.5 | 11 | 7.5 | 11 |  11
  9 |      7 |      9 |    3 |            3 | 11 | 7.5 | 14 | 7.5
 10 |      9 |     10 |    7 |            7 | 14 | 7.5 | 21 | 7.5
 11 |      9 |     11 |  3.5 |          3.5 | 14 | 7.5 | 14 |  11
 12 |     11 |     12 |    4 |            4 | 14 |  11 | 18 |  11
 13 |     12 |     13 |    2 |            2 | 18 |  11 | 18 |  13
 14 |     12 |     14 |    2 |            2 | 18 |  11 | 20 |  11
(10 lignes)
Et pour ce qui est des rond-points qui doivent avoir un cost différent de leur reverse_cost 
georoute=# select * from troncon_route_edges where cost != reverse_cost;
 id | source | target |       cost       |   reverse_cost   |  x1  |  y1  |  x2  |  y2
----+--------+--------+------------------+------------------+------+------+------+------
  5 |      5 |      6 | 2.91421356237309 | 1000002.91421356 |    7 |  7.5 |    5 |    9
 23 |     21 |     22 | 1.70710678118655 | 1000001.70710678 | 15.5 | 15.5 |   14 |   16
  3 |      4 |      2 | 2.91421356237309 | 1000002.91421356 |    3 |  7.5 |    5 |    6
 22 |     20 |     21 | 1.20710678118655 | 1000001.20710678 |   16 | 14.5 | 15.5 | 15.5
 20 |     19 |     16 | 2.91421356237309 | 1000002.91421356 |   12 | 14.5 |   14 |   13
  4 |      2 |      5 | 2.91421356237309 | 1000002.91421356 |    5 |    6 |    7 |  7.5
 21 |     16 |     20 | 2.91421356237309 | 1000002.91421356 |   14 |   13 |   16 | 14.5
  6 |      6 |      4 | 2.91421356237309 | 1000002.91421356 |    5 |    9 |    3 |  7.5
 24 |     22 |     19 | 2.91421356237309 | 1000002.91421356 |   14 |   16 |   12 | 14.5
(9 lignes)
Pour la suite, je vais créer deux tables aller et retour qui sont basées sur troncon_route à  partir duquel j'ai récupéré les gid qui valent respectivement 1 et 27 
begin transaction;
create table aller(gid int4) with oids;SELECT AddGeometryColumn( 'aller', 'the_geom', -1, 'MULTILINESTRING', 2 );
create table retour(gid int4) with oids;SELECT AddGeometryColumn( 'retour', 'the_geom', -1, 'MULTILINESTRING', 2 );


insert into aller(the_geom) (select the_geom from troncon_route where edge_id in (select edge_id from 
shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',(select id 
from troncon_route_vertices where geom_id in (select target_id from troncon_route 
where gid=27 limit 1)),(select id from troncon_route_vertices where geom_id in (select source_id from troncon_route where gid=1 limit 1)),false,false)));

insert into retour(the_geom) (select the_geom from troncon_route where edge_id in (select edge_id from
 shortest_path_astar('select  id,source::int4,target::int4,cost::double precision,
reverse_cost::double precision as reverse_cost, x1::double precision,y1::double precision,
x2::double precision,y2::double precision from troncon_route_edges',(select id 
from troncon_route_vertices where geom_id in (select target_id from troncon_route where gid=27 limit 1)),(select id from troncon_route_vertices where geom_id in (select source_id from troncon_route where gid=1 limit 1)),false,true)));
end transaction;
Au niveau de la visualisation, je vais obtenir

Parcours à  l'aller.

Parcours au retour.
Et pour la route, un autre exemple pour se convaincre que mon graphe orienté marche bien

Parcours à  l'aller.

Parcours au retour.

2 commentaires :: aucun trackback