plotRF {SoPhy} | R Documentation |
the function plots random fields including stones and roots using
the definition of the horizons in xswms2d
.
plotRF(h, col.txt='black', col.stones='white', col.rf=if (is.null(h$col.rf)) rainbow(100) else h$col.rf, titl=TRUE, line = 0.7, lim=1, quadratic=TRUE, cex=1, cex.leg=0.8, pch=5, cex.pch = max(0.05, 3/min(length(h$grid.x), length(h$grid.y))), root.col=grey(1 - (i-1) / (1.6 * length(h$plants))), what=c('Root.RF', 'Stone.RF', 'RF'), transf=c("K", "H", "theta", "none"), legend = TRUE, ylim, zlim)
h |
a list of the
same format as the output of xswms2d |
col.txt |
colour of the title |
col.stones |
colour of the stones and the air |
col.rf |
colour spectrum of the random field; plants are always grey |
titl |
logical. If TRUE a title is plotted using the
colour col.txt |
line |
place where the titel is plottes, see
title . |
lim |
value in [0,1] that is used only if
zlim is missing; quantile for the upper bound
zlim[2] ; the value of lim should be
less than 1 if the random field contains some extreme values, since
the colour scale is linear |
quadratic |
logical. If TRUE the figure matrix is enlarged and
filled with NA symmetrically in x-direction or on the bottom in
y-direction such the matrix (and the figure) become quadratic |
cex |
A numerical value giving the amount by which the title text, the axes and the labels of the axes should be scaled relative to the default |
cex.leg |
A numerical value giving the amount by which the legend text should be scaled relative to the default |
pch |
the symbol for a root segment |
cex.pch |
size of the symbol |
root.col |
the colours for a root segments |
what |
chooses between the random field that includes modifications by roots and stones or by stones allone, or the random field without modifications |
transf |
Kind of random field:
saturated hydraulic conductivity "K" ,
pressure head "H" ,
the water content "theta" ,
or the original Gaussian random field ("none" ) |
legend |
logical. If (TRUE) a legend is added. |
ylim |
missing or vector of two components. If missing the range of the simulation is used. |
zlim |
missing or vector of two components. Values below
zlim[1] are plotted in white,
values above zlim[2] are plotted also in white.
If missing then zlim[1] is the minimal value
and zlim[2] is calculated by means
of lim .
|
If transf="K"
then
the function plots the random field of hydraulic conductivity, i.e.,
the Gaussian random field after
h$miller.link
and h$millerK
have been applied.
the matrix of the random field values, always excluding the plants, if no error has occured and an error message otherwise.
Martin Schlather, martin.schlather@math.uni-goettingen.de http://www.stochastik.math.uni-goettingen.de/institute
## see simulateHorizons