Next Article in Journal
Endogenous Biomarkers for SLC Transporter-Mediated Drug-Drug Interaction Evaluation
Next Article in Special Issue
Impact of Xylose on Dynamics of Water Diffusion in Mesoporous Zeolites Measured by NMR
Previous Article in Journal
Bioactive Electrospun Fibers of Poly(ε-Caprolactone) Incorporating α-Tocopherol for Food Packaging Applications
Previous Article in Special Issue
The Effect of an Accelerator on Cement Paste Capillary Pores: NMR Relaxometry Investigations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI

1
Department of Biomedical Engineering, Kyung Hee University, Seoul 02447, Korea
2
Division of Applied RI, Korea Institute of Radiological and Medical Science, Seoul 01812, Korea
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(18), 5499; https://doi.org/10.3390/molecules26185499
Submission received: 28 June 2021 / Revised: 30 August 2021 / Accepted: 7 September 2021 / Published: 10 September 2021
(This article belongs to the Special Issue Advances in NMR and MRI of Materials)

Abstract

:
Imaging of the electrical conductivity distribution inside the human body has been investigated for numerous clinical applications. The conductivity tensors of biological tissue have been obtained from water diffusion tensors by applying several models, which may not cover the entire phenomenon. Recently, a new conductivity tensor imaging (CTI) method was developed through a combination of B1 mapping, and multi-b diffusion weighted imaging. In this study, we compared the most recent CTI method with the four existing models of conductivity tensors reconstruction. Two conductivity phantoms were designed to evaluate the accuracy of the models. Applied to five human brains, the conductivity tensors using the four existing models and CTI were imaged and compared with the values from the literature. The conductivity image of the phantoms by the CTI method showed relative errors between 1.10% and 5.26%. The images by the four models using DTI could not measure the effects of different ion concentrations subsequently due to prior information of the mean conductivity values. The conductivity tensor images obtained from five human brains through the CTI method were comparable to previously reported literature values. The images by the four methods using DTI were highly correlated with the diffusion tensor images, showing a coefficient of determination (R 2 ) value of 0.65 to 1.00. However, the images by the CTI method were less correlated with the diffusion tensor images and exhibited an averaged R 2 value of 0.51. The CTI method could handle the effects of different ion concentrations as well as mobilities and extracellular volume fractions by collecting and processing additional B1 map data. It is necessary to select an application-specific model taking into account the pros and cons of each model. Future studies are essential to confirm the usefulness of these conductivity tensor imaging methods in clinical applications, such as tumor characterization, EEG source imaging, and treatment planning for electrical stimulation.

1. Introduction

Electrical conductivity of biological tissues is determined by the cell density, extracellular volume fraction, composition and amount of extracellular matrix materials, and membrane characteristics as well as concentrations and mobility of ions in the extracellular and intracellular fluids [1]. The apparent macroscopic conductivity of such a composite material has been studied since the early 1900s and can be expressed as a weighted sum of conductivity values of its components based on the volume fractions and other factors [2,3]. The extracellular and intracellular fluids are conductors through which conductivity values are determined by concentrations and mobilities of ions and other mobile charge carriers. Cells with thin membranes behave like an insulator and lossy dielectric at low and high frequencies, respectively. Extracellular matrix materials are lossy dielectrics. Therefore, the macroscopic tissue conductivity exhibits frequency dependency [4,5].
When elongated cells are aligned towards a certain direction, movements of ions in the extracellular fluid are consequently hindered. Under a low-frequency electric field, the ions in the extracellular space are forced to move along the longitudinal direction, thereby making their mobilities direction-dependent. Therefore, in the white matter and muscle, the conductivity exhibits anisotropic properties at low frequencies. However, at high frequencies, the insulating cell membranes behave like a capacitor, and the anisotropic properties disappear above 1 MHz, for example [1,6]. In this paper, we approximately express the low-frequency conductivity of biological tissue as a tensor that is a symmetric positive definite 3 × 3 matrix [7,8]. At a high-frequency above 1 MHz, the conductivity is expressed as a scalar quantity [9].
Electrical conductivity is a passive material property whose measurement requires a probing current to generate a signal affected by the conductivity. In impedance imaging area, there are two different approaches in conductivity imaging using MRI. Magnetic resonance electrical impedance tomography (MREIT) reconstructs an image of low-frequency isotropic conductivity ( σ L ) distribution by injecting low-frequency currents into a subject and measuring the induced magnetic flux density distributions using an MRI scanner [10,11,12]. Magnetic resonance electrical properties tomography (MREPT) produces high-frequency isotropic conductivity ( σ H ) and permittivity ( ϵ H ) images by generating a radio-frequency (RF) eddy current that is affected by σ H and ϵ H and measuring an induced RF magnetic field using a B1 mapping method [13,14,15].
In the diffusion tensor imaging (DTI) area, conductivity tensor image reconstructions have been investigated based on a physical relationship that conductivity and water diffusion tensors, denoted as C and D , respectively, have the same eigenvectors [16]. Based on this, Tuch et al. derived the first linear model between C and D , which enabled transformation of a diffusion tensor image into a conductivity tensor image [17,18]. Based on the idea of cross-property relation [17], three more conductivity tensor models and image reconstruction methods using DTI were developed [19,20,21,22]. These conductivity tensor models are based on a linear relation of C = η D and the corresponding image reconstruction methods determine the scale factor η .
Although MREIT and MREPT reconstruct images of σ L and σ H , respectively, through measurement of the effects of both ion concentrations and mobilities, neither of these two methods can produce an image of C . On the other hand, the four methods using DTI produce an image of C , but the effects of ion concentrations are not adequately demonstrated. To overcome these limitations, Sajib et al. developed a novel low-frequency conductivity tensor imaging (CTI) method [23]. In the first step of the CTI method, MREPT is used to reconstruct an image of σ H including the effects of both ion concentrations and mobilities at the Larmor frequency. In the second step, they apply a multi-b diffusion weighted imaging method to obtain an image of D and pixel-by-pixel information about extracellular and intracellular spaces. Combining all of these, the scale factor η between C and D is determined for every pixel using the following: (1) a physical relationship between the water diffusivity and ion mobility and (2) a model-based relation between σ H and σ L [6,24].
Recently, Wu et al. reviewed the DTI-based reconstruction models [25]. They reported that the conversion coefficient of water diffusion tensor to electrical conductivity tensor should require the information on ion concentration and extracellular volume fraction [25]. In addition, knowledge of accurate anisotropic conductivity can achieve reliable volume conduction models of electrical brain stimulation and EEG dipole reconstructions [26,27]. Over the years, considerable variation in the conductivity reconstruction models used for such approaches has been observed [27,28,29]. Consequently, an in-depth investigation of conductivity tensor models was required.
In this paper, we compared the accuracy of the five methods in reconstructed conductivity images of two phantoms using a 9.4 T research MRI scanner. The data acquired from five human subjects using a clinical 3 T MRI scanner were used to reconstruct conductivity tensor images of the brains using the five methods. The reconstructed conductivity tensor images of the human brains were analyzed and compared with each other and also with previous literature values.

2. Five Conductivity Tensor Models

2.1. Linear Eigenvalue Model (LEM)

Adopting the physical analysis that C and D have the same eigenvectors [16], Tuch et al. assumed that the conductivity of the intracellular space is negligible, that is, σ i 0 at low frequencies, and as a result derived the following relation between the eigenvalues of C and D [17]:
c m σ e d e d m d i 3 d e + 1 + d m 2 d i 3 d e 2 2 3 d i
where c m and d m for m = 1 , 2 , 3 are the eigenvalues of C and D , respectively; σ e is the extracellular conductivity; and d i and d e are the intracellular and extracellular water diffusion coefficients, respectively. Assuming that d i 0 , the following linear relation between C and D was derived:
C = σ e d e D = η D .
Without measuring σ e and d e in (2), the empirically-calculated scale factor η = 0.844 S·s/mm 3 was applied to all pixels [17].
Since this empirical scale factor does not consider any intra-subject and inter-subject variabilities, a modified approach was proposed where the volume of the conductivity tensor ellipsoid was matched with the cubed value of the isotropic conductivity using the least square method [18]. For the human brain, the scale factor η was determined as
η = d W M σ W M + d G M σ G M d W M 2 + d G M 2
where d W M and d G M are the measured water diffusion coefficient of the white matter (WM) and gray matter (GM), respectively. In this paper, we used (3) with the literature values of σ W M , G M = 0.14 and 0.27 S/m, respectively, for human brain imaging experiments [30,31,32]. For a conductivity phantom, σ W M and σ G M were replaced by measured conductivity values using an impedance analyzer.

2.2. Force Equilibrium Model (FEM)

The relation between the diffusion coefficient d and viscosity ν is given by the Stokes–Einstein relation as follows:
d = k B T 6 π ν r
where k B is Boltzmann’s constant, T is the absolute temperature, and r is the radius of a spherical particle [19]. From the equilibrium condition between the electric and viscous forces, the conductivity σ can be expressed as
σ = J E = q 2 N 6 π r ν
where J and E are the magnitude of the current density and electric field, respectively, q = 1.6 × 10 19 C and N = 2 × 10 25 m 3 . Applying (4) and (5) to the extracellular space, the following relation can be derived:
σ e = 0.76 q 2 N k B T d e
where k B T = 4.1 × 10 21 J. For an anisotropic case, Sekino et al. assumed that the following relation can be inferred from (6) [19]:
C = 0.76 q 2 N k B T D e = η D e .
To measure D e in (7) for each pixel, Sekino et al. [20] adopted the following bi-exponential model for a diffusion-weighted MRI signal S b with a given b value [33,34]:
S i ( b ) S i ( 0 ) = v f , i e b d f , i + v s , i e b d s , i
where S i ( b ) denote the MR signal with a diffusion gradient ( b 0 ) along the ith direction, S i ( 0 ) is the MRI signal without applying a diffusion gradient ( b = 0 ), v f , i and d f , i are the volume fraction and diffusion coefficient of a fast diffusion component along the ith direction, respectively, and v s , i and d s , i are the volume fraction and diffusion coefficient of a slow diffusion component along the ith direction, respectively. Using a curve fitting method, the fast diffusion tensor D f is extracted from the measured data of S i ( b ) and S i ( 0 ) in three orthogonal directions. In (7), D e is replaced by D f assuming that the fast diffusion corresponds to the extracellular diffusion.

2.3. Volume Constraint Model (VCM)

Miranda et al. [21] suggested a method to determine the scaling factor η at each pixel using the measured water diffusion tensor D and the values of equivalent isotropic conductivity σ i s o for different brain tissues from the given literature [30,31,32,35] as follows:
C = 3 σ i s o t r a c e ( D ) D = η D
where trace ( D ) is the sum of the three eigenvalues of D . In this paper, the conductivity phantom, σ i s o was replaced by a measured conductivity value using an impedance analyzer. For the human brain, we used the literature values of σ i s o = 0.14, 0.27, and 1.79 S/m for the WM, GM, and cerebrospinal fluid (CSF), respectively, [30,31,32,35].

2.4. Volume Fraction Model (VFM)

For the brain tissue, Wang et al. assumed that movements of water molecules and ions are constrained by its multi-compartment environment including axons, glial cells, and CSF [22]. The measured diffusion tensor D was expressed as
D = S d i a g ( d l , d t 1 , d t 2 ) S T
where d l , d t 1 , and d t 2 are the eigenvalues of D along the longitudinal and two transversal directions, respectively. Although the WM may include myelinated axons with various directions, the researchers assumed that three groups of the WM exist with their longitudinal directions in parallel to the longitudinal and two transversal directions of the measured diffusion tensor D . In addition, they assumed that the volume fractions of the three WM groups and the remaining isotropic tissues are α l , α t 1 , α t 2 , and α i s o . Therefore, for this multi-compartment model, the following equations can be obtained:
α l e b d l W + α t 1 e b d t 1 W + α t 2 e b d t 2 W + ( 1 α l α t 1 α t 2 ) e b d i s o = e b d l α l e b d t 1 W + α t 1 e b d l W + α t 2 e b d t 2 W + ( 1 α l α t 1 α t 2 ) e b d i s o = e b d t 1 α l e b d t 2 W + α t 1 e b d t 1 W + α t 2 e b d l W + ( 1 α l α t 1 α t 2 ) e b d i s o = e b d t 2
where d i s o is the diffusion coefficient of the isotropic tissues; d l W , d t 1 W , and d t 2 W are the diffusion coefficients of the WM in the longitudinal and two transversal directions, respectively; and α l + α t 1 + α t 2 + α i s o = 1 .
The eigenvalues of the conductivity tensor C are expressed as a weighted sum of the conductivity values of the four compartments as:
σ l = α l σ l W + α t 1 σ t 1 W + α t 2 σ t 2 W + α i s o σ i s o σ t 1 = α l σ t 1 W + α t 1 σ l W + α t 2 σ t 2 W + α i s o σ i s o σ t 2 = α l σ t 2 W + α t 1 σ t 1 W + α t 2 σ l W + α i s o σ i s o
where σ l W , σ t 1 W , and σ t 2 W are the conductivity values of the WM in the longitudinal and two transversal directions, respectively. The values of σ l W , σ t 1 W , σ t 2 W , and σ i s o are not measured; instead, the literature values are adopted [32]. Note that the VFM method does not use the scale factor η and can be used only for the brain tissue including the WM. For the GM and CSF, isotropic conductivity values from existing literature are used. The constraint of σ t 1 = σ t 2 is applied in the implementation of the VFM method.

2.5. Conductivity Tensor Imaging (CTI) Model

The CTI method derives a low-frequency conductivity tensor C by using a high-frequency isotropic conductivity σ H obtained using the MREPT technique and the information about water diffusion obtained by the multi-b diffusion weighted imaging method. Since the details of its basic theory and algorithm are available in [6,23,24], we simply introduce the following CTI formula in this paper:
Here, the low-frequency conductivity tensor is expressed as
C = α c ¯ e D e
Using the reference value of β = 0.41 [23,24], the apparent extracellular ion concentration ( c ¯ e ) of (13) can be estimated as
c ¯ e = σ H α d e w + ( 1 α ) d i w β
Using c ¯ e from (14) in (13), low-frequency conductivity tensor can be expressed as
C = α c ¯ e D e = α σ H α d e w + ( 1 α ) d i w β D e = η D e
where σ H is the high-frequency conductivity at the Larmor frequency, α is the extracellular volume fraction, c ¯ e is apparent extracellular ion concentration, β is the ion concentration ratio of intracellular and extracellular spaces, d e w and d i w are the extracellular and intracellular water diffusion coefficients, respectively, and D e is the extracellular water diffusion tensor.

3. Imaging Experiments and Data Processing

3.1. Phantom Imaging

To compare the accuracy of the reconstructed conductivity images based on five methods described in Section 2, we used two conductivity phantoms, each with three compartments of known conductivity values. The compartments were filled with electrolytes or giant vesicle suspensions. The giant vesicles were cell-like materials with thin insulating membranes [36].
Phantom #1 comprised of two compartments of electrolytes (EL 1 and EL 2 ) and one compartment of a giant vesicle suspension (GVS 1 ) where giant vesicles were suspended in the electrolyte EL 1 . Phantom #2 comprised of two compartments of different electrolytes (EL 3 and EL 4 ) and one compartment of a different giant vesicle suspension (GVS 2 ). Table 1 shows the concentrations of NaCl and CuSO 4 , the extracellular volume fraction, mobility, and low-frequency conductivity values, which are measured by an impedance analyzer (SI1260A, AMETEK, West Sussex, UK) at 10 Hz. For the electrolyte EL 3 , we increased its viscosity by adding 2 g/L of hyaluronic acid and 10 g/L of polyethylene glycol (PEG, average Mv 8000) solution, thereby decreased the ion mobility and also its conductivity value. In the giant vesicle suspension GVS 2 , we used electrolyte EL 3 .
A 9.4 T research MRI scanner (Agilent Technologies, Santa Clara, CA, USA) equipped with a single-channel birdcage coil (Model: V-HQS-094-00638-029, RAPID Biomedical GmbH, Rimpar, Germany) with a cubic voxel having 0.5 mm edge length was used for phantom imaging. For multi-b diffusion weighted imaging, we used the single-shot spin-echo echo-planar imaging (SS-SE-EPI) pulse sequence. The imaging parameters were as follows: repetition time (TR)/echo time (TE) = 2000/70 ms, number of signal acquisitions = 2, field-of-view (FOV) = 65 × 65 mm 2 , slice thickness = 0.5 mm, flip angle = 90 , and image matrix size = 128 × 128. The number of diffusion-weighting gradient directions was 30 with b-values of 50, 150, 300, 500, 700, 1000, 1400, 1800, 2200, 2600, 3000, 3600, 4000, 4500, and 5000 s/mm 2 . For high-frequency conductivity image reconstructions in the CTI method, B1 phase maps were acquired using the multi-slice multi-echo spin-echo (MS-ME-SE) pulse sequence. The imaging parameters were as follows: TR/TE = 2200/22 ms, number of signal acquisitions = 5, FOV = 65 × 65 mm 2 , slice thickness = 0.5 mm, flip angle = 90 , and image matrix size = 128 × 128. More details of the phantom preparation and data acquisition are described in [6,24].

3.2. In Vivo Human Imaging

Five healthy volunteers were recruited based on the protocol approved by the Institutional Review Board (IRB) at Kyung Hee University (KHSIRB-18-073). Informed consent forms were obtained from the volunteers before conducting the studies. In vivo human brain imaging experiments were performed using a clinical 3 T MRI scanner (Magnetom Trio A Tim, Siemens Medical Solution, Erlangen, Germany) with body coil in transmitting mode and an 8-channel head coil in receiving mode (3D head matrix, A Tim Coil, Siemens Medical Solution, Erlangen, Germany).
For multi-b diffusion weighted imaging, we used the SS-SE-EPI pulse sequence. The number of diffusion-weighting gradient directions was 15 with similar b-values used in the phantom experiments. The imaging parameters were as follows: TR/TE = 2000/70 ms, slice thickness = 4 mm, flip angle = 90 , number of averaging = 2, number of slices = 5, and acquisition matrix = 64 × 64. The matrix size of 64 × 64 was extended to 128 × 128 for subsequent data processing steps. The imaging time was 23 min for the multi-b diffusion data acquisition. For B1 mapping, the MS-ME-SE pulse sequence was used. The imaging parameters were as follows: TR/TE = 1500/15 ms, number of echoes = 6, number of averaging = 5, slice thickness = 4 mm, number of slices = 5, acquisition matrix = 128 × 128, and FOV = 240 × 240 mm 2 using a scan duration of 16 min. For anatomical reference, a conventional T 2 -weighted scan was obtained. The total imaging time for each human subject was about 41 min. More details of the in vivo scans are described in [24].

3.3. Data Processing

The acquired multi-b diffusion data were preprocessed using the MRtrix3 (www.mrtrix.org, accessed on 10 March 2021) [37] and FMRIB software library (FSL, www.fmrib.ox.ac.uk/fsl, accessed on 10 March 2021) [38]. The preprocessing steps included MP-PCA denoising [39], Gibbs-ringing correction [40], and eddy current distortion correction [38]. The averaged images at b = 0 were linearly coregistered to the T 2 -weighted images using the FLIRT method, and the affine transformation matrix was used to nonlinearly coregister the diffusion weighted images at b 0 to the T 2 -weighted images using the FNIRT method [38]. The b-value of 700 s/mm 2 was used for diffusion tensor image reconstructions in the LEM, VCM, and VFM methods. In the FEM and CTI methods, we calculated D f and D s using a bi-exponential model [24,34].
For high-frequency conductivity image reconstructions in CTI method, the acquired phase map data were first corrected for Gibbs-ringing artifacts [40], and then multi-channel multiple echoes were combined to achieve the best signal-to-noise ratio (SNR) [41,42]. The high-frequency conductivity images were reconstructed using the method proposed by Gurler et al. [14] in order to suppress boundary artifacts.
For human brain images, T 2 -weighted images were segmented into the WM, GM, and CSF regions using the MICO [43]. The summary of numbers of pixels in the regions of interest (ROIs) are given in Table 2. In each ROI, we excluded the outermost layer of two-pixel width to reduce partial volume effects. All the data processing steps were implemented using the MATLAB software (Mathworks, Natick, MA, USA). For the image reconstructions using the CTI method, we used the MRCI toolbox (https://iirc.khu.ac.kr/toolbox.html, accessed on 10 March 2021) [44].

4. Results

4.1. Two Phantoms

Figure 1 shows the reconstructed conductivity tensor and diffusion tensor images of phantoms #1 and #2. For both tensors, we plotted their longitudinal and two transversal components. The images of the scale factor η between two tensors were also plotted. Note that the VFM method could not be used for the phantoms since it is only applicable to an anisotropic object. The ROIs were defined as shown in Figure 2a,b corresponding to the three different compartments of phantoms #1 and #2, respectively. The mean and standard deviation (SD) of the reconstructed conductivity values for each ROI were calculated and compared with those that were independently measured using the impedance analyzer at 10 Hz as in Figure 2.
In the case of phantom #1 shown in Figure 1a and Figure 2a, the electrolytes EL 1 and EL 2 had different NaCl concentrations as shown in Table 1. The LEM and FEM methods failed to distinguish the difference in concentration, and the VCM and CTI methods were able to differentiate the difference. In phantom #2 shown in Figure 1b and Figure 2b, the conductivity of EL 3 was reduced compared to that of EL 4 due to the reduced mobility in EL 3 as shown in Table 1. Since the diffusion coefficient was altered by this change in mobility, all four methods could distinguish the difference between EL 3 and EL 4 with the same NaCl concentration in phantom #2. Only the CTI method was able to measure the changes in both concentration and mobility without using prior information of the conductivity values measured by the impedance analyzer.
Table 3 shows the errors in the reconstructed conductivity values with respect to the reference values measured by the impedance analyzer (SI1260A, AMETEK Inc., Berwyn, PA, USA). Since the measured conductivity values using the impedance analyzer were used in place of σ W M , G M , C S F in (3), the values of η in the LEM method were different from each other for two phantoms. However, the LEM method could not account for the effects of the concentration difference in phantom #1 despite being able to detect the mobility difference in phantom #2. In the VCM method, the errors in the four electrolyte regions were zero since the measured conductivity values were used in place of σ i s o in (9). The CTI method could recover the conductivity values with an error of 1.10% to 5.26% in all the regions of the electrolytes and giant vesicle suspensions.

4.2. Five Human Brains

Figure 3 shows the images of the longitudinal ( σ l ) and transversal ( σ t 1 and σ t 2 ) components of the reconstructed conductivity tensor images of the brains of the first and second subject using the five methods. Figure 4 shows the mean and SD values of the reconstructed conductivity tensors in the WM, GM, and CSF regions, which are summarized in Table 2 for all five subjects. For the WM region, the mean and SD values of σ l , σ t 1 , and σ t 2 were plotted for the five different methods. For the GM and CSF regions, we plotted the mean and SD values of σ L = σ l + σ t 1 + σ t 2 3 .
In Figure 4, the LEM and FEM methods underestimated the conductivity values of the CSF region compared with the literature values [35]. For the WM and GM regions, the VFM method overestimated conductivity values compared to the literature values [30,31,32,45]. The conductivity values observed in the VCM method, especially for GM and CSF regions, were comparable to the literature [31,32,45]. Nevertheless, in the LEM and VCM method, we used the literature values of σ i s o for the WM, GM, and CSF regions for all five subjects, and this resulted in a small amount of inter-subject variability. The CTI method produced conductivity values that were comparable to the existing literature values without using prior information of the literature values. In the WM region, the value of σ l obtained by the CTI method was between 0.19 and 0.32 S/m and the values of σ t 1 and σ t 2 were between 0.07 and 0.19 S/m, respectively. In the GM region, the value of σ L was between 0.23 and 0.30 S/m. The value of σ L in the CSF region was between 1.59 and 1.82 S/m with the mean value of 1.72 S/m for all the five subjects, consistent with those found in the existing literature [35,45].
We performed a correlation analyses to visualize the simultaneous influence of water diffusion tensor in reconstructed conductivity tensor. Figure 5 shows the plots of the linear regression analyses between σ l and d l , the longitudinal component of both tensors and averaged transversal components σ t = ( σ t 1 + σ t 2 ) / 2 and d t = ( d t 1 + d t 2 ) / 2 for the WM, GM, and CSF regions in all five human subjects (Appendix A). For the LEM and FEM methods, the coefficient of determination ( R 2 ) was 0.99 and 1.00, respectively. Therefore, these two methods, may not provide additional information that is not available in the water diffusion tensor. For the VCM and VFM methods, R 2 ranges from 0.65 to 0.79, respectively, whereas the CTI method showed lesser correlation between σ l , t and d l , t with R 2 ranges from 0.46 and 0.56. In the CTI method, the magnitudes of C and D provided moderately dependent information although they have the same directional property, i.e., the same eigenvectors. This is also because the conductivity tensor in the CTI method has a predominant effect of apparent ion concentrations.
To analyze the directional property, we also computed the anisotropy ratio (AR) in the WM region expressed as A R C = 2 σ l σ t 1 + σ t 2 and A R D = 2 d l d t 1 + d t 2 for the conductivity tensor and diffusion tensor, respectively [46]. The mean and SD value of A R D were 2.52 ± 0.76, whereas the mean and SD values of A R C were 2.53 ± 0.77, 2.52 ± 0.78, 2.50 ± 0.78, 1.73 ± 0.47, and 2.52 ± 0.84 for the LEM, FEM, VCM, VFM, and CTI methods, respectively, for the five human subjects. The lower value of A R C in the VFM model can be attributed to the assumption of σ t 1 = σ t 2 . For the other four methods, the values of A R C and A R D were similar since the structural property primarily determined the anisotropy.
Using the CTI method as a reference, we plotted the images of the relative difference ( r d D I R ), which are defined as
r d D I R = σ D I R , C T I σ D I R , M T H σ D I R , C T I × 100 ( % )
where DIR is l , t 1 , or t 2 and σ D I R , C T I and σ D I R , M T H are the reconstructed σ l , σ t 1 , or σ t 2 using the CTI method and one of the LEM, FEM, VCM, and VFM methods, respectively (Appendix A). Figure 6 shows the images of r d D I R for all five subjects pooled together. Table 4 summarizes the absolute mean values of | r d D I R | for the WM, GM, and CSF regions for each subject. The bar plots at top of Figure 6 show absolute mean values of | r d D I R | for all five subjects for WM, GM, and CSF regions, and images at the bottom of Figure 6 show the relative difference with sign of deviation to CTI method. The mean and SD of the relative difference for the entire brain slice of the five subjects were 69.63 ± 31.11, 104.94 ± 35.41, 53.35 ± 25.15, and 68.83 ± 14.38% for the LEM, FEM, VCM, and VFM methods, respectively.

5. Discussion

In this paper, we compared five conductivity tensor models and image reconstruction methods using two phantoms and five human subjects. The data from the phantom experiments were used to evaluate the accuracy of the reconstructed conductivity images against the measured conductivity values obtained using the impedance analyzer. In vivo human experiment data were used to compare the values of the reconstructed conductivity tensors with existing literature values and also with each other using the CTI method as a reference.
In the phantom experiments, the LEM and FEM methods could not distinguish the electrolytes with different NaCl concentrations. The errors in the reconstructed conductivity images using the LEM and FEM were relatively large. The errors of the VCM method were the smallest since we used the conductivity values measured using the impedance analyzer in place of σ i s o in (9). In most clinical applications where the value of σ i s o is unknown, the VCM method may fail to recover the inhomogeneous conductivity values. In addition, living tissues are heterogeneous, and their conductivity varies with the pathophysiology condition. Thus, a single global value of conductivity to the entire tissue type may not sufficient. The VFM method could not be applied to the phantoms since it was explicitly designed for anisotropic brain tissues. Without relying on prior information of mean isotropic conductivity values, the CTI method recovered conductivity values with an error ranging from 1.10% to 5.26% for all six different compartments in the two phantoms having different ion concentrations and mobilities. Furthermore, the CTI method could properly handle the effects of different cell densities in two giant vesicle suspensions with different extracellular volume fractions.
Using the data from in vivo human brain imaging experiments, we produced conductivity tensor images of the brains as shown in Figure 3 using the five methods. Compared to the conductivity values of the WM, GM, and CSF from existing literature, the values using the LEM and FEM methods had large amounts of bias in the WM and GM regions. Since the global scale factor was used for all pixels in the FEM method, the method’s ability to handle intra-subject and inter-subject variability appeared to be primarily limited. Although the scale factor η in the LEM method contains inter-subject variability, but still lacks intra-subject tissue heterogeneity. Furthermore, the global scale factor does not account for position dependence, and its deviation in the tissues only explains the water diffusion alterations. The conductivity values in CSF from both LEM and FEM methods were significantly underestimated compared to the literature values. The calculated conductivity tensor in both these methods was simply linearly scaled water diffusion tensor, and hence may explain their considerable low conductivity measurements [17,19]. The conductivity values measured in our study were matched with previous studies by Rullman et al. and Sekino et al. [18,20].
In the volume constrained method (VCM), conductivity values of GM and CSF were matched with the literature and also with the CTI method in our study. In contrast, conductivity values of WM were underestimated at longitudinal directions. This might be due to the use of similar literature values of σ i s o in all three fiber directions. This method still lacks the intra-voxel heterogenic property of the tissues, and segmentation inaccuracies may end up assigning incorrect σ i s o values. An introduction of the isotropic extracellular electrical conductivity from MREIT could be helpful to overcome this limitation [12]. In the volume fraction method (VFM), conductivity values in WM and GM were highly overestimated. This stems from the fact that the eigenvalues of conductivity tensor were computed as a weighted sum of conductivity values associated with four compartments of individual fiber bundle direction. VFM method is computationally complicated, and the ill-posed problem may lead to failed fit in noisy voxels [22]. In addition, the assumption of constant σ L and σ t 1 = σ t 2 , which were then set to be literature conductivity value, somehow neglected the intra-voxel heterogeneity, and reduced the anisotropic ratio. Instead of such limitations, VFM method effectively handled the partial volume effect.
Although the scale factor η in the LEM, FEM, and VCM methods was computed for each pixel, their ability to handle the intra-subject and inter-subject variability was also limited as a result of fixed literature conductivity values of the WM, GM, and CSF regions for all subjects. The VFM method manages cross-subject variability but failed with intra-subject flexibility. In contrast to the four methods using DTI only, the CTI method produced conductivity tensor images with values that were comparable to existing literature values. Without using prior information about tissue conductivity values, the CTI method appeared to properly handle the intra-subject and inter-subject variability. However, the intra-subject and inter-subject variabilities, may have been affected in all methods by measurement noise, and systematic artifacts.
The recovered conductivity tensors using the LEM and FEM methods in Figure 5a to Figure 5b shows similar results as the water diffusion tensor in our study. In spite of the calculating the pixel-dependent scale factor in the VCM method, a higher correlation with water diffusion was observed (Figure 5c). Considering that the local concentration of a charge carrier can affect the linearity between C and D , the CTI method appeared to properly reflect the effects of position-dependent concentration differences by incorporating the actually measured high-frequency conductivity σ H into the conductivity tensor image reconstruction. Note that these linear regression analyses should not be considered as a general validity of C = η D . In CTI method, the mutual restriction of both the ionic and the water mobility by the geometry of the brain medium builds the basis for the relationship of C and D . The voxels in conductivity tensor map are expressed as a sum of products of the carrier concentration, extracellular volume fraction, and the mobility tensor. Since all of the five methods utilized the mobility information embedded in the water diffusion tensor D , the anisotropy ratio of C was primarily determined by the anisotropy ratio of D .
Data acquisition and processing steps are more involved in the CTI method than in those other four methods using DTI. Given that there are various sources of error in the data acquisition and processing stages, subsequent studies are necessary to rigorously validate the performance of the CTI method in terms of errors due to measurement noise and artifacts, partial volume effects, and coregistration of different images. In particular, high-frequency conductivity images using MREPT suffer from errors at boundaries of two regions with different conductivity values and also assumes the piecewise constant conductivity [47], which does not hold in practice. Although the MREPT method we adopted in this paper could reduce such boundary artifacts [14] at the expense of smooth images, future improvements in high-frequency conductivity image reconstruction algorithms are needed to enhance the accuracy of the CTI method in human subjects. The parameter β in (15) was estimated from the literature values and assumed to be constant for all pixels. Although the sensitivity of a reconstructed conductivity tensor image to β appeared to be small [24], clinical studies are needed to find out any practical limitations imposed by this assumption. The CTI method can provide clinically useful information about pathological and physiological changes in cells and cellular structures associated with disease progression. It is one of most promising clinical applications is tumor imaging in terms for early detection, better characterization, and monitoring after a treatment [48,49]. The recovered low-frequency conductivity tensor could be also used to predict internal current pathways and electric field distributions subject to externally injected or induced therapeutic currents in electrical stimulation [28,29].
This study has several limitations that should be considered in the future works. The distribution of conductivity within the giant vesicle phantom does not reflect an adequate anisotropic environment. Future studies should include phantoms mimicking microstructural properties similar to the brain [50]. Diffusion MRI suffers from various systematic errors, such as gradient inhomogeneity [51]. Although we corrected the data using commonly used correction methods, a study of systematic error propagation or noise analysis is needed to enhance the accuracy of the conductivity images. In human brain, conventional DTI measures water diffusion assuming that displacement distribution of water molecules in a given time is a Gaussian function. However, this assumption may not be valid in complex biological tissue where water molecules often show non-Gaussian diffusion [52]. DTI also has limitations in characterizing the diffusion process in areas of low anisotropy and complex fiber structure in a voxel. Future studies with diffusion kurtosis imaging (DKI) or high angular methods (HARDI) can provide better characterization of human brain architecture [52,53]. Cell membranes are assumed to resist low-frequency currents in the CTI approach. As a result, the CTI approach may underestimate the low-frequency conductivity value of tissue including cells with leaky membranes. It would be worthwhile to investigate a more sophisticated CTI model involving such cells.

6. Conclusions

In this study, we provided an overview of the current state of the art in MRI-based conductivity tensor reconstruction. The accuracy of five conductivity tensor model was investigated using two phantoms with four electrolytes and two giant vesicles suspensions with known internal conductivity values. The findings show that methods with pixel-dependent scale factors work better than methods with global scale factors. The experimental results showed that the accuracy of the VCM and CTI methods was superior to that of the LEM and FEM methods. Contrary to the four methods using only DTI, the CTI method did not use prior information on mean isotropic conductivity values, and produced conductivity images with errors ranging from 1.10% to 5.26%. From in vivo human brain imaging experiments, the reconstructed conductivity values of the white and gray matter using the LEM, VCM, and CTI methods were comparable with the values available in the literature. Except for LEM and FEM, all methods yielded conductivity values of the CSF similar to those of literature. Methods using only the water diffusion tensor and prior knowledge of the isotropic mean conductivity values varied depending on the parameter value used. The CTI method was able to properly handle the effects of different ion concentrations as well as mobilities and extracellular volume fractions in our study. Although the data processing of the CTI method is more involved, it appeared to be the most accurate among the five methods in this study. Future research is required to confirm the clinical utility of these low-frequency conductivity tensor image reconstruction approaches in diagnostic imaging and bioelectromagnetic modeling.

Author Contributions

N.K. designed the MR experiment, processed and analyzed the data. and drafted the manuscript. B.-K.C. and I.-O.K. prepared the giant vesicle, performed the MR experiment, and analyzed the data. J.-A.P. and H.-J.K. analyzed the results and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF), the Ministry of Health and Welfare of Korea, and Korea Institute of Radiological and Medical Sciences (KIRAMS) grants funded by the Korea government (No. 2019R1A2C2088573, 2020R1A2C200790611, 2021R1A2C2004299, HI18C2435, and 50461-2021).

Institutional Review Board Statement

The protocols are approved by the Institutional Review Board at Kyung Hee University(KHSIRB-18-073).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Some of the data used in this study are available on https://iirc.khu.ac.kr/toolbox.html, (accessed on 10 March 2021) or available from corresponding author upon reasonable request.

Acknowledgments

The authors thank Prof EunAh Lee for the helpful advice on giant vesicle suspension.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

LEMLinear Eigenvalue Model
FEMForce Equilibrium Model
VCMVolume Constraint Model
VFMVolume Fraction Model
CTIConductivity Tensor Imaging
WMWhite Matter
GMGray Matter
CSFCerebrospinal Fluid
MREPTMagnetic Resonance Electrical Properties Tomography

Appendix A

In this section, the main symbols used in this paper can be found with description and units.
σ L Isotropic low-frequency conductivity (S/m)
σ l Longitudinal component of conductivity tensor (S/m)
σ t Transversal component of conductivity tensor (S/m)
d l Longitudinal component of water diffusion tensor ( μ m 2/ ms )
d t Transversal component of water diffusion tensor ( μ m 2/ ms )
σ i s o Isotropic low-frequency conductivity value from the literature (S/m)
d e Isotropic extracellular diffusion coefficient ( μ m 2/ ms )
D e Extracellular water diffusion tensor ( μ m 2/ ms )

References

  1. Grimnes, S.; Martinsen, O.G. Bioimpedance and Bioelectricity Basics; Academic Press: London, UK, 2015. [Google Scholar]
  2. Kerner, E.H. The electrical conductivity of composite materials. Proc. Phys. Soc. B 1956, 69, 802–807. [Google Scholar] [CrossRef]
  3. Lux, F. Models proposed to explain the electrical conductivity of mixtures made of conductive and insulating materials. J. Mater. Sci. 1993, 28, 285–301. [Google Scholar] [CrossRef]
  4. Schwan, H.P. Electrical properties of tissue and cell suspensions. Adv. Biol. Med. Phys. 1957, 5, 147–209. [Google Scholar] [PubMed]
  5. Gabreil, C.; Peyman, A.; Grant, E.H. Electrical conductivity of tissue at frequency below 1 MHz. Phys. Med. Biol. 2009, 54, 4863–4878. [Google Scholar] [CrossRef] [PubMed]
  6. Choi, B.K.; Katoch, N.; Park, J.E.; Ko, I.O.; Kim, H.J.; Kwon, O.I.; Woo, E.J. Validation of conductivity tensor imaging using giant vesicle suspensions with different ion mobilities. Biomed. Eng. OnLine 2020, 19, 1–17. [Google Scholar] [CrossRef] [PubMed]
  7. Basser, P.J.; Mattiello, J.; Le bihan, D. MR diffusion tensor spectroscopy and imaging. Biophys. J. 1994, 66, 259–267. [Google Scholar] [CrossRef] [Green Version]
  8. Kwon, O.I.; Jeong, W.C.; Sajib, S.Z.K.; Kim, H.J.; Woo, E.J. Anisotropic conductivity tensor imaging in MREIT using directional diffusion rate of water molecules. Phys. Med. Biol. 2014, 59, 2955–2974. [Google Scholar] [CrossRef]
  9. Liu, J.; Yang, Y.; Katscher, U.; He, B. Electrical properties tomography based on B1 maps in MRI: Principles, applications, and challenges. IEEE Trans. Biomed. Eng. 2017, 64, 2515–2530. [Google Scholar] [CrossRef]
  10. Seo, J.K.; Woo, E.J. Magnetic resonance electrical impedance tomography (MREIT). SIAM Rev. 2011, 53, 40–68. [Google Scholar] [CrossRef]
  11. Seo, J.K.; Kim, D.H.; Lee, J.; Kwon, O.I.; Sajib, S.Z.K.; Woo, E.J. Electrical tissue property imaging using MRI at dc and Larmor frequency. Inv. Prob. 2012, 28, 084002. [Google Scholar] [CrossRef]
  12. Seo, J.K.; Woo, E.J. Electrical tissue property imaging at low frequency using MREIT. IEEE Trans. Biomed. Eng. 2014, 61, 1390–1399. [Google Scholar] [PubMed]
  13. Katscher, U.; Voigt, T.; Findeklee, C.; Vernickel, P.; Nehrke, K.; Dossel, O. Determination of electrical conductivity and local SAR via B1 mapping. IEEE Trans. Med. Imag. 2009, 28, 1365–1374. [Google Scholar] [CrossRef] [PubMed]
  14. Gurler, N.; Ider, Y.Z. Gradient-based electrical conductivity imaging using MR phase. Mag. Reson. Med. 2016, 77, 137–150. [Google Scholar] [CrossRef] [PubMed]
  15. Leijsen, R.; Brink, W.; van den Berg, C.; Webb, A.; Remis, R. Electrical properties tomography: A methodological review. Diagnostics 2021, 11, 176. [Google Scholar] [CrossRef] [PubMed]
  16. Sen, A.K.; Torquato, S. Effective conductivity of anisotropic two-phase composite media. Phys. Rev. B 1989, 39, 4504–4515. [Google Scholar] [CrossRef] [PubMed]
  17. Tuch, D.S.; Wedeen, V.J.; Dale, A.M.; George, J.S.; Belliveau, J.W. Conductivity tensor mapping of the human brain using diffusion tensor MRI. Proc. Nat. Acad. Sci. USA 2001, 98, 11697–11701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Rullmann, M.; Anwander, A.; Dannhauer, M.; Warfield, S.K.; Duffy, F.H.; Wolters, C.H. EEG source analysis of epileptiform activity using a 1 mm anisotropic hexahedra finite element head model. NeuroImage 2009, 44, 399–410. [Google Scholar] [CrossRef] [Green Version]
  19. Sekino, M.; Yamaguchi, K.; Iriguchi, N.; Ueno, S. Conductivity tensor imaging of the brain using diffusion-weighted magnetic resonance imaging. J. App. Phys. 2003, 93, 6430–6732. [Google Scholar] [CrossRef]
  20. Sekino, M.; Inoue, Y.; Ueno, S. Magnetic resonance imaging of electrical conductivity in the human brain. IEEE Trans. Mag. 2005, 41, 4203–4205. [Google Scholar] [CrossRef]
  21. Miranda, P.C.; Pajevic, S.; Hallett, M.; Basser, P. The distribution of currents induced in the brain by magnetic stimulation: A finite element analysis incorporating DT-MRI-derived conductivity data. Proc. Int. Soc. Mag. Res. Med. 2001, 9, 1540. [Google Scholar]
  22. Wang, K.; Zhu, S.; Mueller, B.A.; Lim, K.O.; He, B. A new method to derive white matter conductivity from diffusion tensor MRI. IEEE Trans. Biomed. Eng. 2008, 55, 2481–2486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sajib, S.Z.K.; Kwon, O.I.; Kim, H.J.; Woo, E.J. Electrodeless conductivity tensor imaging (CTI) using MRI: Basic theory and animal experiments. Biomed. Eng. Lett. 2018, 8, 273–282. [Google Scholar] [CrossRef] [PubMed]
  24. Katoch, N.; Choi, B.K.; Sajib, S.Z.K.; Lee, E.; Kim, H.J.; Kwon, O.I.; Woo, E.J. Conductivity tensor imaging of in vivo human brain and experimental validation using giant vesicle suspension. IEEE Trans. Biomed. Eng. 2018, 38, 1569–1577. [Google Scholar] [CrossRef] [PubMed]
  25. Wu, Z.; Liu, Y.; Hong, M.; Yu, X. A review of anisotropic conductivity models of brain white matter based on diffusion tensor imaging. Med. Biol. Eng. Comput. 2018, 56, 1325–1332. [Google Scholar] [CrossRef] [PubMed]
  26. Wolters, C.H.; Anwander, A.; Tricoche, X.; Weinstein, D.; Koch, M.A.; Macleod, R.S. Influence of tissue conductivity anisotropy on eeg/meg field and return current computation in a realistic head model: A simulation and visualization study using high-resolution finite element modeling. NeuroImage 2006, 30, 813–826. [Google Scholar] [CrossRef] [Green Version]
  27. Vorwerk, J.; Aydin, Ü; Wolters, C.H.; Butson, C.R. Influence of head tissue conductivity uncertainties on EEG dipole reconstruction. Front. Neurosci. 2009, 13, 531. [Google Scholar]
  28. Shahid, S.; Wen, P.; Ahfock, T. Numerical investigation of white matter anisotropic conductivity in defining current distribution under tDCS. Comput. Meth. Prog. Biomed. 2013, 9, 48–64. [Google Scholar] [CrossRef]
  29. Lee, W.H.; Liu, Z.; Mueller, B.A.; Lim, K.O.; He, B. Influence of white matter anisotropic conductivity on EEG source localization: Comparison to fMRI in human primary visual cortex. Clin. Neurophysiol. 2009, 120, 2071–2081. [Google Scholar] [CrossRef] [Green Version]
  30. Nicholson, P.W. Specific impedance of cerebral white matter. Exp. Neurol. 1965, 13, 386–401. [Google Scholar] [CrossRef]
  31. Geddes, L.A.; Baker, L.E. The specific resistance of biological material: A compendium of data for the biomedical engineer and physiologist. Phys. Med. Biol. 1967, 44, 193–271. [Google Scholar] [CrossRef]
  32. Gabriel, C.; Gabriel, S.; Corthout, E. The dielectric properties of biological tissues: I. Literature survey. Phys. Med. Biol. 1996, 41, 2231–2249. [Google Scholar] [CrossRef] [Green Version]
  33. Niendorf, T.; Dijkhuizen, R.M.; Norris, D.G.; Campagne, M.V.K.; Nicolay, K. Biexponential diffusion attenuation in various states of brain tissue: Implications for diffusion-weighted imaging. Mag. Reson. Med. 1996, 36, 847–857. [Google Scholar] [CrossRef]
  34. Clark, C.A.; Hedehus, M.; Moseley, M.E. In vivo mapping of the fast and slow diffusion tensors in human brain. Mag. Reson. Med. 2002, 47, 623–628. [Google Scholar] [CrossRef]
  35. Baumann, S.B.; Wozny, D.R.; Kelly, S.K.; Meno, F.M. The electrical conductivity of human cerebrospinal fluid at body temperature. IEEE Trans. Biomed. Eng. 1997, 44, 220–223. [Google Scholar] [CrossRef]
  36. Moscho, A.; Orwar, O.; Chiu, D.T.; Modi, B.P.; Zare, R.N. Rapid preparation of the giant unilamellar vesicles. Proc. Nat. Acad. Sci. USA 1996, 93, 11443–11447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Tournier, J.D.; Smith, R.; Raffelt, D.; Tabbara, R.; Dhollander, T.; Pietsch, M.; Christiaens, D.; Jeurissen, B.; Yeh, C.H.; Connelly, A. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. NeuroImage 2019, 202, 116137. [Google Scholar] [CrossRef]
  38. Smith, S.M.; Jenkinson, M.; Woolrich, M.W.; Beckmann, C.F.; Behrens, T.E.J.; Johansen-Berg, H.; Bannister, P.R.; De Luca, M.; Drobnjak, I.; Flitney, D.E.; et al. Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage 2004, 23, 208–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Veraart, J.; Novikov, D.S.; Christiaens, D.; Ades-Aron, B.; Sijbers, J.; Fieremans, E. Denoising of diffusion MRI using random matrix theory. NeuroImage 2016, 142, 394–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Kellner, E.; Dhital, B.; Kiselev, V.G.; Reisert, M. Gibbs-ringing artifact removal based on local subvoxel-shifts. Magn. Reson. Med. 2016, 76, 1574–1581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Walsh, D.O.; Gmitro, A.F.; Marcellin, M.W. Adaptive reconstruction of phased array MR imagery. Magn. Reson. Med. 2000, 43, 682–690. [Google Scholar] [CrossRef]
  42. Kwon, O.I.; Jeong, W.C.; Sajib, S.Z.K.; Kim, H.J.; Woo, E.J.; Oh, T.I. Reconstruction of dual-frequency conductivity by optimization of phase map in MREIT and MREPT. Biomed. Eng. Online 2014, 13, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, C.; Liu, Z.; Gore, J.C.; Davatzikos, C. Multiplicative intrinsic component optimization (MICO) for MRI bias field estimation and tissue segmentation. Mag. Res. Imag. 2014, 32, 913–923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sajib, S.Z.K.; Katoch, N.; Kim, H.J.; Kwon, O.I.; Woo, E.J. Software toolbox for low-frequency conductivity and current density imaging using MRI. IEEE Trans. Biomed. Eng. 2017, 64, 2505–2514. [Google Scholar] [PubMed]
  45. Chauhan, M.; Indahlastari, A.; Kasinadhuni, A.K.; Schar, M.; Mareci, T.H.; Sadleir, R.J. Low-frequency conductivity tensor imaging of the human head in vivo using DT-MREIT: First study. IEEE Trans. Med. Imag. 2018, 37, 966–976. [Google Scholar] [CrossRef]
  46. Pierpaoli, C.; Basser, P.J. Toward a quantitative assessment of diffusion anisotropy. Mag. Res. Imag. 1996, 36, 893–906. [Google Scholar] [CrossRef]
  47. Seo, J.K.; Kim, M.O.; Lee, J.S.; Choi, N.; Woo, E.J.; Kim, H.J.; Kwon, O.I.; Kim, D.H. Error analysis of non constant admittivity for MR-based electric property imaging. IEEE Trans. Med. Imag. 2012, 31, 430–437. [Google Scholar]
  48. Tha, K.K.; Katscher, U.; Yamaguchi, S.; Stehning, C.; Terasaka, S.; Fujima, N.; Kudo, K.; Kazumata, K.; Yamamoto, T.; Van Cauteren, M.; et al. Noninvasive electrical conductivity measurement by MRI: A test of its validity and the electrical conductivity characteristics of glioma. Eur. Radiol. 2018, 28, 348–355. [Google Scholar] [CrossRef]
  49. Lesbats, C.; Katoch, N.; Minhas, A.S.; Taylor, A.; Kim, H.J.; Woo, E.J.; Poptani, H. High-frequency electrical properties tomography at 9.4T as a novel contrast mechanism for brain tumors. Mag. Reson. Med. 2021, 86, 382–892. [Google Scholar] [CrossRef] [PubMed]
  50. Fieremans, E.; Lee, H.H. Physical and numerical phantoms for the validation of brain microstructural MRI: A cookbook. Neuroimage 2018, 182, 39–61. [Google Scholar] [CrossRef]
  51. Dowell, N.; Tofts, P. Quality Assurance for Diffusion MRI; Oxford University Press Inc.: Oxford, UK, 2011. [Google Scholar]
  52. Veraart, J.; Poot, D.H.; Van Hecke, W.; Blockx, I.; Van der Linden, A.; Verhoye, M.; Sijbers, J. More accurate estimation of diffusion tensor parameters using diffusion kurtosis imaging. Mag. Reson. Med. 2011, 65, 138–145. [Google Scholar] [CrossRef]
  53. Tuch, D.S.; Reese, T.G.; Wiegell, M.R.; Makris, N.; Belliveau, J.W.; Wedeen, V.J. High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Mag. Reson. Med. 2002, 48, 577–582. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Reconstructed conductivity tensor images of phantoms #1 (a) and #2 (b). Target images are the true conductivity tensor images generated using the conductivity values measured by the impedance analyzer. (c,d) are the diffusion tensor images of phantoms #1 and #2, respectively. The VFM method was not applicable to the phantoms since it was specifically designed for the brain tissues.
Figure 1. Reconstructed conductivity tensor images of phantoms #1 (a) and #2 (b). Target images are the true conductivity tensor images generated using the conductivity values measured by the impedance analyzer. (c,d) are the diffusion tensor images of phantoms #1 and #2, respectively. The VFM method was not applicable to the phantoms since it was specifically designed for the brain tissues.
Molecules 26 05499 g001
Figure 2. Values of σ L = σ l + σ t 1 + σ t 2 3 in the reconstructed conductivity tensor images of phantoms #1 (a) and #2 (b). The values were compared with the conductivity values measured at 10 Hz using the impedance analyzer. The square symbol indicates the mean value and the bar represents the SD.
Figure 2. Values of σ L = σ l + σ t 1 + σ t 2 3 in the reconstructed conductivity tensor images of phantoms #1 (a) and #2 (b). The values were compared with the conductivity values measured at 10 Hz using the impedance analyzer. The square symbol indicates the mean value and the bar represents the SD.
Molecules 26 05499 g002
Figure 3. (a,b) are the longitudinal ( σ l ) and transversal ( σ t 1 and σ t 2 ) components of the reconstructed conductivity tensor images of the human brains from subjects #1 and #2, respectively, using the five methods of LEM, FEM, VCM, VFM, and CTI. (c,d) show the images of the water diffusion tensors.
Figure 3. (a,b) are the longitudinal ( σ l ) and transversal ( σ t 1 and σ t 2 ) components of the reconstructed conductivity tensor images of the human brains from subjects #1 and #2, respectively, using the five methods of LEM, FEM, VCM, VFM, and CTI. (c,d) show the images of the water diffusion tensors.
Molecules 26 05499 g003
Figure 4. Mean and SD values of the reconstructed conductivity tensors using the five methods for all five human subjects are compared with the existing literature values. (a) are conductivity values from WM region, (b,c) are from GM and CSF, respectively. For the WM region, the mean and SD values of σ l , σ t 1 , and σ t 2 are plotted. For the GM and CSF regions, the mean and SD values of σ L = σ l + σ t 1 + σ t 2 3 are plotted. The vertical bar represents the SD.
Figure 4. Mean and SD values of the reconstructed conductivity tensors using the five methods for all five human subjects are compared with the existing literature values. (a) are conductivity values from WM region, (b,c) are from GM and CSF, respectively. For the WM region, the mean and SD values of σ l , σ t 1 , and σ t 2 are plotted. For the GM and CSF regions, the mean and SD values of σ L = σ l + σ t 1 + σ t 2 3 are plotted. The vertical bar represents the SD.
Molecules 26 05499 g004
Figure 5. Results of the linear regression analyses between s l and d l , the longitudinal component of both tensors and averaged transversal components σ t = σ t 1 + σ t 2 2 and d t = d t 1 + d t 2 2 in the reconstructed conductivity and water diffusion tensor images of all five human subjects pooled together.
Figure 5. Results of the linear regression analyses between s l and d l , the longitudinal component of both tensors and averaged transversal components σ t = σ t 1 + σ t 2 2 and d t = d t 1 + d t 2 2 in the reconstructed conductivity and water diffusion tensor images of all five human subjects pooled together.
Molecules 26 05499 g005
Figure 6. Images of the relative differences in the longitudinal ( e l ) and transversal ( e t 1 and e t 2 ) components between one of the LEM, FEM, VCM, and VFM methods and the CTI method. The graphs at the top are the absolute mean and SD values of the relative differences ( | r d D I R | ) in the WM, GM, and CSF regions. The bar represents the SD.
Figure 6. Images of the relative differences in the longitudinal ( e l ) and transversal ( e t 1 and e t 2 ) components between one of the LEM, FEM, VCM, and VFM methods and the CTI method. The graphs at the top are the absolute mean and SD values of the relative differences ( | r d D I R | ) in the WM, GM, and CSF regions. The bar represents the SD.
Molecules 26 05499 g006
Table 1. Compositions of two conductivity phantoms. The electrolyte and giant vesicle suspension are denoted as EL and GVS, respectively. EL 1 , EL 2 , and GVS 1 were used in phantom #1. EL 3 , EL 4 , and GVS 2 were used in phantom #2.
Table 1. Compositions of two conductivity phantoms. The electrolyte and giant vesicle suspension are denoted as EL and GVS, respectively. EL 1 , EL 2 , and GVS 1 were used in phantom #1. EL 3 , EL 4 , and GVS 2 were used in phantom #2.
CompartmentEL1EL2GVS1EL3EL4GVS2
NaCl (g/L)7.53.57.5333
CuSO 4 (g/L)010000
Extracellular1001001010010050
volume fraction
(%)
Mobilityhighhighhighlowhighlow
σ at 10 Hz (S/m)1.560.830.290.550.700.45
Table 2. Numbers of pixels in the WM, GM, and CSF ROI of the five human brains. The reconstructed conductivity values in these ROIs were compared based on the five image reconstruction methods and also with existing literature values.
Table 2. Numbers of pixels in the WM, GM, and CSF ROI of the five human brains. The reconstructed conductivity values in these ROIs were compared based on the five image reconstruction methods and also with existing literature values.
ROISubject
#1#2#3#4#5
WM10931026125710701047
GM10579187859271023
CSF209135180237234
Table 3. Errors of the reconstructed conductivity values σ L = σ l + σ t 1 + σ t 2 3 with respect to the reference values measured by the impedance analyzer at 10 Hz. Note that the reference values were themselves used in the VCM method for the compartments of EL 1 , EL 2 , EL 3 , and EL 4 .
Table 3. Errors of the reconstructed conductivity values σ L = σ l + σ t 1 + σ t 2 3 with respect to the reference values measured by the impedance analyzer at 10 Hz. Note that the reference values were themselves used in the VCM method for the compartments of EL 1 , EL 2 , EL 3 , and EL 4 .
ROILEM (%)FEM (%)VCM (%)CTI (%)
EL 1 29.6086.2401.10
EL 2 28.1475.5704.42
GVS 1 131.1754.833.451.74
EL 3 32.9773.2703.39
EL 4 67.8266.2705.26
GVS 2 28.1674.222.022.13
Table 4. Mean values of the absolute relative differences ( | r d D I R | ) in the WM, GM, and CSF regions for all five human subjects. The relative differences of the LEM, FEM, VCM, and VFM methods were computed with respect to the CTI method.
Table 4. Mean values of the absolute relative differences ( | r d D I R | ) in the WM, GM, and CSF regions for all five human subjects. The relative differences of the LEM, FEM, VCM, and VFM methods were computed with respect to the CTI method.
SubjectLEM and CTI (%)FEM and CTI (%)VCM and CTI (%)VFM and CTI (%)
WMGMCSFWMGMCSFWMGMCSFWMGMCSF
#144.2550.99130.2463.1087.90171.3950.2768.6030.5891.8689.7331.86
#249.9048.29112.2158.9295.20151.7453.9373.2539.0391.6089.0840.66
#340.3348.87121.0763.3599.29165.4651.6366.4735.2392.7190.6137.09
#443.3044.00106.8349.3663.41130.4649.9165.1735.3092.2690.5650.81
#546.9847.66158.9056.63101.28207.5962.1973.1040.6689.0480.1539.43
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Katoch, N.; Choi, B.-K.; Park, J.-A.; Ko, I.-O.; Kim, H.-J. Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI. Molecules 2021, 26, 5499. https://doi.org/10.3390/molecules26185499

AMA Style

Katoch N, Choi B-K, Park J-A, Ko I-O, Kim H-J. Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI. Molecules. 2021; 26(18):5499. https://doi.org/10.3390/molecules26185499

Chicago/Turabian Style

Katoch, Nitish, Bup-Kyung Choi, Ji-Ae Park, In-Ok Ko, and Hyung-Joong Kim. 2021. "Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI" Molecules 26, no. 18: 5499. https://doi.org/10.3390/molecules26185499

APA Style

Katoch, N., Choi, B. -K., Park, J. -A., Ko, I. -O., & Kim, H. -J. (2021). Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI. Molecules, 26(18), 5499. https://doi.org/10.3390/molecules26185499

Article Metrics

Back to TopTop