A Regularization Process for Electrical Impedance Equation Employing Pseudoanalytic Function Theory

The electrical impedance equation is considered an ill-posed problem where the solution to the forward problem is more easy to achieve than the inverse problem.Thiswork tries to improve convergence in the forward problemmethod,where the Pseudoanalytic FunctionTheory by means of the Taylor series in formal powers is used, incorporating a regularization method to make a solution more stable and to obtain better convergence. In addition, we include a comparison between the designed algorithms that perform proposed method with and without a regularization process and the autoadjustment parameter for this regularization process.


Introduction
The electrical impedance tomography (EIT) problem, which is set by employing the electrical impedance equation, was mathematically posed by Calderon [1] in 1980 proving that the solution of this problem exists as being unique and steady.This problem is exemplified by the following equation: where  represents the conductivity and  denotes the electric potential for a domain Ω within a boundary Γ.This equation is considered ill-posed problem, due to its high complexity to find a solution in terms of the initial data; it means that every variation in the input data can approximate in different way the solution to this problem.Therefore, if a measurement presents any variation or noise, this variation presents difficulty to obtain good approximation; for example, when the finite element method is employed in noise presence in the electrodes, we need to introduce an extra process to correct this fault or the approximation could not be reached.Some methods should be employed to correct these perturbations present advantages and disadvantages, such as a regularization method that permits to suppress the noise in the measurements.Due to this problem and the continuous advances in the field, the EIT is not considered a practical medical imaging procedure.
Employing (1), we can approach two types of problems; in the first one, the conductivity  value is known and the task is to approximate the electric potential in the boundary | Γ .This problem is known as forward Dirichlet boundary value problem; it is also known as nonhomogeneous Laplace equation.For the second problem, the value of the electric potential in the boundary | Γ is known, but the conductivity  within the domain Ω is unknown; this problem is formulated as inverse Dirichlet boundary value problem or electrical impedance tomography problem [2].
The EIT was considered very unsteady, as well as ill-posed problem [2], but, in 2009, Kravchenko [3] and, independently in 2006, Astala and Päivärinta [4] noticed for first time that the two-dimensional case was completely equivalent to special case of Vekua equation [5].
There exist several assortment methods that try to solve forward and inverse problem for (1), where the best mathematical tool to approximate the solution for both problems by now is the finite element method (FEM).It proves to be 2 Mathematical Problems in Engineering stable and easy to employ, but it presents a difficulty when it computes the approximation; this difficulty is presented in the initial data and the variations on it [2].
The existence of different techniques and methods to reach a solution of inverse problems, in which the solutions are interpreted in terms of linear operators [6], permits analysing and studying inverse problems in order to find a solution.Other works employ the connection between pseudoholomorphic functions and conjugate Beltrami equations to deduce the well-posedness on smooth domains of the Dirichlet problem for 2D isotropic conductivity equations [7].There exists a theory of conjugate functions to solve Dirichlet and Neumann problems for conductivity equations, in which they consider some density properties to trace solutions with a boundary approximation issues [8].
All these works contribute to analysing and reaching a solution to the forward problem and continuing working on the solution to the inverse problem of the electrical conductivity equation.Employing all these techniques and analysis the study of the forward problem could help to find a solution to the inverse problem, but the study to be done should be extensive and incorporate different techniques and methods proposed by the different authors.
For the correct understanding of the inverse problem, first, we need to analyse and study the forward problem, to determinate how the energy propagates within the domain; with the collected data we can determine that the convergence and the performance will be very important to design a method that can be employed in the medical imaging in a future.
In this study, the main purpose is to analyse the behaviour of the algorithm based on the Pseudoanalytic Function Theory, which applies the Taylor series in formal powers as the support, proposing additionally to employ a regularization process to improve the stability.The results obtained are compared with the results presented in previous paper [9,10], in order to analyse the improvement or decrement in the convergence due to the method designed.
The remainder of this paper is organized as follows.Section 2 shows the mathematical tools for the electrical impedance equation; in Section 3, the methodology that contains the main idea of the algorithm is presented; following, a brief review of regularization procedure employed to design this algorithm is explained; Sections 4 and 5 show several experiments to compare the results and a discussion about this comparison; finally Section 6 expresses the conclusions about the behaviour of the novel method and discusses a future work.
Then, any complex function  can be represented by the linear conjugation of  and : where  and  are indeed real valued functions.Therefore, the pair (, ) is called generating pair.Following the Pseudoanalytic Function Theory posed by Bers [11], it is possible to introduce the derivative and antiderivative form of a complex valued function , which can be reviewed in the mentioned work.
Let us suppose (, ) generating pair with the form where  is a nonvanishing function within a domain Ω(R 2 ).
Then considering this (, ) generating pair with a  separable variable function within the domain, Thereby (, ) is embedded into a periodic sequence, in which for an  even the generating pair are and for the odd generating pair the forms are Consider the formal powers  0  ( 0 ,  0 ; ), associated with a (  ,   ) generating pair, with the formal degree 0, complex constant coefficient  0 , and center  0 , and depending upon  = +, are defined in agreement with the expression where  and  are complex constants which fulfills the condition Now, let us suppose  to be a (  ,   )-pseudoanalytic function.Thus, we can express it in terms of the so-called Taylor series in formal powers: Since every (, )-pseudoanalytic function  accepts this expansion, the last equation is an analytic representation of the general solution for the Vekua equation.Consider the two-dimensional case of the electric impedance equation (1), and suppose that the conductivity  function can be expressed in terms of a separable variable function; it follows that Introducing the notations where   = /  and   = /  , then (1) turns into special case of Vekua equation with the form where   =   −   and  2 = −1 for the standard imaginary unit.
Employing the last statement presented in (13), we have to introduce some expressions that were presented in the Pseudoanalytic Function Theory [3] and in assortment works [9,10,12], explaining the mathematical basis of the Taylor series in formal powers: (  ,  0 ; ) ; Equation ( 14) was proved in [11], showing that any (  ,   )-pseudoanalytic function can be represented in Taylor series in formal powers.So, from this point of view, ( 14) is an analytical representation of the general solution of the Vekua equation (13).
In (14), we used the Pseudoanalytic Function Theory to approximate the special case of Vekua equation in Taylor series in formal powers, where we can read in the last statement that it has a center in   , depending upon  and approximating the coefficient   , like a Taylor series.
For this case, we also remember that Taylor series in formal powers are performed by introducing a generating pair (  ,   ) and the real constants   and   that fulfill the condition presented before in (9) for the expression (8).
The admittance of the integral was also introduced in [11], taking place when the next notation is written according the recursive formula: where ( 15) is employed to approximate higher exponents of the formal powers.
Employing the structure of (8), (9), and (15), the special case of Vekua equation presented in (13) can be approached by the next statement introducing the factors from ( 14) to (16), obtaining the next equation: Equation ( 16) is the definition of a (, )-integral of a complex valued function .Specifically, since (  ,   )integral of the (  ,   )-derivative of  reaches the integral expression (17) can be considered the ( 0 ,  0 )antiderivative of the function  ( 0 , 0 ) .
For the correct approximation of the electric potential | Γ for (1), using the Taylor series in formal powers, let us note that this electric potential can be found by the next statements: In this case, if the boundary condition | Γ is provided by (19), we can always approximate asymptotically the experimental electric potential | app by the next expression: where   and   are real numbers.This procedure has proven its effectiveness in assortment works [9,10,12], where the numerical approximation achieved highly accurate results.
The analysis of the problem can include samples in piecewise and nonpiecewise separable variable functions for the conductivity ; to understand this situation we have to employ the conjecture exposed in [10].

Conjecture 1. Suppose a 𝜎 arbitrary conductivity function defined within a bounded domain Ω(R 2 ). This function can be approximated by means of a piecewise separable variables function in the form
here  is a real constant such that  +  ̸ = 0 :  ∈ Ω(R 2 ) and {  }  =1 are the constructed interpolation functions with a finite number of samples M of the conductivity function , valued in the -axis that is parallel line within the subdomain Ω, created by tracing {}  =0 parallel lines.
These piecewise separable variable functions can be employed for numerically approximating the set of formal powers.

Proposition 2. Consider an arbitrary conductivity function 𝜎 defined within a domain Ω(R 2 ). It can be considered the limiting case of a piecewise separable variables function, expressed in the form of the Conjecture (21), when 𝐾 subdomains and M number of samples are every subdomain tending to infinity:
A regularization process is needed to increase the precision of the algorithm and obtain a better solution to be analysed for further research in the inverse problem; in the next section we introduce a simple idea to employ a regularization process in the Taylor series in formal powers and iterative method using this idea to improve the results, or at least to obtain a better convergence.

Regularization Process.
Regularization methods are employed in many fields, introducing additional information with the purpose to solve ill-posed problems.Also, they are employed to prevent overfitting and are usually presented in form of penalty for complexity, such as restrictions for smoothness or bounds in the space norm for a vector.In particular case, the including of these methods to the algorithm designed to compute a solution for inverse problems is essential to fitting the data in order to reduce a norm of the solution.
For current problem, the imposed electric potential | Γ represents a finite number of current densities   |  =0 for a domain Ω with a boundary Γ and the knowledge of the Neumann-Dirichlet map [13] where Remember the notations introduced before in (1); we know that  represents the conductivity and  denotes the electric potential in the electric conductivity equation, where  Γ  = / Γ , being all the values of the electric potential  in the boundary domain Γ, and  is the current density.
For the forward Dirichlet boundary problem, we know the value of the conductivity  but the electric potential | Γ is unknown.In this work, the electric potential | Γ is imposed, and, for the purpose of this paper, the approximation of an experimental value  is computed comparing the results by means of the Lebesgue measure and determining if the method employed is being corrected [14].
The Tikhonov regularization process can be presented in such form where ‖ ⋅ ‖ denotes the  2 norm,  can be a square matrix or a matrix of  ×  dimension,  denotes a vector with  rows,  is the regularization parameter, and the problem turns into a minimization of the parameter .
Algorithm 1, which we propose in this study, is based on the Pseudoanalytic Function Theory, employing the Taylor series in formal powers, and it approximates the forward Dirichlet boundary value problem for (1), by means of  total number of radii,  maximum number of formal powers, and  total number of points.Thus, (24) has an explicit solution expressed as follows: where  is the identity matrix with the same dimension of matrix , the regularization parameter  > 0, in this case  =  (+1) [] represents the orthonormalized system, which is the result of the formal powers employing the conductivity  and the imposed boundary condition, and  is the electric potential | Γ .The main problem is to choose the correct parameter , which approximates the solution for (25).
To choose a parameter , the next expression should be used: where , as it was mentioned before, is the orthonormalized Taylor series in formal powers and  represents the electric potential [15].We employ this regularization process to understand the behaviour of the problem leaving for the next section the methodology of the main process for approximation of the solution.

Procedure
The main goal is to develop an effective stable method, which let us to obtain better approximations for the forward problem by the employment of a regularization procedure for better convergence.The results obtained will be analysed in order to understand the problem presented in (1).For this analysis, to find the forward Dirichlet boundary value problem using the methodology shown before, we use the algorithm from [9,10]; as a result we obtain an accurate solution with a lower computer cost.
The principal difficulty, when the algorithm is being employed, lies in the instability of the method when the approximations are taking place.So, in order to improve the convergence of the result, a regularization process is included in the approximation process.
(5) Perform the computation of the electric potential.
First, a domain should be chosen, employing Algorithm 1.According to [9], there can be used smooth and nonsmooth domains, below; to analyse the results a unitary disk domain is chosen.Then, a conductivity function has to be chosen, and, in this study, the conductivity function is taken from the mathematical analysis and geometrical distribution; both cases are within the domain Ω.
Once the domain Ω and the conductivity function sigma are selected, the process continues by selecting the  maximum number of formal powers,  number of radii.and  number of points per radius, the approximation of the Taylor series in formal powers is performed, and the result obtained from this process should pass through the algorithm to the Gram-Schmidt orthonormalization method.
In this part of the procedure, the boundary solution can be computed, but we propose to modernize the procedure using a regularization method, in order to obtain better convergence of the method.Figure 1(a) shows graphically the procedure explained before.
The main idea is to confirm that novel regularized algorithm is able to compute more stable solution for the electric potential in the boundary and to express numerically the convergence of the method.
Additionally, we propose a modification of the explained algorithm where an autoadjustment should be performed before the error estimation process.Here, we use iterative scheme that can estimate the regularization parameter improving the convergence property of the algorithm.The modification made to Algorithm 1 presents a condition that should be satisfied, presenting an iterative method to approximate the solution of (1).In this algorithm, the condition is the same as presented before in (26), where the   in the iterative algorithm is the approximation of the electric potential  app , as it is exposed in Algorithm 2. This autoadjustment introduced in the algorithm computes the solution of the forward problem increasing the computational cost but permits obtaining a better solution that is more stable and convergent.
Applying Algorithms 1 and 2 shown above, we perform the numerical experiments using mathematical expressions, such as exponential and sinusoidal function, and geometrical distributions such as disk center, and five disks' structure at the center.
This set of examples gives us a full perspective of the behaviour of the method, because we employed some samples to analyse the convergence of the method and then make experiments with the geometric distributions to analyse the results that are reviewed in the next section.
The examples investigated in the next section are designed to compare the results when a regularization method is used in the process versus the actual method, which does not employ a regularization procedure.

Results
As it is exposed in the works [9,10], the possibility to employ geometrical distributions functions and mathematical expression in the algorithm, which use the Taylor series in formal powers, let us analyse the electrical impedance equation ( 1), emphasizing the possibility to approximate the forward problem to this equation.

Tikhonov regularization
Error estimation (c) Tikhonov and error estimation blocks Applying this analysis, we used an exponential and sinusoidal function coming from the mathematical analysis, circle at center, and the five disks' structure at center geometrical distributions.These cases are computed within the domain, and the results of the approximation are passed to regularization process in order to analyse the approximation and stability of the method, obtaining the electric potential in the boundary.
The solution for error, in order to analyse behaviour of the solution, is used in form of the Lebesgue measure where | Γ is the imposed boundary condition and  app denotes the approximation by using the original algorithm  and the new modification employing the regularization process of everything within the unitary disk.We show the results in tables to compare and analyse the results obtained by these methods.The next sections show the results for analytical and geometrical cases, in which the analytical cases are representatives due to the consideration of nonseparable conductivity functions in which the results can be achieved; for the geometrical case, the conductivity employed does not possess a boundary condition, imposing an artificial condition to compare the results and determine the accuracy of the method.

Exponential Conductivity Function. Let us use conductivity function, which fulfills (12) and possesses the form
and this expression is shown in Figure 3(a), and for the boundary condition the expression to be imposed is where  denotes the coefficient employed to change the behaviour of the function employed.Table 1 expresses the behaviour when the  coefficient increases and demonstrates the stability property of the method in two cases when the regularization method is employed and when this regularization process is not used.
The results in Table 1 show the known behaviour; meanwhile the number of formal powers increases and the error decreases, demonstrating a better convergence when the regularization method is employed.
The comparison between two algorithms that use the regularization procedure with and without the autoadjustment algorithms has been proved, the convergence and stability in the autoadjustment procedure being better than for comparison methods, but the computational cost considerably increases due to the iterative procedure performed.
This comparison is presented in Figure 3(b); in this graphic the three methods show the behaviour of the expression and in this case all methods reached a good approximation; since the three methods possess a low error, the difference can not be appreciated in the graphic matter, but Table 1 resumes correctly the full error analysis between all the methods.

Sinusoidal Conductivity Function. Let us employ a conductivity function, with the form
Figure 4(a) shows the expression within the unit disk, and, imposing a boundary condition, where  denotes the coefficient used to change the behaviour of the function employed.
Table 2 shows the behaviour of the method when the  coefficient increases its value, including the simulation with and without the regularization process.
The results in this case are better with a regularization procedure; due to the conductivity that presents so many vibrations in the domain, the algorithm tries to reach a solution but its behaviour is bad; differently, when the regularization process is introduced, the results expose better performance because the perturbations or vibrations within the domain are smoothed and it is more easy to compute an approximation to (1).
In comparison with the other methods, the behaviour of the proposed algorithms is similar: using more formal powers the error is decreased and by the employing a regularization method the stability and convergence are improved considerably, but due to the iterative method used the computational cost is increased.
In Figure 3(b), the behaviour of all methods is shown; in this case the approximation can be appreciated because of the values of its error; for this case, the error can be analysed in  terms of the method employed, proving that a regularization method is needed to approximate and reach a better solution.

Circle at
Center.For the next case, we employed a disk at the center with conductivity  1 = 100 and the conductivity in the rest of the domain Ω is  2 = 10, as is shown in Figure 5(a).For this case, we imposed the following boundary condition: The results are shown in Table 3, where the behaviour of the method with and without the regularization process is presented.
In this case, the figure inside the domain is a circle with a variable radius (in concrete case, radius is 0.5); we know that if the disk within the domain is bigger, the approximation will be good, but if the inner conductivity is smaller, the algorithm could not detect the conductivity good; differently, when a regularization process is introduced, it warrants that the conductivity will be employed to compute the solution and the regularization process permits smoothing the conductivity to reach results.
This example provides unexpected results, because the convergence of the NPSM-Tikhonov and NPSM-Tikhonov with autoadjustment process does not obtain the convergence that the NPSM method can obtain.This case presents the best convergence and stability obtained by the NPSM algorithm, and the method with regularization process appears to obtain worse convergence but is still acceptable due to the error estimated.
This example gives us different ideas to improve and to analyse the problem for the optimized method.The first idea about the bad convergence in the solution is that the imposed  boundary condition is not the correct condition for this case; another idea is to employ a different approximation process based on mesh and finally to determine if the regularization process can improve; a better solution can be imposed in the boundary when it comes from the approximation in the same method, which comes from the Taylor series.Figure 5(b) shows the behaviour of the solution in terms of the error computed; in this case, due to the low error, the difference in the graphic cannot be observed, but the difference exists and can be appreciated in the analysis of Table 3.For this case, the error increases a little, but it possesses a steady behaviour.

Concentric Disk at
Center.In this experiment, we used a disk with four rings at the center with different radius and conductivities:  1 = 10 for the disk,  2 = 15 for the first ring,  3 = 20 for the second ring,  4 = 30 for the third ring, and  = 100 for the last ring, as is shown in Figure 6(a).
For this case, we imposed the boundary condition (32).The simulation results are shown in Table 4, presenting the behaviour of the method with and without the regularization process.
This case is difficult due to the disk structure inside the domain; both methods proved to be good when they approximate the electric potentials, in which the algorithm without regularization presents better convergence, compared to the one which used a regularization process presenting smoothing within the domain to find an approximation with more stability.
The best approximation is obtained for the NPSM algorithm, and the NPSM-Tikhonov with and without autoadjustment presents worse approximation but is useful for the problem to be solved.The behaviours of the methods are similar: a high number of formal powers give a better convergence, the number of radii and points per radius do not affect the convergence, and the regularization process  with and without an autoadjustment method only affects the convergence and stability in the solution.In this case, the approximation that is obtained by these methods is considered good, proving that the regularization method can be used for this problem.
This example is useful to understand the behaviour and, thanks to this analysis, the main idea is preserved; a complex study needs to be performed paying more attention to the boundary condition and employing the results obtained in the Taylor series instead of employing the imposed boundary condition in order to determine if the boundary condition which is employed is the best condition to fulfill the domain and the conductivity within; this analysis could be demonstrated if the inverse problem can be approached by means of the Taylor series in formal powers with and without a different interpolation procedure including the regularization process.
Figure 6(b) shows the error analysis of the example employed, in which the behaviour is the same; the method proves to be steady and to reach the solution in terms of the error computed.For this case, the solution shown in the graphic does not express a considerable difference and fits correctly, but in terms Table 4 shows that the error increases and decreases depending upon the data introduced in the problem.

Discussion
The examples exposed in Tables 1, 2, 3, and 4 illustrate the behaviour of the algorithms with and without the regularization process.In the tables, the columns with the title E NPSM, proved this behavior; meanwhile the formal powers increase, the convergence of the approximation is better, and the number of points and radii do not affect considerably the approximation in this method.
Then, the regularization process, which is employed in all cases, showed better stability and in some examples a better convergence; this is the proof that a regularization process could be employed in order to approximate forward problems when the conductivity function presents a nonsmooth conductivity function, giving the possibility to approximate the forward Dirichlet boundary value problem for (1) when the regularization is employed.
All the cases analyzed and discussed in the last section show better convergence and stability.The reasons for these approximations are due to the conductivity employed in

Figure 2
Figure2illustrates the diagram of the algorithm proposed.The procedure analyses the convergence and the stability when an automatic adjustment takes place.The modification made to Algorithm 1 presents a condition that should be satisfied, presenting an iterative method to approximate the solution of(1).In this algorithm, the condition is the same as presented before in (26), where the   in the iterative algorithm is the approximation of the electric potential  app , as it is exposed in Algorithm 2. This autoadjustment introduced in the algorithm computes the solution of the forward problem increasing the computational cost but permits obtaining a better solution that is more stable and convergent.Applying Algorithms 1 and 2 shown above, we perform the numerical experiments using mathematical expressions, such as exponential and sinusoidal function, and geometrical distributions such as disk center, and five disks' structure at the center.This set of examples gives us a full perspective of the behaviour of the method, because we employed some samples to analyse the convergence of the method and then make experiments with the geometric distributions to analyse the results that are reviewed in the next section.The examples investigated in the next section are designed to compare the results when a regularization method is used in the process versus the actual method, which does not employ a regularization procedure.

Figure 1 :
Figure 1: Methodology to compute the forward problem approximation.

Figure 2 :
Figure 2: Methodology to compute the forward problem approximation employing autoadjustment regularization.