inst/extdata/gage_01013500_data.R

In dblodgett-usgs/hyRefactor: Hydrologic Network Refactoring Tools Based on HY_Features

# nolint start
# library(terra)
# library(dplyr)
# library(sf)
# library(nhdplusTools)
# 
# extdata   <- system.file("extdata", package = "hyRefactor")
# 
# # gage_01013500.gpkg turned into pre-processed sample data.
# base      = "/Volumes/Transcend/ngen/refactor-tests/base-data/"
# # RDS NHD files (dont keep local for space)
# cats      = readRDS(file.path(base, "catchments_all.rds"))
# fps       = readRDS(file.path(base, "nhdplus_flowline_update_sb.rds"))
# 
# UT_COMIDs = get_UT(st_drop_geometry(fps), 724696)
# flowpaths = make_standalone(filter(fps, COMID %in% UT_COMIDs))
# catchments = filter(cats, FEATUREID %in% UT_COMIDs)
# 
# fac <- terra::rast(file.path(base, "fdrfac/01a_fac.tif"))
# 
# fac = crop(fac, st_transform(catchments, crs(fac)), snap = "out")
# 
# fdr = terra::rast(file.path(base, "fdrfac/01a_fdr.tif")) %>%
#   crop(st_transform(catchments, crs(fac)), snap = "out")
# 
# 
# refactor_nhdplus(nhdplus_flines = flowpaths,
#                              split_flines_meters         = 10000,
#                              collapse_flines_meters      = 1000,
#                              collapse_flines_main_meters = 1000,
#                              split_flines_cores          = 2,
#                              out_refactored              = "refactor.gpkg",
#                              out_reconciled              = "reconcile.gpkg",
#                              three_pass                  = TRUE,
#                              purge_non_dendritic         = FALSE,
#                              warn                        = FALSE)
# 
# fline_ref <-  st_transform(read_sf("refactor.gpkg"), st_crs(fac))
# fline_rec <-  st_transform(read_sf("reconcile.gpkg"), st_crs(fac))
# 
# catchments = st_transform(catchments, st_crs(fac))
# st_precision(catchments) <- terra::res(fdr)[1]
# 
# raw   <- reconcile_catchment_divides(catchment = catchments,
#                                               fline_ref = fline_ref,
#                                               fline_rec = fline_rec,
#                                               fdr       = fdr,
#                                               fac       = fac,
#                                               para      = 1,
#                                               cache     = NULL,
#                                               fix_catchments = FALSE)
# 
# write_sf(raw, file.path(extdata, "gage_01013500.gpkg"), "raw-divides")
# unlink(list("refactor.gpkg", "reconcile.gpkg"))