With the introduction of combined PET/MRI systems, AIF conversion can be made under certain circumstances (see [1]). We propose a model that allows modification of the injection parameters in the AIF fit to account for differences caused by different injection durations [2].
Brain 18F-Choline PET and DSC-MRI data were obtained using Siemens mMR. The MR contrast agent was injected with a rate of 4ml/sec and the PET tracer was injected manually. Perfusion Mismatch Analyzer [3] was used to extract the MRI-AIF. Carotid arteries were segmented on a post contrast MPRAGE image. PET frames were registered onto this MPRAGE image using rigid registration and partial volume correction was done using the iterative Yang method [4]. The AIFs were fitted using a convolution of a ‘double Butterworth’ function, representing the injection, with a tri-exponential function representing the elimination [Eq. 1]. The bolus shape can be adjusted by changing Δτ (τ2 - τ1). This was tested with a population based MRI AIF [5], as well as with clinical data.
1
where
For the population based input function, Figure Figure11 shows that when Δτ was increased, lower and wider peaks were seen, and with decreased Δτ, higher but narrower peaks were observed. Figure Figure22 shows that the function fits both clinical PET and MRI AIFs well. Values of τ1 and τ2 were changed to modify the MRI-AIF and Figure Figure33 shows the modified MRI-AIF together with the original fitted PET-AIF, normalized to their peaks. Two AIFs have similar peak shapes but start to differ at the elimination phase as Gd-DOTA and 18F-Choline have different tissue uptake rates.
Figure 1
Simulated MRI-AIFs using Parker’s population-based input function refitted with the developed function. AIF shapes with different injection durations, Δτ is shown.
Figure 2
The double Butterworth convolution function used to fit (a) DSC-MRI data and (b) 18F-Choline PET data together with a plot where the timescale of PET-AIF was limited to MRI-AIF’s to show different bolus widths.
Figure 3
The MRI-AIF with modified τ1 and τ2 values plotted together with the PET-AIF. The MRI-AIF peak is scaled to PET-AIF’s peak.
This enables conversion of the early part of the AIFs from one modality to another even if different injection protocols are used.
[1]
Luc Tremblay,et al.
Conversion of arterial input functions for dual pharmacokinetic modeling using Gd‐DTPA/MRI and 18F‐FDG/PET
,
2013,
Magnetic resonance in medicine.
[2]
Uulke A. van der Heide,et al.
The impact of overall injection time on the arterial input function and pharmacokinetic analysis using the Tofts model in DCE-MRI for prostate cancer patients
,
2011
.
[3]
Tobias J. Hagge,et al.
Physics
,
1929,
Nature.
[4]
A. Jackson,et al.
Experimentally‐derived functional form for a population‐averaged high‐temporal‐resolution arterial input function for dynamic contrast‐enhanced MRI
,
2006,
Magnetic resonance in medicine.
[5]
I. Buvat,et al.
A review of partial volume correction techniques for emission tomography and their applications in neurology, cardiology and oncology
,
2012,
Physics in medicine and biology.