The leri R package facilitates access to the NOAA Landscape Evaporative Response Index data product.
You can install the development version of leri with devtools:
# install.packages("devtools")
devtools::install_github("earthlab/leri")
Or, you can install the stable version from CRAN:
install.packages("leri")
The LERI product is available from the year 2000 to present at a 1 km spatial resolution over the continental United States, at the following timescales:
If we were interested in medium term drought, we could get 3 month timescale LERI data for a particular date as follows:
library(leri)
r <- get_leri(date = "2018-09-01", product = "3 month")
r
#> class : RasterLayer
#> dimensions : 2844, 6617, 18818748 (nrow, ncol, ncell)
#> resolution : 0.009, 0.009 (x, y)
#> extent : -126.009, -66.456, 23.949, 49.545 (xmin, xmax, ymin, ymax)
#> crs : +proj=longlat +datum=WGS84 +no_defs
#> source : memory
#> names : LERI_03mn_20180901.nc
#> values : 3.125, 96.87499 (min, max)
The get_leri
function will download the NetCDF files from the NOAA ftp
server and read them as RasterStack
objects, which facilitates any
other operations that users may need via the raster
package. For
example, we could plot the data using the raster::plot
function:
library(raster)
#> Loading required package: sp
library(viridis)
#> Loading required package: viridisLite
plot(r, col = cividis(10, direction = -1),
main = "3 month LERI data for 2018-09-01")
More information on the creation and recommended usage for LERI can be found here: https://psl.noaa.gov/leri/
Other potentially useful resources include:
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.