Noise models for low counting rate coherent diffraction imaging

被引:86
作者
Godard, Pierre [1 ,2 ]
Allain, Marc [1 ]
Chamard, Virginie [1 ]
Rodenburg, John [2 ]
机构
[1] Univ Aix Marseille, CNRS, Inst Fresnel, Fac St Jerome, F-13397 Marseille 20, France
[2] Univ Sheffield, Dept Elect & Elect Engn, Sheffield S1 3JD, S Yorkshire, England
基金
英国工程与自然科学研究理事会;
关键词
PHASE-RETRIEVAL; TRANSMISSION MICROSCOPY; ORDERED SUBSETS; RAY; RECONSTRUCTION; TOMOGRAPHY; ALGORITHM; EMISSION;
D O I
10.1364/OE.20.025914
中图分类号
O43 [光学];
学科分类号
070207 ; 0803 ;
摘要
Coherent diffraction imaging (CDI) is a lens-less microscopy method that extracts the complex-valued exit field from intensity measurements alone. It is of particular importance for microscopy imaging with diffraction set-ups where high quality lenses are not available. The inversion scheme allowing the phase retrieval is based on the use of an iterative algorithm. In this work, we address the question of the choice of the iterative process in the case of data corrupted by photon or electron shot noise. Several noise models are presented and further used within two inversion strategies, the ordered subset and the scaled gradient. Based on analytical and numerical analysis together with Monte-Carlo studies, we show that any physical interpretations drawn from a CDI iterative technique require a detailed understanding of the relationship between the noise model and the used inversion method. We observe that iterative algorithms often assume implicitly a noise model. For low counting rates, each noise model behaves differently. Moreover, the used optimization strategy introduces its own artefacts. Based on this analysis, we develop a hybrid strategy which works efficiently in the absence of an informed initial guess. Our work emphasises issues which should be considered carefully when inverting experimental data. (C) 2012 Optical Society of America
引用
收藏
页码:25914 / 25934
页数:21
相关论文
共 46 条
[1]   Globally convergent image reconstruction for emission tomography using relaxed ordered subsets algorithms [J].
Ahn, S ;
Fessler, JA .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2003, 22 (05) :613-626
[2]   Imaging techniques for gamma-ray diffuse emission: application to INTEGRAL/SPI [J].
Allain, M ;
Roques, JP .
ASTRONOMY & ASTROPHYSICS, 2006, 447 (03) :1175-1187
[3]  
[Anonymous], 1999, Athena scientific Belmont
[4]   THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA [J].
ANSCOMBE, FJ .
BIOMETRIKA, 1948, 35 (3-4) :246-254
[5]   Chemical Contrast in Soft X-Ray Ptychography [J].
Beckers, Mike ;
Senkbeil, Tobias ;
Gorniak, Thomas ;
Reese, Michael ;
Giewekemeyer, Klaus ;
Gleber, Sophie-Charlotte ;
Salditt, Tim ;
Rosenhahn, Axel .
PHYSICAL REVIEW LETTERS, 2011, 107 (20)
[6]  
Bertero M., 1998, INTRO INVERSE PROBLE
[7]  
BOUCHET L, 1995, ASTRON ASTROPHYS SUP, V113, P167
[8]   A unified approach to statistical tomography using coordinate descent optimization [J].
Bouman, CA ;
Sauer, K .
IEEE TRANSACTIONS ON IMAGE PROCESSING, 1996, 5 (03) :480-492
[9]   The PILATUS 1M detector [J].
Broennimann, C ;
Eikenberry, EF ;
Henrich, B ;
Horisberger, R ;
Huelsen, G ;
Pohl, E ;
Schmitt, B ;
Schulze-Briese, C ;
Suzuki, M ;
Tomizaki, T ;
Toyokawa, H ;
Wagner, A .
JOURNAL OF SYNCHROTRON RADIATION, 2006, 13 :120-130
[10]   BICAV: A block-iterative parallel algorithm for sparse systems with pixel-related weighting [J].
Censor, Y ;
Gordon, D ;
Gordon, R .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2001, 20 (10) :1050-1060