Sparse Bayesian blind image deconvolution with parameter estimation

EURASIP Journal on Image and Video Processing, Nov 2012

In this article, we propose a novel blind image deconvolution method developed within the Bayesian framework. We concentrate on the restoration of blurred photographs taken by commercial cameras to show its effectiveness. The proposed method is based on a non-convex l p quasi norm with 0<p<1 that is used for the image, and a total variation (TV) based prior that is utilized for the blur. Bayesian inference is carried out by utilizing bounds for both the image and blur priors using a majorization-minimization principle. Maximum a posteriori estimates of the unknown image, blur and model parameters are calculated. Experimental results (i.e., restorations of more than 30 blurred photographs) are presented to demonstrate the advantage of the proposed method compared to existing ones.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1186%2F1687-5281-2012-20.pdf

Sparse Bayesian blind image deconvolution with parameter estimation

EURASIP Journal on Image and Video Processing Sparse Bayesian blind image deconvolution with parameter estimation Bruno Amizic 0 Rafael Molina 1 Aggelos K Katsaggelos 0 0 Department of Electrical Engineering and Computer Science, Northwestern University , Evanston, IL, 60208-3118 , USA 1 Departamento de Ciencias de la Computacio ́ n e I.A., Universidad de Granada , 18071 Granada , Spain In this article, we propose a novel blind image deconvolution method developed within the Bayesian framework. We concentrate on the restoration of blurred photographs taken by commercial cameras to show its effectiveness. The proposed method is based on a non-convex lp quasi norm with 0 < p < 1 that is used for the image, and a total variation (TV) based prior that is utilized for the blur. Bayesian inference is carried out by utilizing bounds for both the image and blur priors using a majorization-minimization principle. Maximum a posteriori estimates of the unknown image, blur and model parameters are calculated. Experimental results (i.e., restorations of more than 30 blurred photographs) are presented to demonstrate the advantage of the proposed method compared to existing ones. 1 Introduction Blind image deconvolution (BID) refers to the process of estimating both the original image and the blur from the degraded noisy image observation by using partial information about the imaging system. Blind image deconvolution algorithms represent a valuable tool that can be used for improving image quality without requiring complicated calibrations of the real-time image acquisition and processing system (i.e., medical imaging, videoconferencing, space exploration, x-ray imaging, etc.). The blind image deconvolution problem is encountered in many different technical areas, such as astronomical imaging, remote sensing, microscopy, medical imaging, optics, super-resolution applications, and motion tracking applications among others (see, for example, [1-9]). Astronomical imaging is one of the primary applications of blind image deconvolution algorithms [1,2]. Ground based imaging systems are subject to blurring due to the rapidly changing index of refractions of the atmosphere. Extraterrestrial observations of the Earth and the planets are degraded by motion blur as a result of slow camera shutter speeds relative to the rapid spacecraft motion. Blind image deconvolution is used for improving the quality of the Poisson distributed film grain noise present in the blurred X-rays, mammograms, and digital angiographic images. In such applications, many times, the blurring is unavoidable because the medical imaging systems limit the intensity (e.g., low X-ray intensity) of the incident radiation in order to protect the patient’s health [10]. In optics, blind image deconvolution is used to restore the original image from the degradation introduced by a microscope or any other optical instrument [6,7]. The Hubble Space Telescope main mirror imperfections have provided an inordinate amount of images for the digital image processing community [1]. As a final example, in tracking applications the object being tracked might be blurred due to its speed or the motion of the camera. As a result, the track is lost with conventional tracking approaches and the application of blind restoration approaches can improve tracking results [9]. The standard formulation of the gray-scale image degradation model is given in matrix-vector form by y = Hx + n, where the N × 1 vectors x, y, and n represent, respectively, the original image, the available noisy and blurred image, and the observation noise, and H represents the blurring matrix created from the blur point spread function h. The images are assumed to be of size m × n = N , and they are lexicographically ordered into N × 1 vectors. Given y, the BID problem calls for finding estimates of x and H using prior knowledge on them. A number of methods have been proposed to address BID (a recent literature review can be found in [11]). The most recent methods are based on a Bayesian framework, and have addressed the removal of the camera motion from blurred photographs [12-19]. In [12] the unknown image and blur were estimated in a two step process. In the first step the blur is estimated from the blurred photograph by regularizing the image gradients with a mixture of Gaussian distributions and by regularizing the blur with a mixture of exponential distributions. In the second step the image is estimated from the blurred photograph and the estimated blur by utilizing the Richardson-Lucy (RL) algorithm. Finally, the restored image is obtained by performing the histogram equalization (to that of the observed image) on the output of RL. A multi-scale approach is utilized for the algorithm implementation. The multi-scale approach consists of down-sampling the blurred photograph to number of low resolution images, and by utilizing the proposed algorithm iteratively to obtain the blur estimate at each resolution. Estimating camera motion has also been investigated in [13,14], where the unknown image and blur were estimated in a simultaneous fashion. Additionally, [14] concentrated on synthetic experiments where the performance of the algorithm was evaluated by the improvement in signal to noise ratio. Also, in [13,14] the regularization parameters are not automatically estimated but rather heuristically chosen by the user at each iteration in order to yield an unknown image estimate with good visual quality. The major disadvantage of the methods proposed in [13,14] compared to the method proposed in [12] is the lack of parameter estimation. In this article, we extend our study in [20] by providing (1) a multi-scale based implementation of the algorithm which improves the quality of the obtained restored images, and (2) a complete comparison with many existing state of the art blind deconvolution methods. The proposed Bayesian algorithm for BID utilizes a variant of the non-convex lp quasi norm based prior as the unknown image prior and the TV prior as the unknown blur prior. Furthermore, we utilize the Bayesian framework to provide the estimates for all model parameters. Finally, we evaluate the performance of the proposed algorithm and provide comparisons with [12-14,16,19] by restoring blurred photographs taken by commercial cameras. This article is organized as follows. In Section 2 we provide the proposed Bayesian modeling of the BID problem. The Bayesian inference is presented in Section 3. In Section 4, we describe implementation details of the proposed algorithm. Experimental results are provided in Section 5 and conclusions are drawn in Section 6. 2 Bayesian modeling As already discussed in the previous section, the observation noise is modeled as a zero mean white Gaussian ranp(y|β, x, h) ∝ βN/2 exp − y − Hx 2 , where β is the precision of the multivariate Gaussian distribution. As the image prior we utilize a variant of the generalized Gaussian distribution, given by d∈D where ZGG(α) is the partition function, 0 < p < 1, α denotes the set {αd} and d ∈ D = {h, v, hh, vv, hv}. ih(x) and iv(x) correspond to, respectively, the horizontal and vertical first order differences, at pixel i, that is, ih(x) = xi − xl(i) and iv(x) = xi − xa(i), where l(i) and a(i) denote the nearest neighbors of i, to the left and above, respectively. The operators ihh(x), ivv(x), ihv(x) correspond to, respectively, horizontal, vertical and horizontal-vertical second order differences, at pixel i. In this study, similarly to [15,21], we utilize a nonconvex lp quasi norm with 0 < p < 1 since the derivatives of blurry photographs are expected to be sparse. The distributions of the image derivatives often have heavier tails that are better modeled with the non-convex lp quasi norm prior with 0 < p < 1 compared to the convex priors modeled with p = 1, 2. For reducing the complexity of the problem we assume that αh = αv = α and αhh = αvv = αhv = α/2. Additionally, similarly to [22], the partition function is approximated as ZGG(α) ∝ α−λ1N/p, where λ1 is a positive real number. We then simplify (3) accordingly to obtain the following image prior p(x|α) ∝ αλ1N/p exp −α 21−o(d) d∈D i=1 where o(d) ∈ {1, 2} denotes the order of the difference operator id(x). For the blur we utilize the total-variation prior given by (see [23] for more details) p(h|γ ) ∝ γ λ2N exp [−γ TV(h)] , TV(h) = In this study, we use flat improper hyperpriors on α, β and γ , that is, we utilize Note that with this choice of the hyperpriors, the observed image y is made solely responsible for the estimation of the image, blur and hyperparameters. 3 Bayesian inference Bayesian inference on the unknown components of the blind image deconvolution problem is based on the estimation of the unknown posterior distribution p(α, β, γ , x, h | y), given by Assuming that x and h are independent, the joint distribution p(α, β, γ , x, h, y) can be factorized in terms of the observation model p(y|β, x, h), the prior distributions p(x|α) and p(h|γ ), and the hyperparameter distributions p(α), p(β) and p(γ ), that is, In this study, we adopt the maximum a posteriori (MAP) approach to obtain a single point estimate, ¯ = (α¯ , β¯, γ¯ , x¯, h¯), that maximizes p(α, β, γ , x, h | y) as follows, ¯ = argmax p(α, β, γ , x, h | y) = min 21−o(d) d∈D i λ1N N p log α − 2 log β − λ2N log γ . As can be seen from (10), obtaining the point estimate that maximizes the posterior distribution p(α, β, γ , x, h | y) is not straightforward since it requires the minimization of a non-convex functional. Maximizing the posterior distribution p(α, β, γ , x, h | y) by iteratively optimizing in one variable while fixing the others (the so called iterated conditional modes (ICM) method [24]) is equivalent to the variational Bayesian based maximization (see [25] for an example derivation) for the special case when all the posterior distributions are assumed to be degenerate. In this article, we apply the majorization-minimization approach twice to bound the non-convex functional to be minimized. We start by bounding the non-convex image prior p(x|α) by the functional M1(α, x, Z), that is The majorization-minimization approach has been utilized in several approaches for image restoration [25,26]. The functional M1(α, x, Z) is derived by considering the relationship between the weighted geometric and arithmetic means, which is given by tp/2z1−p/2 ≤ p2 t + 1 − p2 z, (12) p(x|α) = const · αλ1N/p exp −α 21−o(d) d∈D ≥ const · αλ1N/p exp − α2p d∈D 21−o(d) zd1−,i p/2 Then (11) holds by setting M1(α, x, Z) = αλ1N/p exp − α2p d∈D 21−o(d) zd1−,i p/2 where Z is a matrix with elements zd,i, with d ∈ {h, v, hh, vv, hv} and i = 1, . . . , N . Similarly, the majorization-minimization criterion is used to bound the blur prior p(h|γ ) utilizing the functional M2(γ , h, u). Let us define, for γ and any N dimensional vector u ∈ (R+)N , with components ui, i = 1, . . . , N , the following functional Using the inequality in (13) with p = 1, for t ≥ 0 and z > 0, that is, p t + 2−pp z Using (13) we obtain zd1−,i p/2 Therefore, we have t ≤ The lower bounds of p(x|α) and p(h|γ ) defined above lead to the following lower bound of the distribution p(α, β, γ , x, h, y) Table 1 ISNRxˆ and ISNRhˆ values, for the cameraman, satellite, shepp-logan, and airplane images degraded by five different motion blurs (BSNR = 40 dB) p(α, β, γ , x, h, y) = p(α)p(β)p(γ )p(x|α)p(h|γ )p(y|β, x, h) ≥ const · p(α)p(β)p(γ )M1(α, x, V)M2 Therefore, a single point estimate that maximizes the lower bound of the posterior distribution p(α, β, γ , x, h | y) is found as follows ¯ = min y − Hx 2 + 21−o(d) [ id(x)]2 + 2−pp zd,i zd1−,i p/2 − λ1pN log α − N2 log β − λ2N log γ . As shown in (20), we are effectively replacing the original non-convex minimization problem (10) by a series of convex ones by utilizing the majorization-minimization criteria and introducing the additional variational vectors zd and u. By iteratively solving this convex optimization problem in an alternating fashion with respect to all unknowns, we obtain a sequence of point estimates and derive the proposed algorithm as shown next. 3.1 Algorithm Given α1,β1,γ 1, h1, ui1 =[ ih(h1)]2 +[ iv(h1)]2, and zd1,i. for k = 1, 2, . . . until a stopping criterion is met: 1. Calculate 21−o(d)( d)tWkd( d) −1 Figure 3 Example convergence curves for the finest scale of the proposed algorithm of the airplane-color restoration for five different synthetic non-parametric motion blurs: (a) Blur 1. (b) Blur 2. (c) Blur 3. (d) Blur 4 and (e) Blur 5. Figure 4 Four different original images from [27]: (a) Image 1. (b) Image 2. (c) Image 3. (d) Image 4. Wkd(i, i) = (zdk,i)p/2−1. 2. Calculate where Wkd is a diagonal matrix with entries hk+1 = βk(Xk)t(Xk) + γ k d∈{h,v} 4. Calculate 5. Calculate where Uk is a diagonal matrix with entries Uk(i, i) = (uik)−1/2. 3. For each d ∈ {h, v, hh, vv, hv} calculate zdk+,i 1 =[ id(xk)]2 , uik+1 =[ ih(hk+1)]2 +[ iv(hk+1)]2 , d∈D 21−o(d) y − Hk+1xk 2 In the line of study presented in [21] the parameter p is set to 4/5 (see [21] for a detailed discussion). Additionally, the parameters λ1 and λ2 are needed to approximate the partition functions for prior distributions p(x|α) and p(h|γ ), respectively. Unfortunately, the approximations of partition functions for the distributions p(x|α) and p(h|γ ) are necessary since its corresponding partition functions are analytically intractable. We follow the approaches proposed in [22,23], as already described in the previous section, and determine the values of the parameters λ1 and λ2 experimentally. The values of the parameters p, λ1, and λ2 are therefore set throughout all the experiments that follow. The robustness of the proposed method will be tested and evaluated under various blurring and noisy conditions. Note that if the blur h and the hyperparameters α, β, and γ are assumed to be known, the proposed algorithm coincides with the iteratively re-weighted least squares (IRLS) algorithm presented in [21] (i.e., in this case for both algorithms the image estimate is calculated as shown in (21)). Note, that the lp quasi norm based prior is also utilized in [15], and that this study simplifies the prior used in [21] by omitting the second order derivatives. 4 Multi-scale implementation The restoration results presented in [12], and more recently in [27], showed the effectiveness of the multiscale approach in implementing blind image deconvolution algorithms. Furthermore, it is shown in [27] that the multi-scale approach prevents the algorithm from converging to the unit impulse. Alternatively, the authors in [13,14] introduced heuristic re-weighting of the regularization parameters, at each iteration, to prevent the Table 2 SSExˆ and ERxˆ values, for the images and blurs defined in Figures 4 and 5, respectively Blur Method Image 1 Image 2 Image 3 SSExˆ Blur 1 ALG 29.91 Fergus et al. 39.73 Cho et al. 33.05 Shan et al. 49.79 Levin et al. 44.06 Blur 2 ALG 33.47 Fergus et al. 40.70 Cho et al. 34.13 Shan et al. 38.87 Levin et al. 48.50 Blur 3 ALG 27.45 Fergus et al. 30.05 Cho et al. 31.41 Shan et al. 28.12 Levin et al. 34.93 Blur 4 ALG 51.59 Fergus et al. 125.80 Cho et al. 63.73 Shan et al. 100.43 Levin et al. 95.81 Blur 5 ALG 31.39 Fergus et al. 27.32 Cho et al. 38.59 Shan et al. 30.76 Levin et al. 26.50 Blur 6 ALG 20.72 Fergus et al. 44.02 Cho et al. 42.68 Shan et al. 71.33 Levin et al. 28.47 Blur 7 ALG 38.51 Fergus et al. 206.70 Cho et al. 43.46 Shan et al. 252.56 Levin et al. 45.91 Blur 8 ALG 30.07 Fergus et al. 49.42 Cho et al. 45.48 Shan et al. 158.72 Levin et al. 48.19 Entries of the table in bold denote the lowest values. Figure 6 Comparing the proposed algorithm with the methods proposed in [12,13,16,19] based on the restoration results from Table 2: Percentage of cases for which the restored image yields the smallest SSExˆ values. algorithms from converging to unrealistic blur estimates. In our study, no heuristic adjustment of each parameter is performed but instead we estimate the parameters automatically within the Bayesian framework. To avoid unrealistic blur estimates we adopt here a multi-scale scheme similar in spirit to the one proposed in [12]. Additionally, the proposed multi-scale approach allows user to automatically initialize the proposed algorithm without visually inspecting the observed blurred image for determining the initial blur estimate. By analyzing the observed blurred image it is possible to come up with more informative initial blur estimates; however in this study our focus is to develop a completely automated algorithm once the blur support is provided. The basic idea behind the multi-scale approach is to down-sample the observed blurred image to a number of low resolution images. At the lowest resolution the initial blur estimate (i.e., h1) is set to the uniform blur and the lowest resolution of the down-sampled observed image is utilized as the initial image estimate (i.e., x1). After convergence is achieved at each scale we up-sample the image and the blur estimates to the next higher resolution and re-run the proposed algorithm. This iterative process is repeated until the algorithm converges and image and blur estimates are obtained at their native resolutions. The detailed pseudocode used in the implementation of our multi-scale approach is shown in Appendix. 5 Experimental results In this section, we present the experimental results obtained by the use of the proposed algorithm. As the performance metric, for the experiments in which the original image is known, we utilize the improvement in signal to noise ratio of the restored image (denoted as ISNRxˆ ), which is defined as 10 log10 x − y 2/ x − xˆ 2 , Figure 7 Comparing the cumulative histogram of the proposed algorithm with the cumulative histogram of the methods proposed in [12,13,16,19] based on the restoration results from Table 2 in terms of the ERxˆ comparison metric. as 10 log10 h − hδ 2/ h − hˆ 2 , where h, hδ and hˆ are the original blur, the unit impulse, and the estimated blur, respectively. In addition, after the blur support, greater than the original one, is specified by the user, all unknown parameters and the estimates of the unknown image and blur Figure 8 Example restorations from Table 2: 1st column represents eight different blurred observations, 2nd column represents their respective original undistorted versions, 3rd column represents their respective restorations obtained by the proposed algorithm, 4th column represents their respective restorations obtained by the method proposed in [12], 5th column represents their respective original blurs, 6th column represents their respective restored blurs obtained by the proposed algorithm, 7th column represents their respective restored blurs obtained by the method proposed in [12]. are estimated automatically as described in Sections 3 and 4. Specifying the blur support for the unknown blur is common with the state-of-the-art approaches (see [12-14,16,19]). Similarly to these approaches, the proposed algorithm is very robust when the support of the blur is largely overestimated by the user, as can be seen in all the experiments that follow. Also, for the experiments in which Figure 9 Example restorations from Table 2: 1st column represents eight different blurred observations, 2nd column represents their respective original undistorted versions, 3rd column represents their respective restorations obtained by the proposed algorithm, 4th column represents their respective restorations obtained by the method proposed in [16], 5th column represents their respective original blurs, 6th column represents their respective restored blurs obtained by the proposed algorithm, 7th column represents their respective restored blurs obtained by the method proposed in [16]. blurred colored images are considered, only the luminance component of the observed image is restored while the observed chroma components are used to obtain the restored colored image once the original luminance is estimated. Finally, the proposed algorithm is terminated when the criterion xk − xk−1 / xk−1 < 10−3 is achieved or the number of iterations reaches 100. After each iteration, we enforce the following constraints on the blur estimates: the positivity (blur elements less than zero are set to zero), the support constraint (blur elements Figure 10 Example restorations from Table 2: 1st column represents eight different blurred observations, 2nd column represents their respective original undistorted versions, 3rd column represents their respective restorations obtained by the proposed algorithm, 4th column represents their respective restorations obtained by the method proposed in [13], 5th column represents their respective original blurs, 6th column represents their respective restored blurs obtained by the proposed algorithm, 7th column represents their respective restored blurs obtained by the method proposed in [13]. outside of the blur support estimate are set to zero), and the energy conservation (sum of the blur elements equals one). In the first set of experiments, we evaluate the performance of the proposed method on four standard images (cameraman, satellite, shepp-logan phantom and airplane) which are widely used in image restoration experiments. The original images are then blurred with five different motion blurs, which are shown in Figure 1. Realizations of white Gaussian noise are added to the respective blurred Figure 11 Example restorations from Table 2: 1st column represents eight different blurred observations, 2nd column represents their respective original undistorted versions, 3rd column represents their respective restorations obtained by the proposed algorithm, 4th column represents their respective restorations obtained by the method proposed in [19], 5th column represents their respective original blurs, 6th column represents their respective restored blurs obtained by the proposed algorithm, 7th column represents their respective restored blurs obtained by the method proposed in [19]. Figure 12 Comparing the proposed method with the method proposed in [14]: 1st column represents three different blurred observations, 2nd column represents their respective restorations obtained by the proposed algorithm, 3rd column represents their respective restorations obtained by the method proposed in [14], 4th column represents their respective blurs obtained by the proposed algorithm, 5th column represents their respective blurs obtained by the method proposed in [14]. images in order to obtain degraded images with the blurred signal to noise (BSNR) ratio of 40 dB. The blurred signal to noise ratio is defined as follows Var(Hx) BSNR = 10 log10 Var(n) , where Var(·) denotes the variance of the random sequence. Example restorations obtained by the proposed algorithm are shown in Figure 2. The restoration results in terms of the previously defined ISNRxˆ and ISNRhˆ , metrics are shown in Table 1. It can be observed from Table 1 that the proposed algorithm is very robust and it is capable of restoring the blurred images very successfully under various non-parametric motion blurs. Example convergence curves obtained by the proposed algorithm are shown in Figure 3. In the second set of experiments, we evaluate the performance of the proposed method on a set of 32 blurred test images taken by a commercial camera. The test images are obtained from [27] and they are available online (www.wisdom.weizmann.ac.il/∼levina/ papers/LevinEtalCVPR09Data.zip). The set of 32 blurred images was obtained by taking the original images shown in Figure 4 and by putting them side by side in order to form a calibration image. Once the calibration image was formed, a commercial camera was mounted on a tripod and eight photos of the calibration image were obtained. Note that during the acquisition process, the Z-axis rotation handle was locked in while the X-axis and Y-axis handles were loosened up in order to simulate in-plane camera shake (see [27] for details; resulting blurs are shown in Figure 5). In this study, we consider a comparison with the following methods [12,13,16,19]. For convenience, from now on, the methods proposed in [12,13,16], and the best method from [19] will be denoted, respectively, as Fergus et al., Shan et al., Cho et al., and Levin at al., while the proposed algorithm will be denoted as ALG. The restoration results, for the second set of experiments where the original image and blur are both known, in terms of Sum of Squared Errors (i.e., SSExˆ = x − xˆ 2) and the SSE ratio test (i.e., ERxˆ = SSExˆ /SSEx) defined ˜ in [27] are shown in Table 2. Note that x˜ is an image estimation obtained from the non-blind method from [27]. It can be observed from Table 2 that the proposed algorithm is very robust and that it is capable of restoring blurred images taken by a commercial camera very successfully under various non-parametric motion blurs. In addition, the proposed algorithm is very competitive with the state-of-the-art methods. As can be noted in Table 2, the SSExˆ metric for the proposed algorithm is respectively, on the average, 427, 89, 6, and 21 smaller than the SSExˆ metric for the Fergus at al., Shan at al., Cho at al., and Levin et al. methods. Note that in Fergus’ method, the blur estimation is performed separately from the image estimation. In order to understand the differences in the restoration results we provide some additional information. In Figure 6 it can be seen that for 66 % of tested cases the proposed algorithm yields the smallest SSExˆ values. In addition, Figure 7 shows that for a number of test cases the proposed algorithm is capable of achieving very small ERxˆ values. For example, there are 78 % of the cases for which the proposed algorithm has ERxˆ smaller than 1.5 while at the same time (as the second best) there are 53 % of the test cases for which Cho et al. method achieves such condition. Example restorations and comparison with Fergus at al., Cho at al., Shan at al, and Levin at al. are shown, respectively, in Figures 8, 9, 10, and 11. In the third set of experiments, we compare the performance of the proposed method with the method proposed in [14] by using the same set of blurred photographs as presented in [14]. Note that in the method proposed in [14] the parameters are not estimated but rather they are manually tuned which results in the sequence of numbers for each parameter. As can be seen in Figure 12 the restoration results obtained by the proposed algorithm are very competitive with the method proposed in [14]. It is clear from Figure 12 that the proposed algorithm produces much sharper restoration results with higher visual quality. Since we lack the true knowledge of the scene, the comparison metrics SSExˆ and ERxˆ are undefined for this experiment. 6 Conclusions In this article, a novel blind image deconvolution algorithm is presented. The proposed algorithm was developed within a Bayesian framework utilizing a non-convex lp quasi norm based sparse prior on the image, and a total-variation prior on the unknown blur. The proposed algorithm is completely automated once the blur support is provided. Experimental results demonstrate that using sparse priors and the proposed parameter estimation, both the unknown image and blur can be estimated with very high accuracy. Furthermore, numerous restorations of photographs taken by commercial cameras are provided to demonstrate the robustness and effectiveness of the proposed approach. Finally, it was shown that the performance of the proposed algorithm is competitive to existing state-of-the-art blind image deconvolution algorithms. Competing interests The authors declare that they have no competing interests. Acknowledgements This work was supported in part by the Department of Energy under contract DE-NA0000457 and the “Ministerio de Ciencia e Innovacio´ n” under contract TIN2010-15137. 1. J Krist, in Astronomical Data Analysis Software and Systems IV, ed . by byRA Shaw, HE Payne, and JJE Hayes . Simulation of HST PSFs Using Tiny Tim (Astronomical Society of the Pacific , San Francisco, USA, 1995 ), pp. 349 - 353 2. TJ Schultz, Multiframe blind deconvolution of astronomical images . J. Opt. Soc. Am. A . 10 , 1064 - 1073 ( 1993 ) 3. T Bretschneider, P Bones, S McNeill, D Pairman, in Proceedings of the American Society for Photogrammetry & Remote Sensing. Image-based quality assessment of SPOT data , ( 2001 ). [ Unpaginated CD -ROM] 4. FS Gibson , F Lanni, Experimental test of an analytical model of aberration in an oil-immersion objective lens used in three-dimensional light microscopy . J. Opt. Soc. Am. A . 8 , 1601 - 1613 ( 1991 ) 5. O Michailovich, D Adam, A novel approach to the 2-D blind deconvolution problem in medical ultrasound . IEEE Trans. Med. Imaging . 24 , 86 - 104 ( 2005 ) 6. M Roggemann, Limited degree-of-freedom adaptive optics and image reconstruction . Appl. Opt . 30 , 4227 - 4233 ( 1991 ) 7. P Nisenson, R Barakat, Partial atmospheric correction with adaptive optics . J. Opt. Soc. Am. A . 4 , 2249 - 2253 ( 1991 ) 8. AK RMCA Segall, Katsaggelos, High-resolution images from low-resolution compressed video . IEEE Signal Process. Mag . 20 ( 3 ), 37 - 48 ( 2003 ) 9. S Dai, M Yang, Y Wu, AK Katsaggelos , in Proc. IEEE Int Image Processing Conf . Tracking motion-blurred targets in video , ( 2006 ), pp. 2389 - 2392 10. CJKK Faulkner, M Louka, in Third Int. Conf. on Image Proc. and Its Applications . Veiling glare deconvolution of images produced by x-ray image intensifiers , ( 1989 ), pp. 669 - 673 11. TE Bishop, SD Babacan, B Amizic, AK Katsaggelos, T Chan, R Molina, Blind Image Deconvolution: Problem Formulation and Existing Approaches. (CRC Press , 2007 ) 12. R Fergus, B Singh, A Hertzmann, ST Roweis, WT Freeman, Removing camera shake from a single photograph . ACM Trans Graph . 25 ( 3 ), 787 - 794 ( 2006 ) 13. Q Shan, J Jia, A Agarwala, in SIGGRAPH '08: ACM SIGGRAPH 2008 papers. High-quality motion deblurring from a single image (ACM , New York, 2008 ), pp. 1 - 10 14. M Almeida, L Almeida, Blind and semi-blind deblurring of natural images . IEEE Trans. Image Process . 19 , 36 - 52 ( 2010 ) 15. D Krishnan, R Fergus, in Advances in Neural Information Processing Systems 22 , ed. by Y Bengio, D Schuurmans, J Lafferty, CKI Williams, and A Culotta. Fast image deconvolution using hyper-Laplacian priors , ( 2009 ), pp. 1033 - 1041 16. S Cho, S Lee, Fast motion deblurring . ACM Trans Graph. (SIGGRAPH ASIA 2009 ). 28 ( 5 ). (Article No. 145, 2009 ) 17. TS Cho, N Joshi, CL Zitnick, SB Kang, R Szeliski, WT Freeman, in IEEE Conference on Computer Vision and Pattern Recognition (CVPR). A content-aware image prior , ( 2010 ), pp. 169 - 176 18. T Hou, S Wang, H Qin, Image deconvolution with multi-stage convex relaxation and its perceptual evaluation . IEEE Trans. Image Process . 20 (12), 3383 - 3392 ( 2011 ) 19. A Levin, Y Weiss, F Durand, W Freeman, in IEEE Conference on Computer Vision and Pattern Recognition (CVPR) . Efficient marginal likelihood optimization in blind deconvolution , ( 2011 ), pp. 2657 - 2664 20. B Amizic, SD Babacan, R Molina, AK Katsaggelos, in European Signal Processing Conference Eusipco. Sparse Bayesian Blind Image Deconvolution with Parameter Estimation (Aalborg, Denmark , 2010 ), pp. 626 - 630 21. A Levin, R Fergus, F Durand, WT Freeman, in SIGGRAPH '07: ACM SIGGRAPH 2007 papers. Image and depth from a conventional camera with a coded aperture (ACM , New York, 2007 ), p. 70 22. A Mohammad-Djafari , in Maximum Entropy and Bayesian Methods . A full bayesian approach for inverse problems (Kluwer Academic Publishers, 1996 ), pp. 135 - 143 23. J Bioucas-Dias , M Figueiredo , J Oliveira, in Proceedings of EUSIPCO' 2006 . Adaptive total-variation image deconvolution: a majorization-minimization approach (Florence , Italy, p. 2006 24. J Besag, On the statistical analysis of dirty pictures . J. Royal Stat. Soc. Series B (Methodological) . 48 ( 3 ), 259 - 302 ( 1986 ) 25. S Babacan, R Molina, A Katsaggelos, Parameter estimation in TV image restoration using variational distribution approximation . IEEE Trans. Image Process . 17 ( 3 ), 326 - 339 ( 2008 ) 26. J Bioucas-Dias , M Figueiredo , J Oliveira, in ICASSP, vol. 2. Total-variation image deconvolution: a majorization-minimization approach , ( 2006 ), p. II 27. A Levin, Y Weiss, F Durand, WT Freeman, in Proc. IEEE Conf. Computer Vision and Pattern Recognition CVPR. Understanding and evaluating blind deconvolution algorithms , ( 2009 ), pp. 1964 - 1971


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1186%2F1687-5281-2012-20.pdf

Bruno Amizic, Rafael Molina. Sparse Bayesian blind image deconvolution with parameter estimation, EURASIP Journal on Image and Video Processing, 2012, 20, DOI: 10.1186/1687-5281-2012-20