Arterial input function estimation compensating for inflow and partial voluming in dynamic contrast-enhanced MRI.

Both inflow and the partial volume effect (PVE) are sources of error when measuring the arterial input function (AIF) in dynamic contrast-enhanced (DCE) MRI. This is relevant, as errors in the AIF can propagate into pharmacokinetic parameter estimations from the DCE data. A method was introduced for flow correction by estimating and compensating the number of the perceived pulse of spins during inflow. We hypothesized that the PVE has an impact on concentration-time curves similar to inflow. Therefore, we aimed to study the efficiency of this method to compensate for both effects simultaneously. We first simulated an AIF with different levels of inflow and PVE contamination. The peak, full width at half-maximum (FWHM), and area under curve (AUC) of the reconstructed AIFs were compared with the true (simulated) AIF. In clinical data, the PVE was included in AIFs artificially by averaging the signal in voxels surrounding a manually selected point in an artery. Subsequently, the artificial partial volume AIFs were corrected and compared with the AIF from the selected point. Additionally, corrected AIFs from the internal carotid artery (ICA), the middle cerebral artery (MCA), and the venous output function (VOF) estimated from the superior sagittal sinus (SSS) were compared. As such, we aimed to investigate the effectiveness of the correction method with different levels of inflow and PVE in clinical data. The simulation data demonstrated that the corrected AIFs had only marginal bias in peak value, FWHM, and AUC. Also, the algorithm yielded highly correlated reconstructed curves over increasingly larger neighbourhoods surrounding selected arterial points in clinical data. Furthermore, AIFs measured from the ICA and MCA produced similar peak height and FWHM, whereas a significantly larger peak and lower FWHM was found compared with the VOF. Our findings indicate that the proposed method has high potential to compensate for PVE and inflow simultaneously. The corrected AIFs could thereby provide a stable input source for DCE analysis.

Improved reliability of perfusion estimation in dynamic susceptibility contrast MRI by using the arterial input function from dynamic contrast enhanced MRI.

The arterial input function (AIF) plays a crucial role in estimating quantitative perfusion properties from dynamic susceptibility contrast (DSC) MRI. An important issue, however, is that measuring the AIF in absolute contrast-agent concentrations is challenging, due to uncertainty in relation to the measured R 2 ∗ -weighted signal, signal depletion at high concentration, and partial-volume effects. A potential solution could be to derive the AIF from separately acquired dynamic contrast enhanced (DCE) MRI data. We aim to compare the AIF determined from DCE MRI with the AIF from DSC MRI, and estimated perfusion coefficients derived from DSC data using a DCE-driven AIF with perfusion coefficients determined using a DSC-based AIF. AIFs were manually selected in branches of the middle cerebral artery (MCA) in both DCE and DSC data in each patient. In addition, a semi-automatic AIF-selection algorithm was applied to the DSC data. The amplitude and full width at half-maximum of the AIFs were compared statistically using the Wilcoxon rank-sum test, applying a 0.05 significance level. Cerebral blood flow (CBF) was derived with different AIF approaches and compared further. The results showed that the AIFs extracted from DSC scans yielded highly variable peaks across arteries within the same patient. The semi-automatic DSC-AIF had significantly narrower width compared with the manual AIFs, and a significantly larger peak than the manual DSC-AIF. Additionally, the DCE-based AIF provided a more stable measurement of relative CBF and absolute CBF values estimated with DCE-AIFs that were compatible with previously reported values. In conclusion, DCE-based AIFs were reproduced significantly better across vessels, showed more realistic profiles, and delivered more stable and reasonable CBF measurements. The DCE-AIF can, therefore, be considered as an alternative AIF source for quantitative perfusion estimations in DSC MRI.

Validation of cardiac image-derived input functions for functional PET quantification.

PURPOSE: Functional PET (fPET) is a novel technique for studying dynamic changes in brain metabolism and neurotransmitter signaling. Accurate quantification of fPET relies on measuring the arterial input function (AIF), traditionally achieved through invasive arterial blood sampling. While non-invasive image-derived input functions (IDIF) offer an alternative, they suffer from limited spatial resolution and field of view. To overcome these issues, we developed and validated a scan protocol for brain fPET utilizing cardiac IDIF, aiming to mitigate known IDIF limitations. METHODS: Twenty healthy individuals underwent fPET/MR scans using [18F]FDG or 6-[18F]FDOPA, utilizing bed motion shuttling to capture cardiac IDIF and brain task-induced changes. Arterial and venous blood sampling was used to validate IDIFs. Participants performed a monetary incentive delay task. IDIFs from various blood pools and composites estimated from a linear fit over all IDIF blood pools (3VOI) and further supplemented with venous blood samples (3VOIVB) were compared to the AIF. Quantitative task-specific images from both tracers were compared to assess the performance of each input function to the gold standard. RESULTS: For both radiotracer cohorts, moderate to high agreement (r: 0.60-0.89) between IDIFs and AIF for both radiotracer cohorts was observed, with further improvement (r: 0.87-0.93) for composite IDIFs (3VOI and 3VOIVB). Both methods showed equivalent quantitative values and high agreement (r: 0.975-0.998) with AIF-derived measurements. CONCLUSION: Our proposed protocol enables accurate non-invasive estimation of the input function with full quantification of task-specific changes, addressing the limitations of IDIF for brain imaging by sampling larger blood pools over the thorax. These advancements increase applicability to any PET scanner and clinical research setting by reducing experimental complexity and increasing patient comfort.

Improving measurement of blood-brain barrier permeability with reduced scan time using deep-learning-derived capillary input function.

PURPOSE: In Dynamic contrast-enhanced MRI (DCE-MRI), Arterial Input Function (AIF) has been shown to be a significant contributor to uncertainty in the estimation of kinetic parameters. This study is to assess the feasibility of using a deep learning network to estimate local Capillary Input Function (CIF) to estimate blood-brain barrier (BBB) permeability, while reducing the required scan time. MATERIALS AND METHOD: A total of 13 healthy subjects (younger (<40 y/o): 8, older (> 67 y/o): 5) were recruited and underwent 25-min DCE-MRI scans. The 25 min data were retrospectively truncated to 10 min to simulate a reduced scan time of 10 min. A deep learning network was trained to predict the CIF using simulated tissue contrast dynamics with two vascular transport models. The BBB permeability (PS) was measured using 3 methods: (i) Ca-25min, using DCE-MRI data of 25 min with individually sampled AIF (Ca); (ii) Ca-10min, using truncated 10min data with AIF (Ca); and (iii) Cp-10min, using truncated 10 min data with CIF (Cp). The PS estimates from the Ca-25min method were used as reference standard values to assess the accuracy of the Ca-10min and Cp-10min methods in estimating the PS values. RESULTS: When compared to the reference method(Ca-25min), the Ca-10min and Cp-10min methods resulted in an overestimation of PS by 217 ± 241 % and 48.0 ± 30.2 %, respectively. The Bland Altman analysis showed that the mean difference from the reference was 8.85 ± 1.78 (x10-4 min-1) with the Ca-10min, while it was reduced to 1.63 ± 2.25 (x10-4 min-1) with the Cp-10min, resulting in an average reduction of 81%. The limits of agreement also reduced by up to 39.2% with the Cp-10min. We found a 75% increase of BBB permeability in the gray matter and a 35% increase in the white matter, when comparing the older group to the younger group. CONCLUSIONS: We demonstrated the feasibility of estimating the capillary-level input functions using a deep learning network. We also showed that this method can be used to estimate subtle age-related changes in BBB permeability with reduced scan time, without compromising accuracy. Moreover, the trained deep learning network can automatically select CIF, reducing the potential uncertainty resulting from manual user-intervention.

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.

Non-invasive kinetic modelling approaches for quantitative analysis of brain PET studies.

Pharmacokinetic modelling with arterial sampling is the gold standard for analysing dynamic PET data of the brain. However, the invasive character of arterial sampling prevents its widespread clinical application. Several methods have been developed to avoid arterial sampling, in particular reference region methods. Unfortunately, for some tracers or diseases, no suitable reference region can be defined. For these cases, other potentially non-invasive approaches have been proposed: (1) a population based input function (PBIF), (2) an image derived input function (IDIF), or (3) simultaneous estimation of the input function (SIME). This systematic review aims to assess the correspondence of these non-invasive methods with the gold standard. Studies comparing non-invasive pharmacokinetic modelling methods with the current gold standard methods using an input function derived from arterial blood samples were retrieved from PubMed/MEDLINE (until December 2021). Correlation measurements were extracted from the studies. The search yielded 30 studies that correlated outcome parameters (VT, DVR, or BPND for reversible tracers; Ki or CMRglu for irreversible tracers) from a potentially non-invasive method with those obtained from modelling using an arterial input function. Some studies provided similar results for PBIF, IDIF, and SIME-based methods as for modelling with an arterial input function (R2 = 0.59-1.00, R2 = 0.71-1.00, R2 = 0.56-0.96, respectively), if the non-invasive input curve was calibrated with arterial blood samples. Even when the non-invasive input curve was calibrated with venous blood samples or when no calibration was applied, moderate to good correlations were reported, especially for the IDIF and SIME (R2 = 0.71-1.00 and R2 = 0.36-0.96, respectively). Overall, this systematic review illustrates that non-invasive methods to generate an input function are still in their infancy. Yet, IDIF and SIME performed well, not only with arterial blood calibration, but also with venous or no blood calibration, especially for some tracers without plasma metabolites, which would potentially make these methods better suited for clinical application. However, these methods should still be properly validated for each individual tracer and application before implementation.

Influence of Different Measurement Methods of Arterial Input Function on Quantitative Dynamic Contrast-Enhanced MRI Parameters in Head and Neck Cancer.

BACKGROUND: Head and neck cancer (HNC) is the sixth most prevalent cancer worldwide. Dynamic contrast-enhanced MRI (DCE-MRI) helps in diagnosis and prognosis. Quantitative DCE-MRI requires an arterial input function (AIF), which affects the values of pharmacokinetic parameters (PKP). PURPOSE: To evaluate influence of four individual AIF measurement methods on quantitative DCE-MRI parameters values (Ktrans , ve , kep , and vp ), for HNC and muscle. STUDY TYPE: Prospective. POPULATION: A total of 34 HNC patients (23 males, 11 females, age range 24-91) FIELD STRENGTH/SEQUENCE: A 3 T; 3D SPGR gradient echo sequence with partial saturation of inflowing spins. ASSESSMENT: Four AIF methods were applied: automatic AIF (AIFa) with up to 50 voxels selected from the whole FOV, manual AIF (AIFm) with four voxels selected from the internal carotid artery, both conditions without (Mc-) or with (Mc+) motion correction. Comparison endpoints were peak AIF values, PKP values in tumor and muscle, and tumor/muscle PKP ratios. STATISTICAL TESTS: Nonparametric Friedman test for multiple comparisons. Nonparametric Wilcoxon test, without and with Benjamini Hochberg correction, for pairwise comparison of AIF peak values and PKP values for tumor, muscle and tumor/muscle ratio, P value ≤ 0.05 was considered statistically significant. RESULTS: Peak AIF values differed significantly for all AIF methods, with mean AIFmMc+ peaks being up to 66.4% higher than those for AIFaMc+. Almost all PKP values were significantly higher for AIFa in both, tumor and muscle, up to 76% for mean Ktrans values. Motion correction effect was smaller. Considering tumor/muscle parameter ratios, most differences were not significant (0.068 ≤ Wilcoxon P value ≤ 0.8). DATA CONCLUSION: We observed important differences in PKP values when using either AIFa or AIFm, consequently choice of a standardized AIF method is mandatory for DCE-MRI on HNC. From the study findings, AIFm and inflow compensation are recommended. The use of the tumor/muscle PKP ratio should be of interest for multicenter studies. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

A simplified method to correct saturation of arterial input function for cardiac magnetic resonance first-pass perfusion imaging: validation with simultaneously acquired PET.

BACKGROUND: First-pass perfusion imaging in magnetic resonance imaging (MRI) is an established method to measure myocardial blood flow (MBF). An obstacle for accurate quantification of MBF is the saturation of blood pool signal intensity used for arterial input function (AIF). The objective of this project was to validate a new simplified method for AIF estimation obtained from single-bolus and single sequence perfusion measurements. The reference MBF was measured simultaneously on 13N-ammonia positron emission tomography (PET). METHODS: Sixteen patients with clinically confirmed myocardial ischemia were imaged in a clinical whole-body PET-MRI system. PET perfusion imaging was performed in a 10-min acquisition after the injection of 10 mCi of 13N-ammonia. The MRI perfusion acquisition started simultaneously with the start of the PET acquisition after the injection of a 0.075 mmol/kg gadolinium contrast agent. Cardiac stress imaging was initiated after the administration of regadenoson 20 s prior to PET-MRI scanning. The saturation part of the MRI AIF data was modeled as a gamma variate curve, which was then estimated for a true AIF by minimizing a cost function according to various boundary conditions. A standard AHA 16-segment model was used for comparative analysis of absolute MBF from PET and MRI. RESULTS: Overall, there were 256 segments in 16 patients, mean resting perfusion for PET was 1.06 ± 0.34 ml/min/g and 1.04 ± 0.30 ml/min/g for MRI (P = 0.05), whereas mean stress perfusion for PET was 2.00 ± 0.74 ml/min/g and 2.12 ± 0.76 ml/min/g for MRI (P < 0.01). Linear regression analysis in MBF revealed strong correlation (r = 0.91, slope = 0.96, P < 0.001) between PET and MRI. Myocardial perfusion reserve, calculated from the ratio of stress MBF over resting MBF, also showed a strong correlation between MRI and PET measurements (r = 0.82, slope = 0.81, P < 0.001). CONCLUSION: The results demonstrated the feasibility of the simplified AIF estimation method for the accurate quantification of MBF by MRI with single sequence and single contrast injection. The MRI MBF correlated strongly with PET MBF obtained simultaneously. This post-processing technique will allow easy transformation of clinical perfusion imaging data into quantitative information.

Arterial input functions in dynamic susceptibility contrast MRI (DSC-MRI) in longitudinal evaluation of brain metastases.

BACKGROUND: Dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) could be helpful to separate true disease progression from pseudo-progression in brain metastases when assessing the need for retreatment. However, the selection of arterial input functions (AIFs) is not standardized for analysis, limiting its use for this application. PURPOSE: To compare population-based AIFs, AIFs specific to each patient, and AIFs specific to every visit in the longitudinal follow-up of brain metastases. MATERIAL AND METHODS: Longitudinal data were collected from eight patients before treatment (6 of 8 patients) and after treatment (6-17 visits). Imaging was performed using a 1.5-T MRI system. Lesions were segmented by subtracting precontrast images from postcontrast images. Cerebral blood volume (rCBV) and cerebral blood flow (rCBF) were computed, and Pearson's product moment correlation coefficients were calculated to evaluate similarity of DSC parameters dependent on various AIF choices across time. AIF shape characteristics were compared. Parameter differences between white matter (WM) and gray matter (GM) were obtained to determine which AIF choice maximizes tissue differentiation. RESULTS: Although DSC parameters follow similar patterns in time, the various AIF selections cause large parameter variations with relative standard deviations of up to ±60%. AIFs sampled in one patient across sessions more similar in shape than AIFs sampled across patients. Estimates of rCBV based on scan-specific AIFs differentiated better between perfusion in WM and GM than patient-specific or population-based AIFs (P ≤ 0.02). CONCLUSION: Results indicate that scan-specific AIFs are the best choice for DSC-MRI parameter estimations in the longitudinal follow-up of brain metastases.

Estimation of the capillary level input function for dynamic contrast-enhanced MRI of the breast using a deep learning approach.

PURPOSE: To develop a deep learning approach to estimate the local capillary-level input function (CIF) for pharmacokinetic model analysis of DCE-MRI. METHODS: A deep convolutional network was trained with numerically simulated data to estimate the CIF. The trained network was tested using simulated lesion data and used to estimate voxel-wise CIF for pharmacokinetic model analysis of breast DCE-MRI data using an abbreviated protocol from women with malignant (n = 25) and benign (n = 28) lesions. The estimated parameters were used to build a logistic regression model to detect the malignancy. RESULT: The pharmacokinetic parameters estimated using the network-predicted CIF from our breast DCE data showed significant differences between the malignant and benign groups for all parameters. Testing the diagnostic performance with the estimated parameters, the conventional approach with arterial input function (AIF) showed an area under the curve (AUC) between 0.76 and 0.87, and the proposed approach with CIF demonstrated similar performance with an AUC between 0.79 and 0.81. CONCLUSION: This study shows the feasibility of estimating voxel-wise CIF using a deep neural network. The proposed approach could eliminate the need to measure AIF manually without compromising the diagnostic performance to detect the malignancy in the clinical setting.

Optimal saturation recovery time for minimizing the underestimation of arterial input function in quantitative cardiac perfusion MRI.

PURPOSE: The purpose of this study was to determine an optimal saturation-recovery time (TS) for minimizing the underestimation of arterial input function (AIF) in quantitative cardiac perfusion MRI without multiple gadolinium injections per subject. METHODS: We scanned 18 subjects (mean age = 59 ± 14 years, 9/9 males/females) to acquire resting perfusion data and 1 additional subject (age = 38 years, male) to obtain stress-rest perfusion data using a 5-fold accelerated pulse sequence with radial k-space sampling and applied k-space weighted image contrast (KWIC) filters on the same k-space data to retrospectively reconstruct five AIF images with effective TS ranging from 10 to 21.2 ms (2.8 ms steps). Undersampled images were reconstructed using a compressed sensing framework with temporal-total-variation and temporal-principal-component as 2 orthogonal sparsifying transforms. The image processing steps included, same motion correction across five different AIF images, signal normalization by the proton-density-weighted-image, signal-to-T1 conversion using a Bloch equation, T1 -to-gadolinium-concentration conversion assuming fast water exchange, T2 * correction to the AIF, and gadolinium-concentration to myocardial blood flow (MBF) conversion based on a Fermi model. RESULTS: Among five TS values, the shortest TS (10 ms) produced significantly (P < 0.05) higher peak AIF and lower resting MBF (13.73 mM, 0.73 mL g-1 min-1 ) than 12.8 ms (11.24 mM, 0.89 mL g-1 min-1 ), 15.6 ms (9.56 mM, 1.05 mL g-1 min-1 ), 18.4 ms (8.55 mM, 1.17 mL g-1 min-1 ), and 21.2 ms (7.95 mM, 1.27 mL g-1 min-1 ). Similarly, shorter TS reduced underestimation of AIF (or overestimation of MBF) for both during stress and at rest, but this effect was canceled in myocardial-perfusion-reserve (MPR). CONCLUSION: This study demonstrates that TS of 10 ms reduces the underestimation of AIF and, hence, the overestimation of MBF compared with longer TS values (12.8-21.2 ms).

Dependency of R2 and R2 * relaxation on Gd-DTPA concentration in arterial blood: Influence of hematocrit and magnetic field strength.

Dynamic susceptibility contrast (DSC) MRI is clinically used to measure brain perfusion by monitoring the dynamic passage of a bolus of contrast agent through the brain. For quantitative analysis of the DSC images, the arterial input function is required. It is known that the original assumption of a linear relation between the R2(*) relaxation and the arterial contrast agent concentration is invalid, although the exact relation is as of yet unknown. Studying this relation in vitro is time-consuming, because of the widespread variations in field strengths, MRI sequences, contrast agents, and physiological conditions. This study aims to simulate the R2(*) versus contrast concentration relation under varying physiological and technical conditions using an adapted version of an open-source simulation tool. The approach was validated with previously acquired data in human whole blood at 1.5 T by means of a gradient-echo sequence (proof-of-concept). Subsequently, the impact of hematocrit, field strength, and oxygen saturation on this relation was studied for both gradient-echo and spin-echo sequences. The results show that for both gradient-echo and spin-echo sequences, the relaxivity increases with hematocrit and field strength, while the hematocrit dependency was nonlinear for both types of MRI sequences. By contrast, oxygen saturation has only a minor effect. In conclusion, the simulation setup has proven to be an efficient method to rapidly calibrate and estimate the relation between R2(*) and gadolinium concentration in whole blood. This knowledge will be useful in future clinical work to more accurately retrieve quantitative information on brain perfusion.

Dynamic contrast-enhanced MRI parametric mapping using high spatiotemporal resolution Golden-angle RAdial Sparse Parallel MRI and iterative joint estimation of the arterial input function and pharmacokinetic parameters.

The aim of this work is to develop a data-driven quantitative dynamic contrast-enhanced (DCE) MRI technique using Golden-angle RAdial Sparse Parallel (GRASP) MRI with high spatial resolution and high flexible temporal resolution and pharmacokinetic (PK) analysis with an arterial input function (AIF) estimated directly from the data obtained from each patient. DCE-MRI was performed on 13 patients with gynecological malignancy using a 3-T MRI scanner with a single continuous golden-angle stack-of-stars acquisition and image reconstruction with two temporal resolutions, by exploiting a unique feature in GRASP that reconstructs acquired data with user-defined temporal resolution. Joint estimation of the AIF (both AIF shape and delay) and PK parameters was performed with an iterative algorithm that alternates between AIF and PK estimation. Computer simulations were performed to determine the accuracy (expressed as percentage error [PE]) and precision of the estimated parameters. PK parameters (volume transfer constant [Ktrans ], fractional volume of the extravascular extracellular space [ve ], and blood plasma volume fraction [vp ]) and normalized root-mean-square error [nRMSE] (%) of the fitting errors for the tumor contrast kinetic data were measured both with population-averaged and data-driven AIFs. On patient data, the Wilcoxon signed-rank test was performed to compare nRMSE. Simulations demonstrated that GRASP image reconstruction with a temporal resolution of 1 s/frame for AIF estimation and 5 s/frame for PK analysis resulted in an absolute PE of less than 5% in the estimation of Ktrans and ve , and less than 11% in the estimation of vp . The nRMSE (mean ± SD) for the dual temporal resolution image reconstruction and data-driven AIF was 0.16 ± 0.04 compared with 0.27 ± 0.10 (p < 0.001) with 1 s/frame using population-averaged AIF, and 0.23 ± 0.07 with 5 s/frame using population-averaged AIF (p < 0.001). We conclude that DCE-MRI data acquired and reconstructed with the GRASP technique at dual temporal resolution can successfully be applied to jointly estimate the AIF and PK parameters from a single acquisition resulting in data-driven AIFs and voxelwise PK parametric maps.

Quantitative assessment of myelin density using [11C]MeDAS PET in patients with multiple sclerosis: a first-in-human study.

PURPOSE: Multiple sclerosis (MS) is a disease characterized by inflammatory demyelinated lesions. New treatment strategies are being developed to stimulate myelin repair. Quantitative myelin imaging could facilitate these developments. This first-in-man study aimed to evaluate [11C]MeDAS as a PET tracer for myelin imaging in humans. METHODS: Six healthy controls and 11 MS patients underwent MRI and dynamic [11C]MeDAS PET scanning with arterial sampling. Lesion detection and classification were performed on MRI. [11C]MeDAS time-activity curves of brain regions and MS lesions were fitted with various compartment models for the identification of the best model to describe [11C]MeDAS kinetics. Several simplified methods were compared to the optimal compartment model. RESULTS: Visual analysis of the fits of [11C]MeDAS time-activity curves showed no preference for irreversible (2T3k) or reversible (2T4k) two-tissue compartment model. Both volume of distribution and binding potential estimates showed a high degree of variability. As this was not the case for 2T3k-derived net influx rate (Ki), the 2T3k model was selected as the model of choice. Simplified methods, such as SUV and MLAIR2 correlated well with 2T3k-derived Ki, but SUV showed subject-dependent bias when compared to 2T3k. Both the 2T3k model and the simplified methods were able to differentiate not only between gray and white matter, but also between lesions with different myelin densities. CONCLUSION: [11C]MeDAS PET can be used for quantification of myelin density in MS patients and is able to distinguish differences in myelin density within MS lesions. The 2T3k model is the optimal compartment model and MLAIR2 is the best simplified method for quantification. TRIAL REGISTRATION: NL7262. Registered 18 September 2018.

Operator dependency of arterial input function in dynamic contrast-enhanced MRI.

OBJECTIVE: To investigate the effect of inter-operator variability in arterial input function (AIF) definition on kinetic parameter estimates (KPEs) from dynamic contrast-enhanced (DCE) MRI in patients with high-grade gliomas. METHODS: The study included 118 DCE series from 23 patients. AIFs were measured by three domain experts (DEs), and a population AIF (pop-AIF) was constructed from the measured AIFs. The DE-AIFs, pop-AIF and AUC-normalized DE-AIFs were used for pharmacokinetic analysis with the extended Tofts model. AIF-dependence of KPEs was assessed by intraclass correlation coefficient (ICC) analysis, and the impact on relative longitudinal change in Ktrans was assessed by Fleiss' kappa (κ). RESULTS: There was a moderate to substantial agreement (ICC 0.51-0.76) between KPEs when using DE-AIFs, while AUC-normalized AIFs yielded ICC 0.77-0.95 for Ktrans, kep and ve and ICC 0.70 for vp. Inclusion of the pop-AIF did not reduce agreement. Agreement in relative longitudinal change in Ktrans was moderate (κ = 0.591) using DE-AIFs, while AUC-normalized AIFs gave substantial (κ = 0.809) agreement. DISCUSSION: AUC-normalized AIFs can reduce the variation in kinetic parameter results originating from operator input. The pop-AIF presented in this work may be applied in absence of a satisfactory measurement.

Maximum Entropy Technique and Regularization Functional for Determining the Pharmacokinetic Parameters in DCE-MRI.

This paper aims to solve the arterial input function (AIF) determination in dynamic contrast-enhanced MRI (DCE-MRI), an important linear ill-posed inverse problem, using the maximum entropy technique (MET) and regularization functionals. In addition, estimating the pharmacokinetic parameters from a DCE-MR image investigations is an urgent need to obtain the precise information about the AIF-the concentration of the contrast agent on the left ventricular blood pool measured over time. For this reason, the main idea is to show how to find a unique solution of linear system of equations generally in the form of [Formula: see text] named an ill-conditioned linear system of equations after discretization of the integral equations, which appear in different tomographic image restoration and reconstruction issues. Here, a new algorithm is described to estimate an appropriate probability distribution function for AIF according to the MET and regularization functionals for the contrast agent concentration when applying Bayesian estimation approach to estimate two different pharmacokinetic parameters. Moreover, by using the proposed approach when analyzing simulated and real datasets of the breast tumors according to pharmacokinetic factors, it indicates that using Bayesian inference-that infer the uncertainties of the computed solutions, and specific knowledge of the noise and errors-combined with the regularization functional of the maximum entropy problem, improved the convergence behavior and led to more consistent morphological and functional statistics and results. Finally, in comparison to the proposed exponential distribution based on MET and Newton's method, or Weibull distribution via the MET and teaching-learning-based optimization (MET/TLBO) in the previous studies, the family of Gamma and Erlang distributions estimated by the new algorithm are more appropriate and robust AIFs.

Modified Maximum Entropy Method and Estimating the AIF via DCE-MRI Data Analysis.

Background: For the kinetic models used in contrast-based medical imaging, the assignment of the arterial input function named AIF is essential for the estimation of the physiological parameters of the tissue via solving an optimization problem. Objective: In the current study, we estimate the AIF relayed on the modified maximum entropy method. The effectiveness of several numerical methods to determine kinetic parameters and the AIF is evaluated-in situations where enough information about the AIF is not available. The purpose of this study is to identify an appropriate method for estimating this function. Materials and Methods: The modified algorithm is a mixture of the maximum entropy approach with an optimization method, named the teaching-learning method. In here, we applied this algorithm in a Bayesian framework to estimate the kinetic parameters when specifying the unique form of the AIF by the maximum entropy method. We assessed the proficiency of the proposed method for assigning the kinetic parameters in the dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), when determining AIF with some other parameter-estimation methods and a standard fixed AIF method. A previously analyzed dataset consisting of contrast agent concentrations in tissue and plasma was used. Results and Conclusions: We compared the accuracy of the results for the estimated parameters obtained from the MMEM with those of the empirical method, maximum likelihood method, moment matching ("method of moments"), the least-square method, the modified maximum likelihood approach, and our previous work. Since the current algorithm does not have the problem of starting point in the parameter estimation phase, it could find the best and nearest model to the empirical model of data, and therefore, the results indicated the Weibull distribution as an appropriate and robust AIF and also illustrated the power and effectiveness of the proposed method to estimate the kinetic parameters.

Robust arterial input function surrogate measurement from the superior sagittal sinus complex signal for fast dynamic contrast-enhanced MRI in the brain.

PURPOSE: Accurately estimating the arterial input function for dynamic contrast-enhanced MRI is challenging. An arterial input function is typically determined from signal magnitude changes related to a contrast agent, often leading to underestimation of peak concentrations. Alternatively, signal phase recovers the accurate peak concentration for straight vessels but suffers from high noise. A recent method proposed to fit the signal in the complex plane by combining the advantages of the previous 2 methods. The purpose of this work is to refine this complex-based method to determine the venous output function (VOF), an arterial input function surrogate, from the superior sagittal sinus. METHODS: We propose a state-of-the-art complex-based method that includes direct compensation for blood inflow and signal phase correction accounting for the curvature of the superior sagittal sinus, generally assumed collinear with B0 . We compared the magnitude-, phase-, and complex-based VOF determination methods against various simulated biases as well as for 29 brain metastases patients. RESULTS: Angulation of the superior sagittal sinus relative to B0 varied widely within patients, and its effect on the signal phase caused an underestimation of peak concentrations of up to 65%. Correction significantly increased the VOF peak concentration for the phase- and complex-based VOFs in the cohort. The phase-based method recovered accurate peak concentrations but lacked precision in the tail of the VOF. Our complex-based VOF completely recovered the effect of inflow and resulted in a high-peak concentration with limited noise. CONCLUSION: The new complex-based method resulted in high-quality VOF robust against superior sagittal sinus curvature and variations in patient positioning.

A theoretical framework for retrospective T2∗ correction to the arterial input function in quantitative myocardial perfusion MRI.

PURPOSE: To develop and evaluate a flexible, Bloch-equation based framework for retrospective T2∗ correction to the arterial input function (AIF) obtained with quantitative cardiac perfusion pulse sequences. METHODS: Our framework initially calculates the gadolinium concentration [Gd] based on T1 measurements alone. Next, T2∗ is estimated from this initial calculation of [Gd] while assuming fast water exchange and using the literature native T2 and static magnetic field variation (ΔB0 ) values. Finally, the [Gd] is recalculated after performing T2∗ correction to the Bloch equation signal model. Using this approach, we performed T2∗ correction to historical phantom and in vivo, dual-imaging perfusion data sets from 3 different patient groups obtained using different pulse sequences and imaging parameters. Images were processed to quantify both the AIF and resting myocardial blood flow (MBF). We also performed a sensitivity analysis of our T2∗ correction to ±20% variations in native T2 and ΔB0 . RESULTS: Compared with the ground truth [Gd] of phantom, the normalized root-means-square-error (NRMSE) in measured [Gd] was 5.1%, 1.3%, and 0.6% for uncorrected, our corrected, and Kellman's corrected, respectively. For in vivo data, both the peak AIF (7.0 ± 3.0 mM vs. 8.6 ± 7.1 mM, 7.2 ± 0.9 mM vs. 8.6 ± 1.7 mM, 7.7 ± 1.8 mM vs. 10.3 ± 5.1 mM, P < .001) and resting MBF (1.3 ± 0.1 mL/g/min vs. 1.1 ± 0.1 mL/g/min, 1.3 ± 0.1 mL/g/min vs. 1.1 ± 0.1 mL/g/min, 1.2 ± 0.1 mL/g/min vs. 0.9 ± 0.1 mL/g/min, P < .001) values were significantly different between uncorrected and corrected for all 3 patient groups. Both the peak AIF and resting MBF values varied by <5% over the said variations in native T2 and ΔB0 . CONCLUSION: Our theoretical framework enables retrospective T2∗ correction to the AIF obtained with dual-imaging, cardiac perfusion pulse sequences.

Short-term repeatability of myocardial blood flow using 82Rb PET/CT: The effect of arterial input function position and motion correction.

BACKGROUND: We tested the repeatability of myocardial blood flow (MBF) quantified using 82Rb with and without motion correction (MC) and with arterial input functions estimated from left ventricle (LV) and atrium (LA). METHODS: Twenty-one patients referred for clinical 82Rb PET/CT underwent repeated rest scans in a single imaging session. Global MBF was quantified using three different assessments by two operators: (1) automatic processing without MC and LV arterial input function (AIF), (2) with MC and LV-AIF, and (3) with MC and LA-AIF. Inter-scan and inter-operator repeatability were tested using coefficient of variation (CV). RESULTS: MC with LV-AIF did not change MBF (no MC: 1.01 ± 0.30 mL/min/g vs MC with LV-AIF: 1.01 ± 0.29, P = 0.70), whereas MC with LA-AIF showed significantly lower MBF assessments (0.95 ± 0.28 mL/min/g, P = 0.0006). We report significant improvement for test-retest reproducibility for global MBF following MC (CV; No MC: 16.0, MC (LV-AIF): 9.2, MC (LA-AIF): 8.8). Good inter-operator repeatability was observed for LV-AIF (CV = 4.7) and LA-AIF (CV = 5.6) for global MBF assessments. CONCLUSIONS: MC significantly improved the test-retest repeatability between operators and between scans. MBF obtained after MC with LV-AIF were comparable, whereas MBFs after MC and LA-AIF were significantly reduced.