USGSDataApp/server.R

In cahhansen/Shiny-WaDE-Data-Explorer: Shiny WaDE Data Explorer

# server.R
library(ggplot2)
library(XML)
library(RColorBrewer)
library(RCurl)
library(wadeR)
options(RCurlOptions = list(ssl.verifypeer = FALSE))
shinyServer(
  function(input, output,session) {
    # input <- list(reportingunit = "12040204-TX West Galveston Bay", year = "1990") # test input
    #REPORTING UNIT INFO##############################################################################################################################
    load("data/ReportingUnits.Rdata") # loads RU_df
    
    #CONSUMPTIVE USE###################################################################################################################################    
    #Create reactive function that will fetch consumptive use data when the user changes the inputs (year or location)
    CU_data <- reactive({
      
      wade_url <- paste0('https://wade-development.usgs.chs.ead/WADE/v0.2/',
                         'GetSummary/GetSummary.php?loctype=REPORTUNIT&loctxt=',
                         RU_df[(RU_df$Name_ID == input$reportingunit), "ReportingUnit"],
                         '&orgid=NWUSP&reportid=',
                         input$year,
                         '-CONSUMPTIVEUSE&datatype=ALL')
      
      return(get_wade_data(wade_url = wade_url))
      
    })
    
    #Returns information about the method
    xml.urlCUMethod=getURLContent(paste0('https://wade-development.usgs.chs.ead/WADE/v0.2/GetMethod/GetMethod.php?methodid=NWUSP&methodname=NOT%20KNOWN'))
    CUmethodroot= xmlRoot(xmlParse(xml.urlCUMethod))
    #Extract Report Summary
    CUmethodsummary=CUmethodroot[["Organization"]][["Method"]]
    CUmethodinfo=xmlSApply(CUmethodsummary,function(x) xmlSApply(x,xmlValue))
    CUmethod_df=data.frame(t(CUmethodinfo),row.names=NULL,stringsAsFactors = FALSE)
    
    #Plots the results from the GetSummary call
    output$CUplot <- renderPlot({
      #Get data from reactive function (based on inputs of year and reporting unit)
      CU_df=CU_data()
      title="Comparison of Water Use by Sector"
      ggplot(data=CU_df,environment = environment())+
        geom_bar(aes(x=Sector,y=Amount, fill=Sector),stat="identity")+
        theme_bw()+scale_fill_brewer(palette="Paired")+
        xlab("Sector")+ylab("Water Use (acre-feet/year)")+ggtitle(title)+
        theme(legend.position="none",plot.title = element_text(hjust = 0.5))
    })
    
    #List amounts in table
    output$CUtable <- renderTable({
      CU_df=CU_data()
    },bordered=TRUE, striped=TRUE)
    
    
    #Prints the method information
    output$CUmethod <- renderText(
        paste0("For more information about the methods used, see: ",CUmethod_df$MethodLinkText[[1]])
      )
    
    #DIVERSION##########################################################################################################################################    
    #Create reactive function that will fetch diversion data when the user changes the inputs (year or location)
    Div_data <- reactive({
      
      wade_url <- paste0('https://wade-development.usgs.chs.ead/WADE/v0.2/',
                        'GetSummary/GetSummary.php?loctype=REPORTUNIT&loctxt=',
                        RU_df[(RU_df$Name_ID==input$reportingunit),"ReportingUnit"],
                        '&orgid=NWUSP&reportid=',
                        input$year,
                        '-WITHDRAWALS&datatype=ALL')

      return(get_wade_data(wade_url = wade_url))
    })
    
    #Returns information about the method
    xml.urlDivMethod=getURLContent(paste0('https://wade-development.usgs.chs.ead/WADE/v0.2/GetMethod/GetMethod.php?methodid=NWUSP&methodname=NOT%20KNOWN'))
    Divmethodroot= xmlRoot(xmlParse(xml.urlDivMethod))
    #Extract Report Summary
    Divmethodsummary=Divmethodroot[["Organization"]][["Method"]]
    Divmethodinfo=xmlSApply(Divmethodsummary,function(x) xmlSApply(x,xmlValue))
    Divmethod_df=data.frame(t(Divmethodinfo),row.names=NULL,stringsAsFactors = FALSE)
    
    #Plots results from the GetSummary call
    output$Divplot <- renderPlot({
      Div_df=Div_data()
      title="Comparison of Diversions by Sector"
      ggplot(data=Div_df,environment = environment())+
        geom_bar(aes(x=Sector,y=Amount, fill=SourceType),stat="identity")+
        theme_bw()+scale_fill_brewer(palette="Paired", name="Source")+
        xlab("Sector")+ylab("Diversions (acre-feet/year)")+ggtitle(title)+
        theme(legend.position="bottom",plot.title = element_text(hjust = 0.5))
    })
    
    #List amounts in table
    output$Divtable <-renderTable({
      Div_df=Div_data()
    },bordered=TRUE,striped=TRUE)
    
    
    #Prints the method information
    output$Divmethod <- renderText(
      paste0("For more information about the methods used, see: ",Divmethod_df$MethodLinkText[[1]])
      )
  }
)