Mostrar una manera de trabajar objetos espaciales del tipo vectorial (Shape) y raster (tiff) en R para realizar estudios agroambientales.
Si conocen otra forma de hacerlo seria un gusto para nosotros recibir su valioso aporte y, asà colaborar con nuestra comunidad.
-Emplearemos aplicaciones como R, RStudio y QGis.
-Datos de la cuenca hidrográfica del rÃo Canoabo en Venezuela, tales como partes de agua, modelos digital de elevación, sitios de muestreos de suelos entre otros.
-Los procedimientos que llevaremos a cabo incluyen: Creaciones de objetos R (vectores y raster), ploteo de capas, proyección cartográfica, entre otros que a continuación se describen.
Verificación de directorio de trabajo
getwd()
## [1] "C:/ecuador"
Activamos las librerias o paquetes necesarios:
library(raster)
## Loading required package: sp
Trabajamos con capas espaciales de SIG:
#Creamos un objeto R de un archivo shapefile shp
perimetro <- shapefile("C:/Sevilla/Canoabo1/cuenca1.shp")
#Creamos un objeto R de un archivo formato tiff raster
DEM <- raster("C:/Sevilla/Canoabo1/MDE3.tif")
#Ploteamos juntos los objetos R recien creados (Modelo digital de elevacion y el perimetro)
plot(DEM, main="Modelo Digital de Elevacion (msnm)")
lines(perimetro)
# Si desea ver los objetos R en memoria
ls()
## [1] "DEM" "perimetro"