Noise reduction and convergence of Bayesian algorithms with blobs based on the Huber function and median root prior

被引:39
作者
Chlewicki, W
Hermansen, F
Hansen, SB
机构
[1] Aarhus Univ Hosp, PET Ctr, DK-8000 Aarhus C, Denmark
[2] Univ Patras, Sch Med, Dept Med Phys, GR-26500 Patras, Greece
[3] Tech Univ Szczecin, Inst Elect Telecommun & Informat, PL-70313 Szczecin, Poland
关键词
D O I
10.1088/0031-9155/49/20/004
中图分类号
R318 [生物医学工程];
学科分类号
0831 ;
摘要
Iterative image reconstruction algorithms have the potential to produce low noise images. Early stopping of the iteration process is problematic because some features of the image may converge slowly. On the other hand, there may be noise build-up with increased number of iterations. Therefore, we examined the stabilizing effect of using two different prior functions as well as image representation by blobs so that the number of iterations could be increased without noise build-up. Reconstruction was performed of simulated phantoms and of real data acquired by positron emission tomography. Image quality measures were calculated for images reconstructed with or without priors. Both priors stabilized the iteration process. The first prior based on the Huber function reduced the noise without significant loss of contrast recovery of small spots, but the drawback of the method was the difficulty in finding optimal values of two free parameters. The second method based on a median root prior has only one Bayesian parameter which was easy to set, but it should be taken into account that the image resolution while using that prior has to be chosen sufficiently high not to cause the complete removal of small spots. In conclusion, the Huber penalty function gives accurate and low noise images, but it may be difficult to determine the parameters. The median root prior method is not quite as accurate but may be used if image resolution is increased.
引用
收藏
页码:4717 / 4730
页数:14
相关论文
共 27 条
[1]   Generalization of median root prior reconstruction [J].
Alenius, S ;
Ruotsalainen, U .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2002, 21 (11) :1413-1420
[2]   Bayesian image reconstruction for emission tomography based on median root prior [J].
Sakari Alenius ;
Ulla Ruotsalainen .
European Journal of Nuclear Medicine, 1997, 24 (3) :258-265
[3]   Tomographic image reconstruction based on a content-adaptive mesh model [J].
Brankov, JG ;
Yang, Y ;
Wernick, MN .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2004, 23 (02) :202-212
[4]   FINITE SERIES-EXPANSION RECONSTRUCTION METHODS [J].
CENSOR, Y .
PROCEEDINGS OF THE IEEE, 1983, 71 (03) :409-419
[5]  
Chlewicki W, 2001, MEDICON 2001: PROCEEDINGS OF THE INTERNATIONAL FEDERATION FOR MEDICAL & BIOLOGICAL ENGINEERING, PTS 1 AND 2, P495
[6]   Grouped coordinate descent algorithms for robust edge-preserving image restoration [J].
Fessler, JA .
IMAGE RECONSTRUCTION AND RESTORATION II, 1997, 3170 :184-194
[7]  
GARCIA I, 1997, DEV GLOBAL OPTIMIZAT, P289
[8]   BAYESIAN RECONSTRUCTIONS FROM EMISSION TOMOGRAPHY DATA USING A MODIFIED EM ALGORITHM [J].
GREEN, PJ .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1990, 9 (01) :84-93
[9]   LOCAL BASIS-FUNCTION APPROACH TO COMPUTED-TOMOGRAPHY [J].
HANSON, KM ;
WECKSUNG, GW .
APPLIED OPTICS, 1985, 24 (23) :4028-4039
[10]   A GENERALIZED EM ALGORITHM FOR 3-D BAYESIAN RECONSTRUCTION FROM POISSON DATA USING GIBBS PRIORS [J].
HEBERT, T ;
LEAHY, R .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1989, 8 (02) :194-202