mlayer.vario {RGeostats}R Documentation

Calculate the Variogram in the Multi-Layer framework

Description

Calculate the Variogram in the Multi-Layer framework

Usage

mlayer.vario(dbin, dbout, nlayers = 2, dirvect = NA, lag = NA,
             toldis = 0.5, tolang = NA, cylrad = 0, nlag = NA, breaks = NA,
             opt.code = 0, tolcode = 0, flag.vel = FALSE, flag.ext=FALSE,
             irf.rank=0, match.time = FALSE,
             colrefd = NA, colreft = NA, verbose=FALSE)

Arguments

dbin

The db-class structure containing the 2-D data file. The file must be 2-D and contain one Z-variable corresponding to the sample depth and a locator "Layer".

dbout

The db-class structure containing the 2-D ouput file where the additional variables must be defined.

nlayers

Number of layers

dirvect

The directions in which the variogram must be calculated. For more information, refer to get.directions

lag

Array containing the distance lags for each calculation direction. If no corresponding lag is not defined, a default lag is calculated so that the maximum distance is equal to half of the field diagonal. For one distance, if the lag is not defined, it is set to the default lag.

toldis

Tolerance on the distance, expressed as a proportion of the lag. This parameter is only used in the case of regular lags.

tolang

Array containing the angular tolerance for each calculation direction. If no angular tolerance is not defined, a default tolerance is calculated as one half of the angle. For one direction, if the angular tolerance is not defined, it is set to the default tolerance.

cylrad

Array containing the cylrad tolerance for each calculation direction. When specified and positive, a pair is discarded if the separation distance in the plane orthogonal to the calculation distance is larger than 'cylrad'.

nlag

Array containing the number of lags for each calculation direction. If not defined, the defaulted number of lags is set to 10. For one direction, if the number of lags is not defined, it is set to 10.

breaks

Vector giving the boundaries of irregular distance lags. For multidrectional computations, breaks is either: a vector (The breaks are then the same for all directions) a list of vectors of the breaks for each direction. If provided, nlag is adapted to the breaks. Intervals are open left and closed right (e.g ]0;15]).

opt.code

Array containing the option concerning the use of the code for each direction. This option defines the strategy for the code usage when considering a pair of samples:

  • 0 : a pair of samples is selected whatever their codes (if active)

  • 1 : a pair is retained if the code is active in the data base and the codes of the two samples are closer than tolcode

  • 2 : a pair is retained only if the codes of the two samples are different

tolcode

Array giving the maximum distance between the codes of the two samples for each direction.

flag.vel

Should be TRUE if the estimation must be performed using velocities (the Model should correspond to interval velocities) or FALSE when working with depth (the Model should correspond to interval thickness). If flag.vel is TRUE, the layer time maps should be defined in the grid output file (locator "time" or "f" if match_time is TRUE).

flag.ext

Should be TRUE if External Drift should be used. Then the grid output file should contain the corresponding maps (locator "f").

irf.rank

Rank of the Intrisic Random Function. Should be -1 (for strict stationarity), 0 (for intrinsic) or 1 (for a first order drift).

match.time

When TRUE, the external drift maps and the time maps coincide in the output grid file. They are both defined using the locator "f".

colrefd

Rank of the attribute containing the Reference Depth map. If 0, the reference depth is set to 0

colreft

Rank of the attribute containing the Reference Time map. If 0, the reference time is set to 0.

verbose

Verbose flag

Value

A vario-class structure containing the experimental variogram of the hidden variables.


[Package RGeostats version 11.0.6 Index]