A Coupled Double-Layer Electrical Impedance Tomography-Based Sensing Skin for Pressure and Leak Detection

There is an extensive need for surface sensors for applications such as tactile sensing for robotics, damage and strain detection for structural health monitoring and leak detection for buried structures. One type of surface sensor is electrical impedance tomography (EIT)-based sensing skins, which use electrically conductive coatings applied on the object’s surface to monitor physical or chemical phenomena on the surface. In this article, we propose a sensing skin with two electrically coupled layers separated by an insulator. Based on electrical measurements, the spatial distribution of the electrical coupling between the layers is estimated. This coupling is sensitive to both the pressure distribution on the surface and water entering between the layers through a leak. We present simulations and experimental studies to evaluate the feasibility of the proposed method for pressure sensing and leak detection. The results support the feasibility of the proposed method for both of these applications.

Although sensing skins show promise in several of the fields that they are studied in, they have certain challenges.One of them is that there can be several different physical and chemical factors contributing to the conductivity of the sensing skin, and, therefore, it may be difficult to distinguish the effect of interest.Hence, there is still need for innovative smart surface sensors for, e.g., detecting surface pressure in robotics and detecting leaks in buried structures.
In publication [37], a functional double-layer sensing skin was proposed for SHM: the two layers of the sensing skin were used separately, yet simultaneously, to distinguish different physical/chemical effects on the surface-particularly, one layer was sensitive to cracking while the other one was sensitive to both cracking and strong electrolytes, such as chlorides dissolved in water in contact with the surface.Another type of multilayer sensing skin proposed for tactile sensing is based on having a second more conductive layer that is weakly separated from the original measured surface [7,38,39].Applied pressure introduces a contact between the two surfaces, thus effectively increasing the local conductivity of the measurement layer on the region of pressure.
In this paper, we propose a novel EIT-based sensing skin technique that is based on the use of two parallel electrically conductive layers.Similarly to the sensing skins presented in [38,39], these layers are coupled electrically, causing current to flow between the layers.In addition to the pressure on the surface, the electrical coupling is also sensitive to water entering in between the two layers.Unlike in the previous studies, we use two similar conductive layers and have electrodes attached on both of the conductive surfaces.We also measure currents flowing between electrodes of different surfaces.We formulate the mathematical model for the proposed technique and introduce its computational implementation.We study the feasibility of the proposed technique for pressure sensing and leak detection using numerical simulations and experiments.

Single-Layer EIT-Based Sensing Skin
EIT-based sensing skins consist of a conductive surface with electrodes attached to it, through which electrical currents can be injected into the surface.Traditionally, using a so-called current injection EIT device, the currents of the electrodes are set to different injection patterns and the resulting potentials are measured.In contrast, so-called potential excitation EIT devices set the potentials of the electrodes to different excitation patterns and measure the resulting currents.In many cases, both the theory and applications work similarly independent of which type of device is used.In this study, we use a potential excitation device for measurements, and hence we also formulate the theory accordingly.We use the conventional, single-layer sensing skin model in the reconstruction of the background conductivities of the two layers before they are installed in the coupled double-layer sensor system.

Forward Problem
The EIT measurements (whether potentials or currents are measured) are usually modeled by the complete electrode model (CEM), which consists of the following equations [40][41][42]: where σ and u are the conductivity and potential inside the imaging domain Ω, z l is the contact impedance of the l:th electrode, n is the outward unit normal vector on the boundary ∂Ω, U l and I l are the potential and current on the l:th electrode, ∂Ω e l is the boundary region corresponding to the l:th electrode, x is the spatial variable and N el is the number of electrodes in the system.Either the potentials U l or currents I l can be solved from the CEM as long as the other ones are known.In this article, as we use a potential injection device, correspondingly, we solve the currents from the CEM.In general, the potential u inside the domain has to be solved simultaneously with the electrode potentials or currents.Finding an analytical solution to the CEM is generally infeasible, and hence a numerical approximation is used instead.In this article, we use a finite element method (FEM) with piecewise linear basis functions to approximate the forward model (for details, see [43]).In FEM, the imaging domain is approximated by a mesh consisting of nodes and elements.Spatially distributed quantities inside the domain are then approximated as a linear combination of basis functions ϕ i .For example, the conductivity on the domain is then approximated as where σ i ∈ R are coefficients of the basis functions and N is the number of basis functions.
Using this approximation, we identify the conductivity by the vector of coefficients σ i in ( 5): From the FEM, we obtain the computational model Ĩσ ( σ) mapping the conductivities to a vector containing electrode currents corresponding to all the used excitations.

Inverse Problem
To solve the coefficients of conductivity σ, a variety of methods exist.In this article, we use the Bayesian framework, in which unknown or uncertain parameter values are treated as random variables.The Bayesian framework also allows the use of prior knowledge of the measured system and the use of reference measurements.Further information can be found, e.g., in [44].
In this article, we use the maximum a posteriori (MAP) point as the reconstruction.The MAP point can be found by solving the following optimization problem [44]: where I m contains all the measured currents from different excitations stacked into a single vector, R( σ) is the prior potential function (more details in Section 3.2) and L is the Cholesky factor of the data precision matrix.
The optimization problem (7) can be solved by many optimization algorithms.In this study, we use the Gauss-Newton (GN) method [45,46] with a line search to find the step length.For each iteration σi , the next GN iteration σi+1 with step length τ is computed by where J is the Jacobian of Ĩ( σ) with regard to σ and H R and ∇R are the Hessian and the gradient of R. The step length τ is chosen optimally with respect to the objective function

A Coupled Double-Layer Sensing Skin
The double-layer sensing skin presented in this paper consists of two conductive surfaces separated by an insulating layer, as depicted in Figure 1.Although the middle layer is insulating, some current will pass through it as displacement current since alternating potential difference is applied to the electrodes.Alternatively, in case of leak detection, water can form a conductive pathway between the layers.We use single-electrode excitations, i.e., one electrode at a time is excited to a higher potential, while the others remain grounded.Currents passing through the grounded electrodes are then measured.In this paper, we model the double-layer sensing skin as two two-dimensional conductive surfaces modeled by the CEM and a coupling between the surfaces.More specifically, the coupling is the conductance per area, and thus, together with the local potential difference between the layers, it determines the current density flowing between the layers.This coupling strengthens if water leaks into the space between the layers (a conductive pathway) or if applied pressure pushes the layers closer to each other (increased local capacitance).For the purely capacitive case, the coupling can be computed by assuming that the layers form a parallel plate capacitor with possibly spatially varying distance d(x) between the layers.For a small area A having d(x) approximately constant, the capacitance between the layers for that area is where ϵ is the electrical permittivity of the insulator between the layers.With signal angular frequency ω, this leads to conductance where j is the imaginary unit.From this, we obtain In this article, we neglect the phase of the coupling, i.e., we use only the absolute value of the coupling computed by (12).Due to this approximation, the reconstructed couplings are to be interpreted only qualitatively.
The spatial coordinate system used for the two conductive surfaces consists of the two coordinates of the point along the surface and an additional index distinguishing the two layers.Therefore, the domain for potential u and conductivity σ is Ω ⊂ R 2 × {1, 2}.On the other hand, the coupling between the layers is shared by the two surfaces, and thus has the domain Ω c ⊂ R 2 .We approximate the areas under the attached electrodes to remain in the same potential as the electrode, and therefore do not consider these areas as a part of the domain Ω, as the local potential u is not solved from the CEM there.However, the area under an electrode on the surface that the electrode is not attached on is part of domain Ω.
We notate the area where neither of the surfaces has electrodes as Ω 0 , i.e., Ω 0 = {(x, i) ∈ Ω | (x, j) ∈ Ω ∀ j = 1, 2}.In addition, the area under electrode l on the surface that the electrode is not attached to is notated as Ω e l .Similarly as before, the boundary corresponding to electrode l is denoted by ∂Ω e l .The subdomains are illustrated in Figure 2. It is Adding the coupling into the CEM requires modifying Equation (1), which is the conservation of currents in the continuum.In the coupled double-layer system, the current on each layer does not have to be conserved, but instead, any current sink on one layer has to be matched by a current source on the other layer, and vice versa.Thus, the total current is conserved, although, on the 2D surfaces, it is not.
Adding these sinks and sources to Equation (1) leads to the set of equations where u(x, i) and σ(x, i) are the potential and conductivity, respectively, at point (x, i) and c(x) is the coupling at point x.The coupling c(x), together with the potential difference in the two layers, determines the current density flowing between the surfaces.In Equations ( 13) and ( 14), we consider only the area without electrodes on either surface.Wherever one of the surfaces has an electrode, the potential function u is not defined, and instead we use the potential of the electrode where Ω e l is the area under electrode l on the surface that the electrode is not attached to.
Equations ( 13)-( 15) can be combined into a single equation by first defining a function for the potential on the other layer Hence, we can write the CEM extended for the coupled double-layer sensing skin: Compared to the original CEM, Equation ( 19) is also modified to take into account the current sources and sinks on Ω e l that contribute to the measured currents as well.

Forward Model Solution
In this paper, we express the coupling c(x) via a separate mesh covering the domain Ω c .The coupling has values c i at nodes x i of this mesh.The function c(x) is then approximated in the FEM equations by where the coefficients c ′ i are interpolated linearly from values c i based on the two meshes.Further, we notate that The finite element formulation for the conventional CEM is well known, but in order to take into account the modifications made to model the double-layer case, we show the full derivation in Appendix A. We derive the FEM for the potential injection case, i.e., we solve the currents from the CEM, similarly as in [43].We emphasize that the additions presented here work with few changes also using the more conventional case of solving the potentials from the CEM.We denote the computational model of the measurements obtained from FEM by Ĩ( σ, c, z), which is a mapping from the parameter space to the vectors containing currents corresponding to all used potential excitations stacked into a single vector.
The basic principles of finding the coefficients of coupling c based on the measurements are the same as in Section 2.2.However, as the conductivity distributions of the surfaces and the contact impedances of the electrodes are not known, they are also estimated.Hence, we define a generalized estimate variable Θ, i.e., any one of c, σ, z or a vector obtained by stacking any combination of these.The conductivity estimate σ in the equations presented in Section 2.2 is then replaced by Θ.The reconstruction process and the constituents of Θ in each step are described in more detail in Section 4.3.For details on the simultaneous estimation of contact impedances, the reader is referred to [47,48].

Prior Functions
To mitigate the ill-posedness of the inverse problem, prior knowledge of the estimated distributions is used.In this study, we assume the conductivity distributions to have spatially smooth variations throughout the surfaces caused by inhomogeneities of the paint layers.On the other hand, the coupling is assumed to have a sparsely steep gradient caused by either compression of the insulating region under pressure or by water applied in between the two layers.Furthermore, both conductivity and the coupling are known to always be positive.
To promote smooth variations in the conductivity distribution, a Gaussian smoothness prior is used [49].It has the functional form where σ e is the expected value of σ and Γ sm determines the covariance between conductivities of all pairs of nodes in the FE mesh of each layer of the sensor.Similarly to [49], we use a Gaussian distance covariance, i.e., where x i , x j ∈ R are the positions of nodes i and j, a sm is the strength of the prior, b sm is a constant related to the length scale of the spatial smoothness, c sm is a small stabilizing constant and δ ij is the Kronecker delta.Throughout this study, we use c sm = 10 −4 × a sm .
To facilitate choosing the constant b sm , we define correlation length d cor as the length where the correlation of two nodes given by ( 25) is 0.01.Then, for the constant b sm , we have The total variation prior is used for favouring sparse spatial changes in the coupling.The total variation prior functional is [ where A k is the area of element k, n is the number of elements in the mesh, a TV is a strength parameter and β is a smoothing parameter used to make the functional differentiable.D x and D y are discrete differential operators corresponding to differentiation in x-and y-directions.The strength parameter a TV is chosen based on the procedure outlined in [50].
In short, where p is the probability of having the maximal gradient norm value of g max on the domain.
The positivity of estimates is enforced by a parabolic barrier function where a pos is the strength parameter for the constraint, Θ is the generalized estimate variable and Θ m is a vector containing minimum estimate values that are not penalized.
Here, we allow for different minimum values for different elements of Θ because the scales of σ, c and z, any of which may be in Θ, may be different.

Simulation Setup
We use numerical simulations to study the overall feasibility of the proposed method.The simulations were computed using the finite element (FE) approximation of the coupled CEM detailed in Section 3.1.The imaging domain was chosen to correspond to the experimental setup, i.e., two square surfaces with side lengths of 25 cm and 16 electrodes on each layer.Measurement setups were simulated corresponding to both pressure and leak detection.For all targets, three different measurements were simulated: 1.
Reference measurements of the layers separately, i.e., with zero coupling between the layers.

2.
Reference measurements of the assembled sensing skin, i.e., with homogeneous coupling.

3.
Measurements with the effect (i.e., applied pressure or water) in place.
Gaussian noise with a standard deviation of 0.1% of each measurement value was added to all measurements.
In the simulations, the insulating layer had a thickness of 5 mm at rest and relative permittivity ϵ r = 4.The couplings for simulations were computed using Equation ( 12), except for the water ingress.
For pressure sensing, two targets (targets I and II) were simulated, both modeling pressing the distance between the layers from 5 mm down to 2 mm on a circular region having a diameter of 5 cm.Target I was simulated with the sensing skin having homogeneous conductivity of 2×10 −3 Ω −1 , while target II was simulated with the sensing skin having the conductivity distribution drawn from the smoothness prior distribution Equation (24) with parameter values d cor = 10 cm, a sm = 1.6 × 10 −7 S 2 and σ e a vector containing values 2 × 10 −3 S in each element.For leak detection, water ingress between the layers was simulated by an area of increased coupling (target III).Equation (12) was used for the background coupling, and the peak value on the wet area was chosen to have ten-fold coupling compared to the background.The initial conductivity distribution for this simulation was drawn from the same distribution as for target II but featured an additional increase in conductivity on the region of water ingress.For target I, contact impedances were set to a negligible value of 10 −6 Ω.For targets II and III, the contact impedances were drawn from normal distribution N (1, 1) Ω, with all values less than 10 −6 Ω replaced by the value 10 −6 Ω.
A mesh of 47,639 nodes and 92,490 elements was used for the FEM simulations.The element widths varied between 0.5 mm (near the electrodes) and 16 mm (on the central regions).The conductivity was represented in a mesh of 4122 nodes and 7616 elements (element widths ranging from 3 mm to 19 mm) in total for the two surfaces, and the coupling was represented in a mesh of 790 nodes and 1478 elements (element widths approximately 10 mm).The conductivity and coupling values were linearly interpolated to the nodal points of the FE mesh.

Experimental Setup
For measurements, a commercial potential excitation EIT device manufactured by Rocsole Ltd. (based in Kuopio, Finland) was used.The injection frequency was chosen as 320 kHz and single-electrode excitations were used.
The two conductive surfaces were made of rubber sheets painted with conductive paint.The paint was a mixture of rubber coating (Maston RUBBERcomp) and graphite powder (CRETACOLOR) added for conductivity in a mix ratio of 1:10 by weight.Electrodes were manufactured of copper tape, which was attached to the rubber sheet with the pre-applied adhesive before applying the paint.The pieces of copper tape were longer than the modeled electrode so that the rest of the piece of tape was turned on the back side of the rubber sheet.Short pieces of wire, with SMB connectors on the other end, were soldered on the copper tapes on the back side.The conductive paint was applied by brush so that the layer of paint covers the entire front surface of each rubber sheet, including the parts of the attached electrodes on the front side.Two layers of conductive paint were applied to obtain a more even coating.The first layer was left to dry overnight before applying the second layer of paint.The insulating layer between the conductive surfaces was made of polyester padding when sensing pressure, and polymer foam when sensing water.The sensing skin was assembled by placing the insulating layer between the rubber sheets, both of the rubber sheets having their front sides, i.e., the conductive surfaces, facing the insulating layer.The layers were not attached together but were only placed resting on top of each other in order to facilitate reusability with different insulating materials.
We note that the choice of the materials used in this experiment is not necessarily optimal for the final applications, but they were chosen to facilitate the laboratory experiments.For example, instead of using rubber sheets for the conductive layers for leak detection, a water-transparent material, such as conductive fabric (e.g., the one used in [9]), may be desirable.Otherwise, a rupture in one of the conductive layers is required for the water to enter the insulating layer and be detected.
A similar set of three measurements was taken as described in Section 4.1 for the simulations, except measurement set 1 (layers measured separately), which was taken only once at the beginning and used as the reference data for all the targets.The second set of measurements (assembled sensing skin without pressure or water) was taken separately for each target before applying the effect.
Weights of 11.6 kg and 3.4 kg were placed separately on the sensing skin for the measurements.Target IV was a 11.6 kg weight standing on the sensing skin on a circular surface having a diameter of 4.5 cm.Target V was a 20 cm tall circular cylinder (diameter 10 cm) weighing 3.4 kg resting on its side (see Figure 3).For target VI, 10 ml of water was injected in between the layers for the leak detection measurements.

Reconstruction
The reference measurements can be used in several ways to mitigate some of the modeling errors present in the system.The reconstruction process used in this study is detailed in this section.During the whole process, the noise of the forward model is assumed to consist of two additive Gaussian components.The first component, a baseline noise, has the same standard deviation s 1 for each measured value, with s 1 being 0.01% of the difference in maximum and minimum measured current.The second component has a standard deviation of 3% of each measured value.We note that the actual measured noise is well below this level, but these values try to approximately include the possible modeling errors as well.For further mitigating the effects of modeling errors, the Bayesian approximation error approach [51,52] could be applied, but it is left outside the scope of this paper.
To use the reference data, solving the inverse problem is divided into three steps: 1.
From measurements made of the two layers separately, their initial conductivity distributions σ0 and the contact impedances z0 of electrodes are estimated as in singlelayer EIT described in Section 2.

2.
From measurements, where the layers were combined, but no pressure or water had yet been applied, a homogeneous estimate c 0 of the coupling is computed, using the conductivity distribution σ0 and contact impedances z0 estimated in step 1.Using the estimate c 0 , a correction to the measurements is computed as specified in reference [24].

3.
From measurements after applying pressure or water, the coupling distribution is estimated together with re-estimation of the conductivity distribution.
The conductivity distribution is re-estimated at the final step in addition to the coupling distribution since applying pressure or water is also presumed to change the conductivity distribution.The estimations in steps 1 and 3 are performed using the GN algorithm outlined in Section 2.2.In step 2, the estimation is simple, as the homogeneous estimate is characterized by a single parameter.The estimation of contact impedances is performed as presented in [47,48].
In step 3, the generalized estimate variable Θ = σT zT T and the prior potential function is Here, the parameters used are a sm = (0.01 × 1/n ∑ n i=1 σ0,i ) 2 , d cor = 5 cm, σ e = σ0 , β = 10 −4 .For reconstructions corresponding to measured data, the parameter a TV is chosen by Equation (28), with p = 0.975 and g max = c0 /d elem , where d elem is the average edge length of mesh elements.For reconstructions corresponding to simulated data, the same procedure is used, but the parameter value is multiplied by an additional constant of 10 −4 for pressure-sensing simulations and 10 −2 for leak detection simulations.The positivity constraint parameters are Θ m,k = 10 −4 S and a pos = 10 5 S −2 for conductivity values, and Θ m,k = 10 −2 Sm −2 and a pos = 10 4 S −2 m 4 for coupling values.
The forward problem is solved by FEM as described in Section 3.1 using a mesh with 4122 nodes and 7616 elements (element widths ranging from 3 mm to 19 mm).The coupling is defined in a mesh containing 790 nodes and 1478 elements (element widths approximately 10 mm).

Simulations
The reconstructions computed from the pressure detection simulations (targets I and II) are shown in Figure 4.Both reconstructed locations and shapes of the inclusions are consistent with the targets within an accuracy expected from an EIT-based method.For target I, the value of coupling inside the inclusion (ca.4.5 × 10 −3 Sm −2 ) does not agree with the true value (ca.5.6 × 10 −3 Sm −2 ).The discrepancy of reconstructed and true values is likely caused by the total variation prior, as it penalizes stronger variations more.The homogeneous conductivity distribution of target I is reconstructed well, with the average conductivity of the reconstruction close to the true value of 2 × 10 −3 S and variations in the order of 1% in the conductivity distribution of the reconstruction.For target II, the reconstructed coupling value inside the inclusion (ca.5.2 × 10 −3 Sm −2 ) agrees with the true value (ca.5.6 × 10 −3 Sm −2 ) better than for target I, but the shape of the inclusion is more distorted.The distortion of the shape may be caused by the proximity of the boundary and electrodes, increasing the modeling error caused by using a different FE mesh for the reconstruction than that used for the data simulation.The reconstructed conductivity distribution of target II shows all the features of the target within the expected accuracy.These reconstructions show that, in the absence of large modeling errors, the size and location of single areas of pressure can be reconstructed fairly accurately even with inhomogeneous background conductivity of the sensing skin.The robustness of quantifying the value of coupling inside an inclusion has to be improved for applications requiring the absolute value, e.g., for quantifying the applied pressure on a surface.For many applications, however, the location of the applied pressure is enough, and, for those applications, the proposed method is viable.
In Figure 5, the reconstructions corresponding to water ingress simulation (target III) are shown.The area of water ingress is depicted well in the coupling reconstruction, with the numerical value (ca.2.1 × 10 −2 Sm −2 ) being quite close to correct (ca.2.4 × 10 −2 Sm −2 ).In addition, the conductivities are reconstructed well within the accuracy expected, with also the increased conductivity in the area of water ingress visible.
This evidence supports the feasibility of leak detection by the proposed double-layer sensing skin.The minor discrepancies between the reconstructed coupling and the true distribution are irrelevant for the application of leak detection.

Experimental Results
The experimental results of pressure detection (targets IV and V) are shown in Figure 6.In this figure, only the reconstructions of the coupling distribution are represented.Both reconstructions show an increased coupling on or close to the correct area.Both reconstructions feature some artifacts of decreased coupling as well.As the applied pressure should not cause a decrease in coupling, these artifacts are easily distinguished from the locations of actual applied pressure.
The reconstruction of target IV in Figure 6a shows the elongated shape of the pressure distribution well.However, the area of increased coupling is reconstructed as continuing through the sensing skin, from one boundary to another, whereas the edges of the object resting on the surface are 2.5 cm from both boundaries.This error is likely caused by the surface deforming further than the exact location of applied pressure.From the reconstruction of target V in Figure 6b, the size and shape of the inclusion cannot be inferred.The spreading of the increased coupling to a larger area is likely caused by the surface deforming on a larger area than the applied pressure.This deformation may also include wrinkling of the surface, which may explain some artifacts even having decreased electrical coupling.Another possible cause of significant error is the alteration in contact impedances caused by the stretching and wrinkling of the conductive sheets under the electrodes.
Although the proposed method seems prone to having some artifacts and inaccuracies in the location and shape of the observed distribution, the results support the feasibility of detecting the location of applied pressure using the proposed method.In this study, we used weights in the order of several kilograms, but we emphasize that the weight range that the sensing skin is sensitive to can be tuned by choosing the materials of conductive and insulating layers.
In Figure 7, the reconstructed coupling distribution corresponding to the water ingress (target VI) is shown.The reconstruction features increased coupling coinciding well with the approximate location of the water ingress.However, in addition, there is another region of increased coupling with no clear explanation.This may be caused, for example, by variations in contact impedances or by slight bending of the edges of the conductive layers.Despite the artifact, the location of water ingress is clearly visible, and thus the results support the viability of using the proposed method for leak monitoring.In this appendix, we present the derivation of the finite element method formulation 439 for the coupled double layer sensing skin.The derivation follows along the lines of previous 440

Conclusions
In this article, we studied feasibility of using EIT-based electrically coupled doublelayer sensing skin for pressure sensing and leak detection.The aim of the method is to reconstruct the spatial distribution of the electrical coupling between the two layers.This distribution contains information of possible external pressure pushing the layers closer to each other, or a conductive path of water formed between the layers.
The feasibility of the proposed method for pressure sensing and leak detection was tested on both simulated and measured data.The reconstructions corresponding to simulated data localized both the applied pressure and the ingress of water, with slight distortions in the shape or the reconstructed electrical coupling value.The reconstructions from measured data had some artifacts, but the approximate locations of the applied pressure or water leak could still be inferred from the reconstructions.While the reconstruction from a circular pressure distribution failed to show the shape of the applied pressure, the reconstruction of a line-like target showed the elongated shape of the applied pressure.
The results of both the simulations and the experimental studies support the feasibility of the proposed method for pressure sensing and leak detection.Hence, the proposed method may be used in applications such as tactile sensing in robotics and leak monitoring in buried structures.and thus we have all parts required for calculating the Jacobian with respect to coupling c.

Figure 1 .
Figure 1.(a) A schematic diagram showing the layers of the sensing skin and locations of electrodes.The electrodes are drawn in copper brown.(b) The two conductive layers of the sensing skin without the insulating middle layer.The black side of the layer is the conductive paint and the yellow side shows the rubber base of the layer.(c) The assembled sensing skin with the insulating layer in the middle.

1 Figure 2 .
Figure2.A schematic diagram of the imaging domain of the double-layer sensing skin.Ω 0 is the area common for both layers.Ω e l is the area under electrode e l on the surface that the electrode is not attached to and ∂Ω e l is the domain of the boundary corresponding to electrode e l .Electrodes are drawn in copper color on the layer that they are attached to.Labels for electrodes 1 and 2 are shown as examples, with blue highlights used to differentiate where the boundary and where the area of the electrode are meant.

Figure 3 .
(a) Target IV: 11.6 kg weight and (b) target V: 3.4 kg weight resting on the sensing skin during measurements.

Figure 4 .
Figure 4. (a) The coupling distribution and (b) the conductivity distributions of target I, and (c) The coupling distribution and (d) the conductivity distributions of target II.First row: the true values.Second row: reconstructions.The conductivity distributions of the two layers are shown side by side.Couplings are in units of Sm −2 and conductivities in units of S.

Figure 5 .
Figure 5. (a) The coupling distribution and (b) the conductivity distributions of target III.First row: the true values.Second row: reconstructions.The conductivity distributions of the two layers are shown side by side.Couplings are in units of Sm −2 and conductivities in units of S.

Figure 6 . 20 Figure 7 .
Figure 6.Reconstructed coupling distributions computed from measured data.(a) The black line corresponds to the central line of contact of the 3.4 kg cylinder resting on its side.(b) The black circle corresponds to the bottom surface of the 11.6 kg weight.The units are in Sm −2 .Version June 20, 2024 submitted to Sensors 14 of 20

Figure 7 .
Figure 7. Reconstructed coupling distribution computed from measured data corresponding to water injected between the layers.The black circle marks the approximate location of water injection.The units are in Sm −2 .