Étant donné que vous avez le LIDAR DEM, vous devez utiliser les flux qui en dérivent. Cela garantit un enregistrement parfait.
L'essentiel de l'idée est d'estimer les pentes moyennes en fonction des élévations aux extrémités des segments.
L'une des procédures les plus simples consiste à «exploser» le réseau de flux en ses arcs non ramifiés composants. Convertissez la collection en une couche "route" basée sur la distance, la rendant "mesurable". Il est maintenant simple de générer une collection d '"événements" d'itinéraire basés sur un tableau de jalons (à des intervalles de 100 m par exemple) pour chaque arc et d'extraire les élévations DEM de ces points d'événement. Des différences d'élévation successives le long de chaque arc, divisées par 100 m, estiment les pentes moyennes des segments.
La figure suivante cartographie les arcs de cours d'eau dérivés d'une analyse d'accumulation de flux d'un DEM USGS de 7,5 minutes (partie du comté de Highland, VA). C'est environ 10 km de diamètre (6 mi).
Étant donné que vous recherchez un barrage restant, ce qui pourrait être indiqué par un changement de pente sur quelques dizaines de mètres (pour un très petit barrage), pensez à utiliser des segments encore plus petits . Si le jeu de données est trop approximatif pour fournir des signaux clairs, vous pouvez facilement le filtrer plus tard (au moyen de moyennes mobiles ou autrement, telles que les courbes de splining des élévations et la différenciation de la spline). En effet, cette approche vous place dans le domaine de l'analyse des séries chronologiques où la variable d'intérêt est l' élévation, pas le gradient, et vous recherchez des modèles composés de courtes sections de niveau suivies de changements soudains.
Il s'agit d'un tracé des élévations DEM observées à des intervalles de 100 m le long de la plupart (pas tous) des segments de cours d'eau représentés. (La taille des cellules est de 30 m.) Si nécessaire, les arcs ont été réorientés pour que l'élévation diminue généralement de gauche à droite. (Si vous regardez attentivement, vous pouvez voir où j'en ai manqué un: il monte de gauche à droite.)
Ce détail de l'arc 16 (le long segment en haut de la carte) montre ce que vous pourriez obtenir lorsque les flux ne sont pas parfaitement enregistrés avec le DEM: par endroits, le flux semble couler vers le haut. Néanmoins, les segments suggérant des caractéristiques de bassin et de chute sont facilement identifiables, en particulier après les jalons 1800 (mètres le long du segment), 4000, 4600 et 6500. Cette identification peut être automatisée de diverses manières, en particulier après le nettoyage de la série d'élévation (par lissage il).
Vous pouvez voir que l'intervalle d'échantillonnage de 100 m utilisé ici n'est vraiment pas assez bon pour identifier des entités beaucoup plus petites de 400 à 500 mètres de long. Donc, pour trouver un petit barrage restant, vous voudrez probablement échantillonner autour d'un intervalle de 10-25 m sur votre LIDAR DEM.
BTW, ce qui rend un segment de flux "trop petit" pour ce genre de travail n'est ni une courte longueur ni une grande taille de cellule, bien que les deux jouent un rôle dans la décision. «Trop petit» dépend de la façon dont vous utiliserez les pentes estimées et du degré d'incertitude de ces estimations. Pour certains travaux, il pourrait même être judicieux d'estimer les gradients à des intervalles de 10 m sur une grille de 10 m!