Abstract:Some recent methods for lossy signal and image compression store only a few selected pixels and fill in the missing structures by inpainting with a partial differential equation (PDE). Suitable operators include the Laplacian, the biharmonic operator, and edge-enhancing anisotropic diffusion (EED). The quality of such approaches depends substantially on the selection of the data that is kept. Optimising this data in the domain and codomain gives rise to challenging mathematical problems that shall be addressed in our work. In the 1D case, we prove results that provide insights into the difficulty of this problem, and we give evidence that a splitting into spatial and tonal (i.e. function value) optimisation does hardly deteriorate the results. In the 2D setting, we present generic algorithms that achieve a high reconstruction quality even if the specified data is very sparse. To optimise the spatial data, we use a probabilistic sparsification, followed by a nonlocal pixel exchange that avoids getting trapped in bad local optima. After this spatial optimisation we perform a tonal optimisation that modifies the function values in order to reduce the global reconstruction error. For homogeneous diffusion inpainting, this comes down to a least squares problem for which we prove that it has a unique solution. We demonstrate that it can be found efficiently with a gradient descent approach that is accelerated with fast explicit diffusion (FED) cycles. Our framework allows to specify the desired density of the inpainting mask a priori. Moreover, is more generic than other data optimisation approaches for the sparse inpainting problem, since it can also be extended to nonlinear inpainting operators such as EED. This is exploited to achieve reconstructions with state-of-the-art quality. We also give an extensive literature survey on PDE-based image compression methods.

Compensated Convexity Methods for Approximations and Interpolations of Sampled Functions in Euclidean Spaces: Applications to Contour Lines, Sparse Data and Inpainting

Surface Interpolation Scheme by Distance Blending over Convex Sets

Improved Sampling and Reconstruction in Spline Subspaces

Compensated Convexity on Bounded Domains, Mixed Moreau Envelopes and Computational Methods

Quasi-interpolators with application to postprocessing in finite element methods

Survey on Geometric Iterative Methods and Their Applications.

Optimising Spatial and Tonal Data for PDE-based Inpainting

Interpolation on Convex-set Valued Lebesgue Spaces and its Applications

Quasi-interpolation projectors for Subdivision Surfaces

First Order Methods beyond Convexity and Lipschitz Gradient Continuity with Applications to Quadratic Inverse Problems

Convexity preserving interpolatory subdivision with conic precision

Hybrid hyperinterpolation over general regions

Polyharmonic and Related Kernels on Manifolds: Interpolation and Approximation

A HYBRID ALTERNATING MINIMIZATION ALGORITHM FOR STRUCTURED CONVEX OPTIMIZATION PROBLEMS WITH APPLICATION IN POISSONIAN IMAGE PROCESSING

Superconvergent Non-Polynomial Approximations

A study of convex convex-composite functions via infimal convolution with applications

Curves defined by a class of discrete operators: approximation result and applications

A Hybrid Inpainting Model Combining Diffusion and Enhanced Exemplar Methods

Multivariate Splines and Their Applications

Medical Imaging Applications of Effectively Multi-Dimensional Interpolation