Study and comparison of different Machine Learning-based approaches to solve the inverse problem in Electrical Impedance Tomographies

Electrical Impedance Tomography (EIT) is a non-invasive technique used to obtain the electrical internal conductivity distribution from the interior of bodies. This is a promising method from the manufacturing viewpoint, since it could be used to estimate different physical inner body properties during the production of goods. Nevertheless, this technique requires dealing with an inverse problem that makes its usage in real-time processes challenging. Recently, Machine Learning techniques have been proposed to solve the inverse problem accurately. However, the majority of prior research is focused on qualitative results, and they typically lack a systematic methodology to determine the optimal hyperparameters appropriately. This work presents a systematic comparison of six popular Machine Learning algorithms: Artificial Neural Network, Random Forest, K-Nearest Neighbors, Elastic Net, Ada Boost, and Gradient Boosting. Particularly, the last two algorithms were based on decision tree learners. Furthermore, we studied the relationship between model performance and different EIT configurations. Specifically, we analyzed whether the measurement pattern and the number of used electrodes could increase the model performance. Experiments revealed that tree-based models present high performance, even better than Neural Networks, the most widely-used Machine Learning model to deal with EIT. Experiments also showed a model performance improvement when the EIT configuration was optimized. Most favorable metrics were attained using the tree-based Gradient Boosting model with a combination of both adjacent and mono measurement patterns as well as with 32 electrodes deployed during the tomographic process. With this particular setting, we achieved an accuracy of 99.14% detecting internal artifacts and a Root Mean Square Error of 4.75 predicting internal conductivity distributions.


Introduction
Different industrial scenarios, such as the wood manufacturing process, the inspection of internal structures, or the study of buried pipes, require a non-invasive analysis of target body internal compositions. Typically, bodies are evaluated using methods that qualitatively examine some inner physical property and generate a representative image with the corresponding inner feature distribution. After that, an expert would be able to generate a diagnosis according to the resulting image, which is usually enough for dealing with many of these situations. However, some industrial problems require a more extensive effort since it is necessary to know the specific property value for each target point.
Currently, there are several tomography techniques that allow the internal study of a body. The most well-known method is the Computed Tomography (CT) scan [1], which is widely used in medicine. However, the CT scan implies the use of ionizing radiation, which has many restrictions, such as the need for qualified staff and high-security measures. These constraints make it challenging to deploy this technology in industrial processes. Nevertheless, there are other less limited approaches that could be used not only in medicine but also in industrial environments. One of the emerging alternatives that have experienced a high development in the last years is the Electrical Impedance Tomography (EIT) [2].
EIT is a non-invasive tomography technique that does not involve the utilization of ionizing radiation. This technique is used to obtain information about the electrical conductivity distribution of a body [3]. For that purpose, a set of external electrodes are placed around a target body. A group of these electrodes is fed with small electrical currents in order to send a signal across the body. Then, the signal is altered according to the body's internal composition, and the remaining electrodes record the results. This process will then be repeated with several electrodes configurations following some stimulation pattern. With the recorded voltage values, it is possible to determine the internal electrical conductivity distribution, which can be used to estimate other physical properties such as the internal moisture distribution. Literature shows several EIT applications focused on various fields such as medicine (brain imaging [4][5][6], lung imaging [7,8], etc.), industrial processes [9], geophysical subsurface imaging [10] and a wide variety of engineering applications such as damage detection and strain localization [11][12][13][14].
The procedure of recording the electrode measurements is called the forward problem, and the method to determine the conductivity distribution from the voltages captured by the electrodes is usually known as the inverse problem, which is an ill-posed and highly nonlinear challenge [15].
Two main approaches have been traditionally used to solve the inverse problem. The first one includes algorithms such as back-projection and the one-step Gauss-Newton, that assume a certain linearity in the response of the body [16]. They are able to solve the problem faster at the expense of a worse accuracy. The second option is composed of iterative algorithms such as the Primal Dual Interior Point Method (PDIPM) [17] and the Iterative Gauss-Newton (IGN) [18], which are based on physical differential equations. They start from an initial hypothesis about the internal electrical conductivity distribution and iteratively change the simulation values in order to obtain a convergence between the simulations results and the actual electrode measurements. Although the latter option is more accurate than the first one, it requires both time and large computational resources. These prerequisites make it difficult to apply iterative algorithms in real-time scenarios.
The use of Machine Learning (ML) seems to be an appropriate approach to develop an accurate inverse problem solver. In fact, Literature shows notorious progress in the application of these types of techniques to deal with the inverse problem in the last decade. However, these approaches are usually focused on qualitative results and they typically lack of a systematic methodology to configure and to select the appropriate ML algorithm. Approaches based on Artificial Neuronal Networks (ANN) are largely used without a previous comparative process.
The aim of this article is to provide a systematic study of different ML models to determine what type of model attains better results in solving the inverse problem. In this work, the model hyperparameter optimization step, which is usually poorly addressed in Literature, is also considered. Despite, previously approaches are mainly focused on models based on ANNs, we conclude that tree-based models are able to surpass ANNs from both qualitative and quantitative viewpoints. We also analyzed the impact of the EIT forward problem process configuration from the ML algorithm accuracy viewpoint, which is a subject not enough explored, yet. This configuration is key to develop a ML-based inverse problem solver since forward problem outputs are the solver inputs. Moreover, this configuration has a direct impact on the deployment of a real system. The obtained conclusions give valuable configuration options to be considered in future deployments of ML-based EIT systems.
This paper is structured as follows: the mathematical basis behind EIT is summarized in Sect. 2; Sect. 3 introduces prior ML approaches for solving the EIT inverse problem; All the materials and methods used to develop this study, including software, ML algorithms, and datasets are described in Sect. 4; Sect. 5 describes the developed experiments and the appropriate discussion; and finally, Sect. 6 draws the conclusions and the future work.

Electrical impedance tomography
The EIT process consists of two main steps. On the one hand, the forward problem is the process of recording the electrode measurements, which is characterized through Eq. 1. This shows the governing equation for the voltage field generated by placing a current across a body [19]: where r is the electric impedance of the medium, I is the injected current, x is the frequency, e is the electric permittivity, and / stands for the electric potential. Assuming that a low frequency or direct current is used (x % 0), Eq. 1 can be reduced to Eq. 2, which is known as the standard governing equation for EIT: On the other hand, the inverse problem (also known as the Calderon's problem [20]) is the method used to determine the body conductivity distribution from the voltages captured by the electrodes. Thus, solving the inverse problem consists of determining the unknown impedance distribution r from the known injected current I and measured voltage V. However, this is a highly nonlinear problem, because the potential distribution is a function of the impedance (/ ¼ /ðrÞ). Furthermore, the inverse problem is also an ill-posed problem because of the diffusive nature of electricity.

Related work
Traditionally, both linear and iterative algorithms have been used to deal with the EIT inverse problem. However, in the last few years, Literature shows promising approaches based on ML. Typically, these approaches are composed of the following simplified main phases: 1) the generation of a dataset of target bodies with inner artifacts, which could be obtained using virtual or real data; 2) the recording of electrode measurements (forward problem), which could also be addressed via simulations or real measurements; 3) the development of an ML model to deal with the inverse problem; and 4) the testing phase, where the trained model is finally examined from the accuracy viewpoint. A simplified version of this methodology is shown in [21], where a dataset of 200 virtual body meshes was used to develop a Radial Basis Function (RBF)-ANN model for solving the inverse problem. The electrode measurements (forward problem) were obtained through simulation using the EIDORS software [22]. Half of the meshes were used to train the model and the remainder to test it. However, there is not a validation phase to appropriately configure the model. The model accuracy was not obtained comparing the ML predictions with the original inputs, but with inverse problem solutions from EIDORS algorithms.
Other approaches based on RBF-ANNs are described in [23,24]. Nevertheless, they were focused on detecting the correct position of the artifacts instead of the specific inner feature distribution values. Moreover, their metrics were obtained comparing their results with linear inverse problem solvers, and it is well-known that these types of solvers have poor accuracy.
The development of noise-robust ML algorithms was deeply explored by Sébastien Martin and Charles T.
M. Choi [25,26]. They presented methodologies able to deal with noise in voltage measurements and with target body deformations. Nonetheless, they only evaluated their approaches from a qualitative perspective.
Convolutional Neural Networks (CNN) [27] have also been used by several authors to solve the EIT inverse problem. A hybrid algorithm to reconstruct EIT images is shown in [28], where a D-bar method combined with a CNN model was applied in order to develop a robust approach to minimize the blurred effect. Wei et al. developed a CNN model to reconstruct ''heart and lung'' phantoms [29]. Nonetheless, the artificially-generated dataset used contained a limited number of elements (800 bodies for training and 40 for testing). A more systematic and comprehensive methodology can be found in [30], where a CNN method is proposed for image reconstruction of Electrical Resistance Tomography (ERT), whose inverse problem is mathematically the same as the EIT's one.
Performance comparison of various ML models is shown in [31], where Least Angle Regression (LARS), Elastic Net (EN), and ANN models were evaluated. They used a complex algorithm composed of a set of models, where each specific model instance was trained to predict the value linked to one pixel in the tomography imagery. They trained the ML models using an early stop methodology through a validation set. However, the network hyperparameters were not optimized. Results showed a good algorithm performance, since all the inner artifacts were detected. Nevertheless, a significant amount of noise was also visualized in the resulting images. A real application of the algorithm is described in [32], where the inner moisture of buildings was analyzed. Again, this research was only focused on qualitative results since the main goal was to visually identify the inner moisture areas using a color scale.
Recently, other approaches based on ANN were published [33,34]. The former approach carried out a minimal ANN hyperparameter optimization. However, a validation set was not used to do it. The latter work only used four virtual body meshes to test the models. In both cases, the model performance could have been overestimated.
Some authors have also explored the use of ANN to optimize the position of the electrodes used for measuring the voltages [27]. With this approach, they were able to improve the quality and reliability of the voltage measurements, reducing errors in EIT reconstructions.
To summarize, there are some common features linked to the majority of the available approaches that must be highlighted. On the one hand, the use of ML models based on some kinds of ANN is the most popular approach. Nonetheless, these alternatives usually lack a systematic methodology to both select and optimize the ML model and its hyperparameters. Moreover, there is also an absence of studies about how different EIT configurations affect the model accuracy. On the other hand, most of the approaches are usually focused on solving the inverse problem from a qualitative point of view. Nevertheless, the specific inner feature distribution values are necessary to deal with some industrial problems. In these scenarios, it is required to solve the inverse problem through a quantitative approach.

Materials and methods
In this paper, we developed and compared six popular Machine Learning algorithms in order to determine which one achieved better results for solving the inverse problem. In particular, we contrasted: ANN [35], EN [36], K-Nearest Neighbors (KNN) [37], Random Forest (RF) [38], Ada Boost (AB) [39] and Gradient Boosting (GB) [40]. All these algorithms were tested from both qualitative and quantitative points of view. It must be noted that RF, AB, and GB are ensemble methods, where decision trees were employed as ensemble components (estimators) to carry out all the experiments. Apart from the ML comparative, different forward problem configurations were also tested to analyze their impact on the model performance.
Generally, the methodology followed in this paper to develop each ML model is aligned with the traditional main phases cited in the previous section. In particular, a simulation approach was included to deal with both the dataset generation and the forward problem solution. This is a common approach to fairly test several experiment configurations, as it is shown in the linked Literature. A detailed explanation is presented in [35], which was used as the primary study of reference to develop our work. We highly encourage readers to read it for a deeper context.

Target bodies
A dataset of 10,220 virtual 2D body meshes was used in this work at the starting point. Each body was simulated through a circular mesh with a one-meter radius. Every mesh was composed by a set of triangular elements, and each triangle was set up with an electrical conductivity default value of 1 S/m. The specific number of triangles per mesh is linked to the experiment configuration.
In order to insert inner artifacts in the bodies, we used an algorithm previously developed [35] that allowed us to simulate a random number (from 1 to 3) of artifacts with arbitrary shapes. The artifact electrical conductivity values were set up according to a conductivity gradient value from a maximum of 100 S/m (artifact core) to a minimum of 40 S/m (artifact edge). Figure 1 shows meshes with one, two, and three inner artifacts inside. The color legend symbolizes the conductivity gradient values used to fill the mesh triangles.

Forward problem simulation
As mentioned before, the forward problem consists of recording the electrode measurements during the tomography process. There are different ways of injecting and recording currents via the electrodes placed around the target body. Typically, they are known as stimulation and measurement patterns, respectively [41].
We used EIDORS to simulate the forward problem. In particular, our experiments were performed employing different numbers of electrodes (8,16,32) placed equidistantly around the bodies. All our experiments used the adjacent pattern as the stimulation pattern. With this pattern, a fixed current of 0.010 A was injected iteratively in every pair of adjacent electrodes. It must be remarked that the number of iterations is equal to the number of used electrodes. At iteration n i , the current is injected between the electrodes i and i þ 1. In the final iteration, the pair of electrodes is composed of the electrode n and the electrode 1. At each iteration, a set of voltages is measured using the rest of the electrodes. We used different measurement patterns depending on the experiment. In order to clarify each one, we describe them in a specific scenario where 16 electrodes (16 iterations) were placed around the target body: • Adjacent pattern: at each iteration, the potential was measured in every pair of the remaining electrodes, obtaining 13 voltage values per iteration and 208 values altogether (see Fig. 2-left). • Mono pattern: at each iteration, the potential was individually measured in every one of the remaining electrodes with respect to a common ground, obtaining 14 voltage values per iteration and 224 values altogether (see Fig. 2-center). • Opposite pattern: at each iteration, the potential was recorded between opposite electrodes. The combinations linked to the electrodes used to inject the current are rejected, so 12 voltage measurements were obtained after each iteration. It should be highlighted that using this pattern the last six measurements are the opposite of the first six for each repetition (v 7 ¼ Àv 1 ; v 8 ¼ Àv 2 ; v 9 ¼ Àv 3 ; :::). At the end of the process, a bunch of 192 voltage values was captured (see Fig. 2-right).

Datasets
The datasets used to develop the ML models were composed of both the voltage measurements obtained after applying the forward problem over the virtual bodies (inputs) and the conductivity values of such virtual bodies (outputs).
We named observation to the combination of inputsoutputs linked to each virtual body. Specifically, we used six datasets to carry out the different experiments: the original one, which was previously developed in [35], and five newly generated datasets with different configurations. All of them were composed of 10,220 observations in accordance with the dataset's original size.
The number of voltage values per observation (inputs) is conditioned by the stimulation and measurement patterns as well as the number of electrodes employed. Taking into account that in all our datasets the adjacent pattern has been used as the stimulation pattern, the number of voltage values can be determined through Eq. 3: where n v is the number of voltages values, n e is the number of electrodes employed and m stands for the number of voltage values recorded at every iteration of the stimulation process. This last variable value depends on both the measurement pattern utilized and the number of electrodes used. In the specific case of 16 electrodes, it adopts a value of 13 for the adjacent pattern, a value of 14 for the mono pattern, and 12 for the opposite pattern. However, it must be remarked that we have also developed datasets joining the voltage measurements obtained from two different patterns for some experiments (see Sect. 5.1). In these cases, the individual results obtained via Eq. 3 should be added up. It must be also noted that when using the combination of the adjacent and mono patterns, there exists redundancy in the data since every voltage measure for the adjacent pattern complies with the following expression: In Eq. 4, i simbolizes the electrode number where the voltage is being meassured. This implies that every voltage measure obtained with the adjacent pattern can also be achieved by applying simple transformations to the voltages obtained through the mono pattern. The number of conductivity values per mesh (outputs) is always equal to the mesh density, considering mesh density as the number of triangles that compose the mesh. When using 8 and 16 electrodes, we have employed 844-triangles. However, in the case of using 32 electrodes, we have utilized meshes with 1,128 triangles. Table 1 shows the specific features associated with every employed dataset. These datasets are identified to make the references in the following sections easier. Specifically, D1 identifies the original dataset from [35]. All datasets are available online, so the experiments can be repeated if necessary (https://gitlab.citius.usc.es/cograde/ datasets_eit_ml_comparative).

Experiment results and discussion
In the experimental work to compare the proposed ML algorithms, we used the same methodology to develop all the different models. The dataset D1 was split into three sets: training (70% of the data), validation (15% of the data), and test (15% of the data). Every model was trained and optimized from the hyperparameter viewpoint using both training and validation sets. Particularly, we used an intensive Grid search optimization process [42] based on the model performance, which was determined via the MSE linked to the validation set, in order to select the optimal hyperparameters. The specific evaluated ranges are shown in Table 2. Once the best configuration was selected for every model, they were trained from scratch with the optimal hyperparameters. This was achieved using a fused training set composed of both the original training data and the validation set (85% of the available data). Finally, model metrics were obtained using the test set. The optimal hyperparameters selected for every model are also listed in Table 2.
Our starting point was the aforementioned study [35], and the ANN model was directly trained using the optimal hyperparameters proposed there, since both studies shared the same model development methodology. It must also be noted that the AB and the GB algorithms cannot handle a multioutput problem natively. To solve this inconvenience, we trained a specific model to predict the value of every mesh element. Then, metrics were calculated evaluating the resulting models (844) for each one of these algorithms. This model-point approach was also followed in other works [31], as described in the related work section.
In order to remove noise from predictions, the same post-processing technique proposed in [35] was applied over all the evaluated approaches. Basically, this process transforms the model outputs with values lower than a certain threshold into background. The procedure to select this optimal threshold was based on ROC curves.
Then, results were evaluated from both qualitative and quantitative perspectives according to the following metrics: RMSE, MAE, accuracy, and Cohen's Kappa coefficient (j) [43]. On the one hand, the RMSE and the MAE are quantitative metrics that allow us to analyze the model performance by identifying specific distribution values. On the other hand, the accuracy and the Cohen's Kappa coefficient measure qualitative results in order to evaluate the model performance to correctly place the inner artifacts inside bodies. Since qualitative metrics require two classes, mesh triangles with a background conductivity value of 1 S/m were associated with the negative class (they are not part of an inner artifact), while triangles with higher conductivity values were associated with the positive class (they are part of an inner artifact). Table 3 enumerates the obtained metrics for every model. From a quantitative viewpoint, the GB model achieved the lowest metrics (which implies a better model performance), cutting down the ANN metrics by more than a 50%. However, it is interesting to observe the detailed information about the specific model behavior for every test item, which is shown in Fig. 3. This figure depicts the RMSE boxplot per virtual body (mesh) from the test set. This metric allowed us to evaluate the model variability.
The GB showed the best RMSE with a restrained variability. However, it must also be observed that the KNN had significant variability, despite having a relatively low RMSE. This proves a model's weakness in some specific scenarios, which is undesired in industrial environments.
These outliers are related to meshes with several close artifacts, where the KNN had difficulty separating them in its predictions. Figure 4 shows a mesh with two artifacts and the tomography reconstruction by all the evaluated ML models. It can be noted that, in this case, all the models were able to separate the two inner artifacts. Nonetheless, in those situations where the artifacts were placed too close to each other, the EN and the KNN generated inaccurate reconstructions, joining the artifacts. This scenario can be observed in Fig. 5. From the qualitative perspective, Table 3 also shows that the GB achieved the highest accuracy and j (which implies a better model performance). However, considering only these metrics, the second-best model is the ANN, with very close results to the GB. This means that the ANN is capable of separating inner artifacts with high accuracy, but it has more difficulties than the GB to predict the exact conductivity values. Specific quality results for these two models are shown in Tables 4 and 5 through the respective confusion matrix.
In the first view, it is surprising that a decision-treebased model (GB) obtained the best performance, taking into account that ANN approaches are the most popular solvers in Literature. Nevertheless, it is worth remembering that most of the previously published approaches are only based on qualitative results, where model metrics are very similar. From all the evaluated metrics, we selected the GB as the most appropriate model to deal with the inverse problem, and it was consequently used for the rest of the experiments described in this paper, in particular, for the evaluation of the model performance with several EIT forward problem configurations.

Forward problem configuration: measurement pattern
The GB model performance was evaluated according to different measurement pattern configurations. Specifically, the following ones were studied: adjacent, mono, opposite, adjacent?mono. In order to evaluate them, the GB model was trained with the appropriate dataset (D1 to D4, respectively, see Table 1 for more details). It must be highlighted that the number of model inputs changes according to the pattern used. Thus, the selected model was trained from scratch for each combination. The same methodology described before, including the hyperparameter optimization and the post-processing phase, was used to develop the appropriate GB model for these experiments. Table 6 shows the model metrics obtained in the test phase for every one of the four studied pattern approaches. The lowest RMSE is attained by combining both the adjacent and the mono patterns. Although this is the configuration with the highest number of inputs, the improvement in the prediction is not merely the consequence of a higher input dimension because, as stated in Sect. 4.2, the adjacent data can be deduced from the mono data by applying linear transformations. This means that the adjacent-mono configuration does not contain additional information that did not exist previously in the adjacent or mono configurations. Therefore, the better performance is related to the fact that Gradient Boosting works better with certain sets of features over others representing the same information. Figure 6 shows an image with three artifacts and the reconstruction with the studied measurement patterns. The reconstruction made with the adjacent-mono measurement pattern was clearly the most similar one to the original image. Moreover, the reconstruction of the artifact boundaries was very precise.

Forward problem configuration: number of electrodes
Using the GB model as well as the measurement pattern (adjacent?mono) proposed in subsection 5.1, we carried out a new comparative focused on the number of electrodes involved in the forward problem. In particular, 8,16, and 32 electrodes were tested. Readers should note that the results from 16 electrodes were already obtained in  It should be pointed out that when the number of electrodes is modified, the number of the generated voltage values (inputs) also changes, as in the case of modifying the measurement pattern. Furthermore, when using 32 electrodes, the number of mesh elements (outputs) is also increased, from 844 to 1128. This dataset modifications require a new training phase from scratch. Again, the same methodology and hyperparameter optimization described in previous sections was used to retrain the GB model. Table 7 shows the model metrics according to the number of used electrodes. Globally, Table 7 and Fig. 7 show that best results were achieved with 32 electrodes. Furthermore, it can be noted that the higher the number of electrodes, the lower the RMSE. This is due to the fact that with a higher number of electrodes, the number of voltage values increases. Consequently, the input dimension also growths, so the ML models have more information in every observation to learn.
Electrode experiments must also be carefully analyzed in order to deploy the approach in real environments. Fewer electrodes would be able to simplify and speed up the forward problem. From a qualitative point of view, there is a significant improvement between 8 and 16 electrodes. Between these two combinations, the Kappa coefficient increases in four percentage units, so the probability of agreement by random chance is minimized. Moreover, the statistical dispersion is also cutting down, as is depicted in Fig. 7. Fewer differences were observed between 32 and 16 electrodes, but it must be noted that the 32-electrodes configuration is a more complex configuration that also forces to increase the mesh density. Mesh triangles with a conductivity value of 1 S/m were associated with class-0 (no artifact), while triangles with higher conductivity values were associated with class-1 (artifact). More than the 98% of the mesh elements were correctly classified with the ANN model. However, there are a 1.28% of the false negatives Mesh triangles with a conductivity value of 1 S/m were associated with class-0 (no artifact), while triangles with higher conductivity values were associated with class-1 (artifact). Almost the 99% of the mesh elements were correctly classified with the GB model. In addition, false positives and false negatives are reduced by almost 50% compared to the metrics obtained by the ANN model The column Optimal tuned hyperparameters shows the optimal hyperparameter values selected to configure the model in each scenario. RMSE and MAE symbolize quantitative metrics. Accuracy and j are values related to qualitative metrics

Conclusions and further work
The use of EIT is interesting in industrial environments because it allows inspecting target bodies internally in a non-invasive way. However, the solution of the inverse problem linked to the EIT method is a major limitation due to the fact that traditional techniques are not able to solve it in real-time accurately. Nevertheless, ML models have proved to be significantly superior to traditional solver algorithms, opening up opportunities to use EIT in industrial processes.
In this work, we carried out a systematic comparison of several ML approaches: ANN, EN, RF, KNN, AB, and GB. These models were compared as inverse solvers from both viewpoints: quantitative and qualitative. Although all the evaluated models achieved acceptable results, the GB model based on decision tree learners was the one that obtained the best metrics from both perspectives and with less statistical dispersion. Furthermore, this experiment  All of the evaluated approaches were trained with data generated with the measurement pattern selected in the previous step (adjacent ?mono) showed that all the tree-based models evaluated (RF, GB and AB) are able to achieve comparable and even better results than the ANN. The use of these models to solve the EIT inverse problem remains underinvestigated, so it could be a potential line of research.
In the last part of this work, we carried out different experiments linked to the impact of the forward problem configuration in the ML model accuracy, which is not frequently analyzed in Literature. In the first one, we compared the performance of the GB model using different kinds of measurement patterns: adjacent, opposite, mono, and the combination of the adjacent and mono patterns. Best results were attained with the last configuration, but the performance differences between distinct combinations were minimal.
In the second experiment, we compared different combinations of the number of electrodes during the forward problem. Clearly, a higher number of electrodes offers better results. Nevertheless, the combination with 16 electrodes is the most balanced approach between model performance and complexity. Again, the use of a more complex solution (32 electrodes) should be considered whether the highest accuracy is required.
Generally, configurations that allow increasing the number of inputs improve the model performance. This is mainly due to the EIT inverse problem being a very unbalanced challenge (more outputs than inputs), and this fact makes it difficult to train the models. Thus, more training information is helpful in the model development.
Oppositely, more training information usually means more complex systems to work in real time. Hence, there exists a tradeoff that should be analyzed according to the problem and the required accuracy.
All these experiments were carried out using virtual datasets generated with EIDORS. However, in future works, we plan to study a real use-case related to the wood industry. Particularly, we will use a ML approach to analyze the moisture distribution in wood-based products for detecting anomalies. In order to get an optimum ML model for this real use-case, we also expect to improve the hyperparameter selection phase evaluating larger parameter ranges as well as using a more complex search strategy. CiTIUS -Centro Singular de Investigación en Tecnoloxías Intelixentes da Universidade de Santiago de Compostela as a Research Center of the Galician University System. Computational resources were provided by the Galicia Supercomputing Center (CESGA).
We greatly appreciate the valuable comments from the reviewers that have helped us to improve our work.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.
Data availability The datasets generated during and/or analyzed during the current study are available in the following repository: https://gitlab.citius.usc.es/cograde/datasets_eit_ml_comparative Declarations Conflict of interest There is no potential conflict of interest in our paper. And all authors have seen the manuscript and approved to submit to your journal.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.