La chaîne de traitements pour la génération des axes de ruissellement concentré a été développée par l’association temporaire Alexandre Maugnard - Champs Libres Coopérative. Elle se base sur le MNT LiDAR 2013-2014. Celui-ci a été agrégé en 2m de résolution, gravé, rendu hydrologiquement continu et ensuite transformé en un MNT où les valeurs d’altitude sont effacées au niveau des dépressions, du réseau hydrographique et des surfaces des grands cours d’eau. C’est sur ce dernier MNT qu’est utilisé l’algorithme de direction de flux répartissant le flux d’amont en avant entre les pixels du MNT. Cet algorithme permet de définir un pattern du réseau de drainage en surface. Pour la génération de la couche de variation possible du tracé des axes de concentration du ruissellement, un algorithme de direction de flux multiple a été utilisé.


Le GMPD est un algorithme de direction de flux basé sur le gradient maximal de pente descendante. En effet, pour pallier au problème du coefficient de partition du flux qui est fixe dans la plupart des algorithmes de direction de flux multiple, Qin et al. (2007)* proposent un algorithme dont le coefficient de partition varie en fonction du gradient maximal de pente descendante. L'idée sous-jacente est d'adapter le degré de dispersion du flux au relief rencontré. La fonction développée par Qin et al. (2007) permet de limiter la divergence là où les gradients de pente sont élevés alors que la divergence est de plus en plus importante lorsque le gradient de pente diminue. En d'autres termes, il s'agit d'un algorithme de direction de flux adaptatif qui doit être moins sensible à de soudains changements de relief que les algorithmes de direction d8 classiques.
La couche de variation possible du tracé obtenue est traitée par un filtre. En effet, après classification du raster des écoulements GMPD, celui-ci est « tamisé » (gdal_sieve) ce qui permet d’éliminer des groupes de pixels isolés. Ensuite, la couche finale de variation possible du tracé est une combinaison du raster classifié et du raster tamisé. Les pixels ayant une valeur de flux supérieure au seuil de classification et qui ont changé de valeur par le tamisage sont supprimés dans la couche de variation possible du tracé, tandis que les pixels ayant une valeur inférieure au seuil de classification et qui ont changé de valeur par le tamisage sont conservés. Cela évite que l'étape de tamisage par « gdal_sieve » créée des valeurs de variation possible du tracé égale à 1 là où le raster classifié n'indiquait pas de valeur égale à 1.

* Qin, C. Z., Zhu, A. X., Pei, T., Li, B. L., Scholten, T., Behrens, T., & Zhou, C. H. (2011). An approach to computing topographic wetness index based on maximum downslope gradient. Precision Agriculture, 12(1), 32-43.