Blind Hierarchical Deconvolution

Deconvolution is a fundamental inverse problem in signal processing and the prototypical model for recovering a signal from its noisy measurement. Nevertheless, the majority of model-based inversion techniques require knowledge on the convolution kernel to recover an accurate reconstruction and additionally prior assumptions on the regularity of the signal are needed. To overcome these limitations, we parametrise the convolution kernel and prior length-scales, which are then jointly estimated in the inversion procedure. The proposed framework of blind hierarchical deconvolution enables accurate reconstructions of functions with varying regularity and unknown kernel size and can be solved efficiently with an empirical Bayes two-step procedure, where hyperparameters are first estimated by optimisation and other unknowns then by an analytical formula.

[1]  Christopher K. I. Williams,et al.  Gaussian Processes for Machine Learning (Adaptive Computation and Machine Learning) , 2005 .

[2]  Neil D. Lawrence,et al.  Deep Gaussian Processes , 2012, AISTATS.

[3]  Christopher J Paciorek,et al.  Spatial modelling using a new class of nonstationary covariance functions , 2006, Environmetrics.

[4]  M. Girolami,et al.  Hyperpriors for Matérn fields with applications in Bayesian inversion , 2016, Inverse Problems & Imaging.

[5]  S. Siltanen,et al.  Can one use total variation prior for edge-preserving Bayesian inversion? , 2004 .

[6]  Samuli Siltanen,et al.  Linear and Nonlinear Inverse Problems with Practical Applications , 2012, Computational science and engineering.

[7]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[8]  Jean-Christophe Pesquet,et al.  Euclid in a Taxicab: Sparse Blind Deconvolution with Smoothed ${\ell _1}/{\ell _2}$ Regularization , 2014, IEEE Signal Processing Letters.

[9]  Lassi Roininen,et al.  Cauchy difference priors for edge-preserving Bayesian inversion , 2016, Journal of Inverse and Ill-posed Problems.

[10]  Andrew M. Stuart,et al.  How Deep Are Deep Gaussian Processes? , 2017, J. Mach. Learn. Res..

[11]  Mark Girolami,et al.  Posterior inference for sparse hierarchical non-stationary models , 2018, Comput. Stat. Data Anal..

[12]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.