List-mode OSEM Research Articles

Improvement of double photon emission Compton imaging using angular and polarization correlations in cascade photons

A Compton camera is a gamma ray imaging system that produces images by superposing Compton cones calculated from Compton kinematics. However, the reconstructed images exhibit a high background. In this study, the angular correlation and polarization correlation were used to decrease the background. 111In, which emits 171 keV and 245 keV cascade gamma rays, was imaged using Compton camera modules in the ring structure. The 3D list-mode ordered subset expectation maximization method was used to reconstruct the images. Probability modeling of events reflecting angular and polarization correlations was applied to the image-reconstruction algorithm. The results were compared with those of a typical Compton camera method using only 245 keV Compton cones. The sensitivity of the proposed method was 1.03% of that of the typical method. However, compared to the typical method, the application of angular and polarization correlations improved the signal-to-background ratio. The proposed method is expected to be useful in capturing low-background images wherein a ring-type Compton camera and cascade photon-emitting nuclides can be used.

Super-resolution in brain positron emission tomography using a real-time motion capture system

Super-resolution (SR) is a methodology that seeks to improve image resolution by exploiting the increased spatial sampling information obtained from multiple acquisitions of the same target with accurately known sub-resolution shifts. This work aims to develop and evaluate an SR estimation framework for brain positron emission tomography (PET), taking advantage of a high-resolution infra-red tracking camera to measure shifts precisely and continuously. Moving phantoms and non-human primate (NHP) experiments were performed on a GE Discovery MI PET/CT scanner (GE Healthcare) using an NDI Polaris Vega (Northern Digital Inc), an external optical motion tracking device. To enable SR, a robust temporal and spatial calibration of the two devices was developed as well as a list-mode Ordered Subset Expectation Maximization PET reconstruction algorithm, incorporating the high-resolution tracking data from the Polaris Vega to correct motion for measured line of responses on an event-by-event basis. For both phantoms and NHP studies, the SR reconstruction method yielded PET images with visibly increased spatial resolution compared to standard static acquisitions, allowing improved visualization of small structures. Quantitative analysis in terms of SSIM, CNR and line profiles were conducted and validated our observations. The results demonstrate that SR can be achieved in brain PET by measuring target motion in real-time using a high-resolution infrared tracking camera.

3D Compton image reconstruction method for whole gamma imaging

Compton imaging or Compton camera imaging has been studied well, but its advantages in nuclear medicine and molecular imaging have not been demonstrated yet. Therefore, the aim of this work was to compare Compton imaging with positron emission tomography (PET) by using the same imaging platform of whole gamma imaging (WGI). WGI is a concept that combines PET with Compton imaging by inserting a scatterer ring into a PET ring. This concept utilizes diverse types of gamma rays for 3D tomographic imaging. In this paper, we remodeled our previous WGI prototype for small animal imaging, and we developed an image reconstruction method based on a list-mode ordered subset expectation maximization algorithm incorporating detector response function modeling, random correction and normalization (sensitivity correction) for either PET and Compton imaging. To the best of our knowledge, this is the world’s first realization of a full-ring Compton imaging system. We selected 89Zr as an imaging target because a 89Zr nuclide emits a 909 keV single-gamma ray as well as a positron, and we can directly compare Compton imaging of 909 keV photons with PET, a well-established modality. We measured a cylindrical phantom and a small rod phantom filled with 89Zr solutions of 10.3 MBq and 10.2 MBq activity, respectively, for 1 h each. The uniform radioactivity distribution of the cylindrical phantom was reconstructed with normalization in both PET and Compton imaging. Coefficients of variation for region-of-interest values were 4.2% for Compton imaging and 3.3% for PET; the difference might be explained by the difference in the detected count number. The small rod phantom experiment showed that the WGI Compton imaging had spatial resolution better than 3.0 mm at the peripheral region although the center region had lower resolution. PET resolved 2.2 mm rods clearly at any location. We measured a mouse for 1 h, 1 d after injection of 9.8 MBq 89Zr oxalate. The 89Zr assimilated in the mouse bony structures was clearly depicted, and Compton imaging results agreed well with PET images, especially for the region inside the scatterer ring. In conclusion, we demonstrated the performance of WGI using the developed Compton image reconstruction method. We realized Compton imaging with a quality approaching that of PET, which is supporting a future expectation that Compton imaging outperforms PET.

Whole gamma imaging: a new concept of PET combined with Compton imaging

We proposed a concept of whole gamma imaging (WGI) that utilizes all detectable gamma rays for imaging. An additional detector ring, which is used as the scatterer, is inserted in the field-of-view of a PET ring so that single gamma rays can be detected by the Compton imaging method. In particular, for the non-pure positron emitters which emit an additional gamma ray almost at the same time, triple gamma imaging will be enabled; localization on each line-of-response (LOR) is possible by using the Compton cone of the additional gamma ray. We developed a prototype to show a proof of the WGI concept. The diameters of scatterer ring and PET ring were set as 20 cm and 66 cm, respectively. For Compton imaging of the 662-keV gamma ray from a 137Cs point source, spatial resolution obtained by the list-mode OSEM algorithm was 4.4 mm FWHM at the 8 cm off-center position and 13.1 mm FWHM at the center position. For PET imaging of a 22Na point source, spatial resolution was about 2 mm FWHM at all positions. For the triple gamma imaging, 5.7 mm FWHM (center) and 4.8 mm FWHM (8 cm off-center) were obtained for the 22Na point source just by plotting the intersecting points between each LOR and each Compton cone of the 1275-keV gamma ray. No image reconstruction was applied. Scandium-44 was produced as a practical candidate of the non-pure positron emitters, and 6.6 mm FWHM (center) and 5.8 mm FWHM (8 cm off-center) were obtained in the same manner. This direct imaging approach which neither requires time-consuming event integration nor iterative image reconstruction may allow in vivo real-time tracking of a tiny amount of activity. Our initial results showed the feasibility of the WGI concept, which is a novel combination of PET and Compton imaging.

Total-Body PET and Highly Stable Chelators Together Enable Meaningful 89Zr-Antibody PET Studies up to 30 Days After Injection.

The use of 89Zr-antibody PET imaging to measure antibody biodistribution and tissue pharmacokinetics is well established, but current PET systems lack the sensitivity needed to study 89Zr-labeled antibodies beyond 2-3 isotope half-lives (7-10 d), after which a poor signal-to-noise ratio is problematic. However, studies across many weeks are desirable to better match antibody circulation half-life in human and nonhuman primates. These studies investigated the technical feasibility of using the primate mini-EXPLORER PET scanner, making use of its high sensitivity and 45-cm axial field of view, for total-body imaging of 89Zr-labeled antibodies in rhesus monkeys up to 30 d after injection. Methods: A humanized monoclonal IgG antibody against the herpes simplex viral protein glycoprotein D (gD) was radiolabeled with 89Zr via 1 of 4 chelator-linker combinations (benzyl isothiocyanate-DFO [DFO-Bz-NCS], where DFO is desferrioxamine B; DFO-squaramide; DFO*-Bz-NCS, where DFO* is desferrioxamine*; and DFO*-squaramide). The pharmacokinetics associated with these 4 chelator-linker combinations were compared in 12 healthy young male rhesus monkeys (∼1-2 y old, ∼3 ± 1 kg). Each animal was initially injected intravenously with unlabeled antibody in a peripheral vessel in the right arm (10 mg/kg, providing therapeutic-level antibody concentrations), immediately followed by approximately 40 MBq of one of the 89Zr-labeled antibodies injected intravenously in a peripheral vessel in the left arm. All animals were imaged 6 times over a period of 30 d, with an initial 60-min dynamic scan on day 0 (day of injection) followed by static scans of 30-45 min on approximately days 3, 7, 14, 21, and 30, with all acquired using a single bed position and images reconstructed using time-of-flight list-mode ordered-subsets expectation maximization. Activity concentrations in various organs were extracted from the PET images using manually defined regions of interest. Results: Excellent image quality was obtained, capturing the initial distribution phase in the whole-body scan; later time points showed residual 89Zr mainly in the liver. Even at 30 d after injection, representing approximately 9 half-lives of 89Zr and with a total residual activity of only 20-40 kBq in the animal, the image quality was sufficient to readily identify activity in the liver, kidneys, and upper and lower limb joints. Significant differences were noted in late time point liver uptake, bone uptake, and whole-body clearance between chelator-linker types, whereas little variation (±10%) was observed within each type. Conclusion: These studies demonstrate the ability to image 89Zr-radiolabeled antibodies up to 30 d after injection while maintaining satisfactory image quality, as provided by the primate mini-EXPLORER with high sensitivity and long axial field of view. Quantification demonstrated potentially important differences in the behavior of the 4 chelators. This finding supports further investigation.

Total-Body Dynamic Reconstruction and Parametric Imaging on the uEXPLORER.

The world's first 194-cm-long total-body PET/CT scanner (uEXPLORER) has been built by the EXPLORER Consortium to offer a transformative platform for human molecular imaging in clinical research and health care. Its total-body coverage and ultra-high sensitivity provide opportunities for more accurate tracer kinetic analysis in studies of physiology, biochemistry, and pharmacology. The objective of this study was to demonstrate the capability of total-body parametric imaging and to quantify the improvement in image quality and kinetic parameter estimation by direct and kernel reconstruction of the uEXPLORER data. Methods: We developed quantitative parametric image reconstruction methods for kinetic analysis and used them to analyze the first human dynamic total-body PET study. A healthy female subject was recruited, and a 1-h dynamic scan was acquired during and after an intravenous injection of 256 MBq of 18F-FDG. Dynamic data were reconstructed using a 3-dimensional time-of-flight list-mode ordered-subsets expectation maximization (OSEM) algorithm and a kernel-based algorithm with all quantitative corrections implemented in the forward model. The Patlak graphical model was used to analyze the 18F-FDG kinetics in the whole body. The input function was extracted from a region over the descending aorta. For comparison, indirect Patlak analysis from reconstructed frames and direct reconstruction of parametric images from the list-mode data were obtained for the last 30 min of data. Results: Images reconstructed by OSEM showed good quality with low noise, even for the 1-s frames. The image quality was further improved using the kernel method. Total-body Patlak parametric images were obtained using either indirect estimation or direct reconstruction. The direct reconstruction method improved the parametric image quality, having a better contrast-versus-noise tradeoff than the indirect method, with a 2- to 3-fold variance reduction. The kernel-based indirect Patlak method offered image quality similar to the direct Patlak method, with less computation time and faster convergence. Conclusion: This study demonstrated the capability of total-body parametric imaging using the uEXPLORER. Furthermore, the results showed the benefits of kernel-regularized reconstruction and direct parametric reconstruction. Both can achieve superior image quality for tracer kinetic studies compared with the conventional indirect OSEM for total-body imaging.

PSF reconstruction for Compton-based prompt gamma imaging

Compton-based prompt gamma (PG) imaging has been proposed for in vivo range verification in proton therapy. However, several factors degrade the image quality of PG images, some of which are due to inherent properties of a Compton camera such as spatial resolution and energy resolution. Moreover, Compton-based PG imaging has a spatially variant resolution loss. In this study, we investigate the performance of the list-mode ordered subset expectation maximization algorithm with a shift-variant point spread function (LM-OSEM-SV-PSF) model. We also evaluate how well the PG images reconstructed using an SV-PSF model reproduce the distal falloff of the proton beam. The SV-PSF parameters were estimated from simulation data of point sources at various positions. Simulated PGs were produced in a water phantom irradiated with a proton beam. Compared to the LM-OSEM algorithm, the LM-OSEM-SV-PSF algorithm improved the quality of the reconstructed PG images and the estimation of PG falloff positions. In addition, the 4.44 and 5.25 MeV PG emissions can be accurately reconstructed using the LM-OSEM-SV-PSF algorithm. However, for the 2.31 and 6.13 MeV PG emissions, the LM-OSEM-SV-PSF reconstruction provides limited improvement. We also found that the LM-OSEM algorithm followed by a shift-variant Richardson–Lucy deconvolution could reconstruct images with quality visually similar to the LM-OSEM-SV-PSF-reconstructed images, while requiring shorter computation time.

PET Image Reconstruction and Deformable Motion Correction Using Unorganized Point Clouds.

Quantitative positron emission tomography imaging often requires correcting the image data for deformable motion. With cyclic motion, this is traditionally achieved by separating the coincidence data into a relatively small number of gates, and incorporating the inter-gate image transformation matrices into the reconstruction algorithm. In the presence of non-cyclic deformable motion, this approach may be impractical due to a large number of required gates. In this paper, we propose an alternative approach to iterative image reconstruction with correction for deformable motion, wherein unorganized point clouds are used to model the imaged objects in the image space, and motion is corrected for explicitly by introducing a time-dependence into the point coordinates. The image function is represented using constant basis functions with finite support determined by the boundaries of the Voronoi cells in the point cloud. We validate the quantitative accuracy and stability of the proposed approach by reconstructing noise-free and noisy projection data from digital and physical phantoms. The point-cloud-based maximum likelihood expectation maximization (MLEM) and one-pass list-mode ordered-subset expectation maximization (OSEM) algorithms are validated. The results demonstrate that images reconstructed using the proposed method are quantitatively stable, with noise and convergence properties comparable to image reconstruction based on the use of rectangular and radially-symmetric basis functions.

Impact of detector design on imaging performance of a long axial field-of-view, whole-body PET scanner

Current generation of commercial time-of-flight (TOF) PET scanners utilize 20–25 mm thick LSO or LYSO crystals and have an axial FOV (AFOV) in the range of 16–22 mm. Longer AFOV scanners would provide increased intrinsic sensitivity and require fewer bed positions for whole-body imaging. Recent simulation work has investigated the sensitivity gains that can be achieved with these long AFOV scanners, and has motivated new areas of investigation such as imaging with a very low dose of injected activity as well as providing whole-body dynamic imaging capability in one bed position. In this simulation work we model a 72 cm long scanner and prioritize the detector design choices in terms of timing resolution, crystal size (spatial resolution), crystal thickness (detector sensitivity), and depth-of-interaction (DOI) measurement capability. The generated list data are reconstructed with a list-mode OSEM algorithm using a Gaussian TOF kernel that depends on the timing resolution and blob basis functions for regularization. We use lesion phantoms and clinically relevant metrics for lesion detectability and contrast measurement. The scan time was fixed at 10 min for imaging a 100 cm long object assuming a 50% overlap between adjacent bed positions. Results show that a 72 cm long scanner can provide a factor of ten reduction in injected activity compared to an identical 18 cm long scanner to get equivalent lesion detectability. While improved timing resolution leads to further gains, using 3 mm (as opposed to 4 mm) wide crystals does not show any significant benefits for lesion detectability. A detector providing 2-level DOI information with equal crystal thickness also does not show significant gains. Finally, a 15 mm thick crystal leads to lower lesion detectability than a 20 mm thick crystal when keeping all other detector parameters (crystal width, timing resolution, and DOI capability) the same. However, improved timing performance with 15 mm thick crystals can provide similar or better performance than that achieved by a detector using 20 mm thick crystals.

A novel partial volume effects correction technique integrating deconvolution associated with denoising within an iterative PET image reconstruction.

Partial volume effect (PVE) plays an important role in both qualitative and quantitative PET image accuracy, especially for small structures. A previously proposed voxelwise PVE correction method applied on PET reconstructed images involves the use of Lucy-Richardson deconvolution incorporating wavelet-based denoising to limit the associated propagation of noise. The aim of this study is to incorporate the deconvolution, coupled with the denoising step, directly inside the iterative reconstruction process to further improve PVE correction. The list-mode ordered subset expectation maximization (OSEM) algorithm has been modified accordingly with the application of the Lucy-Richardson deconvolution algorithm to the current estimation of the image, at each reconstruction iteration. Acquisitions of the NEMA NU2-2001 IQ phantom were performed on a GE DRX PET/CT system to study the impact of incorporating the deconvolution inside the reconstruction [with and without the point spread function (PSF) model] in comparison to its application postreconstruction and to standard iterative reconstruction incorporating the PSF model. The impact of the denoising step was also evaluated. Images were semiquantitatively assessed by studying the trade-off between the intensity recovery and the noise level in the background estimated as relative standard deviation. Qualitative assessments of the developed methods were additionally performed on clinical cases. Incorporating the deconvolution without denoising within the reconstruction achieved superior intensity recovery in comparison to both standard OSEM reconstruction integrating a PSF model and application of the deconvolution algorithm in a postreconstruction process. The addition of the denoising step permitted to limit the SNR degradation while preserving the intensity recovery. This study demonstrates the feasibility of incorporating the Lucy-Richardson deconvolution associated with a wavelet-based denoising in the reconstruction process to better correct for PVE. Future work includes further evaluations of the proposed method on clinical datasets and the use of improved PSF models.

Multi-tracer Imaging of a Compton Camera

Since a Compton camera has high detection sensitivity due to electronic collimation and a good energy resolution, it is a potential imaging system for nuclear medicine. In this study, we investigated the feasibility of a Compton camera for multi-tracer imaging and proposed a rotating Compton camera to satisfy Orlov’s condition for 3D imaging. Two software phantoms of 140 and 511 keV radiation sources were used for Monte-Carlo simulation and then the simulation data were reconstructed by listmode ordered subset expectation maximization to evaluate the capability of multi-tracer imaging in a Compton camera. And the Compton camera rotating around the object was proposed and tested with different rotation angle steps for improving the limited coverage of the fixed conventional Compton camera over the field-of-view in terms of histogram of angles in spherical coordinates. The simulation data showed the separate 140 and 511 keV images from simultaneous multi-tracer detection in both 2D and 3D imaging and the number of valid projection li nes on the conical surfaces was inversely proportional to the decrease of rotation angle. Considering computation load and proper number of projection lines on the conical surface, the rotation angle of 30 degree was sufficient for 3D imaging of the Compton camera in terms of 26 min of computation time and 5 million of detected event number and the increased detection time can be solved with multiple Compton camera system. The Compton camera proposed in this study can be effective system for multi-tracer imaging and is a potential system for development of various disease diagnosis and therapy approaches.

Non rigid respiratory motion correction in whole body PET/MR imaging

Respiratory motion in PET/MR imaging leads to reduced quantitative and qualitative image accuracy. Correction methodologies include the use of respiratory synchronized gated frames which lead to low signal to noise ratio (SNR) given that each frame contains only part of the count available throughout an average PET acquisition. In this work, 4D MRI extracted elastic transformations were applied to list-mode data either inside the image reconstruction or to the reconstructed respiratory synchronized images to obtain respiration corrected PET images. Five patient datasets acquired on the SIEMENS mMR PET/MR system were used. T1-weighted 4D MR images were registered to the end expiration MR image using an elastic B-spline registration algorithm to derive deformation matrices used to account for the respiratory motion. The derived matrices were subsequently applied during the List Mode OSEM PET image reconstruction of the original emission list-mode data. The corrected images were compared with those produced by applying the elastic transformation after reconstructing PET frames followed by summing the realigned gated frames, as well as with uncorrected motion averaged images. Significant improvement was obtained by using any of the two respiratory motion correction techniques compared to uncorrected motion average images (>30% differences in tumor position and size and 25% SNR increase). Both correction approaches lead to nearly equivalent improvements (differences of <8% and <6% for the lesion FWHM and the SNR respectively). A list-mode reconstruction based respiratory motion correction for PET has been implemented and its performance evaluated. This approach was based on the use of elastic transformations derived from 4D MRI during PET image reconstrucion. Our results show significant respiratory motion compensation when compared to the motion average PET images, with slightly better contrast and SNR compared to an equivalent 4D PET image space elastic motion correction method.

Optimization, evaluation, and comparison of standard algorithms for image reconstruction with the VIP-PET

A novel positron emission tomography (PET) scanner design based on a room-temperature pixelated CdTe solid-state detector is being developed within the framework of the Voxel Imaging PET (VIP) Pathfinder project [1]. The simulation results show a great potential of the VIP to produce high-resolution images even in extremely challenging conditions such as the screening of a human head [2]. With unprecedented high channel density (450 channels/cm3) image reconstruction is a challenge. Therefore optimization is needed to find the best algorithm in order to exploit correctly the promising detector potential. The following reconstruction algorithms are evaluated: 2-D Filtered Backprojection (FBP), Ordered Subset Expectation Maximization (OSEM), List-Mode OSEM (LM-OSEM), and the Origin Ensemble (OE) algorithm. The evaluation is based on the comparison of a true image phantom with a set of reconstructed images obtained by each algorithm. This is achieved by calculation of image quality merit parameters such as the bias, the variance and the mean square error (MSE). A systematic optimization of each algorithm is performed by varying thereconstruction parameters, such as the cutoff frequency of the noise filters and the number of iterations. The region of interest (ROI) analysis of the reconstructed phantom is also performed for each algorithm and the results are compared. Additionally, the performance of the image reconstruction methods is compared by calculating the modulation transfer function (MTF). The reconstruction time is also taken into account to choose the optimal algorithm. The analysis is based on GAMOS [3] simulation including the expected CdTe and electronic specifics.

Evaluation of list-mode ordered subset expectation maximization image reconstruction for pixelated solid-state compton gamma camera with large number of channels

The Voxel Imaging PET (VIP) Pathfinder project intends to show the advantages of using pixelated solid-state technology for nuclear medicine applications. It proposes designs for Positron Emission Tomography (PET), Positron Emission Mammography (PEM) and Compton gamma camera detectors with a large number of signal channels (of the order of 106). For Compton camera, especially with a large number of readout channels, image reconstruction presents a big challenge. In this work, results are presented for the List-Mode Ordered Subset Expectation Maximization (LM-OSEM) image reconstruction algorithm on simulated data with the VIP Compton camera design. For the simulation, all realistic contributions to the spatial resolution are taken into account, including the Doppler broadening effect. The results show that even with a straightforward implementation of LM-OSEM, good images can be obtained for the proposed Compton camera design. Results are shown for various phantoms, including extended sources and with a distance between the field of view and the first detector plane equal to 100 mm which corresponds to a realistic nuclear medicine environment.

Development of a prototype PET scanner with depth-of-interaction measurement using solid-state photomultiplier arrays and parallel readout electronics

In this study, we developed a prototype animal PET by applying several novel technologies to use solid-state photomultiplier (SSPM) arrays to measure the depth of interaction (DOI) and improve imaging performance. Each PET detector has an 8 × 8 array of about 1.9 × 1.9 × 30.0 mm3 lutetium-yttrium-oxyorthosilicate scintillators, with each end optically connected to an SSPM array (16 channels in a 4 × 4 matrix) through a light guide to enable continuous DOI measurement. Each SSPM has an active area of about 3 × 3 mm2, and its output is read by a custom-developed application-specific integrated circuit to directly convert analogue signals to digital timing pulses that encode the interaction information. These pulses are transferred to and are decoded by a field-programmable gate array-based time-to-digital convertor for coincident event selection and data acquisition. The independent readout of each SSPM and the parallel signal process can significantly improve the signal-to-noise ratio and enable the use of flexible algorithms for different data processes. The prototype PET consists of two rotating detector panels on a portable gantry with four detectors in each panel to provide 16 mm axial and variable transaxial field-of-view (FOV) sizes. List-mode ordered subset expectation maximization image reconstruction was implemented. The measured mean energy, coincidence timing and DOI resolution for a crystal were about 17.6%, 2.8 ns and 5.6 mm, respectively. The measured transaxial resolutions at the center of the FOV were 2.0 mm and 2.3 mm for images reconstructed with and without DOI, respectively. In addition, the resolutions across the FOV with DOI were substantially better than those without DOI. The quality of PET images of both a hot-rod phantom and mouse acquired with DOI was much higher than that of images obtained without DOI. This study demonstrates that SSPM arrays and advanced readout/processing electronics can be used to develop a practical DOI-measureable PET scanner.

Conceptual Design and Simulation Study of an ROI-Focused Panel-PET Scanner

Positron emission tomography (PET) is an important imaging modality for clincial use. Conventionally, the PET scanner is generally built to provide a roomy enough transverse field-of-view (FOV) for imaging most adults’ torsos. However, in many cases, the region-of-interest (ROI) for imaging is usually a small area inside the human body. Therefore, to fulfill a PET system which provides an FOV comparable in size to the target ROI seems appealing and more cost effective. Meanwhile, such a PET system has the potential for portable or bedside application with the reduced system size. In this work, we have investigated the feasibility of using dual-headed panel-detectors to build an ROI-focused PET scanner. A novel windowed list-mode ordered subset expectation maximization method was developed to perform the ROI image reconstruction. With this method, the ROI of the object can be reconstructed from the coincidences whose position determined by time-of-flight (TOF) measurements was inside the ROI. Monte Carlo simulation demonstrates the feasibility of detecting lesions not less than 1 cm in diameter, with a 300 ps full width at half maximum timing resolution. As a critical system performance, the impact of TOF information on image quality has been studied and the required TOF capability was assessed. With enhanced timing resolution, the distortions and artifacts were reduced effectively. The further improved TOF capability also shows a noticeable improvement of detection performance for low uptake lesions, as well as the recovery speed of lesion contrast, which is of practical significance in the lesion detection task.

Resolution recovery reconstruction for a Compton camera

The spatial resolution from Compton cameras suffers from measurement uncertainties in interaction positions and energies. The degree of degradation in spatial resolution is shift-variant (SV) over the field-of-view (FOV) because the imaging principle is based on the conical surface integration. In our study, the shift-variant point spread function (SV-PSF) is derived from point source measurements at various positions in the FOV and is incorporated into the system matrix of a fully three-dimensional, accelerated reconstruction, i.e. the listmode ordered subset expectation maximization (LMOSEM) algorithm, for resolution recovery. Simulation data from point sources were used to estimate SV and asymmetric parameters for Gaussian, Cauchy, and general parametric PSFs. Although little difference in the fitness accuracy between Gaussian and general parametric PSFs was observed, the general parametric model showed greater flexibility over the FOV in shaping the curve between that for Gaussian and Cauchy functions. The estimated asymmetric SV-PSFs were incorporated into the LMOSEM for resolution recovery. For simulation data from a single point source at the origin, all LMOSEM-SV-PSFs improved the spatial resolution by 2.6 times over the standard LMOSEM. For two point-source simulations, reconstructions also gave a two-fold improvement in spatial resolution and resulted in a greater recovered activity ratio at different positions in the FOV.

Massively parallelizable list-mode reconstruction using a Monte Carlo-based elliptical Gaussian model

A fully three-dimensional (3D) massively parallelizable list-mode ordered-subsets expectation-maximization (LM-OSEM) reconstruction algorithm has been developed for high-resolution PET cameras. System response probabilities are calculated online from a set of parameters derived from Monte Carlo simulations. The shape of a system response for a given line of response (LOR) has been shown to be asymmetrical around the LOR. This work has been focused on the development of efficient region-search techniques to sample the system response probabilities, which are suitable for asymmetric kernel models, including elliptical Gaussian models that allow for high accuracy and high parallelization efficiency. The novel region-search scheme using variable kernel models is applied in the proposed PET reconstruction algorithm. A novel region-search technique has been used to sample the probability density function in correspondence with a small dynamic subset of the field of view that constitutes the region of response (ROR). The ROR is identified around the LOR by searching for any voxel within a dynamically calculated contour. The contour condition is currently defined as a fixed threshold over the posterior probability, and arbitrary kernel models can be applied using a numerical approach. The processing of the LORs is distributed in batches among the available computing devices, then, individual LORs are processed within different processing units. In this way, both multicore and multiple many-core processing units can be efficiently exploited. Tests have been conducted with probability models that take into account the noncolinearity, positron range, and crystal penetration effects, that produced tubes of response with varying elliptical sections whose axes were a function of the crystal's thickness and angle of incidence of the given LOR. The algorithm treats the probability model as a 3D scalar field defined within a reference system aligned with the ideal LOR. This new technique provides superior image quality in terms of signal-to-noise ratio as compared with the histogram-mode method based on precomputed system matrices available for a commercial small animal scanner. Reconstruction times can be kept low with the use of multicore, many-core architectures, including multiple graphic processing units. A highly parallelizable LM reconstruction method has been proposed based on Monte Carlo simulations and new parallelization techniques aimed at improving the reconstruction speed and the image signal-to-noise of a given OSEM algorithm. The method has been validated using simulated and real phantoms. A special advantage of the new method is the possibility of defining dynamically the cut-off threshold over the calculated probabilities thus allowing for a direct control on the trade-off between speed and quality during the reconstruction.

Fast GPU-based computation of the sensitivity matrix for a PET list-mode OSEM algorithm

During the last decade, studies have shown that 3D list-mode ordered-subset expectation-maximization (LM-OSEM) algorithms for positron emission tomography (PET) reconstruction could be effectively computed and considerably accelerated by graphics processing unit (GPU) devices. However, most of these studies rely on pre-calculated sensitivity matrices. In many cases, the time required to compute this matrix can be longer than the reconstruction time itself. In fact, the relatively long time required for the calculation of the patient-specific sensitivity matrix is considered as one of the main obstacle in introducing a list-mode PET reconstruction algorithm for routine clinical use. The objective of this work is to accelerate a fully 3D LM-OSEM algorithm, including the calculation of the sensitivity matrix that accounts for the patient-specific attenuation and normalization corrections. For this purpose, sensitivity matrix calculations and list-mode OSEM reconstructions were implemented on GPUs, using the geometry of a commercial PET system. The system matrices were built on-the-fly by using an approach with multiple rays per detector pair. The reconstructions were performed for a volume of 188×188×57 voxels of 2×2×3.15 mm3 and for another volume of 144×144×57 voxels of 4×4×3.15 mm3. The time to compute the sensitivity matrix for the 188×188×57 array was 9 s while the LM-OSEM algorithm performed at a rate of 1.1 millions of events per second. For the 144×144×57 array, the respective numbers are 8 s for the sensitivity matrix and 0.8 million of events per second for the LM-OSEM step. This work lets envision fast reconstructions for advanced PET applications such as real time dynamic studies and parametric image reconstructions.

MRI-Based Nonrigid Motion Correction in Simultaneous PET/MRI

Respiratory and cardiac motion is the most serious limitation to whole-body PET, resulting in spatial resolution close to 1 cm. Furthermore, motion-induced inconsistencies in the attenuation measurements often lead to significant artifacts in the reconstructed images. Gating can remove motion artifacts at the cost of increased noise. This paper presents an approach to respiratory motion correction using simultaneous PET/MRI to demonstrate initial results in phantoms, rabbits, and nonhuman primates and discusses the prospects for clinical application. Studies with a deformable phantom, a free-breathing primate, and rabbits implanted with radioactive beads were performed with simultaneous PET/MRI. Motion fields were estimated from concurrently acquired tagged MR images using 2 B-spline nonrigid image registration methods and incorporated into a PET list-mode ordered-subsets expectation maximization algorithm. Using the measured motion fields to transform both the emission data and the attenuation data, we could use all the coincidence data to reconstruct any phase of the respiratory cycle. We compared the resulting SNR and the channelized Hotelling observer (CHO) detection signal-to-noise ratio (SNR) in the motion-corrected reconstruction with the results obtained from standard gating and uncorrected studies. Motion correction virtually eliminated motion blur without reducing SNR, yielding images with SNR comparable to those obtained by gating with 5-8 times longer acquisitions in all studies. The CHO study in dynamic phantoms demonstrated a significant improvement (166%-276%) in lesion detection SNR with MRI-based motion correction as compared with gating (P < 0.001). This improvement was 43%-92% for large motion compared with lesion detection without motion correction (P < 0.001). CHO SNR in the rabbit studies confirmed these results. Tagged MRI motion correction in simultaneous PET/MRI significantly improves lesion detection compared with respiratory gating and no motion correction while reducing radiation dose. In vivo primate and rabbit studies confirmed the improvement in PET image quality and provide the rationale for evaluation in simultaneous whole-body PET/MRI clinical studies.