meshing {RGeostats} | R Documentation |
Perform Delaunay triangulation
meshing(db=NA, dbaux=NA, triswitch = "Q", gext=NA, verbose = FALSE, flag.plot=FALSE, ...)
db |
First |
dbaux |
Auxilary |
triswitch |
Command line which may contain the following keywords:
|
gext |
The 'dbaux' may be dilated by gext. This argument designates an array, with its dimension equal to the dimension of the space and which contains the extension in each direction. |
verbose |
Verbose flag |
flag.plot |
When TRUE, the constraining vertices, holes and segments are represented. The created triangles are also represented. |
... |
Arguments passed for graphic function such as |
The meshing list with the following contents:
points Matrix of Network vertices (Dimension: nrow=npoint * ncol=ndim)
meshes Matrix of Vertex indices for each mesh of the Network (Dimension: nrow=nmesh * ncol=ncorner)
holes Matrix of Hole coordinates (Dimension: nrow=npoint * ncol=ndim)
segments Matrix of Vertex indices for the segments (Dimension: nrow=nsegments * ncol=2)
inddat which gives the indices of the samples of 'db' in Q
nbnodup which gives the number of samples of Dbin which do not coincide with samples of Dbout.
dupl.data Indices of the data samples close to a grid node
dupl.grid Indices of the target points close to a data sample
inddat which gives the indices of the samples of 'db' in Q
nbin which gives the number of samples in Dbin
nbout which gives the number of samples in Dbout
nbstein which gives the number of Steiner points added by meshing
nbdup which gives the number of samples of Dbin which coincide with samples of Dbout
nbnodup which gives the number of samples of Dbin which do not coincide with samples of Dbout.
Note that the following parameters can be indirectly accessed (they have been suppressed from the output list):
ndim Space dimension. Obtained as dim(points)[2]
npoint Number of Vertices in Network. Obtained as dim(points)[1]
ncorner Number of corners in each meshing element. Obtained as dim(meshes)[2]
nmesh Number of meshes in the Network. Obtained as dim(meshes)[1]
# Load the Data Set from Jonathan Richard Shewchuk web site. rg.load("Demo.double.hex.db","tri.db") # Perform the basic triangulation # res = meshing(db=tri.db,triswitch="Q") # Perform the triangulation with surface constraint # res = meshing(db=tri.db,triswitch="Qa0.005") # Remove file rm(tri.db, tri.res)