¿Cómo identifica las crestas en un ráster de elevación?

Me preguntaba si alguien podría decirme en términos simples cómo extraer o diferenciar entre crestas y cimas de los datos de elevación de 10 m con ArcGIS 10.1.

La definición de las crestas frente a las cimas de las colinas / montañas depende bastante de la escala. Jeff Jenness cubre conceptualmente cómo modelar formas topográficas topográficas en su artículo Algunos pensamientos sobre el análisis de las características del hábitat topográfico . Si hurga en su sitio web, también puede encontrar su póster sobre esto, en Herramientas de ArcGIS> Diseñador de corredor de facetas terrestres. (El enlace está aquí )

Jeff tiene una herramienta para modelar la posición topográfica en su extensión de herramientas de facetas de tierra para ArcMap, y reescribí una versión simple de 4 clases (crestas, cañones, pendientes planas, pendientes pronunciadas) de su herramienta como parte de una caja de herramientas de Python (descargue el CorridorDesigner general usar caja de herramientas ).

Si está haciendo esto desde cero, puede:

1. Ejecute estadísticas focales (media) en un DEM utilizando un tamaño de vecindario que tenga sentido para su análisis (aunque los vecindarios grandes se ahogarán en un DEM de 10 m).
2. Reste las estadísticas focales medias del DEM original. Esta capa (índice de posición topográfica - TPI) le dará una medida de si un píxel está por encima o por debajo de su promedio de vecindad. Los valores negativos indican que un píxel está por debajo de la elevación promedio de su vecindario (valles y fondos de cañones). Los valores positivos indican que un píxel está por encima de la elevación promedio de su vecindario (cumbres y cumbres).
3. Reclasifique la capa TPI para elegir drenajes, crestas y pendientes. Con mi implementación de la herramienta, utilicé valores> = 12 para indicar crestas, <= para indicar fondos de cañones y entre -12 y 12 para indicar pendientes. Lo que use para estos umbrales realmente depende de su paisaje y estudio.
4. Si está interesado en dividir las clases de pendientes (p. Ej., Áreas planas frente a pendientes pronunciadas), puede calcular opcionalmente la pendiente, reclasificarla en pendientes planas / suaves / inclinadas, y luego usar una instrucción CON para reemplazar los valores de pendiente calculados en el paso 3 con Categorías de pendientes.

La creación de cuencas hidrográficas debería ayudarlo a localizar las crestas y la cima de la colina. Entonces, definiría la cima de una colina como un máximo local, mientras que un punto en una cresta no es el máximo (hay otro punto más alto o igual a este punto). Puede identificar máximos locales utilizando la herramienta de estadística focal.

Otra forma de ver el problema es analizar lo opuesto a su DEM (cero menos DEM). Calcula la acumulación de flujo que convierte en flujos. Uno de los extremos de las corrientes será una parte superior, y la corriente será la cresta.

Tenga en cuenta que 10 m DEM es bastante preciso, por lo que para este tipo de análisis probablemente sea mejor suavizar su capa primero para evitar máximos locales.