R Notebook

Errika Aparicio

01/04/2020

1. Introducción

Este es un cuaderno de R Markdown para ilustrar el cálculo de la pendiente y el aspecto.

En matemáticas, la pendiente o gradiente de un plano es un número que describe tanto la dirección como la inclinación de un plano

El aspecto, porsu parte, se refiere a la orientación de la pendiente del plano en relación con los puntos cardinales. En breve, aspecto es la dirección de la brújula que enfrenta una pendiente.

2. Creación del raster

rm(list=ls())

library(raster)

Se va a crear un juguete DEM. Es necesario tener en cuenta que se debe 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)

Pero, ¿cuántas celdas componen el DEM?

ncell(dem)

[1] 9

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

res(dem)

[1] 5 5

Ahora, se debe asignar valores de elevación a DEM:

valores <- c(50, 45, 50, 30, 30, 30, 8, 10, 10)

(values(dem) <- valores)

[1] 50 45 50 30 30 30  8 10 10

Posteriormente se Traza el DEM a lo largo de los valores de elevación:

Ahora, se debe asignar un sistema de referencia de coordenadas a DEM:

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

3. Cálculo de la pendiente y el aspecto

(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)

(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)

4. Asignación

Ahora se creara un nuevo DEM, no olvidar que es necesario tener en cuenta que se debe pasar las coordenadas del cuadro delimitador (mínimo y máximo a lo largo los dos ejes, x e y):

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

Pero, ¿cuántas celdas componen este DEM?

ncell(dem)

[1] 16

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

res(dem)

[1] 5 5

Ahora, se debe asignar valores de elevación a DEM:

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
[14] 23 19 21

Posteriormente se Traza el DEM a lo largo de los valores de elevación:

Asignación del sistema de referencia de coordenadas a DEM:

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)

(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)