fifelse()
as reported in https://github.com/ropensci/GSODR/issues/121NA
when "STP_ATTRIBUTES" are equal to "0" or no observations used in reporting "STP" to see if this overcomes the issue.Added citation information for the data themselves to README.
Adds new function, get_updates()
, which gets the changelog for the GSOD data and returns it sorted with most recent changes first.
This is really a bug in the GSOD data, not {GSODR}, but we do our best to please.
As reported by @geospacedman in #117, STP values above 1000 were not properly reported.
Upon further inspection, it was found that the GSOD data in the CSV files is incorrect and the leading "1" is truncated and the values are reported as "034" for e.g. "1034".
Further, {data.table} drops the leading zeros by default on import with fread()
unless keepLeadingZeros = TRUE
is set, so that is now set and the 1 is appended when the values are >1000 and all should be well with the world again.
Fixes title that ended with a full stop.
nearest_stations()
now returns a data.table
of all station metadata and a value with the distance in kilometres from the user-provided coordinates.
This function previous returned a single character
vector.
To replicate the previous functionality, you can use the following method by calling the STNID
column name.nearest_stations(LAT = 14.16742, LON = 121.255669, distance = 50)$STNID
Use {roxyglobals} for handling global values.
Update internal isd history database.
Tidy up Appendix 2 in vignette, map of GSOD station locations. Add a caption and remove text from .png image itself.
Remove {httr} as an Import. Use base functionality to check existence first and then {curl} will properly error if a download fails, we will use that functionality.
Fix codecov badge to point at proper branch.
Ensure function titles are all proper title case.
Enhanced handling of data requests for incompatible station-year combinations. A warning()
is issued for each bad combo, but if a vector()
of stations is requested, any valid stations will still be downloaded with the bad combos omitted after the warning is emitted.
Fix bug where the isd_history.Rda
data object was not updated using the data_raw/fetch_isd-history.Rmd
.
Fix bug where stations were not available in the internal isd_history
file, this caused NA
values to appear for some stations when run using format_GSOD()
or to not be available for download with get_GSOD()
.
This adds 1,334 new stations that are now available through GSODR.
These were always available through GSOD but this bug prevented them from being accessed through this package.
Fix minor bug where get_inventory()
returned garbage in the header of the object.
Redoc package to correct HTML issues.
Update internal isd-history
database.
Update changed URLs.
Update internal isd-history
database.
Use return(invisible(NULL))
for functions that check user inputs without any returns.
Replace class(years) == "character")
with inherits(years, what = "character")
in an internal function, .validate_years()
.
Updates invalid URLs.
Skip ANY and ALL tests on CRAN. This fixes the "problems" with GSODR failing on a Solaris instance when the server failed to respond.
Update internal isd-history
database.
Use \CRANpkg{}
in place of \pkg{}
in documentation.
NA
value replacements.Simplify NA
value replacement in "PRCP" column.
The PRCP column values are rounded to two decimal places as in original GSOD data, not one.
The TEMP_ATTRIBUTES, DEWP_ATTRIBUTES, SLP_ATTRIBUTES, STP_ATTRIBUTES, VISIB_ATTRIBUTES and WDSP_ATTRIBUTES columns are formatted as an integer not character.
Better tests for the generated weather data data.table
output checking values and formats.
Tests are updated for updated data availability in the GSOD data due to continuous improvements to the data set.
Standardise handling of author/contributor comments. None have a full stop now in the comment.
Use on.exit()
to reset the working directory to the original user-space value after changing the working directory to untar files located in tempdir()
.
Fixes bug reported in #84 in the FRSHTT columns where the values were all reported as NA
even if there were observed values.
Fixes bug where NA values reported as 99.99, 999.9 or 9999.9 were not replaced with NA
.
Fix bug where FRSHTT (Fog, Rain/Drizzle, Snow/Ice, Hail, Tornado, Thunder) column values split into the respective columns only returned NA
, not the proper values as expected.
Reported in #84.
\donttest{}
but use @examplesIf interactive()
instead.Include columns for COUNTRY_NAME (country name in English), ISO2C and ISO3C in the final output from get_GSOD()
and reformat_GSOD()
.
Diffs in the isd_history are now recorded in the /data-raw/fetch_isd-history.md
file and shipped with GSODR as isd_history.rda
, which can be viewed by using load(system.file("extdata", "isd_diff.rda", package = "GSODR"))
.
Update and improve documentation to reflect country name and ISO code columns.
future.apply::future_lapply()
caused the get_GSOD()
and reformat_GSOD()
functions to run without completing or responding was fixed by simply using R's base lapply()
function. If parallel processing is needed, users should implement their own solutions to download and process in parallel.Fix bug that caused the package to run without responding.
Fix test that failed on CRAN's Solaris server for some reason.
Removes a working DOI link from the reference for the equation used because win-builder checks say it doesn't work (even though it does and there's nothing wrong with the link any time I check).
nearest_stations()
did not always return the nearest station as the first value in the vectorUpdate internal isd-history database, adding 11 stations
Fix any links that redirect found in DESCRIPTION, documentation or other materials in the package
Update internal station list to latest
Correct an error in documentation for update_station_list()
Remove spatial vignettes to slim down Suggests and make CI maintenance easier
Alduchov, O.A. and Eskridge, R.E., 1996. Improved Magnus form approximation of saturation vapor pressure. Journal of Applied Meteorology and Climatology, 35(4), pp.601-609.
Update internal station list to latest
Enhanced documentation
Corrects internal bug that provided a warning message when GSOD files were parsed
Fixes bug where not all files downloaded were cleaned up on the function exit when fetching station inventories
Fixes bug where station inventories from get_inventory()
lacked the location metadata, i.e. country and other geographic information
Update vignette to use latest functions from tidyr, i.e. tidyr::pivot_longer()
Update internal station list to latest
Tidy up documentation, mainly fix functions' title capitalisation
get_GSOD()
checks the number of stations being requested.
If the number is >10, the entire annual file will be downloaded and requested stations will then be selected and returned.
This saves time by reducing the number of requests made to the server.
Users should not see any difference other than quicker responses for a large number of requested stations.
get_GSOD()
now uses https rather than FTP server, correcting bug where the data could not be downloaded any longerCorrected elevation values are no longer available from GSODR
Objects are returned as data.table
objects
get_inventory()
now uses https rather than FTP server
update_station_list()
now uses https rather than FTP server
Implement better error handling when attempting to fetch station inventories
Reduced package dependencies
Improved vignettes that are pre-compiled for faster package installation and updated content with linting and error corrections
Users may now specify country by FIPS code when using get_GSOD()
Improved test coverage
Update internal database of station locations
future_apply
in processing files after downloading. This allows for end users to use a parallel process of their choice.Refactor internal functionality to be more clear and efficient in execution
country-list
is not loaded unless user has specified a country in get_GSOD()
An instance where the FIPS code was determined twice was removed
Replace \dontrun{}
with \donttest{}
in documentation examples
Ensure that DESCRIPTION file follows CRAN guidelines
Format help files, fixing errors and formatting for attractiveness
Update internal database of station locations
Store internal database of station locations fields BEGIN
and END
as integer, not double
Clarify code of conduct statement in README that it only applies to this, GSODR, project
Prompt user for input with warning about reproducibility if using the update_station_list()
function
Adds metadata header to the tibble
returned by get_inventory()
Remove start-up message to conform with rOpenSci guidelines
Remove extra code, clean up code-chunks and use hrbrthemes::theme_ipsum()
for data-raw/fetch_isd-history.md
Fix bug in creating isd-history.rda
file where duplicate stations existed in the file distributed with GSODR
but with different corrected elevation values
Repatch bug reported and fixed previously in version 1.2.0 where Windows users could not successfully download files. This somehow snuck back in.
Introduce a message if a station ID is requested but files are not found on the server. This is in response to an inquiry from John Paul Bigouette where a station is reported as having data in the inventory but the files do not exist on the server.
Fix bug that removed a few hundred stations from the internal GSODR
database of stations in the data-raw
files.
Clean documentation, shortening long lines, fixing formatting, incomplete sentences and broken links
Clarify the reasons for errors that a user may encounter
Update internal databases of station metadata
Clean up this file
Remove ability to export files from get_GSOD()
to slim down the package dependencies and this functions parameters.
Examples of how to convert to a spatial object (both sp and sf are shown) and export ESRI Shapefiles and GeoPackage files are now included in the vignette.
As a result of the previous point, the sp and rgdal packages are no longer Imports but are now in Suggests along with sf for examples in the GSOD vignette.
Fix a nasty bug where GSOD files downloaded using Windows would not untar properly.
This caused the get_GSOD()
function to fail.
Thanks to Ross Darnell, CSIRO, for reporting this.
Correct options in "GSODR use case: Specified years/stations vignette" on line 201 where file
was incorrectly used in place of path
.
Thanks to Ross Darnell, CSIRO, for reporting this.
Correct documentation for reformat_GSOD()
Update internal databases of station metadata
Vignettes contain pre-built figures for faster package installation when building vignettes
Fix start-up message formatting
Correct ORCID comment in author field of DESCRIPTION
Update internal databases for country list and isd_history
MAX_FLAG
and MIN_FLAG
columns now report NA
when there is no flagComment for Bob and Hugh in DESCRIPTION now only ORCID url
dplyr version set to >= 0.7.0 not 0.7 as before
Start-up message statement is more clear in relation to WMO resolution 40, that GSODR does not redistribute any weather data itself
Remove unnecessary function, .onLoad(), from zzz.R
Function titles in documentation now in title case
Correct grammar in documentation
The data.table and fields packages are no longer imported. All internal functions now use dplyr or base R functionality, reducing the dependencies of GSODR
Any data frames returned by GSODR functions are returned as a tibble()
object
The YEARMODA
column is now returned as Date
without time, rather than Character
Add new function, get_inventory()
, which downloads the NCEI's station inventory document and returns a tibble()
object of the data
Use larger images and provide a table of contents in vignettes
Updated and enhanced introductory vignette
Update internal stations list
Fix documentation in vignette where first example would not run due to changes in package data formats
Fix bug in GSODR vignette where examples would not run due to libraries not being loaded
Fix bug where prior server queries would be pre/appended to subsequent queries
Fix bug where invalid stations would return an empty data frame, should stop and return message about checking the station
value supplied to get_GSOD()
and check if data are available for the years requested
Update Appendix 2 of GSODR vignette, map of station locations, to be more clear and follow same format as that of bomrang
package
Update example output in GSODR vignette where applicable
Add welcome message on start-up regarding data use and sharing
Update internal stations list
Fix bug where "Error in read_connection_(con):" when writing to CSV occurs
Fix typo in line 160 of get_GSOD()
where "Rda" should be "rda" to properly load internal package files
Data distributed with GSODR are now internal to the package and not externally exposed to the user
Vignettes have been updated and improved with an improved order of information presented and some have been combined for easier use
Data for station locations and unique identifiers is now provided with the package on installation. Previously this was fetched each time from the ftp server.
The station metadata can now be updated if necessary by using update_station_list()
, this change overwrites the internal data that were originally distributed with the package.
This operation will fetch the latest list of stations and corresponding information from the NCEI ftp server.
Any changes will be overwritten when the R package is updated, however, the package update should have the same or newer data included, so this should not be an issue.
Replace plyr functions with purrr, plyr is no longer actively developed
plyr is no longer an import
get_station_list()
is no longer supported.
Instead use the new
update_station_list()
to update the package's internal station database.
Improved documentation (i.e., spelling corrections and more descriptive)
More descriptive vignette for "GSODR use case: Specified years/stations vignette"
Round MAX/MIN temp to one decimal place, not two
Update SRTM elevation data
Update country list data
Fix missing images in README.html on CRAN
Update documentation for get_GSOD()
when using station
parameter
Edit paper.md for submission to JOSS
Remove extra packages listed as dependencies that are no longer necessary
Correct Working_with_spatial_and_climate_data.Rmd where it was missing the first portion of documentation and thus examples did not work
The get_GSOD()
function returns a data.frame
object in the current R session with the option to save data to local disk
Multiple stations can be specified for download rather than just downloading a single station or all stations
A new function, nearest_stations()
is now included to find stations within a user specified radius (in kilometres) of a point given as latitude and longitude in decimal degrees
A general use vignette is now included
New vignette with a detailed use-case
Output files now include fields for State (US only) and Call (International Civil Aviation Organization (ICAO) Airport Code)
Use FIPS codes in place of ISO3c for file name and in output files because some stations do not have an ISO country code
Spatial file output is now in GeoPackage format (GPKG). This results in a single file output unlike shapefile and allows for long field names
Users can specify file name of output
R >= 3.2.0 now required
Field names in output files use "_" in place of "."
Long field names now used in file outputs
Country is specified using FIPS codes in file name and output file contents due to stations occurring in some locales that lack ISO 3166 3 letter country codes
The get_GSOD()
function will retrieve the latest station data from NCDC and automatically merge it with the CGIAR-CSI SRTM elevation values provided by this package.
Previously, the package provided it's own list of station information, which was difficult to keep up-to-date
A new reformat_GSOD()
function reformats station files in "WMO-WBAN-YYYY.op.gz" format that have been downloaded from the United States National Climatic Data Center's (NCDC) FTP server.
A new function, get_station_list()
allows for fetching latest station list from the FTP server and querying by the user for a specified station or location.
New data layers are provided through a separate package, GSODRdata
, which provide climate data formatted for use with GSODR.
CHELSA (climatic surfaces at 1 km resolution),
MODCF * Remotely sensed high-resolution global cloud dynamics for predicting ecosystem and biodiversity distributions,
ESACCI * ESA's CCI-LC snow cover probability and
CRU CL2.0 (climatic surfaces at 10 minute resolution).
Improved file handling for individual station downloads
Missing values are handled as NA
not -9999
Change from GPL >= 3 to MIT licence to bring into line with ropensci packages
Now included in ropensci, ropensci/GSODR
get_GSOD()
function optimised for speed as best possible after FTPing files from NCDC server
All files are downloaded from server and then locally processed, previously these were sequentially downloaded by year and then processed
A progress bar is now shown when processing files locally after downloading
Reduced package dependencies
The get_GSOD()
function now checks stations to see if the years being queried are provided and returns a message alerting user if the station and years requested are not available
When stations are specified for retrieval using the station = ""
parameter, the get_GSOD()
function now checks to see if the file exists on the server, if it does not, a message is returned and all other stations that have files are processed and returned in output
Documentation has been improved throughout package
Better testing of internal functions
Fixed: Remove redundant code in get_GSOD()
function
Fixed: The stations data frame distributed with the package now include stations that are located above 60 latitude and below -60 latitude
Missing values are reported as NA for use in R, not -9999 as previously
The path
parameter is now instead called dsn
to be more inline with other tools like readOGR()
and writeOGR()
Shapefile file out is no longer supported. Use GeoPackage (GPKG) instead
The option to remove stations with too many missing days is now optional, it now defaults to including all stations, the user must specify how many missing stations to check for an exclude
The max_missing
parameter is now user set, defaults to no check, return all stations regardless of missing days
Users can now select to merge output for station queries across multiple years.
Previously one year = one file per station.
Now are set by user, merge_station_years = TRUE
parameter, only one output file is generated
Country list is now included in the package to reduce run time necessary when querying for a specific country. However, this means any time that the country-list.txt file is updated, this package needs to be updated as well
Updated stations
list with latest version from NCDC published 12-07-2016
Country list is now included in the package to reduce run time necessary when querying for a specific country. However, this means any time that the country-list.txt file is updated, this package needs to be updated as well
Country level, agroclimatology and global data query conversions and calculations are processed in parallel now to reduce runtime
Improved documentation with spelling fixes, clarification and updates
Enable ByteCompile
option upon installation for small increase in speed
Use write.csv.raw
from [iotools]("https://cran.r-project.org/web/packages/iotools/index.html")
to greatly improve runtime by decreasing time used to write CSV files to disk
Use writeOGR()
from rgdal
, in place of raster's
shapefile
to improve runtime by decreasing time used to write shapefiles to disk
Country level, agroclimatology and global data query conversions and calculations are processed in parallel now to reduce runtime
Improved documentation with spelling fixes, clarification and updates
Enable ByteCompile
option upon installation for small increase in speed
Use write.csv.raw
from [iotools]("https://cran.r-project.org/web/packages/iotools/index.html")
to greatly improve runtime by decreasing time used to write CSV files to disk
Fix bug with connection timing out for single station queries.
Somehow the previously working function that checked country names broke with the toupper()
function.
A new function from juba fixes this issue and users can now select country again
User entered values for a single station are now checked against actual station values for validity
stations.rda is compressed
stations.rda now includes a field for "corrected" elevation using hole-filled SRTM data from Jarvis et al. 2008
Set NA or missing values in CSV or shapefile to -9999 from -9999.99 to align with other data sources such as WorldClim
Fix issues with MIN/MAX where MIN referred to MAX (Issue 5)
Fix bug where the tf
item was incorrectly set as tf < * "~/tmp/GSOD-2010.tar
, not tf < * tempfile
, in get_GSOD()
(Issue 6)
CITATION file is updated and corrected
User now has the ability to generate a shapefile as well as CSV file output (Issue 3)
Documentation is more complete and easier to use
Fix issue when reading .op files into R where temperature was incorrectly read causing negative values where T >= 100F, this issue caused RH values of >100% and incorrect TEMP values (Issue 1)
Spelling corrections
Include MIN/MAX flag column
Station data is now included in package rather than downloading from NCDC every time get_GSOD() is run, this data has some corrections where stations with missing LAT/LON values or elevation are omitted, this is not the original complete station list provided by NCDC.
Fixed bug where YDAY not correctly calculated and reported in CSV file
CSV files for station only queries now are names with the Station Identifier. Previously named same as global data
Likewise, CSV files for agroclimatology now are names with the Station Identifier. Previously named same as global data.
Set values where MIN > MAX to NA
Set more MIN/MAX/DEWP values to NA. GSOD README indicates that 999 indicates missing values in these columns, this does not appear to always be true. There are instances where 99 is the value recorded for missing data. While 99F is possible, the vast majority of these recorded values are missing data, thus the function now converts them to NA
Bug fix in MIN/MAX with flags. Some columns have differing widths, which caused a flag to be left attached to some values
Correct URL in README.md for CRAN to point to CRAN not GitHub
Bug fix in importing isd-history.csv file. Previous issues caused all lat/lon/elev values to be >0.
Bug fix where WDSP was mistyped as WDPS causing the creation of a new column, rather than the conversion of the existing
Bug fix if Agroclimatology selected. Previously this resulted in no records.
Set the default encoding to UTF8.
Bug fix for country selection. Some countries did not return proper ISO code.
Bug fix where WDSP was mistyped as WDPS causing the creation of a new column, rather than the conversion of the existing
Use write.csv, not readr::write_csv due to issue converting double to string: https://github.com/tidyverse/readr/issues/387
Now available on CRAN
Add single quotes around possibly misspelled words and spell out comma-separated values and geographic information system rather than just using "CSV" or "GIS" in DESCRIPTION.
Add full name of GSOD (Global Surface Summary of the Day) and URL for GSOD, https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ncdc:C00516 to DESCRIPTION as requested by CRAN.
Require user to specify directory for resulting .csv file output so that any files written to disk are interactive and with user's permission