Get the shoreline from a bathymetric raster

Usage

get_shoreline(bathy_raster, dTolerance = NULL)

Arguments

bathy_raster: SpatRaster object with the bathymetry data.
dTolerance: numeric. The tolerance for simplifying the shoreline. Default is the maximum resolution of the bathymetric raster. If set to 0, the shoreline will not be simplified.

Value

sf POLYGON of the shoreline.

Examples

shoreline <- readRDS(system.file("extdata/rotoma_shoreline.rds",
package = "bathytools"))
depth_points <- readRDS(system.file("extdata/depth_points.rds",
package = "bathytools"))
bathy_raster <- rasterise_bathy(shoreline = shoreline,
depth_points = depth_points, crs = 2193, res = 8)
#> No islands found.
#> ℹ Generating depth points for interpolation
#> ✔ Generating depth points for interpolation [177ms]
#> 
#> ℹ Interpolating depth points to raster
#> ℹ Adjusting depths >= 0 to -0.82m
#> ℹ Interpolating depth points to raster


#> ✔ Interpolating depth points to raster [1.2s]
#> 
shoreline2 <- get_shoreline(bathy_raster = bathy_raster)