Recent Progress and Future Challenges in MR Electric Properties Tomography

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 using B 1 maps. The paper discusses the development of EPT reconstructions, EPT sequences, EPT experiments, and challenging issues of EPT.


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 := − ( / ) at an angular frequency (assumed to be below microwave range) can be probed by the time-harmonic magnetic field H = ( , , ) through the following arrangement of time Maxwell equations (so-called Helmholtz equation): where is the magnetic permeability and is assumed to be isotropic. Here, the corresponding time-varying field is Re{H }. 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 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 algorithm, was developed to provide both conductivity image and current density image. After invention of the harmonic 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.
(A1) A locally ("piecewise") constant (r), that is; ∇ (r) = 0. This assumption has severe consequences, which has to be discussed extensively later on.  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 field 0 ) 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.
(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) | + | 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 of + in the numerator and denominator of the expression (2). Given this cancelation, (2) yields absolute values of , even for arbitrary scaling of + . 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 3 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 + = | + | + . Unfortunately, only the magnitude | + | enters the MR signal in a nonlinear way; for each nominal flip angle of the sequence, the following MR signal is measured: where 0 is the MR magnitude image containing relaxation and spin density effects, − = ( − )/2, and 1 , 2 system-dependent constants [21]. Thus, | + | has the chance to be measured exactly, assuming the ideal working of corresponding 1 mapping techniques (see, e.g., [12][13][14][15][16][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 − 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 phasẽ0 = + +̃− of this setup sometimes 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 | + |, 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: where r is an arbitrarily oriented area centered at r with its boundary r , l the line element, and S 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 − and are not measurable directly, − = 0 and = 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 + . For minimal kernel size, the SNR of the reconstructed is far below the SNR of the input + [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: where r is a volume centered at r with its boundary r and 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 4 Computational and Mathematical Methods in Medicine in [23]. Equation (7) can be rewritten with separate real and imaginary part As already indicated in [9], observation (O2), the assumption |∇ 2 + | ≫ 2|∇ ln | + | ⋅ ∇ + | yields phase-based EPT for conductivity imaging and the assumption |∇ 2 | + |/| + || ≫ |∇ + | 2 yields magnitude-based EPT for permittivity imaging 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 since (9) can be based on + , − , or any phase fulfilling Maxwell's equations, leading to the same as long as ∇| + | = 0 and ∇| − | = 0 are fulfilled.
(F2) Skipping the need of mapping (the magnitude of) 1 , 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 0 = 1.5 T, phase-based EPT has been confirmed at 0 = 7T [26]. The impact of 0 on EPT and the related question of 0 optimal for EPT turned out to be a nontrivial task [27]. Obviously, higher SNR can be expected with increasing 0 . This advantage is counterbalanced by the increasing violation of the assumption ∇| + | = 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 0 = 3 T for conductivity imaging. For permittivity imaging, the violation of the assumption ∇ + = 0 does not increase with 0 , and the highest available 0 seems to be optimal. This trend is further emphasized by the different powers of in (3).
The problem of separating + and − from the transceive phase has been solved analytically by an approach sometimes called "Local Maxwell Tomography" (LMT) [32][33][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 0 , 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 and can be compared based on different TX channels and , but same receive channel with phase − , yielding [34] (r) − (r) = 1 This allows, first, first, the determination of the unknown − and, subsequently, the straightforward determination of the unknown + and + . The central idea of comparing reconstruction results of two or more different TX channels can also be utilized to exactly distinguish | − | from the spin magnetization 0 , 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 | − | 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 Computational and Mathematical Methods in Medicine 5 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] (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 − = 0 and = 0 as it was done in connection with (6) [10]. A challenging problem is to find an exact relation between and + from the full equation (1) in such a way that can be computed robustly and efficiently using + 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 v 1 , v 2 , and v 3 (unit vectors) and its corresponding eigenvalues 1 , 2 , and 3 , respectively: 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 in (6) reflects the degree of anisotropy [40,41]. To be precise, let (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 2 [41]: A perfectly isotropic should not depend on the direction n of (n). On the other hand, the reconstructed should show a minimum for (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 − .

EPT Sequence
3.1. Measuring 1 Magnitude. As seen in (4), the transmit magnitude | + |, 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 | + | mapping ( 1 mapping) are published (see, e.g., [12][13][14][15][16][17]). In principle, EPT can be based on any 1 mapping method. The accuracy of EPT depends on the accuracy of this mapping; that is, the most accurate | + | mapping method leads to the most accurate EPT results. Studies looking for the optimum 1 mapping technique, independent of EPT, have been published elsewhere.

Measuring 1 Phase.
As discussed in the previous chapter, determination of the 1 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 from 0 , 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 0 mapping. In the easiest way, the phase can be measured at two different TE and extrapolated back to TE = 0 [26]. Also more sophisticated 0 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 0 inhomogeneities are too small to cause the wellknown 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.

RF Shimming.
More and more MR systems at 0 = 3 T 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 1 inhomogeneities via RF shimming, that is, the patient-individual adjustment of the complex TX channel weights. It is important to recognize that 1 fields obtained by this RF shimming still satisfy the Helmholtz equations (1) and (2). Thus, if | + | 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].
(iii) Having the complex 1 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 1 can be optimized locally prior to local EPT reconstructions, as was pointed out by [33].

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 0 and off-resonance effects; EPT requires the phase component related to 1 , 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].

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 0 , noise in permittivity images is tremendously higher than in conductivity images, reflecting the low bending of | + | for 0 < 3T. For increasing 0 , 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 mm 3 ) 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 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. (ii) As expected from theory, phase-based EPT decreases obtained conductivity values by 10%, and magnitudebased EPT increases obtained permittivity values by 10% (see Table 1).
(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.
The 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 Figure 2. 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 mm 3 . Boundary artifacts were reduced with a median filter, locally restricted by the SSFP magnitude image.
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.

In Vivo Studies: Body Applications
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 1 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]).

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.

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 intrasubject 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].

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, Computational and Mathematical Methods in Medicine the inclusion of | + | could be necessary, when of the whole pelvis is required [60].

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.