AlphaComplex | R Documentation |
AlphaComplex is a simplicial complex constructed from the finite cells of a Delaunay Triangulation.
The filtration value of each simplex is computed as the square of the circumradius of the simplex if the circumsphere is empty (the simplex is then said to be Gabriel), and as the minimum of the filtration values of the codimension 1 cofaces that make it not Gabriel otherwise. All simplices that have a filtration value strictly greater than a given alpha squared value are not inserted into the complex.
rgudhi::PythonClass
-> AlphaComplex
new()
AlphaComplex
constructor.
AlphaComplex$new(points, precision = "safe")
points
Either a n x d
matrix or a length-n
list of
d
-dimensional vectors or a file with extension .off
.
precision
A string specifying the alpha complex precision. Can be
one of "fast"
, "safe"
or "exact"
. Defaults to "safe"
.
A AlphaComplex
object storing the Alpha complex.
create_simplex_tree()
Generates a simplex tree from the Delaunay triangulation.
AlphaComplex$create_simplex_tree( max_alpha_square = Inf, default_filtration_value = FALSE )
max_alpha_square
A numeric value specifying the maximum alpha
square threshold the simplices shall not exceed. Default is set to
Inf
, and there is very little point using anything else since it does
not save time.
default_filtration_value
A boolean specifying whether filtration
values should not be computed and will be set to NaN
(default_filtration_value = TRUE
). Defaults to FALSE
(which means
compute the filtration values).
A SimplexTree
object storing the computed simplex
tree.
get_point()
This function returns the point corresponding to a given vertex from the SimplexTree.
AlphaComplex$get_point(vertex)
vertex
An integer value specifying the desired vertex.
A numeric vector storing the point corresponding to the input vertex.
clone()
The objects of this class are cloneable with this method.
AlphaComplex$clone(deep = FALSE)
deep
Whether to make a deep clone.
Vincent Rouvreau
Other filtrations and reconstructions:
RipsComplex
,
TangentialComplex
,
WitnessComplex
Xl <- seq_circle(10)
Xm <- Reduce(rbind, Xl, init = numeric())
acm <- AlphaComplex$new(points = Xm)
acl <- AlphaComplex$new(points = Xl)
acl
X <- seq_circle(10)
ac <- AlphaComplex$new(points = X)
st <- ac$create_simplex_tree()
X <- seq_circle(10)
ac <- AlphaComplex$new(points = X)
st <- ac$create_simplex_tree()
ac$get_point(1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.