Logistic Regression with Wave Preprocessing to Solve Inverse Problem in Industrial Tomography for Technological Process Control

: The research presented here concerns the analysis and selection of logistic regression with wave preprocessing to solve the inverse problem in industrial tomography. The presented application includes a specialized device for tomographic measurements and dedicated algorithms for image reconstruction. The subject of the research was a model of a tank ﬁlled with tap water and speciﬁc inclusions. The research mainly targeted the study of developing and comparing models and methods for data reconstruction and analysis. The application allows choosing the appropriate method of image reconstruction, knowing the speciﬁcs of the solution. The novelty of the presented solution is the use of original machine learning algorithms to implement electrical impedance tomography. One of the features of the presented solution was the use of many individually trained subsystems, each of which produces a unique pixel of the ﬁnal image. The methods were trained on data sets generated by computer simulation and based on actual laboratory measurements. Conductivity values for individual pixels are the result of the reconstruction of vector images within the tested object. By comparing the results of image reconstruction, the most efﬁcient methods were identiﬁed.


Introduction
Tomography is a method of examining the interior of objects based on measurements taken at the edge of that object. It is a non-invasive method, which can be used to obtain a cross-section of the element-providing 2D images-or spatially-providing 3D images. The device for taking measurements is called a tomograph, and the image obtained from the measurements is called a tomogram. Mathematical operations and techniques called tomographic reconstruction are used to obtain the image. Different physical phenomena are used depending on the type of object under study. The carriers may be ultrasound, electron beams, electric currents or magnetic fields. Depending on the specificity of a given tomographic technique, we can observe both advantages and disadvantages in the areas of accuracy, frequency or resolution of reconstructed images. Familiarization with the characteristics of individual tomographic techniques allows for the proper selection of an image reconstruction method. Tomographic imaging offers a unique opportunity to study the complexity of a structure without interfering with the object [1,2]. There is an increasing demand for information about the behavior of the internal flows in process equipment. This information should be obtained in a non-invasive way with tomographic instrumentation. Typical measuring equipment is usually inadequate for the complicated process internal conditions, and sometimes its existence can interfere with the process. Industrial tomography is used to control measurement data. The facility to be monitored is a tomographically instrumented part of a manufacturing plant, where a device non-invasively collects data using excitatory electrical potentials on its surface. The instrument sends the raw data to a cloud computing system where an inverse problem (image reconstruction) is resolved. As a result, the state of the monitored object is classified using a machine learning algorithm. Figure 1 presents a general scheme of how the system works.
Energies 2021, 14, x FOR PEER REVIEW 2 of 20 opportunity to study the complexity of a structure without interfering with the object [1,2]. There is an increasing demand for information about the behavior of the internal flows in process equipment. This information should be obtained in a non-invasive way with tomographic instrumentation. Typical measuring equipment is usually inadequate for the complicated process internal conditions, and sometimes its existence can interfere with the process. Industrial tomography is used to control measurement data. The facility to be monitored is a tomographically instrumented part of a manufacturing plant, where a device non-invasively collects data using excitatory electrical potentials on its surface. The instrument sends the raw data to a cloud computing system where an inverse problem (image reconstruction) is resolved. As a result, the state of the monitored object is classified using a machine learning algorithm. Figure 1 presents a general scheme of how the system works. The presented solution enables tomographic sensor-based process control, Big Data analysis, multidimensional industrial process control, advanced human-machine interfaces and knowledge-based process monitoring. Sensor technology is mostly related to electrical tomography (ET) [3,4], which includes capacitive tomography (ECT) [5][6][7][8][9] and electrical resistive tomography (ERT) [10,11]. This solution allows the image to be reconstructed as the conductivity or permittivity distribution of the reservoir under study by measuring its boundary. Another method is optical coherence tomography (OCT), which is an optical technique that allows non-invasive imaging of specimens [12], and ultrasound tomography (UST) [13,14], which is a technique that uses the information included in the ultrasound signal after it has passed through the test object.
Electrical impedance tomography (EIT) is a non-destructive method for creating image reconstructions in various application areas. It is suitable for real-time visualization and the analysis of the electrical conductivity distribution inside the analyzed object [15]. EIT is now a widely used tomographic imaging technique and applies to many areas of everyday life. This technique has found applications in fields such as medical diagnostics [4], industrial process monitoring [13] and geophysical surveys [16]. Mathematical reconstruction of conductivity in EIT involves solving a nonlinear and ill-posed inverse problem from noisy data [17,18]. The current development of EIT algorithms is heavily biased towards machine learning methods [3], and for such an operation, the data must be prepared appropriately. It raises the need to check whether particular algorithms are better than classical ones and can be used. Compared to other known imaging modalities, impedance tomography has many advantages. However, EIT reconstruction can be unstable, and the disadvantage of solving a backward problem should be mentioned. Furthermore, the sensitivity of EIT solutions to measurement, numerical and model errors implies the need to adapt model parameters to specific cases. At the same time, it should be remembered that most reconstruction methods in 2D space are equally effective in 3D. The presented solution enables tomographic sensor-based process control, Big Data analysis, multidimensional industrial process control, advanced human-machine interfaces and knowledge-based process monitoring. Sensor technology is mostly related to electrical tomography (ET) [3,4], which includes capacitive tomography (ECT) [5][6][7][8][9] and electrical resistive tomography (ERT) [10,11]. This solution allows the image to be reconstructed as the conductivity or permittivity distribution of the reservoir under study by measuring its boundary. Another method is optical coherence tomography (OCT), which is an optical technique that allows non-invasive imaging of specimens [12], and ultrasound tomography (UST) [13,14], which is a technique that uses the information included in the ultrasound signal after it has passed through the test object.
Electrical impedance tomography (EIT) is a non-destructive method for creating image reconstructions in various application areas. It is suitable for real-time visualization and the analysis of the electrical conductivity distribution inside the analyzed object [15]. EIT is now a widely used tomographic imaging technique and applies to many areas of everyday life. This technique has found applications in fields such as medical diagnostics [4], industrial process monitoring [13] and geophysical surveys [16]. Mathematical reconstruction of conductivity in EIT involves solving a nonlinear and ill-posed inverse problem from noisy data [17,18]. The current development of EIT algorithms is heavily biased towards machine learning methods [3], and for such an operation, the data must be prepared appropriately. It raises the need to check whether particular algorithms are better than classical ones and can be used. Compared to other known imaging modalities, impedance tomography has many advantages. However, EIT reconstruction can be unstable, and the disadvantage of solving a backward problem should be mentioned. Furthermore, the sensitivity of EIT solutions to measurement, numerical and model errors implies the need to adapt model parameters to specific cases. At the same time, it should be remembered that most reconstruction methods in 2D space are equally effective in 3D.
The article presents an improved method of monitoring and optimizing processes in heterogeneous tank reactors in which specific reactions occur. The method used relates to electrical tomography, and the innovation is the original way of using a hybrid module sys-Energies 2021, 14, 8116 3 of 21 tem in parallel, combining logistic regression with wavelet preprocessing. The significant difference of the presented method over other non-invasive solutions is the reliability of disturbances arising during measurements and the accuracy of imaging reconstruction. In addition, our proposed method enables the appropriate selection of wavelet preprocessing algorithms for image pixels. It results in better reconstruction quality and higher image resolution.
This paper consists of four sections. The architecture of the designed system, industrial processes and application platform is presented in Section 2. Machine learning methods used for image reconstruction, numerical models, tomographic devices and laboratory measurement systems are also discussed. The results of the research work are presented in Section 3 in the form of image reconstructions for synthetic and measured data. Section 4 contains a discussion of the results obtained. Finally, Section 5 summarizes the research work.

Materials and Methods
This section shows the tomographic algorithms, system architecture, industrial tomography, measurement hardware, mathematical algorithms and measurement models used to reconstruct images from real and synthetic data.
The research used the SmartEIT 1.0 (Netrix S.A., Lublin, Poland), our electrical impedance tomograph and a specially prepared tank with EIT measuring electrodes.
For preliminary analyses of numerical models, Python 3.6 (Python Software Foundation, Amsterdam, Netherlands) with NumPy, math, SciPy libraries were used. MATLAB 2020B (MathWorks, Natick, MA, USA) and EIDORS (Ver: 3.9) (Sourceforge, San Diego, CA, USA) [19] were implemented to construct the case generator and visualize the reconstructions. Using R software (The R Foundation, Vienna, Austria), structural parameters were estimated, ROC analysis was performed and the relationship between FPR and TPR was graphically represented. For this purpose, the following packages (The R Foundation, Vienna, Austria) were used: R.matlab (Ver: 3.6.

Novelty of the Proposed Solution
The novelty of the solution presented in this paper is the combination of logistic regression with wavelet preprocessing methods to reconstruct the output image in industrial tomography. The approach proposed by the authors consists in implementing an algorithm with multiple trained subsystems. After converting the predictors to a binary grid, image reconstruction is generated. Voltage drop measurements are used as input data.

Industrial Electrical Tomography
Industrial tomography is used to analyze the technological processes inside the studied facility [2,5]. This approach allows for better real-time process control. Data concentration profiles, phases and chemicals can be studied using fast data acquisition and image reconstruction. The resulting data can be used to monitor process response, improve quality, yield and flow rate. The imaging technique presented here exploits the respective electrical properties of various kinds of substances. Electrical tomography can be divided into impedance tomography and capacitive tomography for dielectric systems. In this method, a current or voltage source is connected to the test object, the voltage distribution at the edge of the test object is measured using a measurement system. An image reconstruction algorithm processes the data collected from the above measurements. The reconstructed image is called a tomogram. It is worth noting that tomography is characterized by low image resolution. This has to do with the number of measurements, which is limited by the nonlinear current flow through the tested element but is also caused by the low sensitivity of the measured voltage waveforms, which depends on the conductivity changes inside the area.

Measuring Device
The device's design is based on the idea of electrical tomography, which is based on the non-invasive measurement of voltage drops from electrodes on the tested object ( Figure 2). Based on measurements of voltages on electrodes directly adjacent to a given medium, it is possible to determine impedance's spatial distribution and thus visualize its internal structure. Due to the built-in microcomputer, it is possible to perform EIT measurements and view the reconstructions created on their basis. In addition, the device has a network interface that allows data transfer to an external server. The tomograph measures voltages by switching channels according to the polar method. First, EXC and GND outputs are connected to two opposite electrodes using multiplexers. The intensity of current flow is programmed to a set value. Then, the signal input is connected successively to the remaining electrodes, on which the voltage is measured to the GND electrode. After completing measurements, the measuring information is disconnected, the forcing electrodes are switched to the next pair and the cycle is repeated. the nonlinear current flow through the tested element but is also caused by the low sensitivity of the measured voltage waveforms, which depends on the conductivity changes inside the area.

Measuring Device
The device's design is based on the idea of electrical tomography, which is based on the non-invasive measurement of voltage drops from electrodes on the tested object (Figure 2). Based on measurements of voltages on electrodes directly adjacent to a given medium, it is possible to determine impedance's spatial distribution and thus visualize its internal structure. Due to the built-in microcomputer, it is possible to perform EIT measurements and view the reconstructions created on their basis. In addition, the device has a network interface that allows data transfer to an external server. The tomograph measures voltages by switching channels according to the polar method. First, EXC and GND outputs are connected to two opposite electrodes using multiplexers. The intensity of current flow is programmed to a set value. Then, the signal input is connected successively to the remaining electrodes, on which the voltage is measured to the GND electrode. After completing measurements, the measuring information is disconnected, the forcing electrodes are switched to the next pair and the cycle is repeated. Raspberry Pi sends and receives data from the measurement module using the SPI interface; Raspberry acts as Master by default, module as Slave. On receiving the command to start the sequence-the module starts measuring using the opposite method. Electrodes A and B force the flow of the alternating current of a sinusoidal waveform, where electrode A is connected via a multiplexer to the circuit of current intensity regulation, and electrode B is connected to GND of the circuit. The voltage is measured at successive electrodes to GND ( Figure 3). Raspberry Pi sends and receives data from the measurement module using the SPI interface; Raspberry acts as Master by default, module as Slave. On receiving the command to start the sequence-the module starts measuring using the opposite method. Electrodes A and B force the flow of the alternating current of a sinusoidal waveform, where electrode A is connected via a multiplexer to the circuit of current intensity regulation, and electrode B is connected to GND of the circuit. The voltage is measured at successive electrodes to GND ( Figure 3).

Measurement System
The measurement system consists of a sensor tank; devices with a model for data transmission, collection and analysis are connected to a communication interface whose task is to read the signal from the selected sensor, process it into a coherent form and then send the read and processed data to the acquisition module ( Figure 4). For this solution, an application was prepared with an algorithm for image reconstruction, which, using learning data obtained by computer simulation from real models, was trained to solve the inverse problem. Conductivity values of the individual pixels of the output vector allow the obtaining of images of the interior of the studied objects.

Algorithms and Methods
This section presents the descriptions of the used algorithms. Industrial tomography also belongs to electromagnetic field inverse problems. An inverse problem is an optimization, identification or synthesis process where parameters are determined by domainspecific information. Such problems do not have straightforward solutions and are illconditioned. Finite element methods solved the forward problem.
There are many optimization methods, using both deterministic and statistical algorithms [20][21][22][23][24][25][26][27]. Machine learning is a group of methods that is gaining increasing popularity in various types of tomography, including EIT. Applying the above methods

Measurement System
The measurement system consists of a sensor tank; devices with a model for data transmission, collection and analysis are connected to a communication interface whose task is to read the signal from the selected sensor, process it into a coherent form and then send the read and processed data to the acquisition module ( Figure 4). For this solution, an application was prepared with an algorithm for image reconstruction, which, using learning data obtained by computer simulation from real models, was trained to solve the inverse problem. Conductivity values of the individual pixels of the output vector allow the obtaining of images of the interior of the studied objects.

Measurement System
The measurement system consists of a sensor tank; devices with a model for data transmission, collection and analysis are connected to a communication interface whose task is to read the signal from the selected sensor, process it into a coherent form and then send the read and processed data to the acquisition module ( Figure 4). For this solution, an application was prepared with an algorithm for image reconstruction, which, using learning data obtained by computer simulation from real models, was trained to solve the inverse problem. Conductivity values of the individual pixels of the output vector allow the obtaining of images of the interior of the studied objects.

Algorithms and Methods
This section presents the descriptions of the used algorithms. Industrial tomography also belongs to electromagnetic field inverse problems. An inverse problem is an optimization, identification or synthesis process where parameters are determined by domainspecific information. Such problems do not have straightforward solutions and are illconditioned. Finite element methods solved the forward problem.
There are many optimization methods, using both deterministic and statistical algorithms [20][21][22][23][24][25][26][27]. Machine learning is a group of methods that is gaining increasing popularity in various types of tomography, including EIT. Applying the above methods

Algorithms and Methods
This section presents the descriptions of the used algorithms. Industrial tomography also belongs to electromagnetic field inverse problems. An inverse problem is an optimization, identification or synthesis process where parameters are determined by domain-specific information. Such problems do not have straightforward solutions and are ill-conditioned. Finite element methods solved the forward problem.
There are many optimization methods, using both deterministic and statistical algorithms [20][21][22][23][24][25][26][27]. Machine learning is a group of methods that is gaining increasing popularity in various types of tomography, including EIT. Applying the above methods requires the 2D cross-section of the test object to be divided into pixels by the finite element method. Typically, a single algorithm converts a set of measurements into the pixels that make up a tomographic image. It can also train separate models to convert measurements to a single image pixel. It then needs to train as many models as the resolution of the tomographic image. It was the approach taken in the research presented here.
For logistic regression, the implementation probability of the output variable is calculated relative to the corresponding category, where the inclusion probability is estimated [28][29][30][31][32][33]. Thus, the implementation of logistic regression enables the resolution of the imaging domain to be determined. When creating a reconstruction for EIT, we must estimate conductivity for each finite element. In the presented case, we calculate the probability that finite elements belong to inclusion. For this purpose, we define a logistic regression model [28] for each finite element. Thus, the reconstruction based on logistic regression must be defined as a set of logistic regressions corresponding to finite elements.
For each finite element, we define a learning set where the class membership is represented as y i ∈ {0, 1} for 1 ≤ i ≤ n. In the analyzed case of inclusion for a pixel, we take y i = 1, whereas if the finite element does not contain an inclusion, we take y i = 0. By observing the signal received from the electrodes x ∈ R m , we make a classification for each finite element. Logistic regression has been used to construct a classifier (mapping) f : R m → {0, 1} .
Let (Ω, F , P) be a probabilistic space and Y a random variable with a discrete distribution, where Y : Ω → {0, 1} . We determine P(Y = y|X ), where y ∈ {0,1}, x ∈ R m . Chances are the ratio of the probability of being successful to the probability of failure.
The logistic regression task is to estimate the probability of success P(Y = 1|X = x ) based on the realization x ∈ R m and we assume P(Y = 1|X = x ) = p(x).
We consider the relationship defined by the equation [28,29] ln θ(x) = ln p(β 0 , β, x) where ε is a random variable with normal distribution N 0, σ 2 and β = (β 1 , β 2 , . . . , β m ) ∈ R m . If there is an intercept in model (2), then we have From Formula (2) we obtain To estimate unknown parameter β we apply the Maximum Likelihood Estimation. Let Y ∈ {0, 1} n and X ∈ R n×(m+1) , where x (1) x (2) . . . The estimation of the parameters consists of solving the task max β L(β, Y, X), where the objective function L(β, Y, X) is given by the equation [28,29] Replace task (5) with max β ln L(β, Y, X), whereby applying the Formula (3), we determine that the natural logarithm of the objective function (6) is following A necessary condition for the existence of an extremum is (2) . . .
The matrix of derivatives of order two is equal to where The matrix of second derivatives defined by Equation (3) is negatively defined. For We estimate the values of the β by the formula The main problem in EIT is the problem of collinearity of measurements obtained from electrodes (see, e.g., [3,9,10,12,16,24]). Extracting a set of stochastically independent features is impossible, so techniques such as Tikhonov regularization, LASSO, elastic net are usually used (see, e.g., [28,29]). Another approach can be used to reduce the redundant variables, such as the discrete wavelet transform [34][35][36]. It consists in determining the projection of the signal {x t } 1≤tleqm obtained from the electrodes to an orthogonal basis. In the considered case, we use the elastic net for feature reduction and decomposition using wavelet preprocessing.
Let Ψ(t) be an orthogonal wavelet basis (mother wavelet) (see, e.g., [34]). For j ∈ Z, we define a sequence Ψ jk k∈Z as follows: and sequence φ jk k∈Z where φ(t) denotes scaling. Thus, the time series {x t } 1≤tleqm can be presented as follows [34,36]: where the value c jk denotes the scaling coefficient, but d ik is the complex coefficient. Due to decomposition level j ∈ Z based on Formula (11), the sequence {x t } 1≤tleqm can be expressed in different forms. The functions Ψ jk (t) and φ jk (t) take non-zero values on a bounded interval. From above, the sequence {x t } 1≤t≤m can be presented as follows: where s << m. Accordingly, we define a projection operator for the time series {x t } 1≤t≤m at level j on an orthogonal basis φ jk (t) 0≤k≤s as follows: In EIT, the sequence of scaling coefficients c jk of projection has been used as input variables for Logistic Regression.

Results
A measuring system consisting of a SmartEIT tomograph and a measuring object with appropriate electrodes was used ( Figure 5). A numerical model of the reservoir was built in order to generate a suitable learning data set. The finite element method was used to design the tank's cross-sectional mesh and the sensor system with MATLAB and EIDORS tools. Methods have been developed to generate learning data to solve the forward problem. For each instance, there is a vector of measurements and an image created on a 2D pixel.
Several tens of measurements were made on such a physical model, adjusting the parameters of the mathematical model based on the generalized Laplace equation, which correctly generated measurement values based on a dense finite element mesh. Based on this model, a set of 50,000 cases containing both the measurements and the corresponding conductivity distributions was generated. In addition, Gaussian noise with standard Energies 2021, 14, 8116 9 of 21 deviation was added to each measurement value. Reference measurements were made several times with the tomograph. By analyzing the signals obtained from the electrodes, the noise was estimated, which averaged about 4% of the standard deviation of the reference measurements. Therefore, for reference, Gaussian white noise was added to the sonde generated from the electrodes in EIDORS, and then reconstructions were determined for the noisy signal using logistic regression models. To assess the quality of the visual area reconstruction, below we present some basic Receiver Operating Characteristic (ROC) analyses [27]. For this purpose, the finite element that does not belong to the inclusions is described as a negative case (N) and interpreted as an element that belongs to the background. On the other hand, the finite element that belongs to the inclusions is taken as the positive case (P). Therefore, to determine the basic characteristics, the confusion matrix is first determined in the following way: TP (True Positive)-denotes the finite elements that correctly belong to the inclusion area; TN (True Negative)-the number of finite elements that are correctly recognized as belonging to the background; FP (False Positive)-the number of finite elements belonging to the background that are recognized as having belonged to the inclusion area (false alarm); FN (False Negative)-the number of finite elements belonging to the inclusion area but recognized as background.
We determine the basic coefficients as follows [27]: To determine the ability of a classifier based on the use of logistic regression [15], we determine a curve describing the Receiver Operating Characteristic (ROC) curve [27]. The ROC curve shows the relationship between True Positive Rate (sensitivity on the Y-axis) To assess the quality of the visual area reconstruction, below we present some basic Receiver Operating Characteristic (ROC) analyses [27]. For this purpose, the finite element that does not belong to the inclusions is described as a negative case (N) and interpreted as an element that belongs to the background. On the other hand, the finite element that belongs to the inclusions is taken as the positive case (P). Therefore, to determine the basic characteristics, the confusion matrix is first determined in the following way: TP (True Positive)-denotes the finite elements that correctly belong to the inclusion area; TN (True Negative)-the number of finite elements that are correctly recognized as belonging to the background; FP (False Positive)-the number of finite elements belonging to the background that are recognized as having belonged to the inclusion area (false alarm); FN (False Negative)-the number of finite elements belonging to the inclusion area but recognized as background.
We determine the basic coefficients as follows [27]: To determine the ability of a classifier based on the use of logistic regression [15], we determine a curve describing the Receiver Operating Characteristic (ROC) curve [27]. The ROC curve shows the relationship between True Positive Rate (sensitivity on the Y-axis) and False Positive Rate (1-specificity on the X-axis) for different settings of cut-off levels (levels of class membership discrimination). Using logistic regression in EIT, the ROC curve shows the detection ability of a binary classifier that determines the affiliation of finite elements to the region of the inclusion for different settings of probability levels. The most optimal models can be determined by analyzing the class of models (a set of classifiers) using the ROC analysis (characterized by a higher TPR value for each established FPR). In some cases, when analyzing the ROC curves, determining the optimal classifier is quite difficult, while as the optimal classifier, we choose the one for which the Area under the Curve (AUC) is the greatest. This quantity is also included in Tables 1-3 describing the reconstructions.
To assess the reliability of reconstruction (agreement between pattern and reconstruction), we estimate Cohen's kappa ratio as follows: .
Generally, κ ∈ [0, 1]. When kappa has the higher value, the agreement between pattern and reconstruction is higher, and when there is complete agreement κ = 1.
To verify whether there is a significant disagreement between pattern and reconstruction, we apply McNemar's test. This test compares the sensitivity and specificity of reconstruction. The McNemar's test statistic is defined as follows: and has χ 2 distribution with 1 degree of freedom. Three patterns are shown below with their reconstructions. The patterns were simulated in MATLAB, and signals were determined as voltage vectors from the electrodes using EIDORS. Figures 6-11 summary the image reconstruction for the example 1.  and has distribution with 1 degree of freedom. Three patterns are shown below with their reconstructions. The patterns were simulated in MATLAB, and signals were determined as voltage vectors from the electrodes using EIDORS. Figures 6-11 summary the image reconstruction for the example 1. Figures  6, 8 and 10 show the results of these reconstructions, but Figures 7, 9 and 11 show the ROC curve for reconstructions presented in Figures 6, 8 and 10. Tables 1-3 present the basic coefficient of ROC analysis, Cohen's kappa, McNemar's test statistic and p-value. Figure 6. Image reconstruction for one object-example 1. Figure 6. Image reconstruction for one object-example 1.     Table 2. Summary of reconstruction of pattern presented in Figure 8.  Figure 8. Image reconstruction for one object-example 2.     Table 3. Summary of reconstruction of pattern presented in Figure 10.  Figure 10. Image reconstruction for two objects-example 3. Figure 10. Image reconstruction for two objects-example 3.

Elastic
Energies 2021, 14, x FOR PEER REVIEW 13 of 20 Figure 11. ROC analysis for example 3. Table 3. Summary of reconstruction of pattern presented in Figure 10. In order to verify the performed numerical tests, measurements were made on a laboratory tank. The image reconstruction results based on logistic regression with elastic net and different wavelets decompositions are shown in Figures 12-20.  In order to verify the performed numerical tests, measurements were made on a laboratory tank. The image reconstruction results based on logistic regression with elastic net and different wavelets decompositions are shown in Figures 12-20.

Discussion
The monitoring system presented is designed to automate, analyze and optimize processes using industrial tomography to analyze without disturbing its interior. This solution allows for better monitoring and control of processes. The system was designed using electrical tomography, which is used to study technological processes. An image reconstruction algorithm then processes the collected data. The main challenge in electrical tomography is the construction of precise measurement devices and image reconstruction methods. The idea of the measurement system was based on tomographic sensors for data collection and gathering and through appropriate communication protocols for processing in the computing cloud.
One of the problems in EIT image reconstruction using binary classifiers is the selection of variables that significantly influence the classifier's result. In this study, the elastic net technique and wavelet decomposition were used to reduce the dimensions of the signal obtained from the electrodes. To answer the question of which method gives the best reconstruction, the same patterns were simulated at the beginning. Next, the reconstructed images have been compared with the patterns. Since visual assessment is not objective enough, popular quantitative Receiver Operating Characteristics were estimated: Accuracy, Sensitivity, Specificity, Positive Predictive Value, Negative Predictive Value, Detection Rate and AUC. Additionally, the assessments described the reliability of reconstructions (Cohen's kappa) and the tests of disagreement between pattern and reconstruction (McNemar's test) were provided.
Analyzing such characteristics as accuracy (the proportion of the visual area that the model correctly recognized) and specificity (the proportion of correctly recognized finite elements belonging to the background) is quite accurate. Thus, the quality of reconstruction is quite good. Nevertheless, from sensitivity (the proportion of correctly recognized finite elements belonging to the inclusion) analysis, we see that the best reconstructions were obtained by applying "db1", "db2", "sym1", "sym2" and "bior1.1" wavelet decompositions. An identical result is obtained by analyzing the values of Cohen's kappa. Quite good reliability of reconstructions was obtained for "db1", "db2", "sym1", "sym2" and "bior1.1" wavelet decompositions too. By analyzing the AUC values, the best reconstructions are obtained for the dimension reduction techniques: elastic net and signal decomposition using "db1", "db2", "coif1", "sym1" and "bior1.1".
Based on McNemar's test, we can conclude that for examples 1-2 at a significance level of 0.05, the null hypothesis (pattern and reconstruction are consistent or disagreement is irrelevant) is rejected in favor of the alternative hypothesis (disagreement between pattern and reconstruction is significant) for "db7", "coif2" and "bior1.5" wavelet decompositions, and "db7" and "coif1" wavelet decompositions for example 3.

Conclusions
This study aimed to develop algorithms based on Logistic Regression and Wavelet Preprocessing methods to solve the inverse problem in EIT. The research focused on the development and comparison of algorithms and models for image reconstruction. All the algorithms presented are well suited to practical implementations in EIT. Depending on the number of measuring sensors used and the analyzed patterns, the obtained quality of the reconstructed images varies in different methods. To perform the reconstruction in EIT, we estimate the conductance for each finite element typically. In the presented case, we calculate the probability that a finite element belongs to an inclusion. For this purpose, we create a logistic regression model for each finite element. In this way, a logistic regression-based reconstruction is defined as a set of logistic regression models resting on the finite elements.
The use of logistic regression for image reconstruction in EIT should be preceded by selecting independent variables for each finite element separately. In order to reduce the dimensions of the independent variables, the elastic net technique was used in this work, as well as wavelet decomposition of signals, which is a good alternative tool to the elastic net. However, in the case of several objects in the field of view, objects located in the center are recognized poorly.
Further work could focus on ROC analysis for finite elements located in the center of the field of view, determining the thresholding for binary classifiers for finite elements depending on the position of the element in the field of view, and improving reconstruction in the center of the field of view (improving classifiers for finite elements further away from the electrodes). For further research and analysis, the decompositions "db1", "db2", "coif1", "sym1", "sym2" and "bior1.1" will be used to reduce the over dimensionality of independent variables in logistic regression models. Additionally, we investigated the dependence of the fit quality of logistic models (e.g., using the Hosmer-Lemeshow test) at different locations in the field of view. It could be the subject of a larger study to improve the quality of the reconstruction further.
Thus, the presented research results contain essential information that may accelerate the development of machine learning methods in industrial tomography. Furthermore, the research contributes to improving the efficiency of tomographic imaging with the use of electrical properties.