pRecipe (version 3.0.1-3)

subset_data: Subset data in space and time

Description

The function subset_data subsets the data in space within a bounding box, and/or in time within a year range.

Usage

subset_data(x, box = NULL, yrs = NULL)
# S4 method for Raster
subset_data(x, box = NULL, yrs = NULL)
# S4 method for data.table
subset_data(x, box = NULL, yrs = NULL)
# S4 method for character
subset_data(x, box = NULL, yrs = NULL)

Value

Raster* object; data.table

Arguments

x

Raster* object; data.table (see details); filename (character, see details)

box

numeric. Bounding box in the form: (xmin, xmax, ymin, ymax)

yrs

numeric. Time range in the form: (start_year, end_year)

Details

If `x` is a data.table, its columns should be named: "lon", "lat", "date", and "value"

If `x` is a filename, it should point to a *.nc file.

If subsetting only in space or time then the arguments must be passed by name. I.e., subset_data(x, box = ...) (space) or subset_data(x, yrs = ...) (time)

Examples

if (FALSE) {
download_data("gldas-vic", tempdir(), timestep = "yearly")
r <- raster::brick(paste0(tempdir(),
"/gldas-vic_tp_mm_land_194801_201412_025_yearly.nc"))
sd <- subset_data(r, c(12.24, 18.85, 48.56, 51.12), c(2000, 2010))
ss <- subset_data(r, box = c(12.24, 18.85, 48.56, 51.12))
st <- subset_data(r, yrs = c(2000, 2010))
}