##Angelica Martinez##

PENDIENTE Y ASPECTO

La pendiente o gradiente de un plano es un número que describe la dirección y lainclinación de este.

El aspecto se refiere a la orientación de la pendiente del plano en relación con los puntos cardinales.

En este cuaderno aprenderemos a crear un DEM y tambien a calcular la pendiente y orientación de un plano.

rm(list=ls())

library(raster)

Vamos a crear un DEM de juguete. Hay que tener en cuenta que es necesario pasar las coordenadasdel cuadro delimitador (mínimo y máximo a lo largo de los ejes x e y)

dem <- raster(ncol=4, nrow=4, xmn=100, xmx=120, ymn=100, ymx=120)

Con la funciónn "ncell podemos ver el número de celdas que tiene nuetro DEM

ncell(dem)

[1] 16

Con la función “res()” podemos conocer la resolución espacial, es decir el tamaño de las celdas del DEM

res(dem)

[1] 5 5

Asignaremos valores a cada celda que compone el DEM, para esto se debe teneren cuenta el orden, se empieza desde la superior izquierda se avanza hacia la derecha y se va desciende a la otra fila siempre empezando por la izquierda.

valores <- c(50, 45, 50, 48, 30, 29,30, 29, 10,9,9,10,25,23,19,21)

(values(dem) <- valores)

 [1] 50 45 50 48 30 29 30 29 10  9  9 10 25 23 19 21

Tracemos el DEM con los valores de elevación asignados anteriormente.

plot(dem, main = "DEM")
text(dem)

Asignaremos al DEM un sistema de coodenadas.

crs(dem) <- CRS('+init=epsg:3115')

(slope = terrain(dem, 'slope', unit='degrees', neighbors=8))

class      : RasterLayer 
dimensions : 4, 4, 16  (nrow, ncol, ncell)
resolution : 5, 5  (x, y)
extent     : 100, 120, 100, 120  (xmin, xmax, ymin, ymax)
crs        : +init=epsg:3115 
source     : memory
names      : slope 
values     : 36.05503, 75.62313  (min, max)

Cálculo de la pendiente y el aspecto

plot(slope, main = "Pendiente")
text(slope)

(aspecto = terrain(dem, 'aspect', unit='degrees', neighbors=8))

class      : RasterLayer 
dimensions : 4, 4, 16  (nrow, ncol, ncell)
resolution : 5, 5  (x, y)
extent     : 100, 120, 100, 120  (xmin, xmax, ymin, ymax)
crs        : +init=epsg:3115 
source     : memory
names      : aspect 
values     : 164.0546, 181.4688  (min, max)

plot(aspecto, main = "Aspecto")
text(aspecto)