4268

RESOLVE for simultaneous mapping high-frequency conductivity and micro-structure parameters in conductivity tensor imaging (CTI)
Tong Sun1, Songxiong Wu2, Nashan Wu2, Qingjun Sun2, Haodong Qin3, Guangyao Wu2, Xin Chen1, and Chunqi Chang1
1Shenzhen University, Shenzhen, China, 2Radiology Department, Shenzhen University General Hospital and Shenzhen University Clinical Medical, Shenzhen University General Hospital, Shenzhen, China, 3Siemens Healthineers, Guangzhou, China, Siemens Healthineers, Shenzhen, China

Synopsis

Keywords: Quantitative Imaging, Electromagnetic Tissue Properties, Biomarkers

Motivation: Conventional two sequences acquisition in conductivity tensor imaging (CTI) can cause geometric mismatch between two acquisition data which may influence subsequent reconstruction.

Goal(s): We aim to simultaneously reconstruct the high-frequency conductivity and to fit microstructure parameters with data from single sequence.

Approach: We propose a novel data acquisition strategy which reconstructs the high-frequency conductivity from the phase of non-diffusion weighted data of RESOLVE.

Results: The reconstructed high-frequency conductivity and fitted microstructure parameters has matched geometric images, and the low-frequency conductivity can be successfully reconstructed.

Impact: We show the conventional two sequence acquisition is reduced to the single sequence acquisition which avoids the geometric mismatch and obtains higher precision of low-frequency conductivity.

Introduction

Electrical conductivity characterizes how electrical currents propagate in biological tissues. and it is influenced by the frequency, ion mobility, ion concentration, and cell shapes, and cell membranes [1]. The property of electrical conductivity may be altered in the presence of pathological conditions, such as stroke and tumors [2, 3]. Meanwhile, the low-frequency conductivity tensor of human brain can realize patient-specific volume conductor models for neuroimaging and electrical stimulation. Recently, a novel low-frequency conductivity tensor imaging (CTI) technique [4], which combined MREPT measurement with multi-b-valued DWI data had been proposed. In previous studies [5, 6], the high-frequency conductivity was reconstructed from the transceive phase acquired by fast spin-echo (FSE) sequence, while the micro-structure parameters were fitted to diffusion-weighted data acquired by the single-shot diffusion-weighted EPI sequence. However, it is sometimes difficult to guarantee matched geometric distortion between transceive phase and diffusion-weighted images. We proposed a novel data acquisition strategy which simultaneously mapped the high-frequency conductivity and micro-structure parameters with the RESOLVE sequence in CTI. With proposed method, the geometric mismatch in conventional acquisition strategy could be eliminated.

Methods

High-frequency conductivity reconstruction with cr-MREPT
The phase-based cr-MREPT formula is given in Figure 1.(1) [7], where φtr = φ+ + φ- is the transceive phase. The artificial diffusion term c is added to the convection reaction equation-based MREPT formula.
Multi-compartment microscopic diffusion imaging model
We use the standard three-compartment model [8], whose kernel functional is shown in Figure 1.(2) where fint is the intra-neurite volume fraction, fw is the free water volume fraction, Dint is the intrinsic diffusion coefficient, Δe=Dext,∥-Dext,⊥,Dext,∥ is the extra-neurite parallel diffusivity, Dext,⊥ is the extra-neurite perpendicular diffusivity. ξ2 = cos2θ. Using nonlinear least square fitting (see Figure 1.(3)), micro-structure parameters fint,fw,Dint,Dext,∥ ,Dext,⊥ can be obtained.
Low-frequency mean conductivity
According to the literature [6], the low-frequency mean conductivity can be expressed by Figure 1.(4), where Dext is the extracellular mean diffusivity, β is the ion concentration ratio of intracellular and extracellular spaces.
Low-frequency anisotropic conductivity tensor
Similar to the literature [6], the low-frequency anisotropic conductivity tensor can be expressed as by Figure 1.(5), where Dsh denotes the spherical harmonic (SH) coefficients of ODF. The overview of the proposed algorithm is shown in Figure 2.
Hollow Polypropylene fibre phantom experiment
Hollow Polypropylene yarns were used to fabricate the fibre buddle. Three fibre buddles are immersed in the 0.8% NaCl solution shown in figure 2. Imaging experiment was performed using a Siemens Prisma scanner equipped with 64-channel head-neck coil. The MR data was acquired by RESOLVE with the following parameters: TR/TE = 5640 / 44ms; in-plane FOV = 220 mm × 220 mm; voxel size = 1.7×1.7×3 mm3; number of slices =35; number of shot = 5; 1 b0 image; diffusion-weighting images with b-values of 1000, 2000 s/mm2 and 30 gradient directions per b-value. The amplitude and phase data were reconstructed from vendor retrospective reconstruction software tool.

Results

The amplitude image of b = 0, 1000 s/mm2 is shown in the middle panel of figure 3. while the phase image of b = 0 s/mm2 is shown in right panel of figure 2. High-frequency conductivity, micro-structure parameters, low-frequency mean conductivity are respectively shown in Figure 4. The intra-neurite volume fraction of fibre is slightly larger than the free water, which may lead to the lower low-frequency conductivity than high-frequency conductivity. The conductivity orientation distribution function (ODF) is shown in Figure 5.

Discussion

In the conventional method, multi-echo spin echo and single-shot EPI sequences were respectively used to acquire phase and diffusion weighted data, which may lead to geometric mismatch in the final low-frequency conductivity. Compare to the conventional method, the proposed method only utilized single multi-shot EPI sequence. The phase is reconstructed from the complex image of b0, thus the phase and amplitude images had the matched geometric distortion. Moreover, we designed a fibre phantom aimed to map the anisotropic distribution. Instead of using second order tensor, we used truncated two order SH coefficient which could reconstruct fibre crossing with four or six order.

Conclusion

The performance of a novel data acquisition method to produce the low-frequency conductivity is validated using an anisotropic fibre conductivity phantom. The proposed method has the advantages of simple post-processing and higher potential clinical application.

Acknowledgements

No acknowledgement found.

References

1 Grimnes, S., and Martinsen, Ø.G.: ‘Chapter 7 - Electrodes’, in Grimnes, S., and Martinsen, Ø.G. (Eds.): ‘Bioimpedance and Bioelectricity Basics (Third Edition)’ (Academic Press, 2015), pp. 179-254

2 Shin, J., Kim, M.J., Lee, J., Nam, Y., Kim, M.-o., Choi, N., Kim, S., and Kim, D.-H.: ‘Initial study on in vivo conductivity mapping of breast cancer using MRI’, Journal of Magnetic Resonance Imaging, 2015, 42, (2), pp. 371-378

3 Jensen-Kondering, U., Shu, L., Böhm, R., Jansen, O., and Katscher, U.: ‘In-vivo pilot study at 3 Tesla: Feasibility of Electric Properties Tomography in a rat model of stroke’, Physics in Medicine, 2020, 9, pp. 100024

4 Sajib, S.Z.K., Kwon, O.I., Kim, H.J., and Woo, E.J.: ‘Electrodeless conductivity tensor imaging (CTI) using MRI: basic theory and animal experiments’, Biomedical Engineering Letters, 2018, 8, (3), pp. 273-282

5 Katoch, N., Choi, B.K., Sajib, S.Z.K., Lee, E., Kim, H.J., Kwon, O.I., and Woo, E.J.: ‘Conductivity Tensor Imaging of In Vivo Human Brain and Experimental Validation Using Giant Vesicle Suspension’, IEEE Transactions on Medical Imaging, 2019, 38, (7), pp. 1569-1577

6 Jahng, G.-H., Lee, M.B., Kim, H.J., Je Woo, E., and Kwon, O.-I.: ‘Low-frequency dominant electrical conductivity imaging of in vivo human brain using high-frequency conductivity at Larmor-frequency and spherical mean diffusivity without external injection current’, NeuroImage, 2021, 225, pp. 117466

7 Gurler, N., and Ider, Y.Z.: ‘Gradient-based electrical conductivity imaging using MR phase’, Magnetic Resonance in Medicine, 2017, 77, (1), pp. 137-150

8 Novikov, D.S., Veraart, J., Jelescu, I.O., and Fieremans, E.: ‘Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI’, NeuroImage, 2018, 174, pp. 518-538

Figures

Figure 1. Formula used in this abstract.

Figure 2. Overview of the proposed approach. (a) MR pulse sequence diagram for RESOLVE. (b) Phase reconstruction from complex images of b=0. (c)Amplitude reconstruction from the complex images of b=0, 500, and 1200. Diffusion MRI kernel function includes “stick” model and “zeppelin” model. (d) The high-frequency is reconstructed from the phase using crMREPT. (e) Micro-structure parameters are fitted to the diffusion weighted MRI amplitude data. (f) Low-frequency mean conductivity is computed by using decomposition formula.

Figure 3. The fibre phantom (left panel). Amplitude image of b=0, 1000 (middle panel). Phase image of b=0 (right panel).

Figure 4. High-frequency conductivity (upper left panel). Intra-neurite volume fraction (upper middle panel). Intrinsic diffusion coefficient (upper right panel). Extracellular mean diffusion (lower left panel). Low-frequency mean conductivity (lower right panel).

Figure 5. Conductivity ODF ROI1 denotes two cross-fibers. ROI 2 denotes one single-orientation fiber.

Proc. Intl. Soc. Mag. Reson. Med. 32 (2024)
4268
DOI: https://doi.org/10.58530/2024/4268