Joint EANM/SNMMI guideline on radiomics in nuclear medicine : Jointly supported by the EANM Physics Committee and the SNMMI Physics, Instrumentation and Data Sciences Council.

PURPOSE: The purpose of this guideline is to provide comprehensive information on best practices for robust radiomics analyses for both hand-crafted and deep learning-based approaches. METHODS: In a cooperative effort between the EANM and SNMMI, we agreed upon current best practices and recommendations for relevant aspects of radiomics analyses, including study design, quality assurance, data collection, impact of acquisition and reconstruction, detection and segmentation, feature standardization and implementation, as well as appropriate modelling schemes, model evaluation, and interpretation. We also offer an outlook for future perspectives. CONCLUSION: Radiomics is a very quickly evolving field of research. The present guideline focused on established findings as well as recommendations based on the state of the art. Though this guideline recognizes both hand-crafted and deep learning-based radiomics approaches, it primarily focuses on the former as this field is more mature. This guideline will be updated once more studies and results have contributed to improved consensus regarding the application of deep learning methods for radiomics. Although methodological recommendations in the present document are valid for most medical image modalities, we focus here on nuclear medicine, and specific recommendations when necessary are made for PET/CT, PET/MR, and quantitative SPECT.

Instrumentation factors affecting variance and bias of quantifying tracer uptake with PET/CT.

PURPOSE: The variances and biases inherent in quantifying PET tracer uptake from instrumentation factors are needed to ascertain the significance of any measured differences such as in quantifying response to therapy. The authors studied the repeatability and reproducibility of serial PET measures of activity as a function of object size, acquisition, reconstruction, and analysis method on one scanner and at three PET centers using a single protocol with long half-life phantoms. METHODS: The authors assessed standard deviations (SDs) and mean biases of consecutive measures of PET activity concentrations in a uniform phantom and a NEMA NU-2 image quality (IQ) phantom filled with 9 months half-life 68Ge in an epoxy matrix. Activity measurements were normalized by dividing by a common decay corrected true value and reported as recovery coefficients (RCs). Each experimental set consisted of 20 consecutive PET scans of either a stationary phantom to evaluate repeatability or a repositioned phantom to assess reproducibility. One site conducted a comprehensive series of repeatability and reproducibility experiments, while two other sites repeated the reproducibility experiments using the same IQ phantom. An equation was derived to estimate the SD of a new PET measure from a known SD based on the ratios of available coincident counts between the two PET measures. RESULTS: For stationary uniform phantom scans, the SDs of maximum RCs were three to five times less than predicted for uncorrelated pixels within circular regions of interest (ROIs) with diameters ranging from 1 to 15 cm. For stationary IQ phantom scans from 1 cm diameter ROIs, the average SDs of mean and maximum RCs ranged from 1.4% to 8.0%, depending on the methods of acquisition and reconstruction (coefficients of variation range 2.5% to 9.8%). Similar SDs were observed for both analytic and iterative reconstruction methods (p > or = 0.08). SDs of RCs for 2D acquisitions were significantly higher than for 3D acquisitions (p < or =s 0.008) for same acquisition and processing parameters. SDs of maximum RCs were larger than corresponding mean values for stationary IQ phantom scans ( < or = 0.02), although the magnitude of difference is reduced due to noise correlations in the image. Increased smoothing decreased SDs ( < or =s 0.045) and decreased maximum and mean RCs (p < or = 0.02). Reproducibility of GE DSTE, Philips Gemini TF, and Siemens Biograph Hi-REZ PET/CT scans of the same IQ phantom, with similar acquisition, reconstruction, and repositioning among 20 scans, were, in general, similar (mean and maximum RC SD range 2.5% to 4.8%). CONCLUSIONS: Short-term scanner variability is low compared to other sources of error. There are tradeoffs in noise and bias depending on acquisition, processing, and analysis methods. The SD of a new PET measure can be estimated from a known SD if the ratios of available coincident counts between the two PET scanner acquisitions are known and both employ the same ROI definition. Results suggest it is feasible to use PET/CTs from different vendors and sites in clinical trials if they are properly cross-calibrated.

Standardization in Quantitative Imaging: A Multicenter Comparison of Radiomic Features from Different Software Packages on Digital Reference Objects and Patient Data Sets.

Radiomic features are being increasingly studied for clinical applications. We aimed to assess the agreement among radiomic features when computed by several groups by using different software packages under very tightly controlled conditions, which included standardized feature definitions and common image data sets. Ten sites (9 from the NCI's Quantitative Imaging Network] positron emission tomography-computed tomography working group plus one site from outside that group) participated in this project. Nine common quantitative imaging features were selected for comparison including features that describe morphology, intensity, shape, and texture. The common image data sets were: three 3D digital reference objects (DROs) and 10 patient image scans from the Lung Image Database Consortium data set using a specific lesion in each scan. Each object (DRO or lesion) was accompanied by an already-defined volume of interest, from which the features were calculated. Feature values for each object (DRO or lesion) were reported. The coefficient of variation (CV), expressed as a percentage, was calculated across software packages for each feature on each object. Thirteen sets of results were obtained for the DROs and patient data sets. Five of the 9 features showed excellent agreement with CV < 1%; 1 feature had moderate agreement (CV < 10%), and 3 features had larger variations (CV ≥ 10%) even after attempts at harmonization of feature calculations. This work highlights the value of feature definition standardization as well as the need to further clarify definitions for some features.

Measured count-rate performance of the Discovery STE PET/CT scanner in 2D, 3D and partial collimation acquisition modes.

We measured count rates and scatter fraction on the Discovery STE PET/CT scanner in conventional 2D and 3D acquisition modes, and in a partial collimation mode between 2D and 3D. As part of the evaluation of using partial collimation, we estimated global count rates using a scanner model that combined computer simulations with an empirical live-time function. Our measurements followed the NEMA NU2 count rate and scatter-fraction protocol to obtain true, scattered and random coincidence events, from which noise equivalent count (NEC) rates were calculated. The effect of patient size was considered by using 27 cm and 35 cm diameter phantoms, in addition to the standard 20 cm diameter cylindrical count-rate phantom. Using the scanner model, we evaluated two partial collimation cases: removing half of the septa (2.5D) and removing two-thirds of the septa (2.7D). Based on predictions of the model, a 2.7D collimator was constructed. Count rates and scatter fractions were then measured in 2D, 2.7D and 3D. The scanner model predicted relative NEC variation with activity, as confirmed by measurements. The measured 2.7D NEC was equal or greater than 3D NEC for all activity levels in the 27 cm and 35 cm phantoms. In the 20 cm phantom, 3D NEC was somewhat higher ( approximately 15%) than 2.7D NEC at 100 MBq. For all higher activity concentrations, 2.7D NEC was greater and peaked 26% above the 3D peak NEC. The peak NEC in 2.7D mode occurred at approximately 425 MBq, and was 26-50% greater than the peak 3D NEC, depending on object size. NEC in 2D was considerably lower, except at relatively high activity concentrations. Partial collimation shows promise for improved noise equivalent count rates in clinical imaging without altering other detector parameters.

Multisite Concordance of DSC-MRI Analysis for Brain Tumors: Results of a National Cancer Institute Quantitative Imaging Network Collaborative Project.

BACKGROUND AND PURPOSE: Standard assessment criteria for brain tumors that only include anatomic imaging continue to be insufficient. While numerous studies have demonstrated the value of DSC-MR imaging perfusion metrics for this purpose, they have not been incorporated due to a lack of confidence in the consistency of DSC-MR imaging metrics across sites and platforms. This study addresses this limitation with a comparison of multisite/multiplatform analyses of shared DSC-MR imaging datasets of patients with brain tumors. MATERIALS AND METHODS: DSC-MR imaging data were collected after a preload and during a bolus injection of gadolinium contrast agent using a gradient recalled-echo-EPI sequence (TE/TR = 30/1200 ms; flip angle = 72°). Forty-nine low-grade (n = 13) and high-grade (n = 36) glioma datasets were uploaded to The Cancer Imaging Archive. Datasets included a predetermined arterial input function, enhancing tumor ROIs, and ROIs necessary to create normalized relative CBV and CBF maps. Seven sites computed 20 different perfusion metrics. Pair-wise agreement among sites was assessed with the Lin concordance correlation coefficient. Distinction of low- from high-grade tumors was evaluated with the Wilcoxon rank sum test followed by receiver operating characteristic analysis to identify the optimal thresholds based on sensitivity and specificity. RESULTS: For normalized relative CBV and normalized CBF, 93% and 94% of entries showed good or excellent cross-site agreement (0.8 ≤ Lin concordance correlation coefficient ≤ 1.0). All metrics could distinguish low- from high-grade tumors. Optimum thresholds were determined for pooled data (normalized relative CBV = 1.4, sensitivity/specificity = 90%:77%; normalized CBF = 1.58, sensitivity/specificity = 86%:77%). CONCLUSIONS: By means of DSC-MR imaging data obtained after a preload of contrast agent, substantial consistency resulted across sites for brain tumor perfusion metrics with a common threshold discoverable for distinguishing low- from high-grade tumors.

Evaluation of Multiclass Model Observers in PET LROC Studies.

A localization ROC (LROC) study was conducted to evaluate nonprewhitening matched-filter (NPW) and channelized NPW (CNPW) versions of a multiclass model observer as predictors of human tumor-detection performance with PET images. Target localization is explicitly performed by these model observers. Tumors were placed in the liver, lungs, and background soft tissue of a mathematical phantom, and the data simulation modeled a full-3D acquisition mode. Reconstructions were performed with the FORE+AWOSEM algorithm. The LROC study measured observer performance with 2D images consisting of either coronal, sagittal, or transverse views of the same set of cases. Versions of the CNPW observer based on two previously published difference-of-Gaussian channel models demonstrated good quantitative agreement with human observers. One interpretation of these results treats the CNPW observer as a channelized Hotelling observer with implicit internal noise.

Evaluation of event position reconstruction in monolithic crystals that are optically coupled.

A PET detector featuring a pseudo-monolithic crystal is being developed as a more cost-effective alternative to a full monolithic crystal PET detector. This work evaluates different methods to localize the scintillation events in quartered monolithic crystals that are optically coupled. A semi-monolithic crystal assembly was formed using four 26 × 26 × 10 mm3 LYSO crystals optically coupled together using optical adhesive, to mimic a 52 × 52 × 10 mm3 monolithic crystal detector. The crystal assembly was coupled to a 64-channel multi-anode photomultiplier tube using silicon grease. The detector was calibrated using a 34 × 34 scan grid. Events were first filtered and depth separated using a multi-Lorentzian fit to the collected light distribution. Next, three different techniques were explored to generate the look up tables for the event positioning. The first technique was 'standard interpolation' across the interface. The second technique was 'central extrapolation', where a bin was placed at the midpoint of the interface and events positioned within the interface region were discarded. The third technique used a 'central overlap' method where an extended region was extrapolated at each interface. Events were then positioned using least-squares minimization and maximum likelihood methods. The least-squares minimization applied to the look up table generated with the standard interpolation technique had the best full width at half maximum (FWHM) intrinsic spatial resolution and the lowest bias. However, there were discontinuities in the event positioning that would most likely lead to artifacts in the reconstructed image. The central extrapolation technique also had discontinuities and a 30% sensitivity loss near the crystal-crystal interfaces. The central overlap technique had slightly degraded performance metrics, but it still provided ~2.1 mm intrinsic spatial resolution at the crystal-crystal interface and had a symmetric and continuously varying response function. Results using maximum likelihood positioning were similar to least-squares minimization for the central overlap data.

Imaging and dosimetric errors in 4D PET/CT-guided radiotherapy from patient-specific respiratory patterns: a dynamic motion phantom end-to-end study.

Effective positron emission tomography / computed tomography (PET/CT) guidance in radiotherapy of lung cancer requires estimation and mitigation of errors due to respiratory motion. An end-to-end workflow was developed to measure patient-specific motion-induced uncertainties in imaging, treatment planning, and radiation delivery with respiratory motion phantoms and dosimeters. A custom torso phantom with inserts mimicking normal lung tissue and lung lesion was filled with [(18)F]FDG. The lung lesion insert was driven by six different patient-specific respiratory patterns or kept stationary. PET/CT images were acquired under motionless ground truth, tidal breathing motion-averaged (3D), and respiratory phase-correlated (4D) conditions. Target volumes were estimated by standardized uptake value (SUV) thresholds that accurately defined the ground-truth lesion volume. Non-uniform dose-painting plans using volumetrically modulated arc therapy were optimized for fixed normal lung and spinal cord objectives and variable PET-based target objectives. Resulting plans were delivered to a cylindrical diode array at rest, in motion on a platform driven by the same respiratory patterns (3D), or motion-compensated by a robotic couch with an infrared camera tracking system (4D). Errors were estimated relative to the static ground truth condition for mean target-to-background (T/Bmean) ratios, target volumes, planned equivalent uniform target doses, and 2%-2 mm gamma delivery passing rates. Relative to motionless ground truth conditions, PET/CT imaging errors were on the order of 10-20%, treatment planning errors were 5-10%, and treatment delivery errors were 5-30% without motion compensation. Errors from residual motion following compensation methods were reduced to 5-10% in PET/CT imaging, <5% in treatment planning, and <2% in treatment delivery. We have demonstrated that estimation of respiratory motion uncertainty and its propagation from PET/CT imaging to RT planning, and RT delivery under a dose painting paradigm is feasible within an integrated respiratory motion phantom workflow. For a limited set of cases, the magnitude of errors was comparable during PET/CT imaging and treatment delivery without motion compensation. Errors were moderately mitigated during PET/CT imaging and significantly mitigated during RT delivery with motion compensation. This dynamic motion phantom end-to-end workflow provides a method for quality assurance of 4D PET/CT-guided radiotherapy, including evaluation of respiratory motion compensation methods during imaging and treatment delivery.

Amphetamine-induced dopamine release in human ventral striatum correlates with euphoria.

BACKGROUND: Studies in experimental animals have implicated the mesolimbic dopaminergic projections into the ventral striatum in the neural processes underlying behavioral reinforcement and motivated behavior; however, understanding the relationship between subjective emotional experience and ventral striatal dopamine (DA) release has awaited human studies. Using positron emission tomography (PET), we correlated the change in endogenous dopamine concentrations following dextroamphetamine (AMPH) administration with the associated hedonic response in human subjects and compared the strength of this correlation across striatal subregions. METHODS: We obtained PET measures of [(11)C]raclopride specific binding to DA D2/D3 receptors before and after AMPH injection (0.3 mg/kg IV) in seven healthy subjects. The change in [(11)C]raclopride binding potential (DeltaBP) induced by AMPH pretreatment and the correlation between DeltaBP and the euphoric response to AMPH were compared between the anteroventral striatum (AVS; comprised of accumbens area, ventromedial caudate, and anteroventral putamen) and the dorsal caudate (DCA) using an MRI-based region of interest analysis of the PET data. RESULTS: The mean DeltaBP was greater in the AVS than in the DCA (p <.05). The AMPH-induced changes in euphoria analog scale scores correlated inversely with DeltaBP in the AVS (r = -.95; p <.001), but not in the DCA (r =.30, ns). Post hoc assessments showed that changes in tension-anxiety ratings correlated positively with DeltaBP in the AVS (r =.80; p [uncorrected] <.05) and that similar relationships may exist between DeltaBP and emotion ratings in the ventral putamen (as were found in the AVS). CONCLUSIONS: The preferential sensitivity of the ventral striatum to the DA releasing effects of AMPH previously demonstrated in experimental animals extends to humans. The magnitude of ventral striatal DA release correlates positively with the hedonic response to AMPH.

PET measures of amphetamine-induced dopamine release in ventral versus dorsal striatum.

Regional differences in dextroamphetamine (AMPH)-induced dopamine (DA) release in the baboon striatum were assessed using positron emission tomographic (PET) measures of [11C]raclopride specific binding to DA D2/D3 receptors acquired before and after AMPH administration. The magnitude of the reduction in [11C]raclopride binding, following AMPH administration, was two-fold greater in the anteroventral striatum (comprised of ventral caudate, anteroventral putamen, and nucleus accumbens) than the dorsal striatum (dorsal caudate). A simulation study demonstrated that any potential biases due to resolution (partial volume) and alignment effects were significantly smaller than the magnitude of the observed results. These regional differences in the sensitivity of AMPH are compatible with microdialysis evidence in rats indicating that the magnitude of DA release in response to AMPH concentrations in the range tested is greater in ventral than dorsal striatal regions. Post hoc tests involving measures in other striatal regions showed that the baseline DA D2/D3 binding was highest and the correlation between AMPH dose and change in [11C]raclopride binding most significant in the putamen.

Three-dimensional imaging characteristics of the HEAD PENN-PET scanner.

UNLABELLED: A volume-imaging PET scanner, without interplane septa, for brain imaging has been designed and built to achieve high performance, specifically in spatial resolution and sensitivity. The scanner is unique in its use of a single annular crystal of Nal(Tl), which allows a field of view (FOV) of 25.6 cm in both the transverse and axial directions. Data are reconstructed into an image matrix of 128(3) with (2 mm)3 voxels, using three-dimensional image reconstruction algorithms. METHODS: Point-source measurements are performed to determine spatial resolution over the scanner FOV, and cylindrical phantom distributions are used to determine the sensitivity, scatter fraction and counting rate performance of the system. A three-dimensional brain phantom and 18F-FDG patient studies are used to evaluate image quality with three-dimensional reconstruction algorithms. RESULTS: The system spatial resolution is measured to be 3.5 mm in both the transverse and axial directions, in the center of the FOV. The true sensitivity, using the standard NEMA phantom (6 liter), is 660 kcps/microCi/ml, after subtracting a scatter fraction of 34%. Due to deadtime effects, we measure a peak true counting rate, after scatter and randoms subtraction, of 100 kcps at 0.7 mCi for a smaller brain-sized (1.1 liter) phantom, and 70 kcps for a head-sized (2.5 liter) phantom at the same activity. A typical 18F-FDG clinical brain study requires only 2 mCi to achieve high statistics (100 million true events) with a scan time of 30 min. CONCLUSION: The HEAD PENN-PET scanner is based on a cost-effective design using Nal(Tl) and has been shown to achieve high performance for brain studies and pediatric whole-body studies. As a full-time three-dimensional imaging scanner with a very large axial acceptance angle, high sensitivity is achieved. The system becomes counting-rate limited as the activity is increased, but we achieve high image quality with a small injected dose. This is a significant advantage for clinical imaging, particularly for pediatric patients.

Comparative evaluation of MR-based partial-volume correction schemes for PET.

UNLABELLED: Because of limitations of spatial resolution, quantitative PET measurements of cerebral blood flow, glucose metabolism and neuroreceptor binding are influenced by partial-volume averaging among neighboring tissues with differing tracer concentrations. METHODS: Two MR-based approaches to partial-volume correction of PET images were compared using simulations and a multicompartment phantom. The two-compartment method corrects PET data for the diluting effects of cerebrospinal fluid (CSF) spaces. The more complex three-compartment method also accounts for the effect of partial-volume averaging between gray and white matter. The effects of the most significant sources of error on MR-based partial-volume correction, including misregistration, resolution mismatch, segmentation errors and white matter heterogeneity, were evaluated. We also examined the relative usefulness of both approaches in PET studies of aging and neurodegenerative disease. RESULTS: Although the three-compartment method was highly accurate (with 100% gray matter recovery achieved in simulations), it was also more sensitive to all errors tested, particularly image segmentation and PET-MR registration. CONCLUSION: Based on these data, we conclude that the two-compartment approach is better suited for comparative PET studies, whereas the three-compartment algorithm is capable of greater accuracy for absolute quantitative measures.

A combined PET/CT scanner for clinical oncology.

UNLABELLED: The availability of accurately aligned, whole-body anatomical (CT) and functional (PET) images could have a significant impact on diagnosing and staging malignant disease and on identifying and localizing metastases. Computer algorithms to align CT and PET images acquired on different scanners are generally successful for the brain, whereas image alignment in other regions of the body is more problematic. METHODS: A combined PET/CT tomograph with the unique capability of acquiring accurately aligned functional and anatomical images for any part of the human body has been designed and built. The PET/CT scanner was developed as a combination of a Siemens Somatom AR.SP spiral CT and a partial-ring, rotating ECAT ART PET scanner. All components are mounted on a common rotational support within a single gantry. The PET and CT components can be operated either separately, or in combined mode. In combined mode, the CT images are used to correct the PET data for scatter and attenuation. Fully quantitative whole-body images are obtained for an axial extent of 100 cm in an imaging time of less than 1 h. When operated in PET mode alone, transmission scans are acquired with dual 137Cs sources. RESULTS: The scanner is fully operational and the combined device has been operated successfully in a clinical environment. Over 110 patients have been imaged, covering a range of different cancers, including lung, esophageal, head and neck, melanoma, lymphoma, pancreas, and renal cell. The aligned PET and CT images are used both for diagnosing and staging disease and for evaluating response to therapy. We report the first performance measurements from the scanner and present some illustrative clinical studies acquired in cancer patients. CONCLUSION: A combined PET and CT scanner is a practical and effective approach to acquiring co-registered anatomical and functional images in a single scanning session.

Attenuation correction for a combined 3D PET/CT scanner.

In this work we demonstrate the proof of principle of CT-based attenuation correction of 3D positron emission tomography (PET) data by using scans of bone and soft tissue equivalent phantoms and scans of humans. This method of attenuation correction is intended for use in a single scanner that combines volume-imaging (3D) PET with x-ray computed tomography (CT) for the purpose of providing accurately registered anatomical localization of structures seen in the PET image. The goal of this work is to determine if we can perform attenuation correction of the PET emission data using accurately aligned CT attenuation information. We discuss possible methods of calculating the PET attenuation map at 511 keV based on CT transmission information acquired from 40 keV through 140 keV. Data were acquired on separate CT and PET scanners and were aligned using standard image registration procedures. Results are presented on three of the attenuation calculation methods: segmentation, scaling, and our proposed hybrid segmentation/scaling method. The results are compared with those using the standard 3D PET attenuation correction method as a gold standard. We demonstrate the efficacy of our proposed hybrid method for converting the CT attenuation map from an effective CT photon energy of 70 keV to the PET photon energy of 511 keV. We conclude that using CT information is a feasible way to obtain attenuation correction factors for 3D PET.

[11C]WAY 100635: a radioligand for imaging 5-HT1A receptors with positron emission tomography.

The potent and selective 5-HT1A antagonist WAY 100635 (N-[2-]4-(2-methoxyphenyl)-1-piperazinyl]ethyl]-N-(2- pyridinyl)cyclohexanecarboxamide) was radiolabeled with 11C in high specific activity, and the in vivo properties of this radioligand were assessed in the brains of rats and monkeys. Following i.v. tail vein injection in rats, [11C]WAY 100635 rapidly penetrated into brain tissue and was retained over a 30-90 min time period in a manner consistent with the known distribution of 5-HT1A receptors. Pretreatment of rats with the selective 5-HT1A agonist (+/-)-8-OH-DPAT effectively blocked the retention of radioactivity in brain regions known to contain high densities of 5-HT1A receptors. The hippocampus-to-cerebellum radioactivity concentration ratio reached a maximum of 16:1 at 60 min post injection. Following i.v. injection of [11C]WAY 100635 in rhesus monkeys, the concentrations of radioactivity in brain regions were consistent with the reported distribution of 5-HT1A receptors in primates, and the frontal cortex-to-cerebellum ratio reached 5.5:1 at 80 min post injection. Pretreatment of the monkeys with (+/-)-8-OH-DPAT reduced this ratio to 1.4:1, and injection of (+/-)-8-OH-DPAT 20 min after the injection of [11C]WAY 100635 significantly displaced frontal cortex binding. The in vivo properties of [11C]WAY 100635 in rats and monkeys strongly support the future utility of this radioligand for imaging 5-HT1A receptors using positron emission tomography (PET).

The Theory of Three-Dimensional Image Reconstruction for PET.

A review is made of selected recent publications on three-dimensional image reconstruction for PET. A new algorithm is proposed which is designed to fully utilize all emitted gamma rays which can be detected in a truncated spherical detector. By such full utilization of emitted rays the new algorithm should produce images of better statistical accuracy than could be produced by previously known algorithms.

Positron emission tomography with a large axial acceptance angle: signal-to-noise considerations.

A volume imaging positron emission tomography (PET) scanner with a large acceptance angle, such as the PENN-PET, offers fine spatial sampling and resolution in three dimensions, and a high sensitivity because of the inclusion of all cross-plane rays. The signal-to-noise ratio (SNR) is used to evaluate image quality for different scanning conditions of the PENN-PET using an activated cylindrical phantom with cold spheres of various sizes. Raising the energy threshold to 400 keV improves the SNR by lowering the scatter fraction, though it also reduces the sensitivity. Increasing the axial acceptance angle from +/-1.3 degrees to +/-6.5 degrees improves the SNR by increasing the sensitivity, even with a two-dimensional reconstruction algorithm, which compromises spatial resolution in the axial direction for points at the edge of the radial field of view. Initial results show that a three-dimensional reconstruction offers an improved SNR over a two-dimensional reconstruction that does not use all cross-plane rays.

Effect of increased axial field of view on the performance of a volume PET scanner.

The performance of the PENN-PET 240H scanner from UGM Medical Systems is tested and compared to the prototype PENN-PET scanner built at the University of Pennsylvania. The UGM PENN-PET scanner consists of six continuous position-sensitive NaI(Tl) detectors, which results in a 50 cm transverse field-of-view and a 12.8 cm axial field-of-view. The fine spatial sampling in the axial direction allows the data to be sorted into as many as 64 transverse planes, each 2 mm thick. A large axial acceptance angle, without interplane septa, results in a high-sensitivity and low-randoms fraction, with a low-scatter fraction due to the use of a narrow photopeak energy window. This work emphasizes those performance measurements that illustrate the special characteristics of a volume imaging scanner and how they change as the axial length is increased.

Exact and approximate rebinning algorithms for 3-D PET data.

This paper presents two new rebinning algorithms for the reconstruction of three-dimensional (3-D) positron emission tomography (PET) data. A rebinning algorithm is one that first sorts the 3-D data into an ordinary two-dimensional (2-D) data set containing one sinogram for each transaxial slice to be reconstructed; the 3-D image is then recovered by applying to each slice a 2-D reconstruction method such as filtered-backprojection. This approach allows a significant speedup of 3-D reconstruction, which is particularly useful for applications involving dynamic acquisitions or whole-body imaging. The first new algorithm is obtained by discretizing an exact analytical inversion formula. The second algorithm, called the Fourier rebinning algorithm (FORE), is approximate but allows an efficient implementation based on taking 2-D Fourier transforms of the data. This second algorithm was implemented and applied to data acquired with the new generation of PET systems and also to simulated data for a scanner with an 18 degrees axial aperture. The reconstructed images were compared to those obtained with the 3-D reprojection algorithm (3DRP) which is the standard "exact" 3-D filtered-backprojection method. Results demonstrate that FORE provides a reliable alternative to 3DRP, while at the same time achieving an order of magnitude reduction in processing time.

Fast fully 3-D image reconstruction in PET using planograms.

We present a method of performing fast and accurate three-dimensional (3-D) backprojection using only Fourier transform operations for line-integral data acquired by planar detector arrays in positron emission tomography. This approach is a 3-D extension of the two-dimensional (2-D) linogram technique of Edholm. By using a special choice of parameters to index a line of response (LOR) for a pair of planar detectors, rather than the conventional parameters used to index a LOR for a circular tomograph, all the LORs passing through a point in the field of view (FOV) lie on a 2-D plane in the four-dimensional (4-D) data space. Thus, backprojection of all the LORs passing through a point in the FOV corresponds to integration of a 2-D plane through the 4-D "planogram." The key step is that the integration along a set of parallel 2-D planes through the planogram, that is, backprojection of a plane of points, can be replaced by a 2-D section through the origin of the 4-D Fourier transform of the data. Backprojection can be performed as a sequence of Fourier transform operations, for faster implementation. In addition, we derive the central-section theorem for planogram format data, and also derive a reconstruction filter for both backprojection-filtering and filtered-backprojection reconstruction algorithms. With software-based Fourier transform calculations we provide preliminary comparisons of planogram backprojection to standard 3-D backprojection and demonstrate a reduction in computation time by a factor of approximately 15.