Análisis multicriterio de riesgo de erosión con ArcGIS

En este apartado vamos a hacer una análisis de riesgo de erosión en un terreno. Es un análisis multicriterio, vamos a analizar la pendiente y la proximidad a los ríos.

Podemos dividir el proceso en 3 partes:

Analizar las pendientes del terreno
Analizar la proximidad a los ríos
Realización del análisis multicriterio

1. ANALIZAR LAS PENDIENTES DEL TERRENO

Debemos contar con una capa raster MDT.

Es el paso más sencillo, en primer lugar, usaremos la herramienta de Slope. Dejaremos la opción Degree.

Spatial Analyst Tools > Surface > Slope.

Ahora el problema es que tenemos los datos en una escala que no nos interesa.

Por lo tanto, vamos a reclasificarlo en 4 grados según su inclinación, siendo el número más alto el de mayor inclinación.

Spatial Analyst Tools > Reclass > Reclassify

Hacemos clic en Classify primero para fijar la cantidad de datos por el que queremos partirlo, y podemos definir el valor que queremos meter en cada uno.

Luego hacemos clic en OK, volveremos a la pantalla anterior y darmos los lugares con mayor inclinación un número más alto.

2. ANALIZAR LA PROXIMIDAD A LOS RÍOS

Para ello tendremos que sacar las cuencas de los ríos, seguiremos el siguiente proceso.

Spatial Analyst Tools > Hidrology > Flow accumulation

A continuación, con el raster obtenido iremos a:

TRUCO: Volviéndonos a meter en reaclasiffy y una vez dentro de la herramienta en classify, podemos fijarnos en la gráfica y encontramos el valor de acumulación que podemos decir que es el río, por lo tanto ahora usaremos la calculadora raster y diremos que los valores mayores que ese valor los seleccione.

Spatial Analyst Tools > Hydrology > Stream to Order. Pondremos en la primera ventana el raster que acabamos de conseguir, y en la segunda el resultante de «flow direction«.

Esto lo pasaremos a polilínea y de aquí haremos un buffer de 300m y de ello volveremos a raster.

Geoprocessing > Buffer, meteremos la capa de líneas que acabamos de obtener y seleccionamos 300 metros (sacará el área alrededor de 300m de toda la polilínea), haremos clic en disolver todo.

Finalmente, haremos un ráster.

3. REALIZACIÓN DEL ANÁLISIS MULTICRITERIO

Ahora reclasificaremos la el raster del buffer de forma que tenga valor 4 (ya que todo el raster tendrá valor 0, daremos a este un valor alto para la futura suma).

Y haremos una calculadora raster donde sumaremos la pendiente reclasificada y el buffer reclasificado.

Las zonas que sean 8 serán las de mayor peligro, e igual conforme mas va bajando habrá menos riesgo.

Superponiéndolo al sombreado del MDT y solo marcando en rojo las zonas de mayor riesgo de erosión, podremos apreciar más claramente cuáles son las zonas más susceptibles a la erosión.