Learn R Programming

rLakeHabitat (version 1.0.0)

littoralVol: Calculate Littoral Volume

Description

Calculate littoral and pelagic volume across water levels from a DEM based on estimated photic depth.

Usage

littoralVol(
  DEM,
  photic,
  secchi = NULL,
  DEMunits = "m",
  depthUnits = "ft",
  by = 1
)

Value

data frame of littoral and pelagic volume estimates

Arguments

DEM

SpatRaster object of a given waterbody, rasters can be transformed to SpatRaster via the rast() function in 'terra'

photic

number giving the average photic depth, overwrites Secchi depth

secchi

number giving the average secchi depth, photic zone estimated as 2.6m * secchi

DEMunits

character describing units of raster coordinate system. Can be meters, kilometers, or hectares ("m", "km", "ha"), default = "m"

depthUnits

character describing units of depth measurement. Can be either feet or meters ("ft", "m"), default = "ft"

by

numeric increment per unit by which volumes are calculated. Higher values will result in lower resolution. Default = 1

Author

Tristan Blechinger, Department of Zoology & Physiology, University of Wyoming

Examples

Run this code
#load raster
DEM <- terra::rast(system.file("extdata", "example_raster.tif", package = 'rLakeHabitat'))
#run function
littoralVol(DEM, photic = 2, DEMunits = "m", depthUnits = "m", by = 1)

Run the code above in your browser using DataLab