#install.packages("devtools") #devtools::install_github("anguswg-ucsb/cdssr") library(cdssr) library(tidyverse) library(janitor) #install.packages("janitor") #API <- 'PwevUQJCStcYZfqrOYbuyztmNPlUJWby' STR <- get_structures(div=3) %>% as_tibble %>% clean_names SBD_STR_DATA <- rbind(read_csv("Data/SBD_Data/StructureList_SBD1.csv") %>% mutate(SBD=1), read_csv("Data/SBD_Data/StructureList_SBD2.csv") %>% mutate(SBD=2), read_csv("Data/SBD_Data/StructureList_SBD3.csv") %>% mutate(SBD=3), read_csv("Data/SBD_Data/StructureList_SBD4.csv") %>% mutate(SBD=4), read_csv("Data/SBD_Data/StructureList_SBD5.csv") %>% mutate(SBD=5), read_csv("Data/SBD_Data/StructureList_SBD6.csv") %>% mutate(SBD=6)) %>% clean_names %>% mutate(wdid=as.character(wdid)) SBD_LINK <- SBD_STR_DATA %>% select(wdid,sbd) nrow(SBD_STR_DATA ) STR <- STR %>% left_join(SBD_LINK) %>% mutate(sbd=ifelse(is.na(sbd),"none",sbd)) STR <- STR %>% select(wdid,sbd,ciu_code,county,water_district,sbd,structure_type,management_district_name,latdecdeg,longdecdeg) STR ?get_structures_divrec_ts get_structures_divrec_ts(wdid='2005007') WDID <- 2005007 LINK <- paste0('https://dwr.state.co.us/Rest/GET/api/v2/structures/divrec/divrecyear/?format=tsv&dateFormat=dateOnly&fields=wdid%2CwaterClassNum%2CwcIdentifier%2CdataMeasDate%2CdataValue%2CmeasUnits%2CapprovalStatus&measInterval=Annual&wdid=',WDID,'&apiKey=PwevUQJCStcYZfqrOYbuyztmNPlUJWby') TEMP <- read_tsv(LINK,skip=2 ) TEMP %>% print(n=20) TOTAL_ROWS <- str_detect(TEMP$wcIdentifier,"Total") Total_Values <- TEMP[TOTAL_ROWS,] Total_Values %>% select(wdid,year=dataMeasDate,Volume=dataValue) OTHER <- TEMP[!TOTAL_ROWS,] OTHER$wcIdentifier TEMP %>% filter((wcIdentifier=='Total (Diversion)') %>% filter( STR %>% filter(structure_type=='WELL') get_parceluse_ts(wdid='2005000') STR <- as_tibble(STR)