Determinación del sub factor L, longitud del terreno (USLE-RUSLE), a partir de un DEM con GRASS-QGIS

En el artículo anterior se señaló un procedimiento para la determinación del sub factor S, (topografía del terreno USLE-RUSLE), a partir de un DEM con GRASS-QGIS. En éste, para complementar la determinación del factor topográfico (LS), se plantea el cálculo del subfactor longitud del terreno L. Sin embargo, a diferencia del subfactor pendiente del terreno (S), su estimación puede hacerse empleando varios enfoques y es, en consecuencia, algo más complicado establecer en algunos casos la variabilidad espacial del mismo. Por ello, en esta primera entrega se establecerá el cálculo empleando el caso más sencillo, es decir, asumiendo una longitud del terreno (λ) constante: 50 y 100 m; respectivamente.

Según Barrios (2000), el cálculo individual de L se basa en las siguientes expresiones:

    L = (λ/22.13)m
    m = β/(1 + β)
    β = {(sen θ/0.0896)/[3.0 (sen θ)0.8 + 0.56]}*r

donde:
L subfactor longitud del terreno (adimensional)
λ longitud del terreno (en metros)
m exponente variable según β
β relación erosión en surco a erosión en entre surco
θ ángulo de inclinación del terreno
r coeficiente igual a: 0.5, en tierras forestales o pastizales; 1.0, en terrenos agrícolas; y 2.0, en sitios de construcción (en este caso, como no se dispone de mapa de uso de la tierra se asume el valor de 1 para r).

Por tanto, a partir del mapa de pendientes, mediante álgebra de mapas, sería posible obtener la distribución espacial de L mediante la siguiente serie de comandos ejecutados en cónsola de GRASS-QGIS o de bash.

r.mapcalc beta="(sin(pendiente_deg)/0.0896)/(3*(sin(pendiente_deg))^0.8+0.56)"
r.mapcalc m="beta/(1+beta)"
r.mapcalc L50="(50/22.13)^m"
r.mapcalc L100="(100/22.13)^m"

Este es el aspecto del ráster L100:

A continuación, mediante toda la serie de mapas producidos se verifica, mediante una hoja de cálculo, que el álgebra de mapas se realizó de la manera adecuada; tal como se presenta en la siguiente imagen:

About these ads
Esta entrada fue publicada en GRASS, QGIS, SIG, Software Libre, Teledetección. Guarda el enlace permanente.

4 respuestas a Determinación del sub factor L, longitud del terreno (USLE-RUSLE), a partir de un DEM con GRASS-QGIS

  1. Pingback: Determinación del sub factor L, longitud del terreno (USLE-RUSLE), a partir de un DEM con GRASS-QGIS | #Geoprocessamento em Foco | Scoop.it

  2. Pingback: Estimación del factor topográfico LS (RUSLE) en GRASS-QGIS según la clase de pendiente |

  3. Pingback: Cálculo de la longitud del terreno (λ de la RUSLE) según enfoque CALSITE mediante r.drain en GRASS-QGIS |

  4. Pingback: Cálculo de la longitud del terreno (λ de la RUSLE) según enfoque CALSITE mediante un script de Python con GRASS |

Deja un comentario

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s