ExtractVeg: Extract Vegegation

In essatech/MNAI.CPBT: Municipal Natural Assets Initiative - Coastal Protection Benefit Toolbox

Extract Vegegation

Description

Merges the cleaned cross-shore profile dataset with vegetation data. Note that this function must be run before the wave evolution and erosion models can be run. Even if your project has no vegetation, running this function is still necessary to format the data frame with appropriate columns.

Usage

ExtractVeg(pt_exp = NA, Vegetation = NA)

Arguments

pt_exp	Cross-shore point elevations spatial points object returned from CleanTransect.
Vegetation	A simple feature polygon of class sf and dataframe. Patch attributes must include hc Numeric, Blade height in meters N Numeric, Shoot density as number of shoots per meter squared. d Numeric, Blade width in meters. Type Either 'Eelgrass', 'Kelp' or 'Marsh'. Cd Numeric, Drag coefficent as per Guannel et al 2015. If unsure use a value of 0.1 for Eelgrass

Value

An object of class sf and data.frame updated with vegetation data attributes and ready for the wave evolution and erosion model.

Examples

## Not run: 
library(MNAI.CPBT)
data(Coastline)
# Generate cross-shore profile lines along the coastline.
crossshore_profiles <- samplePoints(
  Coastline = Coastline,
  ShorelinePointDist = 150,
  BufferDist = 50,
  RadLineDist = 1.5
)
crossshore_lines <- crossshore_profiles[[2]]

# Extract elevation values along each profile
rpath <-  system.file("extdata", "TopoBathy.tif", package = "MNAI.CPBT")
TopoBathy <- raster::raster(rpath)
pt_elevs <- ExtractElev(crossshore_lines, TopoBathy)

# Run SignalSmooth function to smooth elevation profiles
pt_elevs <- SignalSmooth(point_elev = pt_elevs,
SmoothParameter = 5)

# Clean the cross-shore profiles with CleanTransect
cleantransect <- CleanTransect(
  point_elev = pt_elevs,
  RadLineDist = 1.5,
  MaxOnshoreDist = 1,
  trimline = NA
)

# Merge vegetation onto lines
data(Vegetation)
dat_veg <- ExtractVeg(pt_exp = cleantransect, Vegetation = Vegetation)

# Cross-shore profiles ready for modeling

## End(Not run)

essatech/MNAI.CPBT documentation built on July 1, 2023, 12:34 p.m.