Our goal in this paper is the estimation of kinetic model parameters for each voxel corresponding to a dense three-dimensional (3-D) positron emission tomography (PET) image. Typically, the activity images are first reconstructed from PET sinogram frames at each measurement time, and then the kinetic parameters are estimated by fitting a model to the reconstructed time-activity response of each voxel. However, this "indirect" approach to kinetic parameter estimation tends to reduce signal-to-noise ratio (SNR) because of the requirement that the sinogram data be divided into individual time frames. In 1985, Carson and Lange proposed, but did not implement, a method based on the expectation-maximization (EM) algorithm for direct parametric reconstruction. The approach is "direct" because it estimates the optimal kinetic parameters directly from the sinogram data, without an intermediate reconstruction step. However, direct voxel-wise parametric reconstruction remained a challenge due to the unsolved complexities of inversion and spatial regularization. In this paper, we demonstrate and evaluate a new and efficient method for direct voxel-wise reconstruction of kinetic parameter images using all frames of the PET data. The direct parametric image reconstruction is formulated in a Bayesian framework, and uses the parametric iterative coordinate descent (PICD) algorithm to solve the resulting optimization problem. The PICD algorithm is computationally efficient and is implemented with spatial regularization in the domain of the physiologically relevant parameters. Our experimental simulations of a rat head imaged in a working small animal scanner indicate that direct parametric reconstruction can substantially reduce root-mean-squared error (RMSE) in the estimation of kinetic parameters, as compared to indirect methods, without appreciably increasing computation.

Direct reconstruction of kinetic parameter images from dynamic PET data

Context: Expectations play a central role in the mechanism of the placebo effect. In Parkinson disease (PD), the placebo effect is associated with release of endogenous dopamine in both nigrostriatal and mesoaccumbens projections, yet the factors that control this dopamine release are undetermined. Objective: To determine how the strength of expectation of clinical improvement influences the degree of striatal dopamine release in response to placebo in patients with moderate PD. Design: Randomized, repeated-measures study with perceived expectation as the independent betweensubjects variable.

/pdf/effects-of-expectation-on-placebo-induced-dopamine-release-1xne1yvpha.pdf

Effects of expectation on placebo-induced dopamine release in Parkinson disease.

The HRRT PET system has the potential to produce human brain images with resolution better than 3 mm. To achieve the best possible accuracy and precision, we have designed MOLAR, a motion-compensation OSEM list-mode algorithm for resolution-recovery reconstruction on a computer cluster with the following features: direct use of list mode data with dynamic motion information (Polaris); exact reprojection of each line-of- response (LOR); system matrix computed from voxel-to-LOR distances (radial and axial); spatially varying resolution model implemented for each event by selection from precomputed line spread functions based on factors including detector obliqueness, crystal layer, and block detector position; distribution of events to processors and to subsets based on order of arrival; removal of voxels and events outside a reduced field-of-view defined by the attenuation map; no pre-corrections to Poisson data, i.e., all physical effects are defined in the model; randoms estimation from singles; model-based scatter simulation incorporated into the iterations; and component-based normalization. Preliminary computation estimates suggest that reconstruction of a single frame in one hour is achievable. Careful evaluation of this system will define which factors play an important role in producing high resolution, low-noise images with quantitative accuracy.

/pdf/design-of-a-motion-compensation-osem-list-mode-algorithm-for-mrxhxxpodd.pdf

Design of a motion-compensation OSEM list-mode algorithm for resolution-recovery reconstruction for the HRRT

PET and MRI provide complementary information in the study of the human brain. Simultaneous PET/MRI data acquisition allows the spatial and temporal correlation of the measured signals, creating opportunities impossible to realize using stand-alone instruments. This paper reviews the methodologic improvements and potential neurologic and psychiatric applications of this novel technology. We first present methods for improving the performance and information content of each modality by using the information provided by the other technique. On the PET side, we discuss methods that use the simultaneously acquired MRI data to improve the PET data quantification. On the MRI side, we present how improved PET quantification can be used to validate several MRI techniques. Finally, we describe promising research, translational, and clinical applications that can benefit from these advanced tools.

https://jnm.snmjournals.org/content/jnumed/53/12/1916.full.pdf

PET/MRI for Neurologic Applications

Head motion is difficult to avoid in long PET studies, degrading the image quality and offsetting the benefit of using a high-resolution scanner. As a potential solution in an integrated MR-PET scanner, the simultaneously acquired MRI data can be used for motion tracking. In this work, a novel algorithm for data processing and rigid-body motion correction (MC) for the MRI-compatible BrainPET prototype scanner is described, and proof-of-principle phantom and human studies are presented. Methods: To account for motion, the PET prompt and random coincidences and sensitivity data for postnormalization were processed in the line-of-response (LOR) space according to the MRI-derived motion estimates. The processing time on the standard BrainPET workstation is approximately 16 s for each motion estimate. After rebinning in the sinogram space, the motion corrected data were summed, and the PET volume was reconstructed using the attenuation and scatter sinograms in the reference position. The accuracy of the MC algorithm was first tested using a Hoffman phantom. Next, human volunteer studies were performed, and motion estimates were obtained using 2 high-temporal-resolution MRI-based motion-tracking techniques. Results: After accounting for the misalignment between the 2 scanners, perfectly coregistered MRI and PET volumes were reproducibly obtained. The MRI output gates inserted into the PET list-mode allow the temporal correlation of the 2 datasets within 0.2 ms. The Hoffman phantom volume reconstructed by processing the PET data in the LOR space was similar to the one obtained by processing the data using the standard methods and applying the MC in the image space, demonstrating the quantitative accuracy of the procedure. In human volunteer studies, motion estimates were obtained from echo planar imaging and cloverleaf navigator sequences every 3 s and 20 ms, respectively. Motion-deblurred PET images, with excellent delineation of specific brain structures, were obtained using these 2 MRI-based estimates. Conclusion: An MRI-based MC algorithm was implemented for an integrated MR-PET scanner. High-temporal-resolution MRI-derived motion estimates (obtained while simultaneously acquiring anatomic or functional MRI data) can be used for PET MC. An MRI-based MC method has the potential to improve PET image quality, increasing its reliability, reproducibility, and quantitative accuracy, and to benefit many neurologic applications.

https://jnm.snmjournals.org/content/jnumed/52/1/154.full.pdf

MRI-Assisted PET Motion Correction for Neurologic Studies in an Integrated MR-PET Scanner

A method is described to monitor the motion of the head during neurological positron emission tomography (PET) acquisitions and to correct the data post acquisition for the recorded motion prior to image reconstruction. The technique uses an optical tracking system, Polaris™, to accurately monitor the position of the head during the PET acquisition. The PET data are acquired in list mode where the events are written directly to disk during acquisition. The motion tracking information is aligned to the PET data using a sequence of pseudo-random numbers, which are inserted into the time tags in the list mode event stream through the gating input interface on the tomograph. The position of the head is monitored during the transmission acquisition, and it is assumed that there is minimal head motion during this measurement. Each event, prompt and delayed, in the list mode event stream is corrected for motion and transformed into the transmission space. For a given line of response, normalization, including corrections for detector efficiency, geometry and crystal interference and dead time are applied prior to motion correction and rebinning in the sinogram. A series of phantom experiments were performed to confirm the accuracy of the method: (a) a point source located in three discrete axial positions in the tomograph field of view, 0 mm, 10 mm and 20 mm from a reference point, (b) a multi-line source phantom rotated in both discrete and gradual rotations through ±5° and ±15°, including a vertical and horizontal movement in the plane. For both phantom experiments images were reconstructed for both the fixed and motion corrected data. Measurements for resolution, full width at half maximum (FWHM) and full width at tenth maximum (FWTM), were calculated from these images and a comparison made between the fixed and motion corrected datasets. From the point source measurements, the FWHM at each axial position was 7.1 mm in the horizontal direction, and increasing from 4.7 mm at the 0 mm position, to 4.8 mm, 20 mm offset, in the vertical direction. The results from the multi-line source phantom with ±5° rotations showed a maximum degradation in FWHM, when compared with the stationary phantom, of 0.6 mm, in the horizontal direction, and 0.3 mm in the vertical direction. The corresponding values for the larger rotation, ±15°, were 0.7 mm and 1.1 mm, respectively. The performance of the method was confirmed with a Hoffman brain phantom moved continuously, and a clinical acquisition using [11C]raclopride (normal volunteer). A visual comparison of both the motion and non-motion corrected images of the Hoffman brain phantom clearly demonstrated the efficacy of the method. A sample time-activity curve extracted from the clinical study showed irregularities prior to motion correction, which were removed after correction. A method has been developed to accurately monitor the motion of the head during a neurological PET acquisition, and correct for this motion prior to image reconstruction. The method has been demonstrated to be accurate and does not add significantly to either the acquisition or the subsequent data processing.

http://iopscience.iop.org/article/10.1088/0031-9155/48/8/301/pdf

Anthony M Westrip

Papers

The design and implementation of a motion correction scheme for neurological PET.