FedData is an R package implementing functions to automate downloading geospatial data available from several federated data sources (mainly sources maintained by the US Federal government).
FedData version 2.5 will be the final minor CRAN release of FedData 2. FedData 3 will be released in the coming months, but some code built on FedData 2 will not be compatible with FedData 3.
Currently, the package enables extraction from seven datasets:
This package is designed with the large-scale geographic information system (GIS) use-case in mind: cases where the use of dynamic web-services is impractical due to the scale (spatial and/or temporal) of analysis. It functions primarily as a means of downloading tiled or otherwise spatially-defined datasets; additionally, it can preprocess those datasets by extracting data within an area of interest (AoI), defined spatially. It relies heavily on the sp, raster, and rgdal packages.
This package has been built and tested on a source (Homebrew) install of R on macOS 10.14 (High Sierra), and has been successfully run on Ubuntu 14.04.5 LTS (Trusty), Ubuntu 16.04.1 LTS (Xenial) and binary installs of R on Mac OS 10.14 and Windows 10.
Development version from GitHub:
Linux (Ubuntu 14.04.5 or 16.04.1):
First, in terminal:
bash sudo add-apt-repository ppa:ubuntugis/ppa -y sudo apt-get update -q sudo apt-get install libssl-dev libcurl4-openssl-dev netcdf-bin libnetcdf-dev gdal-bin libgdal-dev Then, in R:
This demonstration script is available as an R Markdown document in the GitHub repository: https://github.com/ropensci/FedData.
# Get the DAYMET (North America only) # Returns a raster DAYMET <- get_daymet(template = vepPolygon, label = "VEPIIN", elements = c("prcp","tmax"), years = 1980:1985) # Plot with raster::plot raster::plot(DAYMET$tmax$X1985.10.23)
# Get the daily GHCN data (GLOBAL) # Returns a list: the first element is the spatial locations of stations, # and the second is a list of the stations and their daily data GHCN.prcp <- get_ghcn_daily(template = vepPolygon, label = "VEPIIN", elements = c('prcp')) # Plot the NED again raster::plot(NED) # Plot the spatial locations sp::plot(GHCN.prcp$spatial, pch = 1, add = TRUE) legend('bottomleft', pch = 1, legend="GHCN Precipitation Records")
# Elements for which you require the same data # (i.e., minimum and maximum temperature for the same days) # can be standardized using standardize==T GHCN.temp <- get_ghcn_daily(template = vepPolygon, label = "VEPIIN", elements = c('tmin','tmax'), years = 1980:1985, standardize = TRUE) # Plot the NED again raster::plot(NED) # Plot the spatial locations sp::plot(GHCN.temp$spatial, add = TRUE, pch = 1) legend('bottomleft', pch = 1, legend = "GHCN Temperature Records")
# Get the NRCS SSURGO data (USA ONLY) SSURGO.VEPIIN <- get_ssurgo(template = vepPolygon, label = "VEPIIN") # Plot the NED again raster::plot(NED) # Plot the SSURGO mapunit polygons plot(SSURGO.VEPIIN$spatial, lwd = 0.1, add = TRUE)
# Or, download by Soil Survey Area names SSURGO.areas <- get_ssurgo(template = c("CO670","CO075"), label = "CO_TEST") # Let's just look at spatial data for CO675 SSURGO.areas.CO675 <- SSURGO.areas$spatial[SSURGO.areas$spatial$AREASYMBOL=="CO075",] # And get the NED data under them for pretty plotting NED.CO675 <- get_ned(template = SSURGO.areas.CO675, label = "SSURGO_CO675") # Plot the SSURGO mapunit polygons, but only for CO675 plot(NED.CO675) plot(SSURGO.areas.CO675, lwd = 0.1, add = TRUE)
# Get the ITRDB records ITRDB <- get_itrdb(template = vepPolygon, label = "VEPIIN", recon.years = 850:2000, calib.years = 1924:1983, measurement.type = "Ring Width", chronology.type = "ARSTND") #> Warning in eval(jsub, SDenv, parent.frame()): NAs introduced by coercion #> Warning: attribute variables are assumed to be spatially constant #> throughout all geometries # Plot the NED again raster::plot(NED) # Map the locations of the tree ring chronologies plot(ITRDB$metadata$geometry, pch = 1, add = TRUE) legend('bottomleft', pch = 1, legend = "ITRDB chronologies")