Le portail francophone de la géomatique

J'ai trouvé une méthode qui fonctionne:

Couches en entrée:
- Couche de lignes avec un identifiant unique FID
- MNT (couche nommée altitude)

Étapes:
1) Création de points le long des lignes avec le traitement QGIS QGIS Create points along lines, avec une distance entre les points X
2) Calculatrice de champs: création sur la couche de points d'un champ num_point avec la formule suivante:

"Distance"  / X

3) Calculatrice de champs: création sur la couche de points d'un champ id_point avec la formule suivante:

"FID" || '_' || "num_point"

4) Récupération des valeurs d'altitude sur les points avec le traitement SAGA Add grid values to points
5) Calculatrice de champs: calcul du dénivelé par rapport au point précédent dans un nouveau champ deniv avec la formule suivante:

CASE WHEN num_point = 1 THEN 0
ELSE     
attribute( 
     get_feature( 
        @layer_id,
        'id_point',
        "FID" || '_' || ("num_point" -1)
     ),
     'altitude'
 ) - "altitude"
END

6) Calculatrice de champs: création d'un nouveau champ sens_deniv avec la formule suivante:

"deniv" / abs( "deniv" )

7) Group Stat pour calculer les dénivelés positif et négatif cumulés avec les paramètres suivants:
Columns: sens_deniv
Rows: FID
Value: deniv, sum
8) Export du tableau en .csv et jointure avec la couche de ligne sur le champ FID pour récupérer les dénivelés.


Ça fonctionne même si c'est un peu laborieux. J'ai essayé d'automatiser le tout en passant par le modeleur graphique mais je bloque à l'étape 5 car l'appel à la variable @layer_id ne semble pas fonctionner au sein du modeleur. Dans tous les cas l'étape 7 ne peut pas se faire avec le modeleur du fait de l'appel à Group Stats.

Si vous avez des améliorations, je suis preneur.

Dernière modification par gustavecoste (Tue 24 May 2016 13:06)

gustavecoste, merci beaucoup pour ta technique. Aujourd'hui encore, je n'ai pas trouvé d'extention sur QGis permettant de calculer le denivelé positif.

Cependant, dans ta formule de "deniv", le champ me retourne 1 lorsque je perds 1 mètre, et me retourne -1 lorsque j'en gagne 1. Ca ne devrait pas être l'inverse ?

Je suis obligé de faire un "deniv * -1" après, pour que ça soit cohérent.

Merci à toi pour ta réponse.

Malheureusement, je n'arrive pas à trouver l’algorithme "Clim along line". Est-il trouvable dans l'onglet Traitement ?