Progressive Wasserstein Barycenters of Persistence Diagrams

Abstract

This paper presents an efficient algorithm for the progressive approximation of Wasserstein barycenters of persistence diagrams, with applications to the visual analysis of ensemble data. Given a set of scalar fields, our approach enables the computation of a persistence diagram which is representative of the set, and which visually conveys the number, data ranges and saliences of the main features of interest found in the set. Such representative diagrams are obtained by computing explicitly the discrete Wasserstein barycenter of the set of persistence diagrams, a notoriously computationally intensive task. In particular, we revisit efficient algorithms for Wasserstein distance approximation [12,51] to extend previous work on barycenter estimation [94]. We present a new fast algorithm, which progressively approximates the barycenter by iteratively increasing the computation accuracy as well as the number of persistent features in the output diagram. Such a progressivity drastically improves convergence in practice and allows to design an interruptible algorithm, capable of respecting computation time constraints. This enables the approximation of Wasserstein barycenters within interactive times. We present an application to ensemble clustering where we revisit the k-means algorithm to exploit our barycenters and compute, within execution time constraints, meaningful clusters of ensemble data along with their barycenter diagram. Extensive experiments on synthetic and real-life data sets report that our algorithm converges to barycenters that are qualitatively meaningful with regard to the applications, and quantitatively comparable to previous techniques, while offering an order of magnitude speedup when run until convergence (without time constraint). Our algorithm can be trivially parallelized to provide additional speedups in practice on standard workstations. [...]

References in corpus (9)

• A Comparative Study of Efficient Initialization Methods for the K-Means Clustering Algorithm
• The persistent cosmic web and its filamentary structure I: Theory and implementation
• Fast Computation of Wasserstein Barycenters
• The Topology ToolKit
• Sliced Wasserstein Kernel for Persistence Diagrams
• Felix: A Topology based Framework for Visual Exploration of Cosmic Filaments
• Topological Machine Learning with Persistence Indicator Functions
• Large Scale computation of Means and Clusters for Persistence Diagrams using Optimal Transport
• Strong Equivalence of the Interleaving and Functional Distortion Metrics for Reeb Graphs

Cited by in corpus (8)

• Scalar Field Comparison with Topological Descriptors: Properties and Applications for Scientific Visualization
• Fuzzy c-Means Clustering for Persistence Diagrams
• Merge Tree Geodesics and Barycenters with Path Mappings
• Fast Topological Clustering with Wasserstein Distance
• Discrete Morse Sandwich: Fast Computation of Persistence Diagrams for Scalar Data -- An Algorithm and A Benchmark
• Fast Unbalanced Optimal Transport on a Tree
• Estimation and Quantization of Expected Persistence Diagrams
• A Domain-Oblivious Approach for Learning Concise Representations of Filtered Topological Spaces for Clustering