Patlak Images Research Articles

Ultrashort Oncologic Whole-Body [18F]FDG Patlak Imaging Using LAFOV PET.

Methods to shorten [18F]FDG Patlak PET imaging procedures ranging from 65-90 to 20-30 min after injection, using a population-averaged input function (PIF) scaled to patient-specific image-derived input function (IDIF) values, were recently evaluated. The aim of the present study was to explore the feasibility of ultrashort 10-min [18F]FDG Patlak imaging at 55-65 min after injection using a PIF combined with direct Patlak reconstructions to provide reliable quantitative accuracy of lung tumor uptake, compared with a full-duration 65-min acquisition using an IDIF. Methods: Patients underwent a 65-min dynamic PET acquisition on a long-axial-field-of-view (LAFOV) Biograph Vision Quadra PET/CT scanner. Subsequently, direct Patlak reconstructions and image-based (with reconstructed dynamic images) Patlak analyses were performed using both the IDIF (time to relative kinetic equilibrium between blood and tissue concentration (t*) = 30 min) and a scaled PIF at 30-60 min after injection. Next, direct Patlak reconstructions were performed on the system console using only the last 10 min of the acquisition, that is, from 55 to 65 min after injection, and a scaled PIF using maximum crystal ring difference settings of both 85 and 322. Tumor lesion and healthy-tissue uptake was quantified and compared between the differently obtained parametric images to assess quantitative accuracy. Results: Good agreement was obtained between direct- and image-based Patlak analyses using the IDIF (t* = 30 min) and scaled PIF at 30-60 min after injection, performed using the different approaches, with no more than 8.8% deviation in tumor influx rate value (Ki ) (mean difference ranging from -0.0022 to 0.0018 mL/[min × g]). When direct Patlak reconstruction was performed on the system console, excellent agreement was found between the use of a scaled PIF at 30-60 min after injection versus 55-65 min after injection, with 2.4% deviation in tumor Ki (median difference, -0.0018 mL/[min × g]; range, -0.0047 to 0.0036 mL/[min × g]). For different maximum crystal ring difference settings using the scan time interval of 55-65 min after injection, only a 0.5% difference (median difference, 0.0000 mL/[min × g]; range, -0.0004 to 0.0013 mL/[min × g]) in tumor Ki was found. Conclusion: Ultrashort whole-body [18F]FDG Patlak imaging is feasible on an LAFOV Biograph Vision Quadra PET/CT system without loss of quantitative accuracy to assess lung tumor uptake compared with a full-duration 65-min acquisition. The ultrashort 10-min direct Patlak reconstruction with PIF allows for its implementation in clinical practice.

Impact of different parametric Patlak imaging approaches and comparison with a 2-tissue compartment pharmacokinetic model with a long axial field-of-view (LAFOV) PET/CT in oncological patients.

The recently introduced Long-Axial-Field-of-View (LAFOV) PET-CT scanners allow for the first-time whole-body dynamic- and parametric imaging. Primary aim of this study was the comparison of direct and indirect Patlak imaging as well as the comparison of different time frames for Patlak calculation with the LAFOV PET-CT in oncological patients. Secondary aims of the study were lesion detectability and comparison of Patlak analysis with a two-tissue-compartment model (2TCM). 50 oncological patients with 346 tumor lesions were enrolled in the study. All patients underwent [18F]FDG PET/CT (skull to upper thigh). Here, the Image-Derived-Input-Function) (IDIF) from the descending aorta was used as the exclusive input function. Four sets of images have been reviewed visually and evaluated quantitatively using the target-to-background (TBR) and contrast-to-noise ratio (CNR): short-time (30min)-direct (STD) Patlak Ki, short-time (30min)-indirect (STI) Patlak Ki, long-time (59.25min)-indirect (LTI) Patlak Ki, and 50-60min SUV (sumSUV). VOI-based 2TCM was used for the evaluation of tumor lesions and normal tissues and compared with the results of Patlak model. No significant differences were observed between the four approaches regarding the number of tumor lesions. However, we found three discordant results: a true positive liver lesion in all Patlak Ki images, a false positive liver lesion delineated only in LTI Ki which was a hemangioma according to MRI and a true negative example in a patient with an atelectasis next to a lung tumor. STD, STI and LTI Ki images had superior TBR in comparison with sumSUV images (2.9-, 3.3- and 4.3-fold higher respectively). TBR of LTI Ki were significantly higher than STD Ki. VOI-based k3 showed a 21-fold higher TBR than sumSUV. Parameters of different models vary in their differential capability between tumor lesions and normal tissue like Patlak Ki which was better in normal lung and 2TCM k3 which was better in normal liver. 2TCM Ki revealed the highest correlation (r = 0.95) with the LTI Patlak Ki in tumor lesions group and demonstrated the highest correlation with the STD Patlak Ki in all tissues group and normal tissues group (r = 0.93 and r = 0.74 respectively). Dynamic [18F]-FDG with the new LAFOV PET/CT scanner produces Patlak Ki images with better lesion contrast than SUV images, but does not increase the lesion detection rate. The time window used for Patlak imaging plays a more important role than the direct or indirect method. A combination of different models, like Patlak and 2TCM may be helpful in parametric imaging to obtain the best TBR in the whole body in future.

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

BackgroundThe 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.MethodsThirty-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).ResultsThe 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).ConclusionOur 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.

Prognostic value of whole-body dynamic 18F-FDG PET/CT Patlak in diffuse large B-cell lymphoma

ObjectiveThis study aims to investigate the significance of interim whole-body dynamic 18F-fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) Patlak parameters for predicting the prognosis of patients with diffuse large B-cell lymphoma. To estimate the predictive value of the whole-body dynamic 18F-FDG PET/CT Patlak parameter for 2-year progression-free survival (PFS) and 2-year overall survival (OS). MethodsThis study reports the findings of 67 patients with diffuse large B-cell lymphoma (DLBCL). These patients underwent interim whole-body dynamic 18F-FDG PET/CT scans from June 2021 to January 2023 at the Department of Nuclear Medicine, First Affiliated Hospital of Anhui Medical University. The predictive values of maximum standard uptake value (SUVmax), maximum of net glucose uptake rate (Kimax) and the predictive model combining Kimax and interim treatment response on the prognosis of patients was analyzed using receiver operating characteristic (ROC) curves. Kaplan-Meier survival curves and log-rank tests were used for survival analysis. Univariate and multivariate analyses were performed to screen for independent prognostic risk factors. ResultsAfter a median follow-up of 18 months, 21 patients (31.3%) experienced disease recurrence or death. The cut-off values for the SUVmax and the Kimax were 6.1 and 0.13 μmol min−1·ml−1, respectively. Ann Arbor stage, IPI, SUVmax, Kimax and interim treatment response were associated with PFS and OS in the univariate analysis. However, only Kimax and interim treatment response were independent influences on PFS and OS in multivariate analysis. ConclusionInterim whole-body dynamic 18F-FDG PET/CT Patlak imaging has significant prognostic value in patients with DLBCL. Among them, the interim dynamic parameter Kimax showed the best predictive value for prognosis compared with the interim SUVmax and interim treatment response. The predictive model established by Kimax and the interim treatment response allowed for the accurate stratification of the prognostic risk of DLBCL.

Short-time-window Patlak imaging using a population-based arterial input function and optimized Bayesian penalized likelihood reconstruction: a feasibility study

BackgroundTo explore the feasibility of short-time-window Ki imaging using a population-based arterial input function (IF) and optimized Bayesian penalized likelihood (BPL) reconstruction as a practical alternative to long-time-window Ki imaging with an individual patient-based IF. Myocardial Ki images were generated from 73 dynamic 18F-FDG-PET/CT scans of 30 patients with cardiac sarcoidosis. For each dynamic scan, the Ki images were obtained using the IF from each individual patient and a long time window (10–60 min). In addition, Ki images were obtained using the normalized averaged population-based IF and BPL algorithms with different beta values (350, 700, and 1000) with a short time window (40–60 min). The visual quality of each image was visually rated using a 4-point scale (0, not visible; 1, poor; 2, moderate; and 3, good), and the Ki parameters (Ki-max, Ki-mean, Ki-volume) of positive myocardial lesions were measured independently by two readers. Wilcoxon’s rank sum test, McNemar’s test, or linear regression analysis were performed to assess the differences or relationships between two quantitative variables.ResultsBoth readers similarly rated 51 scans as positive (scores = 1–3) and 22 scans as negative (score = 0) for all four Ki images. Among the three types of population-based IF Ki images, the proportion of images with scores of 3 was highest with a beta of 1000 (78.4 and 72.5%, respectively) and lowest with a beta of 350 (33.3 and 23.5%) for both readers (all p < 0.001). The coefficients of determination between the Ki parameters obtained with the individual patient-based IF and those obtained with the population-based IF were highest with a beta of 1000 for both readers (Ki-max, 0.91 and 0.92, respectively; Ki-mean, 0.91 and 0.92, respectively; Ki-volume, 0.75 and 0.60, respectively; and all p < 0.001).ConclusionsShort-time-window Ki images with a population-based IF reconstructed using the BPL algorithm and a high beta value were closely correlated with long-time-window Ki images generated with an individual patient-based IF. Short-time-window Ki images using a population-based IF and BPL reconstruction might represent practical alternatives to long-time-window Ki images generated using an individual patient-based IF.

Rapid high-quality PET Patlak parametric image generation based on direct reconstruction and temporal nonlocal neural network

Parametric imaging based on dynamic positron emission tomography (PET) has wide applications in neurology. Compared to indirect methods, direct reconstruction methods, which reconstruct parametric images directly from the raw PET data, have superior image quality due to better noise modeling and richer information extracted from the PET raw data. For low-dose scenarios, the advantages of direct methods are more obvious. However, the wide adoption of direct reconstruction is inevitably impeded by the excessive computational demand and deficiency of the accessible raw data. In addition, motion modeling inside dynamic PET image reconstruction raises more computational challenges for direct reconstruction methods. In this work, we focused on the 18F-FDG Patlak model, and proposed a data-driven approach which can estimate the motion corrected full-dose direct Patlak images from the dynamic PET reconstruction series, based on a proposed novel temporal non-local convolutional neural network. During network training, direct reconstruction with motion correction based on full-dose dynamic PET sinograms was performed to obtain the training labels. The reconstructed full-dose /low-dose dynamic PET images were supplied as the network input. In addition, a temporal non-local block based on the dynamic PET images was proposed to better recover the structural information and reduce the image noise. During testing, the proposed network can directly output high-quality Patlak parametric images from the full-dose /low-dose dynamic PET images in seconds. Experiments based on 15 full-dose and 15 low-dose 18F-FDG brain datasets were conducted and analyzed to validate the feasibility of the proposed framework. Results show that the proposed framework can generate better image quality than reference methods.

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

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. Clinical data from 24 patients referred for tumor staging were included in this study. The patients underwent a whole-body dynamic PET study, 20min after FDG injection (0.13mCi/kg). The proposed whole-body scanning protocol includes 6 passes with 4-5 bed positions, depending on the size of the patient, with 2min 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. 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. 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

BackgroundArterial 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.MethodsThe 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.ResultsThe 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.ConclusionsBoth 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.

Generalized whole-body Patlak parametric imaging for enhanced quantification in clinical PET

We recently developed a dynamic multi-bed PET data acquisition framework to translate the quantitative benefits of Patlak voxel-wise analysis to the domain of routine clinical whole-body (WB) imaging. The standard Patlak (sPatlak) linear graphical analysis assumes irreversible PET tracer uptake, ignoring the effect of FDG dephosphorylation, which has been suggested by a number of PET studies. In this work: (i) a non-linear generalized Patlak (gPatlak) model is utilized, including a net efflux rate constant kloss, and (ii) a hybrid (s/g)Patlak (hPatlak) imaging technique is introduced to enhance contrast to noise ratios (CNRs) of uptake rate Ki images. Representative set of kinetic parameter values and the XCAT phantom were employed to generate realistic 4D simulation PET data, and the proposed methods were additionally evaluated on 11 WB dynamic PET patient studies. Quantitative analysis on the simulated Ki images over 2 groups of regions-of-interest (ROIs), with low (ROI A) or high (ROI B) true kloss relative to Ki, suggested superior accuracy for gPatlak. Bias of sPatlak was found to be 16–18% and 20–40% poorer than gPatlak for ROIs A and B, respectively. By contrast, gPatlak exhibited, on average, 10% higher noise than sPatlak. Meanwhile, the bias and noise levels for hPatlak always ranged between the other two methods. In general, hPatlak was seen to outperform all methods in terms of target-to-background ratio (TBR) and CNR for all ROIs. Validation on patient datasets demonstrated clinical feasibility for all Patlak methods, while TBR and CNR evaluations confirmed our simulation findings, and suggested presence of non-negligible kloss reversibility in clinical data. As such, we recommend gPatlak for highly quantitative imaging tasks, while, for tasks emphasizing lesion detectability (e.g. TBR, CNR) over quantification, or for high levels of noise, hPatlak is instead preferred. Finally, gPatlak and hPatlak CNR was systematically higher compared to routine SUV values.

The 2D Hotelling filter - a quantitative noise-reducing principal-component filter for dynamic PET data, with applications in patient dose reduction

BackgroundIn this paper we apply the principal-component analysis filter (Hotelling filter) to reduce noise from dynamic positron-emission tomography (PET) patient data, for a number of different radio-tracer molecules. We furthermore show how preprocessing images with this filter improves parametric images created from such dynamic sequence.We use zero-mean unit variance normalization, prior to performing a Hotelling filter on the slices of a dynamic time-series. The Scree-plot technique was used to determine which principal components to be rejected in the filter process. This filter was applied to [11C]-acetate on heart and head-neck tumors, [18F]-FDG on liver tumors and brain, and [11C]-Raclopride on brain. Simulations of blood and tissue regions with noise properties matched to real PET data, was used to analyze how quantitation and resolution is affected by the Hotelling filter. Summing varying parts of a 90-frame [18F]-FDG brain scan, we created 9-frame dynamic scans with image statistics comparable to 20 MBq, 60 MBq and 200 MBq injected activity. Hotelling filter performed on slices (2D) and on volumes (3D) were compared.ResultsThe 2D Hotelling filter reduces noise in the tissue uptake drastically, so that it becomes simple to manually pick out regions-of-interest from noisy data. 2D Hotelling filter introduces less bias than 3D Hotelling filter in focal Raclopride uptake. Simulations show that the Hotelling filter is sensitive to typical blood peak in PET prior to tissue uptake have commenced, introducing a negative bias in early tissue uptake. Quantitation on real dynamic data is reliable. Two examples clearly show that pre-filtering the dynamic sequence with the Hotelling filter prior to Patlak-slope calculations gives clearly improved parametric image quality. We also show that a dramatic dose reduction can be achieved for Patlak slope images without changing image quality or quantitation.ConclusionsThe 2D Hotelling-filtering of dynamic PET data is a computer-efficient method that gives visually improved differentiation of different tissues, which we have observed improve manual or automated region-of-interest delineation of dynamic data. Parametric Patlak images on Hotelling-filtered data display improved clarity, compared to non-filtered Patlak slope images without measurable loss of quantitation, and allow a dramatic decrease in patient injected dose.

Measurement of metabolic tumor volume: static versus dynamic FDG scans

BackgroundMetabolic tumor volume assessment using positron-emission tomography [PET] may be of interest for both target volume definition in radiotherapy and monitoring response to therapy. It has been reported, however, that metabolic volumes derived from images of metabolic rate of glucose (generated using Patlak analysis) are smaller than those derived from standardized uptake value [SUV] images. The purpose of this study was to systematically compare metabolic tumor volume assessments derived from SUV and Patlak images using a variety of (semi-)automatic tumor delineation methods in order to identify methods that can be used reliably on (whole body) SUV images.MethodsDynamic [18F]-fluoro-2-deoxy-D-glucose [FDG] PET data from 10 lung and 8 gastrointestinal cancer patients were analyzed retrospectively. Metabolic tumor volumes were derived from both Patlak and SUV images using five different types of tumor delineation methods, based on various thresholds or on a gradient.ResultsIn general, most tumor delineation methods provided more outliers when metabolic volumes were derived from SUV images rather than Patlak images. Only gradient-based methods showed more outliers for Patlak-based tumor delineation. Median measured metabolic volumes derived from SUV images were larger than those derived from Patlak images (up to 59% difference) when using a fixed percentage threshold method. Tumor volumes agreed reasonably well (< 26% difference) when applying methods that take local signal-to-background ratio [SBR] into account.ConclusionLarge differences may exist in metabolic volumes derived from static and dynamic FDG image data. These differences depend strongly on the delineation method used. Delineation methods that correct for local SBR provide the most consistent results between SUV and Patlak images.

FDG–PET parametric imaging by total variation minimization

Parametric imaging of the cerebral metabolic rate for glucose (CMRGlc) using [ 18F]-fluoro deoxyglucose positron emission tomography is considered. Traditional imaging is hindered due to low signal-to-noise ratios at individual voxels. We propose to minimize the total variation of the tracer uptake rates while requiring good fit of traditional Patlak equations. This minimization guarantees spatial homogeneity within brain regions and good distinction between brain regions. Brain phantom simulations demonstrate significant improvement in quality of images by the proposed method as compared to Patlak images with post-filtering using Gaussian or median filters.

CT derived Patlak images of the human kidney.

This study aimed to produce Patlak images of the kidney from dynamic CT data and to determine whether such images are substantially affected by fluid movement between renal tubular segments. Renal permeability was measured in 31 kidneys by applying Patlak analysis to time-density data from kidney and aorta during dynamic CT. Permeability parameters were correlated against plasma urea. The renal region (cortex or medulla) with the greatest permeability was determined from parametric images generated using pixel by pixel analysis. The mean value for whole kidney permeability was 517.5 microliters min-1 ml-1. A correlation was found between whole kidney permeability and plasma urea (p < 0.01). Permeability values were highest in the renal medulla in 24 (77%) kidneys. The higher medullary values of permeability are artefactual, resulting from movement of fluid and contrast medium between cortex and medulla. Although Patlak images do not reflect true intrarenal permeability values, the apparent medullary permeability may provide diagnostically useful information about the concentrating ability of the kidney. CT measurements of whole kidney permeability reflect filtration function but the apparent intrarenal variations in permeability will result in measurement errors dependent upon the relative amounts of renal cortex and medulla included in the CT slice studied.