calcLittoral: Calculate Littoral Area

In rLakeHabitat: Interpolate Bathymetry and Quantify Physical Aquatic Habitat

Calculate Littoral Area

Description

Calculates littoral surface area (2D) of a given waterbody across water levels based on an average photic depth value.

Usage

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

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 (secchi and DEM). Can be either feet or meters ("ft", "m"), default = "ft"
by	numeric increment per unit depth by which areas are calculated. Higher values will result in lower resolution. Default = 1
stop	optional numeric value specifying depth at which to stop calculations, default = NULL

Value

data frame of areas in specified units for each depth, as well as the littoral percentage of total surface area

Author(s)

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

Examples

#load raster
DEM <- terra::rast(system.file("extdata", "example_raster.tif", package = 'rLakeHabitat'))
#run function
calcLittoral(DEM, secchi = 1, depthUnits = "m", DEMunits = "m")

rLakeHabitat documentation built on April 16, 2025, 1:10 a.m.