Download the A-soil geochemistry data of the Tellus survey (Northern Ireland), and if desired produce a training image of the geochemical subcomposition Mg-Al-Ca-Fe-Rest.
getTellus(wd = ".", destfile = "TellusASoil.RData", TI = FALSE, cleanup = TRUE)
working directory, where intermediate files and output will be produced
file name where to put the Tellus data (including the extension ".Rdata"!)
either a logical (should the training image be created? defaults to FALSE) or else the file name where to put the created training image (including the extension ".Rdata")
should the downloaded excel file be removed? defaults to TRUE
TRUE if everything went OK. The function is actually called for the
side effect of downloading data from the Tellus survey project website, so be careful
to call it specifying the right directory at wd
.
The function will download an excel file to your working directory, do some
manipulations, save the resulting data.frame in a "TellusASoil.RData" file (or any other
name you provide on the destfile
argument) and
eventually, remove the excel file (if you left cleanup=TRUE
). This data set has
6862 obervations and 58 variables:
Sample ID
X coordinate of each point
Y coordinate of each point
a flag marking some data as a study subset, fixed for reproducibility
Silver concentration in ppm
Cadmium concentration in ppm
Indium concentration in ppm
Tin concentration in ppm
Antimony concentration in ppm
Tellurium concentration in ppm
Indim concentration in ppm
Caesium concentration in ppm
Barium concentration in ppm
Lanthanum concentration in ppm
Cerium concentration in ppm
Sodium oxide in percent
Magnesium oxide in percent
Aluminium oxide in percent
Silicium oxide in percent
Phosphorous(V) oxide in percent
Sulphur(III) oxide in percent
Potasium oxide in percent
Calcium oxide in percent
Titanium oxide in percent
Manganese oxide in percent
Total Iron(III) oxide in percent
Chlorine concentration in ppm
Scandium concentration in ppm
Vanadium concentration in ppm
Chromium concentration in ppm
Cobalt concentration in ppm
Nickel concentration in ppm
Copper concentration in ppm
Zinc concentration in ppm
Callium concentration in ppm
Germanium concentration in ppm
Arsenic concentration in ppm
Selenium concentration in ppm
Bromine concentration in ppm
Rubidium concentration in ppm
Strontium concentration in ppm
Yttrium concentration in ppm
Zirconium concentration in ppm
Niobium concentration in ppm
Molybdenum concentration in ppm
Neodymium concentration in ppm
Samarium concentration in ppm
Ytterbium concentration in ppm
Hafnium concentration in ppm
Tantalum concentration in ppm
Tungsten concentration in ppm
Thallium concentration in ppm
Lead concentration in ppm
Bismuth concentration in ppm
Thorium concentration in ppm
Uranium concentration in ppm
soil acidity
loss on ignition in percent
Additionally, if you give TI!=FALSE
,
the function will produce an additional
file "Tellus_TI.RData" (if TI=TRUE
, or any other filename that you
specify on the argument TI
) with a data.frame with 13287 rows and 8 columns:
X coordinate of each cell point
Y coordinate of each cell point
Magnesia proportion
Alumina proportion
Calcium oxide proportion
Iron oxide proportion
Residual complementing the sum to 1
Indicator of grid point outside the boundary of the country.
NOTE: to use it with setMask()
you will need to invert it using !as.logical(TellusTI$Mask)
This is a migrated version of the data set to a regular grid, ideal for illustrating and testing multiple point geostatistical algorithms.
The function is provided due to conflicting licenses. You download the data from the server https://www.bgs.ac.uk/gsni/tellus/index.html at your own accord. Actually a visit to the project data webpage is highly recommended, in particular for learning about the QA/QC process of the data acquisition, other data sources available and how to use this wealth of data.
# NOT RUN {
getwd()
getTellus(TI=TRUE, cleanup=TRUE)
dir(pattern="Tellus")
# }
Run the code above in your browser using DataCamp Workspace