Direct Estimation of Voxel-Wise Neurotransmitter Response Maps From Dynamic PET Data.

Computational methods, such as the linear parametric neurotransmitter PET (lp-ntPET) method, have been developed to characterize the transient changes in radiotracer kinetics in the target tissue during endogenous neurotransmitter release. In this paper, we describe and evaluate a parametric reconstruction algorithm that uses an expectation maximization framework, along with the lp-ntPET model, to estimate the endogenous neurotransmitter response to stimuli directly from the measured PET data. Computer simulations showed that the proposed direct reconstruction method offers improved accuracy and precision for the estimated timing parameters of the neurotransmitter response at the voxel level ( td=1±2 min, for activation onset bias and standard deviation) compared with conventional post reconstruction modeling ( td=4±7 min). In addition, we applied the proposed direct parameter estimation methodology to a [11C]raclopride displacement study of an awake rat and generated parametric maps illustrating the magnitude of ligand displacement from striatum. Although the estimated parametric maps of activation magnitude obtained from both direct and post reconstruction methodologies suffered from false positive activations, the proposed direct reconstruction framework offered more reliable parametric maps when the activation onset parameter was constrained.

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.

List-mode image reconstruction for positron emission tomography using tetrahedral voxels.

Image space decomposition based on tetrahedral voxels are interesting candidates for use in emission tomography. Tetrahedral voxels provide many of the advantages of point clouds with irregular spacing, such as being intrinsically multi-resolution, yet they also serve as a volumetric partition of the image space and so are comparable to more standard cubic voxels. Additionally, non-rigid displacement fields can be applied to the tetrahedral mesh in a straight-forward manner. So far studies incorporating tetrahedral decomposition of the image space have concentrated on pre-calculated, node-based, system matrix elements which reduces the flexibility of the tetrahedral approach and the capacity to accurately define regions of interest. Here, a list-mode on-the-fly calculation of the system matrix elements is described using a tetrahedral decomposition of the image space and volumetric elements-voxels. The algorithm is demonstrated in the context of awake animal PET which may require both rigid and non-rigid motion compensation, as well as quantification within small regions of the brain. This approach allows accurate, event based, motion compensation including non-rigid deformations.

Performance of MACACO Compton telescope for ion-beam therapy monitoring: first test with proton beams.

In order to exploit the advantages of ion-beam therapy in a clinical setting, delivery verification techniques are necessary to detect deviations from the planned treatment. Efforts are currently oriented towards the development of devices for real-time range monitoring. Among the different detector concepts proposed, Compton cameras are employed to detect prompt gammas and represent a valid candidate for real-time range verification. We present the first on-beam test of MACACO, a Compton telescope (multi-layer Compton camera) based on lanthanum bromide crystals and silicon photo-multipliers. The Compton telescope was first characterized through measurements and Monte Carlo simulations. The detector linearity was measured employing (22)Na and Am-Be sources, obtaining about 10% deviation from linearity at 3.44 MeV. A spectral image reconstruction algorithm was tested on synthetic data. Point-like sources emitting gamma rays with energy between 2 and 7 MeV were reconstructed with 3-5 mm resolution. The two-layer Compton telescope was employed to measure radiation emitted from a beam of 150 MeV protons impinging on a cylindrical PMMA target. Bragg-peak shifts were achieved via adjustment of the PMMA target location and the resulting measurements used during image reconstruction. Reconstructed Bragg peak profiles proved sufficient to observe peak-location differences within 10 mm demonstrating the potential of the MACACO Compton Telescope as a monitoring device for ion-beam therapy.

Sensitivity recovery for the AX-PET prototype using inter-crystal scattering events.

The development of novel detection devices and systems such as the AX-positron emission tomography (PET) demonstrator often introduce or increase the measurement of atypical coincidence events such as inter-crystal scattering (ICS). In more standard systems, ICS events often go undetected and the small measured fraction may be ignored. As the measured quantity of such events in the data increases, so too does the importance of considering them during image reconstruction. Generally, treatment of ICS events will attempt to determine which of the possible candidate lines of response (LoRs) correctly determine the annihilation photon trajectory. However, methods of assessment often have low success rates or are computationally demanding. In this investigation alternative approaches are considered. Experimental data was taken using the AX-PET prototype and a NEMA phantom. Three methods of ICS treatment were assessed--each of which considered all possible candidate LoRs during image reconstruction. Maximum likelihood expectation maximization was used in conjunction with both standard (line-like) and novel (V-like in this investigation) detection responses modeled within the system matrix. The investigation assumed that no information other than interaction locations was available to distinguish between candidates, yet the methods assessed all provided means by which such information could be included. In all cases it was shown that the signal to noise ratio is increased using ICS events. However, only one method, which used full modeling of the ICS response in the system matrix--the V-like model--provided enhancement in all figures of merit assessed in this investigation. Finally, the optimal method of ICS incorporation was demonstrated using data from two small animals measured using the AX-PET demonstrator.