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

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

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

Note

x and depth or i0 must cover the same spatial region.

See also

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