R/sites2sp.R

In wasquith-usgs/akqdecay: Many Distributional Statistics of USGS Daily Streamflow Hydrographs and Powerful Utilities

"sites_to_SpatialPointsDataFrame" <-
function(siteNumbers, latlongcrs="+init=epsg:4269", proj4string=NA, usesp=TRUE, silent=FALSE, ...) {
   sitefile <- NULL; haveFirst <- TRUE; i <- 0; n <- length(siteNumbers)
   badsites <- vector("character")
   for(site in siteNumbers) {
      i <- i + 1
      if(! silent) message(" dataRetrieval::readNWISsite() for ",site," (",i,"/",n,")")
      asitefile <- NULL
      try(asitefile <- dataRetrieval::readNWISsite(site), silent=TRUE)
      if(is.null(asitefile)) {
         warning(" site ",site," does not 'exist' in NWISWeb")
         badsites <- c(badsites, site)
         next;
      }
      if(haveFirst) {
         sitefile <- asitefile
         haveFirst <- FALSE
      } else {
         sitefile <- rbind(sitefile, asitefile)
      }
   }
   if(! usesp) {
      attr(sitefile, "bad sites") <- badsites
      return(sitefile)
   }

   LATLONG <- sp::CRS(latlongcrs)

   sitefile$CDA <- pmin(sitefile$drain_area_va,
                        sitefile$contrib_drain_area_va, na.rm=TRUE)

   sitefile$inventory_dt    <- NULL
   sitefile$instruments_cd  <- NULL
   sitefile$construction_dt <- NULL
   sitefile$gw_file_cd      <- NULL
   sitefile$nat_aqfr_cd     <- NULL
   sitefile$aqfr_cd         <- NULL
   sitefile$aqfr_type_cd    <- NULL
   sitefile$well_depth_va   <- NULL
   sitefile$hole_depth_va   <- NULL
   sitefile$depth_src_cd    <- NULL
   sitefile$project_no      <- NULL
   sitefile$basin_cd        <- NULL
   sitefile$topo_cd         <- NULL

   # The sp library is needed here
   coords  <- cbind(sitefile$dec_long_va, sitefile$dec_lat_va)
   spSites <- sp::SpatialPointsDataFrame(coords, sitefile, proj4string=LATLONG)

   if(is.na(proj4string)) {
      ALBEA <- "+init=epsg:5070"
      message("**unprojected data are being returned**\n",
              "**consider using the 'proj4string' argument with a CRS, which could look\n",
              "**something like this:\n'",
              ALBEA,"'\n")
      return(spSites)
   } else {
      return(sp::spTransform(spSites, sp::CRS(proj4string)))
   }
}