coast_rates estimate the main parameters about the coastal variations, as NSM, LRR or WLR.
Arguments
- inter_dist
Shapefile. A point format shapefile with the intersections between each transect and the shoreline.
- normals
Shapefile. A polyline format shapefile with all normal lines (transects) included in the study site.
- table
CSV. A table with the information about the coastlines dates in a column named "Day" in format (YYYYY-mm-dd). If the information about the acquisition hour is available, the user should introduce a column called "Hour" in format (HH:MM:SS) and the associated uncertainty to each coastline in meters. The column names should be "Day", "Hour" and "Uncertainty".
- out_name
Integer. Output name for the resulting shapefile with the rates of each transect.
Value
A polyline shapefile with all rates associated to each transect.
A csv file with the central tendency stats and quantiles .25, .75 and .90 for each parameter.
Details
This function is part of ODSAS method (Gomez-Pazo et al., 2022) and estimates the main five key measures for each individual transect and also produce aggregates
for all transects identically as within the DSAS tool (Himmelstoss et al., 2018).
The rates and their units are explained in detail in CoastCR documentation.
References
Gomez-Pazo, A., Payo, A., Paz-Delgado, M.V., Delgadillo-Calzadilla, M.A. (2022) Open Digital Shoreline Analysis System: ODSAS v1.0 Journal of Marine Science and Engineering, 10, 26.Himmelstoss, E.A., Henderson, R.E., Kratzmann, M.G., Farris, A.S. (2018) Digital Shoreline Analysis System (DSAS) version 5.0 user guide US Geological Survey Open File Report 2018-1179, 110 pp
Examples
library(sf)
library(CoastCR)
#Normal lines shapefile
normals <- st_read(system.file("extdata/normals.shp", package = "CoastCR"))
#> Reading layer `normals' from data source
#> `/home/runner/work/_temp/Library/CoastCR/extdata/normals.shp'
#> using driver `ESRI Shapefile'
#> Simple feature collection with 16 features and 7 fields
#> Geometry type: LINESTRING
#> Dimension: XY
#> Bounding box: xmin: 508009.1 ymin: 4674223 xmax: 508116 ymax: 4674330
#> Projected CRS: ETRS89 / UTM zone 29N
#Table with dates and associated uncertainty
table <- read.csv(system.file("extdata/table_coastlines.csv", package = "CoastCR"))
#Filtered intersections shapefile
inter_dist <- st_read(system.file("extdata/dist.shp", package = "CoastCR"))
#> Reading layer `dist' from data source
#> `/home/runner/work/_temp/Library/CoastCR/extdata/dist.shp'
#> using driver `ESRI Shapefile'
#> Simple feature collection with 128 features and 4 fields
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: 508058.2 ymin: 4674243 xmax: 508101 ymax: 4674329
#> Projected CRS: ETRS89 / UTM zone 29N
#Define output name
out_name <- tempfile("normals_coast_rates", fileext = ".shp")
coast_rates(inter_dist, normals, table, out_name)
#> Writing layer `normals_coast_rates199142846ef2' to data source
#> `/tmp/RtmpwOQn8r/normals_coast_rates199142846ef2.shp' using driver `ESRI Shapefile'
#> Writing 16 features with 9 fields and geometry type Line String.