Feasibility of using abbreviated scan protocols with population-based input functions for accurate kinetic modeling of [18F]-FDG datasets from a long axial FOV PET scanner.

BACKGROUND: Accurate kinetic modeling of 18F-fluorodeoxyglucose ([18F]-FDG) positron emission tomography (PET) data requires accurate knowledge of the available tracer concentration in the plasma during the scan time, known as the arterial input function (AIF). The gold standard method to derive the AIF requires collection of serial arterial blood samples, but the introduction of long axial field of view (LAFOV) PET systems enables the use of non-invasive image-derived input functions (IDIFs) from large blood pools such as the aorta without any need for bed movement. However, such protocols require a prolonged dynamic PET acquisition, which is impractical in a busy clinical setting. Population-based input functions (PBIFs) have previously shown potential in accurate Patlak analysis of [18F]-FDG datasets and can enable the use of shortened dynamic imaging protocols. Here, we exploit the high sensitivity and temporal resolution of a LAFOV PET system and explore the use of PBIF with abbreviated protocols in [18F]-FDG total body kinetic modeling. METHODS: Dynamic PET data were acquired in 24 oncological subjects for 65 min following the administration of [18F]-FDG. IDIFs were extracted from the descending thoracic aorta, and a PBIF was generated from 16 datasets. Five different scaled PBIFs (sPBIFs) were generated by scaling the PBIF with the AUC of IDIF curve tails using various portions of image data (35-65, 40-65, 45-65, 50-65, and 55-65 min post-injection). The sPBIFs were compared with the IDIFs using the AUCs and Patlak Ki estimates in tumor lesions and cerebral gray matter. Patlak plot start time (t*) was also varied to evaluate the performance of shorter acquisitions on the accuracy of Patlak Ki estimates. Patlak Ki estimates with IDIF and t* = 35 min were used as reference, and mean bias and precision (standard deviation of bias) were calculated to assess the relative performance of different sPBIFs. A comparison of parametric images generated using IDIF and sPBIFs was also performed. RESULTS: There was no statistically significant difference between AUCs of the IDIF and sPBIFs (Wilcoxon test: P > 0.05). Excellent agreement was shown between Patlak Ki estimates obtained using sPBIF and IDIF. Using the sPBIF55-65 with the Patlak model, 20 min of PET data (i.e., 45 to 65 min post-injection) achieved < 15% precision error in Ki estimates in tumor lesions compared to the estimates with the IDIF. Parametric images reconstructed using the IDIF and sPBIFs with and without an abbreviated protocol were visually comparable. Using Patlak Ki generated with an IDIF and 30 min of PET data as reference, Patlak Ki images generated using sPBIF55-65 with 20 min of PET data (t* = 45 min) provided excellent image quality with structural similarity index measure > 0.99 and peak signal-to-noise ratio > 55 dB. CONCLUSION: We demonstrate the feasibility of performing accurate [18F]-FDG Patlak analysis using sPBIFs with only 20 min of PET data from a LAFOV PET scanner.

Dynamic whole-body PET imaging: principles, potentials and applications.

PURPOSE: In this article, we discuss dynamic whole-body (DWB) positron emission tomography (PET) as an imaging tool with significant clinical potential, in relation to conventional standard uptake value (SUV) imaging. BACKGROUND: DWB PET involves dynamic data acquisition over an extended axial range, capturing tracer kinetic information that is not available with conventional static acquisition protocols. The method can be performed within reasonable clinical imaging times, and enables generation of multiple types of PET images with complementary information in a single imaging session. Importantly, DWB PET can be used to produce multi-parametric images of (i) Patlak slope (influx rate) and (ii) intercept (referred to sometimes as "distribution volume"), while also providing (iii) a conventional 'SUV-equivalent' image for certain protocols. RESULTS: We provide an overview of ongoing efforts (primarily focused on FDG PET) and discuss potential clinically relevant applications. CONCLUSION: Overall, the framework of DWB imaging [applicable to both PET/CT(computed tomography) and PET/MRI (magnetic resonance imaging)] generates quantitative measures that may add significant value to conventional SUV image-derived measures, with limited pitfalls as we also discuss in this work.

The effect of time-of-flight and point spread function modeling on 82Rb myocardial perfusion imaging of obese patients.

BACKGROUND: The effect of time-of-flight (TOF) and point spread function (PSF) modeling in image reconstruction has not been well studied for cardiac PET. This study assesses their separate and combined influence on 82Rb myocardial perfusion imaging in obese patients. METHODS: Thirty-six obese patients underwent rest-stress 82Rb cardiac PET. Images were reconstructed with and without TOF and PSF modeling. Perfusion was quantitatively compared using the AHA 17-segment model for patients grouped by BMI, cross-sectional body area in the scanner field of view, gender, and left ventricular myocardial volume. Summed rest scores (SRS), summed stress scores (SSS), and summed difference scores (SDS) were compared. RESULTS: TOF improved polar map visual uniformity and increased septal wall perfusion by up to 10%. This increase was greater for larger patients, more evident for patients grouped by cross-sectional area than by BMI, and more prominent for females. PSF modeling increased perfusion by about 1.5% in all cardiac segments. TOF modeling generally decreased SRS and SSS with significant decreases between 2.4 and 3.0 (P < .05), which could affect risk stratification; SDS remained about the same. With PSF modeling, SRS, SSS, and SDS were largely unchanged. CONCLUSION: TOF and PSF modeling affect regional and global perfusion, SRS, and SSS. Clinicians should consider these effects and gender-dependent differences when interpreting 82Rb perfusion studies.

A multi-centre evaluation of eleven clinically feasible brain PET/MRI attenuation correction techniques using a large cohort of patients.

AIM: To accurately quantify the radioactivity concentration measured by PET, emission data need to be corrected for photon attenuation; however, the MRI signal cannot easily be converted into attenuation values, making attenuation correction (AC) in PET/MRI challenging. In order to further improve the current vendor-implemented MR-AC methods for absolute quantification, a number of prototype methods have been proposed in the literature. These can be categorized into three types: template/atlas-based, segmentation-based, and reconstruction-based. These proposed methods in general demonstrated improvements compared to vendor-implemented AC, and many studies report deviations in PET uptake after AC of only a few percent from a gold standard CT-AC. Using a unified quantitative evaluation with identical metrics, subject cohort, and common CT-based reference, the aims of this study were to evaluate a selection of novel methods proposed in the literature, and identify the ones suitable for clinical use. METHODS: In total, 11 AC methods were evaluated: two vendor-implemented (MR-ACDIXON and MR-ACUTE), five based on template/atlas information (MR-ACSEGBONE (Koesters et al., 2016), MR-ACONTARIO (Anazodo et al., 2014), MR-ACBOSTON (Izquierdo-Garcia et al., 2014), MR-ACUCL (Burgos et al., 2014), and MR-ACMAXPROB (Merida et al., 2015)), one based on simultaneous reconstruction of attenuation and emission (MR-ACMLAA (Benoit et al., 2015)), and three based on image-segmentation (MR-ACMUNICH (Cabello et al., 2015), MR-ACCAR-RiDR (Juttukonda et al., 2015), and MR-ACRESOLUTE (Ladefoged et al., 2015)). We selected 359 subjects who were scanned using one of the following radiotracers: [18F]FDG (210), [11C]PiB (51), and [18F]florbetapir (98). The comparison to AC with a gold standard CT was performed both globally and regionally, with a special focus on robustness and outlier analysis. RESULTS: The average performance in PET tracer uptake was within ±5% of CT for all of the proposed methods, with the average±SD global percentage bias in PET FDG uptake for each method being: MR-ACDIXON (-11.3±3.5)%, MR-ACUTE (-5.7±2.0)%, MR-ACONTARIO (-4.3±3.6)%, MR-ACMUNICH (3.7±2.1)%, MR-ACMLAA (-1.9±2.6)%, MR-ACSEGBONE (-1.7±3.6)%, MR-ACUCL (0.8±1.2)%, MR-ACCAR-RiDR (-0.4±1.9)%, MR-ACMAXPROB (-0.4±1.6)%, MR-ACBOSTON (-0.3±1.8)%, and MR-ACRESOLUTE (0.3±1.7)%, ordered by average bias. The overall best performing methods (MR-ACBOSTON, MR-ACMAXPROB, MR-ACRESOLUTE and MR-ACUCL, ordered alphabetically) showed regional average errors within ±3% of PET with CT-AC in all regions of the brain with FDG, and the same four methods, as well as MR-ACCAR-RiDR, showed that for 95% of the patients, 95% of brain voxels had an uptake that deviated by less than 15% from the reference. Comparable performance was obtained with PiB and florbetapir. CONCLUSIONS: All of the proposed novel methods have an average global performance within likely acceptable limits (±5% of CT-based reference), and the main difference among the methods was found in the robustness, outlier analysis, and clinical feasibility. Overall, the best performing methods were MR-ACBOSTON, MR-ACMAXPROB, MR-ACRESOLUTE and MR-ACUCL, ordered alphabetically. These methods all minimized the number of outliers, standard deviation, and average global and local error. The methods MR-ACMUNICH and MR-ACCAR-RiDR were both within acceptable quantitative limits, so these methods should be considered if processing time is a factor. The method MR-ACSEGBONE also demonstrates promising results, and performs well within the likely acceptable quantitative limits. For clinical routine scans where processing time can be a key factor, this vendor-provided solution currently outperforms most methods. With the performance of the methods presented here, it may be concluded that the challenge of improving the accuracy of MR-AC in adult brains with normal anatomy has been solved to a quantitatively acceptable degree, which is smaller than the quantification reproducibility in PET imaging.

Inter-pass motion correction for whole-body dynamic PET and parametric imaging.

Whole-body dynamic FDG-PET imaging through continuous-bed-motion (CBM) mode multi-pass acquisition protocol is a promising metabolism measurement. However, inter-pass misalignment originating from body movement could degrade parametric quantification. We aim to apply a non-rigid registration method for inter-pass motion correction in whole-body dynamic PET. 27 subjects underwent a 90-min whole-body FDG CBM PET scan on a Biograph mCT (Siemens Healthineers), acquiring 9 over-the-heart single-bed passes and subsequently 19 CBM passes (frames). The inter-pass motion correction was executed using non-rigid image registration with multi-resolution, B-spline free-form deformations. The parametric images were then generated by Patlak analysis. The overlaid Patlak slope Ki and y-intercept Vb images were visualized to qualitatively evaluate motion impact and correction effect. The normalized weighted mean squared Patlak fitting errors (NFE) were compared in the whole body, head, and hypermetabolic regions of interest (ROI). In Ki images, ROI statistics were collected and malignancy discrimination capacity was estimated by the area under the receiver operating characteristic curve (AUC). After the inter-pass motion correction was applied, the spatial misalignment appearance between Ki and Vb images was successfully reduced. Voxel-wise normalized fitting error maps showed global error reduction after motion correction. The NFE in the whole body (p = 0.0013), head (p = 0.0021), and ROIs (p = 0.0377) significantly decreased. The visual performance of each hypermetabolic ROI in Ki images was enhanced, while 3.59% and 3.67% average absolute percentage changes were observed in mean and maximum Ki values, respectively, across all evaluated ROIs. The estimated mean Ki values had substantial changes with motion correction (p = 0.0021). The AUC of both mean Ki and maximum Ki after motion correction increased, possibly suggesting the potential of enhancing oncological discrimination capacity through inter-pass motion correction.

Diagnostic Performance of Dynamic Whole-Body Patlak [18F]FDG-PET/CT in Patients with Indeterminate Lung Lesions and Lymph Nodes.

BACKGROUND: Static [18F]FDG-PET/CT is the imaging method of choice for the evaluation of indeterminate lung lesions and NSCLC staging; however, histological confirmation of PET-positive lesions is needed in most cases due to its limited specificity. Therefore, we aimed to evaluate the diagnostic performance of additional dynamic whole-body PET. METHODS: A total of 34 consecutive patients with indeterminate pulmonary lesions were enrolled in this prospective trial. All patients underwent static (60 min p.i.) and dynamic (0-60 min p.i.) whole-body [18F]FDG-PET/CT (300 MBq) using the multi-bed-multi-timepoint technique (Siemens mCT FlowMotion). Histology and follow-up served as ground truth. Kinetic modeling factors were calculated using a two-compartment linear Patlak model (FDG influx rate constant = Ki, metabolic rate = MR-FDG, distribution volume = DV-FDG) and compared to SUV using ROC analysis. RESULTS: MR-FDGmean provided the best discriminatory power between benign and malignant lung lesions with an AUC of 0.887. The AUC of DV-FDGmean (0.818) and SUVmean (0.827) was non-significantly lower. For LNM, the AUCs for MR-FDGmean (0.987) and SUVmean (0.993) were comparable. Moreover, the DV-FDGmean in liver metastases was three times higher than in bone or lung metastases. CONCLUSIONS: Metabolic rate quantification was shown to be a reliable method to detect malignant lung tumors, LNM, and distant metastases at least as accurately as the established SUV or dual-time-point PET scans.

Effect of Scan Time on Oncologic Lesion Detection in Whole-Body PET.

Lesion-detection performance in oncologic PET depends in part upon count statistics, with shorter scans having higher noise and reduced lesion detectability. However, advanced techniques such as time-of-flight (TOF) and point spread function (PSF) modeling can improve lesion detection. This work investigates the relationship between reducing count levels (as a surrogate for scan time) and reconstructing with PSF model and TOF. A series of twenty-four whole-body phantom scans was acquired on a Biograph mCT TOF PET/CT scanner using the experimental methodology prescribed for the Utah PET Lesion Detection Database. Six scans were acquired each day over four days, with up to 23 (68)Ge shell-less lesions (diam. 6, 8, 10, 12, 16 mm) distributed throughout the phantom thorax and pelvis. Each scan acquired 6 bed positions at 240 s/bed in listmode format. The listmode files were then statistically pruned, preserving Poisson statistics, to equivalent count levels for scan times of 180 s, 120 s, 90 s, 60 s, 45 s, 30 s, and 15 s per bed field-of-view, corresponding to whole-body scan times of 1.5-24 min. Each dataset was reconstructed using ordinary Poisson line-of-response (LOR) OSEM, with PSF model, with TOF, and with PSF+TOF. Localization receiver operating characteristics (LROC) analysis was then performed using the channelized non-prewhitened (CNPW) observer. The results were analyzed to delineate the relationship between scan time, reconstruction method, and strength of post-reconstruction filter. Lesion-detection performance degraded as scan time was reduced, and progressively stronger filters were required to maximize performance for the shorter scans. PSF modeling and TOF were found to improve detection performance, but the degree of improvement for TOF was much larger than for PSF for the large phantom used in this study. Notably, the images using TOF provided equivalent lesion-detection performance to the images without TOF for scan durations 40% shorter, suggesting that TOF may offset, at least in part, the need for longer scan times in larger patients.

Unsupervised inter-frame motion correction for whole-body dynamic PET using convolutional long short-term memory in a convolutional neural network.

Subject motion in whole-body dynamic PET introduces inter-frame mismatch and seriously impacts parametric imaging. Traditional non-rigid registration methods are generally computationally intense and time-consuming. Deep learning approaches are promising in achieving high accuracy with fast speed, but have yet been investigated with consideration for tracer distribution changes or in the whole-body scope. In this work, we developed an unsupervised automatic deep learning-based framework to correct inter-frame body motion. The motion estimation network is a convolutional neural network with a combined convolutional long short-term memory layer, fully utilizing dynamic temporal features and spatial information. Our dataset contains 27 subjects each under a 90-min FDG whole-body dynamic PET scan. Evaluating performance in motion simulation studies and a 9-fold cross-validation on the human subject dataset, compared with both traditional and deep learning baselines, we demonstrated that the proposed network achieved the lowest motion prediction error, obtained superior performance in enhanced qualitative and quantitative spatial alignment between parametric Ki and Vb images, and significantly reduced parametric fitting error. We also showed the potential of the proposed motion correction method for impacting downstream analysis of the estimated parametric images, improving the ability to distinguish malignant from benign hypermetabolic regions of interest. Once trained, the motion estimation inference time of our proposed network was around 460 times faster than the conventional registration baseline, showing its potential to be easily applied in clinical settings.

Benefit of Improved Performance with State-of-the Art Digital PET/CT for Lesion Detection in Oncology.

The latest digital whole-body PET scanners provide a combination of higher sensitivity and improved spatial and timing resolution. We performed a lesion detectability study on two generations of Biograph PET/CT scanners, the mCT Flow and the Vision, to study the impact of improved physical performance on clinical performance. Our hypothesis was that the improved performance of the Vision would result in improved lesion detectability, allowing shorter imaging times or, equivalently, a lower injected dose. Methods: Data were acquired with the Society of Nuclear Medicine and Molecular Imaging Clinical Trials Network torso phantom combined with a 20-cm-diameter cylindrical phantom. Spherical lesions were emulated by acquiring sphere-in-air data and combining them with the phantom data to generate combined datasets with embedded lesions of known contrast. Two sphere sizes and uptakes were used: 9.89-mm-diameter spheres with 6:1 (lung) and 3:1 (cylinder) local activity concentration uptakes and 4.95-mm-diameter spheres with 9.6:1 (lung) and 4.5:1 (cylinder) local activity concentration uptakes. Standard image reconstruction was performed: an ordinary Poisson ordered-subsets expectation maximization algorithm with point-spread function and time-of-flight modeling and postreconstruction smoothing with a 5-mm gaussian filter. The Vision images were also generated without any postreconstruction smoothing. Generalized scan statistics methodology was used to estimate the area under the localized receiver-operating-characteristic curve (ALROC). Results: The higher sensitivity and improved time-of-flight performance of the Vision leads to reduced contrast in the background noise nodule distribution. Measured lesion contrast is also higher on the Vision because of its improved spatial resolution. Hence, the ALROC is noticeably higher for the Vision than for the mCT Flow. Conclusion: Improved overall performance of the Vision provides a factor of 4-6 reduction in imaging time (or injected dose) over the mCT Flow when using the ALROC metric for lesions at least 9.89 mm in diameter. Smaller lesions are barely detected in the mCT Flow, leading to even higher ALROC gains with the Vision. The improved spatial resolution of the Vision also leads to a higher measured contrast that is closer to the real uptake, implying improved quantification. Postreconstruction smoothing, however, reduces this improvement in measured contrast, thereby reducing the ALROC for small, high-uptake lesions.

Impact of time-of-flight on PET tumor detection.

UNLABELLED: Time-of-flight (TOF) PET uses very fast detectors to improve localization of events along coincidence lines-of-response. This information is then utilized to improve the tomographic reconstruction. This work evaluates the effect of TOF upon an observer's performance for detecting and localizing focal warm lesions in noisy PET images. METHODS: An advanced anthropomorphic lesion-detection phantom was scanned 12 times over 3 days on a prototype TOF PET/CT scanner (Siemens Medical Solutions). The phantom was devised to mimic whole-body oncologic (18)F-FDG PET imaging, and a number of spheric lesions (diameters 6-16 mm) were distributed throughout the phantom. The data were reconstructed with the baseline line-of-response ordered-subsets expectation-maximization algorithm, with the baseline algorithm plus point spread function model (PSF), baseline plus TOF, and with both PSF+TOF. The lesion-detection performance of each reconstruction was compared and ranked using localization receiver operating characteristics (LROC) analysis with both human and numeric observers. The phantom results were then subjectively compared to 2 illustrative patient scans reconstructed with PSF and with PSF+TOF. RESULTS: Inclusion of TOF information provides a significant improvement in the area under the LROC curve compared to the baseline algorithm without TOF data (P = 0.002), providing a degree of improvement similar to that obtained with the PSF model. Use of both PSF+TOF together provided a cumulative benefit in lesion-detection performance, significantly outperforming either PSF or TOF alone (P < 0.002). Example patient images reflected the same image characteristics that gave rise to improved performance in the phantom data. CONCLUSION: Time-of-flight PET provides a significant improvement in observer performance for detecting focal warm lesions in a noisy background. These improvements in image quality can be expected to improve performance for the clinical tasks of detecting lesions and staging disease. Further study in a large clinical population is warranted to assess the benefit of TOF for various patient sizes and count levels, and to demonstrate effective performance in the clinical environment.

Data-driven voluntary body motion detection and non-rigid event-by-event correction for static and dynamic PET.

PET has the potential to perform absolute in vivo radiotracer quantitation. This potential can be compromised by voluntary body motion (BM), which degrades image resolution, alters apparent tracer uptakes, introduces CT-based attenuation correction mismatch artifacts and causes inaccurate parameter estimates in dynamic studies. Existing body motion correction (BMC) methods include frame-based image-registration (FIR) approaches and real-time motion tracking using external measurement devices. FIR does not correct for motion occurring within a pre-defined frame and the device-based method is generally not practical in routine clinical use, since it requires attaching a tracking device to the patient and additional device set up time. In this paper, we proposed a data-driven algorithm, centroid of distribution (COD), to detect BM. In this algorithm, the central coordinate of the time-of-flight (TOF) bin, which can be used as a reasonable surrogate for the annihilation point, is calculated for every event, and averaged over a certain time interval to generate a COD trace. We hypothesized that abrupt changes on the COD trace in lateral direction represent BMs. After detection, BM is estimated using non-rigid image registrations and corrected through list-mode reconstruction. The COD-based BMC approach was validated using a monkey study and was evaluated against FIR using four human and one dog studies with multiple tracers. The proposed approach successfully detected BMs and yielded superior correction results over conventional FIR approaches.

Respiratory Motion Compensation for PET/CT with Motion Information Derived from Matched Attenuation-Corrected Gated PET Data.

Respiratory motion degrades the detection and quantification capabilities of PET/CT imaging. Moreover, mismatch between a fast helical CT image and a time-averaged PET image due to respiratory motion results in additional attenuation correction artifacts and inaccurate localization. Current motion compensation approaches typically have 3 limitations: the mismatch among respiration-gated PET images and the CT attenuation correction (CTAC) map can introduce artifacts in the gated PET reconstructions that can subsequently affect the accuracy of the motion estimation; sinogram-based correction approaches do not correct for intragate motion due to intracycle and intercycle breathing variations; and the mismatch between the PET motion compensation reference gate and the CT image can cause an additional CT-mismatch artifact. In this study, we established a motion correction framework to address these limitations. Methods: In the proposed framework, the combined emission-transmission reconstruction algorithm was used for phase-matched gated PET reconstructions to facilitate the motion model building. An event-by-event nonrigid respiratory motion compensation method with correlations between internal organ motion and external respiratory signals was used to correct both intracycle and intercycle breathing variations. The PET reference gate was automatically determined by a newly proposed CT-matching algorithm. We applied the new framework to 13 human datasets with 3 different radiotracers and 323 lesions and compared its performance with CTAC and non-attenuation correction (NAC) approaches. Validation using 4-dimensional CT was performed for one lung cancer dataset. Results: For the 10 18F-FDG studies, the proposed method outperformed (P < 0.006) both the CTAC and the NAC methods in terms of region-of-interest-based SUVmean, SUVmax, and SUV ratio improvements over no motion correction (SUVmean: 19.9% vs. 14.0% vs. 13.2%; SUVmax: 15.5% vs. 10.8% vs. 10.6%; SUV ratio: 24.1% vs. 17.6% vs. 16.2%, for the proposed, CTAC, and NAC methods, respectively). The proposed method increased SUV ratios over no motion correction for 94.4% of lesions, compared with 84.8% and 86.4% using the CTAC and NAC methods, respectively. For the 2 18F-fluoropropyl-(+)-dihydrotetrabenazine studies, the proposed method reduced the CT-mismatch artifacts in the lower lung where the CTAC approach failed and maintained the quantification accuracy of bone marrow where the NAC approach failed. For the 18F-FMISO study, the proposed method outperformed both the CTAC and the NAC methods in terms of motion estimation accuracy at 2 lung lesion locations. Conclusion: The proposed PET/CT respiratory event-by-event motion-correction framework with motion information derived from matched attenuation-corrected PET data provides image quality superior to that of the CTAC and NAC methods for multiple tracers.

Development and Evaluation of mini-EXPLORER: A Long Axial Field-of-View PET Scanner for Nonhuman Primate Imaging.

We describe a long axial field-of-view (FOV) PET scanner for high-sensitivity and total-body imaging of nonhuman primates and present the physical performance and first phantom and animal imaging results. Methods: The mini-EXPLORER PET scanner was built using the components of a clinical scanner reconfigured with a detector ring diameter of 43.5 cm and an axial length of 45.7 cm. National Electrical Manufacturers Association (NEMA) NU-2 and NU-4 phantoms were used to measure sensitivity and count rate performance. Reconstructed spatial resolution was investigated by imaging a radially stepped point source and a Derenzo phantom. The effect of the wide acceptance angle was investigated by comparing performance with maximum acceptance angles of 14°-46°. Lastly, an initial assessment of the in vivo performance of the mini-EXPLORER was undertaken with a dynamic 18F-FDG nonhuman primate (rhesus monkey) imaging study. Results: The NU-2 total sensitivity was 5.0%, and the peak noise-equivalent count rate measured with the NU-4 monkey scatter phantom was 1,741 kcps, both obtained using the maximum acceptance angle (46°). The NU-4 scatter fraction was 16.5%, less than 1% higher than with a 14° acceptance angle. The reconstructed spatial resolution was approximately 3.0 mm at the center of the FOV, with a minor loss in axial spatial resolution (0.5 mm) when the acceptance angle increased from 14° to 46°. The rhesus monkey 18F-FDG study demonstrated the benefit of the high sensitivity of the mini-EXPLORER, including fast imaging (1-s early frames), excellent image quality (30-s and 5-min frames), and late-time-point imaging (18 h after injection), all obtained at a single bed position that captured the major organs of the rhesus monkey. Conclusion: This study demonstrated the physical performance and imaging capabilities of a long axial FOV PET scanner designed for high-sensitivity imaging of nonhuman primates. Further, the results of this study suggest that a wide acceptance angle can be used with a long axial FOV scanner to maximize sensitivity while introducing only minor trade-offs such as a small increase in scatter fraction and slightly degraded axial spatial resolution.

Spatiotemporal distribution modeling of PET tracer uptake in solid tumors.

OBJECTIVE: Distribution of PET tracer uptake is elaborately modeled via a general equation used for solute transport modeling. This model can be used to incorporate various transport parameters of a solid tumor such as hydraulic conductivity of the microvessel wall, transvascular permeability as well as interstitial space parameters. This is especially significant because tracer delivery and drug delivery to solid tumors are determined by similar underlying tumor transport phenomena, and quantifying the former can enable enhanced prediction of the latter. METHODS: We focused on the commonly utilized FDG PET tracer. First, based on a mathematical model of angiogenesis, the capillary network of a solid tumor and normal tissues around it were generated. The coupling mathematical method, which simultaneously solves for blood flow in the capillary network as well as fluid flow in the interstitium, is used to calculate pressure and velocity distributions. Subsequently, a comprehensive spatiotemporal distribution model (SDM) is applied to accurately model distribution of PET tracer uptake, specifically FDG in this work, within solid tumors. RESULTS: The different transport mechanisms, namely convention and diffusion from vessel to tissue and in tissue, are elaborately calculated across the domain of interest and effect of each parameter on tracer distribution is investigated. The results show the convection terms to have negligible effect on tracer transport and the SDM can be solved after eliminating these terms. CONCLUSION: The proposed framework of spatiotemporal modeling for PET tracers can be utilized to comprehensively assess the impact of various parameters on the spatiotemporal distribution of PET tracers.

Generalized PSF modeling for optimized quantitation in PET imaging.

Point-spread function (PSF) modeling offers the ability to account for resolution degrading phenomena within the PET image generation framework. PSF modeling improves resolution and enhances contrast, but at the same time significantly alters image noise properties and induces edge overshoot effect. Thus, studying the effect of PSF modeling on quantitation task performance can be very important. Frameworks explored in the past involved a dichotomy of PSF versus no-PSF modeling. By contrast, the present work focuses on quantitative performance evaluation of standard uptake value (SUV) PET images, while incorporating a wide spectrum of PSF models, including those that under- and over-estimate the true PSF, for the potential of enhanced quantitation of SUVs. The developed framework first analytically models the true PSF, considering a range of resolution degradation phenomena (including photon non-collinearity, inter-crystal penetration and scattering) as present in data acquisitions with modern commercial PET systems. In the context of oncologic liver FDG PET imaging, we generated 200 noisy datasets per image-set (with clinically realistic noise levels) using an XCAT anthropomorphic phantom with liver tumours of varying sizes. These were subsequently reconstructed using the OS-EM algorithm with varying PSF modelled kernels. We focused on quantitation of both SUVmean and SUVmax, including assessment of contrast recovery coefficients, as well as noise-bias characteristics (including both image roughness and coefficient of-variability), for different tumours/iterations/PSF kernels. It was observed that overestimated PSF yielded more accurate contrast recovery for a range of tumours, and typically improved quantitative performance. For a clinically reasonable number of iterations, edge enhancement due to PSF modeling (especially due to over-estimated PSF) was in fact seen to lower SUVmean bias in small tumours. Overall, the results indicate that exactly matched PSF modeling does not offer optimized PET quantitation, and that PSF overestimation may provide enhanced SUV quantitation. Furthermore, generalized PSF modeling may provide a valuable approach for quantitative tasks such as treatment-response assessment and prognostication.

Evaluation of MLACF based calculated attenuation brain PET imaging for FDG patient studies.

Calculating attenuation correction for brain PET imaging rather than using CT presents opportunities for low radiation dose applications such as pediatric imaging and serial scans to monitor disease progression. Our goal is to evaluate the iterative time-of-flight based maximum-likelihood activity and attenuation correction factors estimation (MLACF) method for clinical FDG brain PET imaging. FDG PET/CT brain studies were performed in 57 patients using the Biograph mCT (Siemens) four-ring scanner. The time-of-flight PET sinograms were acquired using the standard clinical protocol consisting of a CT scan followed by 10 min of single-bed PET acquisition. Images were reconstructed using CT-based attenuation correction (CTAC) and used as a gold standard for comparison. Two methods were compared with respect to CTAC: a calculated brain attenuation correction (CBAC) and MLACF based PET reconstruction. Plane-by-plane scaling was performed for MLACF images in order to fix the variable axial scaling observed. The noise structure of the MLACF images was different compared to those obtained using CTAC and the reconstruction required a higher number of iterations to obtain comparable image quality. To analyze the pooled data, each dataset was registered to a standard template and standard regions of interest were extracted. An SUVr analysis of the brain regions of interest showed that CBAC and MLACF were each well correlated with CTAC SUVrs. A plane-by-plane error analysis indicated that there were local differences for both CBAC and MLACF images with respect to CTAC. Mean relative error in the standard regions of interest was less than 5% for both methods and the mean absolute relative errors for both methods were similar (3.4% ± 3.1% for CBAC and 3.5% ± 3.1% for MLACF). However, the MLACF method recovered activity adjoining the frontal sinus regions more accurately than CBAC method. The use of plane-by-plane scaling of MLACF images was found to be a crucial step in order to obtain improved activity estimates. Presence of local errors in both MLACF and CBAC based reconstructions would require the use of a normal database for clinical assessment. However, further work is required in order to assess the clinical advantage of MLACF over CBAC based method.

Simultaneous reconstruction of the activity image and registration of the CT image in TOF-PET.

Previously, maximum-likelihood methods have been proposed to jointly estimate the activity image and the attenuation image or the attenuation sinogram from time-of-flight (TOF) positron emission tomography (PET) data. In this contribution, we propose a method that addresses the possible alignment problem of the TOF-PET emission data and the computed tomography (CT) attenuation data, by combining reconstruction and registration. The method, called MLRR, iteratively reconstructs the activity image while registering the available CT-based attenuation image, so that the pair of activity and attenuation images maximise the likelihood of the TOF emission sinogram. The algorithm is slow to converge, but some acceleration could be achieved by using Nesterov's momentum method and by applying a multi-resolution scheme for the non-rigid displacement estimation. The latter also helps to avoid local optima, although convergence to the global optimum cannot be guaranteed. The results are evaluated on 2D and 3D simulations as well as a respiratory gated clinical scan. Our experiments indicate that the proposed method is able to correct for possible misalignment of the CT-based attenuation image, and is therefore a very promising approach to suppressing attenuation artefacts in clinical PET/CT. When applied to respiratory gated data of a patient scan, it produced deformations that are compatible with breathing motion and which reduced the well known attenuation artefact near the dome of the liver. Since the method makes use of the energy-converted CT attenuation image, the scale problem of joint reconstruction is automatically solved.

Whole-body direct 4D parametric PET imaging employing nested generalized Patlak expectation-maximization reconstruction.

Whole-body (WB) dynamic PET has recently demonstrated its potential in translating the quantitative benefits of parametric imaging to the clinic. Post-reconstruction standard Patlak (sPatlak) WB graphical analysis utilizes multi-bed multi-pass PET acquisition to produce quantitative WB images of the tracer influx rate K i as a complimentary metric to the semi-quantitative standardized uptake value (SUV). The resulting K i images may suffer from high noise due to the need for short acquisition frames. Meanwhile, a generalized Patlak (gPatlak) WB post-reconstruction method had been suggested to limit K i bias of sPatlak analysis at regions with non-negligible (18)F-FDG uptake reversibility; however, gPatlak analysis is non-linear and thus can further amplify noise. In the present study, we implemented, within the open-source software for tomographic image reconstruction platform, a clinically adoptable 4D WB reconstruction framework enabling efficient estimation of sPatlak and gPatlak images directly from dynamic multi-bed PET raw data with substantial noise reduction. Furthermore, we employed the optimization transfer methodology to accelerate 4D expectation-maximization (EM) convergence by nesting the fast image-based estimation of Patlak parameters within each iteration cycle of the slower projection-based estimation of dynamic PET images. The novel gPatlak 4D method was initialized from an optimized set of sPatlak ML-EM iterations to facilitate EM convergence. Initially, realistic simulations were conducted utilizing published (18)F-FDG kinetic parameters coupled with the XCAT phantom. Quantitative analyses illustrated enhanced K i target-to-background ratio (TBR) and especially contrast-to-noise ratio (CNR) performance for the 4D versus the indirect methods and static SUV. Furthermore, considerable convergence acceleration was observed for the nested algorithms involving 10-20 sub-iterations. Moreover, systematic reduction in K i % bias and improved TBR were observed for gPatlak versus sPatlak. Finally, validation on clinical WB dynamic data demonstrated the clinical feasibility and superior K i CNR performance for the proposed 4D framework compared to indirect Patlak and SUV imaging.

Optimizing injected dose in clinical PET by accurately modeling the counting-rate response functions specific to individual patient scans.

UNLABELLED: To optimize the injected dose of radiopharmaceutical in PET, one needs to know its relationship to some metric of data quality for individual patient scans, such as noise-equivalent counting rate (NECR). In this paper, we show how one may accurately model the clinical NECR response corresponding to specific patient scans much as if a counting-rate test had been performed on each patient. We apply this technique to patient data and show how it can lead to improved clinical scanning protocols. METHODS: True and random coincidence rates expressed as functions of an appropriate measurable system parameter such as the detector single-event rate have functional forms that are largely independent of the object being scanned. Thus, reference true and random response functions may be scaled directly to the specific counting rates measured on a clinical scan, thereby yielding a curve of NECR versus injected dose. We have applied this technique to 2 groups of 163 clinical (18)F-FDG scans each. One of the groups was obtained on a lutetium oxyorthosilicate PET/CT scanner with conventional front-end electronics, and the other was obtained on a lutetium oxyorthosilicate PET/CT scanner with a new digital data processing system (Pico-3D). RESULTS: At 90%-95% of maximum signal-to-noise ratio (SNR), the mean optimal dose for a 60-min uptake period ranged from 366 to 717 MBq depending on the electronics and randoms processing method. There was only a slight (1 MBq/kg) dependence of optimal dose on patient weight but a larger dependence on position in the body. Pico-3D electronics improved optimal data SNR by 35% for a 70-kg person, but in both cases NECR fell rapidly with increasing weight (1.4%/kg). For an equivalent data SNR, a 120-kg person would have to be scanned 2.3 times longer than a 60-kg person. Over this range of weight, the mean scatter fraction increased by 12% whereas the ratio of mean randoms to trues increased by 48%. CONCLUSION: The methodology developed here allows one to directly estimate the optimal dose to inject for specific clinical scans and permits a detailed analysis of the sources of noise in PET data and of their variation with parameters such as patient weight.

Fourier rebinning of time-of-flight PET data.

This paper explores fast reconstruction strategies for 3D time-of-flight (TOF) positron emission tomography (PET), based on 2D data rebinning. Starting from pre-corrected 3D TOF data, a rebinning algorithm estimates for each transaxial slice the 2D TOF sinogram that would have been acquired by a single-ring scanner. The rebinned sinograms can then be reconstructed using any algorithm for 2D TOF reconstruction. We introduce TOF-FORE, an approximate rebinning algorithm obtained by extending the Fourier rebinning method for non-TOF data. In addition, we identify two partial differential equations that must be satisfied by consistent 3D TOF data, and use them to derive exact rebinning algorithms and to characterize the degree of the approximation in TOF-FORE. Numerical simulations demonstrate that TOF-FORE is more accurate than two different TOF extensions of the single-slice rebinning method, and suggest that TOF-FORE will be a valuable tool for practical TOF PET in the range of axial apertures and time resolutions typical of current scanners.