Develop method to retrieve small sub-watersheds that are custom in VAHydro in order to add FTABLEs for them.

[rburghol]

Create batch script to create all subshed ftables :

• [x] Retrieve all from VAHydro
• [x] Retrieve model data from VAHydro for a single segment with REST
• [x] Is there a variable that marks it as a sub-watershed? (hydrocode DOES, but this is imperfect, tho maybe criteria "has 3 underscores in hydrocode?)
• [x] Calls simple command line based R script ftable_creation.R #331 to call to create a new FTABLE with arguments describing province, DA, and maybe lat/lon? (also river segment name of course), destination folder, model version (cbp-5.3.2 or cbp-6.0)
  • [x] Save FTABLEs to model scenario folder in CBP model tree (get with @glenncampagna and @juliabruneau who have been saving a lot of files)
  • [x] Estimate floodplain from VAHydro data
  • See: https://github.com/HARPgroup/HARParchive/blob/master/HARP-2022-Summer/AutomatedScripts/ftables/ftable_creation.R

[jdkleiner]

@megpritch See example method below for retrieving all rsegs (next step is to filter those down to just the sub-watersheds)

library("hydrotools")
basepath='/var/www/R';
source("/var/www/R/config.R")

##################################################################
# retrieve rsegs
##################################################################
# set up our data source environment
ds <- RomDataSource$new(site, rest_uname = rest_uname)
ds$get_token(rest_pw)

# retrieve all vahydro rseg features (this method is optimal b/c it allows retrieval of >100 records)
rsegs <- RomFeature$new(ds,list(ftype='vahydro',bundle='watershed'),TRUE)

# generate dataframe of rseg hydrocodes
rsegs_hydrocode_df <- data.frame("hydrocode"=character(),stringsAsFactors=FALSE)
for (i in 1:length(rsegs$hydrocode)){
  rsegs_hydrocode_df <- rbind(rsegs_hydrocode_df,rsegs$hydrocode[i])
}
print(head(rsegs_hydrocode_df))
##################################################################

##################################################################
# Next, use sqldf to query the sub-watersheds within the dataframe
# of hydrocodes (those with hydrocode greater than 29 characters)
##################################################################

# sqldf() here 

##################################################################

[megpritch]

Revised:

rsegs <- RomFeature$new(ds,list(ftype='vahydro',bundle='watershed'),TRUE)
  subshed<- subset(rsegs$hydrocode, nchar(rsegs$hydrocode)>29)
  subshed <- data.frame(subshed)