geotopbricks (version 1.3.3)

getvalues.brick.at.depth: Interpolates the values of a 'brick' at a certain depth and returns the map of brick values at the "depth" level

Description

Interpolates the values of a 'brick' at a certain depth and returns the map of brick values at the "depth" level

Usage

getvalues.brick.at.depth(x, depth, layers, i0 = NULL,
    verify = FALSE, ...)

Arguments

x

a 'RasterBrick' or a three-dimensional array

depth

depth map, generally a 'RasterLayer' object

layers

vector of layer thickness

i0

a 'Raster' containing the number of soil laver just over the bedrock. Default is NULL and is then calculated.

verify

logical. Default is FALSE. If it is TRUE, it verifies that function is working correctly.

...

further argument

Value

• a list of 'Raster' maps:

  i0 a 'Raster' containing the number of soil laver just over the bedrock

  val_z0 a 'Raster' containing the values of x at the i0-th layer

  val_z1 a 'Raster' containing the values of x at the (i0+1)-th layer

  z0 a 'Raster' containing the depth of the center of the i0-th layer

  z1 a 'Raster' containing the depth of the center of the (i0+1)-th layer

See Also

code{vertical.aggregate.brick.within.depth}

Examples

library(geotopbricks)
# The examples is the following R script conteined in a 'inst' directory of the package source
f <- system.file("doc/examples/example.getvalues.brick.at.depth.R",package="geotopbricks")
#  source(f) # Uncomment this line to run the example.
# You can copy the example file using file.copy(from=f,to=....,...) See file.copy documentation