Reconstruction algorithm for single photon emission computed tomography and its numerical implementation

被引:28
作者
Fokas, AS [1 ]
Iserles, A [1 ]
Marinakis, V [1 ]
机构
[1] Univ Cambridge, Dept Appl Math & Theoret Phys, Cambridge, England
关键词
single photon emission computed tomography (SPECT); attenuated Radon transform; imaging techniques;
D O I
10.1098/rsif.2005.0061
中图分类号
O [数理科学和化学]; P [天文学、地球科学]; Q [生物科学]; N [自然科学总论];
学科分类号
07 ; 0710 ; 09 ;
摘要
The modern imaging techniques of positron emission tomography and of single photon emission computed tomography are not only two of the most important tools for studying the functional characteristics of the brain, but they now also play a vital role in several areas of clinical medicine, including neurology, oncology and cardiology. The basic mathematical problems associated with these techniques are the construction of the inverse of the Radon transform and of the inverse of the so-called attenuated Radon transform, respectively. An exact formula for the inverse Radon transform is well known, whereas that for the inverse attenuated Radon transform was obtained only recently by R. Novikov. The latter formula was constructed by using a method introduced earlier by R. Novikov and the first author in connection with a novel derivation of the inverse Radon transform. Here, we first show that the appropriate use of that earlier result yields immediately an analytic formula for the inverse attenuated Radon transform. We then present an algorithm for the numerical implementation of this analytic formula, based on approximating the given data in terms of cubic splines. Several numerical tests are presented which suggest that our algorithm is capable of producing accurate reconstruction for realistic phantoms such as the well-known Shepp-Logan phantom.
引用
收藏
页码:45 / 54
页数:10
相关论文
共 38 条
[21]   REDUCED CEREBRAL GLUCOSE-METABOLISM IN ASYMPTOMATIC SUBJECTS AT RISK FOR HUNTINGTONS-DISEASE [J].
MAZZIOTTA, JC ;
PHELPS, ME ;
PAHL, JJ ;
HUANG, SC ;
BAXTER, LR ;
RIEGE, WH ;
HOFFMAN, JM ;
KUHL, DE ;
LANTO, AB ;
WAPENSKI, JA ;
MARKHAM, CH .
NEW ENGLAND JOURNAL OF MEDICINE, 1987, 316 (07) :357-362
[22]  
Mazziotta JC, 1992, CLIN BRAIN IMAGING P, P244
[23]  
MINOSHIMA S, 1995, J NUCL MED, V36, P1238
[24]  
NATTERE F, 1986, MATH COMPUTERIZED TO
[25]   Inversion of the attenuated Radon transform [J].
Natterer, P .
INVERSE PROBLEMS, 2001, 17 (01) :113-119
[26]   An inversion formula for the attenuated X-ray transformation [J].
Novikov, RG .
ARKIV FOR MATEMATIK, 2002, 40 (01) :145-167
[27]   A penalized-likelihood image reconstruction method for emission tomography, compared to postsmoothed maximum-likelihood with matched spatial resolution [J].
Nuyts, J ;
Fessler, JA .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2003, 22 (09) :1042-1052
[28]   The solitary pulmonary nodule [J].
Ost, D ;
Fein, AM ;
Feinsilver, SH .
NEW ENGLAND JOURNAL OF MEDICINE, 2003, 348 (25) :2535-2542
[29]  
PARK MG, 1996, ARCH NEUROL-CHICAGO, V53, P665
[30]  
Press W.H., 1992, NUMERICAL RECIPES FO