From 716278a44c74e6fefa1d2d173f4975fc5114b83d Mon Sep 17 00:00:00 2001
From: Bock <ABOCK@USGS.GOV>
Date: Tue, 28 Apr 2020 10:29:10 -0600
Subject: [PATCH] modified nsegment aggregation to result in correct US-DS node
sequencing
---
workspace/NHD_navigate.Rmd | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/workspace/NHD_navigate.Rmd b/workspace/NHD_navigate.Rmd
index fe2a879..e905aa8 100644
--- a/workspace/NHD_navigate.Rmd
+++ b/workspace/NHD_navigate.Rmd
@@ -352,12 +352,14 @@ if(needs_layer(out_gpkg, nsegment_raw)) {
# Write out nsegments composed of nhdflowlines
write_sf(ncombined, out_gpkg, nsegment_raw)
- # Aggregate flowlines per POI_ID to a single segment, carrying over useful information
- nsegments<-ncombined %>% group_by(POI_ID) %>% mutate(VA_MA = na_if(VA_MA, -9998)) %>%
+ # Dissolve flowlines to aggregated segments
+ nsegments <- ncombined %>% group_by(POI_ID) %>% arrange(desc(LevelPathI), desc(Hydroseq)) %>%
mutate(TT_Hours = (LENGTHKM * 3280.84) / (VA_MA * 3600)) %>%
- summarize(TotalLength = sum(LENGTHKM), TotalDA = max(TotDASqKM), HW = max(StartFlag), TT = sum(TT_Hours)) %>%
+ summarize(TotalLength = sum(LENGTHKM),TotalDA = max(TotDASqKM), HW = max(StartFlag), TT = sum(PathTimeMA),
+ do_union=FALSE) %>% st_cast("MULTILINESTRING") %>%
inner_join(st_drop_geometry(hucgagessegs) %>% select(COMID, Hydroseq, DnHydroseq), by = c("POI_ID" = "COMID"))
-
+
+
# produce a short data frame for populating TO_POI for downstream segment
to_from<-segPOIs %>% inner_join(st_drop_geometry(ncombined), by = c("DnHydroseq" = "Hydroseq")) %>%
select(COMID.x, Hydroseq, DnHydroseq, POI_ID)
--
GitLab