Las herramientas de Spatial Analyst e Hydrology se pueden usar para identificar y extraer líneas de división (de montañas) desde un raster de elevación (DEM), esto al crear una red de flujos invirtiendo los valores de un DEM.
Si en algún momento tienes el interés de buscar una forma de identificar las crestas o filos de una montaña, a continuación se muestra el procedimiento. En esta entrada se usa ArcMap, pero siguiendo el mismo fundamento se puede realizar en otras aplicaciones como ArcGIS Pro o QGIS.
Antes de iniciar asegurarse corregir el DEM con la herramienta Fill ubicada en:
ArcToolbox > Spatial Analyst > Hydrology >
Ahora con ayuda de la herramienta Raster Calculator (Algebra de mapas) invertir los valores del DEM. Es un proceso extremadamente sencillo, simplemente multiplicar por -1. La herramienta Raster Calculator se la puede localizar en:
ArcToolbox > Spatial Analyst > Map Algebra >
Con la herramienta Flow Direction (dirección de flujo) usar el DEM invertido como datos de entrada, activar la casilla «Force all edge cells to flow outward». La herramienta de dirección de flujo se encuentra en:
ArcToolbox > Spatial Analyst > Hydrology >
Usando la misma dirección de ArcToolbox se ejecuta la herramienta Flow Accumulation utilizando como datos de entrada el ráster de acumulación generado anteriormente. Asegúrese de establecer el tipo de datos de salida como Integer (Entero), en lugar de Float (Flotante).
Configure el umbral del ráster creado con Flow Accumulation y cree la red de flujo con Raster Calculator usando la siguiente expresión:
Con(«FlowAcc»,1,0, «Value >= 2500»)
Nota: el valor de 2500 se lo debe reemplazar acorde a la extensión espacial del área de estudio, y la densidad de la red de crestas a identificar. Por ejemplo, es este caso si deseo solamente las crestas principales aumentaría el valor a 10000, pero si deseo una red muy densa bajo el valor a 500. En cada caso se recomienda probar diversos valores, hasta obtener el resultado deseado.
Si todo marcha bien, debes obtener una imagen similar a la siguiente:
Ahora para transformar la imagen ráster en vector es necesario crear los órdenes de flujo. Para ello usar la herramienta Stream Order, usar como datos de entrada el ráster de crestas y el de dirección, esta herramienta se encuentra ubicada en:
ArcToolbox > Spatial Analyst > Hydrology >
Para transformar en vector usar la herramienta Stream To Feature localizada en:
ArcToolbox > Spatial Analyst > Hydrology >
Finalmente, si todo ha salido bien se obtendrá una imagen como la que se muestra seguidamente. También se debe considerar que este proceso usa el mismo principio para calcular una red hídrica, pero en este caso es invertida y con dificultad se conectará toda la red.
Buenas tardes
Tengo un pequeño inconveniente, cuando coloco la expresión en la calculadora raster me emite un mensaje de error el cual dice lo siguiente python syntax error: parsing error syntaxerror: invalid syntax (line 1)
Por ende no puedo terminar el proceso, tu blog me ha sido de gran ayuda, excelente trabajo!
saludos desde Santa Marta – Colombia