Skip to content
Snippets Groups Projects
readNWISqw.Rd 5.23 KiB
Newer Older
  • Learn to ignore specific revisions
  • % Generated by roxygen2 (4.1.0): do not edit by hand
    % Please edit documentation in R/readNWISqw.r
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \name{readNWISqw}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \alias{readNWISqw}
    
    \title{Raw Data Import for USGS NWIS QW Data}
    \usage{
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    readNWISqw(siteNumbers, parameterCd, startDate = "", endDate = "",
    
      expanded = TRUE, reshape = FALSE, tz = "")
    
    \item{siteNumbers}{character of USGS site numbers.  This is usually an 8 digit number}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \item{parameterCd}{character that contains the code for a parameter
    group, or a character vector of 5-digit parameter codes. See \bold{Details}.}
    
    \item{startDate}{character starting date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
    retrieval for the earliest possible record.}
    
    \item{endDate}{character ending date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
    retrieval for the latest possible record.}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \item{expanded}{logical defaults to \code{TRUE}. If \code{TRUE}, retrieves additional information. Expanded data includes
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    remark_cd (remark code), result_va (result value), val_qual_tx (result value qualifier code), meth_cd (method code),
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    dqi_cd (data-quality indicator code), rpt_lev_va (reporting level), and rpt_lev_cd (reporting level type). If \code{FALSE},
    only returns remark_cd (remark code) and result_va (result value). Expanded = \code{FALSE} will not give
    sufficient information for unbiased statistical analysis.}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \item{reshape}{logical, reshape the expanded data. If \code{TRUE}, then return a wide data frame with all water-quality in a single row for each sample.
    If \code{FALSE} (default), then return a long data frame with each water-quality result in a single row. This
    argument is only applicable to expanded data. Data requested using \code{expanded=FALSE} is always returned in the wide format.}
    
    \item{tz}{character to set timezone attribute of datetime. Default is an empty quote, which converts the
    
    datetimes to UTC (properly accounting for daylight savings times based on the data's provided tz_cd column).
    Possible values to provide are "America/New_York","America/Chicago", "America/Denver","America/Los_Angeles",
    "America/Anchorage","America/Honolulu","America/Jamaica","America/Managua","America/Phoenix", and "America/Metlakatla"}
    
    A data frame with the following columns:
    \tabular{lll}{
    Name \tab Type \tab Description \cr
    agency_cd \tab character \tab The NWIS code for the agency reporting the data\cr
    site_no \tab character \tab The USGS site number \cr
    datetime \tab POSIXct \tab The date and time of the value converted to UTC (if asDateTime = TRUE), \cr
    \tab character \tab or raw character string (if asDateTime = FALSE) \cr
    tz_cd \tab character \tab The time zone code for datetime \cr
    code \tab character \tab Any codes that qualify the corresponding value\cr
    value \tab numeric \tab The numeric value for the parameter \cr
    }
    
    There are also several useful attributes attached to the data frame:
    \tabular{lll}{
    Name \tab Type \tab Description \cr
    url \tab character \tab The url used to generate the data \cr
    queryTime \tab POSIXct \tab The time the data was returned \cr
    comment \tab character \tab Header comments from the RDB file \cr
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    siteInfo \tab data.frame \tab A data frame containing information on the requested sites \cr
    variableInfo \tab data.frame \tab A data frame containing information on the requested parameters \cr
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    Imports data from NWIS web service. This function gets the data from here: \url{http://nwis.waterdata.usgs.gov/nwis/qwdata}
    A list of parameter codes can be found here: \url{http://nwis.waterdata.usgs.gov/nwis/pmcodes/}
    A list of statistic codes can be found here: \url{http://nwis.waterdata.usgs.gov/nwis/help/?read_file=stat&format=table}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \details{
    Valid parameter code groups are "All," or group codes:
    \tabular{ll}{
    Code \tab Description\cr
    INF \tab Information \cr
    PHY \tab Physical \cr
    INM \tab Inorganics, Major, Metals (major cations) \cr
    INN \tab Inorganics, Major, Non-metals (major anions) \cr
    NUT \tab Nutrient \cr
    MBI \tab Microbiological \cr
    BIO \tab Biological \cr
    IMN \tab Inorganics, Minor, Non-metals \cr
    IMM \tab Inorganics, Minor, Metals \cr
    TOX \tab Toxicity \cr
    OPE \tab Organics, pesticide \cr
    OPC \tab Organics, PCBs \cr
    OOT \tab Organics, other \cr
    RAD \tab Radiochemical \cr
    SED \tab Sediment \cr
    POP \tab Population/community \cr
    }
    
    If more than one parameter group is requested, only sites that data for all requested groups are returned.
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    }
    
    siteNumbers <- c('04024430','04024000')
    
    startDate <- '2010-01-01'
    endDate <- ''
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    parameterCd <- c('34247','30234','32104','34220')
    
    \dontrun{
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    rawNWISqwData <- readNWISqw(siteNumbers,parameterCd,startDate,endDate)
    rawNWISqwDataReshaped <- readNWISqw(siteNumbers,parameterCd,
    
              startDate,endDate,reshape=TRUE)
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    parameterCd <- "all"
    rawNWISall <- readNWISqw(siteNumbers,parameterCd,
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
              startDate,endDate)
    pgroup <- c("NUT")
    rawNWISNutrients <- readNWISqw(siteNumbers,pgroup,
              startDate,endDate)
    groups <- c("NUT","OPE")
    rawNWISNutOpe <- readNWISqw(siteNumbers,groups,
              startDate,endDate)
    rawNWISOpe <- readNWISqw(siteNumbers,"OPE",
              startDate,endDate)
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \seealso{
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    \code{\link{readWQPdata}}, \code{\link{whatWQPsites}},
    \code{\link{readWQPqw}}, \code{\link{constructNWISURL}}
    
    Laura A DeCicco's avatar
    Laura A DeCicco committed
    }
    
    \keyword{data}
    \keyword{import}
    \keyword{service}
    \keyword{web}