Snippet to calculate persistence diagram from point cloud data

persistence_point_cloud.py

import gtda.homology as hl

# represent data as a point cloud
point_cloud = ...

# define topological features to track
homology_dimensions = [0, 1, 2]

# define simplicial complex to construct
persistence = hl.VietorisRipsPersistence(
    metric="euclidean", homology_dimensions=homology_dimensions
)

# calculate persistence diagram
persistence_diagram = persistence.fit_transform(point_cloud)