An image restoration technique with error estimates

被引:46
作者
Esch, DN
Connors, A
Karovska, M
van Dyk, DA
机构
[1] Dept Vet Affairs, Ctr Hlth Qual Outcomes & Econ Res, Bedford, MA 01730 USA
[2] Eureka Sci, Oakland, CA 94602 USA
[3] Harvard Smithsonian Ctr Astrophys, Cambridge, MA 02138 USA
[4] Univ Calif Irvine, Dept Stat, Irvine, CA 92697 USA
关键词
methods : data analysis; techniques : high angular resolution;
D O I
10.1086/421761
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
Image restoration including deconvolution techniques offers a powerful tool to improve resolution in images and to extract information on the multiscale structure stored in astronomical observations. We present a new method for statistical deconvolution, which we call expectation through Markov Chain Monte Carlo (EMC2). This method is designed to remedy several shortfalls of currently used deconvolution and restoration techniques for Poisson data. We use a wavelet-like multiscale representation of the true image to achieve smoothing at all scales of resolution simultaneously, thus capturing detailed features in the image at the same time as larger scale extended features. Thus, this method smooths the image, while maintaining the ability to effectively reconstruct point sources and sharp features in the image. We use a principled, fully Bayesian model-based analysis, which produces extensive information about the uncertainty in the fitted smooth image, allowing assessment of the errors in the resulting reconstruction. Our method also includes automatic fitting of the multiscale smoothing parameters. We show several examples of application of EMC2 to both simulated data and a real astronomical X-ray source.
引用
收藏
页码:1213 / 1227
页数:15
相关论文
共 32 条
[1]  
[Anonymous], MEMSYS5 QUANTIFIED M
[2]  
Carter C., 2003, ASP C SER, V295, P477
[3]  
COIFMAN R. R., 1995, Wavelets and statistics, P125, DOI DOI 10.1007/978-1-4612-2544-7_9
[4]   MAXIMUM LIKELIHOOD FROM INCOMPLETE DATA VIA EM ALGORITHM [J].
DEMPSTER, AP ;
LAIRD, NM ;
RUBIN, DB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-METHODOLOGICAL, 1977, 39 (01) :1-38
[5]  
ESCH DN, 2003, THESIS HARVARD U
[6]   RESTORING WITH MAXIMUM LIKELIHOOD AND MAXIMUM ENTROPY [J].
FRIEDEN, BR .
JOURNAL OF THE OPTICAL SOCIETY OF AMERICA, 1972, 62 (04) :511-&
[7]  
Gelman A, 1996, STAT SINICA, V6, P733
[8]  
Gelman A., 1992, STAT SCI, V7, P457, DOI DOI 10.1214/SS/1177011136
[9]  
HASTINGS WK, 1970, BIOMETRIKA, V57, P97, DOI 10.1093/biomet/57.1.97
[10]   Deconvolution of sparse spike trains by iterated window maximization [J].
Kaaresen, KF .
IEEE TRANSACTIONS ON SIGNAL PROCESSING, 1997, 45 (05) :1173-1183