1. Creacion de un raster

La siguiente funcion nos permite iniciar con la consola de R en blanco, borrando lo que tenga esta, debido a que la funcion “rm ()” es básicamente ‘remove {base}’, se usa para eliminar objetos de un entorno especificado. “list = ls ()” es la base de este comando, lo que significa que se refiere a todos los objetos presentes en el espacio de trabajo.

rm(list=ls())

Un raster es una estructura o fichero de datos que representa una rejilla rectangular de píxeles o puntos de color, denominada matriz. Para la creacion del raster lo primero que vamos a realizar es la instalacion del paquete estadistico de R, “raster” y posteriormente llamarlo con la funcion “library()”

#install.packages("raster")
library(raster)

Ahora vamos a crear un objeto al que llamaremos “dem”, pondremos las dimensiones del “DEM” utilizando las funciones “xmn”, “xmx”como minimos y maximos del eje x y las funciones “ymn” y “ymx” del eje y. Ademas las funciones “ncol” como el numero de columnas y “nrow” como el numero de filas.

dem <- raster(ncol=4, nrow=4, xmn=100, xmx=120, ymn=100, ymx=120)
ncell(dem) #No. de celdas que compone el dem
## [1] 16
res(dem) #Es la resolucion del raster o el tamaño de la celda
## [1] 5 5
#se asignan los valores a cada celda
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
plot(dem, main="DEM") #La grafica 
text(dem)

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

¿Qué significa el código EPSG ?

EPSG significa European Petroleum Survey Group quienes publican de forma libre una base de datos de informacion del sistema de coordenadas mas algunos documentos relacionados sobre proyecciones de mapas y datos. El conjunto de datos de parametros geodesicos de EPSG es un conjunto de datos estructurado de sistemas de referencia de coordenadas y transformaciones de coordenadas

¿Dónde se pueden encontrar los códigos EPSG ?: El conjunto de datos EPSG primario se mantiene en su registro en linea (http://www.epsg-registry.org/ ), desde el cual se puede acceder a los datos a traves de una interfaz grafica de usuario o una interfaz de servicio. El registro en linea contiene los datos mas actuales. El conjunto de datos de parametros geodesicos de EPSG, ofrecido a traves de las paginas web de IOGP, se puede usar de forma gratuita, pero su uso esta sujeto a la aceptacion de los Terminos de uso .

¿Cuál es el código EPSG para MAGNA Ciudad de Bogotá ?: 3116

2. Calculo de Pendiente y Aspecto

#Pendiente
(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:3116 +proj=tmerc +lat_0=4.596200416666666 +lon_0=-74.07750791666666 +k=1 +x_0=1000000 +y_0=1000000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
## source     : memory
## names      : slope 
## values     : 36.05503, 75.62313  (min, max)
plot(slope, main = "Pendiente")
text(slope)

#Aspecto
(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:3116 +proj=tmerc +lat_0=4.596200416666666 +lon_0=-74.07750791666666 +k=1 +x_0=1000000 +y_0=1000000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
## source     : memory
## names      : aspect 
## values     : 164.0546, 181.4688  (min, max)
plot(aspecto, main = "Aspecto")
text(aspecto)

 #El valor de cada celda indica la direccion a la que apunta la pendiente de la celda, en este caso valores entre 135 y 220 indican direccion norte