ich weiß nicht, ob sich hier jemand mit dem soiltexture Paket auskennt?
Ich habe bisher folgender programmiert:
Code: Alles auswählen
library(soiltexture)
####DDaten holen und Variablen definieren
PredIn <- read.table(file="PredIn_FK_Corg_4.txt", header=T, sep="\t")
#View(PredIn)
hVar <- 10^(seq(-2,6.6,0.1))
AnzahlProben <- ncol(PredIn)
#########VorbereitunG Dataframe
RetDataFramenFK0 <- as.data.frame(PredIn)
#View(RetDataFramenFK0)
#Dreieck holen und plotten
TT.plot( class.sys = "DE.BK94.TT" , tri.data = RetDataFramenFK0,
grid.show = FALSE, class.lab.show = "none" ,z.name = "nFK" ,main = "I" )
#z.name = "nFK" #col = RetDataFrameFK0$colnFK #pch = 20 , cex = c(3,1)
#Values bestimmen
z.cex.range <- TT.get("z.cex.range")
def.pch <- par("pch")
def.col <- par("col")
def.cex <- TT.get("cex")
nFK.str <- TT.str(RetDataFramenFK0[,"nFK"],
z.cex.range[1],
z.cex.range[2])
# The legend:
legend( x = 80, y = 90,
title = expression( bold("nFK [-]")),
legend = formatC(c( min( RetDataFramenFK0[,"nFK"] ),
quantile(RetDataFramenFK0[,"nFK"] ,probs=c(25,50,75)/100),
max( RetDataFramenFK0[,"nFK"]))),
pt.lwd = 4 , col = def.col,
pt.cex = c( min( nFK.str ), quantile(nFK.str ,probs=c(25,50,75)/100), max( nFK.str )),
pch = def.pch, bty = "o",bg = NA,
box.col = NA, # Uncomment this to remove the legend box
text.col = "black", cex = def.cex)
Das soll ungefähr so aussehen, aber mit dem deutschen Dreieck: Später werde ich 9 von diesen Dreiecken haben und da fände ich es eben schön diesen Farbübergang darstellen zu können.
Vielen Dank,
Eva
