# NOT RUN {
sample_data <- system.file("extdata/sample_natseamless.gpkg",
package = "nhdplusTools")
nhdplus_path(sample_data)
staged_nhdplus <- stage_national_data(output_path = tempdir())
sample_flines <- readRDS(staged_nhdplus$flowline)
subset_rpu(sample_flines, rpu = "07b")
# }
Run the code above in your browser using DataLab