光学 精密工程, 2013, 21 (11): 2943, 网络出版: 2014-01-09
结合参数估计的天文图像极大似然恢复
Astronomical image restoration based on maximum-likelihood incorporated parameter estimation
图像恢复 天文图像 极大似然原理 点扩散函数估计 混合高斯泊松噪声 image restoration astronomical image maximum-likelihood Point Spread Function(PSF) estimation mixed Gaussian-Poisson noise
摘要
分析了Benvenuto等针对天文图像恢复提出的基于极大似然(ML)代价函数的有效逼近模型, 由此提出了一种比传统ML收敛更快的图像恢复算法。该算法在未知点扩散函数(PSF)条件下, 通过观测模糊图像, 自适应估计湍流PSF, 使PSF估计更符合成像环境; 然后, 将该算法与混合高斯泊松噪声的ML算法相结合, 形成增强ML迭代算法。在迭代过程中动态更新PSF, 交替执行恢复图像、去除噪声等策略。结果显示:对于点源目标图像, 本文算法恢复图像的质量在峰值信噪比、均方误差以及平均对比度3个指标上分别提高了96.64%, 69.26%和25.6%; 对于真实湍流退化图像, 恢复质量也有一定改善。结论表明:该方法可以使迭代过程收敛更稳定, 图像恢复质量得到明显提高, 非常适用于天文观测图像的高清晰恢复与重建。
Abstract
The effective approximation mode based on Maximum-likelihood (ML) function proposed by Benvenuto was analyzed for astronomy image restoration, then a new image restoration algorithm with convergence faster than that of traditional ML method was proposed. In this algorithm, PSF known a priori was not required. The turbulence PSF was estimated from observed blur images to make the PSF estimation more accordance with an imaging environment. By incorporating adaptive estimation of PSF into ML restoration, an enhanced ML algorithm was presented. Additionally, the PSF was updated successively during iteration, and the ML restoration and denoising were performed alternatively in iteration. The results show that the proposed algorithm works much better than ML does. Taking the point source image for an instance, proposed method improves the image quality by 96.64%, 69.26% and 25.6% respectively on the peak signal to noise ratio, mean square error and the correlation coefficient. In conclusion, the algorithm allow the iterative process in ML algorithm to converge stably and the image quality to be improved. Experiment results show that the presented method can be used routinely in astronomical image restoration.
耿则勋, 魏小峰, 沈忱. 结合参数估计的天文图像极大似然恢复[J]. 光学 精密工程, 2013, 21(11): 2943. GENG Ze-xun, WEI Xiao-feng, SHEN Chen. Astronomical image restoration based on maximum-likelihood incorporated parameter estimation[J]. Optics and Precision Engineering, 2013, 21(11): 2943.