rLakeHabitat (version 1.0.0)

calcLittoral: 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
)

Value

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

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

Author

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")