From fc5fe2b071765ad9691ece32f7570773a08e8034 Mon Sep 17 00:00:00 2001
From: David Blodgett <dblodgett@usgs.gov>
Date: Wed, 21 Jun 2023 20:07:57 -0500
Subject: [PATCH] remove unneeded attributes, fixes #126. add missing flowlines
from burn line event, fixes #125 fixes #123. Remove duplicate comid rows,
fixes #124.
---
workspace/R/utils.R | 10 ++++++++--
1 file changed, 8 insertions(+), 2 deletions(-)
diff --git a/workspace/R/utils.R b/workspace/R/utils.R
index 8c0d6da..e7727d9 100644
--- a/workspace/R/utils.R
+++ b/workspace/R/utils.R
@@ -27,8 +27,8 @@ fix_headwaters <- function(nhd_rds, out_gpkg, new_atts = NULL,
ble <- sf::st_zm(sf::st_as_sf(ble))
nhd <- sf::st_zm(nhd)
- sf::st_geometry(nhd)[!sf::st_is_empty(sf::st_geometry(ble)) & (nhd$StartFlag == 1 | nhd$Divergence == 2)] <-
- sf::st_geometry(ble)[!sf::st_is_empty(sf::st_geometry(ble)) & (nhd$StartFlag == 1 | nhd$Divergence == 2)]
+ sf::st_geometry(nhd)[!sf::st_is_empty(sf::st_geometry(ble)) & (nhd$StartFlag == 1)] <-
+ sf::st_geometry(ble)[!sf::st_is_empty(sf::st_geometry(ble)) & (nhd$StartFlag == 1)]
if(!is.null(new_atts)) {
new_atts <- data.table::fread(new_atts,
@@ -109,6 +109,12 @@ fix_headwaters <- function(nhd_rds, out_gpkg, new_atts = NULL,
nhd <- select(nhd, -override_tocomid)
+ # remove remnant attributes
+ nhd <- select(nhd, -FromNode, -ToNode, -StartFlag -StreamCalc, -Divergence, -DnMinorHyd)
+
+ warning("naive duplication!")
+ nhd <- nhd[!duplicated(nhd$COMID), ]
+
sf::write_sf(nhd, out_gpkg, "reference_flowlines")
}
--
GitLab