Dynamic PET reconstruction using wavelet regularization with adapted basis functions.
Tomographic reconstruction from positron emission tomography (PET) data is an ill-posed problem that requires regularization. An attractive approach is to impose an l(1) -regularization constraint, which favors sparse solutions in the wavelet domain. This can be achieved quite efficiently thanks to the iterative algorithm developed by Daubechies et al., 2004. In this paper, we apply this technique and extend it for the reconstruction of dynamic (spatio-temporal) PET data. Moreover, instead of using classical wavelets in the temporal dimension, we introduce exponential-spline wavelets (E-spline wavelets) that are specially tailored to model time activity curves (TACs) in PET. We show that the exponential-spline wavelets naturally arise from the compartmental description of the dynamics of the tracer distribution. We address the issue of the selection of the "optimal" E-spline parameters (poles and zeros) and we investigate their effect on reconstruction quality. We demonstrate the usefulness of spatio-temporal regularization and the superior performance of E-spline wavelets over conventional Battle-LemariE wavelets in a series of experiments: the 1-D fitting of TACs, and the tomographic reconstruction of both simulated and clinical data. We find that the E-spline wavelets outperform the conventional wavelets in terms of the reconstructed signal-to-noise ratio (SNR) and the sparsity of the wavelet coefficients. Based on our simulations, we conclude that replacing the conventional wavelets with E-spline wavelets leads to equal reconstruction quality for a 40% reduction of detected coincidences, meaning an improved image quality for the same number of counts or equivalently a reduced exposure to the patient for the same image quality.
Verhaeghe
Jeroen
J
Department of Electronics and Information Systems, MEDISIP, Ghent University-IBBT-IBiTech, De Pintelaan 185 block B, B-9000 Ghent, Belgium. jeroen.verhaeghe@ugent.be
Van de Ville
Dimitri
D
Khalidov
Ildar
I
D'Asseler
Yves
Y
Lemahieu
Ignace
I
Unser
Michael
M
D000465:Algorithms; D001185:Artificial Intelligence; D016000:Cluster Analysis; D003198:Computer Simulation; D003627:Data Interpretation, Statistical; D005246:Feedback; D005583:Fourier Analysis; D006801:Humans; D021621:Imaging, Three-Dimensional / Q000379:methods*; D016247:Information Storage and Retrieval / Q000379:methods; D007668:Kidney / Q000000981:diagnostic imaging; D008099:Liver / Q000000981:diagnostic imaging; D010363:Pattern Recognition, Automated; D049268:Positron-Emission Tomography / Q000379:methods*; D012815:Signal Processing, Computer-Assisted; D013997:Time Factors