Exact parameter identification in PET pharmacokinetic modeling using the irreversible two tissue compartment model.

OBJECTIVE: In quantitative dynamic positron emission tomography (PET), time series of images, reflecting the tissue response to the arterial tracer supply, are reconstructed. This response is described by kinetic parameters, which are commonly determined on basis of the tracer concentration in tissue and the arterial input function. In clinical routine the latter is estimated by arterial blood sampling and analysis, which is a challenging process and thus, attempted to be derived directly from reconstructed PET images. However, a mathematical analysis about the necessity of measurements of the common arterial whole blood activity concentration, and the concentration of free non-metabolized tracer in the arterial plasma, for a successful kinetic parameter identification does not exist. Here we aim to address this problem mathematically. Approach: We consider the identification problem in simultaneous pharmacokinetic modeling of multiple regions of interests of dynamic PET data using the irreversible two-tissue compartment model analytically. In addition to this consideration, the situation of noisy measurements is addressed using Tikhonov regularization. Furthermore, numerical simulations with a regularization approach are carried out to illustrate the analytical results in a synthetic application example. Main results: We provide mathematical proofs showing that, under reasonable assumptions, all metabolic tissue parameters can be uniquely identified without requiring additional blood samples to measure the arterial input function. A connection to noisy measurement data is made via a consistency result, showing that exact reconstruction of the ground-truth tissue parameters is stably maintained in the vanishing noise limit. Furthermore, our numerical experiments suggest that an approximate reconstruction of kinetic parameters according to our analytic results is also possible in practice for moderate noise levels. Significance: The analytical result, which holds in the idealized, noiseless scenario, suggests that for irreversible tracers, fully quantitative dynamic PET imaging is in principle possible without costly arterial blood sampling and metabolite analysis.

Triple-source saddle-curve cone-beam photon counting CT image reconstruction: A simulation study.

PURPOSE: The most common detector material in the PC CT system, cannot achieve the best performance at a relatively higher photon flux rate. In the reconstruction view, the most commonly used filtered back projection, is not able to provide sufficient reconstructed image quality in spectral computed tomography (CT). Developing a triple-source saddle-curve cone-beam photon counting CT image reconstruction method can improve the temporal resolution. METHODS: Triple-source saddle-curve cone-beam trajectory was rearranged into four trajectory sets for simulation and reconstruction. Projection images in different energy bins were simulated by forward projection and photon counting CT respond model simulation. After simulation, the object was reconstructed using Katsevich's theory after photon counts correction using the pseudo inverse of photon counting CT response matrix. The material decomposition can be performed based on images in different energy bins. RESULTS: Root mean square error (RMSE) and structural similarity index (SSIM) are calculated to quantify the image quality of reconstruction images. Compared with FDK images, the RMSE for the triple-source image was improved by 27%, 21%, 14%, 8%, and 6% for the reconstrued image of 20-33, 33-47, 47-58, 58-69, 69-80 keV energy bin. The SSIM was improved by 1.031%, 0.665%, 0.396%, 0.235%, 0.174% for corresponding energy bin. The decomposition image based on corrected images shows improved RMSE and SSIM, each by 33.861% and 0.345%. SSIM of corrected decomposition image of iodine reaches 99.415% of the original image. CONCLUSIONS: A new Triple-source saddle-curve cone-beam PC CT image reconstruction method was developed in this work. The exact reconstruction of the triple-source saddle-curve improved both the image quality and temporal resolution.

Quantitative study on exact reconstruction sampling condition by verifying solution uniqueness in limited-view CT.

In limited-view computed tomography reconstruction, iterative image reconstruction with sparsity-exploiting methods, such as total variation (TV) minimization, inspired by compressive sensing, potentially claims large reductions in sampling requirements. However, a quantitative notion of this claim is non-trivial because of the ill-defined reduction in sampling achieved by the sparsity-exploiting method. In this paper, exact reconstruction sampling condition for limited-view problem is studied by verifying the uniqueness of solution in TV minimization model. Uniqueness is tested by solving a convex optimization problem derived from the sufficient and necessary condition of solution uniqueness. Through this method, the sufficient sampling number of exact reconstruction is quantified for any fixed phantom and settled geometrical parameter in the limited-view problem. This paper provides a reference to quantify the sampling condition. Three phantoms are tested to study the sampling condition of limited view exact reconstruction in this paper. The experiment results show the quantified sampling number and indicate that an object would be accurately reconstructed as the scanning range becomes narrower by increasing sampling number. The increased samplings compensate for the deficiency of the projection angle. However, the lower bound of the scanning range corresponding to three different phantoms are presented, in which an exact reconstruction cannot be obtained once the projection angular is narrowed to this extent no matter how to increase sampling.