This function automates downloading, cleaning, reformatting of data from the Global Surface Summary of the Day (GSOD) data provided by the US National Centers for Environmental Information (NCEI), Three additional useful elements: saturation vapour pressure (es), actual vapour pressure (ea) and relative humidity (RH) are calculated and returned in the final data frame.
Parallel processing can be enabled using plan to set
up a parallel backend of your choice, e.g.,
future::plan("multisession").
get_GSOD(years, station = NULL, country = NULL, max_missing = NULL,
agroclimatology = FALSE)Year(s) of weather data to download.
Optional. Specify a station or multiple stations for which to retrieve, check and clean weather data using STATION. The NCEI reports years for which the data are available. This function checks against these years. However, not all cases are properly documented and in some cases files may not exist for download even though it is indicated that data was recorded for the station for a particular year. If a station is specified that does not have an existing file on the server, this function will silently fail and move on to existing files for download and cleaning.
Optional. Specify a country for which to retrieve weather data; full name, 2 or 3 letter ISO or 2 letter FIPS codes can be used. All stations within the specified country will be returned.
Optional. The maximum number of days allowed to be missing from a station's data before it is excluded from final file output.
Optional. Logical. Only clean data for stations
between latitudes 60 and -60 for agroclimatology work, defaults to FALSE.
Set to TRUE to include only stations within the confines of these
latitudes.
A data frame as a data.table object of
GSOD weather data.
All units are converted to International System of Units (SI), e.g., Fahrenheit to Celsius and inches to millimetres.
Data summarise each year by station, which include vapour pressure and relative humidity elements calculated from existing data in GSOD.
All missing values in resulting files are represented as NA regardless
of which field they occur in.
For a complete list of the fields and description of the contents and units,
please refer to Appendix 1 in the GSODR vignette,
vignette("GSODR", package = "GSODR").
For more information see the description of the data provided by NCEI, http://www7.ncdc.noaa.gov/CDO/GSOD_DESC.txt.
# NOT RUN {
# Download weather station for Toowoomba, Queensland for 2010
tbar <- get_GSOD(years = 2010, station = "955510-99999")
tbar
# Download data for Australia from 2010 to 2011
AUS <- get_GSOD(years = 2010:2011, country = "Australia")
AUS
# }
Run the code above in your browser using DataLab