La définition des crêtes par rapport aux sommets des collines / montagnes dépend assez de l'échelle. Jeff Jenness explique de façon conceptuelle comment modéliser les reliefs topographiques dans son article Quelques réflexions sur l'analyse des caractéristiques de l'habitat topographique . Si vous fouillez son site Web, vous pouvez également trouver son affiche à ce sujet, sous Outils ArcGIS> Concepteur de corridor de facette du sol. (Le lien est ici )
Jeff a un outil pour modéliser la position topographique dans son extension d'outils de facettes terrestres pour ArcMap, et j'ai réécrit une version simple de classe 4 (crêtes, canyons, pentes plates, pentes abruptes) de son outil dans le cadre d'une boîte à outils Python (téléchargez le CorridorDesigner général utiliser la boîte à outils ).
Si vous faites cela à partir de zéro, vous pouvez:
- Exécutez des statistiques focales (moyennes) sur un DEM en utilisant une taille de quartier qui a du sens pour votre analyse (les grands quartiers s'étoufferont cependant avec un DEM de 10 m).
- Soustrayez la moyenne des statistiques focales du DEM d'origine. Cette couche (indice de position topographique - TPI) vous donnera une mesure si un pixel est supérieur ou inférieur à sa moyenne de voisinage. Les valeurs négatives indiquent qu'un pixel est en dessous de l'élévation moyenne de son voisinage (vallées et fonds de canyon). Les valeurs positives indiquent qu'un pixel est au-dessus de l'élévation moyenne de son voisinage (sommets des crêtes et sommets des collines).
- Reclassifiez la couche TPI pour sélectionner les drainages, les crêtes et les pentes. Avec mon implémentation de l'outil, j'ai utilisé des valeurs> = 12 pour indiquer les crêtes, <= pour indiquer les fonds de canyon et entre -12 et 12 pour indiquer les pentes. Ce que vous utilisez pour ces seuils dépend vraiment de votre paysage et de votre étude.
- Si vous souhaitez séparer les classes de pente (par exemple, zones plates vs pentes raides), vous pouvez éventuellement calculer la pente, la reclasser en pentes plates / douces / raides, puis utiliser une instruction CON pour remplacer les valeurs de pente calculées à l'étape 3 par catégories de pente.
La création de bassins versants devrait vous aider à localiser les crêtes et le sommet des collines. Ensuite, je définirais un sommet de colline comme un maximum local, tandis qu'un point sur une crête n'est pas le maximum (il y a un autre point plus haut ou égal à ce point). Vous pouvez identifier les maxima locaux à l'aide de l'outil de statistique focale.
une autre façon de voir le problème est d'analyser à l'opposé de votre DEM (zéro moins DEM). Vous calculez l'accumulation de flux que vous convertissez en flux. L'une des extrémités des ruisseaux sera un sommet et le ruisseau sera la crête.
Notez que 10 m DEM est assez précis, donc pour ce type d'analyse, il est probablement préférable de lisser votre couche en premier afin d'éviter les maxima locaux.
la source