PET Image Reconstruction Incorporating 3D Mean-Median Sinogram Filtering

Q.Clear reconstruction for reducing the scanning time for 68 Ga-DOTA-FAPI-04 PET/MR imaging

PURPOSE

This study aims to determine whether Q.Clear positron emission tomography (PET) reconstruction may reduce tracer injection dose or shorten scanning time in 68Gallium-labelled fibroblast activation protein inhibitor (68 Ga-FAPI) PET/magnetic resonance (MR) imaging.

METHODS

We retrospectively collected cases of 68 Ga-FAPI whole-body imaging performed on integrated PET/MR. PET images were reconstructed using three different methods: ordered subset expectation maximization (OSEM) reconstruction with full scanning time, OSEM reconstruction with half scanning time, and Q.Clear reconstruction with half scanning time. We then measured standardized uptake values (SUVs) within and around lesions, alongside their volumes. We also evaluated image quality using lesion-to-background (L/B) ratio and signal-to-noise ratio (SNR). We then compared these metrics across the three reconstruction techniques using statistical methods.

RESULTS

Q.Clear reconstruction significantly increased SUVmax and SUVmean within lesions (more than 30%) and reduced their volumes in comparison with OSEM reconstruction. Background SUVmax also increased significantly, while background SUVmean showed no difference. Average L/B values for Q.Clear reconstruction were only marginally higher than those from OSME reconstruction with half-time. SNR decreased significantly in Q.Clear reconstruction compared with OSEM reconstruction with full time (but not half time). Differences between Q.Clear and OSEM reconstructions in SUVmax and SUVmean values within lesions were significantly correlated with SUVs within lesions.

CONCLUSIONS

Q.Clear reconstruction was useful for reducing PET injection dose or scanning time while maintaining the image quality. Q.Clear may affect PET quantification, and it is necessary to establish diagnostic recommendations based on Q.Clear results for Q.Clear application.

Masked Joint Bilateral Filtering via Deep Image Prior for Digital X-ray Image Denoising

Medical image denoising faces great challenges. Although deep learning methods have shown great potential, their efficiency is severely affected by millions of trainable parameters. The non-linearity of neural networks also makes them difficult to be understood. Therefore, existing deep learning methods have been sparingly applied to clinical tasks. To this end, we integrate known filtering operators into deep learning and propose a novel Masked Joint Bilateral Filtering (MJBF) via deep image prior for digital X-ray image denoising. Specifically, MJBF consists of a deep image prior generator and an iterative filtering block. The deep image prior generator produces plentiful image priors by a multi-scale fusion network. The generated image priors serve as the guidance for the iterative filtering block, which is utilized for the actual edge-preserving denoising. The iterative filtering block contains three trainable Joint Bilateral Filters (JBFs), each with only 18 trainable parameters. Moreover, a masking strategy is introduced to reduce redundancy and improve the understanding of the proposed network. Experimental results on the ChestX-ray14 dataset and real data show that the proposed MJBF has achieved superior performance in terms of noise suppression and edge preservation. Tests on the portability of the proposed method demonstrate that this denoising modality is simple yet effective, and could have a clinical impact on medical imaging in the future.

Comparison of GPU reconstruction based on different symmetries for dual-head PET

Purpose

Dual‐head positron emission tomography (PET) scanners have increasingly attracted the attention of many researchers. However, with the compact geometry, the depth‐of‐interaction blurring will reduce the image resolution considerably. Monte Carlo (MC)‐based system response matrix (SRM) is able to describe the physical process of PET imaging accurately and improve reconstruction quality significantly. The MC‐based SRM is large and precomputed, which leads to a longer image reconstruction time with indexing and retrieving precomputed system matrix elements. In this study, we proposed a GPU acceleration algorithm to accelerate the iterative reconstruction.

Methods

It has been demonstrated that the line‐of‐response (LOR)‐based symmetry and the Graphics Processing Unit (GPU) technology can accelerate the reconstruction tremendously. LOR‐based symmetry is suitable for the forward projection calculation, but not for the backprojection. In this study, we proposed a GPU acceleration algorithm that combined the LOR‐based symmetry and voxel‐based symmetry together, in which the LOR‐based symmetry is responsible for the forward projection, and the voxel‐based symmetry is used for the backprojection.

Results

Simulation and real experiments verify the efficiency of the algorithm. Compared with the CPU‐based calculation, the acceleration ratios of the forward projection and the backprojection operation are 130 and 110, respectively. The total acceleration ratio is 113×. In order to compare the acceleration effect of the different symmetries, we realized the reconstruction with the voxel‐based symmetry and the LOR‐based symmetry strategies. Compared with the LOR‐based GPU reconstruction, the acceleration ratio is 3.5×. Compared with the voxel‐based GPU reconstruction, the acceleration ratio is 12×.

Conclusion

We have proposed a new acceleration algorithm for the dual‐head PET system, in which both the forward and backprojection operations are accelerated by GPU.

Improved Simultaneous Algebraic Reconstruction Technique Algorithm for Positron-Emission Tomography Image Reconstruction via Minimizing the Fast Total Variation

CONTEXT

There has been considerable progress in the instrumentation for data measurement and computer methods for generating images of measured PET data. These computer methods have been developed to solve the inverse problem, also known as the "image reconstruction from projections" problem.

AIM

In this paper, we propose a modified Simultaneous Algebraic Reconstruction Technique (SART) algorithm to improve the quality of image reconstruction by incorporating total variation (TV) minimization into the iterative SART algorithm.

METHODOLOGY

The SART updates the estimated image by forward projecting the initial image onto the sinogram space. Then, the difference between the estimated sinogram and the given sinogram is back-projected onto the image domain. This difference is then subtracted from the initial image to obtain a corrected image. Fast total variation (FTV) minimization is applied to the image obtained in the SART step. The second step is the result obtained from the previous FTV update. The SART and the FTV minimization steps run iteratively in an alternating manner. Fifty iterations were applied to the SART algorithm used in each of the regularization-based methods. In addition to the conventional SART algorithm, spatial smoothing was used to enhance the quality of the image. All images were sized at 128 × 128 pixels.

RESULTS

The proposed algorithm successfully accomplished edge preservation. A detailed scrutiny revealed that the reconstruction algorithms differed; for example, the SART and the proposed FTV-SART algorithm effectively preserved the hot lesion edges, whereas artifacts and deviations were more likely to occur in the ART algorithm than in the other algorithms.

CONCLUSIONS

Compared to the standard SART, the proposed algorithm is more robust in removing background noise while preserving edges to suppress the existent image artifacts. The quality measurements and visual inspections show a significant improvement in image quality compared to the conventional SART and Algebraic Reconstruction Technique (ART) algorithms.

Automatic cumulative sums contour detection of FBP-reconstructed multi-object nuclear medicine images

The problem of determining the contours of objects in nuclear medicine images has been studied extensively in the past, however most of the analysis has focused on a single object as opposed to multiple objects. The aim of this work is to develop an automated method for determining the contour of multiple objects in positron emission tomography (PET) and single photon emission computed tomography (SPECT) filtered backprojection (FBP) reconstructed images. These contours can be used for computing body edges for attenuation correction in PET and SPECT, as well as for eliminating streak artifacts outside the objects, which could be useful in compressive sensing reconstruction. Contour detection has been accomplished by applying a modified cumulative sums (CUSUM) scheme in the sinogram. Our approach automatically detects all objects in the image, without requiring a priori knowledge of the number of distinct objects in the reconstructed image. This method has been tested in simulated phantoms, such as an image-quality (IQ) phantom and two digital multi-object phantoms, as well as a real NEMA phantom and a clinical thoracic study. For this purpose, a GE Discovery PET scanner was employed. The detected contours achieved root mean square accuracy of 1.14 pixels, 1.69 pixels and 3.28 pixels and a Hausdorff distance of 3.13, 3.12 and 4.50 pixels, for the simulated image-quality phantom PET study, the real NEMA phantom and the clinical thoracic study, respectively. These results correspond to a significant improvement over recent results obtained in similar studies. Furthermore, we obtained an optimal sub-pattern assignment (OSPA) localization error of 0.94 and 1.48, for the two-objects and three-objects simulated phantoms, respectively. Our method performs efficiently for sets of convex objects and hence it provides a robust tool for automatic contour determination with precise results.