GSIF (version 0.3-1)

edgeroi: The Edgeroi Data Set

Description

Soil samples and covariate layers for the Edgeroi area in NSW, Australia (ca 1500 square-km).

Usage

data(edgeroi)

Details

The Edgeroi is one of the standard soil data sets used to test soil mapping methods in Australia. Out of 359 profiles, 210 sites were sampled on a systematic, equilateral triangular grid with a spacing of 2.8 km between sites, the other sites are distributed more irregularly or on transects. The data set is described in detail in http://dx.doi.org/10.1016/j.geoderma.2009.10.007{Malone et al. (2010)} and McGarry et al. (1989). The edgeroi contains only a subset of the original http://www.asris.csiro.au/mapping/hyperdocs/NatSoil/{NatSoil} records. Observed soil classes for TAXGAUC are (alphabetically): Alluvial soil ("A"), Brown clay ("BC"), Black earth ("BE"), Earthy sand ("ES"), Grey clay ("GC"), Grey earth ("GE"), No suitable group ("NSG"), Prairie soil ("PS"), Rendzina ("R"), Red-brown earth ("RBE"), Red clay ("RC"), Red earth ("RE"), Red podzolic soil ("RP"), Solodic soil ("SC"), Soloth ("SH"), Solonchak ("SK"), Siliceous sand ("SS"), and Solonetz ("SZ").

References

• Malone, B.P., McBratney, A.B., Minasny, B. (2010)http://dx.doi.org/10.1016/j.geoderma.2009.10.007{Mapping continuous depth functions of soil carbon storage and available water capacity}. Geoderma 154, 138-152.
• McGarry, D., Ward, W.T., McBratney, A.B. (1989) Soil Studies in the Lower Namoi Valley: Methods and Data. The Edgeroi Data Set. (2 vols) (CSIRO Division of Soils: Adelaide).
• Minty, B., Franklin, R., Milligan, P., Richardson, L.M., and Wilford, J., (2009)http://dx.doi.org/10.1071/EG09025{The Radiometric Map of Australia}. Exploration Geophysics, 40(4), 325-333.

Examples

library(rgdal)
library(aqp)
library(sp)

data(edgeroi)
edgeroi$sites[edgeroi$sites$SOURCEID=="399_EDGEROI_ed095_1",]
edgeroi$horizons[edgeroi$horizons$SOURCEID=="399_EDGEROI_ed095_1",]
## spPoints:
sites <- edgeroi$sites
coordinates(sites) <- ~ LONGDA94 + LATGDA94
proj4string(sites) <- CRS("+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs")
## plot points and grids:
pnts <- list("sp.points", spTransform(sites, CRS("+init=epsg:28355")), pch="+", col="black")
## load the 250 m grids:
download.file("http://worldgrids.org/rda/edgeroi.grids.rda", "edgeroi.grids.rda")
load("edgeroi.grids.rda")
gridded(edgeroi.grids) <- ~x+y
proj4string(edgeroi.grids) <- CRS("+init=epsg:28355")
spplot(edgeroi.grids[1], sp.layout=pnts)
## load the 100 m grids:
download.file("http://worldgrids.org/rda/edgeroi.grids100.rda", "edgeroi.grids100.rda")
load("edgeroi.grids100.rda")
gridded(edgeroi.grids100) <- ~x+y
proj4string(edgeroi.grids100) <- CRS("+init=epsg:28355")
spplot(edgeroi.grids100["TI1LAN6"], sp.layout=pnts)