db.section.plot {RGeostats}R Documentation

Display 2-D surfaces along traces

Description

Display 2-D surfaces along traces

Usage

db.section.plot(x, dbgrid, dbpoint = NA, addnames = NA, disc = 1, radius = 1, 
	cols = NA, flag.fill = TRUE, flag.revert = FALSE, flag.up = TRUE, 
	flag.z = TRUE, flag.vertex = TRUE, flag.text=FALSE, iatt.point=NA,
  fg.col = "black", col0 = "black", pos0=1,
	lwd0 = 1, lty0 = 1, title = NA, pos.legend=0,
  add = FALSE, reset = TRUE, ...)

Arguments

x

The trace defined as a matrix with 2 columns (for first two coordinates) and as many rows as vertices.

dbgrid

The db-class which contains the surfaces to be displayed. It must be defined in 2-D. The surfaces are stored in variables corresponding to "z" or "time" locators (according to 'flag.z' argument).

dbpoint

An auxiliary db-class which contains the samples to be overlaid. It must be defined in 2-D. This file must contain a depth variable (locator "z1") and a layer description (locator "layer").

addnames

A set of variables (defined by their attribute ranks) which define additional surfaces to be overlaid.

disc

Discretization step used for interpolating the different surfaces.

radius

Radius of the search distance used for selecting the samples (contained in 'dbpoint') to be represented. A sample is represented only if its distance to an trace interpolated point is smaller than this argument. The sample is attached to the closest trace interpolated point.

cols

List of colors used for representing the surfaces.

flag.fill

When TRUE, a polygon joining the interpolated surface (along the trace) to the two bottom edge corners of the section is represented. When FALSE, the surface is represented as a broken line.

flag.revert

When TRUE, the surfaces are represented starting from the one corresponding to "z1" (resp. "time1") upwards. When FALSE, the surfaces are represented ordered downwards.

flag.up

The vertical axis is represented upwards when the argument is TRUE and downwards otherwise.

flag.z

When TRUE, the variables of 'dbgrid' which correspond to the "z" locator are displayed. When FALSE, the variables corresponding to the "time" locator are represented instead.

flag.vertex

When TRUE, the vertices of thandare displayed (as vertical solid lines)

flag.text

When TRUE, display the label next to the samples located close enough to the trace.

iatt.point

This gives the attribute containing the literal value to be displayed (only used if 'flag.text' is TRUE)

fg.col

(Foreground) color used for representing the samples from 'dbpoint'.

col0

Array of colors used for representing the additional surfaces from 'dbgrid'.

pos0

Position of the label with respect to the location of the sample located close enough to the trace.

lwd0

Array of line widths for representing the additional surfaces from 'dbgrid'.

lty0

Array of line types for representing the additional surfaces from 'dbgrid'.

title

A character string representing the title.

pos.legend

Position of the legend. For more information, see legend.locate.

add

TRUE to overlay on the previous contents of graphic; FALSE otherwise

reset

When TRUE, the screen is erased before the plot is performed.

...

Additional arguments passed to define.scenes, points and abline.

Examples

# Example due to N. Desassis: Perform a Kriging based on few (simulated)
# data and represent the kriged surface along a diagonal
set.seed(23)
nech   = 20
extend = 100
x1     = seq(1,nech) * floor(extend / (nech+1))
ypos   = floor(extend/2)
dat    = db.create(x1=x1,x2=rep(ypos,nech))
grid   = db.create(nx=c(extend,extend))
neigh  = neigh.create(type=0,ndim=2)
model  = model.create(vartype="Cubic",range=extend/5)
dat    = simtub(,dat,model)

res1 = kriging(dat,grid,model,neigh,flag.std=FALSE)
plot(res1,title="Data, Estimation and Section Trace")
plot(dat,name.post=1,add=TRUE)
abline(0,1)

x = matrix(c(0,extend,ypos,ypos),2)
db.section.plot(x,res1,dat,flag.fill=FALSE,radius=5,cols="blue",cex=1.5)

[Package RGeostats version 14.0.10 Index]