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

Improved correction for the tissue fraction effect in lung PET/CT imaging

Recently, there has been an increased interest in imaging different pulmonary disorders using PET techniques. Previous work has shown, for static PET/CT, that air content in the lung influences reconstructed image values and that it is vital to correct for this 'tissue fraction effect' (TFE). In this paper, we extend this work to include the blood component and also investigate the TFE in dynamic imaging. CT imaging and PET kinetic modelling are used to determine fractional air and blood voxel volumes in six patients with idiopathic pulmonary fibrosis. These values are used to illustrate best and worst case scenarios when interpreting images without correcting for the TFE. In addition, the fractional volumes were used to determine correction factors for the SUV and the kinetic parameters. These were then applied to the patient images. The kinetic parameters K1 and Ki along with the static parameter SUV were all found to be affected by the TFE with both air and blood providing a significant contribution to the errors. Without corrections, errors range from 34-80% in the best case and 29-96% in the worst case. In the patient data, without correcting for the TFE, regions of high density (fibrosis) appeared to have a higher uptake than lower density (normal appearing tissue), however this was reversed after air and blood correction. The proposed correction methods are vital for quantitative and relative accuracy. Without these corrections, images may be misinterpreted

Patient-Specific Method of Generating Parametric Maps of Patlak K(i) without Blood Sampling or Metabolite Correction: A Feasibility Study.

Currently, kinetic analyses using dynamic positron emission tomography (PET) experience very limited use despite their potential for improving quantitative accuracy in several clinical and research applications. For targeted volume applications, such as radiation treatment planning, treatment monitoring, and cerebral metabolic studies, the key to implementation of these methods is the determination of an arterial input function, which can include time-consuming analysis of blood samples for metabolite correction. Targeted kinetic applications would become practical for the clinic if blood sampling and metabolite correction could be avoided. To this end, we developed a novel method (Patlak-P) of generating parametric maps that is identical to Patlak K(i) (within a global scalar multiple) but does not require the determination of the arterial input function or metabolite correction. In this initial study, we show that Patlak-P (a) mimics Patlak K(i) images in terms of visual assessment and target-to-background (TB) ratios of regions of elevated uptake, (b) has higher visual contrast and (generally) better image quality than SUV, and (c) may have an important role in improving radiotherapy planning, therapy monitoring, and neurometabolism studies

Direct Reconstruction of Linear Parametric Images from Dynamic PET Using Nonlocal Deep Image Prior

Direct reconstruction methods have been developed to estimate parametric images directly from the measured PET sinograms by combining the PET imaging model and tracer kinetics in an integrated framework. Due to limited counts received, signal-to-noise-ratio (SNR) and resolution of parametric images produced by direct reconstruction frameworks are still limited. Recently supervised deep learning methods have been successfully applied to medical imaging denoising/reconstruction when large number of high-quality training labels are available. For static PET imaging, high-quality training labels can be acquired by extending the scanning time. However, this is not feasible for dynamic PET imaging, where the scanning time is already long enough. In this work, we proposed an unsupervised deep learning framework for direct parametric reconstruction from dynamic PET, which was tested on the Patlak model and the relative equilibrium Logan model. The patient's anatomical prior image, which is readily available from PET/CT or PET/MR scans, was supplied as the network input to provide a manifold constraint, and also utilized to construct a kernel layer to perform non-local feature denoising. The linear kinetic model was embedded in the network structure as a 1x1 convolution layer. The training objective function was based on the PET statistical model. Evaluations based on dynamic datasets of 18F-FDG and 11C-PiB tracers show that the proposed framework can outperform the traditional and the kernel method-based direct reconstruction methods.Comment: 10 pages, 10 figure

PET Studies of Cerebral Levodopa Metabolism: A Review of Clinical Findings and Modeling Approaches

[18F]Fluoro-3,4-dihydroxyphenyl-l-alanine (FDOPA) was one of the first successful tracers for molecular imaging by positron emission tomography (PET), and has proven immensely valuable for studies of Parkinson’s disease. Following intravenous FDOPA injection, the decarboxylated metabolite [18F] fluorodopamine is formed and trapped within terminals of the nigrostriatal dopamine neurons; reduction in the simple ratio between striatum and cerebellum is indicative of nigrostriatal degeneration. However, the kinetic analysis of dynamic FDOPA-PET recordings is formidably complex due to the entry into brain of the plasma metabolite O-methyl-FDOPA and due to the eventual washout of decarboxylated metabolites. Linear graphical analysis relative to a reference tissue input function is popular and convenient for routine clinical studies in which serial arterial blood samples are unavailable. This simplified approach has facilitated longitudinal studies in large patient cohorts. Linear graphical analysis relative to the metabolite-corrected arterial FDOPA input yields a more physiological index of FDOPA utilization, the net blood-brain clearance. Using a constrained compartmental model, FDOPA-PET recordings can be used to calculate the relative activity of the enzyme DOPA decarboxylase in living brain. We have extended this approach so as to obtain an index of steady-state trapping of [18F]fluorodopamine in synaptic vesicles. Although simple methods of image analysis are sufficient for the purposes of routine clinical studies, the more complex approaches have revealed hidden aspects of brain dopamine in personality, healthy aging, and in the pathophysiologies of Parkinson’s disease and schizophrenia

A Semi-parametric Technique for the Quantitative Analysis of Dynamic Contrast-enhanced MR Images Based on Bayesian P-splines

Dynamic Contrast-enhanced Magnetic Resonance Imaging (DCE-MRI) is an important tool for detecting subtle kinetic changes in cancerous tissue. Quantitative analysis of DCE-MRI typically involves the convolution of an arterial input function (AIF) with a nonlinear pharmacokinetic model of the contrast agent concentration. Parameters of the kinetic model are biologically meaningful, but the optimization of the non-linear model has significant computational issues. In practice, convergence of the optimization algorithm is not guaranteed and the accuracy of the model fitting may be compromised. To overcome this problems, this paper proposes a semi-parametric penalized spline smoothing approach, with which the AIF is convolved with a set of B-splines to produce a design matrix using locally adaptive smoothing parameters based on Bayesian penalized spline models (P-splines). It has been shown that kinetic parameter estimation can be obtained from the resulting deconvolved response function, which also includes the onset of contrast enhancement. Detailed validation of the method, both with simulated and in vivo data, is provided