-
Notifications
You must be signed in to change notification settings - Fork 0
Test extract_data() function
Alexandru Dumitrescu edited this page Mar 21, 2019
·
9 revisions
As an example, we will use air temperature data available (E-OBS, MERA, EUROCORDEX) from package interdecisdata; this package can be installed with:
install.packages("http://193.26.129.60/anm_maps/interdecisdata_0.0.001.tar.gz", repos = NULL, type = "source")
Load EOBS air temperature dataset:
tavg.eobs <- raster::stack(system.file("eobs/tg_ens_mean_0.1deg_reg_2004-2005_v19.0e.nc", package = "interdecisdata"))
Load MERA 2m air temperature dataset:
tavg.mera <- raster::stack(system.file("mera/mera_t2m_daily_2004_2005.nc", package = "interdecisdata"))
raster::plot(tavg.mera[[1:2]])
Load EUROCORDEX tas dataset:
tavg.cordex<- raster::stack(system.file("eurocordex/tas_EUR-11_ECMWF-ERAINT_evaluation_domeniu.mic_r1i1p1_SMHI-RCA4_v1_day_20040101-20051231.nc", package = "interdecisdata"))
raster::plot(tavg.cordex[[1:2]])
Extract reanalysis data to point locations (lonlat coordinates):
library(interdecis)
point <- cbind(lon = -8, lat = 52)
tavg.eobs.co <- extract_data(tavg.eobs, point, date1 = "2004-01-01", date2 = "2005-12-31")
tavg.mera.co <- extract_data(tavg.mera, point, date1 = "2004-01-01", date2 = "2005-12-31")
tavg.cordex.co <- extract_data(tavg.cordex, point, date1 = "2004-01-01", date2 = "2005-12-31")