calculate_nominal_resolution: Calculate "nominal resolution" of a grid

In DrylandEcology/rSW2st: Spatiotemporal Tools for SOILWAT2 and STEPWAT2 Simulation Experiments

Calculate "nominal resolution" of a grid

Description

Calculate "nominal resolution" of a grid

Usage

calculate_nominal_resolution(grid, maskvalue = NA)

Arguments

grid	A terra::rast, a raster::Raster or a stars::stars object. Gridcells with values that are not equal to a maskvalue are included in the "nominal resolution" calculation.
maskvalue	A vector. Values to mask out from grid.

Value

A character string with a "nominal resolution" value in kilometers.

References

CMIP6 Global Attributes, DRS, Filenames, Directory Structure, and CV’s, 10 September 2018 (v6.2.7). Appendix 2: Algorithms for Defining the "nominal_resolution" Attribute

Examples

r1 <- terra::rast(
  xmin = -120, xmax = -90,
  ymin = 30, ymax = 50,
  crs = "OGC:CRS84",
  resolution = c(0.5, 0.5)
)
xy <- terra::spatSample(r1, size = 200L, as.points = TRUE)
r1[xy] <- 1

calculate_nominal_resolution(r1)
calculate_nominal_resolution(stars::st_as_stars(r1))

r2 <- terra::rast(
  xmin = -2480000, xmax = 90000,
  ymin = 650000, ymax = 4020000,
  crs = "EPSG:6350",
  resolution = c(10000, 10000)
)
xy <- terra::spatSample(r2, size = 200L, as.points = TRUE)
r2[xy] <- 1

calculate_nominal_resolution(r2)
calculate_nominal_resolution(stars::st_as_stars(r2))

DrylandEcology/rSW2st documentation built on Jan. 10, 2024, 6:22 p.m.