摘要
不论是在日常摄影还是在高分辨率遥感对地观测中,照明条件不佳都将意味着曝光时间的增长,随之而来的问题便是曝光过程中由于手抖动或平台颤振造成的图像模糊,大大损失了图像包含的信息。在数学上,图像模糊可以表述为清晰图像与点扩散函数(PSF)的卷积,再叠加一定噪声的结果,其逆过程称为图像复原或图像反卷积,是一个典型的病态问题,结果极易受到噪声放大和振铃等负面效应的影响,而正则化算法则能在一定程度上克服这些负面效应,提高复原图像的质量。图像复原算法有着极高的应用价值,尤其是在高分辨率遥感对地观测中,可以有效的避免由于设计制造大口径、长焦距成像系统以及相应的稳像设备而带来的成本急剧上升。为此,本文将图像复原算法作为研究重点,详细讨论了与其相关的若干关键问题,概括如下。
首先是PSF的获取。PSF在图像复原问题中扮演着至关重要的角色,传统的获取PSF的方法多种多样,比较经典的有点光源法和刃边法,以及S.K. Navar和M. Ben—Ezra提出的基于运动轨迹的PSF重构方法。针对颤振运动的特点,本文提出了一种基于光纤陀螺颤振探测的PSF重构方法。光纤陀螺是一种角速度敏感器件,具有结构紧凑,探测精度高的特点,将两个敏感轴互相垂直的光纤陀螺与相机光轴进行刚性连接,就能获取相机在曝光时间内的颤振角位置数据,进而根据物像关系得到像面上像的运动轨迹,重构出PSF。此外,本文还对传统的基于运动轨迹的PSF重构方法做了改进,利用设置较高采样频率的方法,将复杂的几何和插值运算转化为简单的统计运算,有效提高了运行效率。
按照PSF是否已知,可将图像复原算法分为非盲复原算法和盲复原算法两大类,本文在前人工作的基础上提出了三种正则化的非盲复原算法,它们分别是分段局部正则化的Richardson—Lucy(RL)算法,基于自然图像梯度先验分布的正则化RL算法和基于混合高斯型马尔可夫专家场(GSM FoE)的非盲复原算法。前两种算法都基于泊松噪声模型,是对标准RL算法的改进。在分段局部正则化的RL算法中,本文设计了一种基于高斯马尔可夫随机场的正则项,并通过一个分段幂函数,对其平滑强度加以调整,从而达到了抑制复原图像中的噪声放大和振铃效应的目的。与其相似,基于自然图像梯度先验分布的正则化RL算法是将自然图像梯度先验分布模型与泊松噪声模型相结合,由于自然图像梯度先验分布模型很好的符合了真实自然图像梯度的稀疏分布规律,因此该算法能够有效的改善标准RL算法的表现,得到高质量的复原图像。GSM FoE是在马尔可夫专家场(FoE)基础上提出的一种新型自然图像概率模型,与传统模型相比,其优势在于构成它的全部滤波器都是通过特殊的最优化方法,采用自然图像库训练得到的,因此它能够对自然图像进行更加准确的建模。本文将其引入非盲复原算法,并采用了Split Bregman算法对问题进行求解,实验结果表明,它是一种优秀的图像复原算法,达到或超过了目前一些主流复原算法的水平。
在盲复原算法方面,本文同样提出了三种方法,它们是基于马尔可夫专家场(FoE)的单幅图像盲复原算法,基于自然图像梯度先验分布的单幅图像盲复原算法和基于多幅图像的盲复原算法。其中,基于马尔可夫专家场的单幅图像盲复原算法的创新主要有两点,首先是将FoE模型引入盲复原算法,与GSM FoE相似,构成它的全部参数和滤波器都是由训练得到的,因此该模型的准确度很高。同时,还采用了基于Student-t函数的概率模型对PSF进行建模。其次,在问题求解上,本文改进了传统的轮换迭代算法,在每次迭代中,都将上次迭代所得的结果作为约束条件,因此克服了传统方法容易收敛到模糊解的缺点,得到了高质量的复原图像。基于自然图像梯度先验分布的单幅图像盲复原算法建立在基于自然图像梯度先验分布的正则化RL算法基础上,并采用了传统的轮换迭代算法求解,其不同之处在于对PSF进行优化时,用高斯概率模型对泊松概率模型做了近似,因而提高了运算速度和结果的准确度。最后,本文还提出了一种基于多幅图像的盲复原算法,它分别采用了自然图像梯度先验分布模型和基于L1范数的概率模型对清晰图像和PSF进行建模,由于多幅模糊图像包含了更多的信息,因此该算法能够得到高质量的复原图像和较为准确的PSF。
最后,本文还搭建了相应的实验系统对上述PSF重构方法和各种复原算法进行了验证,实验结果表明本文提出的基于光纤陀螺颤振探测的图像复原方法具有很强鲁棒性,能够对模糊图像实现有效复原,显著提高图像质量。同时,本文提出的几种图像复原算法也有较好的表现,其中基于混合高斯型马尔可夫专家场的非盲复原算法和基于马尔可夫专家场盲复原算法达到了或超过了目前一些主流图像复原算法的水平。
No matter in daily photography or high-resolution remote sensing, bad lighting condition means the increase of exposure time, the accompanied problem is, during the exposure time, the hand shake or platform vibration results in blurred images, which will greatly reduce the amount of information contained in the images. In mathematics, image blurring is modeled by convolving the clear image with a point spread function (PSF) plus some noise. Its inverse process is called image deconvolution, which is a typical ill-posed problem, the result is usually contaminated by amplified noise and ringing artifacts. However, regularization methods can be used to improve the performance of image deconvolution algorithm to achieve result of high quality. Image deconvolution algorithm is very useful, especially in the application of high-resolution remote sensing, which can effectively avoid the sharp rise of the cost caused by the design and manufacture of imaging systems with large aperture, long focal length and the corresponding stabilization equipments. Therefore, the thesis focuses on the image deconvolution algorithm and detailed discusses some key problems about it, which are summarized as follows.
The first is about how to obtain the PSF. PSF plays a very important role in image deconvolution. There are various methods which can be adopted to achieve it, such as the point lighting source method, the knife edge method and the method based on motion trajectory which is proposed by S.K.Nayar and M.Ben-Ezra. According to the characteristics of vibration, we propose a PSF reconstruction method using the fiber optic gyroscope. The fiber optic gyroscope is a kind of angular velocity sensing device, which has the advantages of compact structure and high detection accuracy. We rigidly connect two fiber optic gyroscopes whose sensitive axes are perpendicular to each other to the optical axis of the camera, such that we can get the angle position of the vibration in two directions during the exposure time, then according to the relationship between the object and the image, we can obtain the motion trajectory of the image and reconstruct the PSF. In addition, we also improve the traditional PSF reconstruction method based on motion trajectory. In the improved version, with the method of setting a high sampling frequency, the complicated geometric and interpolation calculations are converted into a simple statistical calculation, and thus enhance the efficiency.
According to whether the PSF is known, image deconvolution algorithms can be divided into two kinds, i.e., non-blind and blind image deconvolution. The thesis propose three non-blind image deconvolution algorithms, they are the piecewise local regularized Richardson-Lucy (RL) algorithm, the regularized RL algorithm based on natural image gradient prior and the algorithm based on Gaussian Scale Mixture Fields of Experts (GSM FoE) prior. The first two are derived from the model of Poisson noise, which are all improved versions of standard RL algorithm. In the piecewise local regularized RL algorithm, we design a regularization term based on the Gaussian Markov random field, and adopt a piecewise power function to adjust the smoothing strength to suppress the amplified noise and ringing artifacts in the restored image. Similarly, in the regularized RL algorithm based on natural image gradient prior, the model of Poisson noise is combined with the prior of image gradient, because the prior is in good agreement with the sparse probabilistic distribution of the gradient of nature images, the algorithm can efficiently improve the performance of standard RL algorithm and reach result of high quality. GSM FoE which is derived from Fields of Experts (FoE) is a very effective probabilistic model for natural image, in contrast to traditional priors, its advantage is that all the filters used to construct the GSM FoE prior are trained with a database, thus it can seize the characteristics of natural images more accurately. The thesis introduces this prior into non-blind image deconvolution and adopts the Split Bregman method to solve the resulted cost function. Experimental results show that its performance meets or exceeds some state of the art methods.
The thesis also proposes three blind image deconvolution algorithms, i.e., the single image deconvolution algorithm based on Fields of Experts (FoE) prior, the single image deconvolution algorithm based on natural image gradient prior and the deconvolution algorithm with multiple blurred images. There are two innovations in the single image deconvolution algorithm based on FoE prior. Firstly, we introduce the FoE prior into blind image deconvolution. similar to the GSM FoE prior, all its parameters and filters are trained from natural images, thus it is of high accuracy. Meanwhile, a prior based on Student-t function is also used to regularize the PSF.
Secondly, we improve the traditional alternating minimization (AM) algorithm which is usually used to solve the problem of blind image deconvolution. In each iteration of the improved algorithm, the mid-restored image obtained from the former iteration is used as a constraint, and thus it will not converge to the unwanted "blurry" result. The single image deconvolution algorithm based on natural image gradient prior adopts the proposed regularized RL algorithm with natural image gradient prior. The result is reached by the traditional AM algorithm. The innovation is that the Gaussian noise model is used to approximate the Poisson noise model to estimate the PSF, so the computation speed and accuracy are enhanced. We also propose a blind image deconvolution algorithm adopting multiple blurred images, the natural image gradient prior and a L1norm based prior is used to regularize the clear image and the PSF respectively. Since multiple frames contain more information about the clear image, it can achieve results of high quality.
Finally, we also set up experimental system to verify the effectiveness of the proposed PSF reconstruction method and the image deconvolution algorithms. The results show that the image deconvolution method based on vibration detection using fiber optic gyroscope is robust, it can effectively enhance the quality of the obtained image. Meanwhile, the performances of the deconvolution algorithms proposed in this thesis are also good. Especially, the blind image deconvolution algorithm based on the FoE prior and the non-blind image deconvolution algorithm based on the GSM FoE prior are comparable with some state of the art methods.
引文
[1]K. Janschek, V. Tchernykh, Optical correlator for image motion compensation in the focal plane of a satellite camera, Automatic Control in Aerospace,2002, pp.357-361.
[2]S.K. Nayar, M. Ben-Ezra, Motion-based motion deblurring, IEEE Trans. on Pattern Analysis and Machine Intelligence,26 (2004) 689-698.
[3]W.D. Dong, Y.T. Chen, Z.H. Xu, H.J. Feng, Q. Li, Image stabilization with support vector machine, Journal of Zhejiang University-Science C-Computers & Electronics,12 (2011) 478-485.
[4]H.X. Pang, W.D. Dong, Z.H. Xu, H.J. Feng, Q. Li, Y.T. Chen, Novel linear search for support vector machine parameter selection, Journal of Zhejiang University-Science C-Computers & Electronics,12 (2011) 885-896.
[5]J. Lewis, Fast template matching, Vision Interface,1995, pp.15-19.
[6]B.S. Reddy, B. Chatterji, An FFT-based technique for translation, rotation, and scale-invariant image registration, IEEE Trans. on Image Process.,5 (1996) 1266-1271.
[7]L.G. Brown, A survey of image registration techniques, ACM computing surveys (CSUR), 24(1992)325-376.
[8]R.C. Gonzalez, R.E. Woods, Digital image processing, Second Edition, Publishing House of Electronics Industry,2002.
[9]W.H. Richardson, Bayesian-Based Iterative Method of Image Restoration, J. Opt. Soc. Am., 62(1972)55-59.
[10]L.B. Lucy, An iterative technique for the rectification of observed distributions, Astronomical Journal,79 (1974) 745-754.
[11]E. Bratsolis, M. Sigelle, A spatial regularization method preserving local photometry for Richardson-Lucy restoration, Astronomy and Astrophysics,375 (2001) 1120-1128.
[12]J.M. Bioucas-Dias, M.A.T. Figueiredo, J.P. Oliveira, Total Variation-Based Image Deconvolution:a Majorization-Minimization Approach, Proceedings of IEEE International Conference on Acoustics, Speech and Signal Processing, IEEE,2006, pp.861-864.
[13]P. Rodriguez, B. Wohlberg, Efficient Minimization Method for a Generalized Total Variation Functional, IEEE Trans. on Image Process.,18 (2009) 322-332.
[14]Y.L. Wang, J.F. Yang, W.T. Yin, Y. Zhang, A New Alternating Minimization Algorithm for Total Variation Image Reconstruction, SIAM Journal on Imaging Sciences,1 (2008) 248-272.
[15]T. Goldstein, S. Osher, The Split Bregman Method for L1-Regularized Problems, SIAM Journal on Imaging Sciences,2 (2009) 323-343.
[16]A. Levin, R. Fergus, F. Durand, W.T. Freeman, Image and depth from a conventional camera with a coded aperture, Acm Trans. on Graphics,26 (2007).
[17]D. Krishnan, R. Fergus, Fast image deconvolution using hyper-laplacian priors, Advances in Neural Information Processing Systems,22 (2009) 1-9.
[18]L. Yuan, J. Sun, L. Quan, H.Y. Shum, Progressive inter-scale and intra-scale non-blind image deconvolution, Acm Trans. on Graphics,27 (2008).
[19]S.G. Chang, B. Yu, M. Vetterli, Adaptive wavelet thresholding for image denoising and compression, IEEE Trans. on Image Process.,9 (2000) 1532-1546.
[20]S.G. Chang, B. Yu, M. Vetterli, Spatially adaptive wavelet thresholding with context modeling for image denoising, IEEE Trans. on Image Process.,9 (2000) 1522-1531.
[21]J. Portilla, V. Strela, M.J. Wainwright, E.P. Simoncelli, Image denoising using scale mixtures of Gaussians in the wavelet domain, IEEE Trans. on Image Process.,12 (2003) 1338-1351.
[22]R. Neelamani, H. Choi, R. Baraniuk, Forward:Fourier-Wavelet regularized deconvolution for Ill-conditioned systems, IEEE Trans. on Signal Processing,52 (2004) 418-433.
[23]J.M. Bioucas-Dias, Bayesian wavelet-based image deconvolution:A GEM algorithm exploiting a class of heavy-tailed priors. IEEE Trans. on Image Process.,15 (2006) 937-951.
[24]D. Weisheng, Z. Lei, S. Guangming, W. Xiaolin, Image Deblurring and Super-Resolution by Adaptive Sparse Domain Selection and Adaptive Regularization, IEEE Trans. on Image Process.,20 (2011) 1838-1857.
[25]V. Katkovnik, K. Egiazarian, J. Astola, A spatially adaptive nonparametric regression image deblurring, IEEE Trans. on Image Process.,14 (2005) 1469-1478.
[26]G. Chantas, N. Galatsanos, A. Likas. M. Saunders, Variational Bayesian image restoration based on a product of t-distributions image prior, IEEE Trans. on Image Process.,17 (2008) 1795-1805.
[27]J. Ni, P. Turaga, V.M. Patel, R. Chellappa, Example-driven manifold priors for image deconvolution, IEEE Trans. on Image Process.,20 (2011) 3086-3096.
[28]T. Hou, S. Wang, H. Qin, Image deconvolution with multi-stage convex relaxation and its perceptual evaluation, IEEE Trans. on Image Process.,20 (2011) 3383-3392.
[29]G. Chantas, N.P. Galatsanos, R. Molina, A.K. Katsaggelos, Variational Bayesian image restoration with a product of spatially weighted total variation image priors, IEEE Trans, on Image Process.,19 (2010) 351-362.
[30]R. Pan, S.J. Reeves, Efficient Huber-Markov edge-preserving image restoration, IEEE Trans, on Image Process.,15 (2006) 3728-3735.
[31]F. Viallefont-Robinet, Edge method for on-orbit defocus assessment, Opt. Express,18 (2010)20845-20851.
[32]F. Viallefont-Robinet, D. Leger, Improvement of the edge method for on-orbit MTF measurement, Opt. Express,18 (2010) 3531-3545.
[33]H. Hwang, Y.-W. Choi, S. Kwak, M. Kim, W. Park, MTF assessment of high resolution satellite images using ISO 12233 slanted-edge method, International Society for Optics and Photonics,2008, pp.710905-710909.
[34]D. Leger, F. Viallefont, E. Hillairet, A. Meygret, In-flight refocusing and MTF assessment of SPOT5 HRG and HRS cameras, International Symposium on Remote sensing,2003, pp. 224-231.
[35]P. Dherete, B. Rouge, Image de-blurring and application to SPOT5 THR satellite imaging, Proceedings of Geoscience and Remote Sensing Symposium, IEEE,2003, pp.318-320.
[36]K. Kohm, Modulation transfer function measurement method and results for the Orbview-3 high resolution imaging satellite, Proceedings of ISPRS,2004. pp.12-23.
[37]T. Choi, IKONOS satellite on orbit modulation transfer function (MTF) measurement using edge and pulse method, A thesis submitted in partial fulfillment of the requirements for the master of Science, Major in Engineering. South Dakota State University,2002.
[38]P. Schwind, M. Schneider, G. Palubinskas. T. Storch, R. Muller, R. Richter, Processors for ALOS optical data:Deconvolution, DEM generation, orthorectification, and atmospheric correction, IEEE Trans. on Geoscience and Remote Sensing,47 (2009) 4074-4082.
[39]P. Schwind, G. Palubinskas, T. Storch, R. Miiller, Evaluation of deconvolution methods for PRISM images, Proceedings of ALOS PI Symposium,2008.
[40]付中梁,冯华君,徐之海,李奇,毛成军,基于快速CCD位移探测的运动模糊图像的恢复,光电工程,36(2009)69-73.
[41]付中梁,冯华君,徐之海,李奇,陈跃庭,二维任意运动形式模糊图像的恢复,光学学报,2(2010)369-372.
[42]W. Stefan, Image Restoration by blind deconvolution (Diploma Thesis), Mathematics, Technische Universitat Munchen,2003.
[43]I. Kopriva, D.J. Garrood, V. Borjanovic, Single frame blind image deconvolution by non-negative sparse matrix factorization, Optics Communications,266 (2006) 456-464.
[44]D.A. Fish, A.M. Brinicombe, E.R. Pike, J.G. Walker, Blind deconvolution by means of the Richardson-Lucy algorithm, J. Opt. Soc. Am. A,12 (1995) 58-65.
[45]J. Zhang, Q. Zhang, G. He, Blind image deconvolution by means of asymmetric multiplicative iterative algorithm, J. Opt. Soc. Am. A,25 (2008) 710-717.
[46]T.F. Chan, W. Chiu-Kwong, Total variation blind deconvolution, IEEE Trans. on Image Process.,7(1998)370-375.
[47]Y.L. You, M. Kaveh, Blind image restoration by anisotropic regularization, IEEE Trans. on Image Process.,8 (1999) 396-407.
[48]R. Fergus, B. Singh, A. Hertzmann, S.T. Roweis, W.T. Freeman, Removing camera shake from a single photograph, Acm Trans. on Graphics,25 (2006) 787-794.
[49]J. Miskin, D.J.C. MacKay, Ensemble learning for blind image separation and deconvolution, Advances in Independent Component Analysis.2000, pp.123-141.
[50]Q. Shan, J. Jia, A. Agarwala, High-quality motion deblurring from a single image, Acm Trans, on Graphics,27 (2008).
[51]D.L. Donoho, Compressed sensing, IEEE Trans. on Information Theory.,52 (2006) 1289-1306.
[52]D.L. Donoho, M. Elad, Optimally sparse representation in general (nonorthogonal) dictionaries via L1 minimization, Proceedings of the National Academy of Sciences of the United States of America,100 (2003) 2197-2202.
[53]D.L. Donoho, M. Elad, V.N. Temlyakov, Stable recovery of sparse overcomplete representations in the presence of noise, IEEE Trans. on Information Theory.,52 (2006) 6-18.
[54]M.S.C. Almeida, L.B. Almeida, Blind and Semi-Blind Deblurring of Natural Images, IEEE Trans. on Image Process.,19 (2010) 36-52.
[55]Z. Xu, E.Y. Lam, Maximum a posteriori blind image deconvolution with Huber-Markov random-field regularization, Opt. Lett.,34 (2009) 1453-1455.
[56]A.P. Dempster, N.M. Laird, D.B. Rubin, Maximum likelihood from incomplete data via the EM algorithm, Journal of the Royal Statistical Society. Series B (Methodological), (1977) 1-38.
[57]D. Krishnan, T. Tay, R. Fergus, Blind deconvolution using a normalized sparsity measure, Proceedings of IEEE Conference on Computer Vision and Pattern Recognition, IEEE,2011, pp. 233-240.
[58]D.G. Tzikas, A.C. Likas, N.P. Galatsanos, Variational Bayesian Sparse Kernel-Based Blind Image Deconvolution With Student's-t Priors, IEEE Trans. on Image Process.,18 (2009) 753-764.
[59]A.C. Likas, N.P. Galatsanos, A variational approach for Bayesian blind image deconvolution, IEEE Trans, on Signal Processing,52 (2004) 2222-2233.
[60]H. Liao, M.K. Ng, Blind deconvolution using generalized cross-validation approach to regularization parameter estimation, IEEE Trans. on Image Process.,20 (2011) 670-680.
[61]S. Unnikrishna Pillai, B. Liang, Blind image deconvolution using a robust GCD approach, IEEE Trans. on Image Process.,8 (1999) 295-301.
[62]M.K. Ng, R.J. Plemmons, S. Qiao, Regularization of RIF blind image deconvolution, IEEE Trans. on Image Process.,9 (2000) 1130-1134.
[63]S.D. Babacan, R. Molina, A.K. Katsaggelos, Variational Bayesian blind deconvolution using a total variation prior, IEEE Trans. on Image Process.,18 (2009) 12-26.
[64]D. Tzikas, A. Likas, N. Galatsanos, Variational Bayesian blind image deconvolution with student-t priors, IEEE International Conference on Image Processing, IEEE,2007. pp. I-109-I-112.
[65]L. Yan, H. Fang, S. Zhong, Blind image deconvolution with spatially adaptive total variation regularization, Opt. Lett.,37 (2012) 2778-2780.
[66]G. Vicidomini, S.W. Hell, A. Schonle, Automatic deconvolution of 4Pi-microscopy data with arbitrary phase, Opt. Lett.,34 (2009) 3583-3585.
[67]S. Cho, S. Lee, Fast Motion Deblurring, Acm Trans. on Graphics,28 (2009).
[68]J.H. Money, S.H. Kang, Total variation minimizing blind deconvolution with shock filter reference, Image and Vision Computing,26 (2008) 302-314.
[69]L. Xu, J.Y. Jia, Two-Phase Kernel Estimation for Robust Motion Deblurring, Proceedings of European Conference on Computer Vision, ECCV,2010, pp.157-170.
[70]A. Levin, Y. Weiss, F. Durand. W.T. Freeman, Understanding and evaluating blind deconvolution algorithms, Proceedings of IEEE Conference on Computer Vision and Pattern Recognition, IEEE,2009, pp.1964-1971.
[71]A. Levin, Y. Weiss, F. Durand, W.T. Freeman, Efficient marginal likelihood optimization in blind deconvolution, IEEE Conference on Computer Vision and Pattern Recognition, IEEE, 2011, pp.2657-2664.
[72]L. Yuan, J. Sun, L. Quan, H.Y. Shum, Image deblurring with blurred/noisy image pairs, Acm Trans. on Graphics,26 (2007).
[73]S.H. Lee, H.M. Park, S.Y. Hwang, Motion deblurring using edge map with blurred/noisy image pairs, Optics Communications,285 (2012) 1777-1786.
[74]J.-F. Cai, H. Ji, C. Liu, Z. Shen, High-quality curvelet-based motion deblurring from an image pair, IEEE Conference on Computer Vision and Pattern Recognition, IEEE.2009, pp. 1566-1573.
[75]A. Rav-Acha, S. Peleg. Two motion-blurred images are better than one, Pattern Recognition Letters,26 (2005) 311-317.
[76]C. Jia, Y. Lu, T. Chi-Keung, Q. Long, Robust dual motion deblurring. Proceedings of IEEE Conference on Computer Vision and Pattern Recognition, IEEE,2008, pp.1-8.
[77]C.A. Berenstein, E.V. Patrick, Exact deconvolution for multiple convolution operators-an overview, plus performance characterizations for imaging sensors. Proceedings of the IEEE, 1990, pp.723-734.
[78]L.P. Yaroslavsky, H.J. Caulfield, Deconvolution of multiple images of the same object, Appl. Opt.,33 (1994) 2157-2162.
[79]F. Sroubek, J. Flusser, Multichannel blind iterative image restoration, IEEE Trans. on Image Process.,12 (2003) 1094-1106.
[80]F. Sroubek, J. Flusser, Multichannel blind deconvolution of spatially misaligned images, IEEE Trans. on Image Process.,14 (2005) 874-883.
[81]J.-F. Cai, H. Ji, C. Liu, Z. Shen, Blind motion deblurring using multiple images, Journal of computational physics,228 (2009) 5057-5071.
[82]Y.-W. Wen, C. Liu, A.M. Yip, Fast splitting algorithm for multiframe total variation blind video deconvolution, Appl. Opt.,49 (2010) 2761-2768.
[83]D.F. Shi, C.Y. Fan, H. Shen, P.F. Zhang, C.H. Qiao, X.X. Feng, Y.J. Wang, Restoration of atmospheric turbulence degraded images, Optics Communications,284 (2011) 5556-5561.
[84]T. Schulz, B. Stribling, J. Miller, Multiframe blind deconvolution with real data:imagery of the Hubble Space Telescope, Opt. Express,1 (1997) 355-362.
[85]D.A. Hope, S.M. Jefferies, Compact multiframe blind deconvolution, Opt. Lett.,36 (2011) 867-869.
[86]Y.V. Zhulina, Multiframe blind deconvolution of heavily blurred astronomical images, Appl. Opt.,45 (2006) 7342-7352.
[87]M. Piana, M. Bertero, Regularized deconvolution of multiple images of the same object, J. Opt. Soc. Am. A.13(1996) 1516-1523.
[88]M.F. Tappen, B.C. Russell, W.T. Freeman, Exploiting the sparse derivative prior for super-resolution and image demosaicing, Proceedings of IEEE Workshop on Statistical and Computational Theories of Vision, Citeseer,2003.
[89]A. Levin, Y. Weiss, User Assisted Separation of Reflections from a Single Image Using a Sparsity Prior, IEEE Trans. on Pattern Analysis and Machine Intelligence,29 (2007) 1647-1654.
[90]邹谋炎,反卷积和信号复原,国防工业出版社,2001.
[91]M. Bertero, P. Boccacci, Introduction to inverse problems in imaging. Institute of Physics Publishing,1998.
[92]A.N. Tikhonov, A. Goncharsky, V. Stepanov, A. Yagola, Numerical methods for the solution of ill-posed problems, Kluwer Academic Publishers Dordrecht,1995.
[93]A.N. Tikhonov, On the stability of inverse problems, Proceedings of Doklady Akademii Nauk SSSR, Russian Academy of Sciences,1943, pp.195-198.
[94]L.I. Rudin, S. Osher, E. Fatemi, Nonlinear total variation based noise removal algorithms, Physica D,60 (1992) 259-268.
[95]S. Geman, D. Geman, Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images, IEEE Trans. on Pattern Analysis and Machine Intelligence, (1984) 721-741.
[96]D. Geman, Y. Chengda, Nonlinear image recovery with half-quadratic regularization, IEEE Trans. on Image Process.,4 (1995) 932-946.
[97]P. Charbonnier, L. Blanc-Feraud, G. Aubert, M. Barlaud, Deterministic edge-preserving regularization in computed imaging, IEEE Trans. on Image Process.,6 (1997) 298-311.
[98]J. Moussouris, Gibbs and Markov random systems with constraints, Journal of Statistical Physics,10(1974)11-33.
[99]M. Bertero, P. Boccacci, G. Desidera, G. Vicidomini, Image deblurring with Poisson data: from cells to galaxies, Inverse Problems,25 (2009) 123006.
[100]P.J. Green, On use of the EM for penalized likelihood estimation, Journal of the Royal Statistical Society. Series B (Methodological), (1990) 443-452.
[101]P.J. Green, Bayesian reconstructions from emission tomography data using a modified EM algorithm, IEEE Trans. on Medical Imaging,9 (1990) 84-93.
[102]R. Anderson, H. Bilger, G. Stedman, " Sagnac" effect:A century of Earth-rotated interferometers, American journal of physics,62 (1994) 975-985.
[103]张桂才,光纤陀螺原理与技术.国防工业出版社,2008.
[104]N. Dey, L. Blanc-Feraud, C. Zimmer, P. Roux, Z. Kam, J.C. Olivo-Marin, J. Zerubia, Richardson-Lucy algorithm with total variation regularization for 3D confocal microscope deconvolution, Microscopy Research and Technique,69 (2006) 260-266.
[105]Z. Hao, L. Yu, W. Qinzhang, Blind image deconvolution subject to bandwidth and total variation constraints. Opt. Lett.,32 (2007) 2550-2552.
[106]W. Zhou, A.C. Bovik, H.R. Sheikh, E.P. Simoncelli, Image quality assessment:from error visibility to structural similarity, IEEE Trans. on Image Process.,13 (2004) 600-612.
[107]W.D. Dong, H.J. Feng, Z.H. Xu, Q. Li, A piecewise local regularized Richardson-Lucy algorithm for remote sensing image deconvolution, Optics and Laser Technology,43 (2011) 926-933.
[108]Y.P. Wang, H.J. Feng, Z.H. Xu, Q. Li, C.Y. Dai, An improved Richardson-Lucy algorithm based on local prior, Optics and Laser Technology,42 (2010) 845-849.
[109]J. Zhao, H. Feng, Z. Xu, Q. Li, An improved image deconvolution approach using local constraint, Optics and Laser Technology,44 (2012) 421-427.
[110]S. Roth, M.J. Black, Fields of Experts:a framework for learning image priors, Proceedings of IEEE Conference on Computer Vision and Pattern Recognition, IEEE,2005, pp. 860-867.
[111]S. Roth, M.J. Black, Fields of experts, International Journal of Computer Vision,82 (2009) 205-229.
[112]Y. Weiss, W.T. Freeman, What makes a good model of natural images?, IEEE Conference on Computer Vision and Pattern Recognition, IEEE,2007, pp.1-8.
[113]M. Elad, Sparse and redundant representations, Springer,2010.
[114]C. Wang, L.F. Sun, Z.Y Chen, J.W. Zhang, S.Q. Yang, Multi-scale blind motion deblurring using local minimum, Inverse Problems,26 (2010).
[115]D. Sun, W.-K. Cham, Postprocessing of low bit-rate block DCT coded images based on a fields of experts prior, IEEE Trans. on Image Process.,16 (2007) 2743-2751.
[116]T. Scavo, J. Thoo, On the geometry of Halley's method, The American mathematical monthly,102(1995)417-426.