A New High-Resolution Perfusion-Weighted Imaging Technique without exogenous tracer and Its Quantitative Model
Hyunseok Seo1, Dongchan Kim1, Jaejin Cho1, Kinam kwon1, Byungjai Kim1, and HyunWook Park1

1Electrical Engineering, KAIST, Daejeon, Korea, Republic of

Synopsis

The proposed method is based on the perfusion model of intravoxel incoherent motion (IVIM), which assumes that perfusion is a microcirculation of blood in the capillary network. The method introduces two-types of bi-polar gradients in a radial spin-echo sequence, a new ordering of bi-polar gradients for isotropic perfusion weighting, and signal models for quantitative analysis of perfusion. To verify the proposed method, in-vivo imaging in 3 T MRI was performed, and the proposed method offers the high-resolution quantitative perfusion map.

Introduction

The perfusion-weighted imaging (PWI) offers a good indicator to measure the metabolism of the biologic tissues. So it has various applications for the clinical diagnosis and for the neuroscience. There are two typical methods for PWI. One is the dynamic susceptibility contrast-MRI (DSC-MRI) and the other is arterial spin labeling (ASL). While the DSC-MRI method offers the high contrast-to-noise ratio (CNR), it is used reluctantly because it requires an invasively exogenous tracer (1). Alternatively, the non-invasive ASL technique can be more widely applied, but it still has the major drawback of a poor signal-to-noise ratio (SNR) (2). In this abstract, a new high-resolution PWI technique without the exogenous is proposed.

Methods

The proposed method obtains perfusion information by separating diffusion and perfusion information based on the concept of IVIM (3) using the radial spin-echo sequence. When a single bi-polar gradient is applied, as shown in Fig. 1(a), the MR signal of a voxel is modeled by the exponential decay for diffusion and the decay of a sine integral function for perfusion, which considers the laminar flow in the capillary (4). However, the double bi-polar gradients in Fig. 1(a) can compensate the perfusion effects, but not the diffusion effects. Thus, the signal intensity from the single (S single) and double (S double) bi-polar gradients can be respectively described as follows,

\[S^{\tt single}=S_{0}\exp(-bD){\tt Si}(2\alpha v)/(2\alpha v)\\ S^{\tt double}=S_{0}\exp(-bD),\\{\tt where}~{\tt Si}(x)=\int_{0}^{x}\frac{\sin\tau}{\tau}d\tau,{\tt and}~\alpha=|\int_{t}^{}G(\tau)d\tau|.\tag1\]

The S0 is the MR signal independent of diffusion and perfusion terms, D is the apparent diffusion coefficients (ADC), and v is voxel-averaged velocity of the capillaries, which is directly proportional to cerebral blood flow (CBF) and cerebral blood volume (CBV). Since the cerebral perfusion is anisotropic (5), the signal decay from perfusion is affected by the direction of the bi-polar gradients. So the direction of bi-polar gradients is changed for each projection view (radial spoke) to obtain an averaged perfusion weighting (6, 7).

\[S^{\tt single}=\frac{S_{0}}{N}\sum_{n=1}^N\exp(-bD_{n}){\tt Si} (2\alpha v_{n})/(2\alpha v_{n})\\S^{\tt double}=\frac{S_{0}}{N}\sum_{n=1}^N\exp(-bD_{n})\tag2\]

Eq. [2] means the arithmetic averaged signal due to changing the direction of bi-polar gradient along the radial spoke, where N is the total number of bi-polar gradient directions, Dn and vn are the ADC and voxel-averaged capillary velocity, respectively, of the nth bi-polar gradient direction. Then, Eq. [2] is approximated as Eq. [3] with error rates of ε1 and ε2.

\[S^{\tt single}=S_{0}\exp(-b\frac{1}{N}\sum_{n=1}^ND_{n}){\tt si}(2\alpha\frac{1}{N}\sum_{n=1}^Nv_{n})/(2\alpha\frac{1}{N}\sum_{n=1}^Nv_{n})+S_{0}\cdot\epsilon1\\S^{\tt double}=S_{0}\exp(-b\frac{1}{N}\sum_{n=1}^ND_{n})+S_{0}\cdot\epsilon2\tag3\]

It was already proven that ε2 is negligible (7). ε1 can be estimated with the Taylor series expansion of exponential function and the series expression of sine integral function (8) with in-vivo imaging parameters. If ε1 is also negligible, the voxel-averaged velocity of capillaries, which considers the perfusion anisotropy (i.e., mean of vns), can be calculated by the inverse process of (S single / S double). In-vivo experiments were performed using a 3 T MRI scanner (Siemens Magnetom Verio, Erlangen, Germany) to validate the proposed method. Imaging parameters are given as field of view (FOV) = 256×256 mm2, matrix size = 328×328, TR/TE = 2000/52 ms, number of averages = 1, b-values = 52 s/mm2, and number of sampling points for each projection view = 328. A total of 328 projection views with 328 different bi-polar gradient directions were used.

Results

Figure 2 shows that the proposed signal model has the error rate (ε1) less than 1.2 %, when fractional anisotropies (FAs) of diffusion and perfusion do not have extreme value (i.e., FA < 0.6). Therefore, the error term in the approximated S single can be ignored. The proposed method offers the voxel-averaged capillary velocity map with an in-plane resolution of 0.78 × 0.78 mm2, as shown in Fig. 3. Figure 4 shows a fusion image that the voxel-averaged capillary velocity map is overlaid on the anatomical T2-weighted image. The proposed method reflects the general characteristic that perfusion information between the gray matter and white matter is different.

Conclusion & Discussion

As demonstrated by the simulations based on the formula and the in-vivo experiments, the proposed method offers the high-resolution quantitative perfusion map. Thus, the proposed method can be a useful PWI tool.

Acknowledgements

This research was supported by the Brain Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2014M3C7033999).

This work was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry for Health and Welfare, Korea (HI14C1135).

References

1. Magn Reson Med 1996;36:715-725.

2. Magn Reson Med 2009;62:165-173.

3. Radiology 1988;168:497-505.

4. Medical Physics 1987;14:43-18.

5. Magn Reson Med 2008;60:1284-1291.

6. 22th Annual Meeting of ISMRM 2014;1499.

7. Phys Med Biol 2014;59:6289-6303.

8. Princeton University Press 2003:105-106.

9. J Magn Reson Imaging 2014;39:624-632.

Figures

(a) The proposed pulse sequence diagram, and (b) 328 directions of bi-polar gradients used in the experiment.

The simulation results of ε1 with respect to FAs, when the perfusion fraction (f) is 0.3 determined by reference (9).

A high-resolution voxel-averaged capillary velocity map.

A fusion image that Fig. 3 is overlaid on the anatomical T2-weighted image.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3321