- Research
- Open Access

# Conductivity image enhancement in MREIT using adaptively weighted spatial averaging filter

- Tong In Oh
^{1}, - Hyung Joong Kim
^{1}, - Woo Chul Jeong
^{1}, - Hun Wi
^{1}, - Oh In Kwon
^{2}Email author and - Eung Je Woo
^{1}

**13**:87

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

© Oh et al.; licensee BioMed Central Ltd. 2014

**Received:**20 April 2014**Accepted:**24 June 2014**Published:**26 June 2014

## Abstract

### Background

In magnetic resonance electrical impedance tomography (MREIT), we reconstruct conductivity images using magnetic flux density data induced by externally injected currents. Since we extract magnetic flux density data from acquired MR phase images, the amount of measurement noise increases in regions of weak MR signals. Especially for local regions of MR signal void, there may occur excessive amounts of noise to deteriorate the quality of reconstructed conductivity images. In this paper, we propose a new conductivity image enhancement method as a postprocessing technique to improve the image quality.

### Methods

Within a magnetic flux density image, the amount of noise varies depending on the position-dependent MR signal intensity. Using the MR magnitude image which is always available in MREIT, we estimate noise levels of measured magnetic flux density data in local regions. Based on the noise estimates, we adjust the window size and weights of a spatial averaging filter, which is applied to reconstructed conductivity images. Without relying on a partial differential equation, the new method is fast and can be easily implemented.

### Results

Applying the novel conductivity image enhancement method to experimental data, we could improve the image quality to better distinguish local regions with different conductivity contrasts. From phantom experiments, the estimated conductivity values had 80% less variations inside regions of homogeneous objects. Reconstructed conductivity images from upper and lower abdominal regions of animals showed much less artifacts in local regions of weak MR signals.

### Conclusion

We developed the fast and simple method to enhance the conductivity image quality by adaptively adjusting the weights and window size of the spatial averaging filter using MR magnitude images. Since the new method is implemented as a postprocessing step, we suggest adopting it without or with other preprocessing methods for application studies where conductivity contrast is of primary concern.

## Keywords

- Magnetic resonance electrical impedance tomography
- Conductivity image
- Noise estimation
- Denoising method
- Adaptively weighted spatial averaging filter

## Background

Electrical conductivity is a passive material property of a biological tissue or organ providing diagnostic information of its physiological function and pathological state [1, 2]. In magnetic resonance electrical impedance tomography (MREIT), we aim to visualize the internal conductivity distribution of the human body by injecting electrical current and measuring induced magnetic flux density data using an MRI scanner [3, 4]. The MREIT technique uses the *z*-component *B*
_{
z
} of the magnetic flux density **B**=(*B*
_{
x
},*B*
_{
y
},*B*
_{
z
}) to recover conductivity and/or current density images [5–10].

In experimental MREIT studies of animal and human subjects [11–14], the quality of reconstructed conductivity images highly depends on the noise level in measured *B*
_{
z
} data and the injection current amplitude. For a given current amplitude, we can improve the image quality by reducing the noise level. If we reduce the current amplitude to avoid adverse effects such as electrical stimulations of nerve and muscle, the range of *B*
_{
z
} decreases proportionally and the measured *B*
_{
z
} data become more vulnerable to noise. In order not to deteriorate the conductivity image quality, we should reduce the noise level as well.

There have been numerous studies to minimize the noise level in measured *B*
_{
z
} data by optimizing the data collection method including pulse sequences, RF coils, shimming, averaging, and so on. For example, the injected current nonlinear encoding (ICNE) method was introduced to reduce the noise level in *B*
_{
z
} data by extending the current injection time without increasing the scan time [15]. Various multi-echo and multi-coil techniques have been developed together with optimization methods to combine multiple signals from multiple coils and echoes [16].

Once we acquire *k*-space data from a designed MREIT experiment, we should carefully process the data to suppress the noise and enhance the image quality. This requires proper understanding of the noise characteristics. The noise standard deviation of *B*
_{
z
} is inversely proportional to the signal-to-noise ratio (SNR) of the MR magnitude image and the current injection time when the magnitude image SNR is higher than about 3 [17, 18]. Since the magnitude image SNR varies over the image, the noise level of *B*
_{
z
} also varies. If the magnitude image SNR is much lower, for example, less than 3, then we should expect an excessive amount of noise in *B*
_{
z
}. This may happen when we scan an animal or human subject since there exist local regions of weak MR signal such as the lungs, gas-filled organs, and outer layers of bones. In those regions with very low magnitude image SNRs, there occur excessive amounts of noise in measured *B*
_{
z
} data. It is, therefore, important to properly suppress the noise effects in preprocessing or postprocessing steps.

There are several denoising methods to process the acquired *B*
_{
z
} data [19–23]. These methods use either a PDE-based approach or a localized data processing approach and require a certain degree of manual adjustment of parameters. All of these methods are applied to measured *B*
_{
z
} data before conductivity image reconstructions as a preprocessing step.

In this paper, we present a new method to suppress the noise effects in reconstructed conductivity images as a postprocessing step. We suggest incorporating *a prior* information from the MR magnitude image which is always available in MREIT. For current density imaging using measured magnetic flux density data, Joy *et al*. suppressed noise effects by setting some appropriate threshold in the MR magnitude image SNR [24]. A similar approach based on the morphological enabled dipole inversion (MEDI) has been proposed in quantitative susceptibility map (QSM) to overcome the ill-posedness in its inverse problem [25].

The proposed method enhances the quality of reconstructed conductivity images by using an adaptively weighted spatial averaging filter. To determine the weights, we will define a distance function with respect to the noise standard deviation in measured *B*
_{
z
} data. Since the noise is inversely proportional to the MR magnitude image, we will incorporate the magnitude image in the spatial filtering process. After describing the details of the proposed method, we will show how it performs with experimental data from a conductivity phantom and also animal subjects.

## Methods

### Noise in magnetic flux density data

We inject current into an imaging object, of which timing is synchronized with a chosen MR pulse sequence. The usual current injection period is from the end of the exciting RF pulse to the beginning of the readout gradient. One may choose different current injection times and patterns depending on designed pulse sequences. For example, the ICNE-multi-echo pulse sequence extends the duration of current injection until the end of multiple read-out gradients.

*z*-component

*B*

_{ z }in a voxel results in an extra phase in the MR phase image. To remove systematic phase artifacts, we sequentially inject positive

*I*

^{+}and negative

*I*

^{-}current pulses to obtain the following

*k*-space data:

where *M* is the MR magnitude image, *δ* is any systematic phase artifact, *γ*=26.75×10^{7}rad/T·s is the gyromagnetic ratio of the proton, and *Ω* is a field-of-view (FOV). Here, the superscript of *S*^{±} denotes a brief notation for *S*^{+} and *S*^{-}.

*B*

_{ z }by

where ${\Gamma}^{\pm}=M{e}^{\mathrm{i\delta}}{e}^{\pm \mathrm{i\gamma}{B}_{z}{T}_{c}}$.

*B*

_{ z }is inversely proportional to the current injection time

*T*

_{ c }and the SNR of the MR magnitude image

*Υ*as

for *Υ*>2.8 [17, 18]. If we reduce the current amplitude *I*^{±}, the range of *B*
_{
z
} decreases proportionally. To reduce the noise standard deviation, we have to increase the current injection time *T*
_{
c
} and the SNR *Υ* simultaneously. However, this is not possible mainly due to the *T*
_{2} or ${T}_{2}^{\ast}$ decay of the MR signal.

### Conductivity image reconstruction

*I*induces distributions of voltage

*u*, current density

**J**=(

*J*

_{ x },

*J*

_{ y },

*J*

_{ z }), and magnetic flux density

**B**=(

*B*

_{ x },

*B*

_{ y },

*B*

_{ z }) inside the imaging object

*Ω*with its boundary

*∂*

*Ω*. The voltage satisfies the following partial differential equation:

**r**=(

*x*,

*y*,

*z*),

**n**is the outward unit normal vector on

*∂*

*Ω*, and

*g*denotes the Neumann boundary data subject to the injection current. From the Biot-Savart law, the magnetic flux density

*B*

_{ z }is related with the current density

**J**=-

*σ*∇

*u*as

where *μ*
_{0} is the magnetic permeability of the free space.

There exists numerous image reconstruction algorithms to visualize the conductivity *σ* or the current density **J** in the imaging object from measured *B*
_{
z
} data [5, 9, 16, 26]. In this paper, we adopted the transversal *J*-substitution algorithm [27] since it differentiates the noisy *B*
_{
z
} data once and does not propagate the noise effects from one region to another. We first assume the imaging object with a homogeneous conductivity distribution *σ*^{
H
}. Solving (4) with *σ*^{
H
} in place of *σ*, we can compute the voltage, current density, and magnetic flux density, which are denoted as *u*^{
H
}, **J**^{
H
}, and ${B}_{z}^{H}$, respectively.

*B*

_{ z }data subject to two orthogonal injection currents. Denoting them as

*B*

_{ z,1}and

*B*

_{ z,2}, the transversal

*J*-substitution algorithm produces an image of the conductivity

*σ*as

where ${\stackrel{~}{\nabla}}^{\perp}f:=\left(\frac{\mathrm{\partial f}}{\mathrm{\partial y}},-\frac{\mathrm{\partial f}}{\mathrm{\partial x}}\right)$ for a given scalar function *f* and 〈·,·〉 denotes the scalar inner product. We used the computed values of two magnetic flux densities ${B}_{z,1}^{H}$ and ${B}_{z,2}^{H}$ corresponding to two computed voltages ${u}_{1}^{H}$ and ${u}_{2}^{H}$, respectively, for the homogeneous case of *σ*^{
H
}. Note that the noisy data of *B*
_{
z,1} and *B*
_{
z,2} may deteriorate the quality of the reconstructed image of *σ*. More details about the image reconstruction algorithm are described in [27].

### Adaptively weighted spatial averaging

*M*(

**r**) is the MR magnitude image,

*B*

_{ r }(

*η*(

**r**)) is a neighborhood of the point

**r**with a radius

*η*(

**r**), and

*h*(

**r**) is a function of the noise level in measured

*B*

_{ z }data. For each pixel

**s**in the neighborhood of the pixel

**r**, that is,

**s**∈

*B*

_{ r }(

*η*(

**r**)), we define the following weighting factor

*w*

_{ r }(

**s**):

**r**, the distance

*D*(

**r**,

**s**) measures the similarity of

*M*(

**r**) and

*M*(

**s**) in the surrounding region of

**r**. Since the noise standard deviation of

*B*

_{ z }is inversely proportional to the magnitude image SNR, the weights are in effect adjusted by the noise level in measured

*B*

_{ z }data. In (7), two parameters of

*η*(

**r**) and

*h*(

**r**) regulate the extent of the spatial averaging filter. We now design the denominator

*h*(

**r**) and the radius

*η*(

**r**) in (7) as

We can choose the proportionality parameters of *h*(**r**) and *η*(**r**) in (10) depending on the quality of the reconstructed conductivity image to be filtered.

*D*(

**r**,

**s**) preserves the conductivity value where the noise level of

*B*

_{ z }is low. To determine the characteristics of the weighting factor

*w*

_{ r }(

**s**), we decompose the weighted conductivity

*σ*

_{ w }(

**r**) to the true conductivity

*σ*

_{ t }(

**r**) and the noise term:

*σ*

_{ t }(

**s**) and

*N*

_{ r }(

**s**) denote the noiseless true conductivity and the noise term at

**s**, respectively. The weighted conductivity

*σ*

_{ w }(

**r**) and the true conductivity

*σ*

_{ t }(

**r**) satisfy the following relation:

For a small amount of noise, both *h*(**r**) and *η*(**r**) are small and, therefore, the relation (12) implies that both ${\mathcal{E}}_{1}\left(\mathbf{r}\right)$ and ${\mathcal{E}}_{2}\left(\mathbf{r}\right)$ are small. This means that the method does no harm to the conductivity image when the noise level is low.

For a large amount of noise, the first error term ${\mathcal{E}}_{1}\left(\mathbf{r}\right)$ becomes small if there was no significant variations of the conductivity values within the neighboring region *B*
_{
r
}(*η*(**r**)) of the pixel at **r**. This is true only when the pixels with similar magnitude image values also have similar conductivity values. We will discuss implications of this restriction later. The second error term ${\mathcal{E}}_{2}\left(\mathbf{r}\right)$ is considerably reduced by the spatial averaging of the random noise.

### Imaging experiments

#### Phantom imaging

*T*

_{ R }=1200 ms, echo time

*T*

_{ E }=15 ms, number of echoes

*N*

_{ E }=5, field of view (FOV) was 240×240 mm

^{2}, number of excitations (NEX) was 1, imaging matrix was 128×128, and total imaging time was 10.24 min. Figure 1(b) shows the injected current pulses synchronized with RF pulses. We combined the acquired multiple echoes to optimize the SNR in measured

*B*

_{ z }data [16]. Figure 1(c) shows the measured

*B*

_{ z }image subject to the vertical current injection.

#### Animal imaging

^{3}. The procedure was approved by the Institutional Animal Care and Use Committee (IACUC). We injected currents in two mutually orthogonal directions between two pairs of electrodes facing each other. The injection current amplitude ranged from 5 to 10 mA. We adopted the ICNE pulse sequence. The imaging parameters were as follows: repetition time

*T*

_{ R }=1200 ms, echo time

*T*

_{ E }=30 ms, FOV was 280×280 mm

^{2}, slice thickness was 4 mm, number of slice was 8, NEX was 6, imaging matrix was 128×128, and total imaging time was 60 min. Figure 2 shows two data sets we used in the paper.

## Results

*B*

_{ z }data, we determined the denominator

*h*in (7) to compute the distance function

*D*(

**r**,

**s**). Figure 3(b) shows the denominator

*h*in (7). Figure 3(c) shows the reconstructed conductivity image using the algorithm in (6). Due to the short

*T*

_{2}relaxation time of the objects, the reconstructed conductivity values of the objects show strong noise defects. Figure 3(d) shows the filtered conductivity image using the proposed method. In the noisy regions of ${\mathcal{D}}_{i},\phantom{\rule{1em}{0ex}}i=1,2,3$, the values of the denominator

*h*were relatively large to produce smaller distance values in (7). Since the values of

*η*were also large in ${\mathcal{D}}_{i}$, more pixels in the neighboring region were included in the spatial averaging with relatively larger weights.

**Variance of the reconstructed conductivity values in the ROIs using the variance function in (** 15**)**

${V}_{{\mathcal{D}}_{1}}\left(\sigma \right)$ | ${V}_{{\mathcal{D}}_{2}}\left(\sigma \right)$ | ${V}_{{\mathcal{D}}_{3}}\left(\sigma \right)$ | |
---|---|---|---|

Before filtering | 25.4 | 15.7 | 18.4 |

After filtering | 3.39 | 3.89 | 4.81 |

*h*in (7) for the upper abdominal region in Figure 2(a). The areas marked by the arrows have large values of

*h*since their MR magnitude image SNRs were low. The average magnitude intensity in the marked areas was 3.78, where the intensity of the entire image ranged from 0 to 30. Since the noise level of

*B*

_{ z }is inversely proportional to the MR magnitude intensity, the recovered conductivity values of the marked regions included a relatively large amount of noise. Figure 4(b) shows the reconstructed conductivity image without using any filtering method and we can observe severe noise effects in those regions.

where *f* is the reconstructed conductivity image shown in Figure 4(b). Figure 4(c) shows the denoised conductivity image using the iteration method (16). We can see that the iteration method blurred the entire image. Figure 4(d) is the conductivity image after applying the adaptively weighted spatial averaging filter. We can see that the proposed method enhanced the conductivity image in the marked regions without affecting other regions. For the blurred image in Figure 4(c), we used the following parameters: iteration number = 200, *α*=0.1, and fidelity term *β*=0.01. These parameters were chosen to make the variance of the conductivity values in the marked regions of the image in Figure 4(c) to be equal to that of the image in 4(d).

*B*

_{ z }data were noisier than the case in Figure 4. Comparing two filtered conductivity images in 5(c) and (d), we can see that the method proposed in this paper is clearly advantageous.

## Discussion and conclusion

For animal or human subjects, their internal structures are heterogenous and there often exist local regions of weak MR signals. This makes the noise level in measured magnetic flux density data vary significantly for different pixels. Since we use the data to reconstruct conductivity images, noise effects are conveyed to the conductivity images with spatially varying image quality. Conventional filtering methods without considering this property, therefore, unnecessarily blur the entire conductivity image.

In this paper, we proposed a fast and simple method to enhance the conductivity image quality by utilizing the MR magnitude image. Noting that the noise level in measured magnetic flux density data is inversely proportional to the pixel value of the MR magnitude image, we could adaptively adjust the weights and window size of the spatial averaging filter. In designing the filter, we used two parameters of *η*(**r**) and *h*(**r**) in (7) to regulate the extent of the spatial averaging filter. Since we set their values to be proportional to the noise level in measured *B*
_{
z
} data, the proposed method adaptively changes the extent of the spatial averaging depending on the noise level at each pixel. It is important to properly choose the proportionality constants for *η*(**r**) and *h*(**r**) in (10). Though we heuristically chose the constants in this paper, we plan to rigorously investigate their effects on the filtered image and develop an automatic method to determine them.

Since the spatial averaging is performed on neighboring pixels with similar MR magnitude image values, it may remove any useful conductivity contrast among those pixels. Using the proposed method, we can apply the spatial averaging to local regions of severe noise where we can not trust reconstructed conductivity values. We can, therefore, suppress noise for the price of reduced image resolution within local regions with large amounts of noise. The proposed method negligibly influences the conductivity image in other regions with enough SNR.For the phantom image in Figure 3, we found that the spatial averaging recovered the correct conductivity values in the noisy regions since there was no conductivity contrast within each region. For the animal images in Figures 4 and 5, it was difficult to assert that the recovered conductivity values of the noisy regions were correct since there could have been some conductivity contrast among those pixels with similar MR magnitude image values. In application studies, we may analyze both unfiltered and filtered conductivity images if quantitative conductivity values of noisy regions are of primary concern.

While the existing denoising methods in MREIT are applied to measured *B*
_{
z
} data before conductivity image reconstructions [19–23], the new method proposed in this paper is a postprocessing method, which can be applied to reconstructed conductivity images. Without requiring any image segmentation, one can easily implement the method without or with chosen preprocessing steps. We suggest adopting the proposed method in future experimental MREIT studies where conductivity contrast is of primary concern to extract diagnostic information.

## Declarations

### Acknowledgements

This paper was supported the Konkuk University research support program and by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) (2011-0022779, 2012R1A1A2008477, 2013R1A2A2A04016066, 2014R1A2A1A09006320).

## Authors’ Affiliations

## References

- Grimnes S, Martinsen OG:
*Bioimpedance and Bioelectricity Basics, 2nd Ed*. Oxford: Academic Press; 2008.Google Scholar - Seo JK, Kim DH, Lee J, Kwon OI, Sajib SZK, Woo EJ:
**Electrical tissue property imaging using MRI at dc ad Larmor frequency.***Inv Prob*2012,**28**(8):0840022.MathSciNetView ArticleGoogle Scholar - Woo EJ, Seo JK:
**Magnetic Resonance Electrical Impedance Tomography (MREIT) for high-resolution conductivity imaging.***Physiol Meas*2008,**29**(10):R1-R26.View ArticleGoogle Scholar - Seo JK, Woo EJ:
*Nonlinear Inverse Problems in Imaging*. New York: Wiley; 2012.Google Scholar - Seo JK, Yoon JR, Woo EJ, Kwon O:
**Reconstruction of conductivity and current density images using only one component of magnetic field measurements.***IEEE Trans Biomed Eng*2003,**50**(9):1121–1124.View ArticleGoogle Scholar - Birgul O, Hamamura M, Muftuler L, Nalcioglu O:
**Contrast and spatial resolution in MREIT using low amplitude current.***Phys Med Biol*2006,**51**(19):5035–5049.View ArticleGoogle Scholar - Gao N, Zhu SA, He BA:
**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.View ArticleGoogle Scholar - Hamamura M, Muftuler L, Birgul O, Nalcioglu O:
**Measurement of ion diffusion using magnetic resonance electrical impedance tomography.***Phys Med Biol*2006,**51**(11):2753–2762.View ArticleGoogle Scholar - Hamamura M, Muftuler L:
**Fast imaging for magnetic resonance electrical impedance tomography.***Magn Reson Imaging*2008,**26**(6):739–745.View ArticleGoogle Scholar - Muftuler LT, Chen G, Hamamura MJ, Ha SH:
**MREIT with SENSE acceleration using a dedicated RF coil design.***Physiol Meas*2009,**30**(9):913–929.View ArticleGoogle Scholar - Kim HJ, Oh TI, Kim YT, Lee BI, Woo EJ, Seo JK, Lee SY, Kwon O, Park C, Kang BT, Park HM:
**In vivo electrical conductivity imaging of a canine brain using a 3 T MREIT system.***Physiol Meas*2008,**29**(10):1145–1155.View ArticleGoogle Scholar - Kim HJ, Kim YT, Minhas AS, Jeong WC, Woo EJ, Seo JK, Kwon OJ:
**In vivo high-resolution conductivity imaging of the human leg using MREIT: the first human experiment.***IEEE Trans Med Imag*2009,**28**(11):1681–1687.View ArticleGoogle Scholar - Jeon K, Minhas AS, Kim YT, Jeong WC, Kim HJ, Kang BT, Park HM, Lee CO, Seo JK, Woo EJ:
**MREIT conductivity imaging of the postmortem canine abdomen using CoReHA.***Physiol Meas*2009,**30**(9):957–966.View ArticleGoogle Scholar - Oh TI, Jeong WC, McEwan A, Park HM, Kim HJ, Kwon OI, Woo EJ:
**Feasibility of Magnetic Resonance Electrical Impedance Tomography (MREIT) conductivity imaging to evaluate brain abscess lesion: in vivo canine model.***J Mag Reson Imag*2013,**38**(1):189–197.View ArticleGoogle Scholar - Park C, Lee BI, Kwon O, Woo EJ:
**Measurement of induced magnetic flux density using injection current nonlinear encoding (ICNE) in MREIT.***Physiol Meas*2006,**28**(2):117–127.View ArticleGoogle Scholar - Nam 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.View ArticleGoogle Scholar - Scott GC, Joy MLG, Armstrong RL, Henkelmantitle RM:
**Sensitivity of magnetic-resonance current-density imaging.***J Mag Res*1992,**97**(2):235–254.Google Scholar - Sadleir 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.View ArticleGoogle Scholar - Lee BI, Lee SH, Kim TS, Kwon O, Woo EJ, Seo JK:
**Harmonic decomposition in PDE-based denoising technique for magnetic resonance electrical impedance tomography.***IEEE Trans Biomed Eng*2005,**52**(11):1912–1920.View ArticleGoogle Scholar - Lee S, Seo JK, Park C, Lee BI, Woo EJ:
**Conductivity image reconstruction from defective data in MREIT: numerical simulation and animal experiment.***IEEE Trans Med Imag*2006,**25**(2):168–176.View ArticleGoogle Scholar - Jeon K, Kim HJ, Lee CO, Seo JK, Woo EJ:
**Integration of denoising, inpainting and local harmonic Bz algorithm for MREIT imaging of intact animals.***Phys Med Biol*2010,**55**(24):7541–7556.View ArticleGoogle Scholar - Lee CO, Jeon K, Ahn S, Kim HJ, Woo EJ:
**Ramp-preserving denoising for conductivity image reconstruction in magnetic resonance electrical impedance tomography.***IEEE Trans Bio Eng*2011,**58**(7):2038–2050.View ArticleGoogle Scholar - Sajib SZK, Kim HJ, Kwon OI, Woo EJ:
**Regional absolute conductivity reconstruction using projected current density in MREIT.***Physics Med Biol*2012,**57**(18):5841–5859.View ArticleGoogle Scholar - Joy MLG, Lebedev VP, Gati JS:
**Imaging of durrent density and durrent pathways in rabbit brain during transcranial electrostimulation.***IEEE Trans Bio Eng*1999,**46**(7):1139–1149.View ArticleGoogle Scholar - Liu T, Liu J, de Rochefort L, Spincemaille P, Khalidov I, Ledoux JR, Wang Y:
**Morphology enabled dipole inversion (MEDI) from a single-angle acquisition: comparison with COSMOS in human brain imaging.***Magn Reson Med*2011,**66**(3):777–783.View ArticleGoogle Scholar - Ider 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.View ArticleGoogle Scholar - Nam HS, Lee BI, Choi J, Park C, Kwon OI:
**Conductivity imaging with low level current injection using transversal J-substitution algorithm in MREIT.***Phys Med Biol*2007,**52:**6717–6730.View ArticleGoogle Scholar - Kim YT, Yoo PJ, Oh TI, Woo EJ:
**Magnetic flux density measurement in magnetic resonance electrical impedance tomography using a low-noise current source.***Meas Sci Technol*2011,**22**(10):105803.View ArticleGoogle Scholar - Rudin L, Osher S:
**Total variation based image restoration with free local constraints.***Proc 1st IEEE ICIP*1994,**1:**31–35.Google Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.