#remotes::install_github("kwb-r/kwb.dwd@dev")
#remotes::install_github("kwb-r/kwb.impetus")
library(kwb.impetus)
shape_obj <- kwb.dwd:::get_shape_of_german_region(name = "berlin")
#> Loading 'shapes_germany' from 'C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/cache/shapes_germany.RData'... ok.
shape_file <- file.path(getwd(), "berlin.shp")
shape_obj %>%
sf::st_as_sf() %>%
sf::write_sf(shape_file)
# Works
yearmonth_start <- "190601"
yearmonth_end <- "190712"
dwd_monthly_vars <- c("air temperature (mean)" = "air_temperature_mean"#,
#"drought index" = "drought_index",
#"evaporation, potential" = "evapo_p",
#"evaporation, real" = "evapo_r",
#"precipitation" = "precipitation",
#"soil moisture" = "soil_moist",
#"soil temperature (5 cm)" = "soil_temperature_5cm"
)
system.time(
dwd_berlin_monthly_list1 <- stats::setNames(lapply(dwd_monthly_vars, function(dwd_var) {
kwb.dwd::read_monthly_data_over_shape(
file = shape_file,
variable = dwd_var,
from = yearmonth_start,
to = yearmonth_end,
quiet = TRUE
)
}), nm = dwd_monthly_vars))
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_evapo_p_202201.asc.gz
#> Warning in OGRSpatialRef(dsn, layer, morphFromESRI = morphFromESRI, dumpSRS
#> = dumpSRS, : Discarded datum Deutsches_Hauptdreiecksnetz in Proj4 definition:
#> +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel
#> +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs
#> OGR data source with driver: ESRI Shapefile
#> Source: "C:\Users\mrustl\Documents\RProjects\kwb.impetus\R\berlin.shp", layer: "berlin"
#> with 1 features
#> It has 11 fields
#> ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/: ok.
#> 1/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/01_Jan/: ok.
#> 2/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/02_Feb/: ok.
#> 3/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/03_Mar/: ok.
#> 4/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/04_Apr/: ok.
#> 5/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/05_May/: ok.
#> 6/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/06_Jun/: ok.
#> 7/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/07_Jul/: ok.
#> 8/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/08_Aug/: ok.
#> 9/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/09_Sep/: ok.
#> 10/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/10_Oct/: ok.
#> 11/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/11_Nov/: ok.
#> 12/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/12_Dec/: ok.
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190601/grids_germany_monthly_air_temp_mean_190601.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190701/grids_germany_monthly_air_temp_mean_190701.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190602/grids_germany_monthly_air_temp_mean_190602.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190702/grids_germany_monthly_air_temp_mean_190702.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190603/grids_germany_monthly_air_temp_mean_190603.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190703/grids_germany_monthly_air_temp_mean_190703.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190604/grids_germany_monthly_air_temp_mean_190604.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190704/grids_germany_monthly_air_temp_mean_190704.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190605/grids_germany_monthly_air_temp_mean_190605.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190705/grids_germany_monthly_air_temp_mean_190705.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190606/grids_germany_monthly_air_temp_mean_190606.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190706/grids_germany_monthly_air_temp_mean_190706.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190607/grids_germany_monthly_air_temp_mean_190607.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190707/grids_germany_monthly_air_temp_mean_190707.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190608/grids_germany_monthly_air_temp_mean_190608.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190708/grids_germany_monthly_air_temp_mean_190708.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190609/grids_germany_monthly_air_temp_mean_190609.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190709/grids_germany_monthly_air_temp_mean_190709.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190610/grids_germany_monthly_air_temp_mean_190610.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190710/grids_germany_monthly_air_temp_mean_190710.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190611/grids_germany_monthly_air_temp_mean_190611.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190711/grids_germany_monthly_air_temp_mean_190711.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190612/grids_germany_monthly_air_temp_mean_190612.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190712/grids_germany_monthly_air_temp_mean_190712.asc.gz
#> Masking and cropping grids_germany_monthly_air_temp_mean_190601 ... ok. (0.47s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190701 ... ok. (0.56s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190602 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190702 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190603 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190703 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190604 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190704 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190605 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190705 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190606 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190706 ... ok. (0.42s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190607 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190707 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190608 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190708 ... ok. (0.41s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190609 ... ok. (0.43s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190709 ... ok. (0.43s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190610 ... ok. (0.42s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190710 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190611 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190711 ... ok. (0.40s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190612 ... ok. (0.42s)
#> Masking and cropping grids_germany_monthly_air_temp_mean_190712 ... ok. (0.41s)
#> User System verstrichen
#> 11.27 1.79 29.97
# Does not work
yearmonth_start <- "190701"
yearmonth_end <- "190812"
dwd_monthly_vars <- c("air temperature (mean)" = "air_temperature_mean"#,
#"drought index" = "drought_index",
#"evaporation, potential" = "evapo_p",
#"evaporation, real" = "evapo_r",
#"precipitation" = "precipitation",
#"soil moisture" = "soil_moist",
#"soil temperature (5 cm)" = "soil_temperature_5cm"
)
system.time(
dwd_berlin_monthly_list1 <- stats::setNames(lapply(dwd_monthly_vars, function(dwd_var) {
kwb.dwd::read_monthly_data_over_shape(
file = shape_file,
variable = dwd_var,
from = yearmonth_start,
to = yearmonth_end,
quiet = TRUE
)
}), nm = dwd_monthly_vars))
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_evapo_p_202201.asc.gz
#> Warning in OGRSpatialRef(dsn, layer, morphFromESRI = morphFromESRI, dumpSRS
#> = dumpSRS, : Discarded datum Deutsches_Hauptdreiecksnetz in Proj4 definition:
#> +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel
#> +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs
#> OGR data source with driver: ESRI Shapefile
#> Source: "C:\Users\mrustl\Documents\RProjects\kwb.impetus\R\berlin.shp", layer: "berlin"
#> with 1 features
#> It has 11 fields
#> ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/: ok.
#> 1/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/01_Jan/: ok.
#> 2/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/02_Feb/: ok.
#> 3/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/03_Mar/: ok.
#> 4/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/04_Apr/: ok.
#> 5/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/05_May/: ok.
#> 6/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/06_Jun/: ok.
#> 7/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/07_Jul/: ok.
#> 8/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/08_Aug/: ok.
#> 9/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/09_Sep/: ok.
#> 10/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/10_Oct/: ok.
#> 11/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/11_Nov/: ok.
#> 12/12: ftp://opendata.dwd.de/climate_environment/CDC/grids_germany/monthly/air_temperature_mean/12_Dec/: ok.
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190701/grids_germany_monthly_air_temp_mean_190701.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190801/grids_germany_monthly_air_temp_mean_190801.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190702/grids_germany_monthly_air_temp_mean_190702.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190802/grids_germany_monthly_air_temp_mean_190802.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190703/grids_germany_monthly_air_temp_mean_190703.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190803/grids_germany_monthly_air_temp_mean_190803.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190704/grids_germany_monthly_air_temp_mean_190704.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190804/grids_germany_monthly_air_temp_mean_190804.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190705/grids_germany_monthly_air_temp_mean_190705.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190805/grids_germany_monthly_air_temp_mean_190805.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190706/grids_germany_monthly_air_temp_mean_190706.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190806/grids_germany_monthly_air_temp_mean_190806.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190707/grids_germany_monthly_air_temp_mean_190707.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190807/grids_germany_monthly_air_temp_mean_190807.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190708/grids_germany_monthly_air_temp_mean_190708.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190808/grids_germany_monthly_air_temp_mean_190808.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190709/grids_germany_monthly_air_temp_mean_190709.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190809/grids_germany_monthly_air_temp_mean_190809.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190710/grids_germany_monthly_air_temp_mean_190710.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190810/grids_germany_monthly_air_temp_mean_190810.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190711/grids_germany_monthly_air_temp_mean_190711.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190811/grids_germany_monthly_air_temp_mean_190811.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190712/grids_germany_monthly_air_temp_mean_190712.asc.gz
#> File already available: C:\Users\mrustl\AppData\Local\Temp/R_kwb.dwd/grids_germany_monthly_air_temp_mean_190812/grids_germany_monthly_air_temp_mean_190812.asc.gz
#> Error in .local(.Object, ...) :
#> Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : Cannot create a RasterLayer object from this file.
#> Zeitnahme beendet um: 1.11 0.05 4.11
works for time period
190601-190712but not for190701-190812(see code below)Created on 2022-09-20 by the reprex package (v2.0.1)
Session info