R: Estimation or Simulation of Multilayer surfaces

m2d.spde {RGeostats}

R Documentation

Estimation or Simulation of Multilayer surfaces

Description

Estimation or Simulation of Multilayer surfaces

Usage

m2d.spde(dbin, dbout, model, nlayer = 1, niter = 100, seed = 232131, 
				 nbsimu = 1, flag.ed = TRUE, flag.drift = FALSE, flag.ce = FALSE, 
				 flag.cstd = FALSE, iptr.pinch = -1, verbose = FALSE, 
				 radix = "m2D", modify.target = db.locmod())

Arguments

`dbin`	The `db-class` 2-D structure containing the dat information. For each layer, this information may involve: the constraining lower bound (pointer "lower") the constraining upper bound (pointer "upper") the elevation hard data (pointer "z") All these variables are optional. When all the variables are defined, please refer to the 'Details' section for further explanations.
`dbout`	The `db-class` 2-D structure which contains the resulting (estimated or simulated) surfaces.
`model`	The `model-class` structure containing the Model information.
`nlayer`	The number of layers
`niter`	The number of iterations for the Gibbs sampler
`seed`	The seed used for the random number generator.
`nbsimu`	The number of simulations. When set to 0, the procedure only fits the optimal drift (same as when 'flag.drift=TRUE'
`flag.ed`	When TRUE, each layer is estimated (or simulated) taking an external surface as external drift. This surface must have the locator "f" for each layer.
`flag.drift`	When TRUE, the optimal drift is calculated for each layer.
`flag.ce`	When TRUE, the simulations (after they have been calculated in a standard manner) are averaged on a pixel basis in order to produce the conditional expectation.
`flag.cstd`	When TRUE, the simulations (after they have been calculated in a standard manner) are compared on a pixel basis in order to produce the standard deviation of the conditional expectation.
`iptr.pinch`	This corresponds to the rank of an attribute of the 'dbout' file which contains a scaled distance value. This value is used as a multiplicative coefficient applied to the trend of the thickness between layers. Consequently, having this value equal to 0 corresponds to an area where this thickness is null: a pinchout.
`verbose`	Verbose flag
`radix`	Radix of the name given to the variables storing the results in the target Db.
`modify.target`	Decides whether or not the newly created variables will have their locator defined or not. For more information, see `db.locmod`.

Details

When the input db-class contains both lower-upper inequality and hard data, the hard data prevails. A warning is issued if the hard data does not lie within the lower-upper inequality interval.

Value

The db-class output structure which contains:

The optimal drift for each layer (if flag.drift=TRUE)
The conditional expectation for each layer (if flag.ce=TRUE)
The standard deviation of the conditional expectation for each layer (if flag.cstd=TRUE)
otherwise the simulations for each layer

[Package RGeostats version 14.0.10 Index]