-
Notifications
You must be signed in to change notification settings - Fork 1
/
README.Rmd
143 lines (100 loc) · 5.57 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
---
output: github_document
---
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r, include = FALSE}
library(knitr)
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-",
cache = TRUE
)
```
# `easyclimate` <a href='https://verughub.github.io/easyclimate/'><img src="man/figures/easyclimate_logo_small.png" align="right" height="139"/></a>
# Easy access to high-resolution daily climate data for Europe
<!-- badges: start -->
`r badger::badge_cran_release()`
<a href="https://diffify.com/R/easyclimate" target="_blank"><img src="https://diffify.com/diffify-badge.svg" alt="The diffify page for the R package easyclimate" style="width: 100px; max-width: 100%;"></a>
[](https://www.repostatus.org/#active) [](https://github.com/VeruGHub/easyclimate/actions/workflows/R-CMD-check.yaml)
[](https://hits.dwyl.com/VeruGHub/easyclimate)
[](https://hits.dwyl.com/VeruGHub/easyclimate)
[](https://cran.r-project.org/package=easyclimate)
<!-- badges: end -->
Get high-resolution (1 km) daily climate data (precipitation, minimum and maximum temperatures) for Europe from the European climatic database hosted at [University of Natural Resources and Life Sciences, Vienna, Austria](https://boku.ac.at/en/wabo/waldbau/wir-ueber-uns/daten). Data are currently available from 1950 to 2022.
This climatic dataset was originally built by [A. Moreno & H. Hasenauer](https://doi.org/10.1002/joc.4436) and further developed by W. Rammer, C. Pucher & M. Neumann (see [this document](https://github.com/VeruGHub/easyclimate/blob/master/inst/Description_Evaluation_Validation_Downscaled_Climate_Data_v2.pdf) for more details on the development and characteristics of the climatic dataset, and [this document](https://doi.org/10.6084/m9.figshare.22962671.v1) for the updates of the last version - v4).
In this R package we implemented [Cloud-Optimised Geotiffs](https://www.cogeo.org/) so that we can obtain daily climate data for thousands of sites/days within seconds/minutes, without having to download huge rasters. But if you need to obtain data for large areas, please download the rasters from the FTP server (<ftp://palantir.boku.ac.at/Public/ClimateData/>) and extract the values locally rather than using this package, so as not to saturate the file server. For that, you may use a FTP client such as [FileZilla](https://filezilla-project.org/).
For a detailed description of {easyclimate}, please read [this paper](https://doi.org/10.1016/j.envsoft.2023.105627) (open access version [here](https://doi.org/10.32942/osf.io/mc8uj)) or visit the package [website](https://verughub.github.io/easyclimate/).
## Installation
Stable version from CRAN:
``` r
install.packages("easyclimate")
```
Development version from GitHub:
``` r
# install.packages("remotes")
remotes::install_github("VeruGHub/easyclimate")
```
## Examples
### Obtain a data frame of climatic values
To obtain a data frame of daily climatic values for point coordinates:
```{r example, message = FALSE}
library(easyclimate)
coords <- data.frame(lon = -5.36, lat = 37.40)
prec <- get_daily_climate(coords,
period = "2001-01-01:2001-01-03",
climatic_var = "Prcp",
version = 4) # default
```
```{r echo=FALSE}
kable(prec)
```
<br>
### Obtain a raster of climatic values
To obtain a (multi-layer) raster of daily climatic values for an area:
```{r message=FALSE, fig.width = 8, fig.height = 3}
library(terra)
## Download the polygon contour of a region
sobrarbe <- mapSpain::esp_get_comarca(comarca = "Sobrarbe")
## Coordinates must be in lonlat
sobrarbe <- project(vect(sobrarbe), "EPSG:4326")
## Download Tmax values for that region between 1st and 3rd May 2020
sobrarbetemp <- get_daily_climate(
coords = sobrarbe,
climatic_var = "Tmax",
period = "2020-05-01:2020-05-03",
output = "raster"
)
```
The output (`sobrarbetemp`) is a SpatRaster with 3 layers (for each of 3 days):
```{r}
sobrarbetemp
```
Let's make a map. First using terra:
```{r map_terra, eval=FALSE}
plot(sobrarbetemp, col = rev(RColorBrewer::brewer.pal(9, "RdYlBu")),
smooth = TRUE, nc = 3)
```
Now using ggplot2 and tidyterra:
```{r map_ggplot, message=FALSE, fig.width = 8, fig.height = 3}
library(ggplot2)
library(tidyterra)
ggplot() +
geom_spatraster(data = sobrarbetemp) +
facet_wrap(~lyr, ncol = 3) +
scale_fill_distiller(palette = "RdYlBu", na.value = "transparent") +
geom_spatvector(data = sobrarbe, fill = NA) +
labs(fill = "Maximum\ntemperature (ºC)") +
scale_x_continuous(breaks = c(-0.25, 0, 0.25)) +
scale_y_continuous(breaks = seq(42.2, 42.8, by = 0.2)) +
theme_minimal()
```
<br> Visit the articles of the [package website](https://verughub.github.io/easyclimate/) for more extended tutorials!
<br>
## CITATION
If you use easyclimate, please cite both the appropriate data source and the package as:
```{r echo=FALSE, results='asis', cache = FALSE}
print(citation("easyclimate"), style = "text")
```