In FluvialGeomorph/fluvgeo: A Package for Fluvial Geomorphic Analysis

library(fluvgeo)
library(sf)
library(dplyr)
library(knitr)
library(pander)
library(ggplot2)
library(ggrepel)
library(tmap)
library(terra)
library(terrainr)

Study Area Overview

if(aerial) {
  print(fluvgeo::map_reach_overview(flowline = params$flowline_sf, 
                                    cross_section = params$cross_section_sf,
                                    background = "aerial",
                                    xs_label_freq = params$xs_label_freq,
                                    extent_factor = params$extent_factor))
}

if(elevation) {
  print(fluvgeo::map_reach_overview(flowline = params$flowline_sf, 
                                    cross_section = params$cross_section_sf,
                                    background = "elevation",
                                    xs_label_freq = params$xs_label_freq,
                                    exaggeration = params$exaggeration,
                                    extent_factor = params$extent_factor))
}

Longitudinal Profile

print(fluvgeo::compare_long_profile(stream = params$stream,
                            flowline_pts_sf_list = params$flowline_pts_sf_list,
                            features_sf = params$features_sf,
                            profile_units = params$profile_units))

print(fluvgeo::compare_xs_long_profile(stream = params$stream,
                                     xs_pts_sf_list = params$xs_points_sf_list,
                                     features_sf = params$features_sf,
                                     xs_label_freq = params$xs_label_freq,
                                     profile_units = params$profile_units))

Cross Section Metrics

print(fluvgeo::xs_metrics_plot_L1(xs_dims_sf = params$cross_section_sf,
                                  features_sf = params$features_sf,
                                  xs_label_freq = params$xs_label_freq,
                                  profile_units = params$profile_units))

Cross Section Profiles

if(params$show_xs_map) {
  # Get dem raster
  dem_rast <- fluvgeo::gdb_raster2SpatRast(raster_path = params$dem)
}

# Calculate cross section bearings and sequentially order
cross_section <- params$cross_section_sf %>%
  fluvgeo::xs_bearing() %>%
  dplyr::arrange(Seq)

# Iterate through cross sections
for(j in unique(cross_section$Seq)) {
  if(params$show_xs_map) {
    # Create the cross section plot
    print(fluvgeo::map_xs(cross_section = params$cross_section_sf,
                          xs_number = j,
                          dem = dem_rast,
                          extent_factor = params$extent_factor))
  }
  print(fluvgeo::xs_compare_plot_L1(stream = stream,
                                    xs_number = j,
                                    xs_pts_sf_list = xs_pts_sf_list))
  # Insert vertical white space so that next figure in loop is recognized
  cat('\n')
  # Must set chunk option`results='asis'` to ensure table is drawn properly
}