submesh.mesh {excursions} | R Documentation |
Extract a part of a mesh
Description
Extracts a part of a mesh
Usage
submesh.mesh(z, mesh)
Arguments
z |
A matrix with values indicating which nodes that should be present in the submesh. |
mesh |
An |
Value
An fm_mesh_2d
object.
Author(s)
Finn Lindgren finn.lindgren@gmail.com
Examples
## Not run:
if (require(fmesher)) {
nxy <- 30
x <- seq(from = 0, to = 4, length.out = nxy)
lattice <- fm_lattice_2d(x = x, y = x)
mesh <- fm_mesh_2d_inla(lattice = lattice, extend = FALSE, refine = FALSE)
# extract a part of the mesh inside a circle
xy.in <- rowSums((mesh$loc[, 1:2] - 2)^2) < 1
submesh <- excursions:::submesh.mesh(matrix(xy.in, nxy, nxy), mesh)
plot(mesh$loc[, 1:2])
lines(2 + cos(seq(0, 2 * pi, length.out = 100)), 2 + sin(seq(0, 2 * pi, length.out = 100)))
plot(submesh, add = TRUE)
points(mesh$loc[xy.in, 1:2], col = "2")
}
## End(Not run)
[Package excursions version 2.5.8 Index]