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, ...)
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 |
verify |
logical. Default is |
... |
further argument |
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
x and depth or i0 must cover the same spatial region.
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 documentationPlease choose more modern alternatives, such as Google Chrome or Mozilla Firefox.