Cálculo de pendiente.

Sandra Katteryne Rodríguez Hurtado

Objetivo

El objetivo de este trabajo es comprende cómo crear un DEM y calcular la pendiente y orientación.

rm(list=ls())

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

ncell(dem) 

[1] 9

rr res(dem)

[1] 5 5

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

rr (values(dem) <- valores)

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

rr plot(dem, main = ) text(dem)

rr crs(dem) <- CRS(‘+init=epsg:3115’)

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

rr plot(slope, main = ) text(slope)

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

rr plot(aspecto, main = ) text(aspecto)

EJERCICIO

Calcular la pendiente y el aspecto para el siguiente DEM, usando el código R.

#[1] 16
#[1] 5 5

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

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

Ahora, se va a mirar el numero de celdas que tiene el DEM, para esto se hace uso de la funcion “ncell()”:

ncell(dem)

[1] 16

Luego, mirar la resolución espacial (es decir, el tamaño de celda) de DEM, para esto se usa “res()”:

res(dem)

[1] 5 5

Ahora, se va a asignar valore a cada una de las celdas del DEM, para ello se debe tener presente el orden.

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

Ahora, se trazará el DEM con los valores de elevación anteriormente asignados.

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

Asignar un sistema de coordenadas al DEM con el que se está trabajando.

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 pendiente y 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)

Cuestionario

1.¿Qué significa el código EPSG?

Por sus siglas en inglés, European Petroleum Survey Group, es una organización que posee una base de datos en parámetros geodésicos con códigos estándar, códigos EPSG, para datos, coordenadas, entre otros, que permite definir datos de posición en los Sistemas de Información geográfica. Juega un papel importante para realizar actividades que requieran gestionar o manipular datos en ambientes digitales. Este grupo está conformado por topógrafos, geógrafos, y cartógrafos que han recolectado información sobre sistemas (marcos) de referencia y modernos (geocéntricos), proyecciones cartográficas y proyecciones de elipsoides.

3.¿Cuál es el código EPSG para MAGNA Ciudad de Bogotá? Se utiliza MAGNA-SIRGAS, es el Marco Geocéntrico Nacional de Referencia, Densificación del Sistema de Referencia Geocéntrico para las Américas. Este sistema fue propuesto para su uso en el país en reemplazo de DATUM-Bgototá.

REFERENCIAS

• https://mappinggis.com/2016/04/los-codigos-epsg-srid-vinculacion-postgis/
• https://www.igac.gov.co/es/contenido/areas-estrategicas/magna-sirgas