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 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 >

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

Un salon multimédia Linux dans le Nord de la France - 28 et 29 avril

Par david techer, mardi 27 mars 2007 à 15:24 :: Linux

Un salon multimédia dans le Nord de la France à Proville! Purée ça a l'air intéressant. http://linuxcambresis.org/wp-content/uploads/2007/03/dossier-presse-salon2007.pdf. Si j'avais les moyens d'y aller ce serait avec un énorme plaisir mais mon calendrier risque de ne pas me le permettre!

aucun commentaire :: aucun trackback

Petit bidouille éléctrique sur un adaptateur péritel pour avoir la couleur sur une télé, reliée la sortie TV-out/S-Vidéo du PC.

Par david techer, mardi 27 mars 2007 à 03:22 :: Linux

Hier soir, j'avais décidé de brancher mon PC à ma TV. J'ai pour celà pour bien faire utiliser mon câble S-vidéo branché sur la sortie s-video de mon ordi de bureau. La jonction entre le câble svidéo et la prise péritel de TV se fait à l'aide d'un adaptateur péritel. Malheureusement la sortie sur l'écran était en noir et blanc. J'ai cru au début à une erreur Nvidia. En allant rechercher sur des forums dédiés, j'ai trouvé la solution, le super-lien suivant: http://forum.generation-nt.com/faq-materiels/sortie-tv-noir-and-blanc-151616/ . Au lieu de souder ou d'essayer de joindre par un petit fil de cuivre - pas si facile que ça - la borne 15 à la borne 20, j'ai tout simplement pris mon briquet pour brûler une petite portion de la gaine en plastique afin de dénuder une petite portion des deux fils qui sont reliés aux bornes 15 et 20. J'ai ensuite torsadé les petites portions de cuivres entre elles. J'ai ensuite pris une petite portion de scotch afin de faire tenir le tout.

Et ça marche nickel ce petit truc . Vraiment sympa cette petite astuce!!!.

2 commentaires :: aucun trackback