powered by
Calculate littoral and pelagic volume across water levels from a DEM based on estimated photic depth.
littoralVol( DEM, photic, secchi = NULL, DEMunits = "m", depthUnits = "ft", by = 1 )
data frame of littoral and pelagic volume estimates
SpatRaster object of a given waterbody, rasters can be transformed to SpatRaster via the rast() function in 'terra'
number giving the average photic depth, overwrites Secchi depth
number giving the average secchi depth, photic zone estimated as 2.6m * secchi
character describing units of raster coordinate system. Can be meters, kilometers, or hectares ("m", "km", "ha"), default = "m"
character describing units of depth measurement. Can be either feet or meters ("ft", "m"), default = "ft"
numeric increment per unit by which volumes are calculated. Higher values will result in lower resolution. Default = 1
Tristan Blechinger, Department of Zoology & Physiology, University of Wyoming
#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