The blind image deconvolution is a challenging, highly ill-posed nonlinear inverse problem. We introduce a Multiscale Hierarchical Decomposition Method (MHDM) that is iteratively solving variational problems with adaptive data and regularization parameters, towards obtaining finer and finer details of the unknown kernel and image. We establish convergence of the residual in the noise-free data case, and then in the noisy data case when the algorithm is stopped early by means of a discrepancy principle. Fractional Sobolev norms are employed as regularizers for both kernel and image, with the advantage of computing the minimizers explicitly in a pointwise manner. In order to break the notorious symmetry occurring during each minimization step, we enforce a positivity constraint on the Fourier transform of the kernels. Numerical comparisons with a single-step variational method and a non-blind MHDM show that our approach produces comparable results, while less laborious parameter tuning is necessary at the price of more computations. Additionally, the scale decomposition of both reconstructed kernel and image provides a meaningful interpretation of the involved iteration steps.
翻译:盲图像去卷积是一个具有挑战性的高度病态非线性逆问题。本文提出一种多尺度层次分解方法,该方法通过迭代求解具有自适应数据和正则化参数的变分问题,逐步获取未知核与图像的更精细细节。我们在无噪声数据情形下建立了残差的收敛性,进而在有噪声数据情形下通过差异原理提前终止算法时证明了收敛性。采用分数阶Sobolev范数作为核与图像的正则项,其优势在于能够以逐点方式显式计算极小化子。为打破每次极小化步骤中出现的典型对称性问题,我们对核的傅里叶变换施加了非负约束。与单步变分方法及非盲多尺度层次分解方法的数值比较表明,本方法在获得相当结果的同时,虽需更多计算量,但显著减少了繁琐的参数调节工作。此外,重建核与图像的尺度分解为迭代步骤提供了有意义的解释。