Anatomical-based partial volume correction for low-dose dedicated cardiac SPECT/CT

Partial volume correction for Lu-177-PSMA SPECT

BACKGROUND

The limited spatial resolution in SPECT images leads to partial volume effect (PVE), degrading the subsequent dosimetric accuracy. We aim to quantitatively evaluate PVE and partial volume corrections (PVC), i.e., recovery coefficient (RC)-PVC (RC-PVC), reblurred Van-Cittert (RVC) and iterative Yang (IY), in 177Lu-PSMA-617 SPECT images.

METHODS

We employed a geometrical cylindrical phantom containing five spheres (diameters ranging from 20 to 40 mm) and 40 XCAT phantoms with various anatomical variations and activity distributions. SIMIND Monte Carlo code was used to generate realistic noisy projections. In the clinical study, sequential quantitative SPECT/CT imaging at 4 time-points post 177Lu-PSMA-617 injections were analyzed for 10 patients. Iterative statistical reconstruction methods were used for reconstruction with attenuation, scatter and geometrical collimator detector response corrections, followed by post-filters. The RC-curves were fit based on the geometrical phantom study and applied for XCAT phantom and clinical study in RC-PVC. Matched and 0.5-2.0 voxels (2.54-10.16 mm) mismatched sphere masks were deployed in IY. The coefficient of variation (CoV) was measured on a uniform background on the geometrical phantom. RCs of spheres and mean absolute activity error (MAE) of kidneys and tumors were evaluated in simulation data, while the activity difference was evaluated in clinical data before and after PVC.

RESULTS

In the simulation study, the spheres experienced significant PVE, i.e., 0.26 RC and 0.70 RC for the 20 mm and 40 mm spheres, respectively. RVC and IY improved the RC of the 20 mm sphere to 0.37 and 0.75 and RC of the 40 mm sphere to 0.96 and 1.04. Mismatch in mask increased the activity error for all spheres in IY. RVC increased noise and caused Gibbs ringing artifacts. For XCAT phantoms, both RVC and IY performed comparably and were superior to RC-PVC in reducing the MAE of the kidneys. However, IY and RC-PVC outperformed RVC for tumors. The XCAT phantom study and clinical study showed a similar trend in the kidney and tumor activity differences between non-PVC and PVC.

CONCLUSIONS

PVE greatly impacts activity quantification, especially for small objects. All PVC methods improve the quantification accuracy in 177Lu-PSMA SPECT.

Sex differences in 123I-mIBG scintigraphy imaging techniques in patients with heart failure

BACKGROUND

123I-mIBG-scintigraphy could be a useful stratifying tool for patients with heart failure (HF). The purpose of this retrospective study is to evaluate whether there are differences between men and women with HF in terms of the prediction of cardiac arrhythmic events (AE).

RESEARCH AND METHODS

A total of 306 patients, before implantable-cardioverter-defibrillator (ICD) implantation, were evaluated. They underwent 123I-mIBG-scintigraphy and an evaluation of the results was performed after 85 months of follow-up. Early and late planar and SPECT cardiac images were acquired. Heart-to-mediastinum ratio (HM) for planar images and the sum of the segmental scores (SS) for SPECT were calculated.

RESULTS

In the general population, age, early SS (ESS), late SS (LSS), and ejection fraction (EF) were statistically significant for the prediction of AE at Cox regression, while early and late HM (eHM,lHM) were not significative for the prediction of AE. Population was divided into females and males and univariate analysis was conducted separately for the two cohorts: no significant variables for prediction of AE were found in females. For males, ESS, LSS, EF, and late HM were statistically significant predictors of AE. The overall survival was similar in males and females, but the risk of AE is lower in males than in females.

CONCLUSIONS

123I-mIBG represents a more effective tool for the prediction of AE in male patients than in women.

Activity recovery for differently shaped objects in quantitative SPECT

Objective.The aim was to theoretically and experimentally investigate recovery in SPECT images with objects of different shapes. Furthermore, the accuracy of volume estimation by thresholding was studied for those shapes.Approach.Nine spheres, nine oblate spheroids, and nine prolate spheroids phantom inserts were used, of which the six smaller spheres were part of the NEMA IEC body phantom and the rest of the inserts were 3D-printed. The inserts were filled with99mTc and177Lu. When filled with99mTc, SPECT images were acquired in a Siemens Symbia Intevo Bold gamma camera and when filled with177Lu in a General Electric NM/CT 870 DR gamma camera. The signal rate per activity (SRPA) was determined for all inserts and represented as a function of the volume-to-surface ratio and of the volume-equivalent radius using VOIs defined according to the sphere dimensions and VOIs defined using thresholding. Experimental values were compared with theoretical curves obtained analytically (spheres) or numerically (spheroids), starting from the convolution of a source distribution with a point-spread function. Validation of the activity estimation strategy was performed using four 3D-printed ellipsoids. Lastly, the threshold values necessary to determine the volume of each insert were obtained.Main results.Results showed that SRPA values for the oblate spheroids diverted from the other inserts, when SRPA were represented as a function of the volume-equivalent radius. However, SRPA values for all inserts followed a similar behaviour when represented as a function of the volume-to-surface ratio. Results for ellipsoids were in agreement with those results. For the three types of inserts the volume could be accurately estimated using a threshold method for volumes larger than 25 ml.Significance.Determination of SRPA independently of lesion or organ shape should decrease uncertainties in estimated activities and thereby, in the long term, be beneficial to patient care.

Segmentation-Free PVC for Cardiac SPECT Using a Densely-Connected Multi-Dimensional Dynamic Network

In nuclear imaging, limited resolution causes partial volume effects (PVEs) that affect image sharpness and quantitative accuracy. Partial volume correction (PVC) methods incorporating high-resolution anatomical information from CT or MRI have been demonstrated to be effective. However, such anatomical-guided methods typically require tedious image registration and segmentation steps. Accurately segmented organ templates are also hard to obtain, particularly in cardiac SPECT imaging, due to the lack of hybrid SPECT/CT scanners with high-end CT and associated motion artifacts. Slight mis-registration/mis-segmentation would result in severe degradation in image quality after PVC. In this work, we develop a deep-learning-based method for fast cardiac SPECT PVC without anatomical information and associated organ segmentation. The proposed network involves a densely-connected multi-dimensional dynamic mechanism, allowing the convolutional kernels to be adapted based on the input images, even after the network is fully trained. Intramyocardial blood volume (IMBV) is introduced as an additional clinical-relevant loss function for network optimization. The proposed network demonstrated promising performance on 28 canine studies acquired on a GE Discovery NM/CT 570c dedicated cardiac SPECT scanner with a 64-slice CT using Technetium-99m-labeled red blood cells. This work showed that the proposed network with densely-connected dynamic mechanism produced superior results compared with the same network without such mechanism. Results also showed that the proposed network without anatomical information could produce images with statistically comparable IMBV measurements to the images generated by anatomical-guided PVC methods, which could be helpful in clinical translation.

Post-reconstruction attenuation correction for SPECT myocardium perfusion imaging facilitated by deep learning-based attenuation map generation

BACKGROUND

Attenuation correction can improve the quantitative accuracy of single-photon emission computed tomography (SPECT) images. Existing SPECT-only systems normally can only provide non-attenuation corrected (NC) images which are susceptible to attenuation artifacts. In this work, we developed a post-reconstruction attenuation correction (PRAC) approach facilitated by a deep learning-based attenuation map for myocardial perfusion SPECT imaging.

METHODS

In the PRAC method, new projection data were estimated via forwardly projecting the scanner-generated NC image. Then an attenuation map, generated from NC image using a pretrained deep learning (DL) convolutional neural network, was incorporated into an offline reconstruction algorithm to obtain the attenuation-corrected images from the forwardly projected projections. We evaluated the PRAC method using 30 subjects with a DL network trained with 40 subjects, using the vendor-generated AC images and CT-based attenuation maps as the ground truth.

RESULTS

The PRAC methods using DL-generated and CT-based attenuation maps were both highly consistent with the scanner-generated AC image. The globally normalized mean absolute errors were 1.1% ± .6% and .7% ± .4% and the localized absolute percentage errors were 8.9% ± 13.4% and 7.8% ± 11.4% in the left ventricular (LV) blood pool, respectively, and - 1.3% ± 8.0% and - 3.8% ± 4.5% in the LV myocardium for PRAC methods using DL-generated and CT-based attenuation maps, respectively. The summed stress scores after PRAC using both attenuation maps were more consistent with the ground truth than those of the NC images.

CONCLUSION

We developed a PRAC approach facilitated by deep learning-based attenuation maps for SPECT myocardial perfusion imaging. It may be feasible for this approach to provide AC images for SPECT-only scanner data.

Nuclear Cardiac Imaging in the Interventional Suite

PURPOSE OF REVIEW

This review presents the current state of imaging approaches that enable real-time molecular imaging in the interventional suite and discusses the potential future use of integrated nuclear imaging and fluoroscopy for intraprocedural guidance in the evaluation and treatment of both cardiovascular and oncological diseases.

RECENT FINDINGS

Although there are no commercially available real-time hybrid nuclear imaging devices that are approved for use in the interventional suite, prototype open gantry hybrid nuclear imaging and x-ray c-arm imaging systems and theranostic catheter for location radiotracer detection are currently undergoing development and testing by multiple groups. The integration of physiological and molecular targeted nuclear imaging for real-time delivery of targeted theranostics in the interventional laboratory may enable more personalized care for a wide variety of cardiovascular procedures and improve patient outcomes.

Quantification of myocardial blood flow (MBF) and reserve (MFR) incorporated with a novel segmentation approach: Assessments of quantitative precision and the lower limit of normal MBF and MFR in patients

BACKGROUND

Quantification of myocardial blood flow (MBF) and myocardial flow reserve (MFR) has shown diagnostic and prognostic values for the assessment of coronary artery disease (CAD). This study aimed to evaluate in patients a highly automatic Yale-MQ (myocardial blood flow quantification) software incorporated with a novel image segmentation approach for quantification of global and regional MBF and MFR from dynamic 82Rb cardiac positron emission tomography (PET).

METHODS

Global and regional MBFs and MFRs were quantified in 80 patients (18 normal and 62 CAD subjects) by two different observers using the Yale-MQ software. Lower limits of normal (LLN) values and intra- and inter-observer variabilities of MBFs and MFRs were calculated for the assessment of quantitative precision. The Yale-MQ was compared with a commercially available software (Corridor 4DM) being used as a reference.

RESULTS

The Yale-MQ method provided precise assessments of LLNs of MBF and MFR. The global and regional MBFs and MFR quantified via Yale-MQ were correlated strongly with those via Corridor4DM (R ≥ 0.867). The intra- and inter-observer variabilities of MBFs and MFRs quantified via Yale-MQ were small (≤ 7.7% for MBFs and ≤ 10.0% for MFRs) with excellent correlations (R ≥ 0.980 for MBFs and R ≥ 0.976 for MFRs).

CONCLUSIONS

The new Yale-MQ software associated with the automatic processing scheme provides a highly reproducible clinical tool for precise quantification of MBF and MFR in patients with reliable LLN values.

Recent advances in cardiac SPECT instrumentation and imaging methods

Cardiac SPECT continues to play a critical role in detecting and managing cardiovascular disease, in particularly coronary artery disease (CAD) (Jaarsma et al 2012 J. Am. Coll. Cardiol. 59 1719-28), (Agostini et al 2016 Eur. J. Nucl. Med. Mol. Imaging 43 2423-32). While conventional dual-head SPECT scanners using parallel-hole collimators and scintillation crystals with photomultiplier tubes are still the workhorse of cardiac SPECT, they have the limitations of low photon sensitivity (~130 count s^-1 MBq^-1), poor image resolution (~15 mm) (Imbert et al 2012 J. Nucl. Med. 53 1897-903), relatively long acquisition time, inefficient use of the detector, high radiation dose, etc. Recently our field observed an exciting growth of new developments of dedicated cardiac scanners and collimators, as well as novel imaging algorithms for quantitative cardiac SPECT. These developments have opened doors to new applications with potential clinical impact, including ultra-low-dose imaging, absolute quantification of myocardial blood flow (MBF) and coronary flow reserve (CFR), multi-radionuclide imaging, and improved image quality as a result of attenuation, scatter, motion, and partial volume corrections (PVCs). In this article, we review the recent advances in cardiac SPECT instrumentation and imaging methods. This review mainly focuses on the most recent developments published since 2012 and points to the future of cardiac SPECT from an imaging physics perspective.

Cardiac motion correction for improving perfusion defect detection in cardiac SPECT at standard and reduced doses of activity

In cardiac SPECT perfusion imaging, cardiac motion can lead to motion blurring of anatomical detail and perfusion defects in the reconstructed myocardium. In this study, we investigated the potential benefit of cardiac motion correction for improving the detectability of perfusion defects. We considered a post-reconstruction motion correction (PMC) approach in which the image motion between two cardiac gates is obtained with optical flow estimation. In the experiments, we demonstrated the proposed post-reconstruction motion correction with optical flow estimation (PMC-OFE) approach on a set of clinical acquisitions from 194 subjects. We quantified the detectability of perfusion defects in the reconstructed images by using the total perfusion deficit scores, calculated by the clinical software tool QPS, and conducted a receiver-operating-characteristic (ROC) study to obtain the detection performance. Besides imaging with conventional standard dose, we also evaluated the approach for reduced dose SPECT imaging where the imaging dose was retrospectively reduced to 50%, 25%, and 12.5% of the standard dose. The proposed PMC-OFE approach achieved at each dose level higher area-under-the-ROC-curve (AUC) for perfusion defect detection than the traditional approach of using ungated data (Non-MC) (p -value  <  0.05); in particular, with half dose, PMC-OFE achieved AUC  =  0.813, which is comparable to Non-MC with standard dose (AUC  =  0.795). Moreover, the proposed PMC-OFE approach also outperformed the 'Motion Frozen' (MF) method implemented in the clinical quantitative gated SPECT (QGS) software. In particular, at 25% and 12.5% of standard dose, the AUC values obtained by PMC-OFE are 0.788 and 0.779, respectively, compared to 0.758 and 0.731 for MF (p -value  <  0.05).

Non-invasive measures of coronary microcirculation: Taking the long road to the clinic

Although coronary microvascular disease is now a well-recognized entity that is associated with significant morbidity and mortality, current non-invasive strategies cannot differentiate between coronary microvascular disease (CMD) and obstructive epicardial stenosis. While the evaluation of intramyocardial blood volume as a surrogate measure for microvascular health may have limited sensitivity in early-stage disease, this strategy does enable the diagnosis of CMD in the presence of concurrent epicardial disease, bringing us one step further toward improving the management of this disease. Herein, we discuss the advantages and limitations of current non-invasive measures of CMD and the need for further investment in bringing these technologies to the bedside.

Quantification of intramyocardial blood volume with 99mTc-RBC SPECT-CT imaging: A preclinical study

BACKGROUND

Currently, there is no established non-invasive imaging approach to directly evaluate myocardial microcirculatory function in order to diagnose microvascular disease independent of co-existing epicardial disease. In this work, we developed a methodological framework for quantification of intramyocardial blood volume (IMBV) as a novel index of microcirculatory function with SPECT/CT imaging of 99mTc-labeled red blood cells (RBCs).

METHODS

Dual-gated myocardial SPECT/CT equilibrium imaging of 99mTc-RBCs was performed on twelve canines under resting conditions. Five correction schemes were studied: cardiac gating with no other corrections (CG), CG with attenuation correction (CG + AC), CG + AC with scatter correction (CG + AC + SC), dual cardiorespiratory gating with AC + SC (DG + AC + SC), and DG + AC + SC with partial volume correction (DG + AC + SC + PVC). Quantification of IMBV using each approach was evaluated in comparison to those obtained from all corrections. The in vivo SPECT estimates of IMBV values were validated against those obtained from ex vivo microCT imaging of the casted hearts.

RESULTS

The estimated IMBV with all corrections was 0.15 ± 0.03 for the end-diastolic phase and 0.11 ± 0.03 for the end-systolic phase. The cycle-dependent change in IMBV (ΔIMBV) with all corrections was 23.9 ± 8.6%. Schemes that applied no correction or partial correction resulted in significant over-estimation of IMBV and significant under-underestimation of ΔIMBV. Estimates of IMBV and ΔIMBV using all corrections were consistent with values reported in the literature using invasive techniques. In vivo SPECT estimates of IMBV strongly correlated (R2 ≥ 0.70) with ex vivo measures for the various correction schemes, while the fully corrected scheme yielded the smallest bias.

CONCLUSIONS

Non-invasive quantification of IMBV is feasible using 99mTc-RBCs SPECT/CT imaging, however, requires full compensation of physical degradation factors.

Fully automatic multi-atlas segmentation of CTA for partial volume correction in cardiac SPECT/CT

Anatomical-based partial volume correction (PVC) has been shown to improve image quality and quantitative accuracy in cardiac SPECT/CT. However, this method requires manual segmentation of various organs from contrast-enhanced computed tomography angiography (CTA) data. In order to achieve fully automatic CTA segmentation for clinical translation, we investigated the most common multi-atlas segmentation methods. We also modified the multi-atlas segmentation method by introducing a novel label fusion algorithm for multiple organ segmentation to eliminate overlap and gap voxels. To evaluate our proposed automatic segmentation, eight canine ^99mTc-labeled red blood cell SPECT/CT datasets that incorporated PVC were analyzed, using the leave-one-out approach. The Dice similarity coefficient of each organ was computed. Compared to the conventional label fusion method, our proposed label fusion method effectively eliminated gaps and overlaps and improved the CTA segmentation accuracy. The anatomical-based PVC of cardiac SPECT images with automatic multi-atlas segmentation provided consistent image quality and quantitative estimation of intramyocardial blood volume, as compared to those derived using manual segmentation. In conclusion, our proposed automatic multi-atlas segmentation method of CTAs is feasible, practical, and facilitates anatomical-based PVC of cardiac SPECT/CT images.

Noise suppressed partial volume correction for cardiac SPECT/CT

PURPOSE

Partial volume correction (PVC) methods typically improve quantification at the expense of increased image noise and reduced reproducibility. In this study, the authors developed a novel voxel-based PVC method that incorporates anatomical knowledge to improve quantification while suppressing noise for cardiac SPECT/CT imaging.

METHODS

In the proposed method, the SPECT images were first reconstructed using anatomical-based maximum a posteriori (AMAP) with Bowsher's prior to penalize noise while preserving boundaries. A sequential voxel-by-voxel PVC approach (Yang's method) was then applied on the AMAP reconstruction using a template response. This template response was obtained by forward projecting a template derived from a contrast-enhanced CT image, and then reconstructed using AMAP to model the partial volume effects (PVEs) introduced by both the system resolution and the smoothing applied during reconstruction. To evaluate the proposed noise suppressed PVC (NS-PVC), the authors first simulated two types of cardiac SPECT studies: a (99m)Tc-tetrofosmin myocardial perfusion scan and a (99m)Tc-labeled red blood cell (RBC) scan on a dedicated cardiac multiple pinhole SPECT/CT at both high and low count levels. The authors then applied the proposed method on a canine equilibrium blood pool study following injection with (99m)Tc-RBCs at different count levels by rebinning the list-mode data into shorter acquisitions. The proposed method was compared to MLEM reconstruction without PVC, two conventional PVC methods, including Yang's method and multitarget correction (MTC) applied on the MLEM reconstruction, and AMAP reconstruction without PVC.

RESULTS

The results showed that the Yang's method improved quantification, however, yielded increased noise and reduced reproducibility in the regions with higher activity. MTC corrected for PVE on high count data with amplified noise, although yielded the worst performance among all the methods tested on low-count data. AMAP effectively suppressed noise and reduced the spill-in effect in the low activity regions. However it was unable to reduce the spill-out effect in high activity regions. NS-PVC yielded superior performance in terms of both quantitative assessment and visual image quality while improving reproducibility.

CONCLUSIONS

The results suggest that NS-PVC may be a promising PVC algorithm for application in low-dose protocols, and in gated and dynamic cardiac studies with low counts.

The impact of system matrix dimension on small FOV SPECT reconstruction with truncated projections

PURPOSE

A dedicated cardiac hybrid single photon emission computed tomography (SPECT)/CT scanner that uses cadmium zinc telluride detectors and multiple pinhole collimators for stationary acquisition offers many advantages. However, the impact of the reconstruction system matrix (SM) dimension on the reconstructed image quality from truncated projections and 19 angular samples acquired on this scanner has not been extensively investigated. In this study, the authors aimed to investigate the impact of the dimensions of SM and the use of body contour derived from adjunctive CT imaging as an object support in reconstruction on this scanner, in relation to background extracardiac activity.

METHODS

The authors first simulated a generic SPECT/CT system to image four NCAT phantoms with various levels of extracardiac activity and compared the reconstructions using SM in different dimensions and with/without body contour as a support for quantitative evaluations. The authors then compared the reconstructions of 18 patient studies, which were acquired on a GE Discovery NM570c scanner following injection of different radiotracers, including (99m)Tc-Tetrofosmin and (123)I-mIBG, comparing the scanner's default SM that incompletely covers the body with a large SM that incorporates a patient specific full body contour.

RESULTS

The simulation studies showed that the reconstructions using a SM that only partially covers the body yielded artifacts on the edge of the field of view (FOV), overestimation of activity and increased nonuniformity in the blood pool for the phantoms with higher relative levels of extracardiac activity. However, the impact on the quantitative accuracy in the high activity region, such as the myocardium, was subtle. On the other hand, an excessively large SM that enclosed the entire body alleviated the artifacts and reduced overestimation in the blood pool, but yielded slight underestimation in myocardium and defect regions. The reconstruction using the larger SM with body contour yielded the most quantitatively accurate results in all the regions of interest for a range of uptake levels in the extracardiac regions. In patient studies, the SM incorporating patient specific body contour minimized extracardiac artifacts, yielded similar myocardial activity, lower blood pool activity, and subsequently improved myocardium-to-blood pool contrast (p < 0.0001) by an average of 7% (range 0%-18%) across all the patients, compared to the reconstructions using the scanner's default SM.

CONCLUSIONS

Their results demonstrate that using a large SM that incorporates a CT derived body contour in the reconstruction could improve quantitative accuracy within the FOV for clinical studies with high extracardiac activity.