RESUMO
Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.
RESUMO
The risk stratification and long-term survival of patients with orthotopic heart transplantation (OHT) is impacted by the complication of cardiac allograft vasculopathy (CAV). This study evaluates changes in myocardial blood flow (MBF) and myocardial coronary flow reserve (CFR) in a group of long-term OHT patients using quantitative cardiac 82Rb-positron emission tomography (PET). Twenty patients (7 females and 13 males, mean age = 72.7 ± 12.2 years with CAV and 62.9 ± 7.2 years without CAV and post-OHT mean time = 13.9 years), were evaluated retrospectively using dynamic cardiac 82Rb-PET at rest and regadenoson-induced stress. The patients also underwent selective coronary angiography (SCA) for diagnosis and risk stratification. CAV was diagnosed based on SCA findings and maximal intimal thickness greater than 0.5 mm, as defined by International Society of Heart and Lung Transplantation (ISHLT). Global and regional MBFs were estimated in three vascular territories using the standard 1-tissue compartment model for dynamic 82Rb-PET. The myocardial CFR was also calculated as the ratio of peak stress MBF to rest MBF. Among twenty patients, seven had CAV in, at least, one major coronary artery (ISHLT CAV grade 1 or higher) while 13 patients did not have CAV (NonCAV). Mean rate-pressure products (RPP) at rest were significantly elevated in CAV patients compared to those without CAV (P = 0.002) but it was insignificant at stress (P = NS). There was no significant difference in the stress MBFs between CAV and NonCAV patients (P = NS). However, the difference in RPP-normalized stress MBFs was significant (P = 0.045), while RPP-normalized MBFs at rest was not significant (P = NS). Both CFR and RPP-normalized CFR were significantly lower in CAV compared to NonCAV patients (P < 0.001). There were significant correlations between MBFs and RPPs at rest for both CAV (ρ = 0.764, P = 0.047) and NonCAV patients (ρ = 0.641, P = 0.017), while there were no correlations at stress for CAV (ρ = 0.232, P = NS) and NonCAV patients (ρ = 0.068, P = NS). This study indicates that the resting MBF is higher in late-term post-OHT patients. The high resting MBF and reduced CFR suggest an unprecedented demand of blood flow and blunted response to stress due to impaired vasodilatory capacity that is exacerbated by the presence of CAV.