L' aide d'ArcGIS indique:
Les bassins versants peuvent être délimités à partir d'un MNT en calculant la direction du flux et en l'utilisant dans l'outil Bassin versant. Pour déterminer la zone de contribution, un raster représentant la direction du flux doit d'abord être créé avec l'outil Direction du flux.
La direction du flux est calculée à partir du DEM à l'aide de la méthode D8 , où le flux est abstrait en calculant pour chaque cellule, laquelle de ses 8 voisins, l'eau de cette cellule va s'écouler.
Il existe de nombreuses alternatives à D8, telles que Rho8, Froh8 et Stream Tubes, mais la plupart des logiciels SIG, y compris ArcGIS, ont tendance à utiliser D8, car il est plus simple et moins intensif en calcul que les autres.
Il y a quelques années, je travaillais sur un projet de délimitation de bassin versant et nous étions confrontés à plusieurs problèmes dus à l'utilisation d'ArcGIS par la méthode D8. Les deux principaux problèmes étaient
- D8 autorise uniquement l'écoulement unidirectionnel. L'eau ne peut s'écouler que dans une direction à partir d'une cellule.
- Les flux de flux générés avaient un énorme biais le long de l'axe diagonal. Cela a donné naissance à des flux étranges.
D'après nos données, nous savions que ces deux problèmes étaient de gros problèmes, alors j'avais développé des outils pour générer des directions de flux à l'aide de méthodes hybrides.
L'une de mes premières tâches a été de désosser l'outil de calcul du bassin versant. J'ai trouvé que c'était logiquement assez simple. Si vous souhaitez trouver le bassin versant pour un point donné (également appelé point d'écoulement), vous devez d'abord trouver la cellule à laquelle il appartient. Souvent, vous essaierez de l'enclencher au point avec l'accumulation de débit la plus élevée dans une tolérance donnée.
Pour cette cellule, vous trouverez toutes les cellules du quartier qui y contribuent. Pour chacune de ces cellules de voisinage, vous trouvez les cellules qui y contribuent et ainsi de suite. Vous poursuivez ce processus itératif jusqu'à ce que vous ne trouviez aucune nouvelle cellule. C'est alors que vous avez atteint les lignes de crête ou la limite du bassin versant.
J'ai trouvé que mon code simple qui faisait cela pour les rasters ASCII, donnait une sortie presque similaire par rapport à l'outil Watershed d'ArcGIS. Parfois, il y avait une différence de quelques cellules sur la frontière, donc je suis convaincu qu'ArcGIS suit un algorithme D8 non modifié.