Dynamic PET reconstruction using wavelet regularization with adapted basis functions

Memory-efficient model-based deep learning with convergence and robustness guarantees

Computational imaging has been revolutionized by compressed sensing algorithms, which offer guaranteed uniqueness, convergence, and stability properties. Model-based deep learning methods that combine imaging physics with learned regularization priors have emerged as more powerful alternatives for image recovery. The main focus of this paper is to introduce a memory efficient model-based algorithm with similar theoretical guarantees as CS methods. The proposed iterative algorithm alternates between a gradient descent involving the score function and a conjugate gradient algorithm to encourage data consistency. The score function is modeled as a monotone convolutional neural network. Our analysis shows that the monotone constraint is necessary and sufficient to enforce the uniqueness of the fixed point in arbitrary inverse problems. In addition, it also guarantees the convergence to a fixed point, which is robust to input perturbations. We introduce two implementations of the proposed MOL framework, which differ in the way the monotone property is imposed. The first approach enforces a strict monotone constraint, while the second one relies on an approximation. The guarantees are not valid for the second approach in the strict sense. However, our empirical studies show that the convergence and robustness of both approaches are comparable, while the less constrained approximate implementation offers better performance. The proposed deep equilibrium formulation is significantly more memory efficient than unrolled methods, which allows us to apply it to 3D or 2D+time problems that current unrolled algorithms cannot handle.

Denoising non-steady state dynamic PET data using a feed-forward neural network

The quality of reconstructed dynamic PET images, as well as the statistical reliability of the estimated pharmacokinetic parameters is often compromised by high levels of statistical noise, particularly at the voxel level. Many denoising strategies have been proposed, both in the temporal and spatial domain, which substantially improve the signal to noise ratio of the reconstructed dynamic images. However, although most filtering approaches are fairly successful in reducing the spatio-temporal inter-voxel variability, they may also average out or completely eradicate the critically important temporal signature of a transient neurotransmitter activation response that may be present in a non-steady state dynamic PET study. In this work, we explore an approach towards temporal denoising of non-steady state dynamic PET images using an artificial neural network, which was trained to identify the temporal profile of a time-activity curve, while preserving any potential activation response. We evaluated the performance of a feed-forward perceptron neural network to improve the signal to noise ratio of dynamic [¹¹C]raclopride activation studies and compared it with the widely used highly constrained back projection (HYPR) filter. Results on both simulated Geant4 Application for Tomographic Emission data of a realistic rat brain phantom and experimental animal data of a freely moving animal study showed that the proposed neural network can efficiently improve the noise characteristics of dynamic data in the temporal domain, while it can lead to a more reliable estimation of voxel-wise activation response in target region. In addition, improvements in signal-to-noise ratio achieved by denoising the dynamic data using the proposed neural network led to improved accuracy and precision of the estimated model parameters of the lp-ntPET model, compared to the HYPR filter. The performance of the proposed denoising approach strongly depends on the amount of noise in the dynamic PET data, with higher noise leading to substantially higher variability in the estimated parameters of the activation response. Overall, the feed-forward network led to a similar performance as the HYPR filter in terms of spatial denoising, but led to notable improvements in terms of temporal denoising, which in turn improved the estimation activation parameters.

Semi-Linearized Proximal Alternating Minimization for a Discrete Mumford-Shah Model

The Mumford-Shah model is a standard model in image segmentation, and due to its difficulty, many approximations have been proposed. The major interest of this functional is to enable joint image restoration and contour detection. In this work, we propose a general formulation of the discrete counterpart of the Mumford-Shah functional, adapted to nonsmooth penalizations, fitting the assumptions required by the Proximal Alternating Linearized Minimization (PALM), with convergence guarantees. A second contribution aims to relax some assumptions on the involved functionals and derive a novel Semi-Linearized Proximal Alternated Minimization (SL-PAM) algorithm, with proved convergence. We compare the performances of the algorithm with several nonsmooth penalizations, for Gaussian and Poisson denoising, image restoration and RGB-color denoising. We compare the results with state-of-the-art convex relaxations of the Mumford-Shah functional, and a discrete version of the Ambrosio-Tortorelli functional. We show that the SL-PAM algorithm is faster than the original PALM algorithm, and leads to competitive denoising, restoration and segmentation results.

Feasibility and robustness of dynamic 18F-FET PET based tracer kinetic models applied to patients with recurrent high-grade glioma prior to carbon ion irradiation

The aim of this study was to analyze the robustness and diagnostic value of different compartment models for dynamic 18F-FET PET in recurrent high-grade glioma (HGG). Dynamic 18F-FET PET data of patients with recurrent WHO grade III (n:7) and WHO grade IV (n: 9) tumors undergoing re-irradiation with carbon ions were analyzed by voxelwise fitting of the time-activity curves with a simplified and an extended one-tissue compartment model (1TCM) and a two-tissue compartment model (2TCM), respectively. A simulation study was conducted to assess robustness and precision of the 2TCM. Parameter maps showed enhanced detail on tumor substructure. Neglecting the blood volume VB in the 1TCM yields insufficient results. Parameter K1 from both 1TCM and 2TCM showed correlation with overall patient survival after carbon ion irradiation (p = 0.043 and 0.036, respectively). The 2TCM yields realistic estimates for tumor blood volume, which was found to be significantly higher in WHO IV compared to WHO III (p = 0.031). Simulations on the 2TCM showed that K1 yields good accuracy and robustness while k2 showed lowest stability of all parameters. The 1TCM provides the best compromise between parameter stability and model accuracy; however application of the 2TCM is still feasible and provides a more accurate representation of tracer-kinetics at the cost of reduced robustness. Detailed tracer kinetic analysis of 18F-FET PET with compartment models holds valuable information on tumor substructures and provides additional diagnostic and prognostic value.

Sparse view cone beam X-ray luminescence tomography based on truncated singular value decomposition

Cone beam X-ray luminescence computed tomography (CB-XLCT) has been proposed as a promising hybrid imaging technique. Though it has the advantage of fast imaging, the inverse problem of CB-XLCT is seriously ill-conditioned, making the image quality quite poor, especially for imaging multi-targets. To achieve fast imaging of multi-targets, which is essential for in vivo applications, a truncated singular value decomposition (TSVD) based sparse view CB-XLCT reconstruction method is proposed in this study. With the weight matrix of the CB-XLCT system being converted to orthogonal by TSVD, the compressed sensing (CS) based L₁-norm method could be applied for fast reconstruction from fewer projection views. Numerical simulations and phantom experiments demonstrate that by using the proposed method, two targets with different edge-to-edge distances (EEDs) could be resolved effectively. It indicates that the proposed method could improve the imaging quality of multi-targets significantly in terms of localization accuracy, target shape, image contrast, and spatial resolution, when compared with conventional methods.

Enhancement of dynamic myocardial perfusion PET images based on low-rank plus sparse decomposition

BACKGROUND AND OBJECTIVE

The absolute quantification of dynamic myocardial perfusion (MP) PET imaging is challenged by the limited spatial resolution of individual frame images due to division of the data into shorter frames. This study aims to develop a method for restoration and enhancement of dynamic PET images.

METHODS

We propose that the image restoration model should be based on multiple constraints rather than a single constraint, given the fact that the image characteristic is hardly described by a single constraint alone. At the same time, it may be possible, but not optimal, to regularize the image with multiple constraints simultaneously. Fortunately, MP PET images can be decomposed into a superposition of background vs. dynamic components via low-rank plus sparse (L + S) decomposition. Thus, we propose an L + S decomposition based MP PET image restoration model and express it as a convex optimization problem. An iterative soft thresholding algorithm was developed to solve the problem. Using realistic dynamic ⁸²Rb MP PET scan data, we optimized and compared its performance with other restoration methods.

RESULTS

The proposed method resulted in substantial visual as well as quantitative accuracy improvements in terms of noise versus bias performance, as demonstrated in extensive ⁸²Rb MP PET simulations. In particular, the myocardium defect in the MP PET images had improved visual as well as contrast versus noise tradeoff. The proposed algorithm was also applied on an 8-min clinical cardiac ⁸²Rb MP PET study performed on the GE Discovery PET/CT, and demonstrated improved quantitative accuracy (CNR and SNR) compared to other algorithms.

CONCLUSIONS

The proposed method is effective for restoration and enhancement of dynamic PET images.

Advances in PET/MR instrumentation and image reconstruction

The combination of positron emission tomography (PET) and MRI has attracted the attention of researchers in the past approximately 20 years in small-animal imaging and more recently in clinical research. The combination of PET/MRI allows researchers to explore clinical and research questions in a wide number of fields, some of which are briefly mentioned here. An important number of groups have developed different concepts to tackle the problems that PET instrumentation poses to the exposition of electromagnetic fields. We have described most of these research developments in preclinical and clinical experiments, including the few commercial scanners available. From the software perspective, an important number of algorithms have been developed to address the attenuation correction issue and to exploit the possibility that MRI provides for motion correction and quantitative image reconstruction, especially parametric modelling of radiopharmaceutical kinetics. In this work, we give an overview of some exemplar applications of simultaneous PET/MRI, together with technological hardware and software developments.

Deep reconstruction model for dynamic PET images

Accurate and robust tomographic reconstruction from dynamic positron emission tomography (PET) acquired data is a difficult problem. Conventional methods, such as the maximum likelihood expectation maximization (MLEM) algorithm for reconstructing the activity distribution-based on individual frames, may lead to inaccurate results due to the checkerboard effect and limitation of photon counts. In this paper, we propose a stacked sparse auto-encoder based reconstruction framework for dynamic PET imaging. The dynamic reconstruction problem is formulated in a deep learning representation, where the encoding layers extract the prototype features, such as edges, so that, in the decoding layers, the reconstructed results are obtained through a combination of those features. The qualitative and quantitative results of the procedure, including the data based on a Monte Carlo simulation and real patient data demonstrates the effectiveness of our method.

Federated optimisation of kinetic analysis problems

Positron Emission Tomography (PET) data is intrinsically dynamic, and kinetic analysis of dynamic PET data can substantially augment the information provided by static PET reconstructions. Yet despite the insights into disease that kinetic analysis offers, it is not used clinically and seldom used in research beyond the preclinical stage. The utility of PET kinetic analysis is hampered by several factors including spatial inconsistency within regions of homogeneous tissue and relative computational expense when fitting complex models to individual voxels. Even with sophisticated algorithms inconsistencies can arise because local optima frequently have narrow basins of convergence, are surrounded by relatively flat (uninformative) regions, have relatively low-gradient valley floors, or combinations thereof. Based on the observation that cost functions for individual voxels frequently bear some resemblance to each-other, this paper proposes the federated optimisation of the individual kinetic analysis problems within a given image. This approach shares parameters proposed during optimisation with other, similar voxels. Federated optimisation exploits the redundancy typical of large medical images to improve the optimisation residuals, computational efficiency and, to a limited extent, image consistency. This is achieved without restricting the formulation of the kinetic model, resorting to an explicit regularisation parameter, or limiting the resolution at which parameters are computed.

Dynamic PET reconstruction using temporal patch-based low rank penalty for ROI-based brain kinetic analysis

Dynamic positron emission tomography (PET) is widely used to measure changes in the bio-distribution of radiopharmaceuticals within particular organs of interest over time. However, to retain sufficient temporal resolution, the number of photon counts in each time frame must be limited. Therefore, conventional reconstruction algorithms such as the ordered subset expectation maximization (OSEM) produce noisy reconstruction images, thus degrading the quality of the extracted time activity curves (TACs). To address this issue, many advanced reconstruction algorithms have been developed using various spatio-temporal regularizations. In this paper, we extend earlier results and develop a novel temporal regularization, which exploits the self-similarity of patches that are collected in dynamic images. The main contribution of this paper is to demonstrate that the correlation of patches can be exploited using a low-rank constraint that is insensitive to global intensity variations. The resulting optimization framework is, however, non-Lipschitz and nonconvex due to the Poisson log-likelihood and low-rank penalty terms. Direct application of the conventional Poisson image deconvolution by an augmented Lagrangian (PIDAL) algorithm is, however, problematic due to its large memory requirements, which prevents its parallelization. Thus, we propose a novel optimization framework using the concave-convex procedure (CCCP)

Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging

In dynamic Positron Emission Tomography (PET), an estimate of the radio activity concentration is obtained from a series of frames of sinogram data taken at ranging in duration from 10 seconds to minutes under some criteria. So far, all the well-known reconstruction algorithms require known data statistical properties. It limits the speed of data acquisition, besides, it is unable to afford the separated information about the structure and the variation of shape and rate of metabolism which play a major role in improving the visualization of contrast for some requirement of the diagnosing in application. This paper presents a novel low rank-based activity map reconstruction scheme from emission sinograms of dynamic PET, termed as SLCR representing Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging. In this method, the stationary background is formulated as a low rank component while variations between successive frames are abstracted to the sparse. The resulting nuclear norm and l₁ norm related minimization problem can also be efficiently solved by many recently developed numerical methods. In this paper, the linearized alternating direction method is applied. The effectiveness of the proposed scheme is illustrated on three data sets.

Homotopic non-local regularized reconstruction from sparse positron emission tomography measurements

Background

Positron emission tomography scanners collect measurements of a patient’s in vivo radiotracer distribution. The system detects pairs of gamma rays emitted indirectly by a positron-emitting radionuclide (tracer), which is introduced into the body on a biologically active molecule, and the tomograms must be reconstructed from projections. The reconstruction of tomograms from the acquired PET data is an inverse problem that requires regularization. The use of tightly packed discrete detector rings, although improves signal-to-noise ratio, are often associated with high costs of positron emission tomography systems. Thus a sparse reconstruction, which would be capable of overcoming the noise effect while allowing for a reduced number of detectors, would have a great deal to offer.

Methods

In this study, we introduce and investigate the potential of a homotopic non-local regularization reconstruction framework for effectively reconstructing positron emission tomograms from such sparse measurements.

Results

Results obtained using the proposed approach are compared with traditional filtered back-projection as well as expectation maximization reconstruction with total variation regularization.

Conclusions

A new reconstruction method was developed for the purpose of improving the quality of positron emission tomography reconstruction from sparse measurements. We illustrate that promising reconstruction performance can be achieved for the proposed approach even at low sampling fractions, which allows for the use of significantly fewer detectors and have the potential to reduce scanner costs.

Simultaneous Reconstruction and Segmentation of Dynamic PET via Low-Rank and Sparse Matrix Decomposition

Although of great clinical value, accurate and robust reconstruction and segmentation of dynamic positron emission tomography (PET) images are great challenges due to low spatial resolution and high noise. In this paper, we propose a unified framework that exploits temporal correlations and variations within image sequences based on low-rank and sparse matrix decomposition. Thus, the two separate inverse problems, PET image reconstruction and segmentation, are accomplished in a simultaneous fashion. Considering low signal to noise ratio and piece-wise constant assumption of PET images, we also propose to regularize low-rank and sparse matrices with vectorial total variation norm. The resulting optimization problem is solved by augmented Lagrangian multiplier method with variable splitting. The effectiveness of proposed approach is validated on realistic Monte Carlo simulation datasets and the real patient data.

Direct estimation of kinetic parametric images for dynamic PET

Dynamic positron emission tomography (PET) can monitor spatiotemporal distribution of radiotracer in vivo. The spatiotemporal information can be used to estimate parametric images of radiotracer kinetics that are of physiological and biochemical interests. Direct estimation of parametric images from raw projection data allows accurate noise modeling and has been shown to offer better image quality than conventional indirect methods, which reconstruct a sequence of PET images first and then perform tracer kinetic modeling pixel-by-pixel. Direct reconstruction of parametric images has gained increasing interests with the advances in computing hardware. Many direct reconstruction algorithms have been developed for different kinetic models. In this paper we review the recent progress in the development of direct reconstruction algorithms for parametric image estimation. Algorithms for linear and nonlinear kinetic models are described and their properties are discussed.

Direct 4-D PET list mode parametric reconstruction with a novel EM algorithm

The production of images of kinetic parameters is often the ultimate goal of positron emission tomography (PET) imaging. The indirect method of PET parametric imaging, also called the frame-based method (FM), is performed by fitting the time-activity curve (TAC) for each voxel with an appropriate compartment model after image reconstruction. The indirect method is simple and easily implemented, however, it usually leads to some loss of accuracy or precision, due to the use of two separate steps. This paper presents a direct 4-D method for producing 3-D images of kinetic parameters from list mode PET data. In this application, the TAC for each voxel is described by a one-tissue compartment model (1T). Extending previous EM algorithms, a new spatiotemporal complete data space was introduced to optimize the maximum likelihood function. This leads to a straightforward closed-form parametric image update equation. This method was implemented by extending the current list mode platform MOLAR to produce a parametric algorithm PMOLAR-1T. Using an ordered subset approach, qualitative and quantitative evaluations were performed using 2-D (x, t) and 4-D (x, y, z, t) simulated list mode data based on brain receptor tracers and also with a human brain study. Comparisons with the indirect method showed that the proposed direct method can lead to accurate estimation of the parametric image values with reduced variance, especially at low count levels. In the 2-D test, the direct method showed similar bias to the frame-based method but with variance reduction of 23%-60%. In the 4-D test, bias values of both methods were no more than 4% and the direct method had lower variability (coefficient of variation reduction of 0%-64% compared to the frame-based method) at the normal count level. The direct method had a larger reduction in variability (27%-81%) and lower bias (1%-5% for 4-D and 1%-19% for FM) at low count levels. The results in the human brain study are similar with PMOLAR-1T showing lower noise than FM.

Fourier domain closed-form formulas for estimation of kinetic parameters in reversible multi-compartment models

Background

Compared with static imaging, dynamic emission computed tomographic imaging with compartment modeling can quantify in vivo physiologic processes, providing useful information about molecular disease processes. Dynamic imaging involves estimation of kinetic rate parameters. For multi-compartment models, kinetic parameter estimation can be computationally demanding and problematic with local minima.

Methods

This paper offers a new perspective to the compartment model fitting problem where Fourier linear system theory is applied to derive closed-form formulas for estimating kinetic parameters for the two-compartment model. The proposed Fourier domain estimation method provides a unique solution, and offers very different noise response as compared to traditional non-linear chi-squared minimization techniques.

Results

The unique feature of the proposed Fourier domain method is that only low frequency components are used for kinetic parameter estimation, where the DC (i.e., the zero frequency) component in the data is treated as the most important information, and high frequency components that tend to be corrupted by statistical noise are discarded. Computer simulations show that the proposed method is robust without having to specify the initial condition. The resultant solution can be fine tuned using the traditional iterative method.

Conclusions

The proposed Fourier-domain estimation method has closed-form formulas. The proposed Fourier-domain curve-fitting method does not require an initial condition, it minimizes a quadratic objective function and a closed-form solution can be obtained. The noise is easier to control, simply by discarding the high frequency components, and emphasizing the DC component.

Kinetic parameter estimation using a closed-form expression via integration by parts

Dynamic emission computed tomographic imaging with compartment modeling can quantify in vivo physiologic processes, eliciting more information regarding underlying molecular disease processes than is obtained from static imaging. However, estimation of kinetic rate parameters for multi-compartment models can be computationally demanding and problematic due to local minima. A number of techniques for kinetic parameter estimation have been studied and are in use today, generally offering a tradeoff between computation time, robustness of fit and flexibility with differing sets of assumptions. This paper presents a means to eliminate all differential operations by using the integration-by-parts method to provide closed-form formulas, so that the mathematical model is less sensitive to data sampling and noise. A family of closed-form formulas are obtained. Computer simulations show that the proposed method is robust without having to specify the initial condition.

Wavelet-based resolution recovery using an anatomical prior provides quantitative recovery for human population phantom PET [¹¹C]raclopride data

The objective of this study was to evaluate a resolution recovery (RR) method using a variety of simulated human brain [¹¹C]raclopride positron emission tomography (PET) images. Simulated datasets of 15 numerical human phantoms were processed by a wavelet-based RR method using an anatomical prior. The anatomical prior was in the form of a hybrid segmented atlas, which combined an atlas for anatomical labelling and a PET image for functional labelling of each anatomical structure. We applied RR to both 60 min static and dynamic PET images. Recovery was quantified in 84 regions, comparing the typical 'true' value for the simulation, as obtained in normal subjects, simulated and RR PET images. The radioactivity concentration in the white matter, striatum and other cortical regions was successfully recovered for the 60 min static image of all 15 human phantoms; the dependence of the solution on accurate anatomical information was demonstrated by the difficulty of the technique to retrieve the subthalamic nuclei due to mismatch between the two atlases used for data simulation and recovery. Structural and functional synergy for resolution recovery (SFS-RR) improved quantification in the caudate and putamen, the main regions of interest, from -30.1% and -26.2% to -17.6% and -15.1%, respectively, for the 60 min static image and from -51.4% and -38.3% to -27.6% and -20.3% for the binding potential (BP(ND)) image, respectively. The proposed methodology proved effective in the RR of small structures from brain [¹¹C]raclopride PET images. The improvement is consistent across the anatomical variability of a simulated population as long as accurate anatomical segmentations are provided.

Direct 4D parametric imaging for linearized models of reversibly binding PET tracers using generalized AB-EM reconstruction

Due to high noise levels in the voxel kinetics, development of reliable parametric imaging algorithms remains one of most active areas in dynamic brain PET imaging, which in the vast majority of cases involves receptor/transporter studies with reversibly binding tracers. As such, the focus of this work has been to develop a novel direct 4D parametric image reconstruction scheme for such tracers. Based on a relative equilibrium (RE) graphical analysis formulation (Zhou et al 2009b Neuroimage 44 661-70), we developed a closed-form 4D EM algorithm to directly reconstruct distribution volume (DV) parametric images within a plasma input model, as well as DV ratio (DVR) images within a reference tissue model scheme (wherein an initial reconstruction was used to estimate the reference tissue time-activity curves). A particular challenge with the direct 4D EM formulation is that the intercept parameters in graphical (linearized) analysis of reversible tracers (e.g. Logan or RE analysis) are commonly negative (unlike for irreversible tracers, e.g. using Patlak analysis). Subsequently, we focused our attention on the AB-EM algorithm, derived by Byrne (1998, Inverse Problems 14 1455-67) to allow inclusion of prior information about the lower (A) and upper (B) bounds for image values. We then generalized this algorithm to the 4D EM framework, thus allowing negative intercept parameters. Furthermore, our 4D AB-EM algorithm incorporated and emphasized the use of spatially varying lower bounds to achieve enhanced performance. As validation, the means of parameters estimated from 55 human (11)C-raclopride dynamic PET studies were used for extensive simulations using a mathematical brain phantom. Images were reconstructed using conventional indirect as well as proposed direct parametric imaging methods. Noise versus bias quantitative measurements were performed in various regions of the brain. Direct 4D EM reconstruction resulted in notable qualitative and quantitative accuracy improvements (over 35% noise reduction, with matched bias, in both plasma and reference-tissue input models). Similar improvements were also observed in the coefficient of variation of the estimated DV and DVR values even for relatively low uptake cortical regions, suggesting the enhanced ability for robust parameter estimation. The method was also tested on a 90 min (11)C-raclopride patient study performed on the high-resolution research tomograph wherein the proposed method was shown across a variety of regions to outperform the conventional method in the sense that for a given DVR value, improved noise levels were observed.

Spatio-temporal diffusion of dynamic PET images

Positron emission tomography (PET) images are corrupted by noise. This is especially true in dynamic PET imaging where short frames are required to capture the peak of activity concentration after the radiotracer injection. High noise results in a possible bias in quantification, as the compartmental models used to estimate the kinetic parameters are sensitive to noise. This paper describes a new post-reconstruction filter to increase the signal-to-noise ratio in dynamic PET imaging. It consists in a spatio-temporal robust diffusion of the 4D image based on the time activity curve (TAC) in each voxel. It reduces the noise in homogeneous areas while preserving the distinct kinetics in regions of interest corresponding to different underlying physiological processes. Neither anatomical priors nor the kinetic model are required. We propose an automatic selection of the scale parameter involved in the diffusion process based on a robust statistical analysis of the distances between TACs. The method is evaluated using Monte Carlo simulations of brain activity distributions. We demonstrate the usefulness of the method and its superior performance over two other post-reconstruction spatial and temporal filters. Our simulations suggest that the proposed method can be used to significantly increase the signal-to-noise ratio in dynamic PET imaging.

Characterization of the ringing artifacts in rotator-based reconstruction with Monte Carlo-based resolution compensation for PET

PURPOSE

Studies have shown that Monte Carlo-based reconstruction could effectively improve the image quality of positron emission tomography. The authors have previously used a Gaussian rotator-based algorithm to efficiently reduce the computational cost for system matrix (SM) calculation and to meet the large memory requirements for SM storage. However, pronounced ringing artifacts were observed in the reconstructed image. In this article, the authors investigated and characterized these artifacts.

METHODS

The authors proposed an "ideal" rotator and used it as a baseline in the artifacts evaluation. This ideal rotator produces perfectly rotated images. The Gaussian rotator method was evaluated by a full system model and a partial system model without positron range and acolinearity, which could be compensated for by the blurring of the Gaussian rotator for 18F studies. Noiseless data, Monte Carlo simulation data, as well as acquired experimental data were used to quantitatively characterize the behavior of the artifacts.

RESULTS

The study of the noiseless data indicated that the artifacts were mainly attributed to the rotator, which further blurred the simulated system responses. The simulation study suggested that the artifacts become less pronounced and not quantitatively significant in practice. This result is consistent with the experimental data study. Better contrast recovery was achieved with an over-compensated system model. Traditionally, an undercompensated system model was preferred to avoid artifacts. The authors' studies suggest that the Gaussian rotator with the full system model yields the best image quality among the evaluated methods with considerably reduced quantitative error and quantitatively insignificant artifacts in practice.

CONCLUSIONS

The authors' investigation indicated that a moderately overcompensated system model (about 2 mm FWHM in this study) yielded better contrast recovery and quantitatively insignificant artifacts in practice.

Motion compensation for fully 4D PET reconstruction using PET superset data

Fully 4D PET image reconstruction is receiving increasing research interest due to its ability to significantly reduce spatiotemporal noise in dynamic PET imaging. However, thus far in the literature, the important issue of correcting for subject head motion has not been considered. Specifically, as a direct consequence of using temporally extensive basis functions, a single instance of movement propagates to impair the reconstruction of multiple time frames, even if no further movement occurs in those frames. Existing 3D motion compensation strategies have not yet been adapted to 4D reconstruction, and as such the benefits of 4D algorithms have not yet been reaped in a clinical setting where head movement undoubtedly occurs. This work addresses this need, developing a motion compensation method suitable for fully 4D reconstruction methods which exploits an optical tracking system to measure the head motion along with PET superset data to store the motion compensated data. List-mode events are histogrammed as PET superset data according to the measured motion, and a specially devised normalization scheme for motion compensated reconstruction from the superset data is required. This work proceeds to propose the corresponding time-dependent normalization modifications which are required for a major class of fully 4D image reconstruction algorithms (those which use linear combinations of temporal basis functions). Using realistically simulated as well as real high-resolution PET data from the HRRT, we demonstrate both the detrimental impact of subject head motion in fully 4D PET reconstruction and the efficacy of our proposed modifications to 4D algorithms. Benefits are shown both for the individual PET image frames as well as for parametric images of tracer uptake and volume of distribution for (18)F-FDG obtained from Patlak analysis.

Direct 4D List Mode Parametric Reconstruction for PET with a Novel EM Algorithm

We present a direct method for producing images of kinetic parameters from list mode PET data. The time-activity curve for each voxel is described by a one-tissue compartment, 2-parameter model. Extending previous EM algorithms, a new spatiotemporal complete data space was introduced to optimize the maximum likelihood function. This leads to a straightforward parametric image update equation with moderate additional computation requirements compared to the conventional algorithm. Qualitative and quantitative evaluations were performed using 2D (x,t) and 4D (x,y,z,t) simulated list mode data for a brain receptor study. Comparisons with the two-step approach (frame-based reconstruction followed by voxel-by-voxel parameter estimation) show that the proposed method can lead to accurate estimation of the parametric image values with reduced variance, especially for the volume of distribution (V(T)).