R/coordlist2shape.R

In jschap1/routR: Prepares inputs for the UW river routing model

#' coordlist2shape
#' 
#' Make shapefile from list of coordinates
#' @param gages data frame containing lat and lon values for each station
#' @param outname savename for the output shapefile
#' @param r raster or vector object with the desired output coordinate system
#' @export

coordlist2shape <- function(gages, outname, r)
{
  raster::coordinates(gages) <- ~lon+lat
  raster::crs(gages) <- raster::crs(r)
  rgdal::writeOGR(gages, 
                  dsn = outname, 
                  layer = "gages", 
                  driver = "ESRI Shapefile", 
                  overwrite = TRUE
                  )
  print(paste("Saved", outname))
  return(gages)
}