GPU-accelerated Monte Carlo based scatter correction in brain PET/MR

Development of a Monte Carlo-based scatter correction method for total-body PET using the uEXPLORER PET/CT scanner

Objective.This study presents and evaluates a robust Monte Carlo-based scatter correction (SC) method for long axial field of view (FOV) and total-body positron emission tomography (PET) using the uEXPLORER total-body PET/CT scanner.Approach.Our algorithm utilizes the Monte Carlo (MC) tool SimSET to compute SC factors in between individual image reconstruction iterations within our in-house list-mode and time-of-flight-based image reconstruction framework. We also introduced a unique scatter scaling technique at the detector block-level for optimal estimation of the scatter contribution in each line of response. First image evaluations were derived from phantom data spanning the entire axial FOV along with image data from a human subject with a large body mass index. Data was evaluated based on qualitative inspections, and contrast recovery, background variability, residual scatter removal from cold regions, biases and axial uniformity were quantified and compared to non-scatter-corrected images.Main results.All reconstructed images demonstrated qualitative and quantitative improvements compared to non-scatter-corrected images: contrast recovery coefficients improved by up to 17.2% and background variability was reduced by up to 34.3%, and the residual lung error was between 1.26% and 2.08%. Low biases throughout the axial FOV indicate high quantitative accuracy and axial uniformity of the corrections. Up to 99% of residual activity in cold areas in the human subject was removed, and the reliability of the method was demonstrated in challenging body regions like in the proximity of a highly attenuating knee prosthesis.Significance.The MC SC method employed was demonstrated to be accurate and robust in TB-PET. The results of this study can serve as a benchmark for optimizing the quantitative performance of future SC techniques.

PET scatter estimation using deep learning U-Net architecture

Objective.Positron emission tomography (PET) image reconstruction needs to be corrected for scatter in order to produce quantitatively accurate images. Scatter correction is traditionally achieved by incorporating an estimated scatter sinogram into the forward model during image reconstruction. Existing scatter estimated methods compromise between accuracy and computing time. Nowadays scatter estimation is routinely performed using single scatter simulation (SSS), which does not accurately model multiple scatter and scatter from outside the field-of-view, leading to reduced qualitative and quantitative PET reconstructed image accuracy. On the other side, Monte-Carlo (MC) methods provide a high precision, but are computationally expensive and time-consuming, even with recent progress in MC acceleration.Approach.In this work we explore the potential of deep learning (DL) for accurate scatter correction in PET imaging, accounting for all scatter coincidences. We propose a network based on a U-Net convolutional neural network architecture with 5 convolutional layers. The network takes as input the emission and computed tomography (CT)-derived attenuation factor (AF) sinograms and returns the estimated scatter sinogram. The network training was performed using MC simulated PET datasets. Multiple anthropomorphic extended cardiac-torso phantoms of two different regions (lung and pelvis) were created, considering three different body sizes and different levels of statistics. In addition, two patient datasets were used to assess the performance of the method in clinical practice.Main results.Our experiments showed that the accuracy of our method, namely DL-based scatter estimation (DLSE), was independent of the anatomical region (lungs or pelvis). They also showed that the DLSE-corrected images were similar to that reconstructed from scatter-free data and more accurate than SSS-corrected images.Significance.The proposed method is able to estimate scatter sinograms from emission and attenuation data. It has shown a better accuracy than the SSS, while being faster than MC scatter estimation methods.

Influences of 3D PET scanner components on increased scatter evaluated by a Monte Carlo simulation

Monte Carlo simulation is widely applied to evaluate the performance of three-dimensional positron emission tomography (3D-PET). For accurate scatter simulations, all components that generate scatter need to be taken into account. The aim of this work was to identify the components that influence scatter. The simulated geometries of a PET scanner were: a precisely reproduced configuration including all of the components; a configuration with the bed, the tunnel and shields; a configuration with the bed and shields; and the simplest geometry with only the bed. We measured and simulated the scatter fraction using two different set-ups: (1) as prescribed by NEMA-NU 2007 and (2) a similar set-up but with a shorter line source, so that all activity was contained only inside the field-of-view (FOV), in order to reduce influences of components outside the FOV. The scatter fractions for the two experimental set-ups were, respectively, 45% and 38%. Regarding the geometrical configurations, the former two configurations gave simulation results in good agreement with the experimental results, but simulation results of the simplest geometry were significantly different at the edge of the FOV. From the simulation of the precise configuration, the object (scatter phantom) was the source of more than 90% of the scatter. This was also confirmed by visualization of photon trajectories. Then, the bed and the tunnel were mainly the sources of the rest of the scatter. From the simulation results, we concluded that the precise construction was not needed; the shields, the tunnel, the bed and the object were sufficient for accurate scatter simulations.

An energy minimization method for the correction of cupping artifacts in cone-beam CT

The purpose of this study was to reduce cupping artifacts and improve quantitative accuracy of the images in cone-beam CT (CBCT). An energy minimization method (EMM) is proposed to reduce cupping artifacts in reconstructed image of the CBCT. The cupping artifacts are iteratively optimized by using efficient matrix computations, which are verified to be numerically stable by matrix analysis. Moreover, the energy in our formulation is convex in each of its variables, which brings the robustness of the proposed energy minimization algorithm. The cupping artifacts are estimated as a result of minimizing this energy. The results indicate that proposed algorithm is effective for reducing the cupping artifacts and preserving the quality of the reconstructed image. The proposed method focuses on the reconstructed image without requiring any additional physical equipment; it is easily implemented and provides cupping correction using a single scan acquisition. The experimental results demonstrate that this method can successfully reduce the magnitude of cupping artifacts. The correction algorithm reported here may improve the uniformity of the reconstructed images, thus assisting the development of perfect volume visualization and threshold-based visualization techniques for reconstructed images.

Innovations in Small-Animal PET/MR Imaging Instrumentation

Multimodal imaging has led to a more detailed exploration of different physiologic processes with integrated PET/MR imaging being the most recent entry. Although the clinical need is still questioned, it is well recognized that it represents one of the most active and promising fields of medical imaging research in terms of software and hardware. The hardware developments have moved from small detector components to high-performance PET inserts and new concepts in full systems. Conversely, the software focuses on the efficient performance of necessary corrections without the use of CT data. The most recent developments in both directions are reviewed.