Low-rank tensor recovery using sequentially optimal modal projections in iterative hard thresholding (SEMPIHT)

© 2017, Society for Industrial and Applied MathematicsInternational audienceIterative hard thresholding (IHT) is a simple and effective approach to parsimonious data recovery. Its multilinear rank (mrank)-based application to low-rank tensor recovery (LRTR) is especially valuable given the difficulties involved in this problem. In this paper, we propose a novel IHT algorithm for LRTR, choosing sequential per-mode SVD truncation as its thresholding operator. This operator is less costly than those used in existing IHT algorithms for LRTR, and often leads to superior performance. Furthermore, by exploiting the sequential optimality of the employed modal projections, we derive recovery guarantees relying on restricted isometry constants. Though these guarantees are suboptimal, our numerical studies indicate that a quasi-optimal number of Gaussian measurements suffices for perfect data reconstruction. We also investigate a continuation technique which yields a sequence of progressively more complex estimated models until attaining a target mrank. When recovering real-world data, this strategy stabilizes the estimation error and can also accelerate convergence. In tensor completion, in particular, it can cope with nonideal characteristics of the sensed tensors and so is crucial for achieving a satisfactory performance. Extensive numerical experiments are reported, including the completion of hyperspectral imaging data and comparisons with several other existing approaches