Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 941745, 8 pages

http://dx.doi.org/10.1155/2013/941745

Research Article

CoReHA 2.0: A Software Package for In Vivo MREIT Experiments

Kiwan Jeon¹ and Chang-Ock Lee²

¹National Institute for Mathematical Sciences, Daejeon 305-811, Republic of Korea

²Department of Mathematical Sciences, KAIST, Daejeon 305-701, Republic of Korea

Received 12 December 2012; Accepted 20 January 2013

Academic Editor: Jin Keun Seo

Copyright © 2013 Kiwan Jeon and Chang-Ock Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Magnetic resonance electrical impedance tomography (MREIT) is a new medical imaging modality visualizing static conductivity images of electrically conducting subjects. Recently, MREIT has rapidly progressed in its theory, algorithm, and experiment technique and now reached to the stage of in vivo animal experiments. In this paper, we present a software, named CoReHA 2.0 standing for the second version of conductivity reconstructor using harmonic algorithms, to facilitate in vivo MREIT reconstruction of conductivity image. This software offers various computational tools including preprocessing of MREIT data, identification of 2D geometry of the imaging domain and electrode positions, and reconstruction of cross-sectional scaled conductivity images from MREIT data. In particular, in the new version, we added several tools including ramp-preserving denoising, harmonic inpainting, and local harmonic algorithm to deal with data from in vivo experiments. The presented software will be useful to researchers in the field of MREIT for simulation, validation, and further technical development.

1. Introduction

Recently, a new imaging modality called magnetic resonance electrical impedance tomography (MREIT) has been introduced, which allows high resolution imaging of tomographic electrical conductivity distributions of biological objects [1, 2]. The technique involves (i) current injection into an electrically conducting object such as animal or human body through surface electrodes, (ii) measurement of induced internal magnetic flux density using an MRI system, typically only the -component of the induced magnetic flux density , where is the axis parallel to the main magnetic field of the MR scanner, and (iii) conductivity reconstruction by solving nonlinear boundary value problems with given injected currents and measured magnetic flux density, employing finite element methods. Although it looks straight forward to reconstruct conductivity distribution as described in the three steps, computation involves several innovative approaches including magnetic flux density estimation, data verification, segmentation, and solving forward/inverse problems. These specific computations cannot be handled through readily available finite element packages; thus, it is required to develop a user friendly software with graphics user interface (GUI) for those who wish to reconstruct conductivity distributions. Upon these requests, we developed a software package, called CoReHA [3, 4] which stands for conductivity reconstructor using harmonic algorithms, using VC++ MFC 6.0 (Microsoft Foundation Class Library 6.0) and OpenGL under the Microsoft Windows operating system. Based on the harmonic algorithm [5], CoReHA supports all procedures from the preprocessing of raw data to conductivity imaging through the intuitively apprehensible graphic user interface, more specifically, data conversion of raw -space data, data verification, segmentation tools for numerical computation, solvers of forward/inverse problems using finite element methods, and 2D/3D data view as well as histogram. This software has been a major tool to facilitate multilateral studies for MREIT and has brought out successful reconstruction of conductivity imaging in many researches [4, 6–8].

In the stage of in vivo MREIT experiments, we have to address the improvement of the signal-to-noise ratio (SNR) of measured data, since data has weak strength due to the low amount of injected current for the safety guide. Moreover, there may exist MR signal void regions in the animal or human body, where noise levels of are excessively high, and as a result, uncertain effects are produced. For dealing with these technical problems, several methods were developed, called ramp preserving denoising [9], harmonic inpainting [10], and local harmonic algorithm [11]. Incorporating these new features, we release the second version of CoReHA for in vivo MREIT experiments. The presented software will be useful to researchers in the field of MREIT experimental studies as well as in vivo animal/human experiments. CoReHA is available from the website http://iirc.khu.ac.kr/.

This paper is organized as follows. In Section 2, MREIT system is introduced and related works are explained. Brief explanation of the previous version of this software is given for comprehensive understanding of our works in Section 3. Details of new tools of the software for dealing with in vivo stage issues are provided in Section 4. Finally, we close in Section 5 with conclusions and a discussion of future works.

2. The Basic of MREIT and Related Works

We briefly explain how MR scanner is used as a tool to capture internal magnetic flux density images, and the MR scanner has its magnetic field in -direction. Let an electrically conducting subject occupy a three-dimensional domain with its boundary and a conductivity distribution . As shown in Figure 1, we attach pairs of electrodes and along in order to inject a current in a form of pulses whose timing is synchronized with an MR pulse sequence for . The injection current produces current density satisfying the following elliptic equations:

where . Also, magnetic flux density is given by Biot-Savart law as follows:

where is magnetic flux density from electrodes, wires, and others. Note that in . Then MR spectrometer provides the complex -space data that is influenced by in the following way. Given current injection , we have

where is a conventional MR magnitude image, any systematic phase artifact,  rad/T·s the gyromagnetic ratio of hydrogen, and the current pulse width in seconds. By the discrete inverse Fourier transformation of the -space data in (3), we obtain the following complex images:

To extract data stably while eliminating systematic phase artifact , we inject the counter directional current providing the counter part of in (4) so that we get

Division of by its counter part leads to

Figure 1: Domain for MREIT.

MRCDI [12] and early stage MREIT [13–15] use measurements of all three components of which require subject rotations inside MR scanner. However, experiences show that these subject rotations are impractical and also cause other problems such as misalignments of pixels. Hence, in order to make the MREIT technique easily applicable to clinical situations, we should use only data for the conductivity reconstruction.

In the harmonic algorithm [5, 16, 17] which is the first constructive -based MREIT algorithm, we inject two independent electrical currents and through two pairs of surface electrodes and , respectively, to get two raw -space data. The amount of each injection current is determined by the product of current amplitude and time duration. A modified spin-echo pulse sequence synchronized with the current injection is typically used as described in experimental works [5, 18, 19], and and are obtained from the formula (6). Then the reconstruction algorithm is based on the following identity:

where

and and are the and components of induced current density , respectively. We now present the conductivity reconstruction procedure using the harmonic algorithm in CoReHA, which consists of six steps.

Step 1. We impress electrical currents and through pairs of surface electrodes and , respectively, and get the -space data for using an MR scanner.

Step 2. We produce an MR magnitude image and induced magnetic flux densities from the -space data for given by (6).

Step 3. Using the MR magnitude image , we perform segmentation of and . Here, we use level-set-based segmentations.

Step 4. We set the initial guess in .

Step 5. We solve the problem (1) with .

Step 6. We solve

with the boundary condition

where and are two-dimensional Laplacian and gradient operators on the -plane, respectively, a two-dimensional slice of , which is perpendicular to the -axis, and the two-dimensional outward normal vector to .

Note that in the harmonic algorithm in [5, 16], Steps 5 and 6 are repeated to get from until it converges. In order to ensure the convergence of the harmonic algorithm, the conductivity values of the subject’s boundary should be homogeneous and the contrast of the conductivity distribution is assumed small enough [16, 17]. But, in general, the conductivity distribution of animal and human is quite inhomogeneous [11]. Also, for the completion of the iterative procedure, we have to resolve heavy cost computational issues for 3D forward problem (1), related to, for example, 3D segmentation, 3D mesh generation, and 3D forward problem solver. Fortunately, for the clinical purpose, it is enough to find the scaled conductivity images instead of the true ones because the reconstructed scaled ones reflect the fine details of the true conductivity contrast when the conductivity contrast is low [11]. Hence, in CoReHA, we iterate Steps 5 and 6 only once. Due to the same reason, in Step 5, we solve a simplified problem

instead of (1). In addition, in Step 6, we adopt the boundary condition

to replace (10). Note that three-dimensional Laplacian of may not be invalid due to the circumstance of the experiments. Hence, CoReHA supports for both 2D and 3D Laplacian of .

3. CoReHA 1.0: For Multilateral Studies

CoReHA implements Step 2 through 6 of the harmonic algorithm since Step 1 is about the experiment using an MR scanner and EIT equipment. In this section, we briefly describe each step of CoReHA 1.0 for understanding the software package on the whole.

After an MREIT experiment, one can perform the phase extraction process by applying fast Fourier transform in Section 2 as well as obtain the MR imaging by taking magnitude of . For this computation, we use FFTW [20], one of the well-known fast Fourier transform libraries, for obtaining (4) from (3). Since in (6) is wrapped due to the branch cut of argument operator, we apply the Goldstein’s algorithm [21] for two-dimensional phase unwrapping. Note that even if data is continuous in the -plane given by phase unwrapping, we do not guarantee the continuity along the -direction. In order to verify the continuity of along the -direction, CoReHA supports a verification tool. Hence, if the verification fails, one should drop in for preventing the artifact coming from wrong Laplacian calculation.

In order to get the imaging domain, we have employed statistically reinstating method (SRM) [22] on MR image, which is a level-set-based segmentation method. After the boundary of the object including electrodes is segmented by SRM, one can operate additional manual work to configure the object without electrodes or to modify the local geometry. After segmentation, one can generate the triangle mesh for the numerical computation. For the triangulation, we adopt the Triangle [23] which is a well-known open source software for the two-dimensional triangulation.

For the computation of (11) and (9), we implement the standard finite element method and apply the conjugate gradient (CG) method for matrix inversion. The size of generated triangulation is automatically determined by given geometrical information and its size is small enough to cover the pixel size. Of course, one can adjust the size of triangulation for one’s own purpose.

For the details of CoReHA 1.0, see [3].

4. CoReHA 2.0: For Better Imaging Quality

In this section, we describe in detail new features in the second version of CoReHA for dealing with in vivo MREIT experiments.

4.1. Preprocessing: Ramp Preserving Denoising

For the safety of the subject in in vivo experiments, the amount of injected current should be reduced less than a few milliamperes. The SNR of data is also directly affected by the noise of . Therefore, a proper denoising method is required for obtaining the reconstructed image with reduced noise artifact. Let us consider the reconstruction identity (7) for the derivation of denoising algorithm. Apparently, the change of the conductivity distribution is directly proportional to the Laplacian of data if the matrix is invertible. When we apply a denoising algorithm, the structure of ramp should be preserved in order to prevent the change of shape or wrong alignment of anomaly location. For our better understanding, we illustrate a simple example using a modified Shepp-Logan phantom. We compute the current density using (11) with the vertically identical conductivity distribution given in Figure 2(a). Applying the Biot-Savart law (2), we obtain the data as shown in Figure 2(b). We can see that the ramp structure of reflects the change of the conductivity distribution as shown in Figure 2(c). In order to keep the ramp structure during the denoising, we adopt a nonlinear diffusion equation method based on structure tensor [9, 24].

(a)

(b)

(c)

Figure 2: Explanation of ramp structure. (a) Conductivity distribution. (b) Corresponding data generated by the current injection along vertical direction. (c) 1D profiles along the line which is indicated by the red line in (a) and (b). Green dotted line is the conductivity distribution and the red one is the data. Locations of conductivity change match the ramp structure which is characterised by changes of slopes of the data.

For the sake of completeness, we summarise the results in [9, 24]. Let us consider the structure tensor defined by

where

To reduce the noise, we solve the following nonlinear PDE cooperating with the structure tensor :

with

where is a two-dimensional region for the slice of data, is the element of , and are maximum and minimum eigenvalues of , respectively, and are corresponding normalized eigenvectors, and is the (element by element) convolution of with the two-dimensional Gaussian kernel with a standard deviation . In the new version of CoReHA, the user only considers the parameter , the total diffusion time of the denoising algorithm. Through extensive numerical studies [9, 24], we figured out that the denoising results are robust to the tensor regularization time . The noise level of affects only . Consequently, we fix in this version, which seems to be enough time to regularize the diffusion tensor. Note that the denoising algorithm is only valid when the user considers two-dimensional Laplacian of since the diffusion process is on the -plane. For more details, see [9].

4.2. Preprocessing: Harmonic Inpainting Using MR Data

Let us revisit (6). In a local region of , where , is defected severely by the amplification of noise due to the division procedure. In the previous version of CoReHA to overcome this trouble, a manual segmentation tool was provided for the extraction of the problematic region with the harmonic inpainting algorithm [10] to solve the Poisson equation. However, the manual segmentation causes user-dependent reconstruction results. For dealing with this situation, we assume that the defected region has low MR magnitude and its conductivity distribution is homogenous so that . Then, we automatically select the defected region by taking a threshold for the MR data. Figure 3 shows an example for the selection of the defected region by thresholding. In the new version of CoReHA, we consider the of maximum magnitude of MR data as a threshold. Incorporating the ramp-preserving denoising, we additionally apply isotropic diffusion on the defected region within enough duration. The user does not need to consider how well the defected region is segmented; therefore, it guarantees the robustness of the conductivity reconstruction without user dependency.

(a)

(b)

Figure 3: Defected region select using threshold (a) MR image and (b) defected region by the thresholding (purple color).

4.3. Reconstruction: Local Harmonic Algorithm

We support a modified local harmonic algorithm [11] for the successful reconstruction by excluding the defected region directly. Let be a local region without defection which attracts our interests. Then, we locally reconstruct the conductivity by solving (9) and (10) on the local region instead of the whole domain . In order to solve (9) with (10), we use a finite element method on a selected triangle mesh and apply conjugate gradient (CG) method for matrix inversion. CoReHA 2.0 supports a user interface to draw and manipulate a polygonal local region by adding or removing points. It is useful when the user wants to see a small local change of conductivity on the local region such as brain or pelvis. In Figure 4(c), we see that the reconstructed conductivity is severely contaminated by locally magnified noise. But, Figure 4(d) shows the successful reconstruction of local brain region by the local harmonic algorithm.

(a)

(b)

(c)

(d)

Figure 4: Local harmonic algorithm: (a) MR image, (b) selected interested region, (c) conductivity reconstruction on the whole domain, and (d) conductivity reconstruction using the local harmonic algorithm.

4.4. Additional Tools

Due to the different physical properties between water and fat, the measured MR signal can be shifted in the fat region, since the MR scanner assumes that spin isochromats precess at the Larmor frequency of water. Therefore, the MR signal may become void, and overlapped images are produced in the animal/human experiments. To overcome this difficulty, new version of CoReHA contains the chemical shift correction using 3-points Dixon method [25]. In the -space data conversion, one can manipulate the MR signal using the user interface which provides separated water and fat images of MR. For more details, see [25].

Furthermore, we implemented many useful features to help our understanding of MR, , and conductivity images, for example, magnification windows and advanced alignment and resizing of windows. Figure 5 is a snapshot of CoReHA 2.0, that shows newly introduced features.

Figure 5: Snapshot of CoReHA 2.0.

5. Conclusions and Future Works

In this paper, we presented new features in a software package, called CoReHA 2.0, to deal with the data from the in vivo stage of MREIT. In this version, we enhanced the preprocessing procedure based on the ramp preserving denoising and harmonic inpainting with MR data to increase SNR of data in the in vivo experiment. Also, we provided a user interface for the local harmonic algorithm for the reconstruction on the locally interested region. The performance results of suggested algorithms and methods are given in [26]. Note that CoReHA reconstructs the contrast image of conductivity. For the medical diagnosis, the contrast image is enough. However, recently the necessity for the true conductivity came to the fore. For the true conductivity image reconstruction, we need a subject-dependent 3D segmentation and modeling tools, 3D denoising algorithm, and 3D finite element solver and higher order methods for more accurate computation.

Conflict of Interests

The authors declare that this paper has no conflict of interests with CoReHA 2.0 software.

Acknowledgment

This work was supported by NRF-2011-0015399.

References

J. K. Seo and E. J. Woo, “Magnetic resonance electrical impedance tomography (MREIT),” SIAM Review, vol. 53, no. 1, pp. 40–68, 2011.

E. J. Woo and J. K. Seo, “Magnetic resonance electrical impedance tomography (MREIT) for high-resolution conductivity imaging,” Physiological Measurement, vol. 29, no. 10, pp. R1–R26, 2008.

K. Jeon, C.-O. Lee, H. J. Kim, E. J. Woo, and J. K. Seo, “CoReHA: conductivity reconstructor using harmonic algorithms for magnetic resonance electrical impedance tomography (MREIT),” Journal of Biomedical Engineering, vol. 30, no. 5, pp. 279–287, 2009.

K. Jeon, A. S. Minhas, Y. T. Kim et al., “MREIT conductivity imaging of the postmortem canine abdomen using CoReHA,” Physiological Measurement, vol. 30, no. 9, pp. 957–966, 2009.

J. K. Seo, J. R. Yoon, E. J. Woo, and O. Kwon, “Reconstruction of conductivity and current density images using only one component of magnetic field measurements,” IEEE Transactions on Biomedical Engineering, vol. 50, no. 9, pp. 1121–1124, 2003.

W. C. Jeong, Y. T. Kim, A. S. Minhas et al., “In vivo conductivity imaging of human knee using 3 mA injection current in MREIT,” Journal of Physics, vol. 224, no. 1, Article ID 012148, 2010.

H. J. Kim, W. C. Jeong, A. S. Minhas et al., “In vivo conductivity imaging of canine male pelvis using a 3T MREIT system,” Journal of Physics, vol. 224, no. 1, Article ID 012020, 2010.

A. S. Minhas, Y. T. Kim, Y. Han, W. C. Jeong, H. J. Kim, and E. J. Woo, “Experiment performance evaluation of multi-echo ICNE pulse sequence in magnetic resonance electrical impedance tomography,” Magnetic Resonance in Medicine, vol. 66, no. 4, pp. 957–956, 2011.

C.-O. Lee, K. Jeon, S. Ahn, H. J. Kim, and E. J. Woo, “Ramp-preserving denoising for conductivity image reconstruction in magnetic resonance electrical impedance tomography,” IEEE Transactions on Biomedical Engineering, vol. 58, no. 7, pp. 2038–2050, 2011.

B. I. Lee, S. H. Lee, T. S. Kim, O. Kwon, E. J. Woo, and J. K. Seo, “Harmonic decomposition in PDE-based denoising technique for magnetic resonance electrical impedance tomography,” IEEE Transactions on Biomedical Engineering, vol. 52, no. 11, pp. 1912–1920, 2005.

J. K. Seo, S. W. Kim, S. Kim et al., “Local harmonic algorithm with domain decomposition in MREIT: computer simulation study,” IEEE Transactions on Medical Imaging, vol. 27, no. 12, pp. 1754–1761, 2008.

Y. Z. Ider and O. Birgul, “Use of magnetic field generated by the internal distribution of injected currents of electrical impedance tomography,” Elektrik, vol. 6, pp. 215–225, 1998.

O. Kwon, E. J. Woo, J. R. Yoon, and J. K. Seo, “Magnetic resonance electrical impedance tomography (MREIT): simulation study of J-substitution algorithm,” IEEE Transactions on Biomedical Engineering, vol. 49, no. 2, pp. 160–167, 2002.

E. J. Woo, S. Y. Lee, and C. W. Mun, “Impedance tomography using internal current density distribution measured by nuclear magnetic resonance,” in Proceedings of SPIE Mathematical Methods in Medical Imaging III, vol. 2299, pp. 377–385, July 1994.

N. Zhang, Electrical impedance tomography based on current density imaging [M.S. thesis], University of Toronto, Tronto, Canada, 1992.

J. J. Liu, J. K. Seo, M. Sinl, and E. J. Woo, “On the convergence of the harmonic algorithm in magnetic resonance electrical impedance tomography,” SIAM Journal on Applied Mathematics, vol. 67, no. 5, pp. 1259–1282, 2007.

J. J. Liu, J. K. Seo, and E. J. Woo, “A posteriori error estimate and convergence analysis for conductivity image reconstruction in MREIT,” SIAM Journal on Applied Mathematics, vol. 70, no. 8, pp. 2883–2903, 2010.

C. Park, O. Kwon, E. J. Woo, and J. K. Seo, “Electrical conductivity imaging using gradient decomposition algorithm in magnetic resonance electrical impedance tomography (MREIT),” IEEE Transactions on Medical Imaging, vol. 23, no. 3, pp. 388–394, 2004.

E. J. Woo, S. Y. Lee, J. K. Seo, O. Kwon, S. H. Oh, and B. I. Lee, “Conductivity images of biological tissue phantoms using a 3.0 Tesla MREIT system,” in Proceedings of the 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '04), pp. 1287–1289, San Francisco, NC, USA, September 2004.

M. Frigo and S. G. Johnson, 2008, http://www.fftw.org/

D. C. Ghiglia and M. D. Pritt, Two-Dimensional Phase Unwrapping: Theory, Algorithm and Software, Wiley Interscience, New York, NY, USA, 1998.

S. H. Park, C.-O. Lee, and J. Hahn, “Image segmentation based on the statistical variational formulation using the local region information,” Research Report 08-03, The Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea, 2008.

J. R. Shewchuk, “Triangle: a two-dimensional quality mesh generator and Delaunay triangulator,” 2005, http://www.cs.cmu.edu/~quake/triangle.html

J. Hahn and C. O. Lee, “A nonlinear structure tensor with the diffusivity matrix composed of the image gradient,” Journal of Mathematical Imaging and Vision, vol. 34, no. 2, pp. 137–151, 2009.

A. S. Minhas, Y. T. Kim, W. C. Jeong, H. J. Kim, S. Y. Lee, and E. J. Woo, “Chemical shift artifact correction in MREIT,” Journal of Biomedical Engineering, vol. 30, pp. 461–468, 2009.

K. Jeon, H. J. Kim, C. O. Lee, J. K. Seo, and E. J. Woo, “Integration of the denoising, inpainting and local harmonic algorithm for MREIT imaging of intact animals,” Physics in Medicine and Biology, vol. 55, no. 24, pp. 7541–7556, 2010.