Le portail francophone de la géomatique

Bonjour,

Je veux calculer la distance en passant par les routes entre 2 villes j'avais fait une première application qui marchait mais qui était pas très précises.
Pensant comprendre comment ça marchait, j'ai mis en place une procédure pour adapter la base  à la fonction dijkstra_sp,  c'est à dire l'implémentation par pgrouting de l'algorythme de dijkstra, d'une façon que je pensais plus propre...
Mais ça ne marche pas : quand j'utilise dijkstra_sp ça ne me donne pas de résultat :'(
Je sais que mon erreur est limité aux tables routes et vertices_tmp (crée par la fonction  assign_vertex_id() ) j'ai conscience que  c'est une piètre circonscription du problème...
Je vous fais part des successions de transformations que j'ai fait j'espère que vous y verrez plus clair que moi...

voici comment j'importe le fichier france.osm (trouvé ici : http://download.geofabrik.de/osm/europe/) :

E:\bddSpaciale\osm2pgsql\osm2pgsql.exe -d France -U postgres --slim -S "E:\bddSpaciale\osm2pgsql\default.style" -P 5432 "E:\bddSpaciale\données\geofabrique\france.osm"

voici les transformations qui "nous intéressent"

UPDATE planet_osm_roads SET way = ST_Transform(ST_SetSRID(way,900913), 4326);
update geometry_columns set srid = 4326 where f_table_name = 'planet_osm_roads' ;

create table routes as (select postal_code, name, place, way as the_geom from planet_osm_roads);


alter table routes add column gid serial;

ALTER TABLE routes ADD COLUMN source integer;

ALTER TABLE routes ADD COLUMN target integer;

ALTER TABLE routes ADD COLUMN length double precision;

insert into geometry_columns values ('', 'public', 'routes', 'the_geom', 2, 4326, 'LINESTRING');


SELECT assign_vertex_id('routes', 0.001, 'the_geom', 'gid');
UPDATE routes SET length = length(the_geom);


CREATE INDEX source_idx ON routes(source);
CREATE INDEX target_idx ON routes(target);
CREATE INDEX geom_idx ON routes USING GIST(the_geom GIST_GEOMETRY_OPS);



UPDATE routes SET the_geom= ST_Transform(ST_SetSRID(the_geom,4326), 2154);
alter table vertices_tmp drop constraint enforce_srid_the_geom;
UPDATE vertices_tmp SET the_geom = ST_Transform(ST_SetSRID(the_geom,4326), 2154);


update geometry_columns set srid = 2154;

Voilà le problème, si je prend 2 noeuds 8898 et 55665 respéctiement des noeuds proches de Rennes et de Strasbourg tel que :

SELECT sum(length(the_geom)) 
    FROM dijkstra_sp('routes', (8898) , 
                (55665) );

ça me met une réponse vide :'(
Si je prend des nœuds très proches (~1km) j'ai une réponse mais j'ai pas vérifié si elle collait à la vérité...

Je vous remercie d'avance d'essayer ou même simplement me dire d'où peut venir le problème juste pour me donner une piste je suis en train de refaire la base pour vérifier mais ça met beaucoup de temps et j'en ai pas :s

Dernière modification par Syphys (Thu 10 June 2010 11:52)