High-quality triangulations

Create a triangle mesh based on initial point locations, specified or automatic boundaries, and mesh quality parameters.

inla.mesh.2d(loc = NULL,
             loc.domain = NULL,
             offset = NULL,
             n = NULL,
             boundary = NULL,
             interior = NULL,
             max.edge,
             min.angle = NULL,
             cutoff = 1e-12,
             max.n.strict = NULL,
             max.n = NULL,
             plot.delay = NULL,
             crs = NULL)

Arguments

loc

Matrix of point locations to be used as initial triangulation nodes. Can alternatively be a SpatialPoints or SpatialPointsDataFrame object.

loc.domain

Matrix of point locations used to determine the domain extent. Can alternatively be a SpatialPoints or SpatialPointsDataFrame object.

offset

The automatic extension distance. One or two values, for an inner and an optional outer extension. If negative, interpreted as a factor relative to the approximate data diameter (default=-0.10???)

n

The number of initial nodes in the automatic extensions (default=16)

boundary

A list of one or two inla.mesh.segment objects describing domain boundaries.

interior

An inla.mesh.segment object describing desired interior edges.

max.edge

The largest allowed triangle edge length. One or two values.

min.angle

The smallest allowed triangle angle. One or two values. (Default=21)

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.

max.n.strict

The maximum number of vertices allowed, overriding min.angle and max.edge (default=-1, meaning no limit). One or two values, where the second value gives the number of additional vertices allowed for the extension.

max.n

The maximum number of vertices allowed, overriding max.edge only (default=-1, meaning no limit). One or two values, where the second value gives the number of additional vertices allowed for the extension.

plot.delay

On Linux (and Mac if appropriate X11 libraries are installed), specifying a nonnegative numeric value activates a rudimentary plotting system in the underlying fmesher program, showing the triangulation algorithm at work, with waiting time factor plot.delay between each step.

On all systems, specifying any negative value activates displaying the result after each step of the multi-step domain extension algorithm.

crs

An optional CRS or inla.CRS object

Value

An inla.mesh object.

Author

Finn Lindgren finn.lindgren@gmail.com

See also

inla.mesh.create, inla.delaunay, inla.nonconvex.hull

Examples

loc <- matrix(runif(10*2),10,2)

if (require("splancs")) {
  boundary <- list(inla.nonconvex.hull(loc, 0.1, 0.15),
                   inla.nonconvex.hull(loc, 0.2, 0.2))
  offset <- NULL
} else {
  boundary <- NULL
  offset <- c(0.1, 0.2)
}
#> Loading required package: splancs
#> 
#> Spatial Point Pattern Analysis Code in S-Plus
#>  
#>  Version 2 - Spatial and Space-Time analysis
#> Error in inla.has_PROJ6(): inla.require("rgdal") is not TRUE
mesh <- inla.mesh.2d(loc, boundary=boundary, offset=offset, max.edge=c(0.05, 0.1))
#> Error: object 'boundary' not found

plot(mesh)
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'mesh' not found