La chaîne de traitements pour la génération des axes de ruissellement concentré a été développée par "Earth and Life Institute" (UCL - ELIe). Elle se base sur le MNT LiDAR 2013-2014. Celui-ci a été lissé et rendu hydrologiquement continu.

La couche de données se base sur la couche d'accumulation de flux dérivé de l'utilisation de l'algorithme à direction multiple MD8. Le seuillage des valeurs d’aire contributive dérivées du MD8 permet de définir une emprise pour chaque axe de concentration correspondant à l’incertitude sur le positionnement de l’axe mais qui peut également refléter l’étalement potentiel de l’axe. Plus le seuil choisi est élevé, plus la zone d’incertitude est étroite. A contrario, l’utilisation d’une valeur trop faible conduit à définir des zones d’incertitude de taille importante et à mettre en évidence des zones qui ne contiennent pas forcément des axes de ruissellement. Il semble qu’une valeur de 0.5 ha soit un bon compromis car elle permet de définir des zones resserrées là où les axes sont relativement certains et de mettre en évidence les zones où la dispersion intervient.

La couche des incertitudes à été générée à partir des solutions WhiteBox GAT, ArcGIS 10.2 et des scripts python associés (module arcpy).

Deux artéfacts sont à noter :

1. Le premier est lié au fait que la couche d’accumulation de flux dérivée de l’utilisation de l’algorithme à direction de flux multiple MD8 dans WhiteBox GAT ne peut être obtenue qu’au départ d’un MNT et non pas au départ d’une couche de direction de flux. De ce fait, il n’est pas possible d’effacer les directions de flux au niveau des cours d’eau et des zones bâties. On obtient donc des valeurs d’accumulation de flux dérivées du MD8, et par conséquent des zones d’incertitudes, au niveau des cours d’eau et des zones bâties. Ces dernières sont en partie éliminées en supprimant les valeurs d’accumulation de flux dérivées du MD8 au droit des cours d’eau et dans les zones bâties. Toutefois, des zones d’incertitude autour des cours d’eau subsistent et il subsiste des zones d’incertitude en aval de zone bâties et qui ne correspondent à aucun axe;

2. Le deuxième artéfact est lié au fait qu’il n’est pas possible de définir les valeurs de direction de flux dérivées du MD8 sur l’ensemble du territoire couvert par les données du MNT LiDAR en un seul bloc ou même par sous-bassin hydrographique lorsque l’on utilise WhiteBox GAT. Ceci aboutit en effet à une erreur mémoire. Nous avons donc procédé de la manière qui suit : le MNT hydrologiquement continu sur l’ensemble du territoire couvert par les données MNT LIDAR a été découpé selon l’emprise des 4551 mailles avec un buffer de 200 m puis, pour chacune de ces mailles les valeurs d’accumulation de flux dérivées du MD8 ont été calculées dans WhiteBox GAT. Pour finir, les 4551 mailles ont été fusionnées et les valeurs maximales d’accumulation de flux dérivées du MD8 ont été conservées dans les zones de recouvrement. Il est donc probable que les valeurs d’accumulation de flux dérivées du MD8 soient sous-estimées car on ne travaille pas à l’échelle du sous-bassin hydrographique. Cependant, puisque les zones d’incertitude sont définies à partir d’une valeur d’accumulation de flux correspondant à 0.5 ha on peut raisonnablement penser qu’elles sont dans l’ensemble définies de manière correcte.

La démarche méthodologique a exploité les outils suivants :
- ArcGIS 10.2;
- Script python (v. 2.7);
- WhiteBox GAT (v. 3.2.2).

Références : Maugnard, A., Bielders, C., Defourny, P. (Earth and Life Institute - UCL). Mise à jour du produit cartographique "axes de ruissellement concentré sur l'ensemble de la région wallonne sur base du nouveau jeu de données LiDAR et autres données disponibles". Rapport final - Août 2016.