Le portail francophone de la géomatique

Bonjour

il me semble qu'avec une formule de ce type, cela calcul la longueur exact à l'intersection des polygones communes

SELECT commune.nom, sum(st_length(st_intersection(assainissement.geom, commune.geom))) AS ml_cana
FROM administratif.commune, technique.assainissement
WHERE st_intersects(commune.geom,assainissement.geom)
GROUP BY commune.nom

J'ai repris cette requête d'une requête plus générale donc elle est a adapter à vos besoins.
Cordialement.

Bonjour,

Merci pour ces éléments ! La requête donnée par Benoît ne renvoie pas le linéaire exact par commune mais le résultat est légèrement différent.

Pour résumer :


Scénario 1 - Ma requête compte en double les lignes qui chevauchent les contours de communes :

SELECT a.nom, st_length(b.geom) as longueur FROM commune a, reseau b WHERE st_intersects(a.geom, b.geom) ;

Scénario 2 - La requête proposée par Benoît tranche et associe chaque ligne à cheval à une commune en particulier :

SELECT a.nom, sum(st_length(st_intersection(a.geom, b.geom)) AS longueur FROM commune a, reseau b WHERE st_intersects(a.geom, b.geom) GROUP BY a.commune ;

Nicolas, comment t'y prendrais-tu pour intégrer cette découpe dans la requête (avec la structure WITH table AS...) ?


Merci encore pour votre aide !

Thomas

Dernière modification par Thomas Williamson (Thu 17 October 2013 14:47)

Ne serait-ce pas votre group by qui pose problème ?

Vous faites un select a.nom et un group by a.commune

Cordialement.

Par calcul de linéaire, vous souhaitez obtenir la somme des longueurs de réseau contenue dans chaque commune ?

Si c'est le cas, il y a deux options concernant les lignes confondues avec une frontiere de commune:
    • comptabiliser leur longueur pour chaque commune, (la requete de Benoit)
    • comptabiliser leur longueur une seule fois, en attribuant le segment a une seule commune: il faut alors un critère pour choisir a quelle commune attribuer un segment commun.

Nicolas

Je pense que la remarque de Nicolas voulait dire : quel choix faites-vous lorsque la ligne longe (entendre, se superpose à) une portion de votre limite communale? Pour le coup, ce linéaire touche les deux communes. La requête proposée par Benoît compte cette portion dans les deux communes.
Après, peut-être que cette situation est marginale dans vos traitements... Mais c'est à savoir.

Bonjour,

Je me permets de réanimer ce sujet car je suis dans un cas similaire : je souhaite moi aussi obtenir la longueur d'un réseau (cours d'eau) par communes et la fonction que j'utilise est la même que celle plus haut avec

sum(st_length(st_intersection(a.geom,b.geom)))

et elle fonctionne très bien pour les cours d'eau qui "entrent" une fois par communes et la traversent par exemple.

Néanmoins, certains cours d'eau en bordure de communes longent la limite, entrant et sortant plusieurs fois de chaque communes (en zigzag) et là la longueur est nulle.

Auriez-vous une idée ? Je sais que la géométrie "en morceau" de la ligne sera discutable, mais j'aimerais bien obtenir la longueur par communes malgré tout.

Merci à tous

Emma