#' @title Road Proximity Sample
#'
#' @description
#' Road Proximity Sample for Site Type B
#'
#' @details
#' WSEP Tier II Road Proximity Sample for Site Type B
#' \itemize{
#' \item{Apply a buffer (20m for < 3rd order and 40 m for ≥ 3rd order) streams.}
#' \item{Taking the intersection of roads and this buffer}
#' \item{Removing any stream crossings (Site Type A), by excluding any cases within 100 m of a crossing to avoid double counting}
#' \item{Removing any segments < 50m in length unless they are near a switch-back (determined by manual review of the map). Switch-back determination is an optional manual step not automated in R code}
#' \item{Providing the start point of the segment as well as the segment length and strata id (< 3rd order vs. ≥ 3rd order) and mapping the entire segment on the field maps to facilitate sampling}
#' \item{Generate a random sample from the Site Type B list the strata with < 3rd order streams}
#' \item{ Append the complete list of Site Type B from the ≥ 3rd order strata }
#' \item{ Create a field checklist with at least the following fields: unique identifier (e.g., WatershedName_SD_B_001), coordinates of start point and end point, segment length, and Strata }
#' }
#' Note that there is no repeat sampling or replacement. If the sample size `n` exceeds the number of segments available for a stratum, then only the limited number of available locations will be returned.
#'
#' @param n Sample size n for each stratum. Defaults to 20.
#' @param strm A streamline data set of class `sf` from the BCFWA. Run `utm_projection()` first to ensure projection is in meters.
#' @param roads A road network data set of class `sf`. Run `utm_projection()` first to ensure projection is in meters.
#' @param buffer_s1_m Numeric. Buffer size (in meters) for Strata 1 streams < 3rd order. Defaults to 20 m.
#' @param buffer_s2_m Numeric. Buffer size (in meters) for Strata 2 streams >= 3rd order. Defaults to 40 m.
#' @param buffer_crossings_m Numeric. Exclusion buffer size (in meters) to stream-road crossings to exclude sample sites from Site Type A. Defaults to 100 m.
#' @param small_strm_segment_m Numeric. Remove small stream-road buffered segments unless they are near a switch-back (determined manually and added at end). Default small segment size is set to 50 m.
#' @param stream_order String. Column name for stream order. Defaults to `STREAM_ORDER`. Rename as needed.
#' @importFrom magrittr %>%
#'
#' @return
#' A object of class `sf` of sample features representing sample sites.
#'
#'
#' @export
<- (
n = 20,
strm = ,
roads = ,
buffer_s1_m = 20,
buffer_s2_m = 40,
buffer_crossings_m = 100,
small_strm_segment_m = 50,
stream_order = "STREAM_ORDER"
) {
# Bind for non-standard
tmp_group_name <-
L1 <-
# Fix stream order column name
(!(stream_order %in% (strm))) {
(("Could not find a column named ", stream_order, ". Check stream_order argument"))
}
vals <- strm, stream_order]
sf::st_geometry(vals) <-
vals <- vals, 1]
vals <- (vals)
strm$tmp_stream_order <- vals
# Recreate stream road crossings
sc2 <- ({ sf::st_intersection(roads, strm) })
# Apply buffer - group s1
ssub <- strm[which(strm$tmp_stream_order < 3), ]
# make placeholder dataframe
skeep_sb_s1 <- ()
buff <- sf::st_buffer(roads, buffer_s1_m) # 20 m
buff <- buff %>% dplyr::group_by() %>% dplyr::summarise()
# Taking the intersection of roads and this buffer,
("Stratum 1 streams and road buffer...")
rint <- ({ sf::st_intersection(ssub, buff) })
rint <- ({ sf::st_cast(rint, "LINESTRING") })
rint$tmp <- "samp"
rint <- rint, ("tmp")]
# Removing any stream crossings (site Type A)
# by excluding any cases within 100m of a
# crossing to avoid double counting,
sc2_100 <- sf::st_buffer(sc2, buffer_crossings_m)
sc2_100 <- sc2_100 %>% dplyr::group_by() %>% dplyr::summarise()
sdiff <- ({ sf::st_difference(rint, sf::st_union(sc2_100)) })
# Removing any segments <50m in length
# unless they are near a switch-back
# (determined by manual review of the map),
sdiff$length_m <- sf::st_length(sdiff)
sdiff$length_m <- (sdiff$length_m)
skeep <- sdiff[which(sdiff$length_m > small_strm_segment_m), ]
skeep_sb_s1 <- skeep
(((skeep_sb_s1), " stream segments for Site Type B - Stratum 1"))
((skeep_sb_s1) > 0) {
skeep_sb_s1$site_id <- ("B_", 1:(skeep_sb_s1), "_stratum_1")
skeep_sb_s1$strata <- "stratum_1"
skeep_sb_s1$ <- "type_b"
skeep_sb_s1$tmp <-
skeep_sb_s1 <- skeep_sb_s1, ("site_id", "strata", "type", "length_m")]
}
# Apply buffer - group s1
ssub <- strm[which(strm$tmp_stream_order >= 3), ]
skeep_sb_s2 <- ()
# Only run if stream segments available
((ssub) > 0) {
buff <- sf::st_buffer(roads, buffer_s1_m) # 20 m
buff <- buff %>% dplyr::group_by() %>% dplyr::summarise()
# taking the intersection of roads and this buffer,
("Stratum 2 streams and road buffer...")
rint <- ({ sf::st_intersection(ssub, buff) })
rint <- ({ sf::st_cast(rint, "LINESTRING") })
((rint) > 0 ) {
rint$tmp <- "samp"
rint <- rint, ("tmp")]
# (c) removing any stream crossings (site Type A)
# by excluding any cases within 100m of a
# crossing to avoid double counting
sc2_100 <- sf::st_buffer(sc2, buffer_crossings_m)
sc2_100 <- sc2_100 %>% dplyr::group_by() %>% dplyr::summarise()
sdiff <- ({ sf::st_difference(rint, sf::st_union(sc2_100)) })
# removing any segments <50m in length unless
# they are near a switch-back
# (determined by manual review of the map),
sdiff$length_m <- sf::st_length(sdiff)
sdiff$length_m <- (sdiff$length_m)
skeep <- sdiff[which(sdiff$length_m > small_strm_segment_m), ]
skeep_sb_s2 <- skeep
(((skeep_sb_s2), " stream segments for Site Type B - Stratum 2"))
((skeep_sb_s2) > 0) {
skeep_sb_s2$site_id <- ("B_", 1:(skeep_sb_s2), "_stratum_2")
skeep_sb_s2$strata <- "stratum_2"
skeep_sb_s2$ <- "type_b"
skeep_sb_s2$tmp <-
skeep_sb_s2 <- skeep_sb_s2, ("site_id", "strata", "type", "length_m")]
}
}
}
# Merge together
((skeep_sb_s1) == 0 & (skeep_sb_s2) == 0) {
("No stream segments remaining - adjust buffer sizes")
}
((skeep_sb_s1) == 0 & (skeep_sb_s2) > 0) {
skeep_all <- skeep_sb_s2
}
((skeep_sb_s1) > 0 & (skeep_sb_s2) == 0) {
skeep_all <- skeep_sb_s1
}
((skeep_sb_s1) > 0 & (skeep_sb_s2) > 0) {
skeep_all <- (skeep_sb_s1, skeep_sb_s2)
}
skeep_all$length_m <- (skeep_all$length_m, 1)
# Sample by strata
s1 <- skeep_all$site_id[which(skeep_all$strata == "stratum_1")]
((s1) > n) {
s1_samp <- (s1, size = n, = )
} {
s1_samp <- s1
}
s2 <- skeep_all$site_id[which(skeep_all$strata == "stratum_2")]
((s2) > n) {
s2_samp <- (s2, size = n, = )
} {
s2_samp <- s2
}
samp_ids <- (s1_samp, s2_samp)
# Point sample
ssamp <- skeep_all[which(skeep_all$site_id %in% samp_ids), ]
my_geom <- ({ sf::st_cast(ssamp, "LINESTRING") })
<- (sf::st_coordinates(my_geom))
<- ()
first_pt <- %>% dplyr::group_by(L1) %>%
dplyr::slice(1)
first_pt <- (first_pt)
first_pt_sf <- sf::st_as_sf(first_pt, coords = ("X", "Y"))
sf::st_crs(first_pt_sf) <- sf::st_crs(strm)$epsg
((first_pt_sf) != (ssamp)) {
("bad coordinate breakdown...")
}
# Add column names
first_pt_sf$site_id <- ssamp$site_id
first_pt_sf$strata <- ssamp$strata
first_pt_sf$ <- ssamp$
first_pt_sf$length_m <- ssamp$length_m
first_pt_sf$L1 <-
first_pt_sf <- first_pt_sf, ("site_id", "strata", "type", "length_m")]
# Add on stream order
int_t <- ({ sf::st_intersection(strm,
sf::st_buffer(first_pt_sf, "site_id"],
= 10))
})
add_order <- int_t, ("site_id", stream_order)]
sf::st_geometry(add_order) <-
add_order <- add_order!((add_order)), ]
add_order$stream_order <- (add_order, stream_order])
add_order <- add_order, ("site_id", "stream_order")]
first_pt_sf_so <- (first_pt_sf, add_order,
by.x = "site_id",
by.y = "site_id",
all.x = ,
all.y = )
# order by stream order then site ID
first_pt_sf_so <- first_pt_sf_so[order(first_pt_sf_so$stream_order), ]
# Return sample points and sample lines
ret_obj <- ()
ret_obj$ <- first_pt_sf_so
ret_obj$line_segments <- skeep_all
(("Total sample sites: ", (skeep_all)))
(ret_obj)
}
