Display a USDA soil texture triangle with optional grid, labels and soil texture points.
soil.texture(soiltexture=NULL, main="", at=seq(0.1, 0.9, by=0.1),
axis.labels=c("percent sand", "percent silt",
"percent clay"),
tick.labels=list(l=seq(10, 90, by=10), r=seq(10, 90, by=10),
b=seq(10, 90, by=10)),
show.names=TRUE, show.lines=TRUE, col.names="gray",
bg.names=par("bg"), show.grid=FALSE, col.axis="black",
col.lines="gray", col.grid="gray", lty.grid=3,
show.legend=FALSE, label.points=FALSE, point.labels=NULL,
col.symbols="black", pch=par("pch"), ...)
If soiltexture was included, a list of the x,y
positions of the soil types plotted. If not, nil.
Matrix of soil textures where each row is a soil sample and three columns contain the proportions of the components sand, silt and clay in the range 0 to 1 or percentages in the range 0 to 100.
The title of the soil texture plot. Defaults to nothing.
Positions on the three axes where ticks will be drawn.
Labels for the axes.
The tick labels for the three axes.
Logical - whether to show the names of different soil types within the soil triangle.
Logical - whether to show the boundaries of the different soil types within the soil triangle.
Color of the soil names. Defaults to gray.
Color to use when drawing a blank patch for the names of soil types.
Logical - whether to show grid lines at each 10 level of each soil component.
Color of the triangular axes, ticks and labels.
Color of the boundary lines. Defaults to gray.
Color of the grid lines. Defaults to gray.
Type of line for the grid. Defaults to dashed.
Logical - whether to display a legend.
Logical - whether to call thigmophobe.labels to label the points.
Optional labels for the points or legend.
Color of the symbols representing each value.
Symbols to use in plotting values.
Additional arguments passed to triax.points and then points.
Sander Oom, Jim Lemon, and Michael Toews
soil.texture displays a triangular plot area on which soil textures defined as proportions of sand, silt and clay can be plotted. Optional grid, vertex labels, soil type divisions and names may also be displayed. If a matrix of soil textures is present, these will be plotted.
U.S. Department of Agriculture, Natural Resources Conservation Service, 2007. National Soil Survey Handbook, title 430-VI.// https://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/class/taxonomy/?cid=nrcs142p2_053577
U.S. Department of Agriculture, Natural Resources Conservation Service, 2007. Soil Texture Calculator// https://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/class/taxonomy/?cid=nrcs142p2_053577
get.soil.texture, triax.plot
data(soils)
soil.texture(main="NO DATA")
soil.texture(soils, main="DEFAULT", pch=2)
soil.texture(soils, main="LINES AND NAMES", show.lines=TRUE,
show.names=TRUE, pch=3)
soiltex.return<-soil.texture(soils[1:6,], main="GRID AND LEGEND",
show.grid=TRUE, pch=4, col.symbols=1:6, show.legend=TRUE)
par(soiltex.return$oldpar)
Run the code above in your browser using DataLab