Time-of-flight PET image reconstruction using origin ensembles

Time of flight dual photon emission computed tomography

Time-of-flight dual photon emission computed tomography (TOF-DuPECT) is an imaging system that can obtain radionuclide distributions using time information recorded from two cascade-decay photons. The potential decay locations in the image space, a hyperbolic response curve, can be determined via time-difference-of-arrival (TDOA) estimations from two instantaneous coincidence photons. In this feasibility study, Monte Carlo simulations were performed to generate list-mode coincidence data. A full-ring positron emission tomography-like detection system geometry was built in the simulation environment. A contrast phantom and a Jaszczak-like phantom filled with Selenium-75 (Se-75) were used to evaluate the image quality. A TOF-DuPECT system with varying coincidence time resolution (CTR) was then evaluated. We used the stochastic origin ensemble (SOE) algorithm to reconstruct images from the recorded list-mode data. The results indicate that the SOE method can be successfully employed for the TOF-DuPECT system and can achieve acceptable image quality when the CTR is less than 100 ps. Therefore, the TOF-DuPECT imaging system is feasible. With the improvement of the detector with time, future implementations and applications of TOF-DuPECT are promising. Further quantitative imaging techniques such as attenuation and scatter corrections for the TOF-DuPECT system will be developed in future.

A noise smoothing origin ensemble algorithm based on local filtering

An origin ensemble (OE) image reconstruction algorithm can be used for the fast reconstruction of unconventional geometrical images, e.g. in a Compton camera (CC) system. Due to the low-count rate in the emission data, the reconstructed image is often noisy and inhomogeneous in density. In this study, we propose a way to smooth out the noise in the OE algorithm. During the OE reconstruction, the algorithm stochastically modifies the current location to a random new voxel along the probable corresponding curve of each event depending on the relative event density of the new and old locations. In the original OE technique, the event density is simply the number of events in the voxel. In the proposed method, the event density is estimated from the filtering of a kernel window centered on the voxel. Incorporating the regional filtering is similar to performing an OE algorithm on a smoothed image at each iteration and enables the reconstruction of a smoother image. A Flangeless Esser PET phantom and a multi-activity phantom are used to study the property of the new reconstruction algorithm. The results indicate that the proposed method performs better than a conventional OE algorithm in terms of normalized mean square error (NMSE) and structural similarity (SSIM). Both contrast noise ratio (CNR) and reconstruction accuracy of the new method are better than the conventional OE algorithm and their performances improve with the increase of object size. The median-OE possesses the highest overall image quality and recovery rate among the three filter-OE algorithms and is the method of choice for image reconstruction. Comparing to conventional post-smoothing OEs, the NMSE of median-OE improves 57.6% (46.9%) and the SSIM increased by 73.2% (51.1%) for the Esser (multi-activity) phantom. The proposed OE algorithm is simple and efficient for noise smoothing without complex calculations and highly suited for low-count cases.

Evaluation of Single-Chip, Real-Time Tomographic Data Processing on FPGA SoC Devices

A novel approach to tomographic data processing has been developed and evaluated using the Jagiellonian positron emission tomography scanner as an example. We propose a system in which there is no need for powerful, local to the scanner processing facility, capable to reconstruct images on the fly. Instead, we introduce a field programmable gate array system-on-chip platform connected directly to data streams coming from the scanner, which can perform event building, filtering, coincidence search, and region-of-response reconstruction by the programmable logic and visualization by the integrated processors. The platform significantly reduces data volume converting raw data to a list-mode representation, while generating visualization on the fly.

Motion compensation using origin ensembles in awake small animal positron emission tomography

In emission tomographic imaging, the stochastic origin ensembles algorithm provides unique information regarding the detected counts given the measured data. Precision in both voxel and region-wise parameters may be determined for a single data set based on the posterior distribution of the count density allowing uncertainty estimates to be allocated to quantitative measures. Uncertainty estimates are of particular importance in awake animal neurological and behavioral studies for which head motion, unique for each acquired data set, perturbs the measured data. Motion compensation can be conducted when rigid head pose is measured during the scan. However, errors in pose measurements used for compensation can degrade the data and hence quantitative outcomes. In this investigation motion compensation and detector resolution models were incorporated into the basic origin ensembles algorithm and an efficient approach to computation was developed. The approach was validated against maximum liklihood-expectation maximisation and tested using simulated data. The resultant algorithm was then used to analyse quantitative uncertainty in regional activity estimates arising from changes in pose measurement precision. Finally, the posterior covariance acquired from a single data set was used to describe correlations between regions of interest providing information about pose measurement precision that may be useful in system analysis and design. The investigation demonstrates the use of origin ensembles as a powerful framework for evaluating statistical uncertainty of voxel and regional estimates. While in this investigation rigid motion was considered in the context of awake animal PET, the extension to arbitrary motion may provide clinical utility where respiratory or cardiac motion perturb the measured data.

Sensitivity estimation in time-of-flight list-mode positron emission tomography

PURPOSE

An accurate quantification of the images in positron emission tomography (PET) requires knowing the actual sensitivity at each voxel, which represents the probability that a positron emitted in that voxel is finally detected as a coincidence of two gamma rays in a pair of detectors in the PET scanner. This sensitivity depends on the characteristics of the acquisition, as it is affected by the attenuation of the annihilation gamma rays in the body, and possible variations of the sensitivity of the scanner detectors. In this work, the authors propose a new approach to handle time-of-flight (TOF) list-mode PET data, which allows performing either or both, a self-attenuation correction, and self-normalization correction based on emission data only.

METHODS

The authors derive the theory using a fully Bayesian statistical model of complete data. The authors perform an initial evaluation of algorithms derived from that theory and proposed in this work using numerical 2D list-mode simulations with different TOF resolutions and total number of detected coincidences. Effects of randoms and scatter are not simulated.

RESULTS

The authors found that proposed algorithms successfully correct for unknown attenuation and scanner normalization for simulated 2D list-mode TOF-PET data.

CONCLUSIONS

A new method is presented that can be used for corrections for attenuation and normalization (sensitivity) using TOF list-mode data.