High Performance 3D PET Reconstruction Using Spherical Basis Functions on a Polar Grid

Qualitative and Quantitative Evaluation of Blob-Based Time-of-Flight PET Image Reconstruction in Hybrid Brain PET/MR Imaging

Purpose

Many neurological diseases affect small structures in the brain and, as such, reliable visual evaluation and accurate quantification are required. Recent technological developments made the clinical use of hybrid positron emission tomography/magnetic resonance (PET/MR) systems possible, providing both functional and anatomical information in a single imaging session. Nevertheless, there is a trade-off between spatial resolution and image quality (contrast and noise), which is dictated mainly by the chosen acquisition and reconstruction protocols. Image reconstruction algorithms using spherical symmetric basis functions (blobs) for image representation have a number of additional parameters that impact both the qualitative and quantitative image characteristics. Hence, a detailed investigation of the blob-based reconstruction characteristics using different parameters is needed to achieve optimal reconstruction results.

Procedures

This work evaluated the impact of a range of blob parameters on image quality and quantitative accuracy of brain PET images acquired on the Ingenuity Time-of-Flight (TOF) PET/MR system. Two different phantoms were used to simulate brain imaging applications. Image contrast and noise characteristics were assessed using an image quality phantom. Quantitative performance in a clinical setting was investigated using the Hoffman 3D brain phantom at various count levels. Furthermore, the visual quality of four clinical studies was scored blindly by two experienced physicians to qualitatively evaluate the influence of different reconstruction protocols, hereby providing indications on parameters producing the best image quality.

Results

Quantitative evaluation using the image quality phantom showed that larger basis function radii result in lower contrast recovery (∼2 %) and lower variance levels (∼15 %). The brain phantom and clinical studies confirmed these observations since lower contrast was seen between anatomical structures. High and low count statistics gave comparable values. The qualitative evaluation of the clinical studies, based on the assessment performed by the physicians, showed a preference towards lower image variance levels with a slightly lower contrast, favoring higher radii and four iterations.

Conclusion

This study systematically evaluated a number of basis function parameters and their quantitative and qualitative effect within PET image reconstruction in the context of brain imaging. A range of blob parameters can minimize error and provided optimal image quality, where the anatomical structures could be recognized but the exact delineation of these structures is simplified in scans with lower image variance levels and thus, higher radii should be preferred. With the optimization of blob parameters, the reconstructed images were found to be qualitatively improved (optimum parameters {d = 2.0375, alpha = 10.4101, radius = 3.9451}) as assessed by the physicians compared to the current clinical protocol. However, this qualitative improvement is task specific, depending on the desired image characteristics to be extracted. Finally, this work could be used as a guide for application-specific optimal parameter selection.

Electronic supplementary material

The online version of this article (doi:10.1007/s11307-015-0824-x) contains supplementary material, which is available to authorized users.

GPU-based high-performance computing for radiation therapy

Recent developments in radiotherapy therapy demand high computation powers to solve challenging problems in a timely fashion in a clinical environment. The graphics processing unit (GPU), as an emerging high-performance computing platform, has been introduced to radiotherapy. It is particularly attractive due to its high computational power, small size, and low cost for facility deployment and maintenance. Over the past few years, GPU-based high-performance computing in radiotherapy has experienced rapid developments. A tremendous amount of study has been conducted, in which large acceleration factors compared with the conventional CPU platform have been observed. In this paper, we will first give a brief introduction to the GPU hardware structure and programming model. We will then review the current applications of GPU in major imaging-related and therapy-related problems encountered in radiotherapy. A comparison of GPU with other platforms will also be presented.

Medical image processing on the GPU - past, present and future

Graphics processing units (GPUs) are used today in a wide range of applications, mainly because they can dramatically accelerate parallel computing, are affordable and energy efficient. In the field of medical imaging, GPUs are in some cases crucial for enabling practical use of computationally demanding algorithms. This review presents the past and present work on GPU accelerated medical image processing, and is meant to serve as an overview and introduction to existing GPU implementations. The review covers GPU acceleration of basic image processing operations (filtering, interpolation, histogram estimation and distance transforms), the most commonly used algorithms in medical imaging (image registration, image segmentation and image denoising) and algorithms that are specific to individual modalities (CT, PET, SPECT, MRI, fMRI, DTI, ultrasound, optical imaging and microscopy). The review ends by highlighting some future possibilities and challenges.