EPFL | Biomedical Imaging Group

TV-Based Spline Reconstruction with Fourier Measurements: Uniqueness and Convergence of Grid-Based Methods

T. Debarre, Q. Denoyelle, J. Fageot

Journal of Computational and Applied Mathematics, vol. 422, paper no. 114937, 13 p., April 2023.

We study the problem of recovering piecewise-polynomial periodic functions from their low-frequency information. This means that we only have access to possibly corrupted versions of the Fourier samples of the ground truth up to a maximum cutoff frequency K_c . The reconstruction task is specified as an optimization problem with total-variation (TV) regularization (in the sense of measures) involving the Mth order derivative regularization operator L = D^M. The order M ≥ 1 determines the degree of the reconstructed piecewise-polynomial spline, whereas the TV regularization norm, which is known to promote sparsity, guarantees a small number of pieces. We show that the solution of our optimization problem is always unique, which, to the best of our knowledge, is a first for TV-based problems. Moreover, we show that this solution is a periodic spline matched to the regularization operator L whose number of knots is upper-bounded by 2K_c . We then consider the grid-based discretization of our optimization problem in the space of uniform L-splines. On the theoretical side, we show that any sequence of solutions of the discretized problem converges uniformly to the unique solution of the gridless problem as the grid size vanishes. Finally, on the algorithmic side, we propose a B-spline-based algorithm to solve the discretized problem, and we demonstrate its numerical feasibility experimentally. On both of these aspects, we leverage the uniqueness of the solution of the original problem.

@ARTICLE(http://bigwww.epfl.ch/publications/debarre2301.html,
AUTHOR="Debarre, T. and Denoyelle, Q. and Fageot, J.",
TITLE="{TV}-Based Spline Reconstruction with {F}ourier Measurements:
	{U}niqueness and Convergence of Grid-Based Methods",
JOURNAL="Journal of Computational and Applied Mathematics",
YEAR="2023",
volume="422",
number="",
pages="",
month="April",
note="paper no.\ 114937")

© 2023 The Authors. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from The Authors. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.