4394

A simple phantom for extracting Intravoxel Incoherent Motion (IVIM) model parameters
Erick O Buko1, Afis Ajala1, Jiming Zhang2, and Pei Herng Hor1
1Physics, Texas Center for Superconductivity at University of Houston, Houston, TX, United States, 2Diagnostic and Interventional Radiology, Baylor St. Luke's Medical Center, Houston, TX, United States

Synopsis

We present a simple phantom set-up to model diffusion weighted imaging signal arising from intra-voxel incoherent motion (IVIM). The model provides means to independently control both the perfusion fraction (f) as well as the pseudodiffusion coefficient (Df).

Synopsis

We present a simple phantom set-up to model diffusion weighted imaging signal arising from intra-voxel incoherent motion (IVIM). The model provides means to independently control both the perfusion fraction (f) as well as the pseudodiffusion coefficient (Df). Across a range of flow velocities, numerical simulations and phantom experiments confirm that: (a) perfusion fraction (f) is overestimated by traditional analysis approaches such as segmented (S), or oversegmented (OS) methods, but is slightly underestimated by the AS-T2 method; (b) Df increases in proportion to flow velocity, and the estimation of Df is unaffected by the analysis technique.

Introduction

Since the initial description of IVIM phenomena by LeBihan1, several groups have suggested various approaches to create phantom models to simulate physiologic conditions 2-4. This has been a challenging issue as tissue diffusion (Dt) and fluid diffusion (Df) vary by two orders of magnitude or more, and creation of physical models with well-defined perfusion volume fractions (f) is also difficult. In this work, we present a simple phantom model that can be used to synthesize signals mimicking IVIM phenomena which can be used to evaluate the performance of various analysis algorithms in the estimation of parameters such as f, and Df.

METHODS

IVIM Phantom: A cartoon representation of the phantom model used to generate IVIM signal is shown in Figure 1. A small cylinder of known diameter filled with material of desired diffusivity, and desired T2 relaxation time serves as a representation of static tissue (Dt, and T2t), and fluid with desired T2 (T2f) is pumped into another cylindrical container with a spout at a specific height which is connected to a straw with an opening at the bottom of the container, acts as a representation of fluid compartment. The same fluid is pumped into the fluid compartment using a peristaltic pump at a constant flow rate, and a steady state is established when the inflow rate equals the outflow rate at the spout which flows into the reservoir of the pump. Adjusting the flow rate alters the fluid velocity (V) within the fluid compartment, thereby altering Df.
Signal generation: This phantom set up can now be imaged using diffusion weighted imaging with various ‘b’ values, with the static tissue compartment and the fluid compartment within the imaging field of view, and when the flow across the fluid compartment is in steady state.
Controlling perfusion fraction (f) :Now depending on the perfusion fraction desired, the ratio of pixels from the static and fluid compartments can be randomly picked, e.g., if a perfusion fraction of 20% is desired, the ratio of pixels chosen from the static tissue compartment and flowing compartment will be 4:1. The addition of these measured signals will yield an IVIM signal curve that includes the effect of the diffusivities and relaxation rate differences between the two compartments (Figure 2).
Controlling pseudo-diffusion coefficient (Df ): To attain a pre-determined value of Df without the confounding effects of relaxation rate differences, the diffusion weighted signal just from the flowing compartment is fitted with the biexponential IVIM model to extract the value of Df as a function of flow velocity which can be estimated from the input flow rate from the peristaltic pump and the geometry of the fluid compartment (cylinder, in this instance).
Numerical Simulations: MR signal was simulated using Eq. 1 with Df as in Eq. 6. The following were fixed to mimic our phantom: f 0.25, Dt 0.002 mm2/s, l 0.7354 mm, T2f 97.3 ms, T2t 63.5 ms, 10 b values (0,10,25,50,75,100,200,300,500,800 s/mm2). We also varied V from 0 to 0.921 mm/s and 8 TE ranging from 0 to 150 ms. Gaussian noise was added to generate MR signal for each sample at SNR levels of 20 and 50. Each sample was simulated 5000 times.
Phantom: Two phantoms with different T2 values, 63.5 ms and 97.3 ms were made using water doped with Gd to represent tissue and fluid compartments respectively. Fluid compartment was passed through a bottle of diameter 48 mm at 5 different flow rates (0, 40, 60, 80, 100 mL/min) resulting into velocities ranging from 0 to 0.921 mm/s.
MR acquisition : Diffusion Weighted Imaging (DWI) signals were acquired from a 1.5 T scanner using the following parameters: FOV 319 mm X 165 mm , Voxel size 3 mm X 3 mm X 10 mm, TR 2000ms, TE 60 ms and 10 b values as used in the simulation.

RESULTS

(1) Numerical simulations show that IVIM analysis methods (S, OS, AS) which do not address the T2 differences between the tissue and fluid compartments overestimate f by as much as 40% at a fixed TE of 60 ms, but with T2 correction, correct value of f can be obtained at all fluid velocities (Figure 3). These numerical simulations are verified from the experimental measurements at various flow velocities (Figure 4). (2) The estimated Df as a function of V for simulated and measured signals is shown in Figure 5.

CONCLUSION

The proposed phantom model is flexible to generate arbitrary signals with different combinations of f and Df, and is easy to include the effect of tissue relaxation rate differences between the compartments.

Acknowledgements

No acknowledgement found.

References

1. Le Bihan, D., Breton, E., Lallemand, D., Aubin, M. L., Vignaud, J., & Laval-Jeantet, M. (1988). Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology, 168(2), 497-505.

2. Schneider, M., Gaaß, T., Dinkel, J., Ingrisch, M., Reiser, M., & Dietrich, O. (2016). Intravoxel incoherent motion MRI in a 3‐dimensional microvascular flow phantom. In Proc Int Soc Magn Reson Med (Vol. 24, p. 0920).

3. Cho, G. Y., Kim, S., Jensen, J. H., Storey, P., Sodickson, D. K., & Sigmund, E. E. (2012). A versatile flow phantom for intravoxel incoherent motion MRI. Magnetic resonance in medicine, 67(6), 1710-1720.

4. Lorenz, C. H., Pickens III, D. R., Puffer, D. B., & Price, R. R. (1991). Magnetic resonance diffusion/perfusion phantom experiments. Magnetic resonance in medicine, 19(2), 254-260.

Figures

Figure 1: A cartoon representation of the phantom model used to generate IVIM signal. Tissue compartment is a non flowing phantom (Top left) while fluid was flowing at a steady flow rate in the fluid compartment (Top right).

Figure 2: A representation of how IVIM signal was achieved by mixing signals from tissue and fluid compartments with predefined fluid fraction.

Figure 3: Mean and STD of the estimated f from simulated data as a function of V using four methods (S, OS, AS and AS-T2) at different SNR levels a) 20 and b)50 at TE of 60 ms. Ffluid is the true vale of f.

Figure 4: a) Estimated f as a function of V at a TE of 60 ms using four methods (S, OS, AS and AS-T2) . Ffluid is the total signal fraction from the fluid while fflow (expected f) is the product of f measured from pure fluid and Ffluid. b) Percentage differences between the estimated f and fflow at different V values.

Figure 5: Estimated Df as a function of V at a TE of 60 ms using four methods (S, OS, AS and AS-T2). Df was measured from pure fluid compartment using S method from simulated data (left) and experimental data (right) .

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
4394