data-raw/data_raw.R

In MariekeDirk/solarEuR: Interpolation and time serie analysis of solar radation in Europe

library(data.table)
library(raster)
main_dir<-"/nobackup/users/dirksen/data/radiation_europe/"

#Ensembles data and metadata 
ensembles_meta<-fread(paste0(main_dir,"stations/ensembles_qq_stations"),nrows = 102)
ensembles_meta<-ensembles_meta[,2:9]
names(ensembles_meta)<-c("STATION","STATION NAME","COUNTRY NAME","LAT","LON","ELEV","START","STOP")

devtools::use_data(ensembles_meta,overwrite = TRUE)

ensembles<-fread(paste0(main_dir,"stations/ensembles_qq"))
ensembles<-subset(ensembles,select=which(names(ensembles)!="|"))

devtools::use_data(ensembles,overwrite = TRUE)

#Synops data and metadat
synops_meta<-fread(paste0(main_dir,"stations/synops_qq_stations"),nrow=598)
synops_meta<-synops_meta[,2:9]
names(synops_meta)<-c("STATION","STATION NAME","COUNTRY NAME","LAT","LON","ELEV","START","STOP")

devtools::use_data(synops_meta,overwrite = TRUE)

synops<-fread(paste0(main_dir,"stations/synops_qq"))
synops<-subset(synops,select=which(names(synops)!="|"))

devtools::use_data(synops,overwrite = TRUE)

#solar radiation data
sarah.example<-list.files("/nobackup/users/dirksen/data/SARAH_raster/",pattern=".grd",full.names = TRUE)[5000]
sarah.grid<-stack(sarah.example)
names(sarah.grid)<-"Irr"

devtools::use_data(sarah.grid,overwrite = TRUE)

#natural earth dataset
WGS84<-CRS("+init=epsg:4326")
mymap.unpro=readOGR(dsn='/nobackup/users/dirksen/data/NaturalEarthData/ne_10m_admin_0_countries',layer="ne_10m_admin_0_countries") # Read in (unprojected) map data
countries=spTransform(mymap.unpro, WGS84) # Reproject the map

devtools::use_data(countries,overwrite = TRUE)

#DEM and land/sea mask
mask_topo_1km<-stack("/nobackup/users/dirksen/data/radiation_europe/DEM/gtopo30_gis_1km.grd")
devtools::use_data(mask_topo_1km,overwrite = TRUE)

#grid used for the current E-OBS
mask2_topo_1km<-stack("/nobackup/users/dirksen/data/radiation_europe/DEM/gmted2010_gis_1km.grd")
devtools::use_data(mask2_topo_1km,overwrite = TRUE)

#NASA global radiation
global_radiation<-fread("https://eosweb.larc.nasa.gov/sse/global/text/global_radiation")
coordinates(global_radiation)<-~Lon+Lat
crs(global_radiation)<-WGS84
gridded(global_radiation)<-TRUE
global_radiation<-stack(global_radiation)
devtools::use_data(global_radiation,overwrite=TRUE)