- Research
- Open Access
- Published:

# Reconstruction of dual-frequency conductivity by optimization of phase map in MREIT and MREPT

*BioMedical Engineering OnLine*
**volume 13**, Article number: 24 (2014)

## Abstract

### Background

The spectroscopic conductivity distribution of tissue can help to explain physiological and pathological status. Dual frequency conductivity imaging by combining Magnetic Resonance Electrical Property Tomography (MREPT) and Magnetic Resonance Electrical Impedance Tomography (MREIT) has been recently proposed. MREIT can provide internal conductivity distributions at low frequency (below 1 kHz) induced by an external injecting current. While MREPT can provide conductivity at the Larmor frequency related to the strength of the magnetic field. Despite this potential to describe the membrane properties using spectral information, MREPT and MREIT techniques currently suffer from weak signals and noise amplification as they both reply on differentiation of measured phase data.

### Methods

We proposed a method to optimize the measured phase signal by finding weighting factors according to the echo signal for MREPT and MREIT using the ICNE (Injected current nonlinear encoding) multi-echo pulse sequence. Our target weights are chosen to minimize the measured noise. The noise standard deviations were precisely analyzed for the optimally weighted magnetic flux density and the phase term of the positive-rotating magnetic field. To enhance the quality of dual-frequency conductivity images, we applied the denoising method based on the reaction-diffusion equation with the estimated noise standard deviations. A real experiment was performed with a hollow cylindrical object made of thin insulating film with holes to control the apparent conductivity using ion mobility and an agarose gel cylinder wrapped in an insulating film without holes to show different spectroscopic conductivities.

### Results

The ability to image different conductivity characteristics in MREPT and MREIT from a single MR scan was shown by including the two objects with different spectroscopic conductivities. Using the six echo signals, we computed the optimized weighting factors for each echo. The qualities of conductivity images for MREPT and MREIT were improved by optimization of the phase map. The proposed method effectively reduced the random noise artifacts for both MREIT and MREPT.

### Conclusion

We enhanced the dual conductivity images using the optimally weighted magnetic flux density and the phase term of positive-rotating magnetic field based on the analysis of the noise standard deviations and applying the optimization and denoising methods.

## Background

The conductivity spectra of biological tissues can provide diagnostic medical information from the estimation of physiological and pathological conditions of *in-vivo* and *ex-vivo* tissue. However it is difficult to produce high resolution conductivity images inside the human body [1, 2]. The conventional conductivity imaging methods have limited spatial resolution and sensitivity inherited from the ill-posed nature of the problem [3]. In order to achieve sensitive conductivity images with high resolution, electric impedance imaging techniques based on a magnetic resonance imaging (MRI) have been propsed. These include magnetic resonance electrical impedance tomography (MREIT) and magnetic resonance electrical property tomography (MREPT) which are under active investigation [4–17]. Both methods use the internal magnetic field information obtained from the phase data of MRI scanner to reconstruct the internal conductivity image.

MREPT can provide the electrical conductivity information at the Larmor frequency by measuring the phase of positive rotating field due to the applied *B*
_{1} field. It does not require application of external current and therefore directly recovers the conductivity distribution by taking the derivative of the measured phase signals twice. MREIT needs a pair of electrodes to inject current into an imaging object during MRI scan for measuring a magnetic flux density induced by the external injecting current. The MREIT technique used only the *z*-component of magnetic flux density, *B*
_{
z
}, of **B**=(*B*
_{
x
},*B*
_{
y
},*B*
_{
z
}) to reconstruct the cross-sectional apparent conductivity image at a lower frequency range (below 1 kHz) [18–22]. The phase difference approach with an interleaved encoding scheme was adopted to cancel the systematic artifacts accumulated in phase signals and also reduce the random noise artifacts. Recently, a simultaneous conductivity imaging technique using a combination of MREPT and MREIT was proposed to provide the dual-frequency conductivities of tissue from a common MR scan [23]. MREPT using *B*
_{1}-mapping technique visualizes the conductivity and permittivity distributions at the Larmor frequency and MREIT recovers the apparent conductivity distribution when injecting low-frequency external current through the attached electrodes. Since the biological tissues show the frequency dependent conductivity property [1, 2], the simultaneous dual-frequency conductivity imaging using a single MR scan is beneficial to provide distinct electrical features of tissues quantitatively.

Both conductivity imaging techniques use the phase signals of measured MR data. MREIT and MREPT commonly suffer from weak signals and noise amplification by using the derivative of the measured signal. The noise level of *B*
_{
z
} in MREIT is inversely proportional to the signal-to-noise ratio (*Υ*^{j}) of the MR magnitude image and the current injection pulse width. Because of the small amount of injection current and the poor quality of measured *B*
_{
z
} change due to the injected current, it is difficult to perform *in vivo* human experiments using a conventional MR pulse sequence. To enhance the magnetic flux density due to the injected current in MREIT, the injected current nonlinear encoding (ICNE) method was introduced. This extended the duration of injecting current until the end of a readout gradient, and improved the signal by using a multi-echo train MR pulse sequence [24–26]. MREPT also suffers from low sensitivity due to the inherently poor signal to noise ratio and noise sensitive characteristics as it also needs the derivative of measured data, is very sensitive to the measured noise.

In this paper, we adopt a multiple spin echo MREIT pulse sequence based on the ICNE scheme to measure multiple phase data for MREIT and MREPT images. The acquired multiple phase data can be decomposed to the phase term reflecting the magnetic flux density signal induced by the injected current and the other phase term of positive rotating field due to the applied *B*
_{1} field. We analyze the noise level of the two decomposed phase terms and minimize the measured random noise artifacts by applying the optimal combination of multiple phase terms. Also, we apply a denoising technique to the optimized magnetic flux density for MREIT and to the phase signal for MREPT in order to improve the quality of the reconstructed conductivity images. We prepared a conductivity phantom consisted of two different kinds of anomalies to show the difference of MREIT and MREPT. A phantom experiment is conducted to validate that the proposed method is able to improve the qualities of reconstructed dual-frequency conductivity images compared to the results of using conventional MREIT and MREPT reconstruction algorithms.

## Methods

### Governing equation

We denote an imaging object as *Ω*. The admittivity in *Ω* is *κ*=*σ*(*ω*)+*i* *ω* *ε*(*ω*), where *σ*(*ω*) and *ε*(*ω*) are the conductivity and permittivity, respectively, at the angular frequency *ω*. Time-harmonic Maxwell’s equations relate the electric field **E** and magnetic flux field **H**:

where the current **J** and the electric field intensity **E** satisfy the relation **J**=*κ* **E** by Ohm’s law.

In MREIT, we can assume that the injecting current is sufficiently low frequency to meet *i* *ω* *ε*(*ω*)≈0 and neglect *i* *ω* *μ*
_{0}≈0, Maxwell’s equation (1) satisfies

where *σ*
_{
L
}=*σ*(*ω*) denotes the low-frequency electrical conductivity responding to the externally injecting current. Since the electric field intensity **E**=−∇*u* is a gradient form under negligible angular frequency *ω*, we have

where **n** is the outward normal vector on the surface *∂* *Ω* and *g* is the current density on the surface.

In MREPT, by assuming the locally homogeneity, the magnetic field **H** can be expressed as a simple form:

The phase terms retrieved from the measured phase are the combination of the phase terms of positive {H}^{+}=\frac{{H}_{x}+{\mathit{\text{iH}}}_{y}}{2} and negative {H}^{-}=\frac{{H}_{x}-{\mathit{\text{iH}}}_{y}}{2} rotating field. For a restricted situation such as the usage of transmit-receive coil, the phase term of positive rotating field can be retrieved from the measured phase. The transverse field of **H** can be decomposed into the positively rotating field {H}^{+}=\frac{{H}_{x}+{\mathit{\text{iH}}}_{y}}{2} and the negatively rotating field {H}^{-}=\frac{{H}_{x}-{\mathit{\text{iH}}}_{y}}{2}. We denote *φ*^{+} and *φ*^{−} as the phase term of *H*^{+} and *H*^{−}, respectively. When single channel transmit-receive coils are used, the transmit *φ*^{+} can be reasonably similar to *φ*^{−}
[16]. The measurable positively rotating field {H}^{+}=\frac{{H}_{x}+{\mathit{\text{iH}}}_{y}}{2} of **H**=(*H*
_{
x
},*H*
_{
y
},*H*
_{
z
}) also satisfies the relation under the locally homogeneity condition

By separating the real and imaginary parts of *κ* and assuming *φ*^{+}≈*φ*^{−}, the phase term of *H*^{+} can estimate the conductivity *σ*
_{
H
}=*σ*(*ω*) which reflects the high-frequency property responding to the RF pulse [16]:

### ICNE multi-echo pulse sequence and multiple phase data

The ICNE multi-echo pulse sequence is adopted to acquire multiple echoes measured at the echo time {T}_{{E}_{j}},\phantom{\rule{1em}{0ex}}j=1,\cdots \phantom{\rule{0.3em}{0ex}},{N}_{E}, where *N*
_{
E
} is the number of echoes per repetition time *T*
_{
R
}. Figure 1 shows its schematic diagram where *I*^{+} and *I*^{−} are sequentially injected positive and negative currents synchronized with the applied 180° rephasing pulses with *N*
_{
E
} times. The *j*-th complex image, {\zeta}_{j}^{\pm}, corresponding to the *j*-th readout gradient generated with the positive and negative injecting currents *I*^{+} and *I*^{−}, can be expressed as

where {M}_{\mathit{\text{xy}},j}^{\alpha}\left(\mathbf{r}\right)={M}_{\mathit{\text{xy}}}^{\alpha}\left(\mathbf{r}\right){e}^{-{T}_{{E}_{j}}/{T}_{2}\left(\mathbf{r}\right)} is the *j*-th transverse magnetization at a flip angle *α*, *δ*
_{
ε
} is the systematic phase artifact, *γ* is the gyromagnetic ratio of hydrogen, *T*
_{2} is the transverse relaxation time, and {T}_{{c}^{\phantom{\rule{0.3em}{0ex}}j}} is the *j*-th duration of the injecting current.

The *j*-th transverse magnetization {M}_{\mathit{\text{xy}},j}^{\alpha}\left(\mathbf{r}\right) can be rewritten in detail as

where *C*
_{1} and *C*
_{2} are system-dependent constants. Assuming *δ*
_{
ε
}≈0, the phase signal can be separated as

By subtracting and adding the phase signals {\mathcal{P}}_{j}^{\pm}, we can obtain the magnetic flux density *B*
_{
z
} induced by the externally injecting current at the low-frequency and the phase term of *H*^{+} at the Larmor frequency simultaneously.

From the relation 10, we can compute the *j*-th phase signals for MREPT and the *j*-th *z*-component of magnetic flux density for MREIT:

The noise standard deviation of the measured magnetic flux density {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} in 11 is inversely proportional to the current injection time {T}_{c}^{\phantom{\rule{0.3em}{0ex}}j} and the signal-to-noise ratio (*Υ*^{j}) of MR magnitude image following [27, 28]

A similar analysis provides the noise standard deviation of *φ* *j*+ as

### Optimal combination of multiple {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}

In order to enhance the conductivity image at the low frequency using MREIT, we have to improve the quality of the *B*
_{
z
} signal. An alternative to reduce the noise of *B*
_{
z
} is to find an optimal combination of multiple {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}},\phantom{\rule{1em}{0ex}}j=1,\cdots \phantom{\rule{0.3em}{0ex}},{N}_{E}, by determining the weighting factors which satisfied conditions of *ξ*
_{
j
}(**r**)>0 and \sum _{j=1}^{{N}_{E}}{\xi}_{j}\left(\mathbf{r}\right)=1. When we assume that the measured *k*-space signals are contaminated with independent identically distributed (IID) complex Gaussian random noise as same amount of expected noise standard deviation of each {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} in 12, the noise variance of the weighted {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} satisfies the following relation

To determine the optimized weighting factor *ξ*
_{
j
}(**r**), the noise variance of the combined {B}_{z}^{\xi} is required to minimize the following conditions:

where {\mathcal{F}}_{{B}_{z}}({\xi}_{1},\cdots \phantom{\rule{0.3em}{0ex}},{\xi}_{{N}_{E}}):=\sum _{j=1}^{{N}_{E}}\frac{{\xi}_{j}^{2}\left(\mathbf{r}\right)}{{T}_{{c}^{\phantom{\rule{0.3em}{0ex}}j}}^{2}{\left|{\zeta}_{j}^{\pm}\right(\mathbf{r}\left)\right|}^{2}} and {\mathcal{G}}_{{B}_{z}}({\xi}_{1},\cdots \phantom{\rule{0.3em}{0ex}},{\xi}_{{N}_{E}}):=\sum _{j=1}^{{N}_{E}}{\xi}_{j}\left(\mathbf{r}\right). The method of Lagrange multipliers to solve the optimization problem in 15 introduces a new variable *λ* and requires the condition that the gradients of {\mathcal{F}}_{{B}_{z}} and {\mathcal{G}}_{{B}_{z}} are parallel. Using the method of Lagrange multipliers, the optimal weighting factor *ξ*
_{
j
}(**r**) can be obtained as

where {\Psi}_{j}\left(\mathbf{r}\right):={T}_{{c}^{\phantom{\rule{0.3em}{0ex}}j}}^{2}\left|{\zeta}_{j}^{\pm}\right(\mathbf{r}){|}^{2}.

The optimized magnetic flux density to produce a conductivity image in MREIT can be generated by the weighted average of multiple {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} with the determined weighting factor of *ξ*
_{
j
}(**r**) as {B}_{z}^{\xi}\left(\mathbf{r}\right)=\sum _{j=1}^{{N}_{E}}{\xi}_{j}\left(\mathbf{r}\right){B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}\left(\mathbf{r}\right) using the ICNE multi-echo pulse sequence. The noise standard deviation of {B}_{z}^{\xi} can be precisely described by substituting 16 into the noise variance of {B}_{z}^{\xi} in 14:

where {\mathcal{N}}_{{M}_{\mathit{\text{xy}}}} denotes the noise level of magnitude image.

### Optimal combination of multiple *φ* *j*+

The noise variance of the weighted *φ* *j*+ satisfies the following relation similar as 14.

The noise variance of the combined *φ*^{+} with weighting factors is required to minimize the following conditions,

where {\mathcal{F}}_{{\phi}^{+}}({\chi}_{1},\cdots \phantom{\rule{0.3em}{0ex}},{\chi}_{{N}_{E}}):=\sum _{j=1}^{{N}_{E}}\frac{{\chi}_{j}^{2}\left(\mathbf{r}\right)}{\left|{\zeta}_{j}^{\pm}\right(\mathbf{r}){|}^{2}} and {\mathcal{G}}_{{\phi}^{+}}({\chi}_{1},\cdots \phantom{\rule{0.3em}{0ex}},{\chi}_{{N}_{E}}):=\sum _{j=1}^{{N}_{E}}{\chi}_{j}\left(\mathbf{r}\right). The optimal weighting factor, *χ*
_{
j
}(**r**), corresponding to each echo phase can be determined by using the method of the Lagrange multipliers as

where {\Phi}_{j}\left(\mathbf{r}\right):=\left|{\zeta}_{j}^{\pm}\right(\mathbf{r}){|}^{2}.

The optimally weighted average of multiple phase signal *φ* *j*+ with the determined weighting factor of *χ*
_{
j
}(**r**), {\phi}^{+,\chi}\left(\mathbf{r}\right)=\sum _{j=1}^{{N}_{E}}{\chi}_{j}\left(\mathbf{r}\right){\phi}_{j}^{+}\left(\mathbf{r}\right), can produce an enhanced conductivity image in MREPT. The noise standard deviation of the optimized phase signal, *φ*^{+,χ}, is computed by substituting (20) into the noise variance of *φ*^{+,χ} in (18) to estimate pixel-by-pixel noise precisely:

Considering only random noise effects, the optimal weighting factor *χ*
_{
j
}(**r**) reduces the noise level, depending on the echo number and *T*
_{2}-decay rate at each imaging pixel.

### Denoising method using the estimated {\mathit{\text{sd}}}_{{B}_{z}^{\xi}} and {\mathit{\text{sd}}}_{{\phi}^{+,\chi}}

A denoising technique is applied to gradients of the optimized magnetic flux density for MREIT image and the optimized phase signal for MREPT image since the reconstruction procedures only require the differentiated measured data. Also, it is advantageous to denoise \nabla {B}_{z}^{\xi} and ∇*φ*^{+,χ} because the measured magnetic flux density of {B}_{z}^{\xi} and the phase of *H*^{+} are inherently continuous without conventional edge information.

To denoise the optimized \nabla {B}_{z}^{\xi} and ∇*φ*^{+,χ}, we use the following conventional reaction-diffusion equation [29, 30]:

where the parameter *β* is a fidelity term and the function *ψ* satisfies

Here, the initial state *f* denotes each component of the optimized magnetic flux density \nabla {B}_{z}^{\xi} or the phase signal ∇*φ*^{+,χ}, which is to be denoised. The time dependent solution *v*(**r**,*t*) is the denoised image of the initial state image *f*.

To remove the random noise artifact while preserving the edge information, it is important to determine the diffusion function *ψ* and *β*. Using the estimated noise level of {\mathit{\text{sd}}}_{{\phi}^{+,\chi}}, the diffusion function *ψ* and the fidelity function *β* can be specified as

where *ς*
_{
ε
} is a small parameter to guarantee the positive sign of |∇*v*(**r**,*t*)|+*ς*
_{
ε
}. For the magnetic flux density {B}_{z}^{\xi}, the diffusion function *ψ* can be selected similarly.

The quality of magnetic flux density signal in MREIT depends on the magnitude intensity and width of injected current simultaneously, while the quality of phase signal in MREPT only depends on the magnitude intensity. Therefore, there are common effects due to the decay of the magnitude intensity in MREPT and MREIT processing. There are also different effects caused by externally injecting current in MREIT. These occur when we calculate each of the optimal weighting factors and the noise standard deviation in the optimal combination of multiple phase and magnetic flux density images. The proposed optimization method uses the decay rate of magnitude intensities to determine the weighting factors in MREIT and MREPT, which relate with the *T*
_{2} values. However, the determined weighting factors in 16 and 20 only include the measured magnitude intensity at each echo time {T}_{{E}_{j}},\phantom{\rule{0.3em}{0ex}}j=1,\cdots \phantom{\rule{0.3em}{0ex}},{N}_{E}. The proposed method therefore does not need estimates *T*
_{2} values to optimize the multiple phase signals.

### Phantom design and experimental setup

In order to evaluate the proposed optimizing method using ICNE multi-echo MR pulse sequence, we designed a phantom which can show different characteristics of reconstructed conductivity images by MREIT and MREPT. Figure 2(a) and (b) show the phantom configuration at the middle slice of the object and the photo view of the phantom, respectively. We built a cylindrical phantom filled with a saline of 0.2 *Sm*^{−1} conductivity (0.3 g/L NaCl and 1 g/L CuSO_{4}). The diameter and height of the phantom were 11 cm and 14 cm, respectively. We attached four carbon hydrogel electrodes (HUREV Co. Ltd, Korea) on the side of the acrylic container for injecting current.

Two different objects were positioned inside phantom. The left one was a thin hollow cylindrical object with the diameter of 4 cm using an insulating thin film of 0.4 mm thickness. We punched four holes which had 2 mm diameter with equally spaced around the circumference. Because we filled the same saline of 0.2 *Sm*^{−1} inside and outside of the hollow cylindrical object with holes, the apparent electrical conductivities inside and outside of object were determined by the movement of ions through the holes [31]. The right cylindrical object with the diameter of 3 cm was made of an agarose gel (1 g/L CuSO_{4}, 2.1 g/L NaCl, 15 g/L Agar) to generate different conductivity, spin density and *T*
_{2}-decay compared with the background saline. The conductivity of agarose gel was 1.10 *Sm*^{−1}. It was wrapped in a thin insulating film without holes.

A transversal injecting current of 10 mA was introduced into the phantom *via* a pair of recessed carbon hydrogel electrodes attached at the middle of the phantom. Figure 2(c) shows the magnitude image at the middle slice. We used the ICNE multi-echo MR pulse sequence by injecting current of alternating polarity synchronized with the multiple refocusing pulses as shown in Figure 1.

Imaging parameters in a 3T MRI scanner (Achieva, Philips) with birdcage transmit-receive (Tx/Rx) RF head coil were as follows: repetition time *T*
_{
R
}=1200 ms, data acquisition time width *T*
_{
s
}=3.584 ms, echo-spacing {T}_{{E}_{\mathit{\text{sp}}}}=15 ms, total number of echo *N*
_{
E
}=6 and number of averaging = 4. The reconstructed image matrix was 128×128, with a FOV of 180×180 mm^{2}, 5 mm slice thickness. Total scanning time for both vertical and horizontal injection currents was 20 minutes with an interleaved phase encoding acquisition.

## Results

Figure 3 shows the multiple MR magnitude images measured at the echo time {T}_{{E}_{j}}=15\times j ms for *j*=1,⋯,6, where the right anomaly shows rapid decay of the transverse magnetization due to magnetic field non-uniformity and spin-spin transverse relaxation.

Figure 4 shows the phase images of *φ* *j*+ of *H*^{+} map and magnetic flux density {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} induced by the vertically injected current using ICNE multi-echo pulse sequence at the echo time {T}_{{E}_{j}}=15\times j ms, *j*=1,⋯,6. The images in Figure 4(a) and (b) were estimated by adding and subtracting the phase signals of {\zeta}_{j}^{\pm}, respectively. Since the electromagnetic wave at the Larmor frequency of 128 MHz could penetrate the thin insulating film and there was same saline inside and outside of the thin insulating wall of the left hollow cylinder, the measured phase values of *φ* *j*+ for the left object in MREPT could not produce conductivity difference between inside and outside. However, the measured {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}} images at the low frequency in MREIT showed the different value inside the left anomaly because MREIT measures the apparent magnetic flux density corresponding to the ion mobility and intrinsic conductivities of composite materials. Note that the slope of magnetic flux density was changed abruptly due to the high current density around the holes in the same direction of the current injection. Regarding the different conductivity agarose gel cylinder wrapped in the thin insulating film without holes, both images in Figure 4(a) and (b) could present distinctions between inside of the object and the background saline.

Using the proposed optimization method with the estimated noise standard deviations {\mathit{\text{sd}}}_{{\phi}_{j}^{+}} and {\mathit{\text{sd}}}_{{B}_{z}^{\phantom{\rule{0.3em}{0ex}}j}},\phantom{\rule{0.3em}{0ex}}j=1,\cdots \phantom{\rule{0.3em}{0ex}},6, we found the optimal weighting factors for the phase of *H*^{+} and the magnetic flux density by the injected current, respectively. Figure 5(a) and (b) show the determined weighting factors for dual-frequency conductivity images in MREPT and MREIT combined method, respectively. In MREPT, the intensity of weighting factors was monotonically decreased because the noise level of *φ* *j*+ only depended on *T*
_{2}-decay rate of {\zeta}_{j}^{\pm},\phantom{\rule{1em}{0ex}}j=1,\cdots \phantom{\rule{0.3em}{0ex}},6, whereas the weighting factors in MREIT show the different characteristics because they depended on the combination of *T*
_{2}-decay rate and the duration of injected current simultaneously. In Figure 5(b), the intensity of weighting factors in the background area was monotonically increased because the effect of time increasing for the injection current was greater than the magnitude attenuation effect of echo signals. However, the weighting factors in the right agarose gel anomaly region were different due to the combination of both effects. The weighting factor values at the second echo were higher than those of other echoes.

Figure 6 shows the reconstructed dual-frequency conductivity images using the multiple echoes *φ* *j*+ and {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}^{d} for *j*=1,⋯,6 and *d*=1,2, where {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}^{1} and {B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}^{2} were the measured magnetic flux density data by the vertically and horizontally injected currents, respectively. To reconstruct the high-frequency conductivity distribution in Figure 6(a), we used phase for each echo signal as

Here, we could not detect the hollow cylinder in the left side made of the thin insulating film with holes. The high-frequency conductivity inside and outside of the left object was estimated to be the same value. The quality of reconstructed conductivity image for each echo in Figure 6(a) was monotonically degraded corresponding to the weighting factors shown in Figure 5(a). Table 1 shows the values of weighting factor inside and outside of the right agarose object wrapped by insulating film.

We estimated the projected current density **J**^{P}, which was the optimal divergence-free component from the measured *z*-component magnetic flux density in order to reconstruct the low-frequency conductivity images [32]. The projected current density consists of the background current density **J**^{0}=−∇*α* and{\nabla}_{\mathit{\text{xy}}}^{\perp}\eta =\left(\frac{\mathrm{\partial \eta}}{\mathrm{\partial y}},-\frac{\mathrm{\partial \eta}}{\mathrm{\partial x}},0\right), where the potential *η* satisfies at the imaging slice *Ω*
_{
t
}⊂*Ω*

where ∇_{
xy
}denotes the two-dimensional gradient. The absolute conductivity distribution can be written as a matrix form using the two measured{B}_{z}^{d},\phantom{\rule{1em}{0ex}}d=1,2,data

where *Ω*
_{
t
}is an imaging slice, *τ*= log*σ*, and

Figure 6(b) shows the reconstructed low-frequency conductivity images by solving 27 using the *j*-th measured two magnetic flux densities{B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}^{d},\phantom{\rule{1em}{0ex}}d=1,2obtained with the vertically and horizontally injected currents. The reconstructed conductivity values inside the hollow cylinder with holes in the left side reflected the externally injected current flowed through the four holes. Regarding the agarose gel object wrapped in the insulating film without holes, the apparent conductivity reconstructed by solving 27 shows the same anomaly as an insulator since the measured *B*
_{
z
}data only reflected the current density distribution due to the externally injected current through the attached electrodes. The quality of conductivity image using the sixth echo data was poor inside the agarose gel anomaly region due to relatively short *T*
_{2}-decay relaxation time. It was consistent with the corresponding weighting factor in Figure 5(b). Table 2 shows the values of weighting factor inside and outside of the right agarose object wrapped by insulating film.

Figure 7(a) and (e) show the normalized noise standard deviations of optimally weighted *φ*^{+,χ}and{B}_{z}^{\xi}, respectively. In previous approaches without using multi-echo data with optimizing method, the conductivity images for MREPT and MREIT were reconstructed as shown in Figure 7(b) and (f). Due to the different *T*
_{2}-decay rate and the duration of injected current, the intensities of{\mathit{\text{sd}}}_{{\phi}^{+,\chi}}and{\mathit{\text{sd}}}_{{B}_{z}^{\xi}}represented different characteristics. The reconstructed dual-frequency conductivity images with the optimally weighted *φ*^{+,χ}and{B}_{z}^{d,\xi},\phantom{\rule{1em}{0ex}}d=1,2,were displayed in Figure 7(c) and (g). Solving the reaction-diffusion equation 22 with the estimated noise standard deviations{\mathit{\text{sd}}}_{{\phi}^{+,\chi}}and{\mathit{\text{sd}}}_{{B}_{z}^{\xi}}, we denoised ∇*φ*^{+,χ}and\nabla {B}_{z}^{d,\xi}; the iteration number was 400 and the regularization function of *β* was\frac{0.05}{\mathcal{N}\left({\mathit{\text{sd}}}_{{\phi}^{+,\chi}}\right)}for MREPT case, where\mathcal{N}\left({\mathit{\text{sd}}}_{{\phi}^{+,\chi}}\right)was the normalized noise standard deviation in Figure 7(a). The denoising parameters for MREIT were the same as those for MREPT. Figure 7(d) and (h) show the reconstructed dual-frequency conductivity images after denoising of ∇*φ*^{+,χ}and\nabla {B}_{z}^{d,\xi}, respectively.

## Discussion

MREPT and MREIT used an internal magnetic flux density information from a phase imaging can provide conductivity images with high resolution compared to the conventional electrical impedance imaging techniques. Moreover, the interleaved measurement technique with alternating injected current polarity based on a multi-spin-echo pulse sequence, enables simultaneous dual conductivity images. From a single MR scan, we could obtain the low-frequency conductivity images by subtracting and the high-frequency conductivity images by adding the measured phase terms. Despite their advantages and usefulness, MREPT and MREIT techniques commonly suffer from noise amplifications in phase imaging. The proposed method used a ICNE multi-echo pulse sequence based on a spin echo pulse sequence in order to suppress the background field inhomogeneity. This also influenced to the reconstructed high-frequency conductivity using the MREPT technique. Even though the estimated weighting factors for MREIT and MREPT optimally reduce the random noise artifacts, non-uniform noise artifacts from various sources can severely deteriorate the reconstructed conductivity distributions. To be a clinical device used practically, it is important to develop noise reduction techniques taking into account the physical properties and experimental environments of MREIT and MREPT.

The designed phantom consists of two different kinds of anomalies to show the frequency dependent conductivity information. The stable apparent conductivity contrast of hollow cylindrical object made of thin insulating transparency film was only controlled by ion mobility through holes on the film wall. It excluded other effects caused by any ion concentration gradient. Therefore, MREIT images could only present the hollow cylindrical object with holes when comparing the phase terms{\mathcal{P}}^{+}and{\mathcal{P}}^{-}in 10. The apparent conductivity of the region inside cylinder was a non-zero conductivity value because the externally injected current entered the anomaly through the holes. In MREPT, the reconstructed high-frequency conductivity was same inside and outside of the anomaly since the measured phase of *H*^{+}only reflected the material conductivities inside and outside of the cylindrical film. Another agarose gel cylinder wrapped in an insulating film without holes shows different conductivity characteristics in MREPT and MREIT. Due to the short *T*
_{2}-decay relaxation time of agarose gel, the noise levels of{\mathit{\text{sd}}}_{{\phi}^{+,\chi}}and{\mathit{\text{sd}}}_{{B}_{z}^{\xi}}were relatively severe compared to those in the other regions. The noise level of each{\mathit{\text{sd}}}_{{\phi}_{j}^{+}}was strictly increasing with respect to the echo number depending on the *T*
_{2}-decay rate, whereas those of{\mathit{\text{sd}}}_{{B}_{{z}^{\phantom{\rule{0.3em}{0ex}}j}}}had different characteristics depending on{T}_{{c}^{\phantom{\rule{0.3em}{0ex}}j}}and the *T*
_{2}-decay rate simultaneously.

MREIT provides the conductivity distribution at a low frequency, whereas MREPT produces the conductivity at the Larmor frequency of 128 MHz at 3 T. The cell membranes consisted of the phospholipid bilayer with embedded proteins behave as a capacitor or insulating film. It shows a complicated pattern of conductivity depending on the measured frequencies. For the biological tissues, the membranes restrict the flow of current of low frequency. MREIT images reflect the membrane properties quantitatively and MREIT has a potential to visualize the anisotropic conductivity tensor map. On the other hand, the membranes become transparent at high frequency and provide a relatively degraded sensitivity in the conductivity image. MREPT may provide different conductivity characteristics compared to MREIT results. Based on the previous studies related to the complex conductivity spectra measured at multiple frequencies within the range of 10 Hz to several kHz, spectroscopic complex conductivity distribution can contribute to explain the physiological and pathological status of internal tissues [33]. The dual-frequency conductivity imaging by combination of MREPT and MREIT is an advanced technique to show different information in spite of displaying two extreme cases at low and the Larmor frequencies. It is meaningful to develop a method to produce multi-frequency conductivity image spectra with high spatial resolution and sensitivity using MR scanner within the range of DC to the Larmor frequency since significantly distinguishable signal changes in biological tissues depending on the applied frequency range.

We have a plan to support the clinically applicable combined conductivity imaging method based on MR techniques. We may apply this method to detect the destroyed cell membrane or the modification of cell or tissue structure due to necrosis or apoptosis initiated by inflammatory response inside the body. Since RF ablation and cryoablation for cancer treatments are accompanied by the destruction of cell membranes, it will be a good method to diagnose ablated lesions based on physiological status of tissue to improve the safety and predict the local recurrence after ablation.

## Conclusions

Cross-sectional conductivity imaging methods for high spatial resolution and sensitivity inside the human body has been actively investigated. MREIT and MREPT techniques show different electrical properties at low frequency (below 1 kHz) as MREIT uses an externally injecting current while MREPT uses the Larmor frequency of an MRI scanner. Recently, the dual-frequency conductivity imaging from a single MR scan simultaneously was proposed based on combination of MREPT and MREIT. Even though it produces conductivity spectral information, MREPT and MREIT have commonly suffered from weak signals and noise amplification since the procedures to reconstruct the conductivity in MREIT and MREPT need to differentiate the measured phase signals. We suggested the optimization method to find weighting factors according to echo signals for MREPT and MREIT using the ICNE multi-echo pulse sequence, which minimized the noise artifacts in the measured phase data. The noise standard deviations were precisely analyzed for the optimally weighted magnetic flux density and the phase term of positive-rotating magnetic field. We applied the denoising method based on the reaction-diffusion equation with the estimated noise standard deviations to enhance the quality of dual-frequency conductivity images. A real phantom experiment was performed to validate the proposed method using common measured data to reconstruct the dual-frequency conductivity distributions for MREPT and MREIT.

## References

Geddes LA, Baker LE:

**The specific resistance of biological material—a compendium of data for the biomedical engineer and physiologist.***Med Biol Eng*1967,**5**(3):271–293. 10.1007/BF02474537Gabriel S, Lau RW, Gabriel C:

**The dielectric properties of biological tissues: Ii. measurements in the frequency range 10 hz to 20 ghz.***Phys Med Biol*1996,**41**(11):2251–2269. 10.1088/0031-9155/41/11/002Holder DS:

*Electrical Impedance Tomography: Methods, History and Applications*. Bristol, UK: IOP Publishing; 2005.Ider YZ, Birgul O:

**Use of the magnetic field generated by the internal distribution of injected currents for electrical impedance tomography (mr-eit).***Elektrik*1998,**6**(3):215–225.Kwon O, Woo EJ, Yoon J-R, Seo JK:

**Magnetic resonance electrical impedance tomography (mreit): simulation study of j-substitution algorithm.***Biomed Eng, IEEE Trans*2002,**49**(2):160–167. 10.1109/10.979355Birgül Ö, Eyüboglu BM, Ider YZ:

**Current constrained voltage scaled reconstruction (ccvsr) algorithm for mr-eit and its performance with different probing current patterns.***Phys Med Biol*2003,**48**(5):653–671. 10.1088/0031-9155/48/5/307Ider YZ, Onart S, Lionheart WRB:

**Uniqueness and reconstruction in magnetic resonance–electrical impedance tomography (mr–eit).***Physiol Meas*2003,**24**(2):591–604. 10.1088/0967-3334/24/2/368Joy MLG:

**Mr current density and conductivity imaging: the state of the art.**In*Proceedings of the 26th Annual International Conference of the IEEE EMBS*. San Francisco, CA, USA; 1–5 September 2004:5315–5319.Muftuler LT, Hamamura M, Birgul O, Nalcioglu O:

**Resolution and contrast in magnetic resonance electrical impedance tomography (mreit) and its application to cancer imaging.***Technol Cancer Res & Treat*2004,**3**(6):599–609.Ödemir MS, Eyüboğlu BM, Öbek O:

**Equipotential projection-based magnetic resonance electrical impedance tomography and experimental realization.***Phys Med Biol*2004,**49**(20):4765–4783. 10.1088/0031-9155/49/20/008Seo JK, Pyo HC, Park C, Kwon O, Woo EJ:

**Image reconstruction of anisotropic conductivity tensor distribution in mreit: computer simulation study.***Phys Med Biol*2004,**49**(18):4371–4382. 10.1088/0031-9155/49/18/012Gao N, Zhu SA, He B:

**A new magnetic resonance electrical impedance tomography (mreit) algorithm: the rsm-mreit algorithm with applications to estimation of human head conductivity.***Phys Med Biol*2006,**51**(12):3067–3083. 10.1088/0031-9155/51/12/005Birgul O, Hamamura MJ, Muftuler LT, Nalcioglu O:

**Contrast and spatial resolution in mreit using low amplitude current.***Phys Med Biol*2006,**51**(19):5035–5049. 10.1088/0031-9155/51/19/020Hamamura MJ, Muftuler LT, Birgul O, Nalcioglu O:

**Measurement of ion diffusion using magnetic resonance electrical impedance tomography.***Phys Med Biol*2006,**51**(11):2753–2762. 10.1088/0031-9155/51/11/005Katscher U, Voigt T, Findeklee C, Vernickel P, Nehrke K, Dossel O:

**Determination of electric conductivity and local sar via b1 mapping.***Med Imaging, IEEE Trans*2009,**28**(9):1365–1374.Voigt T, Katscher U, Doessel O:

**Quantitative conductivity and permittivity imaging of the human brain using electric properties tomography.***Magn Reson Med*2011,**66**(2):456–466. 10.1002/mrm.22832Seo JK, Kim M-O, Lee J, Choi N, Woo EJ, Kim HJ, Kwon OI, Kim D-H:

**Error analysis of nonconstant admittivity for mr-based electric property imaging.***Med Imaging, IEEE Trans*2012,**31**(2):430–437.Lee BI, Oh SH, Woo EJ, Lee SY, Cho MH, Kwon O, Seo JK, Baek WS:

**Static resistivity image of a cubic saline phantom in magnetic resonance electrical impedance tomography (mreit).***Physiol Meas***24**(2):579–589.Seo JK, Yoon J-R, Woo EJ, Kwon O:

**Reconstruction of conductivity and current density images using only one component of magnetic field measurements.***Biomed Eng, IEEE Trans*2003,**50**(9):1121–1124. 10.1109/TBME.2003.816080Ider YZ, Onart S:

**Algebraic reconstruction for 3d magnetic resonance–electrical impedance tomography (mreit) using one component of magnetic flux density.***Physiol Meas*2004,**25**(1):281–294. 10.1088/0967-3334/25/1/032Hamamura MJ, Muftuler LT:

**Fast imaging for magnetic resonance electrical impedance tomography.***Magn Reson Imaging*2008,**26**(6):739–745. 10.1016/j.mri.2008.01.031Sajib SZK, Kim HJ, Kwon OI, Woo EJ:

**Regional absolute conductivity reconstruction using projected current density in mreit.***Phys Med Biol*2012,**57**(18):5841–5859. 10.1088/0031-9155/57/18/5841Kim HJ, Jeong WC, Sajib SZK, Kim MO, Kwon OI, Woo EJ, Kim DH:

**Simultaneous imaging of dual-frequency electrical conductivity using a combination of mreit and mrept.***Magn Reson Med*2014,**71**(1):200–208. 10.1002/mrm.24642Park C, Lee BI, Kwon O, Woo EJ:

**Measurement of induced magnetic flux density using injection current nonlinear encoding (icne) in mreit.***Physiol Meas*2007,**28**(2):117–127. 10.1088/0967-3334/28/2/001Nam HS, Kwon OI:

**Optimization of multiply acquired magnetic flux density b z using icne-multiecho train in mreit.***Phys Med Biol*2010,**55**(9):2743–2759. 10.1088/0031-9155/55/9/021Minhas AS, Jeong WC, Kim YT, Han Y, Kim HJ, Woo EJ:

**Experimental performance evaluation of multi-echo icne pulse sequence in magnetic resonance electrical impedance tomography.***Magn Reson Med*2011,**66**(4):957–965. 10.1002/mrm.22872Scott GC, Joy MLG, Armstrong RL, Henkelman RM:

**Sensitivity of magnetic-resonance current-density imaging.***J Magn Reson (1969)*1992,**97**(2):235–254. 10.1016/0022-2364(92)90310-4Sadleir R, Grant S, Zhang SU, Lee BI, Pyo HC, Oh SH, Park C, Woo EJ, Lee SY, Kwon O, Seo JK:

**Noise analysis in magnetic resonance electrical impedance tomography at 3 and 11 t field strengths.***Physiol Meas*2005,**26**(5):875–884. 10.1088/0967-3334/26/5/023Rudin LI, Osher S:

**Total variation based image restoration with free local constraints.**In*Proceedings of the First IEEE International Conference on Image Processing*. Austin, Texas; 13–16 November 1994:31–35.Perona P, Malik J:

**Scale-space and edge detection using anisotropic diffusion.***Pattern Anal Machine Intell, IEEE Trans*1990,**12**(7):629–639. 10.1109/34.56205Oh TI, Kim YT, Minhas A, Seo JK, Kwon OI, Woo EJ:

**Ion mobility imaging and contrast mechanism of apparent conductivity in mreit.***Phys Med Biol*2011,**56**(7):2265–2277. 10.1088/0031-9155/56/7/022Park C, Lee BI, Kwon OI:

**Analysis of recoverable current from one component of magnetic flux density in mreit and mrcdi.***Phys Med Biol*2007,**52**(11):3001–3013. 10.1088/0031-9155/52/11/005Oh TI, Koo H, Lee KH, Kim SM, Lee J, Kim SW, Seo JK, Woo EJ:

**Validation of a multi-frequency electrical impedance tomography (mfeit) system khu mark1: impedance spectroscopy and time-difference imaging.***Physiol Meas*2008,**29**(3):295. 10.1088/0967-3334/29/3/002

## Acknowledgements

This paper was supported by Konkuk University in 2013.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

OIK developed the algorithm and wrote the manuscript. WCJ prepared and performed the experiments and wrote up the experimental section. SZKS analyzed the experimental data and confirmed it by simulation. HJK assisted with analysis and processing the data. EJW conceived the idea and gave critical revision for important intellectual content. TIO designed the research topic and experiments, analyzed and drafted the manuscript. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access**
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (
https://creativecommons.org/publicdomain/zero/1.0/
) applies to the data made available in this article, unless otherwise stated.

## About this article

### Cite this article

Kwon, O.I., Jeong, W.C., K Sajib, S.Z. *et al.* Reconstruction of dual-frequency conductivity by optimization of phase map in MREIT and MREPT.
*BioMed Eng OnLine* **13, **24 (2014). https://doi.org/10.1186/1475-925X-13-24

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1475-925X-13-24

### Keywords

- MRI
- MREIT
- MREPT
- Conductivity
- Magnetic flux density
- Optimization