In rOpenSpain/climaemet: Climate AEMET Tools

Since the last release, this package has been integrated into rOpenSpain, a community of R enthusiasts whose ultimate goal is to create high-quality R packages for data mining public Spanish open sources.

From version 1.0.0 onward, we have introduced some improvements and (breaking) changes on the package, in order to provide a smoother interaction with the AEMET API service.

API Key

Get your API Key

To be able to download data from AEMET you will need a free API key which you can get at https://opendata.aemet.es/centrodedescargas/obtencionAPIKey

Once that you have your API Key, you can use any of the following methods:

a. Set API Key with `aemet_api_key()`

This is the recommended option. Just type:

aemet_api_key("YOUR_API_KEY", install = TRUE)

Using install = TRUE ensures that the API key is stored on your local computer and it would be reloaded every time you load the library. From now on you can forget about API keys!

b. Use an environment variable

This is a temporary alternative. You can set your API key as an environment variable

Sys.setenv(AEMET_API_KEY = "YOUR_API_KEY")

Note that this is only valid for the current session. You would need to re-run this command each time you restart your session.

c. Modify your `.Renviron` file

This stores your API key permanently on your machine. You can start editing your .Renviron running this command:

usethis::edit_r_environ()

Now you can add the following line to you .Renviron file:

AEMET_API_KEY = YOUR_API_KEY

New features

`tibble` format

From v1.0.0 onward, climaemet provides its results in tibble format. Also, the functions try to guess the correct format of the fields (i.e. something as a Date/Hour now is an hour, numbers are parsed as double, etc.).

See how a tibble is displayed:

# See a tibble in action

aemet_last_obs("9434")
#> # A tibble: 13 × 25
#>    idema   lon fint                 prec   alt  vmax    vv    dv   lat  dmax ubi    pres    hr stdvv    ts pres_nmar
#>    <chr> <dbl> <dttm>              <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl>     <dbl>
#>  1 9434  -1.00 2025-02-09 23:00:00     0   249   1.7   1.1   142  41.7   133 ZARA…  995.    86   0.2   5.6     1025.
#>  2 9434  -1.00 2025-02-10 00:00:00     0   249   1.9   1.4   174  41.7   175 ZARA…  995.    87   0.2   6.3     1025.
#>  3 9434  -1.00 2025-02-10 01:00:00     0   249   2.1   1.7   124  41.7   130 ZARA…  994.    85   0.2   6.9     1025.
#>  4 9434  -1.00 2025-02-10 02:00:00     0   249   3.1   2     140  41.7   125 ZARA…  994.    84   0.2   6.3     1025.
#>  5 9434  -1.00 2025-02-10 03:00:00     0   249   2.9   1.6   147  41.7   140 ZARA…  994.    86   0.3   6.4     1024.
#>  6 9434  -1.00 2025-02-10 04:00:00     0   249   2.6   0.7    12  41.7   130 ZARA…  993.    86   0.3   6.3     1024.
#>  7 9434  -1.00 2025-02-10 05:00:00     0   249   3.9   0.4   138  41.7    95 ZARA…  993.    88   0.2   5.9     1024.
#>  8 9434  -1.00 2025-02-10 06:00:00     0   249   3.1   2.2   139  41.7   143 ZARA…  993.    87   0.4   6.1     1024.
#>  9 9434  -1.00 2025-02-10 07:00:00     0   249   2.6   1.2   177  41.7   138 ZARA…  993.    88   0.2   4.8     1024 
#> 10 9434  -1.00 2025-02-10 08:00:00     0   249   2.4   1.2   169  41.7   175 ZARA…  994.    88   0.3   5.9     1024.
#> 11 9434  -1.00 2025-02-10 09:00:00     0   249   2.1   1      92  41.7    75 ZARA…  994.    85   0.2   7.6     1024.
#> 12 9434  -1.00 2025-02-10 10:00:00     0   249   2.9   1.8   108  41.7    90 ZARA…  994.    78   0.4   9.5     1024.
#> 13 9434  -1.00 2025-02-10 11:00:00     0   249   4.1   2.6   116  41.7    85 ZARA…  994.    71   0.5  11.3     1024.
#> # ℹ 9 more variables: tamin <dbl>, ta <dbl>, tamax <dbl>, tpr <dbl>, stddv <dbl>, inso <dbl>, tss5cm <dbl>,
#> #   pacutp <dbl>, tss20cm <dbl>

Note that when possible, data representing dates and numbers are converted to the right format.

Spatial objects: sf

Another major change in v1.0.0 is the ability of return information on spatial sf format, using return_sf = TRUE. The coordinate reference system (CRS) used is EPSG 4326, that correspond to the World Geodetic System (WGS) and return coordinates in latitude/longitude (unprojected coordinates):

# You would need to install `sf` if not installed yet
# run install.packages("sf") for installation

library(ggplot2)
library(dplyr)

all_stations <- aemet_daily_clim(
  start = "2021-01-08", end = "2021-01-08",
  return_sf = TRUE
)


ggplot(all_stations) +
  geom_sf(aes(colour = tmed), shape = 19, size = 2, alpha = 0.95) +
  labs(
    title = "Average temperature in Spain",
    subtitle = "8 Jan 2021",
    color = "Max temp.\n(celsius)",
    caption = "Source: AEMET"
  ) +
  scale_colour_gradientn(
    colours = hcl.colors(10, "RdBu", rev = TRUE),
    breaks = c(-10, -5, 0, 5, 10, 15, 20),
    guide = "legend"
  ) +
  theme_bw() +
  theme(
    panel.border = element_blank(),
    plot.title = element_text(face = "bold"),
    plot.subtitle = element_text(face = "italic")
  )

Example: Temperature in Spain

Further enhancements

Other enhancements included on the v1.0.0:

All the functions are now vectorized.
New function get_metadata_aemet().
New function ggclimat_walter_lieth(). This function is now the default for climatogram_* functions . Old behavior can be reproduced with options ggplot2 = FALSE.
Plot functions gains new parameters (verbose and ...). Now it is possible to pass colors to the plotting functions.
New example datasets: climaemet::climaemet_9434_climatogram, climaemet::climaemet_9434_temp and climaemet::climaemet_9434_wind.