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 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 >

dimanche 1 avril 2007

PgRouting - test sur un jeu de données réelles - 3/3

Ici, je vais considérer le réseau suivant que je me suis amusé à extraire de données réelles. Avant de commencer les tests, j'ai du effectuer quelques modifications mineurs des données. Je ne cache pas non plus que j'ai du ici vérifier pour chaque noeud inverse au sens de circulation.

Pré-traitement: tronçons de rond-poind à l'envers.
J'ai eu avec ces données quelques tronçons de rond-points (MULTILINESTRING) dont l'acquisition géométrique a été faite dans le sens opposé au sens de circulation . Rien de bien méchant. J'ai pu résoudre celà grâce à la fonction reverse() de PostGIS :
UPDATE troncon_test SET the_geom = reverse(the_geom)  WHERE gid in (74,75,174,175,186,54);
Mais bon là il s'agit d'un petit jeu. Si j'avais du effectuer celà sur un énorme jeu, bonjour la galère! Mais bon quand on vous livre des données rien ne garantit qu'elles soient déjà configurées comme on voudrait.
Mon réseau et mes noeuds.

Fig 1. Le réseau routier.

Dans QGIS, je vais charger les noeuds que j'ai crée en suivant les mêmes commandes SQL de mon précédent billet

Fig 2. Une portion des noeuds du réseau.

Exemple 1: avec deux rond-poins
Dans un premier temps, on va essayer d'aller du noeud 188 au noeud 143. Là je dois traverser deux rond-points comme le montrent les images suivantes

Fig 3. Parcours aller - image au un niveau du réseau.


Fig 3. Zoom dans le réseau, parcours aller et retour.
Exemple 2: un rond-point et une route à sens unique.
Maintenant je vais aller du noeud 167(=source) au noeud 76(=target) de mon réseau. Les images pour le parcours aller/retour vue à l'échelle du réseau sont présentées sur les figures n° 4.a et 4.b.

Fig 4.a Parcours aller.



Fig 4.b Parcours retour
Maintenant dans le sens du retour, j'ai

Fig 5.a r.Zoom sur le noeud 76. Le rond-point est pris dans le bon sens.



Fig 5.b Zoom sur le noeud 76. Le rond-point est pris dans le bon sens.

Conclusion: La fonction shortest_path_astar a bien été éprouvé - comme on dit - ces deux derniers jours. L'utilisation de cette fonctionnalité de pgRouting s'avère des plus efficaces. Je ne regrette pas de l'avoir tester. A un moment, j'ai eu un doute quand j'ai vu notamment certains rond-points pris en sens inverse. Mais bon celà marche. C'est le plus important .

samedi 31 mars 2007

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

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 .

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

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

jeudi 29 mars 2007

Avoir des noms de fichiers en lettres minuscules

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.