stream <- "Cole Creek R1"
flowline_points_1 <- file.path(system.file("extdata", "y2016_R1.gdb",
package = "fluvgeo"),
"feature_dataset/flowline_points")
flowline_points_2 <- file.path(system.file("extdata", "y2010_R1.gdb",
package = "fluvgeo"),
"feature_dataset/flowline_points")
flowline_points_3 <- file.path(system.file("extdata", "y2006_R1.gdb",
package = "fluvgeo"),
"feature_dataset/flowline_points")
flowline_points_4 <- NULL
survey_name_1 <- "2016"
survey_name_2 <- "2010"
survey_name_3 <- "2006"
survey_name_4 <- NULL
features_fc <- file.path(system.file("extdata", "y2016_R1.gdb",
package = "fluvgeo"),
"feature_dataset/features")
profile_units <- "feet"
# Create list of survey paths
flowline_points_paths <- list(flowline_points_1, flowline_points_2,
flowline_points_3, flowline_points_4)
# Name the survey paths list by their survey names
survey_names <- c(survey_name_1, survey_name_2, survey_name_3, survey_name_4)
flowline_points_paths <- setNames(flowline_points_paths, survey_names)
# Eliminate empty surveys
flowline_points_paths <- purrr::discard(flowline_points_paths, is.null)
# Convert list of survey paths to list of sf objects
flowline_pts_sf_list <- purrr::map(flowline_points_paths, fluvgeo::fc2sf)
# Convert features_fc to an sf object
features_sf <- fluvgeo::fc2sf(features_fc)
# Call the graph function
p <- fluvgeo::compare_long_profile(stream = stream,
flowline_pts_sf_list = flowline_pts_sf_list,
features_sf = features_sf,
profile_units = profile_units)
#print(p)
test_that("compare_long_profile plot exists", {
expect_true("ggplot" %in% class(p))
expect_error(print(p), NA)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.