Performance evaluation of a 90°-rotating dual-head small animal PET system

Hybrid model of photon propagation based on the analytical and Monte Carlo methods for a dual-head PET system

The construction of photon propagation has a close relationship with the quality of reconstructed images. The classical Monte Carlo (MC) based method can model the photon propagation precisely, but it is time-consuming. The analytical method can often quickly construct a model, but its precision is a problem. How to fully exploit the advantages of the MC simulation and analytical model is an open problem. Inspired by the characteristics of the depth of interaction (DOI) detectors, which can help confirm the deposited position of a photon with DOI-encoding technology, we virtually discretize each crystal into several subcrystals to obtain the statistical distribution by MC-based simulation. Then, the statistical distribution is combined with a spatially variant solid-angle model. This combination strategy provides a hybrid model to describe photon propagation with relatively high accuracy and low computational cost. Three discretization schemes are compared to optimize the constructed photon propagation model. Several experiments are carried out to evaluate the performance of the proposed hybrid method. The metrics of full width at half maximum (FWHM), contrast recovery (CR), and coefficient of variation (COV) are adopted to quantitate the imaging results. The classical MC-based method is compared as a gold-standard reference. When a crystal is divided into two discretized positions, the convergent tendencies of CRs and COVs are consistent with that based on MC simulation method, respectively. In terms of FWHMs, the resolutions of using the MC-based model and the proposed hybrid model are 0.71 mm and 0.68 mm in the direction parallel to the detector head, respectively. This indicates the potential of the proposed method in positron emission tomography imaging.

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.

Influence of Rotation Increments on Imaging Performance for a Rotatory Dual-Head PET System

For a rotatory dual-head positron emission tomography (PET) system, how to determine the rotation increments is an open problem. In this study, we simulated the characteristics of a rotatory dual-head PET system. The influences of different rotation increments were compared and analyzed. Based on this simulation, the imaging performance of a prototype system was verified. A reconstruction flowchart was proposed based on a precalculated system response matrix (SRM). The SRM made the relationships between the voxels and lines of response (LORs) fixed; therefore, we added the interpolation method into the flowchart. Five metrics, including spatial resolution, normalized mean squared error (NMSE), peak signal-to-noise ratio (PSNR), contrast-to-noise (CNR), and structure similarity (SSIM), were applied to assess the reconstructed image quality. The results indicated that the 60° rotation increments with the bilinear interpolation had advantages in resolution, PSNR, NMSE, and SSIM. In terms of CNR, the 90° rotation increments were better than other increments. In addition, the reconstructed images of 90° rotation increments were also flatter than that of 60° increments. Therefore, both the 60° and 90° rotation increments could be used in the real experiments, and which one to choose may depend on the application requirement.