AdjSD: Calculates standard deviation of bathymetry (a measure of...

In MultiscaleDTM: Multi-Scale Geomorphometric Terrain Attributes

Calculates standard deviation of bathymetry (a measure of rugosity) adjusted for slope

Description

Calculates standard deviation of bathymetry (a measure of rugosity). Using a sliding rectangular window a plane is fit to the data and the standard deviation of the residuals is calculated (Ilich et al., 2023)

Usage

AdjSD(
  r,
  w = c(3, 3),
  na.rm = FALSE,
  include_scale = FALSE,
  filename = NULL,
  overwrite = FALSE,
  wopt = list()
)

Arguments

r	DTM as a SpatRaster or RasterLayer in a projected coordinate system where map units match elevation/depth units
w	A vector of length 2 specifying the dimensions of the rectangular window to use where the first number is the number of rows and the second number is the number of columns. Window size must be an odd number.
na.rm	A logical indicating whether or not to remove NA values before calculations
include_scale	logical indicating whether to append window size to the layer names (default = FALSE)
filename	character Output filename.
overwrite	logical. If TRUE, filename is overwritten (default is FALSE).
wopt	list with named options for writing files as in writeRaster

Value

a SpatRaster or RasterLayer of adjusted rugosity

References

Ilich, A. R., Misiuk, B., Lecours, V., & Murawski, S. A. (2023). MultiscaleDTM: An open-source R package for multiscale geomorphometric analysis. Transactions in GIS, 27(4). https://doi.org/10.1111/tgis.13067

Examples

r<- rast(volcano, extent= ext(2667400, 2667400 + 
ncol(volcano)*10, 6478700, 6478700 + nrow(volcano)*10), 
crs = "EPSG:27200")
adjsd<- AdjSD(r, w=c(5,5), na.rm = TRUE)
plot(adjsd)

MultiscaleDTM documentation built on April 4, 2025, 5:17 a.m.