English

Sparse Wavelet Representations of Spatially Varying Blurring Operators

Optimization and Control 2015-10-13 v5 Numerical Analysis

Abstract

Restoring images degraded by spatially varying blur is a problem encountered in many disciplines such as astrophysics, computer vision or biomedical imaging. One of the main challenges to perform this task is to design efficient numerical algorithms to approximate integral operators.We introduce a new method based on a sparse approximation of the blurring operator in the wavelet domain. This method requires O(Nϵd/M)\mathcal{O}\left(N \epsilon^{-d/M}\right) operations to provide ϵ\epsilon-approximations, where NN is the number of pixels of a dd-dimensional image and M1M\geq 1 is a scalar describing the regularity of the blur kernel. In addition, we propose original methods to define sparsity patterns when only the operators regularity is known.Numerical experiments reveal that our algorithm provides a significant improvement compared to standard methods based on windowed convolutions.

Keywords

Cite

@article{arxiv.1404.1023,
  title  = {Sparse Wavelet Representations of Spatially Varying Blurring Operators},
  author = {Paul Escande and Pierre Weiss},
  journal= {arXiv preprint arXiv:1404.1023},
  year   = {2015}
}
R2 v1 2026-06-22T03:42:34.315Z