Rahman Ashequr, Zhu Yansong, Clarkson Eric, Kupinski Matthew A, Frey Eric C, Jha Abhinav K
Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA.
Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA.
Inverse Probl. 2020 Aug;36(8). doi: 10.1088/1361-6420/ab958b. Epub 2020 Aug 20.
The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPU) for acceleration. The method was applied to analyze the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2-D digital synthetic and anthropomorphic phantoms. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.
在单光子发射计算机断层扫描(SPECT)成像中,无需单独进行透射扫描即可执行衰减和散射补偿(ASC)的潜力非常巨大。在这种情况下,SPECT中的衰减主要是由于康普顿散射,其中康普顿散射的概率与组织的衰减系数成正比,并且散射光子的能量与散射角相关。基于此前提,我们研究了以列表模式(LM)格式获取并包含能量信息的SPECT散射光子数据是否可用于估计衰减图。为此,我们提出了一种基于费舍尔信息的方法,该方法仅使用SPECT发射数据就能为联合估计活度和衰减分布的任务得出克拉美-罗界(CRB)。在此过程中,提出了一种基于路径的形式主义来处理包括散射光子数据在内的LM SPECT发射数据。费舍尔信息方法在NVIDIA图形处理单元(GPU)上实现以加速计算。该方法被应用于分析SPECT LM发射数据的信息含量,该数据包含高达一阶的散射事件,该分析是在一个模拟SPECT系统中进行的,该系统的参数通过使用二维数字合成和拟人化体模的实际计算研究来模拟临床系统。该方法还应用于包含高达二阶散射的合成体模的LM数据。使用拟人化体模进行的实验模拟了心肌灌注和多巴胺转运体(DaT)-扫描SPECT研究。结果表明,获得的衰减系数和活度系数的CRB通常远低于这些系数的真实值。检测到的光子数量增加会使衰减系数和活度系数的CRB降低。此外,我们观察到能量分辨率更好的系统会使衰减系数的CRB更低。总体而言,结果表明包括散射光子在内的LM SPECT发射数据包含联合估计活度系数和衰减系数的信息。