plotRF {SoPhy} | R Documentation

## Plotting random fields

### Description

the function plots random fields including stones and roots using the definition of the horizons in `xswms2d`.

### Usage

      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)

### Arguments

`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`.

### Details

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.

### Value

the matrix of the random field values, always excluding the plants, if no error has occured and an error message otherwise.

### Author(s)

Martin Schlather, schlath@hsu-hh.de http://www.unibw-hamburg.de/WWEB/math/schlath/schlather.html

### See Also

`plotWater`, `xswms2d`

### Examples

       ## see simulateHorizons

[Package SoPhy version 1.0.28 Index]