Single frame blind image deconvolution by non-negative sparse matrix factorization

Novel approach to single frame multichannel blind image deconvolution has been formulated recently as non-negative matrix factorization problem with sparseness constraints imposed on the unknown mixing vector that accounts for the case of non-sparse source image. Unlike most of the blind image deconvolution algorithms, the novel approach assumed no a priori knowledge about the blurring kernel and original image. Our contributions in this paper are: (i) we have formulated generalized non-negative matrix factorization approach to blind image deconvolution with sparseness constraints imposed on either unknown mixing vector or unknown source image; (ii) the criteria are established to distinguish whether unknown source image was sparse or not as well as to estimate appropriate sparseness constraint from degraded image itself, thus making the proposed approach completely unsupervised; (iii) an extensive experimental performance evaluation of the non-negative matrix factorization algorithm is presented on the images degraded by the blur caused by the photon sieve, out-of-focus blur with sparse and non-sparse images and blur caused by atmospheric turbulence. The algorithm is compared with the state-of-the-art single frame blind image deconvolution algorithms such as blind Richardson-Lucy algorithm and single frame multichannel independent component analysis based algorithm and non-blind image restoration algorithms such as multiplicative algebraic restoration technique and Van-Cittert algorithms. It has been experimentally demonstrated that proposed algorithm outperforms mentioned non-blind and blind image deconvolution methods.     

Application of multi-frame approach in single-frame blind deconvolution

A new single-frame blind deconvolution algorithm for the linear shift-invariant imaging system is presented. The algorithm processes the partial images segmented from one single degraded image by multi-frame approach to recover the point spread function (PSF). Then a deconvolution method is employed to restore the whole image with the recovered PSF. In addition, in order to improve the fidelity and resolution of the recovered PSF, the coprimeness of the partial images is utilized. Results of simulated and real atmospheric turbulence degraded images using the algorithm are reported.     

天文图像多帧盲反卷积收敛性的增强方法

天文图像多帧盲反卷积的收敛性受到初始目标、约束条件和光子噪声等因素的影响.提出了用实际光学成像系统参数确定频率带宽有限约束的方法.用Knox-Thompson 方法重构初始目标相位形成盲反卷积算法的初始目标函数.研究了一种新颖的有效减小光子噪声、边缘效应和振铃现象的方法.根据最大似然估计理论，用期望最大化的优化方法建立了改进的严格约束多帧盲反卷积算法.模拟图像和实际天文图像的复原结果表明，所建立的多帧盲反卷积，可以有效克服大气湍流和减小光子噪声，改善天文观察图像的分辨率，并部分消除光学系统衍射效应对恢复图像的影响. 关键词： 大气光学 天文观测 图像处理和恢复     

利用加权预测的图像迭代盲解卷积

为解决大气湍流造成的图像退化问题,本文鉴于现有的盲解卷积算法收敛性不稳定,计算量大等特点,提出了一种基于加权预测的迭代盲解卷积算法。对目前性能优秀的用迭代实现盲解卷积的L-R算法进行优化,在每次迭代结束后通过加权方法求出预测值,根据预测值计算方向加速算子,从而大大提高算法的收敛速度。实验表明：该算法不仅可对模糊退化图像进行很好的复原,同时与L-R算法相比收敛速度提高约43.8倍,其迭代速度快的特点决定了算法具有较高的工程实用价值。     

Blind image deconvolution with spatially adaptive total variation regularization

A blind deconvolution algorithm with spatially adaptive total variation regularization is introduced. The spatial information in different image regions is incorporated into regularization by using the edge indicator called difference eigenvalue to distinguish edges from flat areas. The proposed algorithm can effectively reduce the noise in flat regions as well as preserve the edge and detailed information. Moreover, it becomes more robust with the change of the regularization parameter. Comparative results on simulated and real degraded images are reported.     

Blind deconvolution of audio-frequency signals using the self-deconvolving data restoration algorithm

A signal processing algorithm has been developed in which a filter function is extracted from degraded data through mathematical operations. The filter function can be used to restore much of the degraded content of the data through use of a deconvolution process. The operation can be performed without prior knowledge of the detection system, a technique known as blind deconvolution. The extraction process, designated self-deconvolving data reconstruction algorithm, is applied here to audio-frequency signals showing significant qualitative improvement. Degradation arising from the process of electronic recording and reproduction is significantly reduced.     

Comparative study on several blind deconvolution algorithms applied to underwater image restoration

Although the use of blind deconvolution of image restoration is a widely known concept, only few reports have discussed in detail its application to solving problem of restoration of underwater range-gated laser images. A comparative study of underwater image restoration using the Richardson-Lucy algorithm, the least-squares algorithm, and the multiplicative iterative algorithm for blind deconvolution is presented. All the deconvolution approaches use denoised underwater images and Wells’ small angle approximation theory of derived point spread function as the initial object and degradation guess, respectively. Owing the underwater no-reference imaging environment, image quality judgment based on the blur metric method is incorporated in our comparison to determine the appropriate deconvolution iteration number for each algorithm, which objectively evaluates the image restoration results. The performance of the three algorithms applied to underwater image restoration is discussed and reported.     

Layered blind deconvolution using subband filter banks

A blind deconvolution method using the concept of subband filter banks has been proposed. This method makes possible the speedy convergence of solution, due to estimation of supports of the reconstructed images by a deconvolution in the lower layer. We apply the method to blind deconvolution by a simulated annealing algorithm and confirm its usefulness. Presented at 1996 International Topical Meeting on Optical Computing (OC’96), April 21–25, Sendai, Japan.     

基于噪声特性的大气湍流退化图像多帧盲反卷积复原

由于大气湍流和噪声的影响,造成观测目标图像的退化.为了目标的精确观测,根据噪声特性,结合符合物理意义的约束条件,提出了新的大气湍流图像盲反卷积复原最小化模型,并以共轭梯度数值优化方法交替迭代求解,复原观测目标图像.为验证提出的算法的有效性,在计算机上模拟参数为望远镜口径为2.0 m,大气相干长度为0.1 m,图像信噪比为10 dB的大气湍流退化和噪声污染的图像,以提出的盲反卷积复原方法复原,实验结果表明,提出的盲反卷积复原算法避免了传统的盲反卷积复原算法的缺陷,有效地克服大气湍流和噪声的影响,复原出了清晰的观测目标图像.该图像盲反卷积复原方法的研究,对地基望远镜的观测有重要的基础性作用.     

基于GPU的多帧盲解卷积图像复原技术并行化实现

多帧盲解卷积图像复原技术能够进一步提高自适应光学图像的分辨力,但其算法比较复杂,处理耗时过长,对序列图像复原经常需要几分钟甚至几十分钟的计算时间,对实际应用造成了极大不便。为了提升算法的运行速度,改善其耗时过长的问题,通过研究和分析盲解卷积算法原理和算法结构,采用目前高速发展的中央处理器（CPU）和图形处理器（GPU）异构加速技术,主要对耗时最长的矩阵卷积运算进行优化,通过使用库函数与算法结构微调相结合的方法并行加速,实现多帧盲解卷积的图像复原算法的并行化。使用并行算法对图像进行复原处理,针对16帧以上分辨率为256256像素的空间目标图像,可以实现17的加速比,为图像复原的实时/准实时提供一种可行的方案。     

利用加权预测的图像迭代盲解卷积

为解决大气湍流造成的图像退化问题,本文鉴于现有的盲解卷积算法收敛性不稳定,计算量大等特点,提出了一种基于加权预测的迭代盲解卷积算法。对目前性能优秀的用迭代实现盲解卷积的L-R算法进行优化,在每次迭代结束后通过加权方法求出预测值,根据预测值计算方向加速算子,从而大大提高算法的收敛速度。实验表明:该算法不仅可对模糊退化图像进行很好的复原,同时与L-R算法相比收敛速度提高约43.8倍,其迭代速度快的特点决定了算法具有较高的工程实用价值。     

Phase diversity imaging through atmospheric turbulence

A method of phase diversity imaging is applied to reconstruct an overly extended object from stellar speckle images. The reconstruction algorithm is based on the iterative blind deconvolution method and the Gerchberg-Saxon phase retrieval procedure. Computer simulations for reconstructing solar granulations are conducted. The simulation results exhibit the effectiveness of the phase diversity method for high-resolution solar imaging.     

基于总变分的湍流退化噪声图像盲反卷积复原

根据总变分的噪声抑制特性和大气湍流成像过程,建立了基于总变分的大气湍流噪声图像多帧盲反卷积复原最小化模型,以基于共轭梯度数值优化方法的交替迭代算法求解,复原出了观测目标的清晰图像。在计算机上模拟了湍流退化和噪声污染图像。实验结果表明,该复原算法能有效地克服大气湍流和噪声的影响,可复原出清晰的原始目标图像。     

Blind image deconvolution using the zero-lag slice of higher-order statistics

A novel algorithm for blind image deconvolution using the zero-lag slice (ZLS) of higher-order statistics only is presented. This method first estimates the point-spread function (PSF) using the ZLS of its third-order moment (TOM) and then uses it with one of the known classical image deconvolution methods. The proposed method has simple computations for PSF estimation because it solves a nonlinear problem by using an iterative method with fast convergence. In each iteration, one need only calculate the ZLS of the TOM and estimate the PSF using simple two-dimensional operations. Furthermore, the method presented achieves good results, since the ZLS estimate obtained from the degraded image exhibits high reliability. The good performance of the proposed algorithm is demonstrated by applying it to synthetic and real data sets.     

一种电磁干扰源成像多分辨率分区算法

在利用抛物反射面对电磁干扰源成像过程中,由于系统衍射受限导致干扰源成像模糊,分辨率低,难以分辨,由于不同频率不同区域干扰源所成图像分辨率不同,具有分区域多分辨率的特征,采用已有超分辨算法难以提高分辨率。利用Mean Shift算法,在原有算法基础上改进使其能够适应多分辨率的电磁干扰源成像,在图像分割的基础上对多分辨率图像进行分块抽离,并采用基于L_R迭代的盲反卷积算法分别对各区域进行分辨率的提高,仿真结果表明算法能够适应对干扰源的多分辨率电磁成像并提高分辨率。     

Restoration of solar and star images with phase diversity-based blind deconvolution

The images recorded by a ground-based telescope are often degraded by atmospheric turbulence and the aberration of the optical system. Phase diversity-based blind deconvolution is an effective post-processing method that can be used to overcome the turbulence-induced degradation. The method uses an ensemble of short-exposure images obtained imultaneously from multiple cameras to jointly estimate the object and the wavefront distribution on pupil. Based on signal estimation theory and optimization theory, we derive the cost function and solve the large-scale optimization problem using a limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method. We apply the method to the urbulence degraded images generated with computer, the solar images acquired with the swedish vacuum solar telescope (SVST, 0.475m) in La Paima and the star images collected with 1.2-m telescope in Yunnan Observatory. In order to avoid edge effect in the restoration of the solar images, a modified Hanning apodized window is adopted.The star image till can be estored when the defocus distance is measured inaccurately. The restored results demonstrate that the method is efficient for removing the effect of turbulence and reconstructing the point-like or extended objects.     

Two-stage blind deconvolution scheme using useful priors

Hand shake blurry image is a common phenomenon in our daily life. In this paper, a novel blind deconvolution scheme is proposed to recover a single hand shake blurry image. The algorithm is subdivided into two main stages, kernel estimation stage and non-blind deconvolution stage. In the kernel estimation stage, we propose a cost function taking a selected map into consideration. In the non-blind decovolution stage, another cost function is designed using image derivatives prior. We also present an adaptive kernel size selection method instead of traditional manual selection. Extensive experiments on real world blurry images are conducted to demonstrate the performance of our algorithm.     

波前解卷积方法中的高频噪声抑制

 介绍了一种波前解卷积中噪声抑制规整化的新方法，并将此方法应用于室内模拟点源实验中。该方法通过在图像复原算法中增加针对图像高频部分的限制条件来抑制高频噪声，以达到对图像复原问题病态特性的规整化。实验结果表明：该规整化方法可以有效地抑制解卷积过程中高频噪声的影响，恢复出达到理论衍射极限分辨率的图像。对于噪声水平较高的降质图像，通过这种解卷积方法可以有效地提高信噪比。同维纳逆滤波方法相比，该方法可以在有效抑制导致病态的高频噪声的基础上充分保持图像的低频；与基于贝叶斯估计的近视解卷积算法相比，该方法不需要知道噪声水平或噪声类型等先验知识，只是从噪声本质出发，通过抑制降质图像高频部分，有效地解决了病态特性问题。     

Real-time blind deconvolution of retinal images in adaptive optics scanning laser ophthalmoscopy

With the use of adaptive optics (AO), the ocular aberrations can be compensated to get high-resolution image of living human retina. However, the wavefront correction is not perfect due to the wavefront measure error and hardware restrictions. Thus, it is necessary to use a deconvolution algorithm to recover the retinal images. In this paper, a blind deconvolution technique called Incremental Wiener filter is used to restore the adaptive optics confocal scanning laser ophthalmoscope (AOSLO) images. The point-spread function (PSF) measured by wavefront sensor is only used as an initial value of our algorithm. We also realize the Incremental Wiener filter on graphics processing unit (GPU) in real-time. When the image size is 512 × 480 pixels, six iterations of our algorithm only spend about 10 ms. Retinal blood vessels as well as cells in retinal images are restored by our algorithm, and the PSFs are also revised. Retinal images with and without adaptive optics are both restored. The results show that Incremental Wiener filter reduces the noises and improve the image quality.     

自适应光学图像非对称图像迭代盲复原算法

 为了提高自适应光学图像复原效果,提出了一种新的多重约束非对称图像迭代盲解卷积算法。首先,在点扩散函数（PSF）频率域引入带宽有限约束来提高迭代盲解卷积算法的可靠性;然后,在PSF空间域引入支持域动态更新的思想以加快迭代盲解卷积算法收敛速度;最后,自动计算迭代盲解卷积算法的非对称因子以提高算法的自适应性。模拟实验结果表明,与RL-IBD算法比较,新算法迭代次数减少22.4%、峰值信噪比提高10.18 dB。在FK5-857和某双星的自适应光学图像复原实验中,也取得很好的复原效果。