USDA Textural Soil Classification
Included in the next version (1.0.3.1) of ggtern, will be a new data set, taken from the USDA Textural Soil Classification[1] which we can use to reproduce the original diagram published by the United Stated Department of Agriculture.
Firstly let us load the packages and data:
# Load the required libraries
library(ggtern)
library(plyr)
library(grid)
# Load the Data. (Available in ggtern 1.0.3.0 next version)
data(USDA)
# Put tile labels at the midpoint of each tile.
USDA.LAB = ddply(USDA, 'Label', function(df) {
apply(df[, 1:3], 2, mean)
})
# Tweak
USDA.LAB$Angle = 0
USDA.LAB$Angle[which(USDA.LAB$Label == 'Loamy Sand')] = -35
Reconstructing the USDA’s Original Plot
Now let us construct the plot from the prepared data
# Construct the plot.
ggplot(data = USDA, aes(y=Clay, x=Sand, z=Silt)) +
coord_tern(L="x",T="y",R="z") +
geom_polygon(aes(fill = Label),
alpha = 0.75, size = 0.5, color = 'black') +
geom_text(data = USDA.LAB,
aes(label = Label, angle = Angle),
color = 'black',
size = 3.5) +
theme_rgbw() +
theme_showsecondary() +
theme_showarrows() +
custom_percent("Percent") +
theme(legend.justification = c(0, 1),
legend.position = c(0, 1)) +
labs(title = 'USDA Textural Classification Chart',
fill = 'Textural Class',
color = 'Textural Class')
Comparison to Existing Format
The above reproduction looks considerably better than the originally published format [1].
Adding Additional Points
A number of people have asked how to add additional points, try the following:
df <- data.frame(SAND=runif(10),SILT=runif(10),CLAY=runif(10)) last_plot() + geom_point(data=df,aes(y=CLAY,x=SAND,z=SILT))
References
[Bibtex]
@Book{usda,
Title = {Soil Mechanics Level 1, Module 3, USDA Textural Classification Study Guide},
Author = {Unknown-Author},
Publisher = {United States Department of Agriculture},
Year = {1987},
Owner = {Nick Hamilton},
Timestamp = {20140116}
}
-
ggtern 2.0 now available
Recently ggplot2 received a severe makeover by releasing version 2.0, and in the spirit -
User Request – Shepards Classification of Sediments
I received a request overnight on how to render the Shepard’s classification diagram, which
About Author
Nick Hamilton
Author of the ggtern package. PhD candidate at the University of New South Wales in the School of Materials Science and Engineering. My field is in the simulations of metallic glasses via several advanced simulation techniques such as QMD and RMC.
Leave a Reply
This site uses Akismet to reduce spam. Learn how your comment data is processed.
I need this diagram in my studies
How can I change the group colors?
thank you,
you need to add a scale control, where you can specify the colours of the polygon fills, try scale_fill_manual(…) or scale_fill_brewer(…)
How can add point (experimental measurements) to identify its texture
Thank¡
Hi Rodrigo,
You can add another layer of points, ie, see geom_point(…), if you need help with that, let me know.
Cheers,
Yes it will still be helpfull for me ^^
I’have a data set with 3 colums SAND, SILT, CLAY.
How can i implement it in your wonderful map ?
I have modified the post to cover your question. In summary, this is how you do it.
df <- data.frame(SAND=runif(10),SILT=runif(10),CLAY=runif(10)) last_plot() + geom_point(data=df,aes(y=CLAY,x=SAND,z=SILT))
A big thanks for your kind help on this question !
(which may be seem casual)
Just in case some people want to plot site in function of a special “variable”
Here what i’ve done:
+guides(fill=FALSE) to the ggplot() in order to erase the current legend
then:
add a site column in the df dataset:
df$site=c(rep(“a”,2),rep(“b”,2),rep(“c”,2),rep(“d”,2),rep(“e”,2));
as.factor(df$site);
Finally:
last_plot() +geom_point(data=df,aes(y=CLAY,x=SAND,z=SILT,colour=site),shape=19,show_guide = TRUE)+scale_color_discrete(as.factor(df$site),name = “site”);
Again many thanks for your website Nicholas Hamilton !
Hey, plot looks great. I can’t seem to get the code working though… do you have any suggestions as to why I’m getting the error message,
Error in (function (el, elname) :
“axis.tern.padding” is not a valid theme element name.
Sorry that theme element was removed some time ago, have updated the code, try again.
Thank you for the speedy response! The new script works. The only problem I’m running into is the fact that my axes are a bit off. I cant seem to get them to match yours or the USDA’s. Im finding clay on the right, silt on the bottom, and sand on the left, though the textural classes are correctly placed… Basically I just need to rotate the image? Sorry for so many questions… Ken
In order to get same the same results with Nicholas, you have to adjust a little bit about the code on the step # Construct the plot. Please see as follows:
# Construct the plot.
ggplot(data = USDA, aes(y=Clay, x=Sand, z=Silt)) +
coord_tern(L=”x”,T=”y”,R=”z”) +
geom_polygon(aes(fill = Label),
alpha = 0.75, size = 0.5, color = ‘black’) +
geom_text(data = USDA.LAB,
aes(label = Label, angle = Angle),
color = ‘black’,
size = 3.5) +
theme_rgbw() +
theme_clockwise() +
theme_showsecondary() +
theme_showarrows() +
custom_percent(“Percent”) +
theme(legend.justification = c(0, 1),
legend.position = c(0, 1)) +
labs(title = ‘USDA Textural Classification Chart’,
fill = ‘Textural Class’,
color = ‘Textural Class’)
Ah! it worked! The vast corners of the internet have rained down upon me the knowledge of conveniently conveying soil textural classes through the wonders of R. Thank you so much (and to Nicholas Hamilton too for designing ggtern!).
Would you mind if I use this colorized textural classification chart in a master’s thesis?
Of course not! I would appreciate a citation though.
Yes of course! Thank you very much, Nicholas.
Hi, Nicholas, will it be possible to plot points with continuous colors (for example, from white to blue) in this USDA chart? Thanks a lot!
Yes of course. Same rules for ggplot apply. This chart has a mapping to ‘fill’ already, but colour remains exposed.
If the plot in the main code is assigned to the ‘base’ variable, try the following:
Thanks a lot. I used colour option, it did work.
Hello Nicholas,
What is an effective way to label in these ternary diagrams?
Thanks for writing this!
It has been so helpful
Thiis is a great ttip particularly to those new to the blogosphere.
Brief but very accurate info� Thanks for sharing
this one. A must read post!
I delight in, result in I found just what I used to be having a look
for. You have ended my 4 day long hunt! God Bless you man. Have a great day.
Bye
Great package!!!!
Question:
I have texture data measured on different sites and at different depths.
I would like to plot these findings on this USDA chart and would like to have different colours per site and different symbols per depth.
This I achieve by doing
base + geom_point(data=df,aes(Sand,Clay,Silt,shape=Depth,colour=Site),size=4)
Is it possible to go further and make sure that each dot is coloured according to ‘Site’ while having a black border/edge? I tried doing this:
base + geom_point(data=df,aes(Sand,Clay,Silt,shape=Depth,fill=Site),colour=”black”,size=5)+
scale_shape_manual(values=c(21,24,22))
But it ends up looking horrible as the colours of the entire chart change.
I realize ‘fill’ is already used to build the chart, but is there a workaround? I tried adding two sets of points, one with size=5 that is completely black and another on top with size =4 that has the right colours, but it just ends up looking weird.
An intriguing discussion is worth comment. I do think that you ought
to publish more about this topic, it might not be a taboo subject but typically folks don’t
talk about these issues. To the next! All the best!!
Not single-handedly as the best Asian gambling operator
in the 2009 and 2010 versions of eGaming review Magazine,
but the online gambling site Sbobet has moreover become the 11th most ranked company in the “50 faculty Company” list in 2011 which is the list of
the most influential companies in the eGaming industry.
This is proof that Sbobet is not an online gambling site.
For starters, Sbobet is the best online bookie to acquit yourself Asian Handicap because here you can choose
both the type of bet and the odds you want. Compared to competitors
such as Dafabet and 188Bet, Sbobet has the
most amazing in-play betting assistance taking into consideration liven up video streaming features for several tournaments and popular matches in the midst of
loyal customers.
suitable and thank you for visiting online casino agents for
Indonesia. Here we are here to manage to pay for entry for fans of sbobet casino who desire to
link the international agent site Promosbobet.
The products offered are very diverse and are
utterly the best choice. Derived from providers that have credibility
and have been proven to have the best predicate in serving their customers.
Providers when Sbobet, Maxbet and many more.
As the best online casino and soccer gambling site in Indonesia.
We manage to pay for maximum minister to gone customer hold that works online 24
full hours all day. appropriately members can deed whenever and wherever
without fearing to incline the possibility of difficulties or confusion.
similar to reality taking into account that, of course you have to be careful in choosing
a trusted gambling site to be a conduct yourself partner.
We are here to back up you incline the difficulties and clip a long way to find the best and most trusted gambling agent and the safest ball in Indonesia.
The Mabosbet site is the most secure ball trusted trusted agent that has been providing the best minister
to to its customers for years. PAGCOR has an qualified license from an authorized body that oversees gambling regulations
in the Asian region. This license guarantees the safest soccer gambling agent that
all games will be provided by Mabosbet guaranteed safe and have A-grade quality.