Effect of Using 2mm Voxels on Observer Performance for PET Lesion Detection.

Positron emission tomography (PET) images are typically reconstructed with an in-plane pixel size of approximately 4mm for cancer imaging. The objective of this work was to evaluate the effect of using smaller pixels on general oncologic lesion-detection. A series of observer studies was performed using experimental phantom data from the Utah PET Lesion Detection Database, which modeled whole-body FDG PET cancer imaging of a 92kg patient. The data comprised 24 scans over 4 days on a Biograph mCT time-of-flight (TOF) PET/CT scanner, with up to 23 lesions (diam. 6-16mm) distributed throughout the phantom each day. Images were reconstructed with 2.036mm and 4.073mm pixels using ordered-subsets expectation-maximization (OSEM) both with and without point spread function (PSF) modeling and TOF. Detection performance was assessed using the channelized non-prewhitened numerical observer with localization receiver operating characteristic (LROC) analysis. Tumor localization performance and the area under the LROC curve were then analyzed as functions of the pixel size. In all cases, the images with ~2mm pixels provided higher detection performance than those with ~4mm pixels. The degree of improvement from the smaller pixels was larger than that offered by PSF modeling for these data, and provided roughly half the benefit of using TOF. Key results were confirmed by two human observers, who read subsets of the test data. This study suggests that a significant improvement in tumor detection performance for PET can be attained by using smaller voxel sizes than commonly used at many centers. The primary drawback is a 4-fold increase in reconstruction time and data storage requirements.

Fourier domain closed-form formulas for estimation of kinetic parameters in reversible multi-compartment models.

BACKGROUND: Compared with static imaging, dynamic emission computed tomographic imaging with compartment modeling can quantify in vivo physiologic processes, providing useful information about molecular disease processes. Dynamic imaging involves estimation of kinetic rate parameters. For multi-compartment models, kinetic parameter estimation can be computationally demanding and problematic with local minima. METHODS: This paper offers a new perspective to the compartment model fitting problem where Fourier linear system theory is applied to derive closed-form formulas for estimating kinetic parameters for the two-compartment model. The proposed Fourier domain estimation method provides a unique solution, and offers very different noise response as compared to traditional non-linear chi-squared minimization techniques. RESULTS: The unique feature of the proposed Fourier domain method is that only low frequency components are used for kinetic parameter estimation, where the DC (i.e., the zero frequency) component in the data is treated as the most important information, and high frequency components that tend to be corrupted by statistical noise are discarded. Computer simulations show that the proposed method is robust without having to specify the initial condition. The resultant solution can be fine tuned using the traditional iterative method. CONCLUSIONS: The proposed Fourier-domain estimation method has closed-form formulas. The proposed Fourier-domain curve-fitting method does not require an initial condition, it minimizes a quadratic objective function and a closed-form solution can be obtained. The noise is easier to control, simply by discarding the high frequency components, and emphasizing the DC component.

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.

Quantitative evaluation of PSMA PET imaging using a realistic anthropomorphic phantom and shell-less radioactive epoxy lesions.

BACKGROUND: Positron emission tomography (PET) with prostate specific membrane antigen (PSMA) have shown superior performance in detecting metastatic prostate cancers. Relative to [18F]fluorodeoxyglucose ([18F]FDG) PET images, PSMA PET images tend to visualize significantly higher-contrast focal lesions. We aim to evaluate segmentation and reconstruction algorithms in this emerging context. Specifically, Bayesian or maximum a posteriori (MAP) image reconstruction, compared to standard ordered subsets expectation maximization (OSEM) reconstruction, has received significant interest for its potential to reach convergence with minimal noise amplifications. However, few phantom studies have evaluated the quantitative accuracy of such reconstructions for high contrast, small lesions (sub-10 mm) that are typically observed in PSMA images. In this study, we cast 3 mm-16-mm spheres using epoxy resin infused with a long half-life positron emitter (sodium-22; 22Na) to simulate prostate cancer metastasis. The anthropomorphic Probe-IQ phantom, which features a liver, bladder, lungs, and ureters, was used to model relevant anatomy. Dynamic PET acquisitions were acquired and images were reconstructed with OSEM (varying subsets and iterations) and BSREM (varying ß parameters), and the effects on lesion quantitation were evaluated. RESULTS: The 22Na lesions were scanned against an aqueous solution containing fluorine-18 (18F) as the background. Regions-of-interest were drawn with MIM Software using 40% fixed threshold (40% FT) and a gradient segmentation algorithm (MIM's PET Edge+). Recovery coefficients (RCs) (max, mean, peak, and newly defined "apex"), metabolic tumour volume (MTV), and total tumour uptake (TTU) were calculated for each sphere. SUVpeak and SUVapex had the most consistent RCs for different lesion-to-background ratios and reconstruction parameters. The gradient-based segmentation algorithm was more accurate than 40% FT for determining MTV and TTU, particularly for lesions [Formula: see text] 6 mm in diameter (R2 = 0.979-0.996 vs. R2 = 0.115-0.527, respectively). CONCLUSION: An anthropomorphic phantom was used to evaluate quantitation for PSMA PET imaging of metastatic prostate cancer lesions. BSREM with ß = 200-400 and OSEM with 2-5 iterations resulted in the most accurate and robust measurements of SUVmean, MTV, and TTU for imaging conditions in 18F-PSMA PET/CT images. SUVapex, a hybrid metric of SUVmax and SUVpeak, was proposed for robust, accurate, and segmentation-free quantitation of lesions for PSMA PET.

Design of an anthropomorphic PET phantom with elastic lungs and respiration modeling.

PURPOSE: Respiratory motion during positron emission tomography (PET) scans can be a major detriment to image quality in oncological imaging. The impact of motion on lesion quantification and detectability can be assessed using phantoms with realistic anatomy representation and motion modeling. In this work, we develop an anthropomorphic phantom for PET imaging that combines anatomic fidelity and a realistic breathing mechanism with deformable lungs. METHODS: We start from a previously developed anatomically accurate but static phantom of a human torso, and add elastic lungs with a highly controllable actuation mechanism which replicates the physics of breathing. The space outside the lungs is filled with a radioactive water solution. To maintain anatomical accuracy and realistic gamma ray attenuation in the torso, all motion mechanisms and actuators are positioned outside of the phantom compartment. The actuation mechanism can produce custom respiratory waveforms with breathing rates up to 25 breaths per minute and tidal volumes up to 1200 mL. RESULTS: Several tests were performed to validate the performance of the phantom assembly, in which the phantom was filled with water and given respiratory waveforms to execute. All parts demonstrated expected performance. Force requirements were not exceeded and no leaks were detected, although continued use of the phantom is required to evaluate wear. The motion of the lungs was determined to be within a reasonable realistic range. CONCLUSIONS: The full mechanical design is described in this paper, as well as a software application with graphical user interface which was developed to plan and visualize respiratory patterns. Both are available online as open source files. The developed phantom will facilitate future work in evaluating the impact of respiratory motion on lesion quantification and detectability in clinical practice.

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.

Rapid Multi-Tracer PET Tumor Imaging With F-FDG and Secondary Shorter-Lived Tracers.

Rapid multi-tracer PET, where two to three PET tracers are rapidly scanned with staggered injections, can recover certain imaging measures for each tracer based on differences in tracer kinetics and decay. We previously showed that single-tracer imaging measures can be recovered to a certain extent from rapid dual-tracer (62)Cu - PTSM (blood flow) + (62)Cu - ATSM (hypoxia) tumor imaging. In this work, the feasibility of rapidly imaging (18)F-FDG plus one or two of these shorter-lived secondary tracers was evaluated in the same tumor model. Dynamic PET imaging was performed in four dogs with pre-existing tumors, and the raw scan data was combined to emulate 60 minute long dual- and triple-tracer scans, using the single-tracer scans as gold standards. The multi-tracer data were processed for static (SUV) and kinetic (K(1), K(net)) endpoints for each tracer, followed by linear regression analysis of multi-tracer versus single-tracer results. Static and quantitative dynamic imaging measures of FDG were both accurately recovered from the multi-tracer scans, closely matching the single-tracer FDG standards (R > 0.99). Quantitative blood flow information, as measured by PTSM K(1) and SUV, was also accurately recovered from the multi-tracer scans (R = 0.97). Recovery of ATSM kinetic parameters proved more difficult, though the ATSM SUV was reasonably well recovered (R = 0.92). We conclude that certain additional information from one to two shorter-lived PET tracers may be measured in a rapid multi-tracer scan alongside FDG without compromising the assessment of glucose metabolism. Such additional and complementary information has the potential to improve tumor characterization in vivo, warranting further investigation of rapid multi-tracer techniques.

Estimating myocardial perfusion from dynamic contrast-enhanced CMR with a model-independent deconvolution method.

BACKGROUND: Model-independent analysis with B-spline regularization has been used to quantify myocardial blood flow (perfusion) in dynamic contrast-enhanced cardiovascular magnetic resonance (CMR) studies. However, the model-independent approach has not been extensively evaluated to determine how the contrast-to-noise ratio between blood and tissue enhancement affects estimates of myocardial perfusion and the degree to which the regularization is dependent on the noise in the measured enhancement data. We investigated these questions with a model-independent analysis method that uses iterative minimization and a temporal smoothness regularizer. Perfusion estimates using this method were compared to results from dynamic 13N-ammonia PET. RESULTS: An iterative model-independent analysis method was developed and tested to estimate regional and pixelwise myocardial perfusion in five normal subjects imaged with a saturation recovery turboFLASH sequence at 3 T CMR. Estimates of myocardial perfusion using model-independent analysis are dependent on the choice of the regularization weight parameter, which increases nonlinearly to handle large decreases in the contrast-to-noise ratio of the measured tissue enhancement data. Quantitative perfusion estimates in five subjects imaged with 3 T CMR were 1.1 +/- 0.8 ml/min/g at rest and 3.1 +/- 1.7 ml/min/g at adenosine stress. The perfusion estimates correlated with dynamic 13N-ammonia PET (y = 0.90x + 0.24, r = 0.85) and were similar to results from other validated CMR studies. CONCLUSION: This work shows that a model-independent analysis method that uses iterative minimization and temporal regularization can be used to quantify myocardial perfusion with dynamic contrast-enhanced perfusion CMR. Results from this method are robust to choices in the regularization weight parameter over relatively large ranges in the contrast-to-noise ratio of the tissue enhancement data.

Evaluation of rapid dual-tracer (62)Cu-PTSM + (62)Cu-ATSM PET in dogs with spontaneously occurring tumors.

We are developing methods for imaging multiple PET tracers in a single scan with staggered injections, where imaging measures for each tracer are separated and recovered using differences in tracer kinetics and radioactive decay. In this work, signal separation performance for rapid dual-tracer (62)Cu-PTSM (blood flow) + (62)Cu-ATSM (hypoxia) tumor imaging was evaluated in a large animal model. Four dogs with pre-existing tumors received a series of dynamic PET scans with (62)Cu-PTSM and (62)Cu-ATSM, permitting evaluation of a rapid dual-tracer protocol designed by previous simulation work. Several imaging measures were computed from the dual-tracer data and compared with those from separate, single-tracer imaging. Static imaging measures (e.g. SUV) for each tracer were accurately recovered from dual-tracer data. The wash-in (k(1)) and wash-out (k(2)) rate parameters for both tracers were likewise well recovered (r = 0.87-0.99), but k(3) was not accurately recovered for PTSM (r = 0.19) and moderately well recovered for ATSM (r = 0.70). Some degree of bias was noted, however, which may potentially be overcome through further refinement of the signal separation algorithms. This work demonstrates that complementary information regarding tumor blood flow and hypoxia can be acquired by a single dual-tracer PET scan, and also that the signal separation procedure works effectively for real physiologic data with realistic levels of kinetic model mismatch. Rapid multi-tracer PET has the potential to improve tumor assessment for image-guide therapy and monitoring, and further investigation with these and other tracers is warranted.

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.

Application of separable parameter space techniques to multi-tracer PET compartment modeling.

Multi-tracer positron emission tomography (PET) can image two or more tracers in a single scan, characterizing multiple aspects of biological functions to provide new insights into many diseases. The technique uses dynamic imaging, resulting in time-activity curves that contain contributions from each tracer present. The process of separating and recovering separate images and/or imaging measures for each tracer requires the application of kinetic constraints, which are most commonly applied by fitting parallel compartment models for all tracers. Such multi-tracer compartment modeling presents challenging nonlinear fits in multiple dimensions. This work extends separable parameter space kinetic modeling techniques, previously developed for fitting single-tracer compartment models, to fitting multi-tracer compartment models. The multi-tracer compartment model solution equations were reformulated to maximally separate the linear and nonlinear aspects of the fitting problem, and separable least-squares techniques were applied to effectively reduce the dimensionality of the nonlinear fit. The benefits of the approach are then explored through a number of illustrative examples, including characterization of separable parameter space multi-tracer objective functions and demonstration of exhaustive search fits which guarantee the true global minimum to within arbitrary search precision. Iterative gradient-descent algorithms using Levenberg-Marquardt were also tested, demonstrating improved fitting speed and robustness as compared to corresponding fits using conventional model formulations. The proposed technique overcomes many of the challenges in fitting simultaneous multi-tracer PET compartment models.

FDG and FLT-PET for Early measurement of response to 37.5 mg daily sunitinib therapy in metastatic renal cell carcinoma.

BACKGROUND: Metastatic renal cell carcinoma has a poor prognosis and an intrinsic resistance to standard treatment. Sunitinib is an oral receptor tyrosine kinase inhibitor that has been used as a first-line targeted therapy in metastatic renal cell carcinoma. While computed tomography (CT) is currently the gold standard for response assessment in oncological trials, numerous studies have shown that positron emission tomography (PET) imaging can provide information predictive of tumor response to treatment earlier than the typical interval for standard of care follow-up CT imaging. In this exploratory study we sought to characterize early tumor response in patients with metastatic renal cell carcinoma treated with continuous daily 37.5 mg sunitinib therapy. METHODS: Twenty patients underwent dynamic acquisition positron emission tomography (PET) imaging using (18) F-fluorodeoxyglucose (FDG) and (18) F-fluorothymidine (FLT) at baseline and early in treatment (after 1, 2, 3 or 4 weeks) with 37.5 mg continuous daily dosing of sunitinib. Semi-quantitative analyses were performed to characterize the tumor metabolic (FDG) and proliferative (FLT) responses to treatment. RESULTS: Proliferative responses were observed in 9/19 patients and occurred in 2 patients at one week (the earliest interval evaluated) after the initiation of therapy. A metabolic response was observed in 5/19 patients, however this was not observed until after two weeks of therapy were completed. Metabolic progression was observed in 2/19 patients and proliferative progression was observed in 1/19 patients. Baseline FDG-PET tumor maximum standardized uptake values correlated inversely with overall survival (p = 0.0036). Conversely, baseline (18) F-fluorothymidine PET imaging did not have prognostic value (p = 0.56) but showed a greater early response rate at 1-2 weeks after initiating therapy. CONCLUSIONS: While preliminary in nature, these results show an immediate and sustained proliferative response followed by a delayed metabolic response beginning after two weeks in metastatic renal cell carcinoma treated with a continuous daily dose of 37.5 mg sunitinib. The results provide evidence of tumor response to lower-dose sunitinib while also supporting the inclusion of PET imaging as a tool for early assessment in oncological clinical trials. TRIAL REGISTRATION: ID: NCT00694096.

Comparative evaluation of lesion detectability for 6 PET imaging platforms using a highly reproducible whole-body phantom with (22)Na lesions and localization ROC analysis.

UNLABELLED: The lesion detectability performance of 6 PET imaging platforms has been compared using a highly reproducible whole-body phantom and localization receiver operating characteristic (LROC) analysis. METHODS: A realistic whole-body phantom consisting of brain, thorax with lungs and liver, and pelvis with bladder was assembled and outfitted with 27 semipermanent (22)Na lesions of various sizes and activity concentrations. The background compartments were reproducibly filled with (18)F solutions. The phantom was imaged under the condition of equal emission scan time on 7 PET platforms: Advance, HR+, HR961, C-PET, IRIX, MCD, and AXIS. Imaging data were processed using manufacturer-supplied software and defaults, and LROC evaluation was performed using 11 human observers. RESULTS: Near-nominal counting rates were obtained for the NaI systems, and the bismuth germanate (BGO) systems were operated well below nominal counting rates. The BGO systems provided the highest lesion detection performance, followed by the large-area dedicated NaI system, and hybrid PET gamma cameras. Lesion detectability was highly dependent on lesion size, with all systems exhibiting similar performance for 16-mm lesions but differentiated performance for lesions < or =12 mm. CONCLUSION: Reconstruction methodology can have a significant effect on lesion detectability. PET lesion detectability performance is correlated with system cost and imaging characteristics. For a particular imaging task, care should be taken to ensure that the scanner being used is appropriate and that the scan time is adjusted accordingly to ensure good lesion detectability.

LOR-OSEM: statistical PET reconstruction from raw line-of-response histograms.

Iterative statistical reconstruction methods are becoming the standard in positron emission tomography (PET). Conventional maximum-likelihood expectation-maximization (MLEM) and ordered-subsets (OSEM) algorithms act on data which have been pre-processed into corrected, evenly-spaced histograms; however, such pre-processing corrupts the Poisson statistics. Recent advances have incorporated attenuation, scatter and randoms compensation into the iterative reconstruction. The objective of this work was to incorporate the remaining pre-processing steps, including arc correction, to reconstruct directly from raw unevenly-spaced line-of-response (LOR) histograms. This exactly preserves Poisson statistics and full spatial information in a manner closely related to listmode ML, making full use of the ML statistical model. The LOR-OSEM algorithm was implemented using a rotation-based projector which maps directly to the unevenly-spaced LOR grid. Simulation and phantom experiments were performed to characterize resolution, contrast and noise properties for 2D PET. LOR-OSEM provided a beneficial noise-resolution tradeoff, outperforming AW-OSEM by about the same margin that AW-OSEM outperformed pre-corrected OSEM. The relationship between LOR-ML and listmode ML algorithms was explored, and implementation differences are discussed. LOR-OSEM is a viable alternative to AW-OSEM for histogram-based reconstruction with improved spatial resolution and noise properties.

An SVD Investigation of Modeling Scatter in Multiple Energy Windows for Improved SPECT Images.

In this work singular value decomposition (SVD) techniques are used to investigate how the use of low energy photons and multiple energy windows affects the noise properties of Tc-99m SPECT imaging. We have previously shown that, when modeling scatter in the projector and backprojector of iterative reconstruction algorithms, simultaneous reconstruction from multiple energy window data can result in very different noise characteristics. Further, the properties depend upon the width and number of energy windows used. To investigate this further, we have generated photon transport matrices using models for scatter, an elliptical phantom containing cold rods of various sizes, and a number of multiple energy window acquisition schemes. Transfer matrices were also generated for the cases of perfect scatter rejection and ideal scatter subtraction. The matrices were decomposed using SVD, and signal power and projection space variance spectra were computed using the basis formed by the left singular vectors. Results indicate very different noise levels for the various energy window combinations. The perfect scatter rejection case resulted in the lowest variance spectrum, and reconstruction-based scatter compensation performed better than the scatter subtraction case. When including lower energy photons in reconstruction-based scatter compensation, using a series of multiple energy windows outperformed a single large energy window. One multiple window combination is presented which achieves a lower variance spectrum than the standard 20% energy window, indicating the potential for using low energy photons to improve the noise characteristics of SPECT images.

Comparison of Static and Dynamic Cardiac Perfusion Thallium-201 SPECT.

Cardiac SPECT is typically performed clinically with static imaging protocols and visually assessed for perfusion defects based upon the relative intensity of myocardial regions. Dynamic imaging, however, has the potential to provide quantitative measures of flow, possibly improving diagnosis. The objective of this study was to compare the information content of dynamic and static thallium SPECT imaging as measures of myocardial perfusion. Studies were performed in four canines, each with an occlusion placed on the left anterior descending coronary artery. Dynamic SPECT imaging was performed at rest and under adenosine stress, and subsets of the data were summed to provide corresponding static datasets for identical physiologic conditions. Microsphere-derived flow measurements were used as the gold standard. The dynamic data were fit to a two-compartment model to provide regional estimates of wash-in rate parameters. Occluded-to-normal ratios were also calculated for each canine study. The results show comparable correlations with microspheres for both wash-in and static scaled image intensities. The dynamic data provided higher defect contrasts, which were more accurate than the static occluded to normal ratios. Preliminary studies were also performed in two patients and the static and dynamic data compared. These results show that dynamic thallium imaging may provide improved diagnostic information compared to static imaging for myocardial perfusion SPECT studies.

Generalized separable parameter space techniques for fitting 1K-5K serial compartment models.

PURPOSE: Kinetic modeling is widely used to analyze dynamic imaging data, estimating kinetic parameters that quantify functional or physiologic processes in vivo. Typical kinetic models give rise to nonlinear solution equations in multiple dimensions, presenting a complex fitting environment. This work generalizes previously described separable nonlinear least-squares techniques for fitting serial compartment models with up to three tissue compartments and five rate parameters. METHODS: The approach maximally separates the linear and nonlinear aspects of the modeling equations, using a formulation modified from previous basis function methods to avoid a potential mathematical degeneracy. A fast and robust algorithm for solving the linear subproblem with full user-defined constraints is also presented. The generalized separable parameter space technique effectively reduces the dimensionality of the nonlinear fitting problem to one dimension for 2K-3K compartment models, and to two dimensions for 4K-5K models. RESULTS: Exhaustive search fits, which guarantee identification of the true global minimum fit, required approximately 10 ms for 2K-3K and 1.1 s for 4K-5K models, respectively. The technique is also amenable to fast gradient-descent iterative fitting algorithms, where the reduced dimensionality offers improved convergence properties. The objective function for the separable parameter space nonlinear subproblem was characterized and found to be generally well-behaved with a well-defined global minimum. Separable parameter space fits with the Levenberg-Marquardt algorithm required fewer iterations than comparable fits for conventional model formulations, averaging 1 and 7 ms for 2K-3K and 4K-5K models, respectively. Sensitivity to initial conditions was likewise reduced. CONCLUSIONS: The separable parameter space techniques described herein generalize previously described techniques to encompass 1K-5K compartment models, enable robust solution of the linear subproblem with full user-defined constraints, and are amenable to rapid and robust fitting using iterative gradient-descent type algorithms.

Effect of varying number of OSEM subsets on PET lesion detectability.

UNLABELLED: Iterative reconstruction has become the standard for routine clinical PET imaging. However, iterative reconstruction is computationally expensive, especially for time-of-flight (TOF) data. Block-iterative algorithms such as ordered-subsets expectation maximization (OSEM) are commonly used to accelerate the reconstruction. There is a tradeoff between the number of subsets and reconstructed image quality. The objective of this work was to evaluate the effect of varying the number of OSEM subsets on lesion detection for general oncologic PET imaging. METHODS: Experimental phantom data were taken from the Utah PET Lesion Detection Database, modeling whole-body oncologic (18)F-FDG PET imaging of a 92-kg patient. The experiment consisted of 24 scans over 4 d on a TOF PET/CT scanner, with up to 23 lesions (diameter, 6-16 mm) distributed throughout the thorax, abdomen, and pelvis. Images were reconstructed with maximum-likelihood expectation maximization (MLEM) and with OSEM using 2-84 subsets. The reconstructions were repeated both with and without TOF. Localization receiver-operating-characteristic (LROC) analysis was applied using the channelized nonprewhitened observer. The observer was first used to optimize the number of iterations and smoothing filter for each case that maximized lesion-detection performance for these data; this was done to ensure that fair comparisons were made with each test case operating near its optimal performance. The probability of correct localization and the area under the LROC curve were then analyzed as functions of the number of subsets to characterize the effect of OSEM on lesion-detection performance. RESULTS: Compared with the baseline MLEM algorithm, lesion-detection performance with OSEM declined as the number of subsets increased. The decline was moderate out to about 12-14 subsets and then became progressively steeper as the number of subsets increased. Comparing TOF with non-TOF results, the magnitude of the performance drop was larger for TOF reconstructions. CONCLUSION: PET lesion-detection performance is degraded when OSEM is used with a large number of subsets. This loss of image quality can be controlled using a moderate number of subsets (e.g., 12-14 or fewer), retaining a large degree of acceleration while maintaining high image quality. The use of more aggressive subsetting can result in image quality degradations that offset the benefits of using TOF or longer scan times.

Methodology for quantitative rapid multi-tracer PET tumor characterizations.

Positron emission tomography (PET) can image a wide variety of functional and physiological parameters in vivo using different radiotracers. As more is learned about the molecular basis for disease and treatment, the potential value of molecular imaging for characterizing and monitoring disease status has increased. Characterizing multiple aspects of tumor physiology by imaging multiple PET tracers in a single patient provides additional complementary information, and there is a significant body of literature supporting the potential value of multi-tracer PET imaging in oncology. However, imaging multiple PET tracers in a single patient presents a number of challenges. A number of techniques are under development for rapidly imaging multiple PET tracers in a single scan, where signal-recovery processing algorithms are employed to recover various imaging endpoints for each tracer. Dynamic imaging is generally used with tracer injections staggered in time, and kinetic constraints are utilized to estimate each tracers' contribution to the multi-tracer imaging signal. This article summarizes past and ongoing work in multi-tracer PET tumor imaging, and then organizes and describes the main algorithmic approaches for achieving multi-tracer PET signal-recovery. While significant advances have been made, the complexity of the approach necessitates protocol design, optimization, and testing for each particular tracer combination and application. Rapid multi-tracer PET techniques have great potential for both research and clinical cancer imaging applications, and continued research in this area is warranted.

Single-scan dual-tracer FLT+FDG PET tumor characterization.

Rapid multi-tracer PET aims to image two or more tracers in a single scan, simultaneously characterizing multiple aspects of physiology and function without the need for repeat imaging visits. Using dynamic imaging with staggered injections, constraints on the kinetic behavior of each tracer are applied to recover individual-tracer measures from the multi-tracer PET signal. The ability to rapidly and reliably image both (18)F-fluorodeoxyglucose (FDG) and (18)F-fluorothymidine (FLT) would provide complementary measures of tumor metabolism and proliferative activity, with important applications in guiding oncologic treatment decisions and assessing response. However, this tracer combination presents one of the most challenging dual-tracer signal-separation problems--both tracers have the same radioactive half-life, and the injection delay is short relative to the half-life and tracer kinetics. This work investigates techniques for single-scan dual-tracer FLT+FDG PET tumor imaging, characterizing the performance of recovering static and dynamic imaging measures for each tracer from dual-tracer datasets. Simulation studies were performed to characterize dual-tracer signal-separation performance for imaging protocols with both injection orders and injection delays of 10-60 min. Better performance was observed when FLT was administered first, and longer delays before administration of FDG provided more robust signal-separation and recovery of the single-tracer imaging measures. An injection delay of 30 min led to good recovery (R > 0.96) of static image values (e.g. SUV), K(net), and K(1) as compared to values from separate, single-tracer time-activity curves. Recovery of higher order rate parameters (k(2), k(3)) was less robust, indicating that information regarding these parameters was harder to recover in the presence of statistical noise and dual-tracer effects. Performance of the dual-tracer FLT(0 min)+FDG(32 min) technique was further evaluated using PET/CT imaging studies in five patients with primary brain tumors where the data from separate scans of each tracer were combined to synthesize dual-tracer scans with known single-tracer components; results demonstrated similar dual-tracer signal recovery performance. We conclude that rapid dual-tracer FLT+FDG tumor imaging is feasible and can provide quantitative tumor imaging measures comparable to those from conventional separate-scan imaging.