R/Al_default_criteria.R

In TravisPritchardODEQ/deqalcrit: Implementation of Oregon DEQ's Aluminum Criteria Procedures

#' Aluminum default criteria Lookup
#'
#' This function looks up default acute or chronic Al criteria values. The function 
#' sends lat/long values via a query to a Oregon DEQ map server to determine deafult DOC values. DOC values are based on 
#' georegion as described in ##insert link to interpretation doc when finalized##. If -9999 is returned, station is 
#' outside of georeferenced area.
#'
#' @param lat Latitude of station
#' @param long Longitude of station
#' @param type If 'Chronic' return CCC, if 'Acute' return CMC 
#' @return default criteria value. If NA is returned, station is outside of georeferenced area.
#' @examples
#' \dontrun{
#'al_default_criteria(45.57770, -122.7476, type = "Acute")
#'al_default_criteria(45.57770, -122.7476, type = "Chronic")
#'
#'#Return chronic criteria
#'aluminum_data_mlocs <- aluminum_data %>%
#'  select(MLocID, Lat_DD, Long_DD) %>%
#'  distinct() %>%
#'  rowwise() %>%
#'  mutate(def_DOC = Al_default_criteria(Lat_DD, Long_DD))
#'}
#'@export


al_default_criteria <- function(lat, long, type = c('Chronic')){
  
  
  if(type == 'Chronic'){
  
  
  path_region <- paste0("https://arcgis.deq.state.or.us/arcgis/rest/services/WQ/OR_Ecoregions_Aluminum/MapServer/0/query?geometry=",
                        long,"%2C+",
                        lat,
                        "&geometryType=esriGeometryPoint&inSR=4269&spatialRel=esriSpatialRelWithin&outFields=Def_CCC+&returnGeometry=false&returnTrueCurves=false&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&historicMoment=&returnDistinctValues=false&resultOffset=&resultRecordCount=&queryByDistance=&returnExtentOnly=false&datumTransformation=&parameterValues=&rangeValues=&quantizationParameters=&featureEncoding=esriDefault&f=geojson")
  
  request_region <- httr::GET(url = path_region)
  
  response_region <- httr::content(request_region, as = "text", encoding = "UTF-8")
  
  region_df<- geojsonsf::geojson_sf(response_region) %>%
    sf::st_drop_geometry()
  
  if(is.null(region_df[1,1])) {
    print("Lat/Long out of area")
    return(NA_real_)
  } else {
    
    return(region_df[1,1])
  }
  
  } else if(type == 'Acute'){
    
    
    path_region <- paste0("https://arcgis.deq.state.or.us/arcgis/rest/services/WQ/OR_Ecoregions_Aluminum/MapServer/0/query?geometry=",
                          long,"%2C+",
                          lat,
                          "&geometryType=esriGeometryPoint&inSR=4269&spatialRel=esriSpatialRelWithin&outFields=Def_CMC+&returnGeometry=false&returnTrueCurves=false&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&historicMoment=&returnDistinctValues=false&resultOffset=&resultRecordCount=&queryByDistance=&returnExtentOnly=false&datumTransformation=&parameterValues=&rangeValues=&quantizationParameters=&featureEncoding=esriDefault&f=geojson")
    
    request_region <- httr::GET(url = path_region)
    
    response_region <- httr::content(request_region, as = "text", encoding = "UTF-8")
    
    region_df<- geojsonsf::geojson_sf(response_region) %>%
      sf::st_drop_geometry()
    
    if(is.null(region_df[1,1])) {
      print("Lat/Long out of area")
      return(NA_real_)
    } else {
      
      return(region_df[1,1])
    }
    
  } else if (type == 'Ecoregion') {
    
    path_region <- paste0("https://arcgis.deq.state.or.us/arcgis/rest/services/WQ/OR_Ecoregions_Aluminum/MapServer/0/query?geometry=",
                          long,"%2C+",
                          lat,
                          "&geometryType=esriGeometryPoint&inSR=4269&spatialRel=esriSpatialRelWithin&outFields=OR_Ecoregi+&returnGeometry=false&returnTrueCurves=false&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&historicMoment=&returnDistinctValues=false&resultOffset=&resultRecordCount=&queryByDistance=&returnExtentOnly=false&datumTransformation=&parameterValues=&rangeValues=&quantizationParameters=&featureEncoding=esriDefault&f=geojson")
    
    request_region <- httr::GET(url = path_region)
    
    response_region <- httr::content(request_region, as = "text", encoding = "UTF-8")
    
    region_df<- geojsonsf::geojson_sf(response_region) %>%
      sf::st_drop_geometry()
    
    if(is.null(region_df[1,1])) {
      print("Lat/Long out of area")
      return(NA_character_)
    } else {
      
      return(region_df[1,1])
    }
  } 
  
}