| getTellus {gmGeostats} | R Documentation |
Download the Tellus survey data set (NI)
Description
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.
Usage
getTellus(wd = ".", destfile = "TellusASoil.RData", TI = FALSE, cleanup = TRUE)
Arguments
wd |
working directory, where intermediate files and output will be produced |
destfile |
file name where to put the Tellus data (including the extension ".Rdata"!) |
TI |
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") |
cleanup |
should the downloaded excel file be removed? defaults to TRUE |
Details
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.
Value
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.
Geochemical soil A survey
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
Sample ID
- EASTING
X coordinate of each point
- NORTHING
Y coordinate of each point
- Flag
a flag marking some data as a study subset, fixed for reproducibility
- Ag
Silver concentration in
ppm- Cd
Cadmium concentration in
ppm- In
Indium concentration in
ppm- Sn
Tin concentration in
ppm- Sb
Antimony concentration in
ppm- Te
Tellurium concentration in
ppm- I
Indim concentration in
ppm- Cs
Caesium concentration in
ppm- Ba
Barium concentration in
ppm- La
Lanthanum concentration in
ppm- Ce
Cerium concentration in
ppm- Na2O
Sodium oxide in
percent- MgO
Magnesium oxide in
percent- Al2O3
Aluminium oxide in
percent- SiO2
Silicium oxide in
percent- P2O5
Phosphorous(V) oxide in
percent- SO3
Sulphur(III) oxide in
percent- K2O
Potasium oxide in
percent- CaO
Calcium oxide in
percent- TiO2
Titanium oxide in
percent- MnO
Manganese oxide in
percent- Fe2O3
Total Iron(III) oxide in
percent- Cl
Chlorine concentration in
ppm- Sc
Scandium concentration in
ppm- V
Vanadium concentration in
ppm- Cr
Chromium concentration in
ppm- Co
Cobalt concentration in
ppm- Ni
Nickel concentration in
ppm- Cu
Copper concentration in
ppm- Zn
Zinc concentration in
ppm- Ga
Callium concentration in
ppm- Ge
Germanium concentration in
ppm- As
Arsenic concentration in
ppm- Se
Selenium concentration in
ppm- Br
Bromine concentration in
ppm- Rb
Rubidium concentration in
ppm- Sr
Strontium concentration in
ppm- Y
Yttrium concentration in
ppm- Zr
Zirconium concentration in
ppm- Nb
Niobium concentration in
ppm- Mo
Molybdenum concentration in
ppm- Nd
Neodymium concentration in
ppm- Sm
Samarium concentration in
ppm- Yb
Ytterbium concentration in
ppm- Hf
Hafnium concentration in
ppm- Ta
Tantalum concentration in
ppm- W
Tungsten concentration in
ppm- Tl
Thallium concentration in
ppm- Pb
Lead concentration in
ppm- Bi
Bismuth concentration in
ppm- Th
Thorium concentration in
ppm- U
Uranium concentration in
ppm- pH
soil acidity
- LOI
loss on ignition in
percent
Training image
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:
- EASTING
X coordinate of each cell point
- NORTHING
Y coordinate of each cell point
- MgO
Magnesia proportion
- Al2O3
Alumina proportion
- CaO
Calcium oxide proportion
- Fe2O3
Iron oxide proportion
- Rest
Residual complementing the sum to 1
- Mask
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.
References
https://www.bgs.ac.uk/gsni/tellus/index.html
Examples
## Not run:
getwd()
getTellus(TI=TRUE, cleanup=TRUE)
dir(pattern="Tellus")
## End(Not run)