Ultra Fast Symmetry and SIMD-Based Projection-Backprojection (SSP) Algorithm for 3-D PET Image Reconstruction

Abstract

Remarkable progress in positron emission tomography (PET) development has occurred in recent years, in hardware, software, and computer implementation of image reconstruction. Recent development in PET scanners such as the high-resolution research tomograph (HRRT) developed by CTI (now Siemens) represents such a case and is capable of greatly enhanced resolution as well as sensitivity. In these PET scanners, the amount of coincidence line data collected contains more than 4.5 x 10(9) coincidence lines of response generated by as many nuclear detectors as 120 000. This formidable amount of data and the reconstruction of this data set pose a real problem in HRRT and have also been of the major bottle neck in further developments of high resolution PET scanners as well as their applications. In these classes of PET scanners, therefore, obtaining one set of reconstructed images often requires many hours of image reconstruction. For example, in HRRT with full data collection in a normal brain scan (using SPAN 3), the image reconstruction time is close to 80 min, making it practically impossible to attempt any list-mode-based dynamic imaging since the image reconstruction time would take many days (as much as 43 h or more for 32-frame dynamic image reconstruction). To remedy this data-handling problem in image reconstruction, we developed a new algorithm based on the symmetry properties of the projection and backprojection processes, especially in the 3-D OSEM algorithm where multiples of projection and back-projection are required. In addition, the single-instruction multiple-data (SIMD) technique also allowed us to successfully incorporate the symmetry properties mentioned above, thereby effectively reducing the total image reconstruction time to a few minutes. We refer to this technique as the symmetry and SIMD-based projection-backprojection (SSP) technique or algorithm and the details of the technique will be discussed and an example of the application of the technique to the HRRT's OSEM algorithm will be presented as a demonstration.

Shine-Through in PET/MR Imaging: Effects of the Magnetic Field on Positron Range and Subsequent Image Artifacts

“…The spatial resolution in the center was 2.3, 3.1, and 2.6 mm in full width at half maximum in the x, y, and z directions, respectively (16). Acquired list-mode data were reconstructed using an ordinary Poisson ordered-subsets expectation maximization 3-dimensional algorithm (6 iterations and 16 subsets) (17). The z-axis of the scanner was defined as aligned with the magnetic field.…”

18F-Florbetapir PET in Patients with Frontotemporal Dementia and Alzheimer Disease

“…PET data were acquired for a total of 60 min after injection in list-mode and were reconstructed using an ultra-fast ordered-subsets expectation maximization algorithm (15) (matrix size, 256 · 256 · 207; voxel size, 1.22 · 1.22 · 1.22 mm). Centroid of motion data were used to correct for head motion by realignment of frames as previously described (16).…”

Measurement of Striatal and Extrastriatal Dopamine Transporter Binding with High-Resolution PET and [11C]PE2I: Quantitative Modeling and Test—Retest Reproducibility

“…Emission data were acquired in list mode and histogrammed using axial compression of span 9, yielding noise -equivalent counts in the range of 3 Â 10 6 to 2.5 Â 10 7 , except in the first frame. Image estimates with isotropic voxel dimensions of 1.22 Â 1.22 Â 1.22 mm 3 were generated using speed-optimized version of OP-OSEM-3D (Ordinary Poisson-OSEM in full 3D) reconstruction (Comtat et al, 2004;Hong et al, 2007) with 16 subsets and 8 full iterations.…”