International Journal of Antennas and Propagation
Volume 2016 (2016), Article ID 9304371, 8 pages
http://dx.doi.org/10.1155/2016/9304371

Research Article

A Banach Space Regularization Approach for Multifrequency Microwave Imaging

Claudio Estatico,¹ Alessandro Fedeli,² Matteo Pastorino,² and Andrea Randazzo²

¹Department of Mathematics, University of Genoa, Via Dodecaneso 35, 16146 Genova, Italy
²Department of Electrical, Electronic, Telecommunications Engineering, and Naval Architecture, University of Genoa, Via Opera Pia 11A, 16145 Genova, Italy

Received 15 October 2015; Accepted 24 February 2016

Academic Editor: Ahmed T. Mobashsher

Copyright © 2016 Claudio Estatico et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A method for microwave imaging of dielectric targets is proposed. It is based on a tomographic approach in which the field scattered by an unknown target (and collected in a proper observation domain) is inverted by using an inexact-Newton method developed in Banach spaces. In particular, the extension of the approach to multifrequency data processing is reported. The mathematical formulation of the new method is described and the results of numerical simulations are reported and discussed, analyzing the behavior of the multifrequency processing technique combined with the Banach spaces reconstruction method.

1. Introduction

The regularization of ill-posed problems in Banach spaces exhibits several advantages over the corresponding classical regularization in Hilbert spaces [1–3]. A significant example is microwave imaging (MI). MI is a well-known technique in which an unknown target is inspected by using an electromagnetic radiation [4]. When illuminated by an incident field, the target scatters the incident radiation depending on its physical properties, geometrical shape, and dimensions. For dielectric targets with dimensions comparable with the wavelength of the incident waves, the scattering phenomena are quite complicated [5]. However, the scattered electric field can be collected around the target and the properties of the object can be retrieved by inverting the equation governing the electromagnetic scattering problem. Several approaches to solve this inverse problem with qualitative [6–8] or quantitative techniques [9–11] have been reported in the scientific literature, with reference to both two- and three-dimensional configurations [12–14]. Among these methods, it has been proved in [15] that, exploiting the properties of the norms in the Banach spaces, it is possible to obtain better reconstructions than developing the corresponding approach in Hilbert spaces. In particular, in [16] an inexact-Newton method has been applied to a tomographic configuration, in which the target is assumed to be a cylindrical one and the imaging process is performed in free space and at a fixed frequency. In this paper, the above formulation is extended to multifrequency imaging, in which the use of more than one incident radiation results in an improved processing, due to the additional information included in the input data (samples of the scattered electric field at different frequencies).

The mathematical formulation of the proposed approach is described in the following with reference to transverse-magnetic illumination conditions. The capabilities and limitations of the proposed extended method are evaluated by means of numerical simulations concerning targets in noisy environments for which the forward problem is solved by using the method of moments. The inverse solution is obtained, as mentioned, by using a two-loop inexact-Newton method, which has been developed in the framework of Banach spaces. In particular, after the discretization of the continuous model (based on the electric field integral equation (EFIE)), the outer loop performs a linearization of the resulting nonlinear algebraic equation, whereas the inner loop solves the obtained (ill-posed) linear equation in a regularized sense by using a truncated Landweber scheme in Banach spaces.

The paper is organized as follows. The mathematical formulation of the developed approach is discussed in Section 2. Section 3 reports some numerical results aimed at validating the inversion procedure. Finally, conclusions are drawn in Section 4.

2. Mathematical Formulation

Let us consider the configuration shown in Figure 1. A cylindrical target, whose cross section is enclosed in an investigation domain , is illuminated by a set of known time-harmonic, TM- electromagnetic fields characterized by angular frequencies . As it is well known, under such hypotheses, the scattering problem turns out to be two-dimensional and scalar [4]. The resulting scattered electric field is collected in a measurement domain .

Figure 1: Configuration of the electromagnetic inverse scattering problem and representation of the measurement setup.

For the sake of simplicity, in the following a single view case is described. Moreover, the dependence is omitted, assuming a frequency domain formulation. The scattered electric field in the measurement domain, ,  , at the angular frequency is related to the dielectric properties of the inspected area by the relationship [4]

where

is a vector containing the relative dielectric permittivity and the electric conductivity ,   ( is the free space dielectric permittivity and the superscript denotes the transposition operator), ,  , is the total electric field inside the investigation domain, and

being the free space Green’s function at the angular frequency [17, 18]. It is worth noting that, assuming nondispersive materials, results to be frequency independent and real-valued. The total electric field inside the investigation area is given by

where is the known incident electric field and is defined as

Combining (1) and (4), the scattering problem at the frequency can be written as

where is defined as

By applying (6) for all the available frequencies and by separating real and imaginary parts of the involved complex functions, the scattering problem can be written as the solution of the following nonlinear system of equations:

where ,  , and .

The inverse problem that must be solved in order to retrieve the dielectric properties of the targets consists in finding the unknown from the nonlinear equation (8), with being the known term. To this end, a Newton-type method [19–21] developed in the framework of Banach spaces [15, 16] is used. The proposed approach performs a minimization of the functional given by

with being the norm of the Banach space . In particular, the method contains two nested loops: the outer one performs a linearization of the scattering operator around the currently reconstructed function , whereas the inner loop finds a suitable regularized solution of the resulting linearized problem by using a truncated Landweber algorithm developed in Banach spaces [15].

The reconstruction method is thus composed by the following steps (Newton iterations are denoted by the index and Landweber iterations by the index ):

(a)Start the algorithm with an initial guess . If no a priori information is available, is assumed.(b) Compute the Newton linearization of the operator around the current value by means of its Fréchet derivative ; that is,

(c) Find a function by using a truncated Landweber regularization loop in Banach spaces, initialized with (the symbol indicates that the function belongs to the dual space of ):

 where is an adequate step length, is the adjoint operator of , and the operators and are the duality maps of the spaces and . We recall that, in Banach spaces, the duality maps are defined as

 with

 The Landweber steps (11) are iterated for .(d) Update the solution with the increment ; that is,

(e) Iterate Newton steps ‎(b)–‎(d) for , or until an a priori stopping criterion is reached.

An important remark has to be devoted to the Fréchet derivative of the operator , necessary for performing steps ‎(b) and ‎(c). In this case, it has the following structure:

where

In (16), the operator is defined as

with being the electric field due to the dielectric profile , and

where is the inhomogeneous Green’s function related to the dielectric profile [4].

3. Numerical Results

The proposed multifrequency imaging method in Banach spaces has been validated by using synthetic data. In the simulations presented in this section, we assume a multiview measurement setup, in which the object is sequentially illuminated by transmitting antennas located in different positions. With reference to the configuration reported in Figure 1, in the th view () the incident electric field is generated by an antenna, modeled as a line-current source, located at

with and . For each view, the field scattered by the target is collected in measurement points on the same circumference at positions

where (with ).

The scattered field data are computed by a numerical code based on the method of moments [22] and then corrupted by a Gaussian noise with zero mean value and variance corresponding to a fixed signal-to-noise ratio (). The considered investigation domain is a square area of side , centered at the origin. It has been discretized into square subdomains for solving the forward electromagnetic problem and cells for the inverse one. The object under test is composed of a known lossless structure (a circular cylinder of radius , centered at the origin, with relative dielectric permittivity ) with three cylindrical inclusions of different shapes and dimensions. The first one is a cylinder with rectangular cross section, characterized by , centered at , with side and side of length and , respectively. In addition, we have two circular cylinders of radii and , centered at and , with relative dielectric permittivities and , respectively. Such values have been chosen since many materials commonly used in industrial processes (e.g., plastic and ceramic compounds, wood, and concrete) have dielectric properties in this range.

The object under test has been illuminated by means of electromagnetic fields at three different frequencies: ,  , and  . The behavior of the inversion method has been evaluated with different numbers of frequencies used by the multifrequency reconstruction procedure (). For instance, with only the data at frequency is considered; for we use the data obtained at and , and for all the available frequency data are used. As regards the inversion method, inner Landweber iterations and outer Gauss-Newton steps have been executed in all cases. The algorithm has been initialized with the outer (known) dielectric cylinder of radius and relative dielectric permittivity as starting guess.

In order to analyze the method’s performance in different operating conditions, numerical simulations with and have been performed. Firstly, let us consider the case with and . The results obtained by the present Banach space inversion technique are shown in Figures 2(a)–2(c).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: Reconstructed distributions of the relative dielectric permittivity of the investigation domain in different operating conditions, with and . Banach space reconstructions: (a) ,  ; (b) ,  ; (c) ,  . Hilbert space reconstructions: (d) ,  ; (e) ,  ; (f) ,  .

For comparison purposes, the results provided by a standard inexact-Newton method working in Hilbert spaces [20] are also provided in Figures 2(d)–2(f). From (a) to (c) and from (d) to (f), the dielectric reconstructions obtained with one, two, and three frequencies are shown. Clearly, the improvement introduced by the Banach space approach is evident, even with a single operating frequency. An increment in the number of considered frequencies produces a further enhancement in the reconstruction quality, with both the Hilbert space and the Banach space techniques. However, as can be seen in Figure 2, the Banach space method outperforms the standard Hilbert space one in all cases. Furthermore, from a practical point of view, the Banach space results obtained with one and three frequencies are comparable. Therefore, we can say that in this case the proposed technique is able to provide very good results even with a reduced amount of data.

The graphs in Figure 3 are about a high-noise case, with . Even in this situation, the benefits of combining the Banach space regularization with the multifrequency processing are clear. In fact, the Hilbert space () reconstruction with frequencies presented in Figure 3(b) suffers from significant ringing effects on the background and on the object area. By lowering the Banach space norm parameter to (Figure 3(c)) it is possible to reduce these unwanted effects, improving the reconstruction performance.

(a)

(b)

(c)

Figure 3: Reconstructed distributions of the relative dielectric permittivity of the investigation domain in different operating conditions, with and : (a) ,  ; (b) , ;  (c) , .

In Figure 4 the inversion results with a low number of views () are presented. Under these operating conditions, due to the small amount of data, the Banach space method combined with the multifrequency approach seems to be useful for reducing the reconstruction artifacts, as can be seen from the comparison of Figure 4(a) (the standard single frequency Hilbert space reconstruction) with the results in Figures 4(b) and 4(c).

(a)

(b)

(c)

Figure 4: Reconstructed distributions of the relative dielectric permittivity of the investigation domain in different operating conditions, with and : (a) ,  ; (b) ,  ; (c) , .

The inversion results have been quantitatively evaluated using the following relative reconstruction error on the whole investigation domain

where denotes the actual value of the relative dielectric permittivity in the th cell () of the investigation domain and is the corresponding reconstructed quantity. The trends of the reconstruction error versus the and the norm parameter are reported in Figure 5. As can be noticed, in the multifrequency results of Figures 5(b) and 5(c) the relative reconstruction error is generally lower than the single frequency case of Figure 5(a). Moreover, the effectiveness of the Banach space reconstruction approach with in the multifrequency case appears to be more evident in presence of low values of . In Figure 6 the behavior of the reconstruction error versus the number of views and the parameter is shown.

(a)

(b)

(c)

Figure 5: Relative reconstruction error on the investigation domain versus the Banach space norm parameter and the , for : (a) one frequency, ; (b) two frequencies, ; (c) three frequencies, .

(a)

(b)

(c)

Figure 6: Relative reconstruction error on the investigation domain versus the Banach space norm parameter and the number of views , for : (a) one frequency, ; (b) two frequencies, ; (c) three frequencies, .

In this case, too, we can observe that on one hand the reconstruction error decreases with an increase in the number of considered frequencies; on the other hand, the use of a Banach space norm leads to significant improvements especially with a reduced number of views (i.e., with a low amount of data).

4. Conclusions

In this paper, a previously developed approach to microwave imaging has been extended to deal with multifrequency processing in free space. The proposed technique is based on an inexact-Newton method and is aimed at inspecting inhomogeneous dielectric cylinders of arbitrary shapes under transverse-magnetic illumination conditions. The use of the inexact-Newton method (modified to take profit of the features of Banach space norms) has been found to represent a very effective regularization tool for solving the electric field integral equations considered in this paper. Numerical results have confirmed the above reconstruction capabilities of the two-loop inversion method with reference to a tomographic imaging configuration. Moreover, the advantages in terms of the reconstruction accuracy of using multifrequency imaging have been proved at least for the considered targets. Future works will be devoted to a more comprehensive numerical assessment of the method, for example, by considering different measurement setups or targets with higher dielectric contrasts. Moreover, experimental data will also be used to further validate the proposed technique in more realistic conditions.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

T. Schuster, Ed., Regularization Methods in Banach Spaces, De Gruyter, Berlin, Germany, 2012.

T. Hein and K. S. Kazimierski, “Accelerated Landweber iteration in Banach spaces,” Inverse Problems, vol. 26, no. 5, Article ID 055002, 2010.

F. Schöpfer, A. K. Louis, and T. Schuster, “Nonlinear iterative methods for linear ill-posed problems in Banach spaces,” Inverse Problems, vol. 22, no. 1, pp. 311–329, 2006.

M. Pastorino, Microwave Imaging, John Wiley & Sons, Hoboken, NJ, USA, 2010.

D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, vol. 93, Springer, New York, NY, USA, 2013.

M. R. Eskandari, R. Safian, and M. Dehmollaian, “Three-dimensional near-field microwave imaging using hybrid linear sampling and level set methods in a medium with compact support,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 10, pp. 5117–5125, 2014.

I. Catapano, A. Randazzo, E. Slob, and R. Solimene, “GPR imaging via qualitative and quantitative approaches,” in Civil Engineering Applications of Ground Penetrating Radar, A. Benedetto and L. Pajewski, Eds., Springer Transactions in Civil and Environmental Engineering, pp. 239–280, Springer, Berlin, Germany, 2015.

M. Fallahpour, J. T. Case, M. T. Ghasr, and R. Zoughi, “Piecewise and Wiener filter-based SAR techniques for monostatic microwave imaging of layered structures,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 1, pp. 282–294, 2014.

M. Pastorino, M. Raffetto, and A. Randazzo, “Electromagnetic inverse scattering of axially moving cylindrical targets,” IEEE Transactions on Geoscience and Remote Sensing, vol. 53, no. 3, pp. 1452–1462, 2015.

W. Zhang and Q. H. Liu, “Three-dimensional scattering and inverse scattering from objects with simultaneous permittivity and permeability contrasts,” IEEE Transactions on Geoscience and Remote Sensing, vol. 53, no. 1, pp. 429–439, 2015.

L. Poli, G. Oliveri, F. Viani, and A. Massa, “MT-BCS-based microwave imaging approach through minimum-norm current expansion,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 9, pp. 4722–4732, 2013.

A. Massa, M. Pastorino, and A. Randazzo, “Reconstruction of two-dimensional buried objects by a differential evolution method,” Inverse Problems, vol. 20, no. 6, pp. S135–S150, 2004.

M. Brignone, G. Bozza, A. Randazzo, M. Piana, and M. Pastorino, “A hybrid approach to 3D microwave imaging by using linear sampling and ACO,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 10, pp. 3224–3232, 2008.

A. Abubakar, T. M. Habashy, G. Pan, and M.-K. Li, “Application of the multiplicative regularized Gauss-Newton algorithm for three-dimensional microwave imaging,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 5, pp. 2431–2441, 2012.

C. Estatico, M. Pastorino, and A. Randazzo, “A novel microwave imaging approach based on regularization in L p Banach spaces,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 7, pp. 3373–3381, 2012.

C. Estatico, A. Fedeli, M. Pastorino, and A. Randazzo, “Microwave imaging of elliptically shaped dielectric cylinders by means of an Lp Banach-space inversion algorithm,” Measurement Science and Technology, vol. 24, no. 7, Article ID 074017, 2013.

C. A. Balanis, Advanced Engineering Electromagnetics, John Wiley & Sons, Hoboken, NJ, USA, 2nd edition, 2012.

M. Pastorino, M. Raffetto, and A. Randazzo, “Two-dimensional Green's function for scattering and radiation problems in elliptically-layered media,” IEEE Transactions on Antennas and Propagation, vol. 62, no. 4, pp. 2071–2080, 2014.

M. Ostadrahimi, P. Mojabi, A. Zakaria, J. LoVetri, and L. Shafai, “Enhancement of Gauss-Newton inversion method for biological tissue imaging,” IEEE Transactions on Microwave Theory and Techniques, vol. 61, no. 9, pp. 3424–3434, 2013.

C. Estatico, G. Bozza, A. Massa, M. Pastorino, and A. Randazzo, “A two-step iterative inexact-Newton method for electromagnetic imaging of dielectric structures from real data,” Inverse Problems, vol. 21, no. 6, pp. S81–S94, 2005.

A. Randazzo, G. Oliveri, A. Massa, and M. Pastorino, “Electromagnetic inversion with the multiscaling inexact Newton method-experimental validation,” Microwave and Optical Technology Letters, vol. 53, no. 12, pp. 2834–2838, 2011.

J. Richmond, “Scattering by a dielectric cylinder of arbitrary cross section shape,” IEEE Transactions on Antennas and Propagation, vol. 13, no. 3, pp. 334–341, 1965.