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

Cuando se evalúa la erosión del suelo a escala de cuencas hidrográficas, usando las tecnologías USLE o RUSLE, se requiere el cálculo del factor topográfico LS que integra, por separado, los subfactores adimensionales longitud (L) y pendiente (S) del terreno. La distribución espacial del subfactor S puede ser obtenida de un mapa de pendientes que a su vez han sido derivados de un DEM a partir del módulo r.slope.aspect en GRASS. Para la obtención de los ráster de pendientes, respectivamente en grados y porcentaje, se tiene:

r.slope.aspect elevation=N09W068_UTM19N_canoa slope=pendiente_deg format=degrees
r.slope.aspect elevation=N09W068_UTM19N_canoa slope=pendiente_per format=percent

donde se utilizó un DEM proyectado en UTM (zona 19N). Los resultados fueron los siguientes:

En varios puntos al azar de ambos rásteres se verificó que la relación matemática era la esperada:

mapa_porcentaje = tan(mapa_grados)*100

Como las expresiones funcionales para la estimación de S son las siguientes:

    S = 10,8.Sen θ + 0,03; p < 9%
    S = 16,8.Sen θ – 0,50; p >= 9%

entonces, mediante álgebra de mapas con el módulo r.mapcalc, se obtiene el ráster subf_S (subfactor pendiente del terreno S):

r.mapcalc subf_S = "if(pendiente_per < 9, 10.8*sin(pendiente_deg)+0.03, 16.8*sin(pendiente_deg)-0.5)"

que se visualiza tal como se presenta en la siguiente imagen:

No obstante, para verificar si el álgebra de mapas se realizó de la manera adecuada para obtener la distribución espacial del subfactor pendiente del terreno, se hizo zoom en una zona dentro del área de trabajo y con la ayuda de los mapas de pendiente, el ráster subf_S y una hoja de cálculo, se verificó para dos celdas al azar con porcentaje de pendiente menor y superior a 9 que r.mapcalc funcionó de la manera esperada; tal como se encuentra en la siguiente imagen:

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

3 respuestas a Determinación del sub factor S, pendiente del terreno (USLE-RUSLE), a partir de un DEM con GRASS-QGIS

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

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

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

Deja un comentario

Introduce tus datos o haz clic en un icono para iniciar sesión:

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