plotPEregs | R Documentation |
Plots the points in and outside of the Delaunay triangles
based on Yp
points which partition
the convex hull of Yp
points and also plots the PE proximity regions
for Xp
points and the Delaunay triangles based on Yp
points.
PE proximity regions are constructed
with respect to the Delaunay triangles with the expansion parameter
r \ge 1
.
Vertex regions in each triangle is
based on the center M=(\alpha,\beta,\gamma)
in barycentric coordinates in the interior of each Delaunay triangle
or based on circumcenter of
each Delaunay triangle (default for M=(1,1,1)
which is the center of mass of the triangle).
See (\insertCiteceyhan:Phd-thesis,ceyhan:arc-density-PE,ceyhan:dom-num-NPE-Spat2011;textualpcds) for more on the PE proximity regions. Also, see (\insertCiteokabe:2000,ceyhan:comp-geo-2010,sinclair:2016;textualpcds) for more on Delaunay triangulation and the corresponding algorithm.
plotPEregs(
Xp,
Yp,
r,
M = c(1, 1, 1),
asp = NA,
main = NULL,
xlab = NULL,
ylab = NULL,
xlim = NULL,
ylim = NULL,
...
)
Xp |
A set of 2D points for which PE proximity regions are constructed. |
Yp |
A set of 2D points which constitute the vertices of the Delaunay triangles. |
r |
A positive real number
which serves as the expansion parameter in PE proximity region;
must be |
M |
A 3D point in barycentric coordinates
which serves as a center in the interior of each Delaunay
triangle or circumcenter of each Delaunay triangle
(for this, argument should be set as |
asp |
A |
main |
An overall title for the plot (default= |
xlab , ylab |
Titles for the |
xlim , ylim |
Two |
... |
Additional |
Plot of the Xp
points,
Delaunay triangles based on Yp
points
and also the PE proximity regions
for Xp
points inside the convex hull of Yp
points
Elvan Ceyhan
plotPEregs.tri
, plotASregs
,
and plotCSregs
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-20; ny<-5; #try also nx<-40; ny<-10 or nx<-1000; ny<-10;
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))
#try also Yp<-cbind(runif(ny,0,1),runif(ny,0,1))
M<-c(1,1,1) #try also M<-c(1,2,3)
r<-1.5 #try also r<-2
plotPEregs(Xp,Yp,r,M,xlab="",ylab="")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.