meshing {RGeostats} | R Documentation |
Perform Meshing
meshing(db=NA, dbaux=NA, triswitch = "Q", gext=NA, verbose = FALSE, flag.plot=FALSE, ...)
db |
First |
dbaux |
Auxilary |
triswitch |
The argument 'triswitch' depends on the type of meshing (i.e triangles, tetrehedra or triangles on the sphere). For Triangulation, please use the following keywords:
For Tetrahedralization, please use the following keywords:
For Triangulation on the Sphere, please use the following keywords (followed by 'xxx') to generate points on the Sphere
|
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 all active samples of input 'Db' in meshing
nbnodup which gives the number of samples of Dbin which do not coincide with samples of Dbout.
dupl.data Indices of the duplicate data samples in the meshing
dupl.grid Indices of the duplicated target points in the meshing
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)