Models for computer simulation studies of input functions for tracer kinetic modeling with positron emission tomography

Unveiling the biological side of PET-derived biomarkers: a simulation-based approach applied to PDAC assessment

PURPOSE

Radiomics has revolutionized clinical research by enabling objective measurements of imaging-derived biomarkers. However, the true potential of radiomics necessitates a comprehensive understanding of the biological basis of extracted features to serve as a clinical decision support. In this work, we propose an end-to-end framework for the in silico simulation of [18F]FLT PET imaging process in Pancreatic Ductal Adenocarcinoma, accounting for the biological characterization of tissues (including perfusion and fibrosis) on tracer delivery. We thus establish a direct association between radiomics features and the underlying biological properties of tissues.

METHODS

We considered 4 immunohistochemically stained Whole Slide Images of pancreatic tissue of one healthy control and three patients with PDAC and/or precursor lesions. From marker-specific images, tissue-depending diffusivity properties were estimated and computational domains were built to simulate the [18F]FLT spatial-temporal uptake exploiting Partial Differential Equations and Finite Elements Method. Consequently, we simulated the imaging process obtaining surrogated PET images for the considered patients, and we performed image-derived features extraction from PET images to be mapped with biological properties via correlation estimation.

RESULTS

The framework captured the phenotypic differences and generated Time Activity Curves reflecting the underlying tissue composition. Image-derived biomarkers were ranked in view of their association with biological characteristics of the tissue, unveiling their molecular correlative. Moreover, we showed that the proposed pipeline could serve as a digital phantom to optimize the image acquisition for lesion detection.

CONCLUSIONS

This innovative framework holds the potential to enhance interpretability and reliability of radiomics, fostering the adoption in personalized nuclear medicine and patient care.

EMATA: a toolbox for the automatic extraction and modeling of arterial inputs for tracer kinetic analysis in [18F]FDG brain studies

PURPOSE

PET imaging is a pivotal tool for biomarker research aimed at personalized medicine. Leveraging the quantitative nature of PET requires knowledge of plasma radiotracer concentration. Typically, the arterial input function (AIF) is obtained through arterial cannulation, an invasive and technically demanding procedure. A less invasive alternative, especially for [18F]FDG, is the image-derived input function (IDIF), which, however, often requires correction for partial volume effect (PVE), usually performed via venous blood samples. The aim of this paper is to present EMATA: Extraction and Modeling of Arterial inputs for Tracer kinetic Analysis, an open-source MATLAB toolbox. EMATA automates IDIF extraction from [18F]FDG brain PET images and additionally includes a PVE correction procedure that does not require any blood sampling.

METHODS

To assess the toolbox generalizability and present example outputs, EMATA was applied to brain [18F]FDG dynamic data of 80 subjects, extracted from two distinct datasets (40 healthy controls, 40 glioma patients). Additionally, to compare with the reference standard, quantification using both IDIF and AIF was carried out on a third open-access dataset of 18 healthy individuals.

RESULTS

EMATA consistently performs IDIF extraction across all datasets, despite differences in scanners and acquisition protocols. Remarkably high agreement is observed when comparing Patlak's Ki between IDIF and AIF (R2: 0.98 ± 0.02).

CONCLUSION

EMATA proved adaptability to different datasets characteristics and the ability to provide arterial input functions that can be used for reliable PET quantitative analysis.

Noninvasive quantification of [18F]SynVesT-1 binding using simplified reference tissue model 2

PURPOSE

[18F]SynVesT-1, a positron emission tomography (PET) radiotracer for the synaptic vesicle glycoprotein 2A (SV2A), demonstrates kinetics similar to [11C]UCB-J, with high brain uptake, fast kinetics fitting well with the one-tissue compartment (1TC) model, and excellent test-retest reproducibility. Challenges arise due to the similarity between k2 and [Formula: see text] (efflux rate of the reference region), when applying the simplified reference tissue model (SRTM) and related methods in [11C]UCB-J studies to accurately estimate [Formula: see text]. This study evaluated the suitability of these methods to estimate [18F]SynVesT-1 binding using centrum semiovale (CS) or cerebellum (CER) as reference regions.

METHOD

Seven healthy participants underwent 120-min PET scans on the HRRT scanner with [18F]SynVesT-1. Six participants underwent test and retest scans. Arterial blood sampling and metabolite analysis provided input functions for the 1TC model, serving as the gold standard for kinetic parameters values. SRTM, coupled SRTM (SRTMC) and SRTM2 estimated were applied to estimate [Formula: see text](ref: CS) and DVRCER(ref: CER) values. For SRTM2, the population average of [Formula: see text] was determined from the 1TC model applied to the reference region. Test-retest variability and minimum scan time were also calculated.

RESULTS

The 1TC k2 (1/min) values for CS and CER were 0.031 ± 0.004 and 0.021 ± 0.002, respectively. Although SRTMC [Formula: see text] was much higher than 1TC [Formula: see text], SRTMC underestimated BPND(ref: CS) and DVRCER by an average of 3% and 1% across regions, respectively, due to similar bias in k2 and [Formula: see text] estimation. SRTM underestimated BPND(ref: CS) by an average of 3%, but with the CER as reference region, SRTM estimation was unstable and DVRCER underestimation varied by region (mean 10%). Using population average [Formula: see text] values, SRTM2 BPND and DVRCER showed the best agreement with 1TC estimates.

CONCLUSION

Our findings support the use of population [Formula: see text] value in SRTM2 with [18F]SynVesT-1 for the estimation of [Formula: see text] or DVRCER, regardless of the choice of reference region.

Numerical simulation method for the assessment of the effect of molar activity on the pharmacokinetics of radioligands in small animals

BACKGROUND

It is well recognized that the molar activity of a radioligand is an important pharmacokinetic parameter, especially in positron emission tomography (PET) of small animals. Occupation of a significant number of binding sites by radioligand molecules results in low radioligand accumulation in a target region (mass effect). Nevertheless, small-animal PET studies have often been performed without consideration of the molar activity or molar dose of radioligands. A simulation study would therefore help to assess the importance of the mass effect in small-animal PET. Here, we introduce a new compartmental model-based numerical method, which runs on commonly used spreadsheet software, to simulate the effect of molar activity or molar dose on the pharmacokinetics of radioligands.

RESULTS

Assuming a two-tissue compartmental model, time-concentration curves of a radioligand were generated using four simulation methods and the well-known Runge-Kutta numerical method. The values were compared with theoretical values obtained under an ultra-high molar activity condition (pseudo-first-order binding kinetics), a steady-state condition and an equilibrium condition (second-order binding kinetics). For all conditions, the simulation method using the simplest calculation yielded values closest to the theoretical values and comparable with those obtained using the Runge-Kutta method. To satisfy a maximum occupancy less than 5%, simulations showed that a molar activity greater than 150 GBq/μmol is required for a model radioligand when 20 MBq is administered to a 250 g rat and when the concentration of binding sites in target regions is greater than 1.25 nM.

CONCLUSIONS

The simulation method used in this study is based on a very simple calculation and runs on widely used spreadsheet software. Therefore, simulation of radioligand pharmacokinetics using this method can be performed on a personal computer and help to assess the importance of the mass effect in small-animal PET. This simulation method also enables the generation of a model time-activity curve for the evaluation of kinetic analysis methods.

A Comparison of Arterial Input Function Interpolation Methods for Patlak Plot Analysis of 68Ga-PSMA-11 PET Prostate Cancer Studies

Positron emission tomography (PET) imaging enables quantitative assessment of tissue physiology. Dynamic pharmacokinetic analysis of PET images requires accurate estimation of the radiotracer plasma input function to derive meaningful parameter estimates, and small discrepancies in parameter estimation can mimic subtle physiologic tissue variation. This study evaluates the impact of input function interpolation method on the accuracy of Patlak kinetic parameter estimation through simulations modeling the pharmacokinetic properties of [68Ga]-PSMA-11. This study evaluated both trained and untrained methods. Although the mean kinetic parameter accuracy was similar across all interpolation models, the trained node weighting interpolation model estimated accurate kinetic parameters with reduced overall variability relative to standard linear interpolation. Trained node weighting interpolation reduced kinetic parameter estimation variance by a magnitude approximating the underlying physiologic differences between normal and diseased prostatic tissue. Overall, this analysis suggests that trained node weighting improves the reliability of Patlak kinetic parameter estimation for [68Ga]-PSMA-11 PET.

Signal separation of simultaneous dual-tracer PET imaging based on global spatial information and channel attention

BACKGROUND

Simultaneous dual-tracer positron emission tomography (PET) imaging efficiently provides more complete information for disease diagnosis. The signal separation has long been a challenge of dual-tracer PET imaging. To predict the single-tracer images, we proposed a separation network based on global spatial information and channel attention, and connected it to FBP-Net to form the FBPnet-Sep model.

RESULTS

Experiments using simulated dynamic PET data were conducted to: (1) compare the proposed FBPnet-Sep model to Sep-FBPnet model and currently existing Multi-task CNN, (2) verify the effectiveness of modules incorporated in FBPnet-Sep model, (3) investigate the generalization of FBPnet-Sep model to low-dose data, and (4) investigate the application of FBPnet-Sep model to multiple tracer combinations with decay corrections. Compared to the Sep-FBPnet model and Multi-task CNN, the FBPnet-Sep model reconstructed single-tracer images with higher structural similarity, peak signal-to-noise ratio and lower mean squared error, and reconstructed time-activity curves with lower bias and variation in most regions. Excluding the Inception or channel attention module resulted in degraded image qualities. The FBPnet-Sep model showed acceptable performance when applied to low-dose data. Additionally, it could deal with multiple tracer combinations. The qualities of predicted images, as well as the accuracy of derived time-activity curves and macro-parameters were slightly improved by incorporating a decay correction module.

CONCLUSIONS

The proposed FBPnet-Sep model was considered a potential method for the reconstruction and signal separation of simultaneous dual-tracer PET imaging.

Machine learning-based approach reveals essential features for simplified TSPO PET quantification in ischemic stroke patients

INTRODUCTION

Neuroinflammation evaluation after acute ischemic stroke is a promising option for selecting an appropriate post-stroke treatment strategy. To assess neuroinflammation in vivo, translocator protein PET (TSPO PET) can be used. However, the gold standard TSPO PET quantification method includes a 90 min scan and continuous arterial blood sampling, which is challenging to perform on a routine basis. In this work, we determine what information is required for a simplified quantification approach using a machine learning algorithm.

MATERIALS AND METHODS

We analyzed data from 18 patients with ischemic stroke who received 0-90 min [18F]GE-180 PET as well as T1-weigted (T1w), FLAIR, and arterial spin labeling (ASL) MRI scans. During PET scans, five manual venous blood samples at 5, 15, 30, 60, and 85 min post injection (p.i.) were drawn, and plasma activity concentration was measured. Total distribution volume (VT) was calculated using Logan plot with the full dynamic PET and an image-derived input function (IDIF) from the carotid arteries. IDIF was scaled by a calibration factor derived from all the measured plasma activity concentrations. The calculated VT values were used for training a random forest regressor. As input features for the model, we used three late PET frames (60-70, 70-80, and 80-90 min p.i.), the ASL image reflecting perfusion, the voxel coordinates, the lesion mask, and the five plasma activity concentrations. The algorithm was validated with the leave-one-out approach. To estimate the impact of the individual features on the algorithm's performance, we used Shapley Additive Explanations (SHAP). Having determined that the three late PET frames and the plasma activity concentrations were the most important features, we tested a simplified quantification approach consisting of dividing a late PET frame by a plasma activity concentration. All the combinations of frames/samples were compared by means of concordance correlation coefficient and Bland-Altman plots.

RESULTS

When using all the input features, the algorithm predicted VT values with high accuracy (87.8 ± 8.3%) for both lesion and non-lesion voxels. The SHAP values demonstrated high impact of the late PET frames (60-70, 70-80, and 80-90 min p.i.) and plasma activity concentrations on the VT prediction, while the influence of the ASL-derived perfusion, voxel coordinates, and the lesion mask was low. Among all the combinations of the late PET frames and plasma activity concentrations, the 70-80 min p.i. frame divided by the 30 min p.i. plasma sample produced the closest VT estimate in the ischemic lesion.

CONCLUSION

Reliable TSPO PET quantification is achievable by using a single late PET frame divided by a late blood sample activity concentration.

Estimation of kinetic parameters in dynamic FDG PET imaging based on shortened protocols: a virtual clinical study

This study investigated the estimation of kinetic parameters and production of related parametric Ki images in FDG PET imaging using the proposed shortened protocol (three 3-min/bed routine static images) by means of the simulated annealing (SA) algorithm. Six realistic heterogeneous tumors and various levels of [18F] FDG uptake were simulated by the XCAT phantom. An irreversible two-tissue compartment model (2TCM) using population-based input function was employed. By keeping two routine clinical scans fixed (60-min and 90-min post injection), the effect of the early scan time on optimizing the estimation of the pharmacokinetic parameters was investigated. The SA optimization algorithm was applied to estimate micro- and macro-parameters (K1, k2, k3, Ki). The minimum bias for most parameters was observed at a scan time of 20-min, which was < 10%. A highly significant correlation (> 0.9) as well as limited bias (< 10%) were observed between kinetic parameters generated from two methods [two-tissue compartment full dynamic scan (2TCM-full) and two-tissue compartment by SA algorithm (2TCM-SA)]. The analysis showed a strong correlation (> 0.8) between (2TCM-SA) Ki and SUV images. In addition, the tumor-to-background ratio (TBR) metric in the parametric (2TCM-SA) Ki images was significantly higher than SUV, although the SUV images provide better Contrast-to-noise ratio relative to parametric (2TCM-SA) Ki images. The proposed shortened protocol by the SA algorithm can estimate the kinetic parameters in FDG PET scan with high accuracy and robustness. It was also concluded that the parametric Ki images obtained from the 2TCM-SA as a complementary image of the SUV possess more quantification information than SUV images and can be used by the nuclear medicine specialist. This method has the potential to be an alternative to a full dynamic PET scan.

Kinetic Modeling of Brain [18-F]FDG Positron Emission Tomography Time Activity Curves with Input Function Recovery (IR) Method

Accurate positron emission tomography (PET) data quantification relies on high-quality input plasma curves, but venous blood sampling may yield poor-quality data, jeopardizing modeling outcomes. In this study, we aimed to recover sub-optimal input functions by using information from the tail (5th-100th min) of curves obtained through the frequent sampling protocol and an input recovery (IR) model trained with reference curves of optimal shape. Initially, we included 170 plasma input curves from eight published studies with clamp [18F]-fluorodeoxyglucose PET exams. Model validation involved 78 brain PET studies for which compartmental model (CM) analysis was feasible (reference (ref) + training sets). Recovered curves were compared with original curves using area under curve (AUC), max peak standardized uptake value (maxSUV). CM parameters (ref + training sets) and fractional uptake rate (FUR) (all sets) were computed. Original and recovered curves from the ref set had comparable AUC (d = 0.02, not significant (NS)), maxSUV (d = 0.05, NS) and comparable brain CM results (NS). Recovered curves from the training set were different from the original according to maxSUV (d = 3) and biologically plausible according to the max theoretical K1 (53//56). Brain CM results were different in the training set (p < 0.05 for all CM parameters and brain regions) but not in the ref set. FUR showed reductions similarly in the recovered curves of the training and test sets compared to the original curves (p < 0.05 for all regions for both sets). The IR method successfully recovered the plasma inputs of poor quality, rescuing cases otherwise excluded from the kinetic modeling results. The validation approach proved useful and can be applied to different tracers and metabolic conditions.

Simultaneous estimation of a model-derived input function for quantifying cerebral glucose metabolism with [18F]FDG PET

BACKGROUND

Quantification of the cerebral metabolic rate of glucose (CMRGlu) by dynamic [18F]FDG PET requires invasive arterial sampling. Alternatives to using an arterial input function (AIF) include the simultaneous estimation (SIME) approach, which models the image-derived input function (IDIF) by a series of exponentials with coefficients obtained by fitting time activity curves (TACs) from multiple volumes-of-interest. A limitation of SIME is the assumption that the input function can be modelled accurately by a series of exponentials. Alternatively, we propose a SIME approach based on the two-tissue compartment model to extract a high signal-to-noise ratio (SNR) model-derived input function (MDIF) from the whole-brain TAC. The purpose of this study is to present the MDIF approach and its implementation in the analysis of animal and human data.

METHODS

Simulations were performed to assess the accuracy of the MDIF approach. Animal experiments were conducted to compare derived MDIFs to measured AIFs (n = 5). Using dynamic [18F]FDG PET data from neurologically healthy volunteers (n = 18), the MDIF method was compared to the original SIME-IDIF. Lastly, the feasibility of extracting parametric images was investigated by implementing a variational Bayesian parameter estimation approach.

RESULTS

Simulations demonstrated that the MDIF can be accurately extracted from a whole-brain TAC. Good agreement between MDIFs and measured AIFs was found in the animal experiments. Similarly, the MDIF-to-IDIF area-under-the-curve ratio from the human data was 1.02 ± 0.08, resulting in good agreement in grey matter CMRGlu: 24.5 ± 3.6 and 23.9 ± 3.2 mL/100 g/min for MDIF and IDIF, respectively. The MDIF method proved superior in characterizing the first pass of [18F]FDG. Groupwise parametric images obtained with the MDIF showed the expected spatial patterns.

CONCLUSIONS

A model-driven SIME method was proposed to derive high SNR input functions. Its potential was demonstrated by the good agreement between MDIFs and AIFs in animal experiments. In addition, CMRGlu estimates obtained in the human study agreed to literature values. The MDIF approach requires fewer fitting parameters than the original SIME method and has the advantage that it can model the shape of any input function. In turn, the high SNR of the MDIFs has the potential to facilitate the extraction of voxelwise parameters when combined with robust parameter estimation methods such as the variational Bayesian approach.

A short 18F-FDG imaging window triple injection neuroimaging protocol for parametric mapping in PET

BACKGROUND

In parametric PET, kinetic parameters are extracted from dynamic PET images. It is not commonly used in clinical practice because of long scan times and the requirement for an arterial input function (AIF). To address these limitations, we designed an 18F-fluorodeoxyglucose (18F-FDG) triple injection dynamic PET protocol for brain imaging with a standard field of view PET scanner using a 24-min imaging window and an input function modeled using measurements from a region of interest placed over the left ventricle.

METHODS

To test the protocol in 6 healthy participants, we examined the quality of voxel-based maps of kinetic parameters in the brain generated using the two-tissue compartment model and compared estimated parameter values with previously published values. We also utilized data from a 36-min validation imaging window to compare (1) the modeled AIF against the input function measured in the validation window; and (2) the net influx rate ([Formula: see text]) computed using parameter estimates from the short imaging window against the net influx rate obtained using Patlak analysis in the validation window.

RESULTS

Compared to the AIF measured in the validation window, the input function estimated from the short imaging window achieved a mean area under the curve error of 9%. The voxel-wise Pearson's correlation between [Formula: see text] estimates from the short imaging window and the validation imaging window exceeded 0.95.

CONCLUSION

The proposed 24-min triple injection protocol enables parametric 18F-FDG neuroimaging with noninvasive estimation of the AIF from cardiac images using a standard field of view PET scanner.

An energy costly architecture of neuromodulators for human brain evolution and cognition

In comparison to other species, the human brain exhibits one of the highest energy demands relative to body metabolism. It remains unclear whether this heightened energy demand uniformly supports an enlarged brain or if specific signaling mechanisms necessitate greater energy. We hypothesized that the regional distribution of energy demands will reveal signaling strategies that have contributed to human cognitive development. We measured the energy distribution within the brain functional connectome using multimodal brain imaging and found that signaling pathways in evolutionarily expanded regions have up to 67% higher energetic costs than those in sensory-motor regions. Additionally, histology, transcriptomic data, and molecular imaging independently reveal an up-regulation of signaling at G-protein-coupled receptors in energy-demanding regions. Our findings indicate that neuromodulator activity is predominantly involved in cognitive functions, such as reading or memory processing. This study suggests that an up-regulation of neuromodulator activity, alongside increased brain size, is a crucial aspect of human brain evolution.

Clinical application of a population-based input function (PBIF) for a shortened dynamic whole-body FDG-PET/CT protocol in patients with metastatic melanoma treated by immunotherapy

BACKGROUND

The aim was to investigate the feasibility of a shortened dynamic whole-body (dWB) FDG-PET/CT protocol and Patlak imaging using a population-based input function (PBIF), instead of an image-derived input function (IDIF) across the 60-min post-injection period, and study its effect on the FDG influx rate (Ki) quantification in patients with metastatic melanoma (MM) undergoing immunotherapy.

METHODS

Thirty-seven patients were enrolled, including a PBIF modeling group (n = 17) and an independent validation cohort (n = 20) of MM from the ongoing prospective IMMUNOPET2 trial. All dWB-PET data were acquired on Vision 600 PET/CT systems. The PBIF was fitted using a Feng's 4-compartments model and scaled to the individual IDIF tail's section within the shortened acquisition time. The area under the curve (AUC) of PBIFs was compared to respective IDIFs AUC within 9 shortened time windows (TW) in terms of linear correlation (R2) and Bland-Altman tests. Ki metrics calculated with PBIF vs IDIF on 8 organs with physiological tracer uptake, 44 tumoral lesions of MM and 11 immune-induced inflammatory sites of pseudo-progression disease were also compared (Mann-Whitney test).

RESULTS

The mean ± SD relative AUC bias was calculated at 0.5 ± 3.8% (R2 = 0.961, AUCPBIF = 1.007 × AUCIDIF). In terms of optimal use in routine practice and statistical results, the 5th-7th pass (R2 = 0.999 for both Ki mean and Ki max) and 5th-8th pass (mean ± SD bias = - 4.9 ± 6.5% for Ki mean and - 4.8% ± 5.6% for Ki max) windows were selected. There was no significant difference in Ki values from PBIF5_7 vs IDIF5_7 for physiological uptakes (p > 0.05) as well as for tumor lesions (mean ± SD Ki IDIF5_7 3.07 ± 3.27 vs Ki PBIF5_7 2.86 ± 2.96 100ml/ml/min, p = 0.586) and for inflammatory sites (mean ± SD Ki IDIF5_7 1.13 ± 0.59 vs Ki PBIF5_7 1.13 ± 0.55 100ml/ml/min, p = 0.98).

CONCLUSION

Our study showed the feasibility of a shortened dWB-PET imaging protocol with a PBIF approach, allowing to reduce acquisition duration from 70 to 20 min with reasonable bias. These findings open perspectives for its clinical use in routine practice such as treatment response assessment in oncology.

Single time point quantitation of cerebral glucose metabolism by FDG-PET without arterial sampling

BACKGROUND

Until recently, quantitation of the net influx of 2-[18F]fluorodeoxyglucose (FDG) to brain (Ki) and the cerebrometabolic rate for glucose (CMRglc) required serial arterial blood sampling in conjunction with dynamic positron emission tomography (PET) recordings. Recent technical innovations enable the identification of an image-derived input function (IDIF) from vascular structures, but are frequently still encumbered by the need for interrupted sequences or prolonged recordings that are seldom available outside of a research setting. In this study, we tested simplified methods for quantitation of FDG-Ki by linear graphic analysis relative to the descending aorta IDIF in oncology patients examined using a Biograph Vision 600 PET/CT with continuous bed motion (Aarhus) or using a recently installed Biograph Vision Quadra long-axial field-of-view (FOV) scanner (Bern).

RESULTS

Correlation analysis of the coefficients of a tri-exponential decomposition of the IDIFs measured during 67 min revealed strong relationships among the total area under the curve (AUC), the terminal normalized arterial integral (theta(52-67 min)), and the terminal image-derived arterial FDG concentration (Ca(52-67 min)). These relationships enabled estimation of the missing AUC from late recordings of the IDIF, from which we then calculated FDG-Ki in brain by two-point linear graphic analysis using a population mean ordinate intercept and the single late frame. Furthermore, certain aspects of the IDIF data from Aarhus showed a marked age-dependence, which was not hitherto reported for the case of FDG pharmacokinetics.

CONCLUSIONS

The observed interrelationships between pharmacokinetic parameters in the IDIF measured during the PET recording support quantitation of FDG-Ki in brain using a single averaged frame from the interval 52-67 min post-injection, with minimal error relative to calculation from the complete dynamic sequences.

Quantitation of dynamic total-body PET imaging: recent developments and future perspectives

BACKGROUND

Positron emission tomography (PET) scanning is an important diagnostic imaging technique used in disease diagnosis, therapy planning, treatment monitoring, and medical research. The standardized uptake value (SUV) obtained at a single time frame has been widely employed in clinical practice. Well beyond this simple static measure, more detailed metabolic information can be recovered from dynamic PET scans, followed by the recovery of arterial input function and application of appropriate tracer kinetic models. Many efforts have been devoted to the development of quantitative techniques over the last couple of decades.

CHALLENGES

The advent of new-generation total-body PET scanners characterized by ultra-high sensitivity and long axial field of view, i.e., uEXPLORER (United Imaging Healthcare), PennPET Explorer (University of Pennsylvania), and Biograph Vision Quadra (Siemens Healthineers), further stimulates valuable inspiration to derive kinetics for multiple organs simultaneously. But some emerging issues also need to be addressed, e.g., the large-scale data size and organ-specific physiology. The direct implementation of classical methods for total-body PET imaging without proper validation may lead to less accurate results.

CONCLUSIONS

In this contribution, the published dynamic total-body PET datasets are outlined, and several challenges/opportunities for quantitation of such types of studies are presented. An overview of the basic equation, calculation of input function (based on blood sampling, image, population or mathematical model), and kinetic analysis encompassing parametric (compartmental model, graphical plot and spectral analysis) and non-parametric (B-spline and piece-wise basis elements) approaches is provided. The discussion mainly focuses on the feasibilities, recent developments, and future perspectives of these methodologies for a diverse-tissue environment.

Dispersion-corrected extracorporeal arterial input functions in PET studies of mice: a comparison to intracorporeal microprobe measurements

BACKGROUND

Kinetic modelling of dynamic PET typically requires knowledge of the arterial radiotracer concentration (arterial input function, AIF). Its accurate determination is very difficult in mice. AIF measurements in an extracorporeal shunt can be performed; however, this introduces catheter dispersion. We propose a framework for extracorporeal dispersion correction and validated it by comparison to invasively determined intracorporeal AIFs using implanted microprobes.

RESULTS

The response of an extracorporeal radiation detector to radioactivity boxcar functions, characterised by a convolution-based dispersion model, gave best fits using double-gamma variate and single-gamma variate kernels compared to mono-exponential kernels for the investigated range of flow rates. Parametric deconvolution with the optimal kernels was performed on 9 mice that were injected with a bolus of 39 ± 25 MBq [18F]F-PSMA-1007 after application of an extracorporeal circulation for three different flow rates in order to correct for dispersion. Comparison with synchronous implantation of microprobes for invasive aortic AIF recordings showed favourable correspondence, with no significant difference in terms of area-under-curve after 300 s and 5000 s. One-tissue and two-tissue compartment model simulations were performed to investigate differences in kinetic parameters between intra- and extracorporeally measured AIFs. Results of the modelling study revealed kinetic parameters close to the chosen simulated values in all compartment models.

CONCLUSION

The high correspondence of simultaneously intra- and extracorporeally determined AIFs and resulting model parameters establishes a feasible framework for extracorporeal dispersion correction. This should allow more precise and accurate kinetic modelling in small animal experiments.

PET AIF estimation when available ROI data is impacted by dispersive and/or background effects

Objective.Blood pool region of interest (ROI) data extracted from the field of view of a PET scanner can be impacted by both dispersive and background effects. This circumstance compromises the ability to correctly extract the arterial input function (AIF) signal. The paper explores a novel approach to addressing this difficulty.Approach.The method involves representing the AIF in terms of the whole-body impulse response (IR) to the injection profile. Analysis of a collection/population of directly sampled arterial data sets allows the statistical behaviour of the tracer's impulse response to be evaluated. It is proposed that this information be used to develop a penalty term for construction of a data-adaptive method of regularisation estimator of the AIF when dispersive and/or background effects maybe impacting the blood pool ROI data.Main results.Computational efficiency of the approach derives from the linearity of the impulse response representation of the AIF and the ability to substantially rely on quadratic programming techniques for numerical implementation. Data from eight different tracers, used in PET cancer imaging studies, are considered. Sample image-based AIF extractions for brain studies with:18F-labeled fluoro-deoxyglucose and fluoro-thymidine (FLT),11C-labeled carbon dioxide (CO2) and15O-labeled water (H2O) are presented. Results are compared to the true AIF based on direct arterial sampling. Formal numerical simulations are used to evaluate the performance of the AIF extraction method when the ROI data has varying amounts of contamination, in comparison to a direct approach that ignores such effects. It is found that even with quite small amounts of contamination, the mean squared error of the regularised AIF is significantly better than the error associated with direct use of the ROI data.Significance.The proposed IR-based AIF extraction scheme offers a practical methodological approach for situations where the available image ROI data may be contaminated by background and/or dispersion effects.

Comparison of cardiac image-derived input functions for quantitative whole body [18F]FDG imaging with arterial blood sampling

Introduction: Dynamic positron emission tomography (PET) and the application of kinetic models can provide important quantitative information based on its temporal information. This however requires arterial blood sampling, which can be challenging to acquire. Nowadays, state-of-the-art PET/CT systems offer fully automated, whole-body (WB) kinetic modelling protocols using image-derived input functions (IDIF) to replace arterial blood sampling. Here, we compared the validity of an automatic WB kinetic model protocol to the reference standard arterial input function (AIF) for both clinical and research settings. Methods: Sixteen healthy participants underwent dynamic WB [18F]FDG scans using a continuous bed motion PET/CT system with simultaneous arterial blood sampling. Multiple processing pipelines that included automatic and manually generated IDIFs derived from the aorta and left ventricle, with and without motion correction were compared to the AIF. Subsequently generated quantitative images of glucose metabolism were compared to evaluate performance of the different input functions. Results: We observed moderate to high correlations between IDIFs and the AIF regarding area under the curve (r = 0.49-0.89) as well as for the cerebral metabolic rate of glucose (CMRGlu) (r = 0.68-0.95). Manual placing of IDIFs and motion correction further improved their similarity to the AIF. Discussion: In general, the automatic vendor protocol is a feasible approach for the quantification of CMRGlu for both, clinical and research settings where expertise or time is not available. However, we advise on a rigorous inspection of the placement of the volume of interest, the resulting IDIF, and the quantitative values to ensure valid interpretations. In protocols requiring longer scan times or where cohorts are prone to involuntary movement, manual IDIF definition with additional motion correction is recommended, as this has greater accuracy and reliability.

Direct reconstruction for simultaneous dual-tracer PET imaging based on multi-task learning

BACKGROUND

Simultaneous dual-tracer positron emission tomography (PET) imaging can observe two molecular targets in a single scan, which is conducive to disease diagnosis and tracking. Since the signals emitted by different tracers are the same, it is crucial to separate each single tracer from the mixed signals. The current study proposed a novel deep learning-based method to reconstruct single-tracer activity distributions from the dual-tracer sinogram.

METHODS

We proposed the Multi-task CNN, a three-dimensional convolutional neural network (CNN) based on a framework of multi-task learning. One common encoder extracted features from the dual-tracer dynamic sinogram, followed by two distinct and parallel decoders which reconstructed the single-tracer dynamic images of two tracers separately. The model was evaluated by mean squared error (MSE), multiscale structural similarity (MS-SSIM) index and peak signal-to-noise ratio (PSNR) on simulated data and real animal data, and compared to the filtered back-projection method based on deep learning (FBP-CNN).

RESULTS

In the simulation experiments, the Multi-task CNN reconstructed single-tracer images with lower MSE, higher MS-SSIM and PSNR than FBP-CNN, and was more robust to the changes in individual difference, tracer combination and scanning protocol. In the experiment of rats with an orthotopic xenograft glioma model, the Multi-task CNN reconstructions also showed higher qualities than FBP-CNN reconstructions.

CONCLUSIONS

The proposed Multi-task CNN could effectively reconstruct the dynamic activity images of two single tracers from the dual-tracer dynamic sinogram, which was potential in the direct reconstruction for real simultaneous dual-tracer PET imaging data in future.

Patient-Adaptive Population-Based Modeling of Arterial Input Functions

Kinetic modeling of dynamic PET data requires knowledge of tracer concentration in blood plasma, described by the arterial input function (AIF). Arterial blood sampling is the gold standard for AIF measurement, but is invasive and labour intensive. A number of methods have been proposed to accurately estimate the AIF directly from blood sampling and/or imaging data. Here we consider fitting a patient-adaptive mixture of historical population time course profiles to estimate individual AIFs. Travel time of a tracer atom from the injection site to the right ventricle of the heart is modeled as a realization from a Gamma distribution, and the time this atom spends in circulation before being sampled is represented by a subject-specific linear mixture of population profiles. These functions are estimated from independent population data. Individual AIFs are obtained by projection onto this basis of population profile components. The model incorporates knowledge of injection duration into the fit, allowing for varying injection protocols. Analyses of arterial sampling data from 18F-FDG, 15O-H2O and 18F-FLT clinical studies show that the proposed model can outperform reference techniques. The statistically significant gain achieved by using population data to train the basis components, instead of fitting these from the single individual sampling data, is measured on the FDG cohort. Kinetic analyses of simulated data demonstrate the reliability and potential benefit of this approach in estimating physiological parameters. These results are further supported by numerical simulations that demonstrate convergence and stability of the proposed technique under varying training population sizes and noise levels.

Clinical validation of a population-based input function for 20-min dynamic whole-body 18F-FDG multiparametric PET imaging

Purpose

Contemporary PET/CT scanners can use 70-min dynamic whole-body (D-WB) PET to generate more quantitative information about FDG uptake than just the SUV by generating parametric images of FDG metabolic rate (MR_FDG). The analysis requires the late (50–70 min) D-WB tissue data combined with the full (0–70 min) arterial input function (AIF). Our aim was to assess whether the use of a scaled population-based input function (sPBIF) obviates the need for the early D-WB PET acquisition and allows for a clinically feasible 20-min D-WB PET examination.

Methods

A PBIF was calculated based on AIFs from 20 patients that were D-WB PET scanned for 120 min with simultaneous arterial blood sampling. MR_FDG imaging using PBIF requires that the area under the curve (AUC) of the sPBIF is equal to the AUC of the individual patient’s input function because sPBIF AUC bias translates into MR_FDG bias. Special patient characteristics could affect the shape of their AIF. Thus, we validated the use of PBIF in 171 patients that were divided into 12 subgroups according to the following characteristics: diabetes, cardiac ejection fraction, blood pressure, weight, eGFR and age. For each patient, the PBIF was scaled to the aorta image-derived input function (IDIF) to calculate a sPBIF, and the AUC bias was calculated.

Results

We found excellent agreement between the AIF and IDIF at all times. For the clinical validation, the use of sPBIF led to an acceptable AUC bias of 1–5% in most subgroups except for patients with diabetes or patients with low eGFR, where the biases were marginally higher at 7%. Multiparametric MR_FDG images based on a short 20-min D-WB PET and sPBIF were visually indistinguishable from images produced by the full 70-min D-WB PET and individual IDIF.

Conclusions

A short 20-min D-WB PET examination using PBIF can be used for multiparametric imaging without compromising the image quality or precision of MR_FDG. The D-WB PET examination may therefore be used in clinical routine for a wide range of patients, potentially allowing for more precise quantification in e.g. treatment response imaging.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-022-00490-y.

Efficient Delay Correction for Total-Body PET Kinetic Modeling Using Pulse Timing Methods

Quantitative kinetic modeling requires an input function. A noninvasive image-derived input function (IDIF) can be obtained from dynamic PET images. However, a robust IDIF location (e.g., aorta) may be far from a tissue of interest, particularly in total-body PET, introducing a time delay between the IDIF and the tissue. The standard practice of joint estimation (JE) of delay, along with model fitting, is computationally expensive. To improve the efficiency of delay correction for total-body PET parametric imaging, this study investigated the use of pulse timing methods to estimate and correct for delay. Methods: Simulation studies were performed with a range of delay values, frame lengths, and noise levels to test the tolerance of 2 pulse timing methods-leading edge (LE) and constant fraction discrimination and their thresholds. The methods were then applied to data from 21 subjects (14 healthy volunteers, 7 cancer patients) who underwent a 60-min dynamic total-body ¹⁸F-FDG PET acquisition. Region-of-interest kinetic analysis was performed and parametric images were generated to compare LE and JE methods of delay correction, as well as no delay correction. Results: Simulations demonstrated that a 10% LE threshold resulted in biases and SDs at tolerable levels for all noise levels tested, with 2-s frames. Pooled region-of-interest-based results (n = 154) showed strong agreement between LE (10% threshold) and JE methods in estimating delay (Pearson r = 0.96, P < 0.001) and the kinetic parameters v_b (r = 0.96, P < 0.001), K_i (r = 1.00, P < 0.001), and K₁ (r = 0.97, P < 0.001). When tissues with minimal delay were excluded from pooled analyses, there were reductions in v_b (69.4%) and K₁ (4.8%) when delay correction was not performed. Similar results were obtained for parametric images; additionally, lesion K_i contrast was improved overall with LE and JE delay correction compared with no delay correction and Patlak analysis. Conclusion: This study demonstrated the importance of delay correction in total-body PET. LE delay correction can be an efficient surrogate for JE, requiring a fraction of the computational time and allowing for rapid delay correction across more than 10⁶ voxels in total-body PET datasets.

Radiogenomic analysis of primary breast cancer reveals [18F]-fluorodeoxglucose dynamic flux-constants are positively associated with immune pathways and outperform static uptake measures in associating with glucose metabolism

BACKGROUND

PET imaging of 18F-fluorodeoxygucose (FDG) is used widely for tumour staging and assessment of treatment response, but the biology associated with FDG uptake is still not fully elucidated. We therefore carried out gene set enrichment analyses (GSEA) of RNA sequencing data to find KEGG pathways associated with FDG uptake in primary breast cancers.

METHODS

Pre-treatment data were analysed from a window-of-opportunity study in which 30 patients underwent static and dynamic FDG-PET and tumour biopsy. Kinetic models were fitted to dynamic images, and GSEA was performed for enrichment scores reflecting Pearson and Spearman coefficients of correlations between gene expression and imaging.

RESULTS

A total of 38 pathways were associated with kinetic model flux-constants or static measures of FDG uptake, all positively. The associated pathways included glycolysis/gluconeogenesis ('GLYC-GLUC') which mediates FDG uptake and was associated with model flux-constants but not with static uptake measures, and 28 pathways related to immune-response or inflammation. More pathways, 32, were associated with the flux-constant K of the simple Patlak model than with any other imaging index. Numbers of pathways categorised as being associated with individual micro-parameters of the kinetic models were substantially fewer than numbers associated with flux-constants, and lay around levels expected by chance.

CONCLUSIONS

In pre-treatment images GLYC-GLUC was associated with FDG kinetic flux-constants including Patlak K, but not with static uptake measures. Immune-related pathways were associated with flux-constants and static uptake. Patlak K was associated with more pathways than were the flux-constants of more complex kinetic models. On the basis of these results Patlak analysis of dynamic FDG-PET scans is advantageous, compared to other kinetic analyses or static imaging, in studies seeking to infer tumour-to-tumour differences in biology from differences in imaging. Trial registration NCT01266486, December 24th 2010.

Temporal information guided dynamic dual-tracer PET signal separation network

PURPOSE

The difficulty of dynamic dual-tracer positron emission tomography (PET) technology is to separate the complete single-tracer information from mixed dual-tracer. Traditional methods cannot separate single injection single-scan dynamic dual-tracer PET images. In this paper, we propose a deep learning framework based on gated recurrent unit (GRU) network and evaluate its performance with simulation experiments and realistic monkey data.

METHODS

The proposed single-scan dynamic dual-tracer PET image separation network consists of three parts, including encoder, separation and decoder module. Encoder part is to map the mixed time activity curves (TACs) from the low-dimensional space to the high-dimensional space to get mixed weight vector matrix. Separation part is to capture the temporal information of mixed weight vector matrix using bi-directional GRU (bi-GRU) layer to obtain the single-tracer masks, and the decoding part remaps the high-dimensional single-tracer weight vector matrix to the low-dimensional space to obtain two separated single tracers.

RESULTS

In the simulation experiments under different tracers, phantoms, noise levels, arterial input function (AIF) and k-parameter with Gaussian random, compared to the stacked auto encoder (SAE) network and traditional background subtraction method, GRU-based network has better performance with low bias and mean squared error (MSE). In the realistic study, the image results of GRU network have higher mean structural similarity (MSSIM), and peak signal to noise ratio (PSNR).

CONCLUSIONS

This study demonstrates the feasibility of temporal information guided neural network in single-injection single-scan dynamic dual-tracer PET images separation. The GRU-based network uses TAC temporal information without AIFs to make the separation results more robust and accurate, which significantly outperforms state-of-the-art method qualitatively and quantitatively. This article is protected by copyright. All rights reserved.

Reduced Acquisition Time [18F]GE-180 PET Scanning Protocol Replaces Gold-Standard Dynamic Acquisition in a Mouse Ischemic Stroke Model

Aim

Understanding neuroinflammation after acute ischemic stroke is a crucial step on the way to an individualized post-stroke treatment. Microglia activation, an essential part of neuroinflammation, can be assessed using [18F]GE-180 18 kDa translocator protein positron emission tomography (TSPO-PET). However, the commonly used 60–90 min post-injection (p.i.) time window was not yet proven to be suitable for post-stroke neuroinflammation assessment. In this study, we compare semi-quantitative estimates derived from late time frames to quantitative estimates calculated using a full 0–90 min dynamic scan in a mouse photothrombotic stroke (PT) model.

Materials and Methods

Six mice after PT and six sham mice were included in the study. For a half of the mice, we acquired four serial 0–90 min scans per mouse (analysis cohort) and calculated standardized uptake value ratios (SUVRs; cerebellar reference) for the PT volume of interest (VOI) in five late 10 min time frames as well as distribution volume ratios (DVRs) for the same VOI. We compared late static 10 min SUVRs and the 60–90 min time frame of the analysis cohort to the corresponding DVRs by linear fitting. The other half of the animals received a static 60–90 min scan and was used as a validation cohort. We extrapolated DVRs by using the static 60–90 min p.i. time window, which were compared to the DVRs of the analysis cohort.

Results

We found high linear correlations between SUVRs and DVRs in the analysis cohort for all studied 10 min time frames, while the fits of the 60–70, 70–80, and 80–90 min p.i. time frames were the ones closest to the line of identity. For the 60–90 min time window, we observed an excellent linear correlation between SUVR and DVR regardless of the phenotype (PT vs. sham). The extrapolated DVRs of the validation cohort were not significantly different from the DVRs of the analysis group.

Conclusion

Simplified quantification by a reference tissue ratio of the late 60–90 min p.i. [18F]GE-180 PET image can replace full quantification of a dynamic scan for assessment of microglial activation in the mouse PT model.

[11C]PBR28 radiotracer kinetics are not driven by alterations in cerebral blood flow

The positron emission tomography (PET) radiotracer [¹¹C]PBR28 has been increasingly used to image the translocator protein (TSPO) as a marker of neuroinflammation in a variety of brain disorders. Interrelatedly, similar clinical populations can also exhibit altered brain perfusion, as has been shown using arterial spin labelling in magnetic resonance imaging (MRI) studies. Hence, an unsolved debate has revolved around whether changes in perfusion could alter delivery, uptake, or washout of the radiotracer [¹¹C]PBR28, and thereby influence outcome measures that affect interpretation of TSPO upregulation. In this simultaneous PET/MRI study, we demonstrate that [¹¹C]PBR28 signal elevations in chronic low back pain patients are not accompanied, in the same regions, by increases in cerebral blood flow (CBF) compared to healthy controls, and that areas of marginal hypoperfusion are not accompanied by decreases in [¹¹C]PBR28 signal. In non-human primates, we show that hypercapnia-induced increases in CBF during radiotracer delivery or washout do not alter [¹¹C]PBR28 outcome measures. The combined results from two methodologically distinct experiments provide support from human data and direct experimental evidence from non-human primates that changes in CBF do not influence outcome measures reported by [¹¹C]PBR28 PET imaging studies and corresponding interpretations of the biological meaning of TSPO upregulation.

Dynamic FDG-PET in localization of focal epilepsy: A pilot study

Epilepsy surgery remains underutilized, in part because non-invasive methods of potential seizure foci localization are inadequate. We used high-resolution, parametric quantification from dynamic 2-[18F] fluoro-2-deoxy-d-glucose positron emission tomography (dFDG-PET) imaging to locate hypometabolic foci in patients whose standard clinical static PET images were normal. We obtained dFDG-PET brain images with simultaneous EEG in a one-hour acquisition on seven patients with no MRI evidence of focal epilepsy to record uptake and focal radiation decay. Images were attenuation- and motion-corrected and co-registered with high-resolution T1-weighted patient MRI and segmented into 18 regions of interest (ROI) per hemisphere. Tracer uptake was calibrated with a model corrected blood input function with partial volume (PV) corrections to generate tracer parametric maps compared between mean radiation values between hemispheres with z-scores. We identified ROI with the lowest negative z scores (<-1.65 SD) as hypometabolic. Dynamic 2-[18F] fluoro-2-deoxy-d-glucose positron emission tomography ( found focal regions of altered metabolism in all cases in which standard clinical FDG-PET found no abnormalities. This pilot study of dynamic FDG-PET suggests that further research is merited to evaluate whether glucose dynamics offer improved clinical utility for localization of epileptic foci over standard static techniques.

Association between cannabinoid 1 receptor availability and glutamate levels in healthy controls and drug-free patients with first episode psychosis: a multi-modal PET and 1H-MRS study

Cannabinoid 1 receptor and glutamatergic dysfunction have both been implicated in the pathophysiology of schizophrenia. However, it remains unclear if cannabinoid 1 receptor alterations shown in drug-naïve/free patients with first episode psychosis may be linked to glutamatergic alterations in the illness. We aimed to investigate glutamate levels and cannabinoid 1 receptor levels in the same region in patients with first episode psychosis. Forty volunteers (20 healthy volunteers, 20 drug-naïve/free patients with first episode psychosis diagnosed with schizophrenia/schizoaffective disorder) were included in the study. Glutamate levels were measured using proton magnetic resonance spectroscopy. CB1R availability was indexed using the distribution volume (VT (ml/cm3)) of [11C]MePPEP using arterial blood sampling. There were no significant associations between ACC CB1R levels and ACC glutamate levels in controls (R = - 0.24, p = 0.32) or patients (R = - 0.10, p = 0.25). However, ACC glutamate levels were negatively associated with CB1R availability in the striatum (R = - 0.50, p = 0.02) and hippocampus (R = - 0.50, p = 0.042) in controls, but these associations were not observed in patients (p > 0.05). Our findings extend our previous work in an overlapping sample to show, for the first time as far as we're aware, that cannabinoid 1 receptor alterations in the anterior cingulate cortex are shown in the absence of glutamatergic dysfunction in the same region, and indicate potential interactions between glutamatergic signalling in the anterior cingulate cortex and the endocannabinoid system in the striatum and hippocampus.

Validation of a combined image derived input function and venous sampling approach for the quantification of [18F]GE-179 PET binding in the brain

Blood-based kinetic analysis of PET data relies on an accurate estimate of the arterial plasma input function (PIF). An alternative to invasive measurements from arterial sampling is an image-derived input function (IDIF). However, an IDIF provides the whole blood radioactivity concentration, rather than the required free tracer radioactivity concentration in plasma. To estimate the tracer PIF, we corrected an IDIF from the carotid artery with estimates of plasma parent fraction (PF) and plasma-to-whole blood (PWB) ratio obtained from five venous samples. We compared the combined IDIF+venous approach to gold standard data from arterial sampling in 10 healthy volunteers undergoing [¹⁸F]GE-179 brain PET imaging of the NMDA receptor. Arterial and venous PF and PWB ratio estimates determined from 7 patients with traumatic brain injury (TBI) were also compared to assess the potential effect of medication. There was high agreement between areas under the curves of the estimates of PF (r = 0.99, p<0.001), PWB ratio (r = 0.93, p<0.001), and the PIF (r = 0.92, p<0.001) as well as total distribution volume (V_T) in 11 regions across the brain (r = 0.95, p<0.001). IDIF+venous V_T had a mean bias of -1.7% and a comparable regional coefficient of variation (arterial: 21.3 ± 2.5%, IDIF+venous: 21.5 ± 2.0%). Simplification of the IDIF+venous method to use only one venous sample provided less accurate V_T estimates (mean bias 9.9%; r = 0.71, p<0.001). A version of the method that avoids the need for blood sampling by combining the IDIF with population-based PF and PWB ratio estimates systematically underestimated V_T (mean bias -20.9%), and produced V_T estimates with a poor correlation to those obtained using arterial data (r = 0.45, p<0.001). Arterial and venous blood data from 7 TBI patients showed high correlations for PF (r = 0.92, p = 0.003) and PWB ratio (r = 0.93, p = 0.003). In conclusion, the IDIF+venous method with five venous samples provides a viable alternative to arterial sampling for quantification of [¹⁸F]GE-179 V_T.

Model Corrected Blood Input Function to Compute Cerebral FDG Uptake Rates From Dynamic Total-Body PET Images of Rats in vivo

Recently, we developed a three-compartment dual-output model that incorporates spillover (SP) and partial volume (PV) corrections to simultaneously estimate the kinetic parameters and model-corrected blood input function (MCIF) from dynamic 2-[18F] fluoro-2-deoxy-D-glucose positron emission tomography (FDG PET) images of mouse heart in vivo. In this study, we further optimized this model and utilized the estimated MCIF to compute cerebral FDG uptake rates, K_i, from dynamic total-body FDG PET images of control Wistar–Kyoto (WKY) rats and compared to those derived from arterial blood sampling in vivo. Dynamic FDG PET scans of WKY rats (n = 5), fasted for 6 h, were performed using the Albira Si Trimodal PET/SPECT/CT imager for 60 min. Arterial blood samples were collected for the entire imaging duration and then fitted to a seven-parameter function. The 60-min list mode PET data, corrected for attenuation, scatter, randoms, and decay, were reconstructed into 23 time bins. A 15-parameter dual-output model with SP and PV corrections was optimized with two cost functions to compute MCIF. A four-parameter compartment model was then used to compute cerebral Ki. The computed area under the curve (AUC) and K_i were compared to that derived from arterial blood samples. Experimental and computed AUCs were 1,893.53 ± 195.39 kBq min/cc and 1,792.65 ± 155.84 kBq min/cc, respectively (p = 0.76). Bland–Altman analysis of experimental vs. computed K_i for 35 cerebral regions in WKY rats revealed a mean difference of 0.0029 min⁻¹ (~13.5%). Direct (AUC) and indirect (Ki) comparisons of model computations with arterial blood sampling were performed in WKY rats. AUC and the downstream cerebral FDG uptake rates compared well with that obtained using arterial blood samples. Experimental vs. computed cerebral K_i for the four super regions including cerebellum, frontal cortex, hippocampus, and striatum indicated no significant differences.

Simplified protocol for whole-body Patlak parametric imaging with 18 F-FDG PET/CT: Feasibility and error analysis

PURPOSE

Parametric imaging using the Patlak model has been shown to provide improved lesion detectability and specificity. The Patlak model requires both tissue time-activity curves (TACs) after equilibrium and knowledge of the input function from the start of injection. Therefore, the conventional dynamic scanning protocol typically starts from the radiotracer injection all the way to equilibrium. In this paper, we propose the use of hybrid population-based and model-based input function estimation and evaluate its use for whole-body Patlak analysis, in order to reduce the total scan time and simplify clinical Patlak parametric imaging protocols. Possible quantitative errors caused by the simplified scanning protocol were also analyzed both theoretically and with the use of clinical data.

MATERIALS AND METHODS

Clinical data from 24 patients referred for tumor staging were included in this study. The patients underwent a whole-body dynamic PET study, 20 min after FDG injection (0.13 mCi/kg). The proposed whole-body scanning protocol includes 6 passes with 4-5 bed positions, depending on the size of the patient, with 2 min for each bed position. An input function from the literature was selected as the shape of the population-based input function. The descending aorta from the corresponding CT image was segmented and applied on the reconstructed dynamic PET images to acquire an image-based input function, which was later fitted using an exponential model. Due to the late scan time, only the later portion of the input function was available, which was used to scale the population-based input function. The hybrid input function was used to derive the whole-body Patlak images. Assuming a given error in the population-based input function, its influence on the final Patlak images were also derived theoretically and verified using the clinical data sets. Finally, the image quality of the reconstructed Patlak slope image was evaluated by an experienced radiologist in four different aspects: image artifacts, image noise, lesion sharpness, and lesion detectability.

RESULTS

It was found that errors in the population-based input function only affect the absolute scale of the Patlak slope image. The induced error is proportional to the percentage area-under-curve (AUC) error in the input function. These findings were also confirmed by numerical analysis. The predicted global scale was in good agreement with results from both image-based Patlak and direct Patlak approach. The fractions of the AUC from the early portion population-based input function were also found to be around 18% of the total AUC of the input function, further limiting the propagation of quantitation error from population-based input function to the final Patlak slope image. The reconstructed Patlak images were also found by the radiologist to provide excellent confidence in lesion detection tasks.

CONCLUSIONS

We have proposed a simplified whole-body scanning protocol that utilizes both population-based input function and model-based input function. The error from the population-based function was found to only affect the global scale and the overall quantitative impact can be predicted using our proposed formulas.

Assessment of population-based input functions for Patlak imaging of whole body dynamic 18F-FDG PET

BACKGROUND

Arterial blood sampling is the gold standard method to obtain the arterial input function (AIF) for quantification of whole body (WB) dynamic 18F-FDG PET imaging. However, this procedure is invasive and not typically available in clinical environments. As an alternative, we compared AIFs to population-based input functions (PBIFs) using two normalization methods: area under the curve (AUC) and extrapolated initial plasma concentration (CP*(0)). To scale the PBIFs, we tested two methods: (1) the AUC of the image-derived input function (IDIF) and (2) the estimated CP*(0). The aim of this study was to validate IDIF and PBIF for FDG oncological WB PET studies by comparing to the gold standard arterial blood sampling.

METHODS

The Feng 18F-FDG plasma concentration model was applied to estimate AIF parameters (n = 23). AIF normalization used either AUC(0-60 min) or CP*(0), estimated from an exponential fit. CP*(0) is also described as the ratio of the injected dose (ID) to initial distribution volume (iDV). iDV was modeled using the subject height and weight, with coefficients that were estimated in 23 subjects. In 12 oncological patients, we computed IDIF (from the aorta) and PBIFs with scaling by the AUC of the IDIF from 4 time windows (15-45, 30-60, 45-75, 60-90 min) (PBIFAUC) and estimated CP*(0) (PBIFiDV). The IDIF and PBIFs were compared with the gold standard AIF, using AUC values and Patlak Ki values.

RESULTS

The IDIF underestimated the AIF at early times and overestimated it at later times. Thus, based on the AUC and Ki comparison, 30-60 min was the most accurate time window for PBIFAUC; later time windows for scaling underestimated Ki (- 6 ± 8 to - 13 ± 9%). Correlations of AUC between AIF and IDIF, PBIFAUC(30-60), and PBIFiDV were 0.91, 0.94, and 0.90, respectively. The bias of Ki was - 9 ± 10%, - 1 ± 8%, and 3 ± 9%, respectively.

CONCLUSIONS

Both PBIF scaling methods provided good mean performance with moderate variation. Improved performance can be obtained by refining IDIF methods and by evaluating PBIFs with test-retest data.

Plasma radio-metabolite analysis of PET tracers for dynamic PET imaging: TLC and autoradiography

BACKGROUND

In molecular imaging with dynamic PET, the binding and dissociation of a targeted tracer is characterized by kinetics modeling which requires the arterial concentration of the tracer to be measured accurately. Once in the body the radiolabeled parent tracer may be subjected to hydrolysis, demethylation/dealkylation and other biochemical processes, resulting in the production and accumulation of different metabolites in blood which can be labeled with the same PET radionuclide as the parent. Since these radio-metabolites cannot be distinguished by PET scanning from the parent tracer, their contribution to the arterial concentration curve has to be removed for the accurate estimation of kinetic parameters from kinetic analysis of dynamic PET. High-performance liquid chromatography has been used to separate and measure radio-metabolites in blood plasma; however, the method is labor intensive and remains a challenge to implement for each individual patient. The purpose of this study is to develop an alternate technique based on thin layer chromatography (TLC) and a sensitive commercial autoradiography system (Beaver, Ai4R, Nantes, France) to measure radio-metabolites in blood plasma of two targeted tracers-[¹⁸F]FAZA and [¹⁸F]FEPPA, for imaging hypoxia and inflammation, respectively.

RESULTS

Radioactivity as low as 17 Bq in 2 µL of pig's plasma can be detected on the TLC plate using autoradiography. Peaks corresponding to the parent tracer and radio-metabolites could be distinguished in the line profile through each sample (n = 8) in the autoradiographic image. Significant intersubject and intra-subject variability in radio-metabolites production could be observed with both tracers. For [¹⁸F]FEPPA, 50% of plasma activity was from radio-metabolites as early as 5-min post injection, while for [¹⁸F]FAZA, significant metabolites did not appear until 50-min post. Simulation study investigating the effect of radio-metabolite in the estimation of kinetic parameters indicated that 32-400% parameter error can result without radio-metabolites correction.

CONCLUSION

TLC coupled with autoradiography is a good alternative to high-performance liquid chromatography for radio-metabolite correction. The advantages of requiring only small blood samples (~ 100 μL) and of analyzing multiple samples simultaneously, make the method suitable for individual dynamic PET studies.

FBP-Net for direct reconstruction of dynamic PET images

Dynamic positron emission tomography (PET) imaging can provide information about metabolic changes over time, used for kinetic analysis and auxiliary diagnosis. Existing deep learning-based reconstruction methods have too many trainable parameters and poor generalization, and require mass data to train the neural network. However, obtaining large amounts of medical data is expensive and time-consuming. To reduce the need for data and improve the generalization of network, we combined the filtered back-projection (FBP) algorithm with neural network, and proposed FBP-Net which could directly reconstruct PET images from sinograms instead of post-processing the rough reconstruction images obtained by traditional methods. The FBP-Net contained two parts: the FBP part and the denoiser part. The FBP part adaptively learned the frequency filter to realize the transformation from the detector domain to the image domain, and normalized the coarse reconstruction images obtained. The denoiser part merged the information of all time frames to improve the quality of dynamic PET reconstruction images, especially the early time frames. The proposed FBP-Net was performed on simulation and real dataset, and the results were compared with the state-of-art U-net and DeepPET. The results showed that FBP-Net did not tend to overfit the training set and had a stronger generalization.

Hybrid PET- and MR-driven attenuation correction for enhanced 18F-NaF and 18F-FDG quantification in cardiovascular PET/MR imaging

BACKGROUND

The standard MR Dixon-based attenuation correction (AC) method in positron emission tomography/magnetic resonance (PET/MR) imaging segments only the air, lung, fat and soft-tissues (4-class), thus neglecting the highly attenuating bone tissues and affecting quantification in bones and adjacent vessels. We sought to address this limitation by utilizing the distinctively high bone uptake rate constant Ki expected from 18F-Sodium Fluoride (18F-NaF) to segment bones from PET data and support 5-class hybrid PET/MR-driven AC for 18F-NaF and 18F-Fluorodeoxyglucose (18F-FDG) PET/MR cardiovascular imaging.

METHODS

We introduce 5-class Ki/MR-AC for (i) 18F-NaF studies where the bones are segmented from Patlak Ki images and added as the 5th tissue class to the MR Dixon 4-class AC map. Furthermore, we propose two alternative dual-tracer protocols to permit 5-class Ki/MR-AC for (ii) 18F-FDG-only data, with a streamlined simultaneous administration of 18F-FDG and 18F-NaF at 4:1 ratio (R4:1), or (iii) for 18F-FDG-only or both 18F-FDG and 18F-NaF dual-tracer data, by administering 18F-NaF 90 minutes after an equal 18F-FDG dosage (R1:1). The Ki-driven bone segmentation was validated against computed tomography (CT)-based segmentation in rabbits, followed by PET/MR validation on 108 vertebral bone and carotid wall regions in 16 human volunteers with and without prior indication of carotid atherosclerosis disease (CAD).

RESULTS

In rabbits, we observed similar (< 1.2% mean difference) vertebral bone 18F-NaF SUVmean scores when applying 5-class AC with Ki-segmented bone (5-class Ki/CT-AC) vs CT-segmented bone (5-class CT-AC) tissue. Considering the PET data corrected with continuous CT-AC maps as gold-standard, the percentage SUVmean bias was reduced by 17.6% (18F-NaF) and 15.4% (R4:1) with 5-class Ki/CT-AC vs 4-class CT-AC. In humans without prior CAD indication, we reported 17.7% and 20% higher 18F-NaF target-to-background ratio (TBR) at carotid bifurcations wall and vertebral bones, respectively, with 5- vs 4-class AC. In the R4:1 human cohort, the mean 18F-FDG:18F-NaF TBR increased by 12.2% at carotid bifurcations wall and 19.9% at vertebral bones. For the R1:1 cohort of subjects without CAD indication, mean TBR increased by 15.3% (18F-FDG) and 15.5% (18F-NaF) at carotid bifurcations and 21.6% (18F-FDG) and 22.5% (18F-NaF) at vertebral bones. Similar TBR enhancements were observed when applying the proposed AC method to human subjects with prior CAD indication.

CONCLUSIONS

Ki-driven bone segmentation and 5-class hybrid PET/MR-driven AC is feasible and can significantly enhance 18F-NaF and 18F-FDG contrast and quantification in bone tissues and carotid walls.

Robust nonlinear parameter estimation in tracer kinetic analysis using infinity norm regularization and particle swarm optimization

In positron emission tomography (PET) studies, the voxel-wise calculation of individual rate constants describing the tracer kinetics is quite challenging because of the nonlinear relationship between the rate constants and PET data and the high noise level in voxel data. Based on preliminary simulations using a standard two-tissue compartment model, we can hypothesize that it is possible to reduce errors in the rate constant estimates when constraining the overestimation of the larger of two exponents in the model equation. We thus propose a novel approach based on infinity-norm regularization for limiting this exponent. Owing to the non-smooth cost function of this regularization scheme, which prevents the use of conventional Jacobian-based optimization methods, we examined a proximal gradient algorithm and the particle swarm optimization (PSO) through a simulation study. Because it exploits multiple initial values, the PSO method shows much better convergence than the proximal gradient algorithm, which is susceptible to the initial values. In the implementation of PSO, the use of a Gamma distribution to govern random movements was shown to improve the convergence rate and stability compared to a uniform distribution. Consequently, Gamma-based PSO with regularization was shown to outperform all other methods tested, including the conventional basis function method and Levenberg-Marquardt algorithm, in terms of its statistical properties.

Probabilistic Graphical Models for Dynamic PET: A Novel Approach to Direct Parametric Map Estimation and Image Reconstruction

In the context of dynamic emission tomography, the conventional processing pipeline consists of independent image reconstruction of single-time frames, followed by the application of a suitable kinetic model to time-activity curves (TACs) at the voxel or region-of-interest level. Direct 4D positron emission tomography (PET) reconstruction, by contrast, seeks to move beyond this scheme and incorporate information from multiple time frames within the reconstruction task. Established direct methods are based on a deterministic description of voxelwise TACs, captured by the chosen kinetic model, considering the photon counting process the only source of uncertainty. In this paper, we introduce a new probabilistic modeling strategy based on the key assumption that activity time course would be subject to uncertainty even if the parameters of the underlying dynamic process are known. This leads to a hierarchical model that we formulate using the formalism of probabilistic graphical modeling. The inference is addressed using a new iterative algorithm, in which kinetic modeling results are treated as prior expectation of activity time course, rather than as a deterministic match, making it possible to control the trade-off between a data-driven and a model-driven reconstruction. The proposed method is flexible to an arbitrary choice of (linear and nonlinear) kinetic models, it enables the inclusion of arbitrary (sub)differentiable priors for parametric maps, and it is simple to implement. Computer simulations and an application to a real-patient scan show how the proposed method is able to generalize over conventional indirect and direct approaches, providing a bridge between them by properly tuning the impact of the kinetic modeling step on image reconstruction.

Direct Parametric Maps Estimation from Dynamic PET Data: An Iterated Conditional Modes Approach

We propose and test a novel approach for direct parametric image reconstruction of dynamic PET data. We present a theoretical description of the problem of PET direct parametric maps estimation as an inference problem, from a probabilistic point of view, and we derive a simple iterative algorithm, based on the Iterated Conditional Mode (ICM) framework, which exploits the simplicity of a two-step optimization and the efficiency of an analytic method for estimating kinetic parameters from a nonlinear compartmental model. The resulting method is general enough to be flexible to an arbitrary choice of the kinetic model, and unlike many other solutions, it is capable to deal with nonlinear compartmental models without the need for linearization. We tested its performance on a two-tissue compartment model, including an analytical solution to the kinetic parameters evaluation, based on an auxiliary parameter set, with the aim of reducing computation errors and approximations. The new method is tested on simulated and clinical data. Simulation analysis led to the conclusion that the proposed algorithm gives a good estimation of the kinetic parameters in any noise condition. Furthermore, the application of the proposed method to clinical data gave promising results for further studies.

Indirect methods for improving parameter estimation of PET kinetic models

PURPOSE

Parametric images obtained from kinetic modeling of dynamic positron emission tomography (PET) data provide a new way of visualizing quantitative parameters of the tracer kinetics. However, due to the high noise level in pixel-wise image-driven time-activity curves, parametric images often suffer from poor quality and accuracy. In this study, we propose an indirect parameter estimation framework which aims to improve the quality and quantitative accuracy of parametric images.

METHODS

Three different approaches related to noise reduction and advanced curve fitting algorithm are used in the proposed framework. First, dynamic PET images are denoised using a kernel-based denoising method and the highly constrained backprojection technique. Second, gradient-free curve fitting algorithms are exploited to improve the accuracy and precision of parameter estimates. Third, a kernel-based post-filtering method is applied to parametric images to further improve the quality of parametric images. Computer simulations were performed to evaluate the performance of the proposed framework.

RESULTS AND CONCLUSIONS

The simulation results showed that when compared to the Gaussian filtering, the proposed denoising method could provide better PET image quality, and consequentially improve the quality and quantitative accuracy of parametric images. In addition, gradient-free optimization algorithms (i.e., pattern search) can result in better parametric images than the gradient-based curve fitting algorithm (i.e., trust-region-reflective). Finally, our results showed that the proposed kernel-based post-filtering method could further improve the precision of parameter estimates while maintaining the accuracy of parameter estimates.

Comparative assessment of linear least-squares, nonlinear least-squares, and Patlak graphical method for regional and local quantitative tracer kinetic modeling in cerebral dynamic 18 F-FDG PET

PURPOSE

Dynamic ¹⁸ F-FDG PET allows quantitative estimation of cerebral glucose metabolism both at the regional and local (voxel) level. Although sensitive to noise and highly computationally expensive, nonlinear least-squares (NLS) optimization stands as the reference approach for the estimation of the kinetic model parameters. Nevertheless, faster techniques, including linear least-squares (LLS) and Patlak graphical method, have been proposed to deal with high resolution noisy data, representing a more adaptable solution for routine clinical implementation. Former research investigating the relative performance of the available algorithms lack precise evaluation of kinetic parameter estimates under realistic acquisition conditions.

METHODS

The present study aims at the systematic comparison of the feasibility and pertinence of kinetic modeling of dynamic cerebral ¹⁸ F-FDG PET using NLS, LLS, and Patlak method, based on numerical simulations and patient data. Numerical simulations were used to study the bias and variance of K₁ and K_i parameters estimation under representative noise levels. Patient data allowed to assess the concordance between the three methods at the regional and voxel scale, and to evaluate the robustness of the estimations with respect to patient head motion.

RESULTS AND CONCLUSIONS

Our findings indicate that at the regional level NLS and LLS provide kinetic parameter estimates (K₁ and K_i ) with similar bias and variance characteristics (K₁ bias ± relative standard deviation [RSD] 0.0 ± 5.1% and 0.1% ± 4.9% for NLS and LLS respectively, K_i bias ± RSD 0.1% ± 4.5% and -0.7% ± 4.4% for NLS and LLS respectively). NLS estimates appear, however, to be slightly less sensitive to patient motion. At the voxel level, provided that patient motion is negligible or corrected, LLS offers an appealing alternative solution for local K₁ mapping. It yields K1 estimates that are highly correlated, with high correlation with NLS values (Pearson's r = 0.95 on actual data) within computations times less than two orders of magnitude lower. Last, Patlak method appears as the most robust and accurate technique for the estimation of K_i values at the regional and voxel scale, with or without head motion. It provides low bias/low variance K_i quantification (bias ± RSD -1.5 ± 9.5% and -4.1 ± 19.7% for Patlak and NLS respectively) as well as smooth parametric images suitable for visual assessment.

Dynamic whole-body PET imaging: principles, potentials and applications

PURPOSE

In this article, we discuss dynamic whole-body (DWB) positron emission tomography (PET) as an imaging tool with significant clinical potential, in relation to conventional standard uptake value (SUV) imaging.

BACKGROUND

DWB PET involves dynamic data acquisition over an extended axial range, capturing tracer kinetic information that is not available with conventional static acquisition protocols. The method can be performed within reasonable clinical imaging times, and enables generation of multiple types of PET images with complementary information in a single imaging session. Importantly, DWB PET can be used to produce multi-parametric images of (i) Patlak slope (influx rate) and (ii) intercept (referred to sometimes as "distribution volume"), while also providing (iii) a conventional 'SUV-equivalent' image for certain protocols.

RESULTS

We provide an overview of ongoing efforts (primarily focused on FDG PET) and discuss potential clinically relevant applications.

CONCLUSION

Overall, the framework of DWB imaging [applicable to both PET/CT(computed tomography) and PET/MRI (magnetic resonance imaging)] generates quantitative measures that may add significant value to conventional SUV image-derived measures, with limited pitfalls as we also discuss in this work.

Relative Patlak plot for dynamic PET parametric imaging without the need for early-time input function

The Patlak graphical method is widely used in parametric imaging for modeling irreversible radiotracer kinetics in dynamic PET. The net influx rate of radiotracer can be determined from the slope of the Patlak plot. The implementation of the standard Patlak method requires the knowledge of full-time input function from the injection time until the scan end time, which presents a challenge for use in the clinic. This paper proposes a new relative Patlak plot method that does not require early-time input function and therefore can be more efficient for parametric imaging. Theoretical analysis proves that the effect of early-time input function is a constant scaling factor on the Patlak slope estimation. Thus, the parametric image of the slope of the relative Patlak plot is related to the parametric image of standard Patlak slope by a global scaling factor. This theoretical finding has been further demonstrated by computer simulation and real patient data. The study indicates that parametric imaging of the relative Patlak slope can be used as a substitute of parametric imaging of standard Patlak slope for tasks that do not require absolute quantification, such as lesion detection and tumor volume segmentation.

Whole tumor kinetics analysis of 18F-fluoromisonidazole dynamic PET scans of non-small cell lung cancer patients, and correlations with perfusion CT blood flow

BACKGROUND

To determine the relative abilities of compartment models to describe time-courses of 18F-fluoromisonidazole (FMISO) tumor uptake in patients with advanced stage non-small cell lung cancer (NSCLC) imaged using dynamic positron emission tomography (dPET), and study correlations between values of the blood flow-related parameter K1 obtained from fits of the models and an independent blood flow measure obtained from perfusion CT (pCT). NSCLC patients had a 45-min dynamic FMISO PET/CT scan followed by two static PET/CT acquisitions at 2 and 4-h post-injection. Perfusion CT scanning was then performed consisting of a 45-s cine CT. Reversible and irreversible two-, three- and four-tissue compartment models were fitted to 30 time-activity-curves (TACs) obtained for 15 whole tumor structures in 9 patients, each imaged twice. Descriptions of the TACs provided by the models were compared using the Akaike and Bayesian information criteria (AIC and BIC) and leave-one-out cross-validation. The precision with which fitted model parameters estimated ground-truth uptake kinetics was determined using statistical simulation techniques. Blood flow from pCT was correlated with K1 from PET kinetic models in addition to FMISO uptake levels.

RESULTS

An irreversible three-tissue compartment model provided the best description of whole tumor FMISO uptake time-courses according to AIC, BIC, and cross-validation scores totaled across the TACs. The simulation study indicated that this model also provided more precise estimates of FMISO uptake kinetics than other two- and three-tissue models. The K1 values obtained from fits of the irreversible three-tissue model correlated strongly with independent blood flow measurements obtained from pCT (Pearson r coefficient = 0.81). The correlation from the irreversible three-tissue model (r = 0.81) was stronger than that from than K1 values obtained from fits of a two-tissue compartment model (r = 0.68), or FMISO uptake levels in static images taken at time-points from tracer injection through to 4 h later (maximum at 2 min, r = 0.70).

CONCLUSIONS

Time-courses of whole tumor FMISO uptake by advanced stage NSCLC are described best by an irreversible three-tissue compartment model. The K1 values obtained from fits of the irreversible three-tissue model correlated strongly with independent blood flow measurements obtained from perfusion CT (r = 0.81).

Accelerated PET kinetic maps estimation by analytic fitting method

In this work, we propose and test a new approach for non-linear kinetic parameters' estimation from dynamic PET data. A technique is discussed, to derive an analytical closed-form expression of the compartmental model used for kinetic parameters' evaluation, using an auxiliary parameter set, with the aim of reducing the computational burden and speeding up the fitting of these complex mathematical expressions to noisy TACs. Two alternative algorithms based on numeric calculations are considered and compared to the new proposal. We perform a simulation study aimed at (i) assessing agreement between the proposed method and other conventional ways of implementing compartmental model fitting, and (ii) quantifying the reduction in computational time required for convergence. It results in a speed-up factor of ∼120 when compared to a fully numeric version, or ∼38, with respect to a more conventional implementation, while converging to very similar values for the estimated model parameters. The proposed method is also tested on dynamic 3D PET clinical data of four control subjects. The results obtained supported those of the simulation study, and provided input and promising perspectives for the application of the proposed technique in clinical practice.