Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging.
A penalized-likelihood (PL) SPECT reconstruction method using a modified regularizer that accounts for anatomical boundary side information was implemented to achieve accurate estimates of both the total target activity and the activity distribution within targets. In both simulations and experimental I-131 phantom studies, reconstructions from (1) penalized likelihood employing CT-side information-based regularization (PL-CT), (2) penalized likelihood with edge preserving regularization (no CT) and (3) penalized likelihood with conventional spatially invariant quadratic regularization (no CT) were compared with (4) ordered subset expectation maximization (OSEM), which is the iterative algorithm conventionally used in clinics for quantitative SPECT. Evaluations included phantom studies with perfect and imperfect side information and studies with uniform and non-uniform activity distributions in the target. For targets with uniform activity, the PL-CT images and profiles were closest to the 'truth', avoided the edge offshoots evident with OSEM and minimized the blurring across boundaries evident with regularization without CT information. Apart from visual comparison, reconstruction accuracy was evaluated using the bias and standard deviation (STD) of the total target activity estimate and the root mean square error (RMSE) of the activity distribution within the target. PL-CT reconstruction reduced both bias and RMSE compared with regularization without side information. When compared with unregularized OSEM, PL-CT reduced RMSE and STD while bias was comparable. For targets with non-uniform activity, these improvements with PL-CT were observed only when the change in activity was matched by a change in the anatomical image and the corresponding inner boundary was also used to control the regularization. In summary, the present work demonstrates the potential of using CT side information to obtain improved estimates of the activity distribution in targets without sacrificing the accuracy of total target activity estimation. The method is best suited for data acquired on hybrid systems where SPECT-CT misregistration is minimized. To demonstrate clinical application, the PL reconstruction with CT-based regularization was applied to data from a patient who underwent SPECT/CT imaging for tumor dosimetry following I-131 radioimmunotherapy.
Dewaraja
Yuni K
YK
Department of Radiology, University of Michigan, Ann Arbor, MI 48109, USA. yuni@umich.edu
Koral
Kenneth F
KF
Fessler
Jeffrey A
JA
D000911:Antibodies, Monoclonal / Q000627:therapeutic use; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*; D008228:Lymphoma, Non-Hodgkin / Q000000981:diagnostic imaging / Q000532:radiotherapy; D019047:Phantoms, Imaging; D016499:Radioimmunotherapy; D015899:Tomography, Emission-Computed, Single-Photon / Q000379:methods*; D014057:Tomography, X-Ray Computed
D000911:Antibodies, Monoclonal; C119496:tositumomab I-131
IEEE Trans Med Imaging. 1987;6(3):228-38
18244025
IEEE Trans Med Imaging. 2005 Aug;24(8):969-76
16092329
IEEE Trans Image Process. 1996;5(9):1346-58
18285223
Phys Med Biol. 2007 Feb 7;52(3):777-90
17228120
IEEE Trans Med Imaging. 1993;12(4):670-80
18218461
IEEE Trans Image Process. 2007 Oct;16(10):2411-22
17926925
IEEE Trans Med Imaging. 1996;15(5):673-86
18215949
Phys Med Biol. 1999 Nov;44(11):2835-51
10588288
IEEE Trans Med Imaging. 2005 May;24(5):667-75
15889553
Cancer Biother Radiopharm. 2009 Aug;24(4):417-26
19694576
J Nucl Med. 2005 May;46(5):840-9
15872359
J Nucl Med. 2010 Apr;51(4):654-9
20237032
Phys Med Biol. 2002 Jan 7;47(1):1-20
11814220
Phys Med Biol. 2008 Sep 7;53(17):4595-604
18678930
Phys Med Biol. 2008 Jun 21;53(12):3099-112
18495976
Phys Med Biol. 1994 Mar;39(3):509-30
15551595
J Nucl Med. 2007 Jun;48(6):1008-16
17504874
IEEE Trans Med Imaging. 1994;13(2):290-300
18218505