View source: R/ArcSliceUGFunctions.R
inci.mat.undAS | R Documentation |
Returns the incidence matrix
for the underlying or reflexivity graph of the AS-PCD
whose vertices are the data points in Xp
in the multiple triangle case.
AS proximity regions are defined
with respect to the Delaunay triangles based on Yp
points and
vertex regions are based on the center M="CC"
for circumcenter of each Delaunay triangle or M=(\alpha,\beta,\gamma)
in barycentric coordinates in the
interior of each Delaunay triangle;
default is M="CC"
, i.e., circumcenter of each triangle.
Loops are allowed, so the diagonal entries are all equal to 1.
Each Delaunay triangle is first converted to
an (nonscaled) basic triangle so that M
will be the same
type of center for each Delaunay triangle
(this conversion is not necessary when M
is CM
).
Convex hull of Yp
is partitioned
by the Delaunay triangles based on Yp
points
(i.e., multiple triangles are the set of these Delaunay triangles
whose union constitutes the
convex hull of Yp
points).
For the incidence matrix loops are allowed,
so the diagonal entries are all equal to 1.
See (\insertCiteceyhan:Phd-thesis,ceyhan:stamet2016;textualpcds.ugraph) for more on the AS-PCDs. Also, see (\insertCiteokabe:2000,ceyhan:comp-geo-2010,sinclair:2016;textualpcds.ugraph) for more on Delaunay triangulation and the corresponding algorithm.
inci.mat.undAS(Xp, Yp, M = "CC", ugraph = c("underlying", "reflexivity"))
Xp |
A set of 2D points which constitute the vertices of the underlying or reflexivity graph of the AS-PCD. |
Yp |
A set of 2D points which constitute the vertices of the Delaunay triangles. |
M |
The center of each triangle.
|
ugraph |
The type of the graph based on AS-PCDs,
|
Incidence matrix for the underlying or reflexivity graph
of the AS-PCD whose vertices are the 2D data set, Xp
.
AS proximity regions are constructed
with respect to the Delaunay triangles and M
-vertex regions.
Elvan Ceyhan
inci.mat.undAStri
, inci.mat.undPE
,
inci.mat.undCS
, and inci.matAS
#\donttest{
nx<-20; ny<-5;
set.seed(1)
Xp<-cbind(runif(nx,0,1),runif(nx,0,1))
Yp<-cbind(runif(ny,0,.25),
runif(ny,0,.25))+cbind(c(0,0,0.5,1,1),c(0,1,.5,0,1))
M<-c(1,1,1)
IM<-inci.mat.undAS(Xp,Yp,M)
IM
pcds::dom.num.greedy(IM)
#}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.