Conductivity image enhancement in MREIT using adaptively weighted spatial averaging filter
© Oh et al.; licensee BioMed Central Ltd. 2014
Received: 20 April 2014
Accepted: 24 June 2014
Published: 26 June 2014
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.
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.
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.
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.
KeywordsMagnetic resonance electrical impedance tomography Conductivity image Noise estimation Denoising method Adaptively weighted spatial averaging filter
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 . Various multi-echo and multi-coil techniques have been developed together with optimization methods to combine multiple signals from multiple coils and echoes .
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 . 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 .
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.
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.
where M is the MR magnitude image, δ is any systematic phase artifact, γ=26.75×107rad/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-.
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 decay of the MR signal.
Conductivity image reconstruction
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  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 , respectively.
where for a given scalar function f and 〈·,·〉 denotes the scalar inner product. We used the computed values of two magnetic flux densities and corresponding to two computed voltages and , 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 .
Adaptively weighted spatial averaging
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.
For a small amount of noise, both h(r) and η(r) are small and, therefore, the relation (12) implies that both and 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 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 is considerably reduced by the spatial averaging of the random noise.
Variance of the reconstructed conductivity values in the ROIs using the variance function in ( 15)
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).
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.
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).
- 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
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.