Respiratory motion in emission tomography leads to reduced image quality. Developed correction methodology has been concentrating on the use of respiratory synchronized acquisitions leading to gated frames. Such frames, however, are of low signal-to-noise ratio as a result of containing reduced statistics. In this work, we describe the implementation of an elastic transformation within a list-mode-based reconstruction for the correction of respiratory motion over the thorax, allowing the use of all data available throughout a respiratory motion average acquisition. The developed algorithm was evaluated using datasets of the NCAT phantom generated at different points throughout the respiratory cycle. List-mode-data-based PET-simulated frames were subsequently produced by combining the NCAT datasets with Monte Carlo simulation. A non-rigid registration algorithm based on B-spline basis functions was employed to derive transformation parameters accounting for the respiratory motion using the NCAT dynamic CT images. The displacement matrices derived were subsequently applied during the image reconstruction of the original emission list mode data. Two different implementations for the incorporation of the elastic transformations within the one-pass list mode EM (OPL-EM) algorithm were developed and evaluated. The corrected images were compared with those produced using an affine transformation of list mode data prior to reconstruction, as well as with uncorrected respiratory motion average images. Results demonstrate that although both correction techniques considered lead to significant improvements in accounting for respiratory motion artefacts in the lung fields, the elastic-transformation-based correction leads to a more uniform improvement across the lungs for different lesion sizes and locations.

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.

To evaluate the carbon 11 ethanol kinetics with positron emission tomography after intratumoral injection of the tracer and assess its redistribution and dilution in patients who have hepatocellular carcinomas and who were scheduled for treatment with percutaneous ethanol injection. The study included eight patients with hepatocellular carcinomas. 11C ethanol was administered via a puncture needle positioned with ultrasonographic guidance. Parametric images based on the Fourier transformation were created for further analysis of the local distribution patterns of the tracer. The ratio of the 45-minute postinjection standardized uptake value to the 5-minute postinjection standardized uptake value was used for the evaluation of ethanol dilution. Five of eight tumors demonstrated almost constant uptake values after the initial distribution phase. In contrast, a rapid elimination of the 11C ethanol from the tumor was documented in three of eight tumors. The 45 minute-to-5 minute ratio was 0.18-0.67 (median value, 0.56) in the tumors. The time-activity curves of the normal liver parenchyma increased slowly but steadily with time owing to a low ethanol elimination from the tumor. Fourier transformation demonstrated inhomogeneous parts on the amplitude images in seven of eight tumors and random redistribution on the phase images in six of eight tumors. Inhomogeneous drug distribution and drug dilution in the target area are likely to be the major limiting parameters for therapy response.