examples/rest_get-set.R

In HARPgroup/hydro-tools: Data processing and accessors for VAHydro

library(httr);
save_directory <- "/var/www/html/files/fe/plots"
#----------------------------------------------
site <- "http://deq2.bse.vt.edu/d.bet"    #Specify the site of interest, either d.bet OR d.dh
#----------------------------------------------
# Load Libraries
basepath='/var/www/R';
source(paste(basepath,'config.local.private',sep='/'));
source(paste(hydro_tools,"VAHydro-2.0/rest_functions.R", sep = "/")); 
source(paste(hydro_tools,"VAHydro-1.0/fn_vahydro-1.0.R", sep = "/"));  
source(paste(hydro_tools,"LowFlow/fn_iha.R", sep = "/"));  
#retrieve rest token
source(paste(hydro_tools,"auth.private", sep = "/"));#load rest username and password, contained in auth.private file
token <- rest_token(site, token, rest_uname, rest_pw);
options(timeout=120); # set timeout to twice default level to avoid abort due to high traffic

riverseg = "OR2_8130_7900";
hydrocode = paste("vahydrosw_wshed_",riverseg,sep="");
ftype = 'vahydro'; # nhd_huc8, nhd_huc10, vahydro
inputs <- list (
  hydrocode = hydrocode,
  bundle = 'watershed',
  ftype = 'vahydro'
)
#property dataframe returned
feature = FALSE;
odata <- getFeature(inputs, token, site, feature);
hydroid <- odata[1,"hydroid"];
fname <- as.character(odata[1,]$name );
print(paste("Retrieved hydroid",hydroid,"for", fname,riverseg, sep=' '));

# get the p5.3.2, scenario  model segment attached to this river feature
inputs <- list(
  varkey = "om_model_element",
  featureid = hydroid,
  entity_type = "dh_feature",
  propcode = "p532cal_062211"
)
model <- getProperty(inputs, site, model)

# now, retrieve august low flow property if set
alfinfo <- list(
  varkey = "auglowflow",
  featureid = as.integer(as.character(model$pid)),
  entity_type = "dh_properties"
)
alfprop <- getProperty(alfinfo, site, alfprop)
# this just sets our property to the 
if (identical(alfprop, FALSE)) {
  # create
  alfprop = alfinfo
}
alfprop$propname = "August Low Flow"
alfprop$propvalue = 155
alfinfo$startdate = format(as.POSIXlt('1984-01-01'),"%s") 
alfinfo$enddate = format(as.POSIXlt('2005-12-31'),"%s") 
postProperty(alfprop,fxn_locations,base_url = site,alfprop) ;