Improved Positron Emission Tomography Quantification: Evaluation of a Maximum-Likelihood Scatter Scaling Algorithm.

Incorrect scatter scaling of positron emission tomography (PET) images can lead to halo artifacts, quantitative bias, or reconstruction failure. Tail-fitted scatter scaling (TFSS) possesses performance limitations in multiple cases. This study aims to investigate a novel method for scatter scaling: maximum-likelihood scatter scaling (MLSS) in scenarios where TFSS tends to induce artifacts or are observed to cause reconstruction abortion. [68Ga]Ga-RGD PET scans of nine patients were included in cohort 1 in the scope of investigating the reduction of halo artifacts relative to the scatter estimation method. PET scans of 30 patients administrated with [68Ga]Ga-uPAR were included in cohort 2, used for an evaluation of the robustness of MLSS in cases where TFSS-integrated reconstructions are observed to fail. A visual inspection of MLSS-corrected images scored higher than TFSS-corrected reconstructions of cohort 1. The quantitative investigation near the bladder showed a relative difference in tracer uptake of up to 94.7%. A reconstruction of scans included in cohort 2 resulted in failure in 23 cases when TFSS was used. The lesion uptake values of cohort 2 showed no significant difference. MLSS is suggested as an alternative scatter-scaling method relative to TFSS with the aim of reducing halo artifacts and a robust reconstruction process.

Population-based deep image prior for dynamic PET denoising: A data-driven approach to improve parametric quantification.

The high noise level of dynamic Positron Emission Tomography (PET) images degrades the quality of parametric images. In this study, we aim to improve the quality and quantitative accuracy of Ki images by utilizing deep learning techniques to reduce the noise in dynamic PET images. We propose a novel denoising technique, Population-based Deep Image Prior (PDIP), which integrates population-based prior information into the optimization process of Deep Image Prior (DIP). Specifically, the population-based prior image is generated from a supervised denoising model that is trained on a prompts-matched static PET dataset comprising 100 clinical studies. The 3D U-Net architecture is employed for both the supervised model and the following DIP optimization process. We evaluated the efficacy of PDIP for noise reduction in 25%-count and 100%-count dynamic PET images from 23 patients by comparing with two other baseline techniques: the Prompts-matched Supervised model (PS) and a conditional DIP (CDIP) model that employs the mean static PET image as the prior. Both the PS and CDIP models show effective noise reduction but result in smoothing and removal of small lesions. In addition, the utilization of a single static image as the prior in the CDIP model also introduces a similar tracer distribution to the denoised dynamic frames, leading to lower Ki in general as well as incorrect Ki in the descending aorta. By contrast, as the proposed PDIP model utilizes intrinsic image features from the dynamic dataset and a large clinical static dataset, it not only achieves comparable noise reduction as the supervised and CDIP models but also improves lesion Ki predictions.

Positron range in combination with point-spread-function correction: an evaluation of different implementations for [124I]-PET imaging.

AIM: To evaluate the effect of combining positron range correction (PRC) with point-spread-function (PSF) correction and to compare different methods of implementation into iterative image reconstruction for 124I-PET imaging. MATERIALS AND METHODS: Uniform PR blurring kernels of 124I were generated using the GATE (GEANT4) framework in various material environments (lung, water, and bone) and matched to a 3D matrix. The kernels size was set to 11 × 11 × 11 based on the maximum PR in water and the voxel size of the PET system. PET image reconstruction was performed using the standard OSEM algorithm, OSEM with PRC implemented before the forward projection (OSEM+PRC simplified) and OSEM with PRC implemented in both forward- and back-projection steps (full implementation) (OSEM+PRC). Reconstructions were repeated with resolution recovery, point-spread function (PSF) included. The methods and kernel variation were validated using different phantoms filled with 124I acquired on a Siemens mCT PET/CT system. The data was evaluated for contrast recovery and image noise. RESULTS: Contrast recovery improved by 2-10% and 4-37% with OSEM+PRC simplified and OSEM+PRC, respectively, depending on the sphere size of the NEMA IQ phantom. Including PSF in the reconstructions further improved contrast by 4-19% and 3-16% with the PSF+PRC simplified and PSF+PRC, respectively. The benefit of PRC was more pronounced within low-density material. OSEM-PRC and OSEM-PSF as well as OSEM-PSF+PRC in its full- and simplified implementation showed comparable noise and convergence. OSEM-PRC simplified showed comparably faster convergence but at the cost of increased image noise. CONCLUSIONS: The combination of the PSF and PRC leads to increased contrast recovery with reduced image noise compared to stand-alone PSF or PRC reconstruction. For OSEM-PRC reconstructions, a full implementation in the reconstruction is necessary to handle image noise. For the combination of PRC with PSF, a simplified PRC implementation can be used to reduce reconstruction times.

Quantitative evaluation of a deep learning-based framework to generate whole-body attenuation maps using LSO background radiation in long axial FOV PET scanners.

PURPOSE: Attenuation correction is a critically important step in data correction in positron emission tomography (PET) image formation. The current standard method involves conversion of Hounsfield units from a computed tomography (CT) image to construct attenuation maps (µ-maps) at 511 keV. In this work, the increased sensitivity of long axial field-of-view (LAFOV) PET scanners was exploited to develop and evaluate a deep learning (DL) and joint reconstruction-based method to generate µ-maps utilizing background radiation from lutetium-based (LSO) scintillators. METHODS: Data from 18 subjects were used to train convolutional neural networks to enhance initial µ-maps generated using joint activity and attenuation reconstruction algorithm (MLACF) with transmission data from LSO background radiation acquired before and after the administration of 18F-fluorodeoxyglucose (18F-FDG) (µ-mapMLACF-PRE and µ-mapMLACF-POST respectively). The deep learning-enhanced µ-maps (µ-mapDL-MLACF-PRE and µ-mapDL-MLACF-POST) were compared against MLACF-derived and CT-based maps (µ-mapCT). The performance of the method was also evaluated by assessing PET images reconstructed using each µ-map and computing volume-of-interest based standard uptake value measurements and percentage relative mean error (rME) and relative mean absolute error (rMAE) relative to CT-based method. RESULTS: No statistically significant difference was observed in rME values for µ-mapDL-MLACF-PRE and µ-mapDL-MLACF-POST both in fat-based and water-based soft tissue as well as bones, suggesting that presence of the radiopharmaceutical activity in the body had negligible effects on the resulting µ-maps. The rMAE values µ-mapDL-MLACF-POST were reduced by a factor of 3.3 in average compared to the rMAE of µ-mapMLACF-POST. Similarly, the average rMAE values of PET images reconstructed using µ-mapDL-MLACF-POST (PETDL-MLACF-POST) were 2.6 times smaller than the average rMAE values of PET images reconstructed using µ-mapMLACF-POST. The mean absolute errors in SUV values of PETDL-MLACF-POST compared to PETCT were less than 5% in healthy organs, less than 7% in brain grey matter and 4.3% for all tumours combined. CONCLUSION: We describe a deep learning-based method to accurately generate µ-maps from PET emission data and LSO background radiation, enabling CT-free attenuation and scatter correction in LAFOV PET scanners.

Implementation of a Spatially-Variant and Tissue-Dependent Positron Range Correction for PET/CT Imaging.

Aim: To develop and evaluate a new approach for spatially variant and tissue-dependent positron range (PR) correction (PRC) during the iterative PET image reconstruction. Materials and Methods: The PR distributions of three radionuclides (18F, 68Ga, and 124I) were simulated using the GATE (GEANT4) framework in different material compositions (lung, water, and bone). For every radionuclide, the uniform PR kernel was created by mapping the simulated 3D PR point cloud to a 3D matrix with its size defined by the maximum PR in lung (18F) or water (68Ga and 124I) and the PET voxel size. The spatially variant kernels were composed from the uniform PR kernels by analyzing the material composition of the surrounding medium for each voxel before implementation as tissue-dependent, point-spread functions into the iterative image reconstruction. The proposed PRC method was evaluated using the NEMA image quality phantom (18F, 68Ga, and 124I); two unique PR phantoms were scanned and evaluated following OSEM reconstruction with and without PRC using different metrics, such as contrast recovery, contrast-to-noise ratio, image noise and the resolution evaluated in terms of full width at half maximum (FWHM). Results: The effect of PRC on 18F-imaging was negligible. In contrast, PRC improved image contrast for the 10-mm sphere of the NEMA image quality phantom filled with 68Ga and 124I by 33 and 24%, respectively. While the effect of PRC was less noticeable for the larger spheres, contrast recovery still improved by 5%. The spatial resolution was improved by 26% for 124I (FWHM of 4.9 vs. 3.7 mm). Conclusion: For high energy positron-emitting radionuclides, the proposed PRC method helped recover image contrast with reduced noise levels and with improved spatial resolution. As such, the PRC approach proposed here can help improve the quality of PET data in clinical practice and research.

First results on kinetic modelling and parametric imaging of dynamic 18F-FDG datasets from a long axial FOV PET scanner in oncological patients.

PURPOSE: To investigate the kinetics of 18F-fluorodeoxyglucose (18F-FDG) by positron emission tomography (PET) in multiple organs and test the feasibility of total-body parametric imaging using an image-derived input function (IDIF). METHODS: Twenty-four oncological patients underwent dynamic 18F-FDG scans lasting 65 min using a long  axial FOV (LAFOV) PET/CT system. Time activity curves (TAC) were extracted from semi-automated segmentations of multiple organs, cerebral grey and white matter, and from vascular structures. The tissue and tumor lesion TACs were fitted using an irreversible two-tissue compartment (2TC) and a Patlak model. Parametric images were also generated using direct and indirect Patlak methods and their performances were evaluated. RESULTS: We report estimates of kinetic parameters and metabolic rate of glucose consumption (MRFDG) for different organs and tumor lesions. In some organs, there were significant differences between MRFDG values estimated using 2TC and Patlak models. No statistically significant difference was seen between MRFDG values estimated using 2TC and Patlak methods in tumor lesions (paired t-test, P = 0.65). Parametric imaging showed that net influx (Ki) images generated using direct and indirect Patlak methods had superior tumor-to-background ratio (TBR) to standard uptake value (SUV) images (3.1- and 3.0-fold mean increases in TBRmean, respectively). Influx images generated using the direct Patlak method had twofold higher contrast-to-noise ratio in tumor lesions compared to images generated using the indirect Patlak method. CONCLUSION: We performed pharmacokinetic modelling of multiple organs using linear and non-linear models using dynamic total-body 18F-FDG images. Although parametric images did not reveal more tumors than SUV images, the results confirmed that parametric imaging furnishes improved tumor contrast. We thus demonstrate the feasibility of total-body kinetic modelling and parametric imaging in basic research and oncological studies. LAFOV PET can enhance dynamic imaging capabilities by providing high sensitivity parametric images and allowing total-body pharmacokinetic analysis.

A CT-less approach to quantitative PET imaging using the LSO intrinsic radiation for long-axial FOV PET scanners.

PURPOSE: Long-axial field-of-view (FOV) positron emission tomography (PET) scanners have gained a lot of interest in the recent years. Such scanners provide increased sensitivity and enable unique imaging opportunities that were not previously feasible. Benefiting from the high sensitivity of a long-axial FOV PET scanner, we studied a computed tomography (CT)-less reconstruction algorithm for the Siemens Biograph Vision Quadra with an axial FOV of 106 cm. METHODS: In this work, the background radiation from radioisotope lutetium-176 in the scintillators was used to create an initial estimate of the attenuation maps. Then, joint activity and attenuation reconstruction algorithms were used to create an improved attenuation map of the object. The final attenuation maps were then used to reconstruct quantitative PET images, which were compared against CT-based PET images. The proposed method was evaluated on data from three patients who underwent a flurodeoxyglucouse PET scan. RESULTS: Segmentation of the PET images of the three studied patients showed an average quantitative error of 6.5%-8.3% across all studied organs when using attenuation maps from maximum likelihood estimation of attenuation and activity and 5.3%-6.6% when using attenuation maps from maximum likelihood estimation of activity and attenuation correction coefficients. CONCLUSIONS: Benefiting from the background radiation of lutetium-based scintillators, a quantitative CT-less PET imaging technique was evaluated in this work.

Three-dimensional Fourier-based reprojection analytic reconstruction from histoprojections for high-resolution time-of-flight positron emission tomography scanners.

Purpose: In spite of the general acceptance of iterative reconstruction for clinical use, analytic algorithms provide an important alternative tool due to their linearity, unbiased performance, and predictability for quantitative imaging and quality control studies. On modern time-of-flight (TOF) positron emission tomography scanners with excellent timing resolution, substantial angular compression of (histoprojection) data is possible without loss of resolution, but this also brings challenges for analytical algorithms. We propose TOF and non-TOF Fourier-based analytic approaches that appropriately handle the data sparsity on modern TOF systems. Approach: The proposed TOF algorithm (3D-DIFTOF-direct inversion Fourier transform for TOF) works directly on histoprojection data. The proposed Fourier-based approaches for histoprojection data are further extended to include non-TOF reconstruction (TOF-binned 3D-DIFT), which is particularly useful in time calibration procedures due to its insensitivity to time calibration errors. TOF information is used here to extend available histoprojection data to a larger number of views, essential for artifact-free non-TOF reconstruction. The proposed algorithms are compared with standard analytic techniques on Siemens scanners-space-based confidence-weighted TOF FBP and non-TOF DIFT. Results: 3D-DIFTOF reconstruction demonstrates both improved NEMA-based resolution and contrast versus background variability trade-offs. Similarly, the TOF-binned 3D-DIFT approach shows improved contrast-noise trade-offs over the standard non-TOF approach and is well suited for timing calibration. Conclusions: Our results demonstrate that the proposed 3D-DIFTOF technique provides an improved and more faithful characterization of image resolution compared with standard space-based analytic reconstructions. The proposed tools also provide accurate translation of sparse TOF data available on clinical scanners to upsampled data for non-TOF algorithms.

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.

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.

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.

Plane-dependent ML scatter scaling: 3D extension of the 2D simulated single scatter (SSS) estimate.

Scatter correction is typically done using a simulation of the single scatter, which is then scaled to account for multiple scatters and other possible model mismatches. This scaling factor is determined by fitting the simulated scatter sinogram to the measured sinogram, using only counts measured along LORs that do not intersect the patient body, i.e. 'scatter-tails'. Extending previous work, we propose to scale the scatter with a plane dependent factor, which is determined as an additional unknown in the maximum likelihood (ML) reconstructions, using counts in the entire sinogram rather than only the 'scatter-tails'. The ML-scaled scatter estimates are validated using a Monte-Carlo simulation of a NEMA-like phantom, a phantom scan with typical contrast ratios of a 68Ga-PSMA scan, and 23 whole-body 18F-FDG patient scans. On average, we observe a 12.2% change in the total amount of tracer activity of the MLEM reconstructions of our whole-body patient database when the proposed ML scatter scales are used. Furthermore, reconstructions using the ML-scaled scatter estimates are found to eliminate the typical 'halo' artifacts that are often observed in the vicinity of high focal uptake regions.

Metal Artifact Reduction of CT Scans to Improve PET/CT.

In recent years, different metal artifact reduction methods have been developed for CT. These methods have only recently been introduced for PET/CT even though they could be beneficial for interpretation, segmentation, and quantification of the PET/CT images. In this study, phantom and patient scans were analyzed visually and quantitatively to measure the effect on PET images of iterative metal artifact reduction (iMAR) of CT data. Methods: The phantom consisted of 2 types of hip prostheses in a solution of 18F-FDG and water. 18F-FDG PET/CT scans of 14 patients with metal implants (either dental implants, hip prostheses, shoulder prostheses, or pedicle screws) and 68Ga-labeled prostate-specific membrane antigen (68Ga-PSMA) PET/CT scans of 7 patients with hip prostheses were scored by 2 experienced nuclear medicine physicians to analyze clinical relevance. For all patients, a lesion was located in the field of view of the metal implant. Phantom and patients were scanned in a PET/CT scanner. The standard low-dose CT scans were processed with the iMAR algorithm. The PET data were reconstructed using attenuation correction provided by both standard CT and iMAR-processed CT. Results: For the phantom scans, cold artifacts were visible on the PET image. There was a 30% deficit in 18F-FDG concentration, which was restored by iMAR processing, indicating that metal artifacts on CT images induce quantification errors in PET data. The iMAR algorithm was useful for most patients. When iMAR was used, the confidence in interpretation increased or stayed the same, with an average improvement of 28% ± 20% (scored on a scale of 0%-100% confidence). The SUV increase or decrease depended on the type of metal artifact. The mean difference in absolute values of SUVmean of the lesions was 3.5% ± 3.3%. Conclusion: The iMAR algorithm increases the confidence of the interpretation of the PET/CT scan and influences the SUV. The added value of iMAR depends on the indication for the PET/CT scan, location and size/type of the prosthesis, and location and extent of the disease.

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.

LSO background radiation as a transmission source using time of flight.

LSO scintillators (Lu2Sio5:Ce) have a background radiation which originates from the isotope Lu-176 that is present in natural occurring lutetium. The decay that occurs in this isotope is a beta decay that is in coincidence with cascade gamma emissions with energies of 307,202 and 88 keV. The coincidental nature of the beta decay with the gamma emissions allow for separation of emission data originating from a positron annihilation event from transmission type data from the Lu-176 beta decay. By using the time of flight information, and information of the chord length between two LSO pixels in coincidence as a result of a beta emission and emitted gamma, a second time window can be set to observe transmission events simultaneously to emission events. Using the time when the PET scanner is not actively acquiring positron emission data, a continuous blank can be acquired and used to reconstruct a transmission image. With this blank and the measured transmission data, a transmission image can be reconstructed. This reconstructed transmission image can be used to perform emission data corrections such as attenuation correction and scatter corrections or starting images for algorithms that estimate emission and attenuation simultaneously. It is observed that the flux of the background activity is high enough to create useful transmission images with an acquisition time of 10 min.

List-mode reconstruction for the Biograph mCT with physics modeling and event-by-event motion correction.

Whole-body PET/CT scanners are important clinical and research tools to study tracer distribution throughout the body. In whole-body studies, respiratory motion results in image artifacts. We have previously demonstrated for brain imaging that, when provided with accurate motion data, event-by-event correction has better accuracy than frame-based methods. Therefore, the goal of this work was to develop a list-mode reconstruction with novel physics modeling for the Siemens Biograph mCT with event-by-event motion correction, based on the MOLAR platform (Motion-compensation OSEM List-mode Algorithm for Resolution-Recovery Reconstruction). Application of MOLAR for the mCT required two algorithmic developments. First, in routine studies, the mCT collects list-mode data in 32 bit packets, where averaging of lines-of-response (LORs) by axial span and angular mashing reduced the number of LORs so that 32 bits are sufficient to address all sinogram bins. This degrades spatial resolution. In this work, we proposed a probabilistic LOR (pLOR) position technique that addresses axial and transaxial LOR grouping in 32 bit data. Second, two simplified approaches for 3D time-of-flight (TOF) scatter estimation were developed to accelerate the computationally intensive calculation without compromising accuracy. The proposed list-mode reconstruction algorithm was compared to the manufacturer's point spread function + TOF (PSF+TOF) algorithm. Phantom, animal, and human studies demonstrated that MOLAR with pLOR gives slightly faster contrast recovery than the PSF+TOF algorithm that uses the average 32 bit LOR sinogram positioning. Moving phantom and a whole-body human study suggested that event-by-event motion correction reduces image blurring caused by respiratory motion. We conclude that list-mode reconstruction with pLOR positioning provides a platform to generate high quality images for the mCT, and to recover fine structures in whole-body PET scans through event-by-event motion correction.

An assessment of the impact of incorporating time-of-flight information into clinical PET/CT imaging.

The introduction of fast scintillators with good stopping power for 511-keV photons has renewed interest in time-of-flight (TOF) PET. The ability to measure the difference between the arrival times of a pair of photons originating from positron annihilation improves the image signal-to-noise ratio (SNR). The level of improvement depends upon the extent and distribution of the positron activity and the time resolution of the PET scanner. While specific estimates can be made for phantom imaging, the impact of TOF PET is more difficult to quantify in clinical situations. The results presented here quantify the benefit of TOF in a challenging phantom experiment and then assess both qualitatively and quantitatively the impact of incorporating TOF information into the reconstruction of clinical studies. A clear correlation between patient body mass index and gain in SNR was observed in this study involving 100 oncology patient studies, with a gain due to TOF ranging from 1.1 to 1.8, which is consistent with the 590-ps time resolution of the TOF PET scanner. The visual comparison of TOF and non-TOF images performed by two nuclear medicine physicians confirmed the advantages of incorporating TOF into the reconstruction, advantages that include better definition of small lesions and image details, improved uniformity, and noise reduction.

Experimental comparison of lesion detectability for four fully-3D PET reconstruction schemes.

The objective of this work was to evaluate the lesion detection performance of four fully-3D positron emission tomography (PET) reconstruction schemes using experimentally acquired data. A multi-compartment anthropomorphic phantom was set up to mimic whole-body (18)F-fluorodeoxyglucose (FDG) cancer imaging and scanned 12 times in 3D mode, obtaining count levels typical of noisy clinical scans. Eight of the scans had 26 (68)Ge "shell-less" lesions (6, 8-, 10-, 12-, 16-mm diameter) placed throughout the phantom with various target:background ratios. This provided lesion-present and lesion-absent datasets with known truth appropriate for evaluating lesion detectability by localization receiver operating characteristic (LROC) methods. Four reconstruction schemes were studied: 1) Fourier rebinning (FORE) followed by 2D attenuation-weighted ordered-subsets expectation-maximization, 2) fully-3D AW-OSEM, 3) fully-3D ordinary-Poisson line-of-response (LOR-)OSEM; and 4) fully-3D LOR-OSEM with an accurate point-spread function (PSF) model. Two forms of LROC analysis were performed. First, a channelized nonprewhitened (CNPW) observer was used to optimize processing parameters (number of iterations, post-reconstruction filter) for the human observer study. Human observers then rated each image and selected the most-likely lesion location. The area under the LROC curve ( A(LROC)) and the probability of correct localization were used as figures-of-merit. The results of the human observer study found no statistically significant difference between FORE and AW-OSEM3D ( A(LROC)=0.41 and 0.36, respectively), an increase in lesion detection performance for LOR-OSEM3D ( A(LROC)=0.45, p=0.076), and additional improvement with the use of the PSF model ( A(LROC)=0.55, p=0.024). The numerical CNPW observer provided the same rankings among algorithms, but obtained different values of A(LROC). These results show improved lesion detection performance for the reconstruction algorithms with more sophisticated statistical and imaging models as compared to the previous-generation algorithms.

Continuous and discrete data rebinning in time-of-flight PET.

This paper investigates data compression methods for time-of-flight (TOF) positron emission tomography (PET), which rebin the 3-D TOF measurements into a set of 2-D TOF data for a stack of transaxial slices. The goal of this work is to develop rebinning algorithms that are more accurate than the TOF single-slice-rebinning (TOF-SSRB) method proposed by Mullani in 1982. Two approaches are explored. The first one is based on a partial differential equation, which expresses a consistency condition for TOF-PET data with a Gaussian TOF profile. From this equation we derive an analytical rebinning algorithm, which is unbiased in the limit of continuous sampling. The second approach is discrete: each 2-D rebinned data sample is calculated as a linear combination of the 3-D TOF samples in the same axial plane parallel to the axis of the scanner. The coefficients of the linear combination are precomputed by optimizing a cost function which enforces both accuracy and good variance reduction, models the TOF profile, the axial PSF of the LORs, and the specific sampling scheme of the scanner. Measurements of a thorax phantom on a prototype TOF-PET scanner with a resolution of 550 ps show that the proposed discrete method improves the bias-variance trade-off and is a promising alternative to TOF-SSRB when data compression is required to achieve clinically acceptable reconstruction time.

Fully 3-D PET reconstruction with system matrix derived from point source measurements.

The quality of images reconstructed by statistical iterative methods depends on an accurate model of the relationship between image space and projection space through the system matrix The elements of the system matrix for the clinical Hi-Rez scanner were derived by processing the data measured for a point source at different positions in a portion of the field of view. These measured data included axial compression and azimuthal interleaving of adjacent projections. Measured data were corrected for crystal and geometrical efficiency. Then, a whole system matrix was derived by processing the responses in projection space. Such responses included both geometrical and detection physics components of the system matrix. The response was parameterized to correct for point source location and to smooth for projection noise. The model also accounts for axial compression (span) used on the scanner. The forward projector for iterative reconstruction was constructed using the estimated response parameters. This paper extends our previous work to fully three-dimensional. Experimental data were used to compare images reconstructed by the standard iterative reconstruction software and the one modeling the response function. The results showed that the modeling of the response function improves both spatial resolution and noise properties.