remove unneeded attributes, fixes #126. add missing flowlines from burn line...

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

@@ -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")
   
   }