Le portail francophone de la géomatique

-- La requete recursive (on met les resultats dans une table)
-- Cette requete recherche tous les chemins possibles a partir d'un troncon
-- donne en suivant les troncons dans leur sens naturel (de la source a la target).
drop table rec_res;
create table 
    rec_res as 
with recursive 
    -- la CTE
    search_graph(id, target, depth, path, cycle) as (
        -- Initialisation
        -- on part d'un troncon specifique
        select 
            g.id, g.target, 1 as depth, ARRAY[g.id] as path, false as cycle 
        from 
            edge_table as g 
        where 
        -- mettre ici les conditions des troncons initiaux desquels on part
            source = 6
        union all
        -- Partie récursive
        select
            g.id
            -- on a besoin de garder la target pour pouvoir jointurer a la prochaine iteration
            , g.target
            -- on incremente la profondeur a chaque iteration
            , sg.depth + 1 as depth
            -- on ajoute le troncon en cours dans le chemin deja parcouru
            , path || g.id as path -- l'opérateur || permet d'ajouter la valeur au tableau
            -- on calcule la longueur du chemin au fur et a mesure
            --, sg.length + g.length as length
            -- on verifie si on est deja passe par ce troncon
            , g.id = ANY(path) as cycle
        from
            -- la table qu'on jointure : c'est le graphe (au format pgrouting)
            edge_table as g
        join
            -- la CTE
            search_graph as sg
        on
            -- on jointure la CTE et edge_table pour se déplacer dans le graphe
            -- nos chemins sont crees en avancant dans le sens naturel des troncons
            -- on ne revient donc pas en arriere
            sg.target = g.source
        where
            -- on ne prend pas en compte les eventuels cycles
            not cycle
        -- si on veut s'arrêter dans la propagation il faut ajouter des conditions ici. Par exemple sur l'attribut length ( voir dans la présentation)
)
select 
    -- tous les champs de la CTE
    sg.*
    -- on renvoie egalement la geometrie pour pouvoir visualiser dans qgis
    -- pour cela on jointure avec la table du graphe
    , edge_table.the_geom
from
    -- La CTE
    search_graph as sg
join
    edge_table
on
    sg.id = edge_table.id
-- on limite toujours les resultats pour ne pas avoir de boucle infinie meme
-- si on s'est trompe dans les criteres d'arret
limit 10000;