Deep learning based automated delineation of the intraprostatic gross tumour volume in PSMA-PET for patients with primary prostate cancer.

PURPOSE: With the increased use of focal radiation dose escalation for primary prostate cancer (PCa), accurate delineation of gross tumor volume (GTV) in prostate-specific membrane antigen PET (PSMA-PET) becomes crucial. Manual approaches are time-consuming and observer dependent. The purpose of this study was to create a deep learning model for the accurate delineation of the intraprostatic GTV in PSMA-PET. METHODS: A 3D U-Net was trained on 128 different 18F-PSMA-1007 PET images from three different institutions. Testing was done on 52 patients including one independent internal cohort (Freiburg: n = 19) and three independent external cohorts (Dresden: n = 14 18F-PSMA-1007, Boston: Massachusetts General Hospital (MGH): n = 9 18F-DCFPyL-PSMA and Dana-Farber Cancer Institute (DFCI): n = 10 68Ga-PSMA-11). Expert contours were generated in consensus using a validated technique. CNN predictions were compared to expert contours using Dice similarity coefficient (DSC). Co-registered whole-mount histology was used for the internal testing cohort to assess sensitivity/specificity. RESULTS: Median DSCs were Freiburg: 0.82 (IQR: 0.73-0.88), Dresden: 0.71 (IQR: 0.53-0.75), MGH: 0.80 (IQR: 0.64-0.83) and DFCI: 0.80 (IQR: 0.67-0.84), respectively. Median sensitivity for CNN and expert contours were 0.88 (IQR: 0.68-0.97) and 0.85 (IQR: 0.75-0.88) (p = 0.40), respectively. GTV volumes did not differ significantly (p > 0.1 for all comparisons). Median specificity of 0.83 (IQR: 0.57-0.97) and 0.88 (IQR: 0.69-0.98) were observed for CNN and expert contours (p = 0.014), respectively. CNN prediction took 3.81 seconds on average per patient. CONCLUSION: The CNN was trained and tested on internal and external datasets as well as histopathology reference, achieving a fast GTV segmentation for three PSMA-PET tracers with high diagnostic accuracy comparable to manual experts.

Observer studies of image quality of denoising reduced-count cardiac single photon emission computed tomography myocardial perfusion imaging by three-dimensional Gaussian post-reconstruction filtering and deep learning.

BACKGROUND: The aim of this research was to asses perfusion-defect detection-accuracy by human observers as a function of reduced-counts for 3D Gaussian post-reconstruction filtering vs deep learning (DL) denoising to determine if there was improved performance with DL. METHODS: SPECT projection data of 156 normally interpreted patients were used for these studies. Half were altered to include hybrid perfusion defects with defect presence and location known. Ordered-subset expectation-maximization (OSEM) reconstruction was employed with the optional correction of attenuation (AC) and scatter (SC) in addition to distance-dependent resolution (RC). Count levels varied from full-counts (100%) to 6.25% of full-counts. The denoising strategies were previously optimized for defect detection using total perfusion deficit (TPD). Four medical physicist (PhD) and six physician (MD) observers rated the slices using a graphical user interface. Observer ratings were analyzed using the LABMRMC multi-reader, multi-case receiver-operating-characteristic (ROC) software to calculate and compare statistically the area-under-the-ROC-curves (AUCs). RESULTS: For the same count-level no statistically significant increase in AUCs for DL over Gaussian denoising was determined when counts were reduced to either the 25% or 12.5% of full-counts. The average AUC for full-count OSEM with solely RC and Gaussian filtering was lower than for the strategies with AC and SC, except for a reduction to 6.25% of full-counts, thus verifying the utility of employing AC and SC with RC. CONCLUSION: We did not find any indication that at the dose levels investigated and with the DL network employed, that DL denoising was superior in AUC to optimized 3D post-reconstruction Gaussian filtering.

Mesh modeling of system geometry and anatomy phantoms for realistic GATE simulations and their inclusion in SPECT reconstruction.

Objective.Monte-Carlo simulation studies have been essential for advancing various developments in single photon emission computed tomography (SPECT) imaging, such as system design and accurate image reconstruction. Among the simulation software available, Geant4 application for tomographic emission (GATE) is one of the most used simulation toolkits in nuclear medicine, which allows building systems and attenuation phantom geometries based on the combination of idealized volumes. However, these idealized volumes are inadequate for modeling free-form shape components of such geometries. Recent GATE versions alleviate these major limitations by allowing users to import triangulated surface meshes.Approach.In this study, we describe our mesh-based simulations of a next-generation multi-pinhole SPECT system dedicated to clinical brain imaging, called AdaptiSPECT-C. To simulate realistic imaging data, we incorporated in our simulation the XCAT phantom, which provides an advanced anatomical description of the human body. An additional challenge with the AdaptiSPECT-C geometry is that the default voxelized XCAT attenuation phantom was not usable in our simulation due to intersection of objects of dissimilar materials caused by overlap of the air containing regions of the XCAT beyond the surface of the phantom and the components of the imaging system.Main results.We validated our mesh-based modeling against the one constructed by idealized volumes for a simplified single vertex configuration of AdaptiSPECT-C through simulated projection data of123I-activity distributions. We resolved the overlap conflict by creating and incorporating a mesh-based attenuation phantom following a volume hierarchy. We then evaluated our reconstructions with attenuation and scatter correction for projections obtained from simulation consisting of mesh-based modeling of the system and the attenuation phantom for brain imaging. Our approach demonstrated similar performance as the reference scheme simulated in air for uniform and clinical-like123I-IMP brain perfusion source distributions.Significance.This work enables the simulation of complex SPECT acquisitions and reconstructions for emulating realistic imaging data close to those of actual patients.

Improved Performance of a Multipinhole SPECT for DAT Imaging by Increasing Number of Pinholes at the Expense of Increased Multiplexing.

SPECT imaging of dopamine transporters (DAT) in the brain is a widely utilized study to improve the diagnosis of Parkinsonian syndromes, where conventional (parallel-hole and fan-beam) collimators on dual-head scanners are commonly employed. We have designed a multi-pinhole (MPH) collimator to improve the performance of DAT imaging. The MPH collimator focuses on the striatum and hence offers a better trade-off for sensitivity and spatial resolution than the conventional collimators within this clinically most relevant region for DAT imaging. Our original MPH design consisted of 9 pinholes with a background-to-striatal (Bkg/Str) projection multiplexing of 1% only. In this simulation study, we investigated whether further improvements in the performance of MPH imaging could be obtained by increasing the number of pinholes, hence by enhancing the sensitivity and sampling, despite the ambiguity in reconstructing images due to increased multiplexing. We performed analytic simulations of the MPH configurations with 9, 13, and 16 pinholes (aperture diameters: 4-6mm) using a digital phantom modeling DAT imaging. Our quantitative analyses indicated that using 13 (Bkg/Str: 12%) and 16 (Bkg/Str: 22%) pinholes provided better performance than the original 9-pinhole configuration for the acquisition with 2 or 4 angular views, but a similar performance with 8 and 16 views.

Robustness and performance of radiomic features in diagnosing cystic renal masses.

PURPOSE: We study the inter-reader variability in manual delineation of cystic renal masses (CRMs) presented in computerized tomography (CT) images and its effect on the classification performance of a machine learning algorithm in distinguishing benign from potentially malignant CRMs. In addition, we assessed whether the inclusion of higher-order robust radiomic features improves the classification performance over the use of first-order features. METHODS: 230 CRMs were independently delineated by two radiologists. Through a combination of random fluctuations, dilation, and erosion operations over the original region of interests (ROIs), we generated four additional sets of synthetic ROIs to capture the inter-reader variability realistically, as confirmed by dice coefficient measurements and visual assessment. We then identified the robust features based on the intra-class coefficient (ICC > 0.85) across these datasets. We applied a tenfold stratified cross-validation (CV) to train and test the performance of the random forest model for the classification of CRMs into benign and potentially malignant. RESULTS: The mean area under the curve (AUC), sensitivity, specificity, positive predictive value, and negative predictive value were 0.87, 0.82, 0.90, 0.85, and 0.93, respectively. With the usage of first-order features alone, the corresponding values were nearly identical. CONCLUSION: AUC ranged for the robust and uncorrelated features from 0.83 ± 0.09 to 0.93 ± 0.04 and for the first-order features from 0.84 ± 0.09 to 0.91 ± 0.04. Our study indicates that the first-order features alone are sufficient for the classification of CRMs, and that inclusion of higher-order features does not necessarily improve performance.

Retrospective fractional dose reduction in Tc-99m cardiac perfusion SPECT/CT patients: A human and model observer study.

BACKGROUND: In the ongoing efforts to reduce cardiac perfusion dose (injected radioactivity) for conventional SPECT/CT systems, we performed a human observer study to confirm our clinical model observer findings that iterative reconstruction employing OSEM (ordered-subset expectation-maximization) at 25% of the full dose (quarter-dose) has a similar performance for detection of hybrid cardiac perfusion defects as FBP at full dose. METHODS: One hundred and sixty-six patients, who underwent routine rest-stress Tc-99m sestamibi cardiac perfusion SPECT/CT imaging and clinically read as normally perfused, were included in the study. Ground truth was established by the normal read and the insertion of hybrid defects. In addition to the reconstruction of the 25% of full-dose data using OSEM with attenuation (AC), scatter (SC), and spatial resolution correction (RC), FBP and OSEM (with AC, SC, and RC) both at full dose (100%) were done. Both human observer and clinical model observer confidence scores were obtained to generate receiver operating characteristics (ROC) curves in a task-based image quality assessment. RESULTS: Average human observer AUC (area under the ROC curve) values of 0.725, 0.876, and 0.890 were obtained for FBP at full dose, OSEM at 25% of full dose, and OSEM at full dose, respectively. Both OSEM strategies were significantly better than FBP with P values of 0.003 and 0.01 respectively, while no significant difference was recorded between OSEM methods (P = 0.48). The clinical model observer results were 0.791, 0.822, and 0.879, respectively, for the same patient cases and processing strategies used in the human observer study. CONCLUSIONS: Cardiac perfusion SPECT/CT using OSEM reconstruction at 25% of full dose has AUCs larger than FBP and closer to those of full-dose OSEM when read by human observers, potentially replacing the higher dose studies during clinical reading.

Theranostics: The Role of Quantitative Nuclear Medicine Imaging.

Theranostics is a precision medicine discipline that integrates diagnostic nuclear medicine imaging with radionuclide therapy in a manner that provides both a tumor phenotype and personalized therapy to patients with cancer using radiopharmaceuticals aimed at the same target-specific biological pathway or receptor. The aim of quantitative nuclear medicine imaging is to plan the alpha or beta-emitting therapy based on an accurate 3-dimensional representation of the in-vivo distribution of radioactivity concentration within the tumor and normal organs/tissues in a noninvasive manner. In general, imaging may be either based on positron emission tomography (PET) or single photon emission computed tomography (SPECT) invariably in combination with X-ray CT (PET/CT; SPECT/CT) or, to a much lesser extent, MRI. PET and SPECT differ in terms of the radionuclides and physical processes that give rise to the emission of high energy photons, as well as the sets of technologies involved in their detection. Using a variety of standardized quantitative parameters, system calibration, patient preparation, imaging acquisition and reconstruction protocols, and image analysis protocols, an accurate quantification of the tracer distribution can be obtained, which helps prescribe the therapeutic dose for each patient.

Inclusion of quasi-vertex views in a brain-dedicated multi-pinhole SPECT system for improved imaging performance.

With brain-dedicated multi-detector systems employing pinhole apertures the usage of detectors facing the top of the patient's head (i.e. quasi-vertex (QV) views) can provide the advantage of additional viewing from close to the brain for improved detector coverage. In this paper, we report the results of simulation and reconstruction studies to investigate the impact of the QV views on the imaging performance of AdaptiSPECT-C, a brain-dedicated stationary SPECT system under development. In this design, both primary and scatter photons from regions located inferior to the brain can contribute to SPECT projections acquired by the QV views, and thus degrade AdaptiSPECT-C imaging performance. In this work, we determined the proportion, origin, and nature (i.e. primary, scatter, and multiple-scatter) of counts emitted from structures within the head and throughout the body contributing to projections from the different AdaptiSPECT-C detector rings, as well as from a true vertex view detector. We simulated phantoms used to assess different aspects of image quality (i.e. uniform activity concentration sphere, and Derenzo), as well as anthropomorphic phantoms with different count levels emulating clinical 123I activity distributions (i.e. DaTscan and perfusion). We determined that attenuation and scatter in the patient's body greatly diminish the probability of the photons emitted outside the volume of interest reaching to detectors and being recorded within the 15% photopeak energy window. In addition, we demonstrated that the inclusion of the residual of such counts in the system acquisition does not degrade visual interpretation or quantitative analysis. The addition of the QV detectors improves volumetric sensitivity, angular sampling, and spatial resolution leading to significant enhancement in image quality, especially in the striato-thalamic and superior regions of the brain. Besides, the use of QV detectors improves the recovery of clinically relevant metrics such as the striatal binding ratio and mean activity in selected cerebral structures. Our findings proving the usefulness of the QV ring for brain imaging with 123I agents can be generalized to other commonly used SPECT imaging agents labelled with isotopes, such as 99mTc and likely 111In.

The clinical utilities of multi-pinhole single photon emission computed tomography.

Single photon emission computed tomography (SPECT) is an important imaging modality for various applications in nuclear medicine. The use of multi-pinhole (MPH) collimators can provide superior resolution-sensitivity trade-off when imaging small field-of-view compared to conventional parallel-hole and fan-beam collimators. Besides the very successful application in small animal imaging, there has been a resurgence of the use of MPH collimators for clinical cardiac and brain studies, as well as other small field-of-view applications. This article reviews the basic principles of MPH collimators and introduces currently available and proposed clinical MPH SPECT systems.

Investigation of Axial and Angular Sampling in Multi-Detector Pinhole-SPECT Brain Imaging.

We designed a dedicated multi-detector multi-pinhole brain SPECT scanner to generate images of higher quality compared to general-purpose systems. The system, AdaptiSPECT-C, is intended to adapt its sensitivity-resolution trade-off by varying its aperture configurations allowing both high-sensitivity dynamic and high-spatial-resolution static imaging. The current system design consists of 23 detector heads arranged in a truncated spherical geometry. In this work, we investigated the axial and angular sampling capability of the current stationary system design. Two data acquisition schemes using limited rotation of the gantry and two others using axial translation of the imaging bed were also evaluated concerning their impact on image quality through improved sampling. Increasing both angular and axial sampling in the current prototype system resulted in quantitative improvements in image quality metrics and qualitative appearance of the images as determined in studies with specifically selected phantoms. Visual improvements for the brain phantoms with clinical distributions were less pronounced but presented quantitative improvements in the fidelity (normalized root-mean-square error (NRMSE)) and striatal specific binding ratio (SBR) for a dopamine transporter (DAT) distribution, and in NRMSE and activity recovery for a brain perfusion distribution. More pronounced improvements with increased sampling were seen in contrast recovery coefficient, bias, and coefficient of variation for a lesion in the brain perfusion distribution. The negligible impact of the most cranial ring of detectors on axial sampling, but its significant impact on sensitivity and angular sampling in the cranial portion of the imaging volume-of-interest were also determined.

Primary, scatter, and penetration characterizations of parallel-hole and pinhole collimators for I-123 SPECT.

Multi-pinhole (MPH) collimators are known to provide better trade-off between sensitivity and resolution for preclinical, as well as for smaller regions in clinical SPECT imaging compared to conventional collimators. In addition to this geometric advantage, MPH plates typically offer better stopping power for penetration than the conventional collimators, which is especially relevant for I-123 imaging. The I-123 emits a series of high-energy (>300 keV, ~2.5% abundance) gamma photons in addition to the primary emission (159 keV, 83% abundance). Despite their low abundance, high-energy photons penetrate through a low-energy parallel-hole (LEHR) collimator much more readily than the 159 keV photons, resulting in large downscatter in the photopeak window. In this work, we investigate the primary, scatter, and penetration characteristics of a single pinhole collimator that is commonly used for I-123 thyroid imaging and our two MPH collimators designed for I-123 DaTscan imaging for Parkinson's Disease, in comparison to three different parallel-hole collimators through a series of experiments and Monte Carlo simulations. The simulations of a point source and a digital human phantom with DaTscan activity distribution showed that our MPH collimators provide superior count performance in terms of high primary counts, low penetration, and low scatter counts compared to the parallel-hole and single pinhole collimators. For example, total scatter, multiple scatter, and collimator penetration events for the LEHR were 2.5, 7.6 and 14 times more than that of MPH within the 15% photopeak window. The total scatter fraction for LEHR was 56% where the largest contribution came from the high-energy scatter from the back compartments (31%). For the same energy window, the total scatter for MPH was 21% with only 1% scatter from the back compartments. We therefore anticipate that using MPH collimators, higher quality reconstructions can be obtained in a substantially shorter acquisition time for I-123 DaTscan and thyroid imaging.

Performance analysis of a high-sensitivity multi-pinhole cardiac SPECT system with hemi-ellipsoid detectors.

PURPOSE: Single-photon emission computed tomography (SPECT) is a noninvasive imaging modality, used in myocardial perfusion imaging. The challenges facing the majority of clinical SPECT systems are low sensitivity, poor resolution, and the relatively high radiation dose to the patient. New generation systems (GE Discovery, DSPECT) dedicated to cardiac imaging improve sensitivity by a factor of 5-8. This improvement can be used to decrease acquisition time and/or dose. However, in the case of ultra-low dose (~3 mCi) injections, acquisition times are still significantly long, taking 10-12 min. The purpose of this work is to investigate a new gamma camera design with 21 hemi-ellipsoid detectors each with a pinhole collimator for cardiac SPECT for further improvement in sensitivity and resolution and reduced patient exposures and imaging times. METHODS: To evaluate the resolution of our hemi-ellipsoid system, GATE Monte-Carlo simulations were performed on point-sources, rod-sources, and NCAT phantoms. For average full-width-half-maximum (FWHM) equivalence with base flat-detector, the pinhole-diameter for the curved hemi-ellipsoid detector was found to be 8.68 mm, an operating pinhole-diameter nominally expected to be ~3 times more sensitive than state-of-the-art systems. Rod-sources equally spaced within the region of interest were acquired with a 21-detector system and reconstructed with our multi-pinhole (MPH) iterative OSEM algorithm with collimator resolution recovery. The results were compared with the results of a state-of-the-art system (GE Discovery) available in the literature. The system was also evaluated using the mathematical anthropomorphic NCAT (NURBS-based Cardiac Torso; Segars et al. IEEE Trans Nucl Sci. 1999;46:503-506) phantom with a full (clinical)-dose acquisition (25 mCi) for 2 min and an ultra-low dose acquisition of 3 mCi for 5.44 min. The estimated left ventricle (LV) counts were compared with the available literature on a state-of-the-art system (DSPECT). FWHM of the LV wall on MPH-OSEM-reconstructed images with collimator resolution recovery was estimated. RESULTS: On acquired rod-sources, the average resolution (FWHM) after reconstruction with resolution recovery in the entire region of interest (ROI) for cardiac imaging was on the average 4.44 mm (±2.84), compared to 6.9 mm (±1 mm) reported for GE Discovery (Kennedy et al., J Nucl Cardiol. 2014:21:443-452). For NCAT studies, improved sensitivity allowed a full-dose (25 mCi) 2-min acquisition (Ell8.68mmFD) which yielded 3.79 M LV counts. This is ~3.35 times higher compared to 1.13 M LV counts acquired in 2 min for clinical full dose for state-of-the-art DSPECT. The increased sensitivity also allowed an ultra-low dose acquisition protocol (Ell8.68 mmULD), 3 mCi (eight times less injected dose) in 5.44 min. This ultra-low dose protocol yielded ~1.23 M LV counts which was comparable to the full-dose 2-min acquisition for DSPECT. The estimated NCAT average FWHM at the LV wall after 12 iterations of the OSEM reconstruction was 4.95 and 5.66 mm around the mid-short-axis slices for Ell8.68mmFD and Ell8.68mmULD, respectively. CONCLUSION: Our Monte-Carlo simulation studies and reconstruction suggest using (inverted wineglass sized) hemi-ellipsoid detectors with pinhole collimators can increase the sensitivity ~3.35 times over the new generation of dedicated cardiac SPECT systems, while also improving the reconstructed resolution for rod-sources with an average of 4.44 mm in region of interest. The extra sensitivity may be used for ultra-low dose imaging (3 mCi) at ~5.44 min for comparable clinical counts as state-of-the-art systems.

Simulations of a Multi-Pinhole SPECT Collimator for Clinical Dopamine Transporter (DAT) Imaging.

SPECT imaging of the dopamine transporter (DAT) is used for diagnosis and monitoring progression of Parkinson's Disease (PD), and differentiation of PD from other neurological disorders. The diagnosis is based on the DAT binding in the caudate and putamen structures in the striatum. We previously proposed a relatively inexpensive method to improve the detection and quantification of these structures for dual-head SPECT by replacing one of the fan-beam collimators with a specially designed multi-pinhole (MPH) collimator. In this work, we developed a realistic model of the proposed MPH system using the GATE simulation package and verified the geometry with an analytic simulator. Point source projections from these simulations closely matched confirming the accuracy of the pinhole geometries. The reconstruction of a hot-rod phantom showed that 4.8 mm resolution is achievable. The reconstructions of the XCAT brain phantom showed clear separation of the putamen and caudate, which is expected to improve the quantification of DAT imaging and PD diagnosis. Using this GATE model, point spread functions modeling physical factors will be generated for use in reconstruction. Also, further improvements in geometry are being investigated to increase the sensitivity of this base system while maintaining a target spatial resolution of 4.5-5 mm.

Correction of hysteretic respiratory motion in SPECT myocardial perfusion imaging: Simulation and patient studies.

PURPOSE: Amplitude-based respiratory gating is known to capture the extent of respiratory motion (RM) accurately but results in residual motion in the presence of respiratory hysteresis. In our previous study, we proposed and developed a novel approach to account for respiratory hysteresis by applying the Bouc-Wen (BW) model of hysteresis to external surrogate signals of anterior/posterior motion of the abdomen and chest with respiration. In this work, using simulated and clinical SPECT myocardial perfusion imaging (MPI) studies, we investigate the effects of respiratory hysteresis and evaluate the benefit of correcting it using the proposed BW model in comparison with the abdomen signal typically employed clinically. METHODS: The MRI navigator data acquired in free-breathing human volunteers were used in the specially modified 4D NCAT phantoms to allow simulating three types of respiratory patterns: monotonic, mild hysteresis, and strong hysteresis with normal myocardial uptake, and perfusion defects in the anterior, lateral, inferior, and septal locations of the mid-ventricular wall. Clinical scans were performed using a Tc-99m sestamibi MPI protocol while recording respiratory signals from thoracic and abdomen regions using a visual tracking system (VTS). The performance of the correction using the respiratory signals was assessed through polar map analysis in phantom and 10 clinical studies selected on the basis of having substantial RM. RESULTS: In phantom studies, simulations illustrating normal myocardial uptake showed significant differences (P < 0.001) in the uniformity of the polar maps between the RM uncorrected and corrected. No significant differences were seen in the polar map uniformity across the RM corrections. Studies simulating perfusion defects showed significantly decreased errors (P < 0.001) in defect severity and extent for the RM corrected compared to the uncorrected. Only for the strong hysteretic pattern, there was a significant difference (P < 0.001) among the RM corrections. The errors in defect severity and extent for the RM correction using abdomen signal were significantly higher compared to that of the BW (severity = -4.0%, P < 0.001; extent = -65.4%, P < 0.01) and chest (severity = -4.1%, P < 0.001; extent = -52.5%, P < 0.01) signals. In clinical studies, the quantitative analysis of the polar maps demonstrated qualitative and quantitative but not statistically significant differences (P = 0.73) between the correction methods that used the BW signal and the abdominal signal. CONCLUSIONS: This study shows that hysteresis in respiration affects the extent of residual motion left in the RM-binned data, which can impact wall uniformity and the visualization of defects. Thus, there appears to be the potential for improved accuracy in reconstruction in the presence of hysteretic RM with the BW model method providing a possible step in the direction of improvement.

Design of a Multi-Pinhole Collimator for I-123 DaTscan Imaging on Dual-Headed SPECT Systems in Combination with a Fan-Beam Collimator.

For the 2011 FDA approved Parkinson's Disease (PD) SPECT imaging agent I-123 labeled DaTscan, the volume of interest (VOI) is the interior portion of the brain. However imaging of the occipital lobe is also required with PD for calculation of the striatal binding ratio (SBR), a parameter of significance in early diagnosis, differentiation of PD from other disorders with similar clinical presentations, and monitoring progression. Thus we propose the usage of a combination of a multi-pinhole (MPH) collimator on one head of the SPECT system and a fan-beam on the other. The MPH would be designed to provide high resolution and sensitivity for imaging of the interior portion of the brain. The fan-beam collimator would provide lower resolution but complete sampling of the brain addressing data sufficiency and allowing a volume-of-interest to be defined over the occipital lobe for calculation of SBR's. Herein we focus on the design of the MPH component of the combined system. Combined reconstruction will be addressed in a subsequent publication. An analysis of 46 clinical DaTscan studies was performed to provide information to define the VOI, and design of a MPH collimator to image this VOI. The system spatial resolution for the MPH was set to 4.7 mm, which is comparable to that of clinical PET systems, and significantly smaller than that of fan-beam collimators employed in SPECT. With this set, we compared system sensitivities for three aperture array designs, and selected the 3 × 3 array due to it being the highest of the three. The combined sensitivity of the apertures for it was similar to that of an ultra-high resolution fan-beam (LEUHRF) collimator, but smaller than that of a high-resolution fan-beam collimator (LEHRF). On the basis of these results we propose the further exploration of this design through simulations, and the development of combined MPH and fan-beam reconstruction.

Comparison of the scanning linear estimator (SLE) and ROI methods for quantitative SPECT imaging.

In quantitative emission tomography, tumor activity is typically estimated from calculations on a region of interest (ROI) identified in the reconstructed slices. In these calculations, unpredictable bias arising from the null functions of the imaging system affects ROI estimates. The magnitude of this bias depends upon the tumor size and location. In prior work it has been shown that the scanning linear estimator (SLE), which operates on the raw projection data, is an unbiased estimator of activity when the size and location of the tumor are known. In this work, we performed analytic simulation of SPECT imaging with a parallel-hole medium-energy collimator. Distance-dependent system spatial resolution and non-uniform attenuation were included in the imaging simulation. We compared the task of activity estimation by the ROI and SLE methods for a range of tumor sizes (diameter: 1-3 cm) and activities (contrast ratio: 1-10) added to uniform and non-uniform liver backgrounds. Using the correct value for the tumor shape and location is an idealized approximation to how task estimation would occur clinically. Thus we determined how perturbing this idealized prior knowledge impacted the performance of both techniques. To implement the SLE for the non-uniform background, we used a novel iterative algorithm for pre-whitening stationary noise within a compact region. Estimation task performance was compared using the ensemble mean-squared error (EMSE) as the criterion. The SLE method performed substantially better than the ROI method (i.e. EMSE(SLE) was 23-174 times lower) when the background is uniform and tumor location and size are known accurately. The variance of the SLE increased when a non-uniform liver texture was introduced but the EMSE(SLE) continued to be 5-20 times lower than the ROI method. In summary, SLE outperformed ROI under almost all conditions that we tested.

Adaptation of the modified Bouc-Wen model to compensate for hysteresis in respiratory motion for the list-mode binning of cardiac SPECT and PET acquisitions: testing using MRI.

PURPOSE: Binning list-mode acquisitions as a function of a surrogate signal related to respiration has been employed to reduce the impact of respiratory motion on image quality in cardiac emission tomography (SPECT and PET). Inherent in amplitude binning is the assumption that there is a monotonic relationship between the amplitude of the surrogate signal and respiratory motion of the heart. This assumption is not valid in the presence of hysteresis when heart motion exhibits a different relationship with the surrogate during inspiration and expiration. The purpose of this study was to investigate the novel approach of using the Bouc-Wen (BW) model to provide a signal accounting for hysteresis when binning list-mode data with the goal of thereby improving motion correction. The study is based on the authors' previous observations that hysteresis between chest and abdomen markers was indicative of hysteresis between abdomen markers and the internal motion of the heart. METHODS: In 19 healthy volunteers, they determined the internal motion of the heart and diaphragm in the superior-inferior direction during free breathing using MRI navigators. A visual tracking system (vts) synchronized with MRI acquisition tracked the anterior-posterior motions of external markers placed on the chest and abdomen. These data were employed to develop and test the Bouc-Wen model by inputting the vts derived chest and abdomen motions into it and using the resulting output signals as surrogates for cardiac motion. The data of the volunteers were divided into training and testing sets. The training set was used to obtain initial values for the model parameters for all of the volunteers in the set, and for set members based on whether they were or were not classified as exhibiting hysteresis using a metric derived from the markers. These initial parameters were then employed with the testing set to estimate output signals. Pearson's linear correlation coefficient between the abdomen, chest, average of chest and abdomen markers, and Bouc-Wen derived signals versus the true internal motion of the heart from MRI was used to judge the signals match to the heart motion. RESULTS: The results show that the Bouc-Wen model generated signals demonstrated strong correlation with the heart motion. This correlation was slightly larger on average than that of the external surrogate signals derived from the abdomen marker, and average of the abdomen and chest markers, but was not statistically significantly different from them. CONCLUSIONS: The results suggest that the proposed model has the potential to be a unified framework for modeling hysteresis in respiratory motion in cardiac perfusion studies and beyond.

Digital anthropomorphic phantoms of non-rigid human respiratory and voluntary body motion for investigating motion correction in emission imaging.

The development of methods for correcting patient motion in emission tomography has been receiving increased attention. Often the performance of these methods is evaluated through simulations using digital anthropomorphic phantoms, such as the commonly used extended cardiac torso (XCAT) phantom, which models both respiratory and cardiac motion based on human studies. However, non-rigid body motion, which is frequently seen in clinical studies, is not present in the standard XCAT phantom. In addition, respiratory motion in the standard phantom is limited to a single generic trend. In this work, to obtain a more realistic representation of motion, we developed a series of individual-specific XCAT phantoms, modeling non-rigid respiratory and non-rigid body motions derived from the magnetic resonance imaging (MRI) acquisitions of volunteers. Acquisitions were performed in the sagittal orientation using the Navigator methodology. Baseline (no motion) acquisitions at end-expiration were obtained at the beginning of each imaging session for each volunteer. For the body motion studies, MRI was again acquired only at end-expiration for five body motion poses (shoulder stretch, shoulder twist, lateral bend, side roll, and axial slide). For the respiratory motion studies, an MRI was acquired during free/regular breathing. The magnetic resonance slices were then retrospectively sorted into 14 amplitude-binned respiratory states, end-expiration, end-inspiration, six intermediary states during inspiration, and six during expiration using the recorded Navigator signal. XCAT phantoms were then generated based on these MRI data by interactive alignment of the organ contours of the XCAT with the MRI slices using a graphical user interface. Thus far we have created five body motion and five respiratory motion XCAT phantoms from the MRI acquisitions of six healthy volunteers (three males and three females). Non-rigid motion exhibited by the volunteers was reflected in both respiratory and body motion phantoms with a varying extent and character for each individual. In addition to these phantoms, we recorded the position of markers placed on the chest of the volunteers for the body motion studies, which could be used as external motion measurement. Using these phantoms and external motion data, investigators will be able to test their motion correction approaches for realistic motion obtained from different individuals. The non-uniform rational B-spline data and the parameter files for these phantoms are freely available for downloading and can be used with the XCAT license.

Use of MRI to assess the prediction of heart motion with gross body motion in myocardial perfusion imaging by stereotracking of markers on the body surface.

PURPOSE: The aim of this study is to determine using MRI in volunteers whether the rigid-body-motion (RBM) model can be approximately used to estimate the gross body-motion of the heart from that of external markers on patient's chest. Our target clinical application is to use a visual-tracking-system (VTS) which employs stereoimaging to estimate heart motion during SPECT/CT and PET∕CT myocardial perfusion imaging. METHODS: To investigate body-motion separate from the respiration the authors had the volunteers hold their breath during the acquisition of a sequence of two sets of EKG-triggered MRI sagittal slices. The first set was acquired pre-motion, and the second postmotion. The motion of the heart within each breath-hold set of slices was estimated by registration to the semiautomatic 3D segmentation of the heart region in a baseline set acquired using the Navigator technique. The motion of the heart between the pre- and postmotion sets was then determined as the difference in the individual motions in comparison to the Navigator sets. An analysis of the combined motion of the individual markers on the chest was used to obtain an estimate of the six-degree-of-freedom RBM from the VTS system. The metric for judging agreement between the motion estimated by MRI and the VTS was the average error. This was defined as the average of the magnitudes of the differences in the vector displacements of all voxels in the heart region. Studies with the Data Spectrum Anthropomorphic Phantom and "No-Motion" studies in which the volunteer did not intentionally move were used to establish a baseline for agreement. With volunteer studies a t-test was employed to determine when statistically significant differences in Average Errors occurred compared to the No-motion studies. RESULTS: For phantom acquisitions, the Average Error when the motion was just translation was 0.1 mm. With complex motions, which included a combination of rotations and translations, the Average Error increased to 3.6 mm. In the volunteers the Average Error averaged over all No-Motion acquisitions was 1.0 mm. For the case of translational motion, which might be expected to be RBM, the Average Error averaged over all volunteer studies increased to 2.6 mm, which was statistically different from the No-Motion studies. For the case of bends and twists of the torso, which would be expected to challenge the RBM model, the Average Error averaged over all such volunteer studies was 4.9 mm and was again statistically different. Investigations of motion of the arm including just bending at the elbow and leg motion resulted in Average Errors which were not statistically different from the No-Motion studies. However, when shoulder movement was included with arm motion the Average Error was near that of torso bends and twists, and statistically different. CONCLUSIONS: Use of the RBM model with VTS predictions of heart motion during reconstruction should decrease the extent of artifacts for the types of patient motion studied. The impact of correction would be less for torso bends and twists, and arm motion which includes the shoulders.

Comparison of methods of acquiring attenuation maps for cardiac SPECT in the presence of respiratory motion.

BACKGROUND: We investigated the effect of the respiratory motion on attenuation-corrected (AC) SPECT images for three different SPECT systems, each using a different approach in obtaining attenuation maps: scanning-line sources (SLS) acquired simultaneous with emission; slow cone-beam CT (CBCT) acquired sequentially to emission; and fast helical CT (HCT) acquired sequentially to emission. METHODS: A torso phantom filled with (99m)Tc was used to model a cardiac perfusion study. Stationary baseline acquisitions were followed by the acquisitions with the phantom moving axially using a computer-controlled motion platform to simulate breathing. In addition, HCT acquisitions were made simulating breath-hold at different extents of misalignment between CT and emission. HCT images were also used to simulate the Average-CT method. Acquisitions were repeated with added breast attachments, and the heart insert in two different orientations. Visual comparison was made of AC maps, AC emission slices and polar maps. Quantitative comparisons were made of global uniformity based on the percent fractional standard deviation (%FSD) of the polar map segment values, and the ratio of the segment values in the Anterior and Inferior walls divided by that of the Lateral and Septal walls (AI/LS ratio). RESULTS: The AC maps for the SLS were inferior to the CT's, and most impacted by added large breast attachment. Motion artifacts seen on CBCT slices were minimized in the derived attenuation maps. AC maps obtained from HCT showed inconsistent organ sizes depending on the direction of respiration at the time of acquisition. Both visually and quantitatively CBCT resulted in the best uniformity (up to 3.4 % lower in %FSD) for all the stationary acquisitions, and for the motion acquisition of the female phantom with large breast attachment (up to 4.0 % lower). For the motion acquisition of the male phantoms, HCT resulted in slightly better uniformity (<0.5 % lower) than CBCT. Breath-hold at end-expiration slightly improved (up to 1.1 %) the uniformity over the HCT acquired during regular breathing. Further improvement was achieved with the Average-CT method. For all the systems, phantom respiratory motion reduced the AI/LS ratio compared to when the phantoms were stationary. CONCLUSIONS: The CBCT approach resulted in the best uniformity of the AC emission images. For the female phantom with larger breast attachment, HCT and SLS were truncated at some projection angles introducing artifacts into the AC emission images. The emission image artifacts observed with HCT could be mitigated by performing breath-hold acquisition at end-expiration or Average-CT type acquisitions.