Robust shape inference from a sparse approximation of the Gaussian trimmed loglikelihood

Claire Brécheteau (1, 2, 3, 4, 5)

1 ECN - École Centrale de Nantes
2 LMJL - Laboratoire de Mathématiques Jean Leray
3 LMO - Laboratoire de Mathématiques d'Orsay
4 DATASHAPE - Understanding the Shape of Data
5 SELECT - Model selection in statistical learning

Given a noisy sample of points lying around some shape M, with possibly outliers or clutter noise, we focus on the question of recovering M, or at least geometric and topological information about M. Often, such inference is based on the sublevel sets of distance-like functions such as the function distance to M, the distance-to-measure (DTM) or the k-witnessed distance. In this paper, we firstly widespread the concept of trimmed log-likelihood to probability distributions. This trimmed log-likelihood can be considered as a generalisation of the DTM. A sparse approximation of the DTM, the m-power distance-to-measure (m-PDTM), has been introduced and studied by Brécheteau and Levrard in 2017. Its sublevel sets are unions of m balls, with m possibly much smaller than the sample size. By miming the construction of the m-PDTM from the DTM, we propose an approximation of the trimmed log-likelihood associated to the family of Gaussian distributions on R^d. This approximation is sparse is the sense that its sublevel sets are unions of m ellipsoids. We provide a Lloyd-type algorithm to compute the centers and covariance matrices associated to the ellipsoids. We improve our algorithm by allowing an additional noise parameter to wipe out some points, just as the trimmed m-means algorithm of Cuesta-Albertos et al. Our algorithm comes together with a heuristic to select this parameter. Some illustrations on different examples enhance that our algorithm is efficient in wiping out clutter noise, recovering the shape and recovering the homology of M; this requiring a storage of only m points and covariance matrices.

Mots clés: ellipsoids; denoising; shape approximation; trimmed log-likelihood; Lloyd’s type algorithm

Domaines: Géométrie algorithmique [cs.CG]; Théorie [stat.TH]; Statistiques [math.ST]

Claire Brécheteau. Robust shape inference from a sparse approximation of the Gaussian trimmed loglikelihood. 2018. ⟨hal-01947424⟩