plotFlow {SoPhy}R Documentation

plotting of the simulated infiltration pattern

Description

plotFlow2d plots two dimensional profiles of simulated stochastic flow patterns; plotFlow3d plots a 3d, perspective figure of simulated stochastic flow patterns.

Usage

plotFlow2d(coord, pointradius=1, slice=2 * pointradius,
           full.size=TRUE, Profiles=1, dev=1, ps="",
           height=4, unit="cm", cex=2, correction=1.2, col=1,
           rl = function(x) readline(paste(x, ": press return")))

plotFlow3d(paths, horizons=c("no", "absorbing", "breakthrough"),
           drop.distr, n.balls=1, pointradius=1,
           dev=1, ps="3d.dye.pattern", ps.background=FALSE,
           profileheight=4, unit="cm", unit.scale=1, inf, sun,
           rl = function(x) readline(paste(x, ": press return")),
           low.resolution=TRUE,
           col=grey(pmin(1, pmax(0, seq(0.95, 0, -0.001)))))

Arguments

coord output of plotFlow3d: n x 3 matrix; 3d coordinates for the centre of the balls that approximate the tube around each path
pointradius radius of the dyed tube around each path; the tube is approximated by a sequence of balls
slice thickness of the slice for which the paths are visible
full.size logical. If FALSE then the projected radii of the sliced balls are plotted. If TRUE then instead of the projected radius the pointradius is used.
Profiles number of the in the y-direction equally spaced profiles of the 3d flow pattern
dev plotting device, see Dev
ps name of the postscript file, see Dev
height height of the plotted figure, see Dev
unit string; if empty then labels are not given. Otherwise unit denotes the units of the plot
cex used for cex.axis and cex.lab, see par
correction positive real number. Do not change its value. This constant takes into account that the plotting procedure used in the function does not plot a circle of exactly the given radius.
col colour of the points; plotFlow3d needs a vector of colours starting with the colour of the nearest points
paths list as returned by flowpattern
horizons Temptative parameter. If "absorbing" the path will be ‘infinitely’ long in a layer at depth 50 about, so that fluxes will stop in this layer If "breakthrough" the path will be very long in a layer at depth 50 about, so that only a few paths are percolated.
The values "absorbing" and "breakthrough" are used to create the figures in the Discussion of Schlather and Huwe (2004a).
drop.distr function or missing. The stained path length reached by a single drop is calculated as follows. First a Gaussian random field is simulated with the parameters given below. Then the Gaussian random field is marginally transformed to uniform distribution by the Gaussian distribution function. Then the inverse of the required distribution function is applied. drop.distr gives this inverse. It should allow for vectors as input variables. The drop.distr used to create paths is overwritten if drop.distr is not missing, but the realisation of the paths is always kept.
n.balls the tube for each path segment is approximated by n.balls balls
ps.background logical. If dev=TRUE then a postscript file is created, which is, in general, huge. Therefore, any postscript file created on a unix system is transformed into a tiff-file by `convert' and back again to a postscript file. These transformation are time consuming. If ps.background=TRUE then these transformations are started in the background.
profileheight essentially the height of the figure; the value is used to calculated the width of the plot. The final height of the plot is higher then profileheight because additional space needed for the perspective plot
unit.scale scaling factor for the numerical values of all three-dimensional coordinates
inf Vector of three components. The first two give the position of the infinitely far point; the third component gives the extension factor for the second dimension that is orthogonal to the screen; inf is of order 1/100
sun the three dimensional coordinates of the fictive sun
rl function. Called, after a figure, but the last, is plot on the screen.
low.resolution logical. If TRUE and dev=FALSE, i.e. postscript file will be produced, then on unix systems the postscript file will be rewritten to have much smaller size, but also lower resolution.

Value

plotFlow2d returns NULL. plotFlow3d returns invisibly a matrix for the centers of the approximating balls

Author(s)

Martin Schlather, martin.schlather@math.uni-goettingen.de http://www.stochastik.math.uni-goettingen.de/institute

References

See Also

flowpattern

Examples

 ## see flowpattern

[Package SoPhy version 1.0.34 Index]