inla.mesh.create {INLA} | R Documentation |
Create a constrained refined Delaunay triangulation (CRDT) for a set of spatial locations.
inla.mesh.create( loc = NULL, tv = NULL, boundary = NULL, interior = NULL, extend = (missing(tv) || is.null(tv)), refine = FALSE, lattice = NULL, globe = NULL, cutoff = 1e-12, plot.delay = NULL, data.dir, keep = (!missing(data.dir) && !is.null(data.dir)), timings = FALSE, quality.spec = NULL, crs = NULL ) inla.delaunay(loc, ...)
loc |
Matrix of point locations. Can alternatively be a
|
tv |
A triangle-vertex index matrix, specifying an existing triangulation. |
boundary |
A list of |
interior |
A list of |
extend |
Setting to |
refine |
|
lattice |
An |
globe |
Subdivision resolution for a semi-regular spherical triangulation with equidistant points along equidistant latitude bands. |
cutoff |
The minimum allowed distance between points. Point at most as far apart as this are replaced by a single vertex prior to the mesh refinement step. |
plot.delay |
On Linux (and Mac if appropriate X11 libraries are
installed), specifying a numeric value activates a rudimentary plotting
system in the underlying |
data.dir |
Where to store the |
keep |
|
timings |
If |
quality.spec |
List of vectors of per vertex |
crs |
An optional |
... |
Optional parameters passed on to |
inla.mesh.create
generates triangular meshes on subsets of
R^2 and S^2. Use the higher level wrapper function
inla.mesh.2d()
for greater control over mesh resolution and
coarser domain extensions.
inla.delaunay
is a wrapper function for obtaining the convex hull of
a point set and calling inla.mesh.create
to generate the classical
Delaunay tringulation.
An inla.mesh
object.
Finn Lindgren finn.lindgren@gmail.com
inla.mesh.2d()
, inla.mesh.1d()
,
inla.mesh.segment()
, inla.mesh.lattice()
,
inla.mesh.query()
loc <- matrix(runif(10 * 2), 10, 2) mesh <- inla.delaunay(loc) plot(mesh) mesh <- inla.mesh.create(loc, interior = inla.mesh.segment(idx = 1:2), extend = TRUE, refine = list(max.edge = 0.1) ) plot(mesh) loc2 <- matrix(c(0, 1, 1, 0, 0, 0, 1, 1), 4, 2) mesh2 <- inla.mesh.create( loc = loc, boundary = inla.mesh.segment(loc2), interior = inla.mesh.segment(idx = 1:2), quality.spec = list(segm = 0.2, loc = 0.05), refine = list(min.angle = 26) ) plot(mesh2)