From 2ddad90d529f20b7fc23b0ac9884914b17887932 Mon Sep 17 00:00:00 2001 From: Ben Varick Date: Mon, 16 Oct 2023 12:53:23 -0500 Subject: [PATCH] overwrite geopackage when combining data --- 01_download_data.R | 18 ++++++++++++++---- 02_process_data.R | 11 ++--------- extent.csv | 4 ++-- 3 files changed, 18 insertions(+), 15 deletions(-) diff --git a/01_download_data.R b/01_download_data.R index 1e2fc02..38ce473 100644 --- a/01_download_data.R +++ b/01_download_data.R @@ -2,6 +2,7 @@ library(nhdplusTools) library(sf) +# set data dir nhdplusTools_data_dir(dir = "data") # load extent of map @@ -28,9 +29,18 @@ extent <- list(longitude_max = max(extent$longitude), extent_huc <- get_huc(AOI = extent_poly, buffer = 0, type = "huc04") # download data -download_nhdplushr("data", extent_huc$huc4, download_files = TRUE) +download_nhdplushr(nhdplusTools_data_dir(), extent_huc$huc4, download_files = TRUE) + +# get vaa +get_vaa() # combine data into geopackage -data <- get_nhdplushr(hr_dir = "data", - out_gpkg = "data/data.gpkg", - layers = NULL) \ No newline at end of file +layers <- c("NHDArea", + "NHDFlowline", + "NHDWaterbody", + "NHDPlusLandSea") + +data <- get_nhdplushr(hr_dir = nhdplusTools_data_dir(), + out_gpkg = paste0(nhdplusTools_data_dir(),"/data.gpkg"), + layers = layers, + overwrite = TRUE) diff --git a/02_process_data.R b/02_process_data.R index 8c877cc..9f8a990 100644 --- a/02_process_data.R +++ b/02_process_data.R @@ -1,16 +1,9 @@ # load libraries library(sf) -# load data - - -layers <- c("NHDArea", - "NHDFlowline", - "NHDWaterbody", - "NHDPlusLandSea") - +# load and crop data data <- list(NULL) sf_use_s2(FALSE) for (layer in layers){ - data[[layer]] <- st_crop(st_read("data/data.gpkg", layer = layer), y = extent_bbox) + data[[layer]] <- st_crop(st_read(paste0(nhdplusTools_data_dir(),"/data.gpkg"), layer = layer), y = extent_bbox) } diff --git a/extent.csv b/extent.csv index 2c3d01e..f693b3c 100644 --- a/extent.csv +++ b/extent.csv @@ -1,3 +1,3 @@ latitude,longitude -43.63839,-89.97779 -42.52677,-88.52104 +47.33402,-93.10419 +42.43754,-86.832