Introducción

Este es un cuaderno R Markdown (http://rmarkdown.rstudio.com) en el que se muestra el cálculo de pendiente y orientación.

Creación de Raster

rm(list = ls())
library(raster)
## Loading required package: sp

Creemos un DEM de juguete. Tenga en cuenta que es necesario pasar las coordenadas del cuadro delimitador (mínimo y máximo a lo largo los dos ejes, x e y):

dem <- raster(ncol=3, nrow=3, xmn=100, xmx=115, ymn=100, ymx=115)

¿Cuántas celdas componen DEM?

ncell(dem)
## [1] 9

¿Cuál es la resolución espacial(tamaño de celda) de DEM?

res(dem)
## [1] 5 5

Para asignar valores de elevación a DEM escribimos:

valores <- c(50, 45, 50, 30, 30, 30, 8, 10, 10)
(values(dem) <- valores)
## [1] 50 45 50 30 30 30  8 10 10

Tracemos el DEM a lo largo de los valores de elevación:

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

Para asignar un sistema de referencia de coordenadas a DEM usamos:

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

Cálculo de pendiente y orientación

(slope =  terrain(dem, 'slope', unit = 'degrees', neighbors = 8))
## class      : RasterLayer 
## dimensions : 3, 3, 9  (nrow, ncol, ncell)
## resolution : 5, 5  (x, y)
## extent     : 100, 115, 100, 115  (xmin, xmax, ymin, ymax)
## crs        : +init=epsg:3115 
## source     : memory
## names      : slope 
## values     : 75.25766, 75.25766  (min, max)
plot(slope, main = "Pendiente")
text(slope)

(aspecto = terrain(dem, 'aspect', unit = 'degrees' , neighbors = 8))
## class      : RasterLayer 
## dimensions : 3, 3, 9  (nrow, ncol, ncell)
## resolution : 5, 5  (x, y)
## extent     : 100, 115, 100, 115  (xmin, xmax, ymin, ymax)
## crs        : +init=epsg:3115 
## source     : memory
## names      : aspect 
## values     : 180.7538, 180.7538  (min, max)
plot(aspecto, main = "Aspecto")
text(aspecto)

Asignación

Escriba el código R para calcular la pendiente y el aspecto para el siguiente DEM: Para este nuevo juguete modificamos el nombre del objeto “dem” a “ndem” solo para distinguirle del primer juguete creado.

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

Como antes, vamos a insertar el comando que nos informa de cuántas celdas se compone el DEM:

ncell(ndem)
## [1] 16

¿Cuál es la resolución espacial (tamaño de celda) de DEM?

res(ndem)
## [1] 5 5

Asignar valores de elevación a DEM:

valoresn <- c(50, 45, 50, 48, 30, 29, 30, 29, 10, 9, 9, 10, 25, 23, 19, 21)
(values(ndem) <- valoresn)
##  [1] 50 45 50 48 30 29 30 29 10  9  9 10 25 23 19 21

Tracemos el DEM a lo largo de los valores de elevación:

plot(ndem, main = "Cálculo de pendiente y aspecto para este DEM")
text(ndem)

Asignemos un sistema de referencia de coordenadas a DEM:

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

Cálculo de pendiente y orientación para DEM

(nslope =  terrain(ndem, '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)
plot(nslope, main = "Pendiente")
text(nslope)

(naspecto = terrain(ndem, '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(naspecto, main = "Aspecto")
text(naspecto)