In KWB-R/kwb.satellite: R Package for Working with Satellite Data from Various Providers (Copernicus, GoogleEarthEngine)

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

1. Register for CDS services

https://cds.climate.copernicus.eu/user/register

2. Setting Access Rights

Copy UUID and API KEY from: https://cds.climate.copernicus.eu/user

Add these to your R environment file by calling usethis::edit_r_environ() and add the following two key value pairs in R the environment file:

"COPERNICUS_CLIMATE_UUID" = "myuuid"
"COPERNICUS_CLIMATE_APIKEY" = "myapikey"

Restart R(Studio) now and set your access details as shown below:

ecmwfr::wf_set_key(user = Sys.getenv("COPERNICUS_CLIMATE_UUID"),
                   key = Sys.getenv("COPERNICUS_CLIMATE_APIKEY"),
                   service = "cds")

4. Create a query

Go to: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form

to create a query.

Copy to Rstudio -> mark all -> goto Tools -> Addins -> Browse_Addins -> ecmwfr -> Python to Mars -> in order to convert query to R list (see: https://cran.r-project.org/web/packages/ecmwfr/vignettes/cds_vignette.html for details)

5. Run your query

request <- list(
    product_type = "reanalysis",
    dataset_short_name = "reanalysis-era5-single-levels",
    file_format = "grib"
    variable = "2m_temperature",
    year = c("2010", "2011", "2012", "2013", "2014", 
             "2015", "2016", "2017", "2018", "2019", "2020"),
    month = c("01", "02", "03", "04", "05", "06", 
              "07", "08", "09", "10", "11", "12"),
    day = c("01", "02", "03", "04", "05", "06", "07", "08", "09", "10",
            "11", "12", "13", "14", "15", "16", "17", "18", "19", "20",
            "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31"),
    time = c("00:00", "01:00", "02:00", "03:00", "04:00", "05:00", "06:00",
             "07:00", "08:00", "09:00", "10:00", "11:00", "12:00", "13:00",
             "14:00", "15:00", "16:00", "17:00", "18:00", "19:00", "20:00",
             "21:00", "22:00", "23:00"),
    area = c(40, 116, 39, 117),
    target =  "download.grib")


  ecmwfr::wf_request(request = request,
                     transfer = TRUE,
                     path = export_dir,
                     verbose = FALSE)

Take a break and wait until request is finished this may take a few hours. To check the progress of the request visit: https://cds.climate.copernicus.eu/cdsapp#!/yourrequests

6. Load raster data

grib_file <- "download.grib"

meta <- kwb.satellite::get_metadata_era5(grib_file)

temperature <- raster::raster(grib_file)

KWB-R/kwb.satellite documentation built on May 29, 2022, 9:08 p.m.