Fan Yi, Lu Hongbing, Hao Chongyang, Liang Zhengrong, Zhou Zhiming
Department of Electronic Information, NorthWestern Polytechnical University, Xi'an, Shaanxi 710071, P. R. China, Email:
Int J Image Graph. 2007 Jan;7(1):87-104. doi: 10.1142/S0219467807002556.
Conventionally, the inverse problem of gated cardiac SPECT is solved by reconstructing the images frame-by-frame, ignoring the inter-frame correlation along the time dimension. To compensate for the non-uniform attenuation for quantitative cardiac imaging, iterative image reconstruction has been a choice which could utilize an a priori constraint on the inter-frame correlation for a penalized maximum likelihood (ML) solution. However, iterative image reconstruction in the 4D space involves intensive computations. In this paper, an efficient method for 4D gated SPECT reconstruction is developed based on Karhune-Loève (KL) transform and Novikov's inverse formula. The temporal KL transform is first applied on the data sequence to de-correlate the inter-frame correlation and then the 3D principal components in the KL domain are reconstructed frame-by-frame using Novikov's inverse formula with non-uniform attenuation compensation. Finally an inverse KL transform is performed to obtain quantitatively-reconstructed 4D images in the original space. With the proposed method, 4D reconstruction can be achieved at a reasonable computational cost. The results from computer simulations are very encouraging as compared to conventional frame-by-frame filtered back-projection and iterative ordered-subsets ML reconstructions. By discarding high-order KL components for further noise reduction, the computation time could be further reduced.
传统上,门控心脏单光子发射计算机断层扫描(SPECT)的逆问题是通过逐帧重建图像来解决的,忽略了时间维度上的帧间相关性。为了补偿定量心脏成像中的非均匀衰减,迭代图像重建是一种选择,它可以利用对帧间相关性的先验约束来进行惩罚最大似然(ML)求解。然而,四维空间中的迭代图像重建涉及大量计算。本文基于卡尔胡宁-勒夫(KL)变换和诺维科夫逆公式,开发了一种有效的四维门控SPECT重建方法。首先对数据序列应用时间KL变换以消除帧间相关性,然后在KL域中使用具有非均匀衰减补偿的诺维科夫逆公式逐帧重建三维主成分。最后进行逆KL变换,以在原始空间中获得定量重建的四维图像。使用所提出的方法,可以以合理的计算成本实现四维重建。与传统的逐帧滤波反投影和迭代有序子集ML重建相比,计算机模拟结果非常令人鼓舞。通过丢弃高阶KL成分以进一步降低噪声,可以进一步减少计算时间。