Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 421619, 12 pages

http://dx.doi.org/10.1155/2013/421619

Research Article

Noninvasive Measurement of Conductivity Anisotropy at Larmor Frequency Using MRI

Joonsung Lee,¹ Yizhuang Song,² Narae Choi,¹ Sungmin Cho,¹ Jin Keun Seo,² and Dong-Hyun Kim¹

¹Department of Electrical and Electronic Engineering, Yonsei University, Seoul 120-749, Republic of Korea

²Department of Computational Science & Engineering, Yonsei University, Seoul 120-749, Republic of Korea

Received 15 October 2012; Revised 2 January 2013; Accepted 18 January 2013

Academic Editor: Ulrich Katscher

Copyright © 2013 Joonsung Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Anisotropic electrical properties can be found in biological tissues such as muscles and nerves. Conductivity tensor is a simplified model to express the effective electrical anisotropic information and depends on the imaging resolution. The determination of the conductivity tensor should be based on Ohm's law. In other words, the measurement of partial information of current density and the electric fields should be made. Since the direct measurements of the electric field and the current density are difficult, we use MRI to measure their partial information such as B1 map; it measures circulating current density and circulating electric field. In this work, the ratio of the two circulating fields, termed circulating admittivity, is proposed as measures of the conductivity anisotropy at Larmor frequency. Given eigenvectors of the conductivity tensor, quantitative measurement of the eigenvalues can be achieved from circulating admittivity for special tissue models. Without eigenvectors, qualitative information of anisotropy still can be acquired from circulating admittivity. The limitation of the circulating admittivity is that at least two components of the magnetic fields should be measured to capture anisotropic information.

1. Introduction

Noninvasive measurement of electrical properties for biological tissues can be useful in EEG/MEG and electromagnetic source imaging [1] and in providing diagnostics information about the physiological and pathological states of the tissues [2–5]. For isotropic conductivity, many approaches have been developed to measure the conductivity at low frequencies and at Larmor frequencies [6]. At low frequencies below 1 kHz, Magnetic Resonance Electrical Impedance Tomography (MREIT) [7] can probe the conductivity distribution. At Larmor frequencies of about 100 MHz, Magnetic Resonance Electrical Property Tomography (MREPT) [8, 9] measures both electric conductivity and permittivity distributions using measurements of positively rotating magnetic fields generated by transmit RF coil, + maps, from MRI.

Microscopically, the conductivity of the biological tissues could be isotropic. However, depending on the imaging resolution, the conductivity of an imaging voxel can be anisotropic. Macroscopically, in other words, if several tissues with different electrical properties are combined in the imaging voxel, the conductivity of the imaging voxel differs when measured in different directions so that it becomes anisotropic. Especially in biological tissues, anisotropic electrical conductivity can be found in muscles and nerves [2–5]. The conductivity tensor is a simplified model with three eigenvectors and three eigenvalues which can include these anisotropic cases.

The eigenvectors of the conductivity tensor at low frequencies (<1 kHz) can be inferred from a prior knowledge of the object or the diffusion tensor imaging [1, 10]. However, at Larmor frequency of about , there have been a few studies on measuring three eigenvectors and eigenvalues of the conductivity tensor. Recently, Katscher et al. [11] proposed a way to estimate partial information of conductivity anisotropy especially in the special case where two minimum eigenvalues are almost equal to zero. In this work, based on Katscher’s approach [11], we generalized and considered more possible cases. Based on numerical phantom simulations and phantom experiments, the performance was evaluated and the limitations and future directions were proposed.

2. Materials and Methods

2.1. Admittivity Tensor Model

Admittivity tensor denoted by is a simplified model for the electrical anisotropic information at the angular frequency , where and are conductivity and permittivity tensors, respectively.

The admittivity tensor can be represented with six parameters:

Expressing its eigenvectors (unit vectors) and its corresponding eigenvalues , the admittivity tensor can be also expressed as

We should note that three equations of Ohm’s law alone are insufficient to identify six unknown components of . However, the eigenvectors of the conductivity tensor could be estimated from prior knowledge of the object or can be determined by measuring the diffusion tensors [1] using MRI. Under the assumption that the eigenvectors of the admittivity tensor are known a priori, using the eigenvectors of the matrix , the conductivity tensor can be decomposed as

where and . In other words, .

2.2. Reconstruction of Admittivity Anisotropy Using Measured Magnetic Fields: Circulating Admittivity

Assume that the three eigenvectors are known. The effective admittivity, , in a voxel can be determined from Ohm’s law as follows:

However in MRI, and are hard to measure. Instead, partial knowledge of the magnetic fields can be acquired from mapping techniques [12–15]. Katscher et al. [11] extended the direct inversion method in MREPT [9] and proposed a way to estimate anisotropy of using the relationship between circulating currents and circulating electric fields over a surface, which can be estimated from the measured magnetic fields.

From time-harmonic Maxwell Equations, the relationship among the current density, electric fields, and the magnetic fields can be expressed as

where is a surface whose normal vector is .

Based on the work by Katscher et al. [11], we define the circulating admittivity, , as the ratio of the circulating currents and the circulating electric fields over the surface with rotating the normal vector :

where is parallel to the normal vector . In (6), the relationship between the curl integral of and the surface integral of hold for homogeneous region of admittivity and can generate artifacts at tissue boundaries [16].

In this work, we investigated the relationship between the circulating admittivity and the admittivity tensor for simple cases. First, if the admittivity tensor is isotropic, the circulating admittivity is also isotropic and equal to the isotropic admittivity. Second, if two eigenvalues of the admittivity tensor are the same, , the eigenvalue can be determined from the circulating admittivity over a surface with a normal vector, , perpendicular to the eigenvectors corresponding to the two eigenvalues, . That is

since the normal vector is perpendicular to the vector over a line and the two eigenvalues are the same .

In addition, as it is considered in [9], if two smaller eigenvalues of the admittivity tensor, , are equal to zero, the largest eigenvalue, , can be determined from the circulating admittivity directly:

As a combination of last two cases, if the largest eigenvalue is much bigger than the two smaller eigenvalues of the admittivity tensor and the two smaller eigenvalues are the same, by measuring the circulating admittivity for several directions of the normal vector, all three eigenvalues could be estimated.

2.3. Unknown Directions of Eigenvectors: Effective Admittivity Map (EAM) and Circulating Admittivity Map (CAM)

In the previous section, we determined the admittivity tensor under the assumption that the eigenvectors of the admittivity tensor were known a priori. Even without the prior knowledge of the eigenvectors, we can still provide a qualitative measurement of anisotropy by computing the dependency on the normal vector, , in (6). For a qualitative measurement of the anisotropy, we define the effective admittivity map, , and the circulating admittivity map, , that describe the distributions of effective admittivity and circulating admittivity over the normal vector, respectively:

As shown in Figure 7, in a later section, the effective admittivity map (EAM) and the circulating admittivity map (CAM) can be drawn with the use of two angles, and which describe the direction of the normal vector. That is, the normal vector was initially located at the positive -axis, was rotated along -axis by , and then was rotated along -axis by .

2.4. Numerical Simulation: Numerical Phantom Model with Anisotropic Effective Admittivity

For numerical evaluation, a numerical phantom with anisotropic effective admittivity can be generated using periodic binary medium. According to homogenization theory, anisotropy can be derived from pointwise admittivity, , that is distributed periodically with respect to the -variable:

where is a constant depending on the binary medium, is a large positive integer, and is the largest integer not greater than .

In this numerical experiment the imaging subject is the box . We divided the domain into two subdomains and . In the admittivity is homogeneous with the value , where the permittivity with the permittivity in the free space. In 33 layers were stacked alternatively with the thickness of 2 mm, the admittivity value and the thickness of 1 mm, and the admittivity value in (10) with setting to be . Figure 1(a) shows the construction of the imaging object.

(a)

(b)

Figure 1: (a) Numerical phantom with anisotropic effective admittivity by stacking periodic binary medium, (b) placement of the imaging object inside the RF coil.

Driven by a birdcage coil at 3T ( MHz) as shown in Figure 1(b), the electric fields, the magnetic fields, and the current densities in microscale were calculated using finite-difference time domain (FDTD) numerical simulations using REMCOM (REMCOM, State College, PA) with the resolution of . Then, to determine the effective admittivity in macroscale, that is, the ratio of the ensemble mean current density to the ensemble mean electrical field, three-dimensional Gaussian filter with the size of and the standard deviation of , which increases the effective voxel size from to , was applied to the simulated fields.

2.5. MRI Experiments

Two phantoms with anisotropic admittivity were generated using straws as shown in Figure 2. As a comparison, one water phantom with isotropic admittivity was made without straws. Three phantoms were cylindrical with the height of and the radius of . The diameters of the straws are for straw and for straw . All three phantoms were filled with the saline water of 0.35 M NaCl concentration as shown in Figure 2.

(a)

(b)

Figure 2: (a) Two straw phantoms before filling saline water, (b) three phantoms filled with saline water.

Using a single-channel transreceive head coil, MR images were measured. The phantoms were located at the isocenter of the coil with the straw orientation of left-right. Only the phase of was measured and the circulating conductivity was determined by the phase-based approximation in MREPT [17]. Three-dimensional balanced steady-state free procession (bSSFP) was acquired with resolution of , field of view (FOV) of , and image size of . The other imaging parameters were the flip angle of 30 degrees, TE of 1.8 ms, TR of 3.6 ms, and the scan time of 5 minutes with 27 averages. All measurements were performed on a 3T Siemens Tim Trio scanner. The phase of was estimated as the half of the measured phase of the image [17].

3. Results

3.1. The Effective Admittivity of the Numerical Phantom

Based on (4), by dividing the filtered current densities and filtered electric fields, the effective anisotropic material can be acquired. Figures 3(a) and 3(d) illustrate the conductivity and the relative permittivity in microscale of resolution in the slice at that we set in this simulation, respectively. The effective conductivities in macroscale , are shown in Figures 3(b) and 3(c) and the effective relative permittivities are shown in Figures 3(e) and 3(f). For the subdomain of homogeneous tissue, , the conductivity and relative permittivity are constant and the same microscopically and macroscopically. The effective conductivity and relative permittivity are almost constant for the subdomain of the alternating layers of tissues, except some distortions, we think, due to simulation errors. Since in this experiment is the same as , was not shown here.

(a) True in microscale

(b) Effective in macroscale

(c) Effective in macroscale

(d) True in microscale

(e) Effective in macroscale

(f) Effective in macroscale

Figure 3: (a) and (d) depict the true conductivity and the relative permittivity values in the microscopic scale. (b)–(f) illustrate the true effective , , , and respectively.

3.2. Observation of Anisotropy Using Circulating Admittivity (6)

In MRI, only partial information of the magnetic fields can be measured. Using conventional single-transmit channel MR scanner, the circularly polarized component of the magnetic fields, , can be measured but the other two components, , are hard to measure. Using a specialized scanner, parallel transmit system, the anticircularly polarized component could be measured [15]. Here, we considered two cases: using and , using only . For the computation of the circulating admittivity using a partial information of the magnetic fields, the unmeasured magnetic fields were assumed to be zero.

Figure 4 illustrates the values of circulating conductivity and the relative permittivity, and , with and at the slice, . The surface for the integration (6) was chosen as a plane with the size of and the normal vector, , of , . In this case, the two components of the magnetic fields, and , were assumed to be known and used to reconstruct the circulating admittivity. As derived in (7), the circulating admittivity with the normal vector of , which is perpendicular on the two eigenvectors of the effective admittivity tensor with the same eigenvalues, is close to the effective admittivity in direction except at the tissue boundaries. For the normal vector of , the circulating admittivity is a weighted average of the effective admittivities, , , . Thus, the circulating conductivity with the normal vector of is bigger than the effective conductivity corresponding to the smallest eigenvalue of the effective admittivity tensor, , shown in Figure 3(c).

(a)

(b)

(c)

(d)

Figure 4: Circulating conductivity and relative permittivity derived from and : (a) , (b) , (c) , and (d) .

However, as shown in Figures 5 and 6, if only one component, , is available, the dependency on the normal vector was lost in the circulating admittivity.

(a)

(b)

(c)

(d)

Figure 5: Circulating conductivity and relative permittivity derived from only (a) , (b) , (c) , and (d) .

(a)

(b)

Figure 6: Circulating conductivity derived with only phase of , (a) , (b) .

(a) Saline water:

(b) Saline water:

(c) Straw phantom 1:

(d) Straw phantom 1:

(e) Straw phantom 2:

(f) Straw phantom 2:

Figure 7: Experiment results: circulating conductivity of phantoms with measured , Coronal Slice, (a) of saline water, (b) of saline water, (c) of straw phantom 1, (d) of straw phantom 1, (e) of straw phantom 2, (f) of straw phantom 2.

The circulating conductivities were determined from phantom experiments in which only the phase of is measurable. As shown in Figure 7, the dependency on the normal vector was also lost in the experimental results. However, the decrease of the conductivity due to the plastic straws was observed.

3.3. Distribution of the Admittivity: Circulating Admittivity Map (CAM)

Using simulated magnetic fields, and , the circulating admittivity map (CAM) was computed for the numerical phantom. As a comparison, the effective conductivity map was also computed using simulated current density and electric fields. The and , which determined the direction of the normal vector , varied from and by one degree. Figure 8 shows the effective conductivity maps, the real part of the EAM, and the circulating conductivity maps, the real part of the CAM, of one anisotropic voxel located at , , and and one isotropic voxel located at , , and .

(a)

(b)

(c)

(d)

Figure 8: Effective conductivity maps and circulating conductivity maps of one anisotropic voxel and one isotropic voxel: (a) effective conductivity map of anisotropic voxel, , , , (b) circulating conductivity map of anisotropic voxel, , , , (c) effective conductivity map of isotropic voxel, , , , (d) circulating conductivity map of isotropic voxel, , , .

For the isotropic voxel, the values of the effective conductivity map and circulating conductivity map were equal to the conductivity of the tissue. For the anisotropic voxel, the circulating conductivity map is also uniform along the direction since only and were used. In both effective conductivity maps and circulating conductivity maps, the direction that maximizes or minimizes the conductivity value does not match with any eigenvector of the admittivity tensor. Thus, with CAM alone, the eigenvectors of the admittivity tensor may be hard to determine and thus a quantitative measurement of the eigenvalues may be hard.

For a qualitative analysis, at each voxel, the maximum value, the minimum value, and the ratio of the maximum to the minimum of the circulating conductivity maps were computed with and . As shown in Figure 9, in this case, the maximum values were almost constant over anisotropic tissues, but the minimum values were not constant. The ratio of the maximum to minimum, which could be used as a qualitative measurement of the anisotropy, was not constant over anisotropic tissues and was smaller than the ratio of the maximum eigenvalue, , to the minimum eigenvalue, , of the conductivity tensor; that is, using CAM, the contrast between isotropic and anisotropic tissues was reduced. However, CAM still can separate anisotropic tissues from isotropic tissues without knowing the eigenvectors of the admittivity tensor.

(a) Maximum

(b) Minimum

(c) Max/Min

(d)

Figure 9: The distribution of the circulating conductivity: (a) the maximum values for the circulating conductivity maps, (b) the minimum values for the circulating conductivity maps, (c) the ratio of the maximum to the minimum, (d) the ratio of two eigenvalues in the conductivity tensor.

4. Discussion

Conductivity tensor is a simplified anisotropy model. Given three eigenvectors, the tensor can be estimated if the electric current densities and electric fields can be measured. In MRI, however, electric fields are hard to be measure without knowing or estimating the conductivity and permittivity of tissues. In this work, using MREPT formulae, the circulating admittivity is proposed as a measure to analyze the anisotropy of the tissues. Circulating admittivity was defined as the ratio of circulating current densities to the circulating electric fields, which can be determined from the magnetic fields. We did not fully investigate, but we derived the relationship between the admittivity tensor for special cases. Using numerical phantom simulations, we verified the relationship for the first two cases: isotropic tissues and two eigenvalues of the admittivity tensor are the same. As a future work, more realistic cases would be considered.

In this work, to deal with unknown eigenvectors, the circulating admittivity map (CAM) was proposed as a qualitative measure. The ratio of the maximum to the minimum conductivity was reduced but still anisotropic tissues can be separated from isotropic tissues.

In the conventional single-transmit channel MR scanner, the circularly polarized magnetic field, , can be measured by mapping methods, but the other two components are hard to measure. If only one measurement of magnetic fields, is available, the anisotropic information is lost in the estimate of the circulating admittivity. Even if only one component of the magnetic fields, , can be measured, the anisotropic information can be acquired by measuring several by rotating the object with respect to transmit coil. In Figures 10 and 11, the circulating conductivities only with for two positions of the object, that is, initial position and rotated, are shown. For simulation data shown in Figure 10, the circulating conductivities were computed with both magnitude and phase of or with only phase of . For experimental data shown in Figure 11, only the phase of the was used. For isotropic tissues, circulating conductivity was not related to the position, but for anisotropic tissues, the circulating conductivity at the center of the phantom was changed. At the boundary of the phantom, very high or negative, especially at the top and the bottom of the phantom after rotating , conductivity values were observed. We think that boundary artifacts [16] created at the air-water boundary were spread inside the phantom due to the spatial filtering used to reduce the noise in the conductivity estimates and thus our conductivity estimates at the boundary of the phantom were not reliable.

(a)

(b)

(c)

(d)

Figure 10: Simulation results: circulating conductivity of simulated phantoms with measured , coronal slice: (a) initial position using both magnitude and phase of , (b) rotated using both magnitude and phase of , (c) initial position using only the phase of , and (d) rotated using only the phase of .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: Experiment results: circulating conductivity of phantoms with measured , coronal slice, (a) saline water: initial position, (b) saline water: rotated, (c) straw phantom 1: initial position, (d) straw phantom 1: rotated, (e) straw phantom 2: initial position, (f) straw phantom 2: rotated.

5. Conclusions

Noninvasive measurement of conductivity tensor at Larmor frequency could be achieved using MRI. Using measured maps from MRI, circulating current density and circulating electric fields can be estimated. In this work, the ratio of the two, called circulating admittivity, was proposed as measure of the conductivity anisotropy at Larmor frequency. Given eigenvectors of the conductivity tensor, quantitative measurement of the eigenvalues can be achieved from circulating admittivity for special tissue models. Without eigenvectors, qualitative information of anisotropy still can be acquired from the distribution of the circulating admittivity. The limitation of the circulating admittivity is that the anisotropic information is lost if only one component of the magnetic field is available. At least, an additional acquisition, either by rotating the object or some other scheme, needs to be performed for anisotropic information.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) Grant funded by Korea Government (MEST) (no. 2012-009903 to J. Lee, N. Choi, and D. H. Kim and no. 2011-8-1782 to Y. Song and J. K. Seo). The authors acknowledge Ulrich Katscher for beneficial discussions.

References

D. S. Tuch, V. J. Wedeen, A. M. Dale, J. S. George, and J. W. Belliveau, “Conductivity tensor mapping of the human brain using diffusion tensor MRI,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 20, pp. 11697–11701, 2001.

C. Gabriel, S. Gabriel, and E. Corthout, “The dielectric properties of biological tissues: I. Literature survey,” Physics in Medicine and Biology, vol. 41, no. 11, pp. 2231–2249, 1996.

S. Gabriel, R. W. Lau, and C. Gabriel, “The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz,” Physics in Medicine and Biology, vol. 41, no. 11, pp. 2251–2269, 1996.

L. A. Geddes and L. E. Baker, “The specific resistance of biological material—a compendium of data for the biomedical engineer and physiologist,” Medical and Biological Engineering, vol. 5, no. 3, pp. 271–293, 1967.

S. Grimnes and O. G. Martinsen, Bioimpedance and Bioelectricity Basics, Academic Press, London, UK, 2000.

J. K. Seo, D. H. Kim, J. Lee, O. I. Kwon, S. Sajib, and E. J. Woo, “Electrical tissue property imaging using MRI at dc and Larmor frequency,” Inverse Problems, vol. 28, Article ID 084002, 2012.

H. S. Khang, B. I. Lee, S. H. Oh et al., “J-substitution algorithm in magnetic resonance electrical impedance tomography (MREIT): phantom experiments for static resistivity images,” IEEE Transactions on Medical Imaging, vol. 21, no. 6, pp. 695–702, 2002.

E. M. Haacke, L. S. Pepropoulos, E. W. Nilges, and D. H. Wu, “Extraction of conductivity and permittivity using magnetic resonance imaging,” Physics in Medicine and Biology, vol. 36, no. 6, pp. 723–733, 1991.

U. Katscher, T. Voigt, C. Findeklee, P. Vernickel, K. Nehrke, and O. Dossel, “Determination of electrical conductivity and local SAR via B1 mapping,” IEEE Transactions on Medical Imaging, vol. 28, pp. 1365–1374, 2009.

A. K. Sen and S. Torquato, “Effective conductivity of anisotropic two-phase composite media,” Physical Review B, vol. 39, no. 7, pp. 4504–4515, 1989.

U. Katscher, T. Voigt, and C. Findeklee, “Estimation of anisotropy of electric conductivity via B1 mapping,” in Proceedings of the 18th Annual Meeting of ISMRM, p. 2866, Stockholm, Sweden, 2010.

S. Akoka, F. Franconi, F. Seguin, and A. Le Pape, “Radiofrequency map of an NMR coil by imaging,” Magnetic Resonance Imaging, vol. 11, no. 3, pp. 437–441, 1993.

R. Stollberger and P. Wach, “Imaging of the active B1 eld in vivo,” Magnetic Resonance in Medicine, vol. 35, pp. 246–251, 1996.

V. L. Yarnykh, “Actual flip-angle imaging in the pulsed steady state: a method for rapid three-dimensional mapping of the transmitted radiofrequency field,” Magnetic Resonance in Medicine, vol. 57, no. 1, pp. 192–200, 2007.

X. Zhang, P. F. V. de Moortele, S. Schmitter, and B. He, “Complex B1 mapping and electrical properties imaging of the human brain using a 16-channel transceiver coil at 7T,” Magnetic Resonance in Medicine, 2012.

J. Seo, M. Kim, J. Lee et al., “Error analysis of non-constant admittivity for MR-based electric property imaging,” IEEE Transactions on Medical Imaging, vol. 31, pp. 430–437, 2011.

T. Voigt, U. Katscher, and O. Doessel, “Quantitative conductivity and permittivity imaging of the human brain using electric properties tomography,” Magnetic Resonance in Medicine, vol. 66, no. 2, pp. 456–466, 2011.