Le portail francophone de la géomatique

Merci Nicolas pour cette remarque. En effet, l'utilisation de st_dwithin simplifie grandement la première requette et va m'être bien utile ! Très intéressant.

Par contre, concernant la méthodologie et mes interrogations, cela ne change rien au résultat obtenu.

Dernière modification par jerome_m (Thu 23 October 2014 10:54)

Une question:

Comment créez-vous les tronçons manquants ? par dessin ou par jointure entre les segments consécutifs ?

Vous pouvez choisir, pour chaque troncon de ligne ordonné, le segment de route le plus proche.
Ensuite, vous reconstruisez la ligne finale, basée sur bd_topo, order by troncon_ligne.

Nico

Bonjour,

Postgis sait aussi faire de l'accrochage auto entre geometries, dans certaines conditions.

La methode est fonction du type de données sur lesquelles vous travaillez, d'ou ma question sur un echantillon test:
Sont-ce des linestring simples, des multi, comment se recoupent ou se chevauchent-elles ?

Vous pouvez nous envoyer qq données, eventuellement anonymisées si nécessaire ?

Nico

C'est un problème qui revient assez souvent : comment faire pour "coller" deux couches qui vues de loin se superposent mais, en fait, ne se supperposent pas du tout; et  où l'une est + précise, l'autre contient des info intéressantes....

Ca fait également penser à la BDTOPAGE qui tente de caler la BD CARTHAGE sur la BDTOPO. Ou à des pistes cyclables numérisées par un cycliste amateur qui ne se superpose pas vraiment avec le référentiel voirie de commune !

Dès lors qu'on souhaite faire des manip automatisée, il faut choisir des seuils : comment les définir ? ST_DWITHIN ( geom , ????) ?

Avec PostGIS, il y a comme le dit Nicolas, un tas de solutions, et ca dépend bien souvent du jeu de données, d'arriver à formuler tous les cas possibles, à éventuellement à créer un indice de 0 à 10 qui tente d'indiquer le % d'erreur à avoir caler une ligne sur la ligne référentielle (pour ensuite faire des corrections manuelles éventuelles).

Mon post ne fait pas vraiment avancer les choses mais la question m'intéresse !

Bon courage,

Dernière modification par ppluvinet (Thu 23 October 2014 18:59)

Merci,

Vous auriez aussi les routes BD TOPO qui vont avec ?

Pour des tracés compliqués, comme le disait Pascal, ca peut vite devenir très complexe de coller à la couche de référence: que faire aux intersections, comment choisir le bon segment qd plusieurs segments de référence sont éligibles, etc.

Nicolas

La solution avec PgROUTING, je n'y avais pas pensé mais y'a de l'idée !

Une autre solution (à tester car je ne l'ai pas fait) pour tenter de repérer les erreurs de mauvais raccordement..
Vous découpez votre linestrind en tronçons de 50m (la distance est un seuil au pif donc à recaler peut-être) et vous numérotez bien les tronçons de 1 à N. Appelons ces toncons "TRONCLINE".
Puis vous créez sur chaque TRONCLINE, 10 points espacés de 5 m .
Puis vous cherchez pour chaque point le tronçons de BDTOPO le + proche.
Par agrégation vous pouvez connaitre le nombre de point qui se sont accrochés au même TRONCLINE. Cela donne un indice de confiance. Si j'ai 9/10 points qui se sont accrochés au même tronçon j'ai plus de chance que si j'en ai que 5. Si l'indice est faible, alors c'est que c'est le cas d'une intersection ou qu'il faut une reprise manuelle.


Sur Postgis, je pense que ca doit être un peu prise de tête mais carrément faisable.

Bon courage,

Dernière modification par ppluvinet (Thu 23 October 2014 18:58)

@Nicolas Ribot : avec ST_dwithin(), si je met un buffer de recherche de quelques mètres, cela fonctionne très bien : en effet, mon tracé initial est très proche de ma couche de référence ce qui simplifie quand même cette étape. Pour les zones sans correspondance, comme précisé au début, je créerai un nouveau tronçon perso.

@Nicolas Granier : votre approche est très intéressante, et elle fait ressortir la notion d'itinéraire, reliant des points ordonnés par des tronçons de route. L'algorithme me semble correct, et si je comprends bien, il me suffirait de créer quelques points de passage pour "spécifier" le sens de parcours pour les zones litigieuse (boucles ou zones de passages dans les deux sens). Ensuite, il sera possible de créer une LineString correcte.


@Pascal PLUVINET : c'est une solution à tester dans le cas où le tracé initial est assez éloigné de la couche référence. Le gros point positif est l'indice de confiance qui permet de choisir dans le cas où plusieurs tronçons sont identifiés.

@tous: suite à ces discussions, je me pose une question concernant la visualisation des LineString et des itinéraires avec QGis. Est il possible de les visualiser en représentant le sens de passage et l'ordre des tronçons ? (avec des petites flèches par exemple, un peu comme dans les logiciels permettant de visualiser des itinéraires).

Merci à tous pour votre participation

Bonjour,

Merci pour les données routes: J'etais parti avec des données OSM, non calées sur la ligne => me suis embarqué dans de la reprojection de points pour d'abord caler la ligne sur OSM, puis ensuite retrouver l'ordre des troncons OSM: bien tricky mais faisable

En voyant vos données, on constate que la ligne a ete numérisée sur bd_topo => les sommets de la ligne sont déja calés sur bd_route, les deux tracés étant superposés exactement (sauf en un point ou il y a un décalage de 10cm).
On constate aussi que BD route est topologiquement correcte: segments individuels connectés aux intersections.

Du coup, le processus pour retrouver les segments BD TOPO dans le bon ordre et suivant la ligne est plus facile: on va se servir de cette caractéristique.

Au passage, je vous ai dit une "bêtise": dans votre cas, st_buffer + st_contains pour trouver toutes les routes qui suivent la ligne est la bonne méthode: st_dwithin avec deux lignes peut renvoyer des lignes proches mais pas completement incluses dans le buffer: ca génère trop de faux positifs

Le principe:

• Explosion de la ligne en segments individuels numérotés, a partir des sommets de la ligne: ca va servir a ordonner bd route à la fin:
=> table vertex (sommets de la ligne) et segments (2 points consécutifs dans la ligne)

• Petit Buffer de la ligne (1m par ex, les données étant calées) pour retrouver tous les segments BD Route + la ligne de départ + la ligne d'arrivée (qui ne sont pas contenues dans le buffer, mais qui l'intersectent, sauf si le trajet ligne commence au début d'une route BD Topo).
=> table route_tmp: les routes bd topo suivant la ligne

• Ensuite, on associe a chaque segment le troncon bd_route qui le contient: on a avoir des doublons: plusieurs segments de la ligne appartiennent parfois au meme troncon bd_route.

• Enfin, on utilise une CTE recursive pour identifier les segments bd routes successifs, non répétitifs, qui forment la ligne. Avec cette méthode, on ne detecte pas un demi tour fait sur le meme troncon BD Route, mais on détecte les boucles (plusieurs passages sur le meme troncon).

En SQL et en image:

creation de la table des sommets de la ligne, ordonnables par le champ idx

drop table IF EXISTS vertex;
create table vertex as (
    -- multilinestring => 2eme element de path a extraire pour avoir l'index du sommet de la ligne
    select gid, (st_dumppoints(geom)).path[2] as idx, (st_dumppoints(geom)).geom
  from ligne
);

self join sur la table des vertex pour construire un segment entre les points n et n+1:

create table segments as (
  select t.idx, st_collect(t.geom, t1.geom) as geom
  from vertex t, vertex t1
  where t.idx <> t1.idx
  and t1.idx = t.idx + 1
);

La ligne de départ
[img]http://i59.tinypic.com/11snkeg.png[/img]

table des routes bd_topo suivant la ligne:
le UNION permet d'ajouter les troncons de départ et d'arrivée a cette table, qui contient pour chaque segment
de ligne le segment bd route qui lui correspond

drop table route_tmp;
create table route_tmp as (
  SELECT
    r.gid AS route_gid,
    r.geom
  FROM route r, ligne l
  WHERE st_contains(st_buffer(l.geom, 1), r.geom)
  UNION
  -- route de depart
    select r.gid, r.geom 
    FROM route r, vertex v 
    where v.idx = 1 and st_dwithin(v.geom, r.geom, 0.01)
  UNION 
  --route d'arrivée
    select r.gid, r.geom 
    FROM route r, vertex v 
    where v.idx = (select max(idx) from vertex) 
    and st_dwithin(v.geom, r.geom, 0.01)
    
);

Routes contenant chaque segment: on crée un identifiant croissant avec row_number, différent de s.idx
qui peut ne pas etre continu si on ne retrouve pas un segment de bd topo:
Cet id servira dans la fonction récursive pour avancer dans le sens de la ligne de départ

drop table if exists seg_route;
create table seg_route as (
    SELECT
      s.idx,
      r.route_gid,
      row_number() OVER (ORDER BY s.idx) as id,
      r.geom
    FROM segments s, route_tmp r
    WHERE st_dwithin(r.geom, st_geometryN(s.geom, 1), 1)
          AND st_dwithin(r.geom, st_geometryN(s.geom, 2), 1)
);

Enfin, la fonction récursive: pour chaque tronçon consécutifs, on repère quand on a affaire à un nouveau troncon et on note
son identifiant, ou null si l'identifiant est le meme que le précédent: on est tjs sur le meme troncon de bd route.
A la fin de la recursion, on ne garde que les routes ayant un identifiant non null et on regénère un id croissant, qui sera le numero de segment
ordonné par rapport à la ligne initiale

with recursive tmp as (
  select idx, id, route_gid, route_gid as newid
  from seg_route
  where id = 1
  
  UNION 
      
  select s.idx, s.id, s.route_gid,
    case when t.route_gid <> s.route_gid then s.route_gid else null end as newid
  from seg_route s, tmp t
  where s.id = t.id + 1
  
) select row_number() over (ORDER BY id) as path, r.geom 
  from tmp t join route_tmp r on t.newid = r.route_gid
where newid is not null;

On obtient une série de segments de bd route ordonnés dans le sens de parcours de la ligne:
En noir la ligne initiale, en bleu les segments BD route retrouvés avec leur ordre

[img]http://i60.tinypic.com/2zhlefr.png[/img]

Les boucles sont gérées: les troncons à plusieurs passages sont dupliqués dans le résultat final:

[img]http://i60.tinypic.com/2nsrvv4.png[/img]

Nico

Enorme merci Nicolas pour ce travail ! Je pense qu'il va servir à de nombreuses personnes vu le peu d'info disponibles sur internet par rapport à cette problématique !
Je reprendrais tout cela la semaine prochaine pour bien l'annalyser et l'appliquer à toutes nos lignes, mais dans tous les cas, j'ai énormément appris avec ces discussions sur les LineString et cela va nous permettre de créer et d'organiser nos données d'une façon efficace.

Pour les "trous" observés dans la BD route, il n'y en a pas dans cet échantillon (volontaire). Cela se produit dans le cas de troncons manquants, (principalement des nouvelles routes par rapport au référentiel utilisé). C'est pour cela qu'il n'est pas nécessaire d'augmenter le buffer de recherche, cela ne servirait à rien.
Je vais essayer de tester rapidement comment se comporte cette méthode lorsqu'il y a des tronçons manquants, et trouver une solution pour répondre entièrement au besoin.

Merci aussi Pascal pour les informations sur le style d'affichage des LineString !

Dernière modification par jerome_m (Fri 24 October 2014 14:37)