From eb0690b00326f30fecf1c12b5651a2842e5d01bb Mon Sep 17 00:00:00 2001
From: Bock <abock@usgs.gov>
Date: Thu, 30 Jun 2022 10:28:41 -0600
Subject: [PATCH] Added wb events to split events
---
workspace/05_hyRefactor_flines.Rmd | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/workspace/05_hyRefactor_flines.Rmd b/workspace/05_hyRefactor_flines.Rmd
index 6efa5e6..1078df2 100644
--- a/workspace/05_hyRefactor_flines.Rmd
+++ b/workspace/05_hyRefactor_flines.Rmd
@@ -72,7 +72,14 @@ events <- read_sf("cache/gages_info.gpkg", "Gages") %>%
filter(reach_meas - FromMeas > 10 & AreaSqKM > 2 &
ToMeas - reach_meas > 10 & LENGTHKM > (combine_meters / 1000)) %>%
select(COMID = comid, REACHCODE = reachcode, REACH_meas = reach_meas, site_no) %>%
- filter(!COMID %in% avoid$COMID)
+ filter(!COMID %in% avoid$COMID) %>%
+ mutate(event_type = "streamgage")
+
+if(!needs_layer(nav_gpkg, "WB_events")){
+ wb_events <- read_sf(nav_gpkg, "WB_events") %>%
+ filter(COMID %in% nhd$COMID)
+ events <- rbind(events, select(wb_events, COMID, REACHCODE, REACH_meas, site_no = WBAREACOMI, event_type))
+}
# Write out events and outelts
write_sf(events, out_refac_gpkg, split_layer)
--
GitLab