Next Article in Journal
A Deep Learning Based Data Recovery Approach for Missing and Erroneous Data of IoT Nodes
Previous Article in Journal
A Hybrid Prognostic Method for Proton-Exchange-Membrane Fuel Cell with Decomposition Forecasting Framework Based on AEKF and LSTM
Previous Article in Special Issue
Eddy Current Measurement for Planar Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microwave-Based Subsurface Characterization through a Combined Finite Element and Variable Exponent Spaces Technique

by
Valentina Schenone
1,
Claudio Estatico
2,
Gian Luigi Gragnani
1,*,
Matteo Pastorino
1,
Andrea Randazzo
1 and
Alessandro Fedeli
1,*
1
Department of Electrical, Electronic, Telecommunications Engineering, and Naval Architecture, University of Genoa, 16145 Genoa, Italy
2
Department of Mathematics, University of Genoa, 16146 Genoa, Italy
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(1), 167; https://doi.org/10.3390/s23010167
Submission received: 29 October 2022 / Revised: 20 December 2022 / Accepted: 21 December 2022 / Published: 24 December 2022
(This article belongs to the Special Issue Electromagnetic Structures for Sensing Applications)

Abstract

:
A microwave characterization technique to inspect subsurface scenarios is proposed and numerically assessed in this paper. The approach is based on a combination of finite element electromagnetic modeling and an inversion procedure in Lebesgue spaces with variable exponents. The former allows for description of the measurement system and subsurface scenario with high accuracy, while the latter exploits the adaptive definition of exponent function to achieve improved results in the regularized solution of the inverse scattering problem. The method has been assessed with numerical simulations regarding two-layered environments with both planar and non-planar air–soil interfaces. The results show the capabilities of the method of detecting buried objects in different operative conditions.

1. Introduction

The ability to investigate subsurface scenarios is needed for a variety of applications, including geophysical reconnaissance, archaeology, mine location and detection, environmental monitoring, and ground mapping in civil applications [1,2,3,4,5,6,7,8]. In this area, ground penetrating radars (GPRs) are one of the most widely adopted tools. However, GPR data are usually represented through B-scans, which are a time-domain representation of GPR measurements. Even though some post-processing techniques can be applied to these data, such as migration techniques [1,9], classical GPRs often require experienced users to interpret data. Moreover, they only provide a qualitative reconstruction of subsurface scenarios, without allowing for a precise dielectric characterization of the inspected targets [6].
Among the various possible approaches to inspect underground structures, microwave imaging provides the ability to characterize subsoil regions based on scattering measurements collected by antennas placed over the area of interest. These techniques can be used to non-invasively determine the presence and the approximate shape of buried objects (in the case of qualitative methods [10,11,12,13,14,15,16]) or the distribution of their dielectric properties (in the case of quantitative strategies [17,18,19,20,21,22,23,24,25]).
The scenario in which prospecting is performed makes the imaging problem particularly challenging, especially focusing on quantitative inverse scattering methods. Indeed, in addition to the general issues of the inverse scattering problem, which is inherently ill-posed and nonlinear [4], there are other issues that are typical of this configuration. First, unlike other applications where it is possible to place antennas around the entire survey area, here, antennas can only be placed on the survey line above the ground or in boreholes, reducing the available information to solve the problem [26,27]. In addition, the layering of the environment is another critical element; for instance, probes are positioned at a certain distance from the ground, and the presence of discontinuities in the dielectric properties of the scattering environment further adds a level of complexity due to the insurgence of reflection phenomena [28].
In order to reconstruct accurate images of objects in the subsurface environments, imaging techniques that integrate electromagnetic models in these settings have been proposed in the literature. Indeed, the ability to accurately model the real environment in which the diagnosis is made is crucial to obtain good results [29]. Moreover, the exploitation of a suitable model to solve the forward problem enables synthetic data to be obtained and testing methods in a simulated environment. In more detail, two main classes of mathematical models can be distinguished: analytical and numerical ones. The former are characterized by good robustness, but assuming simplified hypotheses [30,31,32,33]. The second class includes methods based on finite differences in the time domain, the finite element method, and the method of moments [34,35,36]. These classes of approaches, although more computationally expensive, are particularly advantageous when certain geometric configurations and specific soils are present in the domain. In addition, several radar-like techniques [37,38,39,40] have been recently proposed in this research area to allow soil profiling, so that retrieved soil configuration can be accurately integrated into the electromagnetic model as a priori information [41]. Finally, hybrid solutions combining analytical and numerical models have also been proposed [11].
In this paper, a quantitative imaging method based on the inversion of scattering S -parameters combined with an electromagnetic model of finite element (FE) type is proposed and applied to the diagnosis of shallow subsurface structures. In the literature, several approaches have been described to solve the nonlinear inverse scattering problem by deterministic, stochastic, and neural networks techniques [42,43,44,45,46,47]. Among the deterministic approaches, an interesting class is represented by Newton-type schemes [48,49,50,51]. In the present work, the nonlinear and ill-posed inversion problem is treated using a Newton-type technique formulated in L p ( ) Lebesgue spaces with variable exponent [52,53,54]. This is an iterative regularization procedure which is able to provide accurate diagnostic results due to an adaptive definition of the exponent function. Moreover, the FE-based electromagnetic model combined with the formulation of S -parameters [8] allows for a precise modelling of the electromagnetic problem. The approach, which was first developed for stroke imaging [54], has now been expanded to buried target detection. The main plus of the proposed approach for shallow subsurface inspection is the capability of suitably taking into account the measurements and background configuration inside the inversion procedure. Indeed, in this way, the structure of the measurement system and environments can be considered in the electromagnetic model providing valuable reconstruction results even in this challenging scenario.
The method has been validated by means of numerical data obtained with a two-layer configuration illuminated by open waveguide probes. In detail, the variation in target size, depth, and some other relevant parameters has been evaluated. Finally, a study of the method in the case of non-planar air-soil interface has been conducted; the results of inversion with a priori knowledge of the exact interface have been compared with those achieved by replacing the actual surface profile with a planar one in the embedded model.
The paper is organized as follows: in Section 2, the inverse scattering formulation along with the finite element approach and the inversion procedure are reported. The results are then presented in Section 3. Finally, the Conclusions follow in Section 4.

2. Formulation of the Problem

The shallow subsurface scenario in which detection is performed is shown in Figure 1. A multistatic and multiview system with W waveguide antennas positioned at height h from air–soil interface Λ on a measurement line of length l m parallel to x -axis is used to illuminate the soil. The air is modelled as vacuum (i.e., ε 0 8.85 × 10 12 F/m) and the soil is characterized by a complex dielectric permittivity ε b .
The antennas are activated one by one, and measurements of transmission S-parameters are completed between the antenna ports. In addition, z -polarized fields with no components of the electric field parallel to the x y plane at an angular frequency ω are considered, and the target properties are hypothesized invariant along the z axis. Therefore, a two-dimensional environment is assumed [55].
The distribution of the dielectric permittivity of the region R containing the cross-section of the objects under test is retrieved starting from scattering S -parameters between antenna ports i and w , δ S ^ i w = S t o t i w S i n c i w , where S t o t i w is the S -parameter measured in the presence of the target and S i n c i w is the corresponding quantity related to an empty investigation domain (assumed to be available or estimated).
Moreover, the dielectric properties of the investigation domain R can be described by the contrast function
τ ( x , y ) = ε ( x , y ) ε b ε b
where ε ( x , y ) is the complex dielectric permittivity of the configuration under test and ε b is the background value. In a two-dimensional setting, the scattering parameters are related to the contrast function by means of the following data equation [8]
j ω ε b 2 c i c w R E i n c i ( x , y ) E t o t w ( x , y ) τ ( x , y ) d x d y = δ S ^ i w
where c m is the incoming wave amplitude on m -th antenna port, E i n c i is the z -component of the incident electric field when the i -th antenna acts as a transmitter and E t o t w is the z -component of the total electric field when the w -th antenna transmits, i w . In detail, E t o t w represents the field collected in the presence of the unknown target; therefore, it is a function of the unknown contrast function itself, and this results in a nonlinear relation linking the contrast function with the scattering parameters.
Considering the measurements collected by each antenna for each view, it follows:
[ j ω ε b 2 c 1 c 2 R E i n c 1 ( x , y ) E t o t 2 ( x , y ) τ ( x , y ) d x d y j ω ε b 2 c 1 c W R E i n c 1 ( x , y ) E t o t W ( x , y ) τ ( x , y ) d x d y j ω ε b 2 c W c W 1 R E i n c W ( x , y ) E t o t W 1 ( x , y ) τ ( x , y ) d x d y ] = [ δ S ^ 12 δ S ^ 1 W δ S ^ W ( W 1 ) ]
This is the nonlinear system to be inverted starting from scattering parameters measurements to retrieve the contrast function.
To handle this problem, an FE approach is applied and integrated inside the inversion performed in variable exponent Lebesgue spaces. In detail, the electromagnetic model expressed through the FE method as applied to the present inverse scattering problem is described in Section 2.1. Then, in Section 2.2, the inversion approach is presented.

2.1. FE Approach

A fundamental step for the development of the proposed imaging method is to define a procedure to model the electromagnetic problem, exploiting an FE formulation to compute the electric fields and S -parameters inside the inversion procedure.
In more detail, the analyzed measurement system, shown in Figure 2, is composed of a set of W open waveguides of width a and length b , filled with material of complex dielectric permittivity ε w g . The waveguides are terminated by PEC boundaries on the lateral sides, Π P E C , and a waveguide port on the top side, Π w . In each antenna, a reference system ( x ( w ) , y ( w ) ) is defined at the waveguide port centered at
( x ( w ) , y ( w ) ) = ( l m a 2 + ( w 1 ) ( s + a ) ,   h + b ) , w = 1 , , W
where s is the mutual distance between antenna waveguides.
Assuming the i -th waveguide excited with an incoming wave of the T E 10 fundamental mode, to compute the electric field, the Helmholtz equation solution is necessary along with boundary conditions in order to take into account the exact structure of the environment and measurement system:
t 2 E i n c / t o t i ( x , y ) + ω 2 μ o ε ( x , y ) E i n c / t o t i ( x , y ) = 0
where non-magnetic materials are considered. The following boundary conditions must be imposed: on Π P E C , Dirichlet conditions state that E i n c / t o t i | Π P E C = 0 ; on waveguide ports, the Dirichlet and Neumann conditions are E i n c / t o t i | Π w = E w g ( w , i ) and E i n c / t o t i / y | Π w = E w g , i n c / t o t ( w , i ) / y ( w ) , with w = 1 , , W where E w g , i n c / t o t ( w , i ) , i.e., the z -component of the electric field tangent to the w -th port when the i -th port is in transmitting mode, is:
E w g , i n c / t o t ( i , w ) ( x ( w ) , y ( w ) ) = δ i w e 1 ( i ) ( x ( i ) ) e j β 1 y ( i ) + m = 1 M A m , i n c / t o t ( i , w ) e m ( w ) ( x ( w ) ) e j β m y ( w )
with δ i w the Kronecker delta, A m , i n c / t o t ( i , w ) the amplitude of m -th mode in w th port, e m ( w ) the orthonormal modal function of T E m 0 modes in the w th waveguide, and β m the propagation constant of T E m 0 modes inside the waveguide [56].
The first term of Equation (6) represents the incoming T E 10 mode feeding the i -th port and the second one represents the outcoming T E m 0 modes inside the w th waveguide. For the solution of the forward problem, coefficients A 1 , i n c / t o t ( i , w ) are left as unknowns, and the S -parameters are then retrieved as S i n c / t o t i w = A 1 , i n c / t o t ( i , w ) . Moreover, in the analysed configuration, absorbing boundary conditions (ABC) are required to terminate the simulation domain outside the antenna ports. To this end, the perfectly matched anisotropic absorber (PMA) has been selected [56,57].
The simulation domain includes the above-described shallow subsurface scenario comprised air and soil layers and limited by waveguide ports, waveguide PEC walls, and the PMA layer. In order to solve the forward problem with FE formulation, the domain has been partitioned in a mesh of N triangles of dimension Σ n with frontier Π ( n ) . Then, the electric field E i n c / t o t i ( x , y ) in each n -th element of the mesh is formulated through first-order basis functions. It is worth highlighting that such a triangular form of subdomains gives the possibility of mapping complex structures in an accurate way, making it suitable to describe the problem at hand. In particular, the first-order triangular basis functions for each n -th triangle, Ψ t ( n ) , satisfy properties Ψ t ( n ) ( x u ( n ) , y u ( n ) ) = δ u t and t = 1 3 Ψ t ( n ) ( x , y ) = 1 ,   ( x , y ) Σ n , where ( x u ( n ) , y u ( n ) ) are the coordinates of the u -th node of each n -th triangular element and are defined as follows:
[ Ψ 1 ( n ) ( x , y ) Ψ 2 ( n ) ( x , y ) Ψ 3 ( n ) ( x , y ) ] = [ 1 1 1 x 1 ( n ) x 2 ( n ) x 3 ( n ) y 1 ( n ) y 2 ( n ) y 3 ( n ) ] 1 [ 1 x y ] , ( x , y ) Σ n
The electric field in each n -th triangle can be written as:
E i n c / t o t i ( x , y ) = t = 1 E E i n c / t o t , t i ( n ) Ψ t ( n ) ( x , y ) ,     ( x , y ) Σ n
By considering the FE approximation of field along with the Helmholtz equation, the electric field is numerically computed following the approach in [56].
Then, in the soil region, a rectangular investigation domain R of dimension l R × h R centered at ( x R , y R ) has been defined (Figure 3), which, in its discrete representation, is a subregion of the simulation domain composed of N R triangles. Therefore, we can introduce the vector ε _ [ ε ( 1 ) , , ε ( N R ) ] T , which is obtained by approximating the dielectric properties as constant inside each triangle of the domain R , thus, the corresponding contrast vector results τ _ [ τ ( 1 ) , , τ ( N R ) ] T .
In this way, starting from Equation (2), considering the FE formulation, we obtain the following system:
D _ ( τ _ ) [ j ω ε b 2 b 1 b 2 n = 1 N R E _ i n c 1 ( n ) T [ T ( n ) ] E _ t o t 2 ( n ) τ ( n ) j ω ε b 2 b 1 b W n = 1 N R E _ i n c 1 ( n ) T [ T ( n ) ] E _ t o t W ( n ) τ ( n ) j ω ε b 2 b W b W 1 n = 1 N R E _ W , i n c ( n ) T [ T ( n ) ] E _ t o t W ( n ) τ ( n ) ] = [ δ S ^ 12 δ S ^ 1 W δ S ^ W ( W 1 ) ] δ S ^ _
where E _ i n c / t o t i ( n ) = [ E 1 , i n c / t o t i ( n ) , E 2 , i n c / t o t i ( n ) , E 3 , i n c / t o t i ( n ) ] T is the vector with the nodal values coefficients of the field in each triangle computed by means of the FE method and [ T ( n ) ] is a 3 × 3 matrix of coefficients:
T u t ( n ) = Σ n Ψ u ( n ) Ψ t ( n ) d x d y , u , t = 1 , , 3
By means of Equation (8), a discrete nonlinear operator D _ ( τ _ ) linking the contrast vector with the measurements in δ S ^ _ is defined.

2.2. Inversion Procedure

The imaging problem is solved by inverting Equation (9). In order to address the inversion, an inexact Newton procedure has been exploited where the space X of the unknowns is a variable exponent Lebesgue space L p ( ) (i.e., τ _ X L p ( ) ) and the space Y of the data is a space L p a with constant exponent (i.e., δ S ^ _ Y L p a ) [53]. In particular, since the discrete case is considered, the exponent function of unknowns’ space X turns out to be a vector p _ = [ p ( 1 ) , , p ( N R ) ] T , p ( n ) being the value of the exponent in each n -th triangle. The constant exponent of data space Y is set equal to its spatial average value in R , i.e., p a = 1 N R n = 1 N R p ( n ) .
The inexact Newton method solves the inverse problem in a regularized way by minimizing the following residual function:
Ω ( τ _ ) = 1 2 δ S ^ _ D _ ( τ _ ) Y 2
with · Y norm in the space Y , by moving along a nonstandard gradient direction in the dual space X * of X . In detail, the method consists of two nested iteration cycles as summarized in the flow chart in Figure 4 [54]. At first, at each k -th iteration of the external cycle, Equation (9) is linearized around its current solution (denoted as τ _ k ) with null initial value τ _ 0 = 0 _ . Then, a Landweber-like approach in Lebesgue space L p ( ) is adopted to solve the resulting linearized problem (as shown in the red box of the flowchart). With reference to Figure 4, J X , J X * , and J Y represent the duality maps of X ,   X * , and Y , respectively. ( X * = L p ( ) * is the dual space of X , with p ( ) * as the Hölder conjugate of p ( ) ), D _ k is the Fréchet derivative of D _ , and α = D _ k 2 is the step width [52,54].
Moreover, the exponent function is adaptively modified during the inversion to achieve accurate reconstruction results. In the first iteration of the external cycle, the exponent vector is set to constant values p _ 0 = [ p s t a r t , , p s t a r t ] T since no a priori information is available. Then, at each outer iteration, the current solution τ _ k is exploited for updating p _ , with the l -th component
p k + 1 ( l ) = p m i n + ( p M A X p m i n ) τ k ( l ) / m a x n ϵ N R | τ k ( n ) | , 1 < p m i n p M A X
where [ p m i n , p M A X ] denotes the range of p . This way, values of p close to p m i n are assigned to the portions of R without targets leading to a reduction in ringing effects; values of p close to the maximum are in target regions where a smooth reconstruction is favored.
Finally, the two loops are stopped when the proper convergence criteria are fulfilled. In this case, the number of inner, I m a x , and outer iterations, K m a x , and the relative variation of the residual between two consequent steps, Δ , are considered as stopping rules.

3. Numerical Validation

In this section, numerical validation of the proposed imaging method is presented. In particular, a shallow subsurface scenario has been simulated by the FE forward solver to produce synthetic data and a reference case is described in Section 3.1. Then, the behavior of the method was investigated when some variations in the various parameters were applied. In particular, in Section 3.2 and Section 3.3, the effects and limitations on the performance of the method versus a variation in target size and target depth have been analysed, respectively. Then, in Section 3.4, the method has been studied considering different geometries and the number of targets, while in Section 3.5 the effect of the background uncertainties has been explored. Finally, in Section 3.6, some numerical tests in the presence of a non-planar air–soil interface are reported; the results achieved with exact interface surface known a priori or replaced with planar surfaces in the inversion procedure have been compared.

3.1. Reference Case

In numerical simulations, synthetic S-parameter data were provided as input to the proposed method to determine the distribution of complex dielectric permittivity within the investigation area. The FE forward solver has been exploited for the generation of numerical data.
Concerning the measurement parameters, a set of W = 10 antennas was placed along the measurement line of length l = 1.095   m at a distance s = 55 mm from eachother and located at h = 50 mm over a planar air–soil interface (Figure 1). The waveguides are filled with a material with a dielectric permittivity ε wg = 25 ε 0 and the dimensions a = 60 mm in width and b = 80 mm in length. The soil is characterized as dry sand with a permittivity ε b = ( 3 j 10 3 ) ε 0 .
Configurations with and without a target were simulated at frequency f = 550   MHz and M = 1 mode is considered. The simulation domain is discretized by Gmsh [58] using the frontal Delaunay algorithm with maximum edge length at the antenna ports and elsewhere with a size of s wg = 2 mm and s q = 4 mm, respectively. In this way, total and incident S-parameters have been generated. Total S-parameters were corrupted with a multiplicative Gaussian noise of 3%.
Once synthetic data have been generated, a rectangular investigation domain R of dimension l R = 60   cm and h R = 50 cm centered in ( x R , y R ) = ( 0 , 30 ) cm is considered. Inside R , a coarser discretization of edges s i n v = 7 mm is considered for the inversion.
The results have been evaluated using the following error metric:
e R e g = 1 N R e g n = 1 N R e g | ε ( n ) ε ¯ ( n ) | | ε b |
where R e g is the analysed region inside the investigation domain composed of N R e g elements, with R e g = { R , t a r } (i.e., whole domain R or target region t a r are considered), ε ( n ) is the reconstructed dielectric permittivity in the n -th element, and ε ¯ ( n ) is its reference value in the same element.
Initially, the target under consideration has a square cross-section of side l t = 15   cm , it is centered at ( x t , y t ) = ( 15 , 15 )   cm and characterized by a dielectric permittivity ε t = ε 0 . The method has been run under the following setting: K m a x = I m a x = 100 , p s t a r t = p m i n , and p M A X = 2 . Furthermore, solution loops are terminated when a threshold Δ = 1 % is reached. Concerning the range of variation in the exponent function, a study has been conducted to find its optimum value. Specifically, p m i n = [ 1.2 , 1.9 ] has been considered with step 0.1 .
Figure 5 shows the errors in the whole investigation domain, e R , and in the target region e t a r . As can be noticed, concerning the error in the investigation domain, an increment versus p m i n can be observed whereas the error inside the targets has a minimum in p m i n = 1.3 . Therefore, p m i n = 1.3 has been selected since it offers the best target reconstruction and high-quality results in the investigation domain. In Figure 6, the real part of relative dielectric permittivity is reported along with the actual target shape. As can be noticed, the localization of the target is adequate, and its dielectric permittivity has been reconstructed quite well.

3.2. Variation in Target Size

The size of the target has been changed in the first set of numerical simulations. Specifically, the values l t [ 3 , 18 ] cm with a step of 3 cm have been considered for the side length (i.e., l t has been varied between [ 0.095 λ b , 0.572 λ b ] with λ b the wavelength in the soil). Figure 7 shows the reconstructed distributions of the real part of the dielectric permittivity for l t = { 3 , 12 , 18 } cm. In configurations with l t = 18 cm and l t = 12 cm [Figure 7a,b], the reconstruction of the dielectric permittivity is close to the actual value. A good reconstruction is obtained, and the dielectric target is clearly visible and adequately localized. Conversely, considering a smaller target with l t = 3 cm, the object is barely visible in the reconstruction [Figure 7c]. In Figure 8, the reconstruction errors computed in the investigation domain and inside the target region are reported. As can be noticed, the whole domain error e R tends to increase with the target size. Instead, the error on the object e t a r has a parabolic-like behavior with a minimum at l t = 9 cm ( 0.286 λ b ). In particular, very small targets are weak scatterers and become difficult to detect. Indeed, an underestimation of the dielectric properties can be observed by reducing the target’s size until the method is not able to reconstruct the object. On the contrary, slight degradation of the background reconstruction can be seen by increasing the size of the target [Figure 7a].
In addition, the method was compared with the inversion formulated in the classical Hilbertian space (i.e., X L 2 ,   Y L 2 ). By comparing the errors achieved with the proposed variable exponent space method, it can be observed that the latter approach allows an accuracy improvement both in the investigation domain and in the target under test. This proves that the inversion performed in variable exponent Lebesgue spaces leads to more accurate reconstructions. Indeed, by defining a proper map of the exponent [Equation (12)], low values are assumed in the background and values close to p M A X are in the region where the inhomogeneities are localized. In this way, low values better control sparsity in the background whilst a better estimation of smoothness of the object is reached by assigning relatively high values of exponent in that region.

3.3. Variation in Target Depth

As a further analysis, the depth of the buried target has been varied, i.e., y t [ 15 , 35 ] cm with a step of 2.5 cm (i.e., y t [ 0.477 λ b , 1.112 λ b ] ) has been considered in order to assess the effectiveness and limitations of the method depending on the depth. The reconstructed distributions of the real part of the relative dielectric permittivity, R e { ε r } , are shown in Figure 9 for the cases y t = 32.5 cm, y t = 25 cm, and y t = 20 cm. Moreover, in Figure 10, the scattering S-parameters simulated when the first antenna acts as the source are compared with those computed inside the inversion procedure from the reconstructed dielectric permittivity for cases y t = 25 cm and y t = 20 cm. As can be observed, the reconstructed data match with the simulated ones.
By comparing the results, as the depth increases [Figure 6 and Figure 9a,b], the target reconstruction shows a decrease in the accuracy of the dielectric properties up to the point where it is difficult to locate the object [Figure 9c].
The reconstruction errors in the whole investigation domain and target region are shown in Figure 11. For low values of depth, the best target and whole domain error are achieved; then, by deepening the buried target, both the error in the target region and in the whole investigation domain has an upward trend. In more detail, beyond the depth of y t = 25 cm ( y t = 0.794 λ b ), the detection of the object is more difficult, and the quality of the reconstruction deteriorates as evidenced by the error trend, which for y t < 25 cm shows e R > 0.6 .

3.4. Effect of Different Geometries and a Different Number of Targets

In this section, the behavior of the method has been investigated considering different geometries and a different number of targets.
To this end, the target with a square-cross section of side l t = 12 cm analysed in the previous section has been considered and a second target with a circular-cross section has been introduced in the investigation domain. In detail, the second target is “vacuum-filled”, has a radius r = 6 cm, and is located at ( x t 2 , y t 2 ) = ( 17 , 20 )   cm .
Figure 12 shows the reconstructed distribution of the real part of the dielectric permittivity together with the actual shapes of the targets. As can be noticed, both targets with different geometries are correctly localized, and although in the case of the object with a circular-cross section a slight underestimation can be observed, a quite good reconstruction is achieved. Moreover, by comparing the reconstruction of the square-cross section target alone [Figure 7b] and in the presence of a second object, a slight deterioration can be observed, due to the interaction between the objects. This is confirmed by the computed reconstruction errors in the target region and in the whole domain (i.e., e t a r = 0.553 and e R = 0.074 ), which are higher than in the single-target case (i.e., e t a r = 0.521 and e R = 0.035 ).

3.5. Effect of Uncertainties in Soil Dielectric Properties

In order to assess the robustness of the method versus an inexact knowledge of the dielectric properties of the soil, the background dielectric permittivity used inside the inversion procedure has been set to ε r , b * = ε r , b ± 0.5 . All the other parameters are the same as those described in Section 3.1. In the first case, an underestimation of the background dielectric permittivity is considered whereas in the second case an overestimation is assumed. Figure 13 shows the reconstructed distributions of the real part of the dielectric permittivity. The reconstruction errors in the whole investigation domain and in the target region are shown in Table 1 (for the sake of comparison, we also reported the errors when the exact background properties are known in the inversion analysed in Section 3.1). The error in the whole investigation domain increases when an inaccurate background permittivity is considered. Moreover, a decrement of the target error can be noticed in the first case because the initial estimation of background permittivity is closer to the properties of the object.

3.6. Variation in Air–Soil Interface Roughness

A further study has been conducted to inspect the robustness of the method in the presence of a non-planar air–soil interface. Moreover, the results have been compared when:
  • the model endowed in the inversion takes into account the exact interface profile (i.e., assuming that the surface is known a priori);
  • a planar interface is considered instead (i.e., no a priori information is available).
In the analysed measurement setting, the antennas are located at h = 7 cm over the average soil level. To simulate a non-planar interface, in this set of tests, the profile of Λ is modelled with Catmull–Rom splines [59]. In particular, Λ is based on C = 7 equally spaced control points located to have a no-planar air-soil interface as shown in Figure 14. The root mean square height of the interface has been varied between h r m s = λ 10 and h r m s = λ 30 .
A “vacuum-filled” target (i.e., characterized by a permittivity ε t = ε 0 ) centered at ( x t , y t ) = ( 20 , 20 )   cm with l t = 12 cm has been investigated. The investigation domain is centered at ( x R , y R ) = ( 0 , 32 ) cm. All the other parameters of the forward and inverse procedure are the same as those in the previous sections.
In Figure 15, the reconstruction errors in the target region and in whole investigation domain are reported. Figure 16a–c shows the real part of the reconstructed dielectric permittivity when the surface profile is considered as a priori information in the inversion and Figure 16d–f reports the reconstruction when planar interfaces are adopted in the FE-model inside the inversion procedure.
As can be noticed from the reconstruction of the real part of the permittivity, although the reconstruction is quite successful in both cases, the artifacts in the background appear more pronounced for greater interface roughness when no a priori information about the interface is assumed in the inversion. This is confirmed by the trend of the errors; the error on the object is comparable in the two cases for h r m s λ 20 , while with h r m s = λ 10 , the knowledge of the exact interface brings a significant improvement in the reconstruction; the error on the whole domain is better in all cases with a valuable improvement as the roughness increases.
This proves the potentiality of the model embedded in the inversion procedure combined with S-parameters formulation in various operating conditions. Indeed, an accurate structural description of the involved environments can be integrated and taken into account in the inversion procedure, leading to an enhancement in reconstruction results.

4. Conclusions

In this paper, a microwave imaging technique for shallow subsurface prospecting is proposed, whose aim is to provide the distribution of dielectric permittivity from scattering S-parameters measurements. The method combines an FE approach with an inversion procedure in Lebesgue spaces with variable exponent. In this way, the structure of the measurement system together with the environment can be accurately considered together in the electromagnetic model by the FE formulation that is incorporated into the inversion procedure. Furthermore, the adopted inversion technique, which operates in non-Hilbertian Lebesgue spaces, exploits the adaptive update of the exponent function during the inversion procedure leading to good results even in this challenging scenario.
At first, the method is tested by considering a set of waveguide probes to illuminate a “vacuum-filled” target buried in a two-layered configuration with a planar interface. The behavior and limitations of the method such as the size and depth of the target change have been studied. Moreover, the effects of geometries and the number of targets as well as the uncertainties in soil dielectric permittivity have been investigated. In addition, the behavior of the method in the case of the non-planar air–soil interface has been analyzed, as well as the effect of the a priori knowledge of the interface profile inside the inversion procedure. The potentialities of this method for shallow subsurface prospection are shown by the achieved results. Future developments include both the validation with experimentally measured data and the extension to three-dimensional settings.

Author Contributions

Conceptualization, V.S., C.E., G.L.G., M.P., A.R. and A.F.; Methodology, C.E. and A.F.; Software, V.S.; Formal analysis, V.S. and C.E.; Data curation, A.R.; Writing—original draft, V.S. and A.F.; Writing—review & editing, G.L.G., A.R. and A.F.; Supervision, G.L.G., M.P., A.R. and A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The numerical data presented in this study are available from the corresponding authors on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Persico, R. Introduction to Ground Penetrating Radar: Inverse Scattering and Data Processing; John Wiley & Sons: Hoboken, NJ, USA, 2014; ISBN 978-1-118-30500-3. [Google Scholar]
  2. Benedetto, A.; Pajewski, L. Civil Engineering Applications of Ground Penetrating Radar; Springer: Cham, Switzerland, 2015; ISBN 978-3-319-04813-0. [Google Scholar]
  3. Frigui, H.; Zhang, L.; Gader, P.D. Context-Dependent Multisensor Fusion and Its Application to Land Mine Detection. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2528–2543. [Google Scholar] [CrossRef]
  4. Pastorino, M.; Randazzo, A. Microwave Imaging Methods and Applications; Artech House: London, UK, 2018; ISBN 978-1-63081-526-4. [Google Scholar]
  5. Salucci, M.; Oliveri, G.; Massa, A. GPR Prospecting through an Inverse-Scattering Frequency-Hopping Multifocusing Approach. IEEE Trans. Geosci. Remote Sens. 2015, 53, 6573–6592. [Google Scholar] [CrossRef]
  6. Catapano, I.; Randazzo, A.; Slob, E.; Solimene, R. GPR Imaging via Qualitative and Quantitative Approaches. In Civil Engineering Applications of Ground Penetrating Radar; Benedetto, A., Pajewski, L., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 239–280. ISBN 978-3-319-04812-3. [Google Scholar]
  7. Castaldo, R.; Crocco, L.; Fedi, M.; Garofalo, B.; Persico, R.; Rossi, A.; Soldovieri, F. GPR Microwave Tomography for Diagnostic Analysis of Archaeological Sites: The Case of a Highway Construction in Pontecagnano (Southern Italy). Archaeol. Prospect. 2009, 16, 203–217. [Google Scholar] [CrossRef]
  8. Nikolova, N.K. Introduction to Microwave Imaging; EuMA High Frequency Technologies Series; Cambridge University Press: Cambridge, UK, 2017; ISBN 978-1-107-08556-5. [Google Scholar]
  9. Özdemir, C.; Demirci, Ş.; Yiğit, E.; Yilmaz, B. A Review on Migration Methods in B-Scan Ground Penetrating Radar Imaging. Math. Probl. Eng. 2014, 2014, 280738. Available online: https://www.hindawi.com/journals/mpe/2014/280738/ (accessed on 4 January 2018). [CrossRef] [Green Version]
  10. Fallahpour, M.; Case, J.T.; Ghasr, M.T.; Zoughi, R. Piecewise and Wiener Filter-Based SAR Techniques for Monostatic Microwave Imaging of Layered Structures. IEEE Trans. Antennas Propag. 2014, 62, 282–294. [Google Scholar] [CrossRef]
  11. Firoozabadi, R.; Miller, E.L.; Rappaport, C.M.; Morgenthaler, A.W. Subsurface Sensing of Buried Objects Under a Randomly Rough Surface Using Scattered Electromagnetic Field Data. IEEE Trans. Geosci. Remote Sens. 2007, 45, 104–117. [Google Scholar] [CrossRef]
  12. Soldovieri, F.; Solimene, R. Ground Penetrating Radar Subsurface Imaging of Buried Objects. In Radar Technology; Kouemou, G., Ed.; InTech: London, UK, 2010; ISBN 978-953-307-029-2. [Google Scholar]
  13. Ludeno, G.; Gennarelli, G.; Lambot, S.; Soldovieri, F.; Catapano, I. A Comparison of Linear Inverse Scattering Models for Contactless GPR Imaging. IEEE Trans. Geosci. Remote Sens. 2020, 58, 7305–7316. [Google Scholar] [CrossRef]
  14. Prokopovich, I.; Popov, A.; Pajewski, L.; Marciniak, M. Application of Coupled-Wave Wentzel-Kramers-Brillouin Approximation to Ground Penetrating Radar. Remote Sens. 2018, 10, 22. [Google Scholar] [CrossRef] [Green Version]
  15. Bevacqua, M.T.; Isernia, T. Boundary Indicator for Aspect Limited Sensing of Hidden Dielectric Objects. IEEE Geosci. Remote Sens. Lett. 2018, 15, 838–842. [Google Scholar] [CrossRef]
  16. Ambrosanio, M.; Bevacqua, M.T.; Isernia, T.; Pascazio, V. Performance Analysis of Tomographic Methods against Experimental Contactless Multistatic Ground Penetrating Radar. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 1171–1183. [Google Scholar] [CrossRef]
  17. Gurbuz, T.U.; Aslanyurek, B.; Karabulut, E.P.; Akduman, I. An Efficient Nonlinear Imaging Approach for Dielectric Objects Buried under a Rough Surface. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3013–3022. [Google Scholar] [CrossRef]
  18. Mojabi, P.; LoVetri, J. A Prescaled Multiplicative Regularized Gauss-Newton Inversion. IEEE Trans. Antennas Propag. 2011, 59, 2954–2963. [Google Scholar] [CrossRef]
  19. Estatico, C.; Fedeli, A.; Pastorino, M.; Randazzo, A. Buried Object Detection by Means of a Lp Banach-Space Inversion Procedure. Radio Sci. 2015, 50, 41–51. [Google Scholar] [CrossRef]
  20. Catapano, I.; Crocco, L.; Isernia, T. A Simple Two-Dimensional Inversion Technique for Imaging Homogeneous Targets in Stratified Media. Radio Sci. 2004, 39, RS1012. [Google Scholar] [CrossRef]
  21. Di Donato, L.; Crocco, L. Model-Based Quantitative Cross-Borehole GPR Imaging via Virtual Experiments. IEEE Trans. Geosci. Remote Sens. 2015, 53, 4178–4185. [Google Scholar] [CrossRef]
  22. Tetik, E.; Akduman, I. 3D Imaging of Dielectric Objects Buried under a Rough Surface by Using CSI. Int. J. Antennas Propag. 2015, 2015, 179304. [Google Scholar] [CrossRef] [Green Version]
  23. Schenone, V.; Fedeli, A.; Estatico, C.; Pastorino, M.; Randazzo, A. Experimental Assessment of a Novel Hybrid Scheme for Quantitative GPR Imaging. IEEE Geosci. Remote Sens. Lett. 2022, 19, 3510005. [Google Scholar] [CrossRef]
  24. Beilina, L. Application of the Finite Element Method in a Quantitative Imaging Technique. J. Comput. Methods Sci. Eng. 2016, 16, 755–771. [Google Scholar] [CrossRef] [Green Version]
  25. Beilina, L.; Lindström, E. An Adaptive Finite Element/Finite Difference Domain Decomposition Method for Applications in Microwave Imaging. Electronics 2022, 11, 1359. [Google Scholar] [CrossRef]
  26. Bucci, O.M.; Isernia, T. Electromagnetic Inverse Scattering: Retrievable Information and Measurement Strategies. Radio Sci. 1997, 32, 2123–2137. [Google Scholar] [CrossRef]
  27. Bucci, O.M.; Crocco, L.; Isernia, T.; Pascazio, V. Subsurface Inverse Scattering Problems: Quantifying, Qualifying, and Achieving the Available Information. IEEE Trans. Geosci. Remote Sens. 2001, 39, 2527–2538. [Google Scholar] [CrossRef]
  28. De Coster, A.; Lambot, S. Full-Wave Removal of Internal Antenna Effects and Antenna–Medium Interactions for Improved Ground-Penetrating Radar Imaging. IEEE Trans. Geosci. Remote Sens. 2019, 57, 93–103. [Google Scholar] [CrossRef]
  29. Imperatore, P.; Iodice, A.; Pastorino, M.; Pinel, N. Modelling Scattering of Electromagnetic Waves in Layered Media: An Up-to-Date Perspective. Int. J. Antennas Propag. 2017, 2017, 7513239. [Google Scholar] [CrossRef] [Green Version]
  30. Das, N.K.; Pozar, D.M. A Generalized Spectral-Domain Green’s Function for Multilayer Dielectric Substrates with Application to Multilayer Transmission Lines. IEEE Trans. Microw. Theory Tech. 1987, 35, 326–335. [Google Scholar] [CrossRef]
  31. Xiong, J.L.; Chew, W.C. A Newly Developed Formulation Suitable for Matrix Manipulation of Layered Medium Green’s Functions. IEEE Trans. Antennas Propag. 2010, 58, 868–875. [Google Scholar] [CrossRef]
  32. Imperatore, P.; Iodice, A.; Riccio, D. Electromagnetic Wave Scattering From Layered Structures with an Arbitrary Number of Rough Interfaces. IEEE Trans. Geosci. Remote Sens. 2009, 47, 1056–1072. [Google Scholar] [CrossRef]
  33. Fedeli, A.; Pastorino, M.; Ponti, C.; Randazzo, A.; Schettini, G. Through-the-Wall Microwave Imaging: Forward and Inverse Scattering Modeling. Sensors 2020, 20, 2865. [Google Scholar] [CrossRef]
  34. Bourlier, C.; Pinel, N.; Kubické, G. Method of Moments for 2D Scattering Problems: Basic Concepts and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2013; ISBN 978-1-118-64868-1. [Google Scholar]
  35. Warren, C.; Giannopoulos, A.; Giannakis, I. GprMax: Open Source Software to Simulate Electromagnetic Wave Propagation for Ground Penetrating Radar. Comput. Phys. Commun. 2016, 209, 163–170. [Google Scholar] [CrossRef] [Green Version]
  36. Chew, W.C. Waves and Fields in Inhomogeneous Media; IEEE Press Series on Electromagnetic Wawes; IEEE Press: Piscataway, NJ, USA, 1995; ISBN 978-0-7803-4749-6. [Google Scholar]
  37. Hajebi, M.; Hoorfar, A. A Multiscale Optimization Technique for Large-Scale Subsurface Profiling. IEEE Geosci. Remote Sens. Lett. 2020, 18, 1706–1710. [Google Scholar] [CrossRef]
  38. Liu, Y.; Guo, L. A Time-Delay Calibration Method for Profile Estimation of Two-Layered Rough Surfaces. IEEE Access 2019, 7, 101575–101582. [Google Scholar] [CrossRef]
  39. Sefer, A.; Yapar, A. An Iterative Algorithm for Imaging of Rough Surfaces Separating Two Dielectric Media. IEEE Trans. Geosci. Remote Sens. 2021, 59, 1041–1051. [Google Scholar] [CrossRef]
  40. Song, T.; Kuang, L.; Han, L.; Wang, Y.; Liu, Q.H. Inversion of Rough Surface Parameters From SAR Images Using Simulation-Trained Convolutional Neural Networks. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1130–1134. [Google Scholar] [CrossRef]
  41. Cmielewski, O.; Tortel, H.; Litman, A.; Saillard, M. A Two-Step Procedure for Characterizing Obstacles Under a Rough Surface From Bistatic Measurements. IEEE Trans. Geosci. Remote Sens. 2007, 45, 2850–2858. [Google Scholar] [CrossRef]
  42. Pastorino, M. Stochastic Optimization Methods Applied to Microwave Imaging: A Review. IEEE Trans. Antennas Propag. 2007, 55, 538–548. [Google Scholar] [CrossRef]
  43. Wei, Z.; Chen, X. Deep-Learning Schemes for Full-Wave Nonlinear Inverse Scattering Problems. IEEE Trans. Geosci. Remote Sens. 2019, 57, 1849–1860. [Google Scholar] [CrossRef]
  44. Li, L.; Wang, L.G.; Teixeira, F.L.; Liu, C.; Nehorai, A.; Cui, T.J. DeepNIS: Deep Neural Network for Nonlinear Electromagnetic Inverse Scattering. IEEE Trans. Antennas Propag. 2019, 67, 1819–1825. [Google Scholar] [CrossRef] [Green Version]
  45. Bakhtiari, S.; Ganchev, S.I.; Zoughi, R. Open-Ended Rectangular Waveguide for Nondestructive Thickness Measurement and Variation Detection of Lossy Dielectric Slabs Backed by a Conducting Plate. IEEE Trans. Instrum. Meas. 1993, 42, 19–24. [Google Scholar] [CrossRef]
  46. Chen, X. Computational Methods for Electromagnetic Inverse Scattering, 1st ed.; John Wiley & Sons: Hoboken, NJ, USA; IEEE Press: Singapore, 2018; ISBN 978-1-119-31198-0. [Google Scholar]
  47. Taşkın, U.; Özdemir, Ö. Sparsity Regularized Nonlinear Inversion for Microwave Imaging. IEEE Geosci. Remote Sens. Lett. 2017, 14, 2220–2224. [Google Scholar] [CrossRef]
  48. de Zaeytijd, J.; Franchois, A.; Eyraud, C.; Geffrin, J.-M. Full-Wave Three-Dimensional Microwave Imaging with a Regularized Gauss-Newton Method—Theory and Experiment. IEEE Trans. Antennas Propag. 2007, 55, 3279–3292. [Google Scholar] [CrossRef]
  49. Estatico, C.; Pastorino, M.; Randazzo, A. A Novel Microwave Imaging Approach Based on Regularization in Lp Banach Spaces. IEEE Trans. Antennas Propag. 2012, 60, 3373–3381. [Google Scholar] [CrossRef]
  50. Abubakar, A.; Habashy, T.M.; Pan, G.; Li, M.-K. Application of the Multiplicative Regularized Gauss-Newton Algorithm for Three-Dimensional Microwave Imaging. IEEE Trans. Antennas Propag. 2012, 60, 2431–2441. [Google Scholar] [CrossRef]
  51. Estatico, C.; Fedeli, A.; Pastorino, M.; Randazzo, A.; Tavanti, E. A Phaseless Microwave Imaging Approach Based on a Lebesgue-Space Inversion Algorithm. IEEE Trans. Antennas Propag. 2020, 68, 8091–8103. [Google Scholar] [CrossRef]
  52. Bisio, I.; Estatico, C.; Fedeli, A.; Lavagetto, F.; Pastorino, M.; Randazzo, A.; Sciarrone, A. Variable-Exponent Lebesgue-Space Inversion for Brain Stroke Microwave Imaging. IEEE Trans. Microw. Theory Tech. 2020, 68, 1882–1895. [Google Scholar] [CrossRef]
  53. Estatico, C.; Fedeli, A.; Pastorino, M.; Randazzo, A. Microwave Imaging by Means of Lebesgue-Space Inversion: An Overview. Electronics 2019, 8, 945. [Google Scholar] [CrossRef] [Green Version]
  54. Fedeli, A.; Schenone, V.; Randazzo, A.; Pastorino, M.; Henriksson, T.; Semenov, S. Nonlinear S-Parameters Inversion for Stroke Imaging. IEEE Trans. Microw. Theory Tech. 2021, 69, 1760–1771. [Google Scholar] [CrossRef]
  55. Pastorino, M. Microwave Imaging; John Wiley & Sons: Hoboken, NJ, USA, 2010; ISBN 978-0-470-27800-0. [Google Scholar]
  56. Pelosi, G.; Coccioli, R.; Selleri, S. Quick Finite Elements for Electromagnetic Waves, 2nd ed.; Artech House: Norwood, MA, USA, 2009; ISBN 978-1-4020-9181-0. [Google Scholar]
  57. Berenger, J.-P. A Perfectly Matched Layer for the Absorption of Electromagnetic Waves. J. Comput. Phys. 1994, 114, 185–200. [Google Scholar] [CrossRef]
  58. Geuzaine, C.; Remacle, J.-F. Gmsh: A 3-D Finite Element Mesh Generator with Built-in Pre- and Post-Processing Facilities. Int. J. Numer. Methods Eng. 2009, 79, 1309–1331. [Google Scholar] [CrossRef]
  59. Catmull, E.; Rom, R. A Class of Local Interpolating Splines. In Computer Aided Geometric Design; Barnhill, R.E., Riesenfeld, R.F., Eds.; Academic Press: Cambridge, MA, USA, 1974; pp. 317–326. ISBN 978-0-12-079050-0. [Google Scholar]
Figure 1. Geometry of the shallow subsurface scenario.
Figure 1. Geometry of the shallow subsurface scenario.
Sensors 23 00167 g001
Figure 2. Measurement system.
Figure 2. Measurement system.
Sensors 23 00167 g002
Figure 3. Sketch of simulation and investigation domains.
Figure 3. Sketch of simulation and investigation domains.
Sensors 23 00167 g003
Figure 4. Flowchart of the inversion procedure.
Figure 4. Flowchart of the inversion procedure.
Sensors 23 00167 g004
Figure 5. Relative reconstruction errors in all of the investigation domain R and inside the target region versus the exponent range, defined by the parameter p m i n .
Figure 5. Relative reconstruction errors in all of the investigation domain R and inside the target region versus the exponent range, defined by the parameter p m i n .
Sensors 23 00167 g005
Figure 6. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Reference case.
Figure 6. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Reference case.
Sensors 23 00167 g006
Figure 7. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Target with side length (a) l t = 18 cm, (b) l t = 12 cm, and (c) l t = 3 cm.
Figure 7. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Target with side length (a) l t = 18 cm, (b) l t = 12 cm, and (c) l t = 3 cm.
Sensors 23 00167 g007
Figure 8. Reconstruction errors in the whole investigation domain and inside the target region versus target size, l t with the variable exponent approach and classical Hilbertian space.
Figure 8. Reconstruction errors in the whole investigation domain and inside the target region versus target size, l t with the variable exponent approach and classical Hilbertian space.
Sensors 23 00167 g008
Figure 9. Reconstructed distribution of the real part of the reconstructed relative dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Target centered at (a) y t = 20   cm , (b) y t = 25   cm , and (c) y t = 32.5   cm .
Figure 9. Reconstructed distribution of the real part of the reconstructed relative dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Target centered at (a) y t = 20   cm , (b) y t = 25   cm , and (c) y t = 32.5   cm .
Sensors 23 00167 g009
Figure 10. Simulated and reconstructed scattering S-parameters (real and imaginary part) when the first antenna acts as the source. Target centered at (a) y t = 20   cm and (b) y t = 25   cm .
Figure 10. Simulated and reconstructed scattering S-parameters (real and imaginary part) when the first antenna acts as the source. Target centered at (a) y t = 20   cm and (b) y t = 25   cm .
Sensors 23 00167 g010
Figure 11. Reconstruction errors in the whole domain and inside the target versus the target size, y t .
Figure 11. Reconstruction errors in the whole domain and inside the target versus the target size, y t .
Sensors 23 00167 g011
Figure 12. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. The case of two buried targets with different geometry.
Figure 12. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. The case of two buried targets with different geometry.
Sensors 23 00167 g012
Figure 13. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Background dielectric permittivity in the inversion procedure: (a) underestimation ( ε r , b * = ε r , b 0.5 ) and (b) overestimation ( ε r , b * = ε r , b + 0.5 ).
Figure 13. Reconstructed distribution of the real part of the dielectric permittivity, R e { ε r } , in the investigation domain R together with a square box indicating the shape of the actual target. Background dielectric permittivity in the inversion procedure: (a) underestimation ( ε r , b * = ε r , b 0.5 ) and (b) overestimation ( ε r , b * = ε r , b + 0.5 ).
Sensors 23 00167 g013
Figure 14. Measurement and subsurface model in the forward solver.
Figure 14. Measurement and subsurface model in the forward solver.
Sensors 23 00167 g014
Figure 15. Reconstruction errors in the whole domain and inside the target versus the RMS height of the air–soil interface, h r m s .
Figure 15. Reconstruction errors in the whole domain and inside the target versus the RMS height of the air–soil interface, h r m s .
Sensors 23 00167 g015
Figure 16. Reconstructed distribution of the real part of the complex dielectric permittivity, R e { ε r } , in the investigation domain together with a square box indicating the shape of the actual target. Cases of surface with h r m s = λ 30 and inversion with (a) exact interfaces and (b) planar interfaces; h r m s = λ 20 and inversion with (c) exact interfaces and (d) planar interfaces; and h r m s = λ 10 inversion with (e) exact interfaces and (f) planar interfaces.
Figure 16. Reconstructed distribution of the real part of the complex dielectric permittivity, R e { ε r } , in the investigation domain together with a square box indicating the shape of the actual target. Cases of surface with h r m s = λ 30 and inversion with (a) exact interfaces and (b) planar interfaces; h r m s = λ 20 and inversion with (c) exact interfaces and (d) planar interfaces; and h r m s = λ 10 inversion with (e) exact interfaces and (f) planar interfaces.
Sensors 23 00167 g016
Table 1. Reconstruction errors in the whole domain and inside the target versus estimation of background dielectric permittivity.
Table 1. Reconstruction errors in the whole domain and inside the target versus estimation of background dielectric permittivity.
Background Whole   Domain   Error ,   e R Target   Error ,   e t a r
Exact estimation0.0540.563
Underestimation0.1870.378
Overestimation0.2150.766
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Schenone, V.; Estatico, C.; Gragnani, G.L.; Pastorino, M.; Randazzo, A.; Fedeli, A. Microwave-Based Subsurface Characterization through a Combined Finite Element and Variable Exponent Spaces Technique. Sensors 2023, 23, 167. https://doi.org/10.3390/s23010167

AMA Style

Schenone V, Estatico C, Gragnani GL, Pastorino M, Randazzo A, Fedeli A. Microwave-Based Subsurface Characterization through a Combined Finite Element and Variable Exponent Spaces Technique. Sensors. 2023; 23(1):167. https://doi.org/10.3390/s23010167

Chicago/Turabian Style

Schenone, Valentina, Claudio Estatico, Gian Luigi Gragnani, Matteo Pastorino, Andrea Randazzo, and Alessandro Fedeli. 2023. "Microwave-Based Subsurface Characterization through a Combined Finite Element and Variable Exponent Spaces Technique" Sensors 23, no. 1: 167. https://doi.org/10.3390/s23010167

APA Style

Schenone, V., Estatico, C., Gragnani, G. L., Pastorino, M., Randazzo, A., & Fedeli, A. (2023). Microwave-Based Subsurface Characterization through a Combined Finite Element and Variable Exponent Spaces Technique. Sensors, 23(1), 167. https://doi.org/10.3390/s23010167

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop