Extraer datos de Chelsa

Instalar y cargar paquetes

#install.packages(pkgs=c("sp", "rgdal", "raster", "RStoolbox"))
library(sp)
library(rgdal)
library(raster)
library(RStoolbox)
library(tidyverse)

Hacer una lista de los datos raster

list_dirs <- list.files(path = "DATA/raster/", pattern = glob2rx("*tif$"), all.files=T, full.names=T)
names(list_dirs)<- basename(list_dirs)

Apilar todos los raster en un solo objeto

timeStart<- proc.time()
chelsa_data<- stack(list_dirs)
proc.time() - timeStart

   user  system elapsed 
   0.85    0.06    1.00 

Leer y crear el archivo de coordenadas

pointCoordinates<- read.csv(file = "DATA/coord.csv", sep = ";", dec = ".")
coordinates(pointCoordinates)= ~ Long+ Lat

Extraer los valores del la imagen multibanda a cada uno de los puntos seleccionados

timeStart<- proc.time()
rasValue<- raster::extract(chelsa_data, pointCoordinates)
proc.time() - timeStart

   user  system elapsed 
   4.86    0.00    4.86 

Asignarle los valores a cada localizacion y crear el data frame

combinePointValue<- cbind(pointCoordinates,rasValue)
datos<- as.data.frame(combinePointValue)
datos$ID<- as.factor(datos$ID)

Exportar

write.table(datos, file = "data_Miquel.txt", append = FALSE, sep = ",", row.names = FALSE, col.names = TRUE)