Hypergraph $p$-Laplacian regularization on point clouds for data interpolation
Kehan Shi, Martin Burger
TL;DR
The paper develops a hypergraph $p$-Laplacian framework for data interpolation on point clouds lacking explicit structure, introducing $oldsymbol{ ext{$\\varepsilon_n$-ball}}$ and $oldsymbol{ ext{$k_n$-NN}}$ hypergraphs and establishing variational consistency with the continuum $p$-Laplacian. It proves Γ-convergence of the discrete energies to the continuum energies under mild scaling ($\delta_n\ll\\varepsilon_n\ll1$ and $p>d$), and extends the results to the $k_n$-NN case with weight adjustments, using TL$^p$ topology and transportation maps. The numerical algorithm leverages stochastic primal-dual hybrid gradient to solve the large-scale convex, non-differentiable problem, with convergence guarantees. Empirically, hypergraph HpL regularization suppresses spiky artifacts better than graph-based regulators, improving data interpolation in 1D tests, semi-supervised learning on MNIST, and image inpainting, albeit at higher computational cost. These findings highlight the benefit of higher-order hypergraph connections for robust SSL and image restoration when explicit structure is scarce.
Abstract
As a generalization of graphs, hypergraphs are widely used to model higher-order relations in data. This paper explores the benefit of the hypergraph structure for the interpolation of point cloud data that contain no explicit structural information. We define the $\varepsilon_n$-ball hypergraph and the $k_n$-nearest neighbor hypergraph on a point cloud and study the $p$-Laplacian regularization on the hypergraphs. We prove the variational consistency between the hypergraph $p$-Laplacian regularization and the continuum $p$-Laplacian regularization in a semisupervised setting when the number of points $n$ goes to infinity while the number of labeled points remains fixed. A key improvement compared to the graph case is that the results rely on weaker assumptions on the upper bound of $\varepsilon_n$ and $k_n$. To solve the convex but non-differentiable large-scale optimization problem, we utilize the stochastic primal-dual hybrid gradient algorithm. Numerical experiments on data interpolation verify that the hypergraph $p$-Laplacian regularization outperforms the graph $p$-Laplacian regularization in preventing the development of spikes at the labeled points.
