MR-Based Electrical Conductivity Imaging Using Second-Order Total Generalized Variation Regularization

: Electrical properties provide essential information for cancer detection and specific absorption rate (SAR) estimation. Magnetic resonance electrical properties tomography (MREPT) is an approach to retrieve the distribution of electrical properties. The conventional method suffers from the locally homogeneous assumption and amplification of noise. In this study, a novel approach was introduced to improve the accuracy and the noise robustness of conductivity imaging. The proposed approach reformulated the central equation of the gradient-based method to avoid the calculation of the Laplacian operator. The equation was regularized using the second-order total generalized variation, which formulates an objective function. The optimization problem was solved by the alternating direction method of multipliers (ADMM) method. The proposed method was validated by the simulation data of the cylindrical phantom and Ella head model


Introduction
The electrical properties (EPs) of biological tissues, which include electrical conductivity σ and permittivity ε , are fundamental properties that quantify the ability to transfer electrical current inside the medium and characterize the effect of electric polarization.EPs rely on the frequency of the applied electromagnetic field and are related to multiple tissue characteristics, which involve ion concentration, volume fraction, cellular membrane permeability, and pathological condition [1][2][3].It has been reported that cancerous tissues show a significantly higher value of EPs than normal tissues according to the ex vivo experimental results; for example, the EP value of breast cancer is approximately >200% higher than normal tissue and that of bladder cancer is approximately >100% higher [4][5][6].Thus, mapping EPs has the potential for cancer detection and diagnosis.Moreover, tissue EPs are critical for the application of transcranial magnetic stimulation and transcranial direct current stimulation, which provides potential therapeutic applications for a wide variety of disorders, including Parkinson's disease, stroke, depression, and schizophrenia [7][8][9].Furthermore, the estimation of the specific absorption rate (SAR), which is a crucial problem in high and ultra-high field magnetic resonance imaging (MRI) systems, requires accurate conductivity distribution [10,11].This evidence suggests that imaging the distribution of EPs can provide essential information for both research and medical purposes.
Many efforts have been made to image in vivo EPs based on several distinct measuring techniques [12][13][14].Among them, magnetic resonance electrical properties tomography (MREPT) is a non-invasive method to simultaneously map conductivity and permittivity using the MRI system [15,16].MREPT exploits the function between EPs and electromagnetic field produced by the radiofrequency (RF) coil system, the 1 B + field.The 1 B + field can be measured by MR scanning.
Therefore, MREPT reconstructs the distribution of tissue EPs from the 1 B + measurement [17].
MREPT reconstruction algorithms can be divided into three categories, which are direct methods, inverse methods, and learning-based approaches.
The direct methods reconstruct EPs by applying Maxwell's equations to the measured RF field.The frequently used method of direct methods is the conventional method (referred to as CM); it depends on the assumption that the EPs distribution is locally homogeneous.This method has reconstructed EPs successfully [15], but it includes certain disadvantages.Firstly, the true EPs distribution in the human body, which is spatially heterogeneous, violates the above-mentioned assumption.Thus, the boundary artifacts are brought around the boundary between different tissues [18].Several studies have attempted to address this problem by introducing the gradient information of EPs to the reconstruction method [19,20].Another problem of the CM originates from obtaining phase information.MREPT reconstruction requires the magnitude and phase of 1 B + .The magnitude information can be obtained using the B1 mapping methods [21][22][23], but phase information cannot be explicitly determined using existing technology.The transceive phase assumption (TPA) is used [24].TPA approximates the 1 B + phase as half of the transceive phase, namely the sum of the transmit and receive phase.Van lier et al. found that the TPA is valid in low field strength, low permittivity, and symmetric objects [25].Lastly, the Laplacian operator is sensitive to noise and then introduces numerical errors.Lee et al. provided the theoretical relationship between the signal to noise ratio (SNR) of MREPT and the input 1 B + data [26].
Apart from calculating the EPs distribution directly from the equation derived from Maxwell's equation, inverse methods leverage a forward model for calculating B1 maps and a fitting method for evaluating the EPs distribution [27][28][29].These reconstruction approaches avoid the locally homogeneous assumption and alleviate the noise effect by avoiding the differentiation of measured data.For example, Serralles et al. proposed global Maxwell tomography on the basis of the volume integral equation formulations [30].Balidemaj et al. introduced contrast source inversion EPT (CSI-EPT) [31].However, these methods require extensive numerical efforts and have not been investigated in vivo experiments.
The recent applications of machine learning algorithms in EPT show promising results [32,33].These methods retrieve EPs from large datasets of transmit RF fields and their corresponding EPs via machine learning.Nils Hamp et al. investigated the challenges of deep learning-based methods and found that the networks trained on realistic simulation introduce artifacts when applied to in vivo data [34].
In contrast to acquiring both the conductivity and permittivity from MREPT algorithms, the conductivity reconstruction can be conducted solely by phase information based on the fact that the phase map primarily dominates the conductivity calculation [35].The transceive phase map is used to reconstruct the conductivity.Thereby, it bypasses the TPA and decreases the scan time without B1 mapping scanning.Gurler et al. introduced an electrical conductivity imaging method, which includes the gradient information of conductivity, to improve the reconstruction accuracy.However, this method involves setting the boundary value and the diffusion parameter [36].In our previous study, we proposed a double regularization method that contains total variation and wavelet terms [37].This method suppresses spurious oscillation and enhances the accuracy of reconstruction.However, the method needs the second spatial derivatives of the transceive phase, and staircase effects occur due to the total variation regularization.
In this study, we proposed a novel conductivity imaging method on the basis of the gradient-based conductivity imaging method.We reformulated the central equation through divergence properties, leading to the free measurements of the Laplacian operator.With the finite-difference method, the optimization problem was obtained, and a second-order total generalized variation (TGV) term was used to remove the staircase effect [38].The problem was solved using the alternating direction method of multipliers (ADMM) method [39].Numerical simulation and experimental were applied to validate the performance of the proposed method.The reconstruction results show that the proposed method accurately determines the distribution of conductivity and is robust in terms of noise contamination compared with the conventional and double regularization methods.

The Double Regularization Conductivity Imaging Method
The central equation used in the gradient-based electrical conductivity imaging method is as follows [36]: where μ is the free space permeability, tr ϕ is the transceive phase, and ω is the Larmor frequency.Several assumptions are applied to derive Equation (1), which assumes σ ωε ?
and neglects spatial derivatives of z B .With the first-order central differential scheme, the discretized Equation (1) in an image matrix M N L × × can be written as [37]: where Δ are the spatial resolution along the x , y , z directions, respectively.Equation ( 2) is applied to each point of the imaging region, and then an inverse problem is obtained: where A is the sparse matrix with size According to Sun [37], Equation ( 3) is constructed into an optimization problem with additional double regularization terms of total variation (TV) and wavelet transform.This method is abbreviated as "DRM" in this study.Thus, the corresponding optimization problem is written as: where W is a wavelet transform, , , , λ and 2 λ are the parameters of TV and wavelet transform, respectively.This problem is solved using the Split Bregman method [40].However, ADMM was used in this study and is introduced in detail in a later section.

The Proposed Electrical Conductivity Imaging Method
Inspired by the modified finite difference scheme [41], we can use the divergence properties, which is (5) Equation ( 5) is free of the calculation of the Laplacian operator of the transceive phase.Given the central difference method, the discretized Equation ( 5) in an image region is: tr tr tr tr Equation ( 6) can be transformed into the same linear inverse problem as Equation ( 3) inside the imaging region.The formula A b ρ= is used to represent the procedure result of the proposed method.Given that the diagonal of Matrix A is 0, the results cannot be obtained by the inverse matrix.Subsequently, the problem is expressed as an optimization problem, and second-order TGV regularization is adapted to stabilize the result and improve the accuracy of construction.The TGV regularization can preserve edges while suppressing the staircase effect and is proven as a favorable regularization for medical imaging problems [42][43][44].As previously noted [43], the discrete ( ) TGV α ρ can be written as the following equivalent form: where ( ) ) denotes symmetrized spatial derivatives.
The use of TGV as a regularization term leads to the following reconstruction model: ( ) where 1 α and 2 α are the regularization parameters.
The reconstruction problem is non-smooth and convex, and thus, it is difficult to solve through common methods such as the conjugate gradient (CG) method.The ADMM method is used to solve the optimization problem [39].The key idea of ADMM is to derive an equivalent problem with separable sub-problems by introducing auxiliary variables.The ADMM method does not depend on the smoothness of the problem and converges quickly to a reasonable accuracy result.( ) ( )

Solutions to the Optimization Problem Using ADMM
The corresponding augmented Lagrangian function of Equation ( 9) is: where 1 β , 2 β are penalty parameters, and z % , % y are Lagrange multipliers.Then, Equation ( 10) can be separated into sub-problems according to the variables, respectively.The sub-problems are as follows: ( ) % ( ) arg min 2 2 Given that Equation ( 11) is differentiable, then the optimality condition for Equation ( 15) can be calculated through the CG method or Gauss-Seidel method.The y , z sub-problems are similar, and the solutions are given explicitly by a generalized shrinkage formula [45].Specifically, the solution to the y sub-problem is: Likewise, we obtain the solution to the z sub-problem as: ( ) ( ) ( ) For the sub-problem v, we use the first-order necessary conditions for optimality as follows: ( ) ( ) ( ) Equations ( 18)-( 20) are solved by using the CG method or Gauss-Seidel method.
Correspondingly, the updated procedure of z % and % y is as follows: ( ) Based on the above-described equations, the distribution of conductivity can be calculated by inversing the ρ point-wise.

Electromagnetic Simulations
Figure 1 shows the electromagnetic simulations conducted using a Finite-Difference Time-Domain (FDTD)-based software (SEMCAD X. 14.6, ZMT, Zurich, Switzerland).A high-pass shielded birdcage coil model was used for the numerical simulation.The length of the leg was 300 mm.The diameter of the RF shield and rings were 330 and 300 mm, respectively.By adjusting the capacitor in the rings, the quadrature coil was tuned at 128 MHz (3T) and was driven in a quadrature mode by two ports separated by 90°.The simulation model included a cylindrical phantom and Ella head model (Ella Model, the Virtual family [46]).The diameter and height of the cylindrical phantom were both 160 mm. Figure 1b illustrates that the phantom contained six same size cylinders with small diameters; different EP values were assigned to these cylinders.Figure 1c shows the Ella head model, and the conductivity value of the head tissues was appointed, according to Gabriel et al. [3].The analysis of reconstruction results included three primary tissues of the head, namely gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF).
After the phase information of 1 B + and 1 B − was obtained from the electromagnetic simulation, these data were added to acquire the transceive phase information.Given that the simulation grid was not uniform, the non-uniformly spaced grid was re-gridded into a regular grid by using a linear interpolation method to allow a comparison on the voxel level.The proposed method was validated for noisy simulated data.The noise distribution in MR phase images is assumed as Gaussian noise.The standard deviation (SD) of Gaussian noise in phase images was ( ) 2 SN R , in which SNR is the signal-to-noise ratio in magnetic resonance (MR) magnitude images [47,48].In the simulations, the different MR magnitude SNR values, such as 200, 100, and 50, were set to alter the SD of the Gaussian noise.Then, the noise was added to the transceive phase image to simulate the noisy phase image [36].For the simulations, all methods were repeated for 100 implementations at each noise level to calculate the root mean square error (RMSE) over the reconstructed conductivity profile; the RMSE equation was: where i r is the reconstructed conductivity of pixel i and it r is the target value.

Experimental Setup
A cylindrical phantom was used to verify the proposed approach.The structure of the cylindrical phantom contained inner and outer compartments that were separated by a tube.The materials of the phantom were distilled water, NaCl, agarose, and CuSO4•5H2O.Agarose was used to solidify the phantom, whereas CuSO4•5H2O was used to decrease the T1 values and then reduce the scanning time.The contrast of the phantom was prepared by adjusting the concentrations of saline.The mass ratio of distilled water, agarose, and CuSO4•5H2O was 100:1:0.1.The saline concentrations were 1.5% and 0.5% for the inner and outer compartments, respectively.Conductivity values were measured a by network analyzer (Agilent 4395A, Santa Clara, CA, USA) associated with an S-parameter test set (Agilent 87511A).The conductivity values for the inner and outer compartments were 2.01 and 0.97 s/m, respectively.
The in vivo experiment included a healthy 26-year-old male volunteer.The subject gave their informed consent for inclusion before they participated in the study.The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of Southern Medical University (2018B030333001).
The phantom and in vivo experiments were carried out on a 3.0T MRI scanner (Philips Achieva, Best, The Netherlands) installed in NanFang Hospital.We used the standard body coil as transmitting and receiving and applied the spin-echo (SE) sequence to obtain the transceive phase information.The 2D SE sequence parameters in the phantom experiment were TR/TE = 1000/10 ms, with 5 averages in 20 min.The acquisitions were performed in a total of five slices with a resolution of 1.17 × 1.17 × 3 mm 3 and a field of view (FOV) of 150 ×150 mm 2 .For the in vivo experiment, the parameters of the 2D SE sequence were TR/TE = 1200/7 ms; 10 slices centered on the head center were scanned, with a resolution of 1.88 × 1.88 × 5 mm 3 and a FOV of 240 × 240 mm 2 .This scan took 3 min.

Image Reconstruction
The accuracy of reconstruction was enhanced using a Gaussian kernel to filter the noisy simulation data and experimental data, the kernel size of the filter was 7 × 7 × 3 voxels, and the standard deviation was 1.2 × 1.2 × 1.2 pixels.The proposed method was compared with the CM and the double regularization method (DRM).All calculations were conducted in MATLAB (The MathWorks, Natick, MA, USA).
For the DRM and proposed method, several parameters need to be set.The RMSE was used to estimate the performance of reconstruction for different parameters, and then the proper parameters were selected.Table 1 presents the actual regularization parameters.
Table 1.Parameters of the double regularization method (DRM) and the proposed method.

DRM
Proposed Method

Phantom Simulation Results
The reconstructed conductivity images of the cylindrical phantom were compared with the corresponding target image, as shown in Figure 2.Under the noise-free condition, the conductivity distribution in the homogeneous region is accurate for all the methods in Figure 2a.The distinct difference between these methods depends on the boundary artifacts in the transition region between different contrast components.The boundary artifacts are obviously around the small cylindrical for the CM, with origins from the locally homogeneous assumption.In the DRM conductivity image, the artifacts are reduced due to the adoption of gradient information in the critical equation.However, several artifacts remain around the transition region because of the Laplacian operator.By contrast, the proposed method does not require the calculation of the Laplacian operator, yielding a more precise result than the DRM. Figure 2b depicts the conductivity profiles along the red line in Figure 2a.The proposed method represents a smoother result that is close to the target values.Figure 3 shows the mean and SD of the reconstructed conductivity in each region for different methods.In most areas, the proposed method leads to accurate values and lower SD values than other methods.Figure 4 illustrates the conductivity images under different noise conditions by using three reconstruction methods.The CM conductivity results are severely affected by noise because the Laplacian operator is sensitive to noise.The use of a larger Laplacian kernel can alleviate this phenomenon, but it can also over smoothen the results.The DRM results demonstrate better suppression of the noise effect.In contrast to DRM, the proposed method shows superior conductivity images.The reconstructed results display a clear phantom structure, and the small inner regions retain a sharper edge than the other methods as the noise increase.Thus, the proposed method can alleviate noise effects.The RMSE was used to assess the accuracy of three reconstruction methods with various noise levels.Figure 5 shows the results; the x-axis presents the SD of zero-mean Gaussian noise, in which the SD is ( ) 2 S N R ; the SNR is the signal to noise ratio in the MR magnitude image.Notably, the RMSE of the CM increases with noise and is much larger than that of other methods.For DRM and the proposed method, the change of RMSE in low noise stays small, and the DRM RMSE is larger than the proposed method.These two methods similarly show an increasing tendency of RMSE as the noise increases.The magnitude of the DRM is larger than that of the proposed method.The RMSE of the proposed method stays at a low level, even in severe noise corruption.Thus, the proposed method is robust for noisy phase data.
Figure 5.The root mean square error (RMSE) of reconstructed conductivity using the conventional method (CM), double regularization method (DRM), and the proposed method (PM).

Ella Head Simulation Results
Figure 6 depicts the reconstruction of the Ella head model simulation.Conductivity images are shown from different head positions under noise-free conditions.All the results display the basic structure of the head, which includes WM and GM.Similar to the phantom simulation, the CM conductivity images yield artifacts in the boundary, especially around the CSF tissue.For the DRM, the reconstructed conductivity shows precise results compared with target images.However, the dispersal in the CSF boundary emerges evidently in positions 2 and 3. Compared with DRM, the proposed approach shows a more accurate structure, and the edge is sharper and precise.The results of all methods get worse as the noise increase because the structure of the head is more complicated than the phantom.In addition, the larger Laplacian kernel and the first-order spatial derivative operator dismiss several small structures.The proposed method preserves more details than the results of DRM, such as the WM region in position 2. Figure 8 depicts the RMSE of the reconstructed conductivity images under noise conditions.The RMSE increases as the noise increases and the proposed method shows the lowest values, which display the advantage of the proposed method.

Experimental Results
Figure 9 describes the reconstructed conductivity image of the phantom experiment.Figure 9a and 9b represent the SE magnitude and phase images, respectively.The SNR of magnitude is 132.In the magnitude image, the region between the blue circles is region 1, while the red circle describes region 2. The black circle in the magnitude image is the material used to separate these two compartments.Figure 9b illustrates the conductivity using the CM, DRM, and proposed method.The conductivity of the proposed method exhibits more homogeneity in two regions than CM and DRM.Table 2 shows the calculated mean and SD of these two regions.The proposed method has more accurate mean values than CM and DRM, indicating its feasibility to reconstruct accurate conductivity.

Discussion
In the present study, a novel approach was proposed for imaging the conductivity distribution from the MR transceive phase.On the basis of the gradient-based method, the proposed method reformulates the central equation through divergence properties and then leverages the finite difference method to construct an inverse problem.The resultant optimization problem is regularized by TGV and solved by ADMM.The proposed method does not require the Laplacian operator of the transceive phase, which is the origin of the dominant errors of general approaches.A systematic investigation using numerical simulations and experiments at 3T illustrate that the proposed method can improve the reconstruction accuracy and is more robust to the noise than the CM and DRM.
The CM, which relies on the locally homogeneous assumption, is applicable for those objects with a simple structure but not for complicated objects such as the brain, as shown in Figure 5.The pre-segmentation of objects can partially mitigate these adverse effects before reconstruction.However, this method introduces extra work, and the assumption remains.Therefore, it cannot solve the problem fundamentally.By contrast, the proposed method uses conductivity gradient information and thus in principle is immune to the effects of the above-mentioned assumption.The Ella model simulation results indicate that the proposed approach is favorable for objects with complicated structures.
The CM and DRM require the calculation of the Laplacian operator, which is sensitive to noise.The results of the CM show the effect of noise, while those of the DRM show the impact of the Laplacian operator despite applying conductivity gradient information.A noise-robust Laplacian kernel was used to alleviate the effect of noise [49].However, despite the noise suppression, this leads to over-smoothened reconstruction results.Another method reconstructs conductivity profiles by using the first-order differencing of the RF filed data [50].However, the induced current density image must be determined to obtain the EPs.By contrast, the proposed method transforms the equation to remove the necessity of the Laplacian operator of the transceive phase and directly retrieve the conductivity.The simulation and experiments show that the proposed method decreases the impact of noise.
The proposed method depends on the central equation of the gradient-based method [36].The following three assumptions are made.The first assumption is that the spatial derivatives of 1 B + and 1 B − are negligible.This assumption is valid in the quadrature birdcage coil central region when the coil is used to transmit and receive.However, this assumption appears untrue at high field strengths (>3T).The RF shimming technique may provide a resolution to improve the homogeneity of the RF field.The second assumption is ( ) , which is sufficient for common tissues at field strengths ≤ 3T.The reconstructed conductivity value shows overestimation compared to the real value when violating the assumption, such as region 2 of the phantom simulation.This assumption is close to the assumption used in the conventional approach, which is σ ωε >> [25].Third, the spatial derivatives of z B are negligible compare to 1 B + and 1 B − , which appear rational in the central region of the birdcage coil from the simulation.From the perspective of clinical application, the z B information cannot be directly acquired.
The optimization problem of the proposed method includes a second-order TGV regularization term that can suppress the noise while omitting the staircase effect.ADMM was applied to address the problem and then obtain accurate conductivity images.Compared with DRM, the proposed method eliminates the staircase effect and shows better results in terms of noise suppression.The implementation of the proposed method requires four parameters, which are determined by using RMSE as a criterion.The optimal parameters for different regions vary.Thus, selecting spatially dependent parameters provides the potential to obtain a more accurate result than before, as observed in other applications [51].During the calculation, there are a few iterations to acquire moderate conductivity results, indicating that the time cost is short.The phantom and brain reconstructions take approximately 2 s and 10 s, respectively, on a standard PC with an Intel i5 (2.9 GHz) and 16 Gb of RAM.
The proposed method reformulates the central equation to remove the Laplacian operator through divergence properties.In the discretization for the imaging region, different points are used to construct an inverse problem compared with DRM.Inspired by the recent research [52], we can also use the divergence theorem to transform the equation into an integral formula that can suppress noise, further improving the reconstruction The distribution of conductivity is reconstructed using only transceive phase information.Thus, any sequences can be used to obtain the transceive phase, theoretically.However, other factors, such as the homogeneity of the main field strength, eddy currents, and the phase contribution from flow and motion contaminate the phase data [16].In this study, SE was chosen to acquire the phase.In other studies, the steady-state free precession (SSFP) sequence shows advantages for scanning because of its short time [53].During the present experiment, a volume coil was used for transmitting and receiving, which is a typical configuration in clinical application.The multi-channel coil as the receive coil improves the SNR of the MR image and saves time [54].It would be exploring in future research.In addition, when the phase image is achieved, the magnitude image is also obtained.The magnitude image includes the structural information of the object, the incorporation of which should be explored to get better results.

3 D
are the finite difference operators along the x-axis, y-axis, and z-axis, respectively; v represents an approximation of ρ ∇ and is a matrix field with the same size of ρ ∇ ; and ( ) ( )(

First
and then Equation (8) is transformed into a constrained problem:

Figure 1 .
Figure 1.Simulation setup.(a) The birdcage coil loaded with a numerical phantom, (b) seven regions of phantom, (c) the birdcage coil loaded with the Ella model.The radiofrequency (RF) shield is not shown in (c).

Figure 2 .
Figure 2. Comparison of reconstructed conductivity images using different methods under the noise-free condition.(a) The target image and reconstructed conductivity images using the conventional method (CM), double regularization method (DRM), and proposed method (PM).(b) The reconstructed conductivity profiles along the red line in the target image in (a).

Figure 3 .
Figure 3. Bar graph of reconstructed conductivity value using conventional method (CM), double regularization method (DRM), and proposed method (PM).It corresponds to seven regions.

Figure 4 .
Figure 4. Reconstructed conductivity images under different noisy transceive phase data using the conventional method (CM), double regularization method (DRM), and the proposed method (PM).

Figure 6 .
Figure 6.The reconstructed conductivity images of the Ella head model under noise-free conditions using the conventional method (CM), double regularization method (DRM), and proposed method (PM).

Figure 7
Figure 7 describes the head model results under various noise conditions, which are SNR = 200 and SNR = 100.The results of all methods get worse as the noise increase because the structure of the head is more complicated than the phantom.In addition, the larger Laplacian kernel and the first-order spatial derivative operator dismiss several small structures.The proposed method preserves more details than the results of DRM, such as the WM region in position 2. Figure8depicts the RMSE of the reconstructed conductivity images under noise conditions.The RMSE increases as the noise increases and the proposed method shows the lowest values, which display the advantage of the proposed method.

Figure 7 .
Figure 7. Reconstructed conductivity images of the Ella head model using the conventional method (CM), double regularization method (DRM) and proposed method (PM) under a signal to noise ratio (SNR) = 200 and SNR = 100.

Figure 8 .
Figure 8.The RMSE of reconstructed conductivity images under noise conditions using the conventional method (CM), double regularization method (DRM), and proposed method (PM).

Figure 9 .
Figure 9. Phantom reconstruction result.(a) Magnitude and phase images of the phantom; the SNR of the magnitude image is 132.(b) The reconstructed conductivity images using the conventional method (CM), double regularization method (DRM), and the proposed method (PM).

Figure 10
Figure 10 shows the magnitude, phase, and reconstructed conductivity images of the in vivo experiment.The SNR of the magnitude image is 75.The proposed method shows a clearer structure than the results of the CM and DRM.The CSF region offers more errors, which are possibly caused by low-phase data.

Figure 10 .
Figure 10.Magnitude, phase, and reconstructed conductivity images of head experiment.(a) Spin-echo (SE) magnitude and phase images; the SNR of the magnitude image is 75.(b) Conductivity images using the conventional method (CM), the double regularization method (DRM), and the proposed method (PM).

Table 2 .
Reconstructed conductivity values of phantom experiment