R/ecor_map.R

In USEPA/LakeTrophicModelling: Package to reproduce Hollister et al. (2015) Modeling Lake Trophic State: A Randrom Forest Approach.

#' Create map of Ecoregions
#' 
#' Function to generate Map 1 in hkm2014ESA poster
#' 
#' @param states data for states, as data.frame
#' @param lakes point locations for lake samples as data.frame
#' @param myColor vector of length 3 with colors for fill, lines, and points in that order
#' 
#' @examples
#' wsa9<-readOGR("../inst/extdata/","wsa9_low48")
#' ecor_map(wsa9)
#' @export
#' @import ggplot2
ecor_map<-function(ecor){
  #borrowed from: https://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles#workaround-forced-layering-of-geom_plot
  ecor@data$id <- rownames(ecor@data)
  ecor_f <- fortify(ecor, region="WSA_9")
  names(ecor_f)[which(names(ecor_f)=="id")] = "WSA_9"
  ecor_df <- plyr::join(ecor_f, ecor@data, by="WSA_9")
  gmap<-ggplot(ecor_df, aes(long,lat,group=group,fill=WSA_9)) +
    geom_polygon(data=subset(ecor_df,WSA_9!="Western\nMountains (WMT)")) +
    geom_polygon(data=subset(ecor_df,WSA_9=="Western\nMountains (WMT)")) +
    geom_path(color="lightgrey") +
    coord_equal() +
    coord_map("albers", lat2 = 45.5, lat1 = 29.5)+
    theme(panel.background = element_rect(fill="white"), panel.grid = element_blank(),
          panel.border = element_blank(), 
          axis.text = element_blank(),axis.ticks = element_blank(),
          legend.text = element_text(size=10),
          legend.position = c(0.96,0.5),
          #legend.justification = c(0,0),
          legend.key.width = unit(0.5, "line"),
          legend.key.height = unit(1.6, "line"),
          legend.margin = unit(0,"line"),
          plot.margin = unit(c(1,3,0.5,0.5),"line")
          ) +
    ylab("") +
    xlab("") +
    scale_fill_manual(name="",values=viridis(10))
    
  return(gmap)
}