From 5b15745cf8189772b9563c12793b54baf7cbce4c Mon Sep 17 00:00:00 2001
From: David Blodgett <dblodgett@usgs.gov>
Date: Thu, 7 Mar 2024 15:38:26 -0600
Subject: [PATCH] tidy up get data
---
workspace/00_get_data.Rmd | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/workspace/00_get_data.Rmd b/workspace/00_get_data.Rmd
index ab03de5..183f7a5 100644
--- a/workspace/00_get_data.Rmd
+++ b/workspace/00_get_data.Rmd
@@ -378,10 +378,12 @@ VAA_rds_path <- file.path(nhdplus_dir, "nhd_vaa.rds")
if(!file.exists(VAA_rds_path)) {
message("downloading NHDPlusVAA...")
- get_vaa(path = nhdplus_dir)
+ download_vaa(VAA_fst_path)
+
# Extract other attributes
nhdpv2_vaa <- read_sf(nhdplus_gdb, "NHDFlowline_Network") %>%
st_drop_geometry()
+
saveRDS(nhdpv2_vaa, file.path(nhdplus_dir, "nhd_vaa.rds"))
}
@@ -587,9 +589,9 @@ if(!dir.exists(merit_dir)) {
merit_fac = get_sbfile(file.path(merit_dir,
"ak_merit_fac.zip"),
"64ff628ed34ed30c2057b430"))
- out <- list(merit_catchments = file.path(merit_dir,
+ out <- list(merit_catchments = file.path(merit_dir, "merged_AK_MERIT_Hydro",
"cat_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp"),
- merit_rivers = file.path(merit_dir,
+ merit_rivers = file.path(merit_dir, "merged_AK_MERIT_Hydro",
"riv_pfaf_78_81_82_MERIT_Hydro_v07_Basins_v01.shp"),
aster_dem = file.path(merit_dir, "dem.tif"),
merit_dem = file.path(merit_dir, "ak_merit_dem.tif"),
--
GitLab