# Recent Progress and Future Challenges in MR Electric Properties Tomography

Authors: Ulrich Katscher; Dong-Hyun Kim; Jin Keun Seo
Journal: Computational and Mathematical Methods in Medicine
 (2013)

---

## Abstract

MR Electric Properties Tomography (EPT) is a lately developed medical imaging modality capable of visualizing both conductivity and permittivity of the patient at the Larmor frequency usingB1 maps. The paper discusses the development
of EPT reconstructions, EPT sequences, EPT experiments, and challenging issues of EPT.


---

## Body

## 1. Introduction

The knowledge of electrical tissue properties is expected to be beneficial for clinical diagnosis, therapy monitoring, and RF patient safety. Electrical tissue properties can be described by conductivityσ and permittivity ϵ, and they exhibit frequency-dependent behavior since tissues are heterogenous substances comprising insulating cell membranes and conducting electrolytes. Visualization of frequency-dependent conductivity and permittivity distribution in the range from almost dc to hundreds of MHz may expand our ability to provide diagnostic information about the physiological and pathological state of tissues and organs [1].Due to its electromagnetic background, MRI would be a top candidate for delivering this desired knowledge of electrical tissue properties. The complex permittivityκ:=ϵ-i(σ/ω) at an angular frequency ω (assumed to be below microwave range) can be probed by the time-harmonic magnetic field H=(Hx,Hy,Hz) through the following arrangement of time-Maxwell equations (so-called Helmholtz equation):
(1)-∇2H=∇κκ×[∇×H]+ω2μκH,
where μ is the magnetic permeability and κ is assumed to be isotropic. Here, the corresponding time-varying field is Re{Heiωt}.At frequencies below 1 kHz, Joy et al. in 1989 [2] introduced MR current density imaging (MRCDI) which aims to provide noninvasive visualization of current density J=∇×H inside a body by externally injecting dc current using a pair of surface electrodes and measuring the induced magnetic field H using MRI. In MRCDI, the induced current density J produces a change of the main dc magnetic field, and Hz is a measurable quantity by MRI since it alters the MR phase image. Hence, from MRCDI, obtaining an image of J=∇×H requires mechanical rotation of the subject inside MRI to measure all three components of H [3]. In 1994 [4], MR electrical impedance tomography (MREIT) was proposed to perform the conductivity imaging at dc using the MRCDI technique. In 2001 [5], an imaging technique of MREIT without mechanical rotation, called harmonic Bz algorithm, was developed to provide both conductivity image and current density image. After invention of the harmonic Bz algorithm, MREIT has advanced rapidly [6, 7]. However, it still remains a technical problem to reduce the injection current down to a level for routine clinical use while maintaining the spatial resolution of the resulting conductivity images.At frequencies above 1 MHz, the currents required to image electric properties have not necessarily to be injected by external surface electrodes as in MREIT. Alternatively, eddy currents can be induced by applying magnetic RF fields, avoiding the sensation of pain frequently connected with external current injection. Since magnetic RF fields are an inherent component of MRI, the desired currents can be created by standard MR systems and standard MR sequences. The resulting imaging technique, called Electric Properties Tomography (EPT) [8–11], is the subject of this paper. The basic idea of EPT is that  the electric properties of the patient distort B1, the component of the magnetic RF field responsible for spin excitation. Measuring this distorted B1 by B1 mapping techniques (see, e.g., [12–17]) allows to reconstruct the electric properties causing the observed distortions. This basic idea of EPT is illustrated in Figure 1, depicting the change of the phase of the magnetic RF field due to a brain tumor with a diameter of 1 cm. This phase change increases with the applied frequency, that is, the main field B0 of the MR system used. This phase change also increases with the conductivity of the tumor. The obtained phase changes of several degrees are in a measurable range. Phase changes further increase with increasing tumor size.Figure 1
Change of magnetic RF field phase due to a brain tumor with a diameter of 1 cm. This phase change increases with the applied frequency (corresponding to the system’s main fieldB0) as well as with the conductivity of the tumor. A mean conductivity of 0.4 S/m was assumed for the healthy parts of the brain, and thus, negative phase changes appear for assumed tumor conductivities below 0.4 S/m.RF currents induced in the tissue cause not only a distortion ofB1, which is utilized for EPT as discussed above, but also off-resonance effects. As in MRCDI, these off-resonance effects can be utilized to measure the current density along the direction of B0, called RF-CDI. RF-CDI is not part of this paper, and interested readers might be referred to the corresponding literature [18–20]. The following chapters review the development of EPT reconstructions, EPT sequences, and EPT experiments.


## 2. Development of EPT Reconstruction

Without calling the approach EPT, the first mentioning of EPT was in the early nineties by Haacke et al. [8]. It was suggested to calculate both conductivity σ and permittivity ϵ via the homogeneous Helmholtz equation
(2)κ(r)=-1ω2μ0∇2H+(r)H+(r),
where r:=(x,y,z) and H+:=(Hx+iHy)/2, the positive circularly polarized component of the magnetic field corresponding to the RF transmit field. Here, the main magnetic field is B0=-B0z^, and the corresponding complex rotating vector is a+=x^-iy^, where x^=(1,0,0), y^=(0,1,0),  z^=(0,0,1), and B0>0. From (2), σ and ϵ can be expressed as
(3)σ=1ωμ0Im{∇2H+H+},ϵ=-1ω2μ0Re{∇2H+H+}.Equation (2) is derived from (1) with the following assumptions. (A1)
A locally (“piecewise”) constantκ(r), that is; ∇κ(r)=0. This assumption has severe consequences, which has to be discussed extensively later on.(A2)
A constantμ(r)=μ0=4π×10-7, the magnetic permeability of the free space. This assumption is fairly fulfilled in the human body and does not require further discussion.(A3)
Isotropicκ. At the Larmor frequency of 128 MHz at 3 T, anisotropy is small in most tissues but yields an interesting niche application of EPT.(A4)
|H+| is assumed to be larger than zero to avoid singularities in (2). This is obviously the case in areas of nonzero MR signal.An invaluable advantage is the cancelation of the scaling ofH+ in the numerator and denominator of the expression (2). Given this cancelation, (2) yields absolute values of κ, even for arbitrary scaling of H+. This feature ennobles EPT to the class of quantitative MR methods, opening the chance to directly compare κ between different patients and different lesions. However, in the named publication [8], EPT was not pursued further due to “spurious phase effects unrelated to RF penetration which makes a simple extraction difficult.” Although the mentioned spurious phase effects predominantly belonged to imperfections of the MR systems back in the early nineties, which greatly reduced since then, phase effects unrelated to RF penetration are still one of the major issues for EPT, particularly for in vivo measurements. Instead of the mentioned “simple extraction” of κ, Haacke et al. developed a heterogeneous layer model as a workaround for the observed spurious phase effects [8].The first successful application of EPT (still not called EPT) described by Wen [9] is dated more than 10 years after Haacke’s initial article. In this conference abstract, the expressions (3) are used again. Two further observations are mentioned by Wen, which later on turned out to be of central importance for EPT.(O1)
Computations of (3) require both magnitude and phase of H+=|H+|eiϕ+. Unfortunately, only the magnitude |H+| enters the MR signal in a nonlinear way; for each nominal flip angle α of the sequence, the following MR signal is measured:
(4)S(r)=V1M0(r)H-(r)exp(iϕ+(r))sin(V2α|H+(r)|),whereM0 is the MR magnitude image containing relaxation and spin density effects, H-=(Hx-iHy)/2, and V1,V2 system-dependent constants [21]. Thus, |H+| has the chance to be measured exactly, assuming the ideal working of corresponding B1 mapping techniques (see, e.g., [12–17]). However, the phase ϕ+ is difficult to be determined exactly. The phase of a standard MR image is always the superposition of ϕ+ with its counterpart of the RF reception, ϕ- from H- and, thus, is called “transceive phase” ϕ0=ϕ++ϕ-. In a standard MR system with a quadrature body coil (QBC), the polarization of this coil is switched from RF transmission to RF reception for the sake of optimizing SNR. Wen observed that the resulting ϕ~- of the switched QBC closely resembles ϕ+ [9]. Consequently, a rough approximation of ϕ+ can be obtained by the transceive phase ϕ~0=ϕ++ϕ~- of this setup
(5)ϕ+≈ϕ~0/2=(ϕ++ϕ~-)/2sometimes called “transceive phase assumption.”(O2)
To the leading order, the conductivity response affects the phase of the RF field, while the permittivity response affects the magnitude of the field. Thus,σ can be estimated by applying (2) only to ϕ+, now called “phase-based EPT”. Accordingly, ϵ can be estimated by applying (2) only to |H+|, now called “magnitude-based EPT.”In principle, these two observations pave the way to a clinically feasible EPT. Moreover, the first successful phantom and ex vivo experiments are presented in [9], which will be topic of a later section dedicated to experimental EPT results. After this publication, Wen left the topic of EPT, and again EPT was not further pursued for years.Systematic research on EPT started in 2009 with [10]. This publication is based on the following expression which comes from a modified Helmholtz equation:
(6)κ(r)=∮∂Ar‍∇×H(r′)·dlμ0ω2∫Ar‍H(r′)·dS,
where Ar is an arbitrarily oriented area centered at r with its boundary ∂Ar, dl the line element, and dS the surface element. Equation (6) can be viewed as dividing Ampere’s law by Faraday’s law after suitably integrating these laws. In contrast to (2), the expression (6) requires all three components of the magnetic field H. Since H- and Hz are not measurable directly, H-=0 and Hz=0 were assumed for the applied QBC [10]. Nevertheless, a couple of basic findings studied with this modified Helmholtz equation is valid for all kinds of EPT reconstructions. First of these basic findings, the violation of assumption (A1), the locally constant κ, leads to severe artifacts along boundaries between compartments of different κ. These artifacts are typically strong oscillations (under/overshooting) of the reconstructed κ as shown in Figure 3. This is predominantly a question of the numerical implementation of the calculus operations of the EPT equation applied, which always involves a number of voxels in the neighborhood of the target voxel to be reconstructed. This so-called kernel of involved voxels, regardless of the actual calculus operation, causes the mentioned oscillations as soon as it contains voxels of different κ. Thus, as shown in [10], lowering the kernel size narrows the oscillations, and a kernel of minimal size seems to be optimal. As shown in [10], lowering the kernel size also degrades the noise figure. This touches another basic finding: the second derivative, as enclosed explicitly or implicitly in all EPT equations, tends to significantly enhance the noise in the measured H+. For minimal kernel size, the SNR of the reconstructed κ is far below the SNR of the input H+ [10, 22]. This initiated a bunch of activities to find the optimal tradeoff between artifacts and SNR or to find suitable workarounds, as discussed later.Besides, for the phantom used in [10], the violation of the transceive phase assumption (O1) was proven to be far below 1°, and thus errors arising from this violation are expected to be lower than errors from other sources. Last but not least, in [10], EPT was first applied to estimate local SAR.Another milestone of EPT reconstruction is given by [11]. In contrast to the factual, “physical” modification of (2) in [10], the modification in [11] is the result of reformatting the original equation (2) by suitably integrating its numerator and denominator:
(7)κ(r)=-∫∂Vr‍∇[|H+(r′)|eiϕ+(r′)]·dSμ0ω2∫Vr‍[|H+(r′)|eiϕ+(r′)]dV,
where Vr is a volume centered at r with its boundary ∂Vr and dV is the volume element. The denominator of (7) averages the denominator of (2) over a certain volume, while the numerator integrates the normal derivative to the surface of this volume. From a numerical point of view, the minimal kernel size of (7) is larger than the minimal kernel size of (2), looking for a better solution of the above-mentioned tradeoff between boundary artifacts and noise level.  A systematic comparison of the behavior of (2) and (7) is conducted in [23]. Equation (7) can be rewritten with separate real and imaginary part
(8)κ=-1μ0ω2[(∇2|H+||H+|-|∇ϕ+|2)+i(2∇ln|H+|·∇ϕ++∇2ϕ+)(∇2|H+||H+|-|∇ϕ+|2)].As already indicated in [9], observation (O2), the assumption |∇2ϕ+|≫2|∇ln|H+|·∇ϕ+| yields phase-based EPT for conductivity imaging
(9)σ(r)≈∇2ϕ+(r)μ0ω,
and the assumption |∇2|H+|/|H+||≫|∇ϕ+|2 yields magnitude-based EPT for permittivity imaging
(10)ϵ(r)≈-∇2|H+(r)|μ0ω2|H+(r)|.Reference [11] derives analytically the related errors of these approaches and investigates systematically their feasibility. For typical κ of human tissue at 1.5 T or 3 T, the error introduced by the expressions (9)-(10) is only of the order of 10%. As can be seen directly from the corresponding error terms, phase-based EPT always yields too high conductivities, and magnitude-based EPT always yields too low permittivities. Thus, the discussed errors can be hidden by mapping not σ or ϵ, but |κ| as shown in [24].Phase-based EPT reveals two features invaluable for clinical applications.(F1)
The linearity of the expression (9) supersedes the QBC transceive-phase assumption (5), allowing arbitrary combinations of RF transmit and receive coils. The resulting transceive phase, containing ϕ+ and ϕ- from different RF fields, still yields σ via
(11)σ=∇2ϕ0/2μ0ω=∇2(ϕ++ϕ-)2μ0ω=12(∇2ϕ+μ0ω+∇2ϕ-μ0ω)=2σ2since (9) can be based on ϕ+,  ϕ-, or any phase fulfilling Maxwell’s equations, leading to the same σ as long as ∇|H+|=0 and ∇|H-|=0 are fulfilled.(F2)
Skipping the need of mapping (the magnitude of)B1, which is typically a rather lengthy scan, significantly speeds up the scan time required for EPT. As discussed below, even real-time conductivity measurements seem to be possible [25]. Moreover, it opens the chance that the conductivity can be obtained via sequences which are not primarily driven for EPT, just reusing the transceive phase which usually comes for free with every MR sequence.Shortly after [11], which was performed at B0=1.5T, phase-based EPT has been confirmed at B0=7T [26]. The impact of B0 on EPT and the related question of B0 optimal for EPT turned out to be a nontrivial task [27]. Obviously, higher SNR can be expected with increasing B0. This advantage is counterbalanced by the increasing violation of the assumption ∇|H+|=0 for (9) or the QBC transceive phase assumption (5), respectively. Although not explicitly stated by the authors of [27], the optimal tradeoff between SNR and reconstruction accuracy seems to be given at B0=3T for conductivity imaging. For permittivity imaging, the violation of the assumption ∇ϕ+=0 does not increase with B0, and the highest available B0 seems to be optimal. This trend is further emphasized by the different powers of ω in (3).Since 2009, research on EPT spread out, and more and more groups started to investigate different aspects of EPT [26–31]. Typically, the original expression (2) has been used in these studies.The problem of separatingϕ+ and ϕ- from the transceive phase has been solved analytically by an approach sometimes called “Local Maxwell Tomography” (LMT) [32–34]. LMT is based on the insight that the reconstructed κ must not depend on the applied RF coil [31]. This is a particularly useful insight given a system with multiple, independent RF transmit channels (see, e.g., [35, 36]). Such multitransmit systems, designed primarily for RF shimming at high B0, offer the chance to determine κ separately by each single TX channel or any arbitrary combination of TX channels. Differences of the reconstructed κ based on different RF excitations can be related to a violated transceive phase assumption. For instance, two EPT reconstructions σn and σm can be compared based on different TX channels n and m, but same receive channel with phase ϕ-, yielding [34]
(12)σn(r)-σm(r)=1μ0ω∇ϕ-(r)·∇ln|Hn+(r)||Hm+(r)|.This allows, first, the determination of the unknownϕ- and, subsequently, the straightforward determination of the unknown ϕn+ and ϕm+. The central idea of comparing reconstruction results of two or more different TX channels can also be utilized to exactly distinguish |H-| from the spin magnetization M0, another hitherto unexplored possibility [33]. Unfortunately, the numerical effort to solve the related equations is high. It shall be clarified in future studies, if the accuracy of the obtained ϕ+ is high enough to improve reconstruction results hitherto obtained in the framework of the discussed phase assumptions. The explicit knowledge of ϕ+,  ϕ-, and |H-| seems to be more urgent for the determination of local SAR [34].Thus, at least on a theoretical basis, the issue of phase determination is figured out, and the last remaining issue for the EPT reconstruction is the treatment of nonconstantκ (see Figure 3). Equation (2) can be viewed as a simplified version of (1) using the assumption (A1) of ∇κ=0 so that the effect of (∇κ/κ)×[∇×H] in (1) is neglected. Here, κ might change continuously or discontinuously across boundaries of compartments with different values of κ. Both types of changing κ and the related errors have been analyzed thoroughly in [37]. In (1), the occurring partial derivatives of κ act as additional unknowns. It was suggested that these additional unknowns can be solved using the described comparison of different RF excitations [31]. Comparing (simulations of) two different RF excitations and therefore calling the algorithm “dual excitation algorithm,” the typical boundary artifacts were significantly reduced. Alternatively, it has been proposed to multiply (1) with ∇×H yielding [38]
(13)κ(r)=-∇2H(r)·(∇×H(r))ω2μH(r)·(∇×H(r)).
This equation has the big advantage of removing the term involving ∇κ in (1). However, computation of (13) requires all spatial components of H and |H·∇×H|>0. The authors of [31, 38] propose to assume H-=0 and Hz=0 as it was done in connection with (6) [10]. A challenging problem is to find an exact relation between κ and H+ from the full equation (1) in such a way that κ can be computed robustly and efficiently using H+ only.From a numerical point of view, the occurring oscillations along compartment boundaries are a question of finite kernel size as explained above. Thus, instead of solving (1), boundary artifacts can be avoided by image segmentation prior to reconstruction and performing separate reconstructions on the different compartments. This was demonstrated in [39] in the framework of breast EPT, where different conductivities of highly nested fatty and ductile tissue spoil a standard EPT reconstruction completely. The image segmentation can be used for shaping the applied kernel locally to the current tissue type, as well as for locally restricting a subsequently applied smoothing filter [39]. The mentioned image segmentation can be based on standard T1/T2 weighted images, which implies that same T1/T2 coincides with same σ. This of course is not automatically fulfilled; however, violations of this assumption are expected to occur much less than discontinuous σ across T1/T2 boundaries.This chapter ends with a brief discussion of anisotropicκ violating (A3). Anisotropic κ (a rank-2 tensor) can be characterized by its eigenvectors v1,v2, and v3 (unit vectors) and its corresponding eigenvalues κ1,κ2, and κ3, respectively:
(14)κ=(v1v2v3)(κ1000κ2000κ3)(v1v2v3).
Measuring anisotropy of the tissue conductivity, characterizing the underlying cell structure, might increase diagnostic information. In vivo, anisotropic conductivities can be found in tissue with preferred cell direction, for example, in muscles and nerves. However, one has to keep in mind that anisotropy of κ is expected to decrease with increasing ω, and anisotropy at Larmor frequency could be negligible. Nevertheless, it was pointed out that varying the orientation of the integration area A in (6) reflects the degree of anisotropy [40, 41]. To be precise, let A(n) be a disk with its unit normal vector n. Some anisotropic structure of κ can be observed by displaying the following quantity on the sphere S2 [41]:
(15)Θ(n):=∮∂A(n)‍∇×H·dlμ0ω2∫A(n)‍H·dS,n∈S2.
A perfectly isotropic κ should not depend on the direction n of A(n). On the other hand, the reconstructed κ should show a minimum for A(n) perpendicular to the (main) direction of a (strongly) anisotropic κ. Experimental results of a straw phantom confirmed this concept [40]. However, according to corresponding simulations [41], the concept shall not work without proper knowledge of H-.


## 3. EPT Sequence

### 3.1. MeasuringB1 Magnitude

As seen in (4), the transmit magnitude |H+|, required for EPT, can be measured in a straightforward manner due to its nonlinear impact on the MR signal. In this framework, numerous techniques for |H+| mapping (B1 mapping) are published (see, e.g., [12–17]). In principle, EPT can be based on any B1 mapping method. The accuracy of EPT depends on the accuracy of this mapping; that is, the most accurate |H+| mapping method leads to the most accurate EPT results. Studies looking for the optimum B1 mapping technique, independent of EPT, have been published elsewhere.


### 3.2. MeasuringB1 Phase

As discussed in the previous chapter, determination of theB1 phase ϕ+ always starts with the measurement of the transceive phase. One of the main issues of EPT, as already pointed out in [8], is the contamination of the transceive phase by unwanted phase contributions unrelated to RF penetration. The following steps from (S1) to (S4) describes how to get a transceive phase usable for EPT reconstruction. (S1)
The transceive must not contain any contributions fromB0, that is, any off-resonance effects. The easiest way to exclude off-resonance effects is the use of refocusing pulses, that is, sequences based on spin echoes (SE), like fast spin echo sequences, turbo spin echo sequences, and so on. In contrast, the transceive phase of field-echo based sequences includes off-resonance effects. In this case, these effects can be removed by any kind of B0 mapping. In the easiest way, the phase can be measured at two different TE and extrapolated back to TE = 0 [26]. Also more sophisticated B0 maps can be applied, for example, obtained in the framework of the Dixon techniques (see, e.g., [42]). On the other hand, sequences with balanced gradients (steady-state free-precession, SSFP) are known to have benign off-resonance behavior [25]. As long as B0 inhomogeneities are too small to cause the well-known banding artifacts, which lead to phase jumps of 180°, the SSFP transceive phase fairly resembles the transceive phase of spin echo sequences [25].(S2)
The transceive phase must not contain any contributions from eddy currents in the tissue induced by gradient switching. This can be obtained by averaging two separate measurements with inverted gradient polarization [43]. Alternatively, the balanced gradients of SSFP sequences eliminate this unwanted phase contribution automatically [25].(S3)
Phase contributions from flow and motion should be removed or at least suppressed as much as possible. This job can be done, for example, by double spin echo sequences [44] or, again, by SSFP [25].(S4)
If the “full” complex EPT equation has to be solved and not the phase-based version of EPT, the transceive phase has to be unwrapped before being divided by two, (5). This unwrapping in the three spatial dimensions can be facilitated by performing it separately for each differentiation.Thus, if applicable, SSFP sequences seem to be the sequence of choice for EPT transceive phase measurements. Due to its high efficiency, SSFP is also a good candidate for real-time phase-based EPT [25] as discussed below.


### 3.3. RF Shimming

More and more MR systems atB0=3T and above are equipped with multiple, independent RF transmit channels (see, e.g., [35, 36]). The primary goal of this technique is the compensation of B1 inhomogeneities via RF shimming, that is, the patient-individual adjustment of the complex TX channel weights. It is important to recognize that B1 fields obtained by this RF shimming still satisfy the Helmholtz equations (1) and (2). Thus, if |H+| and ϕ+ of the shimmed system are measured correctly, the accuracy of EPT is not affected. However, the following topics have to be kept in mind regarding EPT in combination with RF shimming. (i)
RF shimming allows RF excitation far from quadrature excitation. Thus, one has to be aware that the transceive phase assumption (5), derived for a QBC with switched polarization, can be significantly violated [45].(ii)
Phase-based EPT assumes|H+| = const, which can be supported by RF shimming [46]. Obviously, RF shimming impairs magnitude-based EPT. Magnitude-based EPT would benefit from RF phase shimming, as long as the resulting B1 inhomogeneities do not cause signal voids, violating assumption (A4).(iii)
Having the complexB1 maps of all TX channels at hand, RF shimming can be performed a posteriori, adjusting the channel weights sequentially for certain ROIs or even voxelwise. Thus, the resulting total B1 can be optimized locally prior to local EPT reconstructions, as was pointed out by [33].


### 3.4. Hybrid Sequences

Since phase-based EPT in principle just requires a standard image’s transceive phase, it is more or less straightforward to combine phase-based EPT with all kinds of other sequences. Below, two possible combinations are depicted. A lot more are expected to come up in the near future. One fruitful combination is obtained by performing EPT and MREIT simultaneously [37]. By injecting external current in the patient for a few ms between RF pulse and signal sampling, MREIT is able to extract tissue conductivity at frequencies below 1 kHz corresponding to the duration of the current injection. Conductivity determined with EPT is related to the Larmor frequency. Thus, the EPT/MREIT hybrid sequence yields a minimal form of conductivity “spectrum.” This frequency dependence can be strikingly proved by wrapping (parts of) phantoms in insulating wrap. The EPT result is not at all affected by the wrapping. For MREIT, the apparent conductivity of the wrapped part is completely suppressed [37]. We refer to [47] that explains how high-frequency current can penetrate the thin membrane and also how this is linked to the frequency-dependent behavior of the complex potential.Another hybrid sequence is described in [48], where EPT is combined with quantitative mapping of the susceptibility χ (QSM). This combination is driven by the idea to get a complete electromagnetic description of the tissue. EPT and QSM are based on two superposed components of the phase of a gradient echo image. QSM requires the phase component related to B0 and off-resonance effects; EPT requires the phase component related to B1, that is, RF transmission and reception as described in Section 3.2. The separation of these two components, for example, by extrapolating the phase of a multiecho sequence to TE = 0 as described above, is the basis of both EPT and QSM. Interested only in EPT or QSM, the noninteresting phase component is just dropped. Thus, the combination of EPT and QSM is simply done by using both phase components as inputs for the respective approaches [48].


## 3.1. MeasuringB1 Magnitude

As seen in (4), the transmit magnitude |H+|, required for EPT, can be measured in a straightforward manner due to its nonlinear impact on the MR signal. In this framework, numerous techniques for |H+| mapping (B1 mapping) are published (see, e.g., [12–17]). In principle, EPT can be based on any B1 mapping method. The accuracy of EPT depends on the accuracy of this mapping; that is, the most accurate |H+| mapping method leads to the most accurate EPT results. Studies looking for the optimum B1 mapping technique, independent of EPT, have been published elsewhere.


## 3.2. MeasuringB1 Phase

As discussed in the previous chapter, determination of theB1 phase ϕ+ always starts with the measurement of the transceive phase. One of the main issues of EPT, as already pointed out in [8], is the contamination of the transceive phase by unwanted phase contributions unrelated to RF penetration. The following steps from (S1) to (S4) describes how to get a transceive phase usable for EPT reconstruction. (S1)
The transceive must not contain any contributions fromB0, that is, any off-resonance effects. The easiest way to exclude off-resonance effects is the use of refocusing pulses, that is, sequences based on spin echoes (SE), like fast spin echo sequences, turbo spin echo sequences, and so on. In contrast, the transceive phase of field-echo based sequences includes off-resonance effects. In this case, these effects can be removed by any kind of B0 mapping. In the easiest way, the phase can be measured at two different TE and extrapolated back to TE = 0 [26]. Also more sophisticated B0 maps can be applied, for example, obtained in the framework of the Dixon techniques (see, e.g., [42]). On the other hand, sequences with balanced gradients (steady-state free-precession, SSFP) are known to have benign off-resonance behavior [25]. As long as B0 inhomogeneities are too small to cause the well-known banding artifacts, which lead to phase jumps of 180°, the SSFP transceive phase fairly resembles the transceive phase of spin echo sequences [25].(S2)
The transceive phase must not contain any contributions from eddy currents in the tissue induced by gradient switching. This can be obtained by averaging two separate measurements with inverted gradient polarization [43]. Alternatively, the balanced gradients of SSFP sequences eliminate this unwanted phase contribution automatically [25].(S3)
Phase contributions from flow and motion should be removed or at least suppressed as much as possible. This job can be done, for example, by double spin echo sequences [44] or, again, by SSFP [25].(S4)
If the “full” complex EPT equation has to be solved and not the phase-based version of EPT, the transceive phase has to be unwrapped before being divided by two, (5). This unwrapping in the three spatial dimensions can be facilitated by performing it separately for each differentiation.Thus, if applicable, SSFP sequences seem to be the sequence of choice for EPT transceive phase measurements. Due to its high efficiency, SSFP is also a good candidate for real-time phase-based EPT [25] as discussed below.


## 3.3. RF Shimming

More and more MR systems atB0=3T and above are equipped with multiple, independent RF transmit channels (see, e.g., [35, 36]). The primary goal of this technique is the compensation of B1 inhomogeneities via RF shimming, that is, the patient-individual adjustment of the complex TX channel weights. It is important to recognize that B1 fields obtained by this RF shimming still satisfy the Helmholtz equations (1) and (2). Thus, if |H+| and ϕ+ of the shimmed system are measured correctly, the accuracy of EPT is not affected. However, the following topics have to be kept in mind regarding EPT in combination with RF shimming. (i)
RF shimming allows RF excitation far from quadrature excitation. Thus, one has to be aware that the transceive phase assumption (5), derived for a QBC with switched polarization, can be significantly violated [45].(ii)
Phase-based EPT assumes|H+| = const, which can be supported by RF shimming [46]. Obviously, RF shimming impairs magnitude-based EPT. Magnitude-based EPT would benefit from RF phase shimming, as long as the resulting B1 inhomogeneities do not cause signal voids, violating assumption (A4).(iii)
Having the complexB1 maps of all TX channels at hand, RF shimming can be performed a posteriori, adjusting the channel weights sequentially for certain ROIs or even voxelwise. Thus, the resulting total B1 can be optimized locally prior to local EPT reconstructions, as was pointed out by [33].


## 3.4. Hybrid Sequences

Since phase-based EPT in principle just requires a standard image’s transceive phase, it is more or less straightforward to combine phase-based EPT with all kinds of other sequences. Below, two possible combinations are depicted. A lot more are expected to come up in the near future. One fruitful combination is obtained by performing EPT and MREIT simultaneously [37]. By injecting external current in the patient for a few ms between RF pulse and signal sampling, MREIT is able to extract tissue conductivity at frequencies below 1 kHz corresponding to the duration of the current injection. Conductivity determined with EPT is related to the Larmor frequency. Thus, the EPT/MREIT hybrid sequence yields a minimal form of conductivity “spectrum.” This frequency dependence can be strikingly proved by wrapping (parts of) phantoms in insulating wrap. The EPT result is not at all affected by the wrapping. For MREIT, the apparent conductivity of the wrapped part is completely suppressed [37]. We refer to [47] that explains how high-frequency current can penetrate the thin membrane and also how this is linked to the frequency-dependent behavior of the complex potential.Another hybrid sequence is described in [48], where EPT is combined with quantitative mapping of the susceptibility χ (QSM). This combination is driven by the idea to get a complete electromagnetic description of the tissue. EPT and QSM are based on two superposed components of the phase of a gradient echo image. QSM requires the phase component related to B0 and off-resonance effects; EPT requires the phase component related to B1, that is, RF transmission and reception as described in Section 3.2. The separation of these two components, for example, by extrapolating the phase of a multiecho sequence to TE = 0 as described above, is the basis of both EPT and QSM. Interested only in EPT or QSM, the noninteresting phase component is just dropped. Thus, the combination of EPT and QSM is simply done by using both phase components as inputs for the respective approaches [48].


## 4. EPT Experiments

### 4.1. Phantom Studies

The principle feasibility of EPT has first been proven using phantoms with different saline concentrations covering roughly the physiologic range, showing a correlation of more than 99% between expected and obtained conductivities [9, 10]. A corresponding permittivity study was published in [49]. For low B0, noise in permittivity images is tremendously higher than in conductivity images, reflecting the low bending of |H+| for B0<3T. For increasing B0, the SNR gain is higher for permittivity than for conductivity, which arises from different powers of ω in (3) [24, 49].The possibility of fast phase-based conductivity imaging has been tested with adding NaCl to a tap water phantom during scanning [25]. The applied SSFP sequence measured images of the whole phantom (size 1 liter, resolution 2.5 × 2.5 × 2.5 mm3) within 4 seconds, visualizing the formation of laminar layers with different salt contents and conductivities at the bottom of the phantom.Conductivity increases with temperature by roughly 2% per Kelvin. This feature has been confirmed with EPT by heating/cooling different phantoms with water and biologic substances (muscle sample, tomato and apple puree) [50]. For the obtained measurement accuracy, no difference in the conductivity’s temperature behavior between the different phantoms has been found. This approach might be helpful in the framework of hyperthermia and thermoablative interventions.


### 4.2. In Vivo Studies: Brain Applications

First EPT in vivo studies have been reported by Voigt et al., investigating different volunteers’ brains [11]. His results reveal several central features of in vivo EPT.(i)
The conductivity and permittivity of healthy grey matter, white matter, and cerebrospinal fluid (CSF) measured with EPT coincide with the literature values [51].(ii)
As expected from theory, phase-based EPT decreases obtained conductivity values by 10%, and magnitude-based EPT increases obtained permittivity values by 10% (see Table1).(iii)
Comparing different volunteers, results vary by roughly 15% for gray and white matter and 20% for CSF. These data comprise measurement uncertainties as well as physiologic intersubject differences.Table 1
Electric properties of different compartments of the brain of a healthy volunteer, measured with EPT at 1.5T [11]. The measured values agree with values expected from the literature [51]. Phase-based EPT yields slightly increased conductivities; magnitude-based EPT yields slightly decreased permittivities.
Conductivity [S/m]
Relative permittivity
Full EPT
Phase-based
Literature
Full EPT
Magnitude-based
Literature
(7)
EPT (9)
[51]
(7)
EPT (10)
[51]
Gray matter
0.69
±
0.14
0.72
±
0.15
0.51
103
±
69
91
±
70
97.4
White matter
0.39
±
0.15
0.43
±
0.15
0.29
72
±
64
63
±
66
67.8
Cerebrospinal fluid
1.75
±
0.34
1.82
±
0.37
2.07
104
±
21
98
±
20
97.3The intersubject variability can be compared with the intra-subject variability, which has been investigated in [52]. Scanning a single volunteer 10 times over a period of three weeks, conductivity variations of roughly 10% across the gray/white matter was obtained. Since this value is close to the intersubject variability reported in [11], physiologic differences between volunteers appear to be rather small. On the other hand, intra-object conductivity variations of phantoms are far below 1%, making system instabilities an unlikely reason for the observed intra-subject variability of 10%. Alternative reasons could be subject motion, spontaneous conductivity fluctuations caused by nutrition or state of health, or others, yet unidentified subject-induced phase instabilities (“unrelated to RF penetration,” see [8]).The conductivity map of a brain of a healthy volunteer is shown in Figure2. EPT is based on the transceive phase of an SSFP scan using a head coil at 1.5 T and a scan time of 3 minutes for an isotropic resolution of 1 × 1 × 1 mm3. Boundary artifacts were reduced with a median filter, locally restricted by the SSFP magnitude image.Figure 2
Volunteer brain study. Left: SSFP image (magnitude); center: SSFP image (phase); right: reconstructed conductivity based on the Laplacian of the SSFP phase and locally restricting the subsequent filtering using the SSFP magnitude.Local inhomogeneity effects using the Helmholtz equation. The reconstructed conductivity and permittivity show strong oscillations. (Simulations were done for a cylinder object of size 100 mm × 100 mm × 120 mm with 3 different values of conductivity (2 S/m, 1 S/m, and 0.5 S/m) and relative permittivity (80, 60, and 40).)
(a)
(b)
(c)
(d)
(e)
(f)For most clinical EPT studies, the focus of interest is on oncology, particularly brain and breast tumors. Initial, single cases of brain tumors have been reported for 1.5 T [53] and 7 T [54]. All cases show a tumor conductivity increased roughly by a factor of two compared with the surrounding white matter, as expected from [55]. In [54], the hypothesis is raised that the tumor conductivity correlates with its sodium content. Systematic studies on a larger number of brain tumor patients, trying to understand the biochemical reasons for enhanced tumor conductivity and to classify different kinds of tumors, are on the way. A single stroke patient was described by [56]. Again, a clear increase of the conductivity within the stroke area has been observed.


### 4.3. In Vivo Studies: Body Applications

#### 4.3.1. Breast

The main components of breast tissue, gland and fat, differ significantly inσ. These components are typically highly nested, leading to high amount of significant conductivity discontinuities throughout the breast. These discontinuities are able to completely spoil any standard EPT reconstruction based on (2) or (9) and (10). Instead of handling conductivity discontinuities on a physical basis via (1), boundary artifacts can be avoided by image segmentation prior to reconstruction and performing separate reconstructions on the different compartments as demonstrated in [39]. Here, the image segmentation was used for shaping the applied kernel locally to the current tissue type, as well as for locally restricting the subsequently applied median filter. The image segmentation was based on the magnitude image of the 3D turbo spin echo performed to obtain the transceive phase for the EPT reconstruction. The resulting conductivity map is more or less free of boundary artifacts. The tumor clearly shows the highest conductivity throughout the breast σ ~2.0–2.5 S/m. Several cysts show an intermediate conductivity σ ~1.0–1.5 S/m. The conductivity of the surrounding fatty tissue is around zero σ ~±0.5 S/m. Remaining inhomogeneities particularly in the fat compartments might arise from insufficient B1 homogeneity or, more likely, from tissue eddy currents, since no second scan with inverted gradient polarization has been performed ([43], step (S2) in the above list). Figure 4 shows breast tumor patient study. Further breast studies are on the way (see, e.g., the breast permittivity study [57]).Figure 4
Breast tumor patient study. Top: turbo spin echo (magnitude); center: reconstructed conductivity image; bottom: contrast enhanced dynamic image. The dynamic image shows the region of the tumor.


#### 4.3.2. Heart

Two isolated, perfused pig hearts where scanned with a gated SSFP sequence [58]. Conductivity values from phase-based EPT of normally perfused heart tissue were compared with the values in ischemic regions after a blockade of the left anterior descending artery. Normal conductivity values turned out to be in agreement with the literature values, while conductivity in ischemic/infarcted areas was 60% lower than that in remote myocardium [58]. These findings are very encouraging for future, challenging in vivo cardiac experiments.


#### 4.3.3. Liver

By analogy to the above-described brain studies, the intra-subject variability of the liver conductivity of healthy volunteers has been investigated [59]. The high efficiency of the applied SSFP sequence allows scanning the whole liver within a single breathhold. The obtained intra-subject and intersubject variability was of the same range as for the brain [11, 52]. A higher artifact level was observed in the liver than in the brain, presumably arising from cardiac motion transferred to the liver and only incompletely suppressed by the applied SSFP sequence. Artifacts are more pronounced for expiration than inspiration breathhold, which can be explained by the contact of heart and liver closer in expiration than in inspiration [59].


#### 4.3.4. Pelvis

In preparation of future pelvis tumor studies, the applicability of phase-based EPT was checked with a pelvis-sized phantom at 3 T [60]. Due to the larger dimensions of the pelvis, as compared to the head, the phase error has to be reinvestigated for this particular anatomy. According to this study, phase-based EPT seems to be sufficient when focusing only on the tumors; however, the inclusion of |H+| could be necessary, when σ of the whole pelvis is required [60].


## 4.1. Phantom Studies

The principle feasibility of EPT has first been proven using phantoms with different saline concentrations covering roughly the physiologic range, showing a correlation of more than 99% between expected and obtained conductivities [9, 10]. A corresponding permittivity study was published in [49]. For low B0, noise in permittivity images is tremendously higher than in conductivity images, reflecting the low bending of |H+| for B0<3T. For increasing B0, the SNR gain is higher for permittivity than for conductivity, which arises from different powers of ω in (3) [24, 49].The possibility of fast phase-based conductivity imaging has been tested with adding NaCl to a tap water phantom during scanning [25]. The applied SSFP sequence measured images of the whole phantom (size 1 liter, resolution 2.5 × 2.5 × 2.5 mm3) within 4 seconds, visualizing the formation of laminar layers with different salt contents and conductivities at the bottom of the phantom.Conductivity increases with temperature by roughly 2% per Kelvin. This feature has been confirmed with EPT by heating/cooling different phantoms with water and biologic substances (muscle sample, tomato and apple puree) [50]. For the obtained measurement accuracy, no difference in the conductivity’s temperature behavior between the different phantoms has been found. This approach might be helpful in the framework of hyperthermia and thermoablative interventions.


## 4.2. In Vivo Studies: Brain Applications

First EPT in vivo studies have been reported by Voigt et al., investigating different volunteers’ brains [11]. His results reveal several central features of in vivo EPT.(i)
The conductivity and permittivity of healthy grey matter, white matter, and cerebrospinal fluid (CSF) measured with EPT coincide with the literature values [51].(ii)
As expected from theory, phase-based EPT decreases obtained conductivity values by 10%, and magnitude-based EPT increases obtained permittivity values by 10% (see Table1).(iii)
Comparing different volunteers, results vary by roughly 15% for gray and white matter and 20% for CSF. These data comprise measurement uncertainties as well as physiologic intersubject differences.Table 1
Electric properties of different compartments of the brain of a healthy volunteer, measured with EPT at 1.5T [11]. The measured values agree with values expected from the literature [51]. Phase-based EPT yields slightly increased conductivities; magnitude-based EPT yields slightly decreased permittivities.
Conductivity [S/m]
Relative permittivity
Full EPT
Phase-based
Literature
Full EPT
Magnitude-based
Literature
(7)
EPT (9)
[51]
(7)
EPT (10)
[51]
Gray matter
0.69
±
0.14
0.72
±
0.15
0.51
103
±
69
91
±
70
97.4
White matter
0.39
±
0.15
0.43
±
0.15
0.29
72
±
64
63
±
66
67.8
Cerebrospinal fluid
1.75
±
0.34
1.82
±
0.37
2.07
104
±
21
98
±
20
97.3The intersubject variability can be compared with the intra-subject variability, which has been investigated in [52]. Scanning a single volunteer 10 times over a period of three weeks, conductivity variations of roughly 10% across the gray/white matter was obtained. Since this value is close to the intersubject variability reported in [11], physiologic differences between volunteers appear to be rather small. On the other hand, intra-object conductivity variations of phantoms are far below 1%, making system instabilities an unlikely reason for the observed intra-subject variability of 10%. Alternative reasons could be subject motion, spontaneous conductivity fluctuations caused by nutrition or state of health, or others, yet unidentified subject-induced phase instabilities (“unrelated to RF penetration,” see [8]).The conductivity map of a brain of a healthy volunteer is shown in Figure2. EPT is based on the transceive phase of an SSFP scan using a head coil at 1.5 T and a scan time of 3 minutes for an isotropic resolution of 1 × 1 × 1 mm3. Boundary artifacts were reduced with a median filter, locally restricted by the SSFP magnitude image.Figure 2
Volunteer brain study. Left: SSFP image (magnitude); center: SSFP image (phase); right: reconstructed conductivity based on the Laplacian of the SSFP phase and locally restricting the subsequent filtering using the SSFP magnitude.Local inhomogeneity effects using the Helmholtz equation. The reconstructed conductivity and permittivity show strong oscillations. (Simulations were done for a cylinder object of size 100 mm × 100 mm × 120 mm with 3 different values of conductivity (2 S/m, 1 S/m, and 0.5 S/m) and relative permittivity (80, 60, and 40).)
(a)
(b)
(c)
(d)
(e)
(f)For most clinical EPT studies, the focus of interest is on oncology, particularly brain and breast tumors. Initial, single cases of brain tumors have been reported for 1.5 T [53] and 7 T [54]. All cases show a tumor conductivity increased roughly by a factor of two compared with the surrounding white matter, as expected from [55]. In [54], the hypothesis is raised that the tumor conductivity correlates with its sodium content. Systematic studies on a larger number of brain tumor patients, trying to understand the biochemical reasons for enhanced tumor conductivity and to classify different kinds of tumors, are on the way. A single stroke patient was described by [56]. Again, a clear increase of the conductivity within the stroke area has been observed.


## 4.3. In Vivo Studies: Body Applications

### 4.3.1. Breast

The main components of breast tissue, gland and fat, differ significantly inσ. These components are typically highly nested, leading to high amount of significant conductivity discontinuities throughout the breast. These discontinuities are able to completely spoil any standard EPT reconstruction based on (2) or (9) and (10). Instead of handling conductivity discontinuities on a physical basis via (1), boundary artifacts can be avoided by image segmentation prior to reconstruction and performing separate reconstructions on the different compartments as demonstrated in [39]. Here, the image segmentation was used for shaping the applied kernel locally to the current tissue type, as well as for locally restricting the subsequently applied median filter. The image segmentation was based on the magnitude image of the 3D turbo spin echo performed to obtain the transceive phase for the EPT reconstruction. The resulting conductivity map is more or less free of boundary artifacts. The tumor clearly shows the highest conductivity throughout the breast σ ~2.0–2.5 S/m. Several cysts show an intermediate conductivity σ ~1.0–1.5 S/m. The conductivity of the surrounding fatty tissue is around zero σ ~±0.5 S/m. Remaining inhomogeneities particularly in the fat compartments might arise from insufficient B1 homogeneity or, more likely, from tissue eddy currents, since no second scan with inverted gradient polarization has been performed ([43], step (S2) in the above list). Figure 4 shows breast tumor patient study. Further breast studies are on the way (see, e.g., the breast permittivity study [57]).Figure 4
Breast tumor patient study. Top: turbo spin echo (magnitude); center: reconstructed conductivity image; bottom: contrast enhanced dynamic image. The dynamic image shows the region of the tumor.


### 4.3.2. Heart

Two isolated, perfused pig hearts where scanned with a gated SSFP sequence [58]. Conductivity values from phase-based EPT of normally perfused heart tissue were compared with the values in ischemic regions after a blockade of the left anterior descending artery. Normal conductivity values turned out to be in agreement with the literature values, while conductivity in ischemic/infarcted areas was 60% lower than that in remote myocardium [58]. These findings are very encouraging for future, challenging in vivo cardiac experiments.


### 4.3.3. Liver

By analogy to the above-described brain studies, the intra-subject variability of the liver conductivity of healthy volunteers has been investigated [59]. The high efficiency of the applied SSFP sequence allows scanning the whole liver within a single breathhold. The obtained intra-subject and intersubject variability was of the same range as for the brain [11, 52]. A higher artifact level was observed in the liver than in the brain, presumably arising from cardiac motion transferred to the liver and only incompletely suppressed by the applied SSFP sequence. Artifacts are more pronounced for expiration than inspiration breathhold, which can be explained by the contact of heart and liver closer in expiration than in inspiration [59].


### 4.3.4. Pelvis

In preparation of future pelvis tumor studies, the applicability of phase-based EPT was checked with a pelvis-sized phantom at 3 T [60]. Due to the larger dimensions of the pelvis, as compared to the head, the phase error has to be reinvestigated for this particular anatomy. According to this study, phase-based EPT seems to be sufficient when focusing only on the tumors; however, the inclusion of |H+| could be necessary, when σ of the whole pelvis is required [60].


## 4.3.1. Breast

The main components of breast tissue, gland and fat, differ significantly inσ. These components are typically highly nested, leading to high amount of significant conductivity discontinuities throughout the breast. These discontinuities are able to completely spoil any standard EPT reconstruction based on (2) or (9) and (10). Instead of handling conductivity discontinuities on a physical basis via (1), boundary artifacts can be avoided by image segmentation prior to reconstruction and performing separate reconstructions on the different compartments as demonstrated in [39]. Here, the image segmentation was used for shaping the applied kernel locally to the current tissue type, as well as for locally restricting the subsequently applied median filter. The image segmentation was based on the magnitude image of the 3D turbo spin echo performed to obtain the transceive phase for the EPT reconstruction. The resulting conductivity map is more or less free of boundary artifacts. The tumor clearly shows the highest conductivity throughout the breast σ ~2.0–2.5 S/m. Several cysts show an intermediate conductivity σ ~1.0–1.5 S/m. The conductivity of the surrounding fatty tissue is around zero σ ~±0.5 S/m. Remaining inhomogeneities particularly in the fat compartments might arise from insufficient B1 homogeneity or, more likely, from tissue eddy currents, since no second scan with inverted gradient polarization has been performed ([43], step (S2) in the above list). Figure 4 shows breast tumor patient study. Further breast studies are on the way (see, e.g., the breast permittivity study [57]).Figure 4
Breast tumor patient study. Top: turbo spin echo (magnitude); center: reconstructed conductivity image; bottom: contrast enhanced dynamic image. The dynamic image shows the region of the tumor.


## 4.3.2. Heart

Two isolated, perfused pig hearts where scanned with a gated SSFP sequence [58]. Conductivity values from phase-based EPT of normally perfused heart tissue were compared with the values in ischemic regions after a blockade of the left anterior descending artery. Normal conductivity values turned out to be in agreement with the literature values, while conductivity in ischemic/infarcted areas was 60% lower than that in remote myocardium [58]. These findings are very encouraging for future, challenging in vivo cardiac experiments.


## 4.3.3. Liver

By analogy to the above-described brain studies, the intra-subject variability of the liver conductivity of healthy volunteers has been investigated [59]. The high efficiency of the applied SSFP sequence allows scanning the whole liver within a single breathhold. The obtained intra-subject and intersubject variability was of the same range as for the brain [11, 52]. A higher artifact level was observed in the liver than in the brain, presumably arising from cardiac motion transferred to the liver and only incompletely suppressed by the applied SSFP sequence. Artifacts are more pronounced for expiration than inspiration breathhold, which can be explained by the contact of heart and liver closer in expiration than in inspiration [59].


## 4.3.4. Pelvis

In preparation of future pelvis tumor studies, the applicability of phase-based EPT was checked with a pelvis-sized phantom at 3 T [60]. Due to the larger dimensions of the pelvis, as compared to the head, the phase error has to be reinvestigated for this particular anatomy. According to this study, phase-based EPT seems to be sufficient when focusing only on the tumors; however, the inclusion of |H+| could be necessary, when σ of the whole pelvis is required [60].


## 5. Conclusion

Tomographic imaging of the conductivity and permittivity distributions inside the human body has been an active research topic in the field of Electrical Impedance Tomography (EIT) measuring boundary current-voltage data [61]. However, wide experience in EIT for more than three decades has shown its methodological limitation in achieving robust reconstructions of static conductivity and permittivity images; EIT is insensitive to any perturbation of internal admittivity whereas it is very sensitive to forward modeling errors. It seems that boundary measurements are insufficient for robust reconstruction of the admittivity distribution inside the subject. For a robust reconstruction, we need an internal measurement that is capable of by MRI scanners as discussed in this paper.Using standard MR systems and standard MR sequences, mapping of the electric properties seems to be clinically feasible, particularly phase-based conductivity imaging. The rapidly evolving field will certainly afford further improved measurement and reconstruction techniques in the near future. The broad spectrum of started preclinical and clinical studies raise hope that answers will soon be available concerning potential diagnostic benefits of EPT.



---

*Source: 546562-2013-03-07.xml*