Open access peer-reviewed chapter

Non-Invasive Microwave Characterization of Dielectric Scatterers

By Sandra Costanzo, Giuseppe Di Massa, Matteo Pastorino, Andrea Randazzo and Antonio Borgia

Submitted: November 26th 2011Reviewed: June 20th 2012Published: November 14th 2012

DOI: 10.5772/50842

Downloaded: 1267

1. Introduction

Microwave tomography is a technique aimed at inspecting unknown bodies by using an incident radiation generated at microwave frequencies (Ali &Moghaddam, 2010; Bellizzi, Bucci, &Catapano, 2011; Catapano, Crocco, &Isernia, 2007; Chen, 2008; Ferraye, Dauvignac, &Pichot, 2003; Gilmore, Mojabi, &LoVetri, 2009; Habashy&Abubakar, 2004; Isernia, Pascazio, &Pierri, 2001; Kharkovsky&Zoughi, 2007; Lesselier& Bowler, 2002; Litman, Lesselier, &Santosa, 1998; Oliveri, Lizzi, Pastorino, & Massa, 2012; Pastorino, 2010; Rekanos, 2008; Schilz&Schiek, 1981; Shea, Kosmas, Van Veen, &Hagness, 2010; Zhang & Liu, 2004; Zhou, Takenaka, Johnson, & Tanaka, 2009). An illuminating system(Costanzo & Di Massa, 2011; Paulsen, Poplack, Li, Fanning, &Meaney, 2000; Zoughi, 2000) is used to produce the incident waves that interact with the body to produce a scattered electromagnetic field. Another system is used to acquire the measurements that are used as input values for the reconstruction procedures. These values are the field samples resulting from the sum of the incident and scattered waves. Since the incident field (i.e., the field produced by the illuminating system when the object is not present) is a known quantity, the scattered electric field can be obtained by a direct subtraction. Moreover, the scattering field is related to the properties of the unknown body by well-known key relationships. In particular, both position, shape and dielectric parameters of the target affect the scattered field. In this Chapter, we consider the inspection of (possibly inhomogeneous) dielectric targets, which are characterized by the distributions of the dielectric permittivity and electric conductivity, whereas magnetic materials (e.g., materials for which the magnetic permittivity is different from the vacuum one) are not considered(El-Shenawee, Dorn, &Moscoso, 2009; Franchois&Pichot, 1997). The relationship between the target properties and the sampled scattered electric field is, in integral form, a Fredholm equation of the first kind, usually indicated as the “data equation”(Bucci, Cardace, Crocco, &Isernia, 2001; Rocca, Benedetti, Donelli, Franceschini, & Massa, 2009; van den Berg &Abubakar, 2001). The kernel of this equation is the Green’s function for two dimensional geometries in free space (in this Chapter we consider imaging configurations in free space, although other configurations - e.g., half space imaging - could be assumed by introducing the proper Green’s function for the specific configuration). The considered problem belongs to the class of inverse problems, which are usually ill-posed, in the sense that the solution can be not unique and unstable. To face the ill-posedness of the problem, the “data equation” is often solved together with the so-called “state equation”, relating the incident electric field inside the inspected object to the problem unknowns. In particular, in the developed approach, the two equations are combined together and a single nonlinear equation is obtained. In order to numerically solve the inverse problem, a discretization is usually necessary. A pixelated image of the scattering cross section can be obtained by using square pulse basis functions. The discrete nature of the measurements (we assume that each measurement antenna is able to collect the field at a given point inside a fixed observation domain) is equivalent to consider Dirac delta functions as testing function. The result of the discretization is a (nonlinear) system of equations to be solved, usually very ill-conditioned. In order to solve, in a regularized sense, the inverse scattering problem in the discrete setting, an iterative algorithm based on an inexact-Newton method is applied (Bozza, Estatico, Pastorino, &Randazzo, 2006; Estatico, Bozza, Massa, Pastorino, &Randazzo, 2005).

The reconstruction method proposed in this Chapter can be, in principle, applied to a large variety of dielectric objects, having homogeneous or multilayer cross-sections with arbitrary shape. Only for demonstration purpose, a simple homogeneous reference target of known dielectric properties is assumed in the following. The scattered field is acquired, both in amplitude and phase, on a square investigation domain around the target, sufficiently extent to be within the radiating near-field region (Costanzo & Di Massa, 2011). The incident field, oriented along the cylindrical target axis, is produced by a standard horn antenna, and a probe of the same kind is used to collect the field on the acquisition domain, for different positions of the illuminating horn. The measured scattered field data are subsequently processed to solve the inverse scattering problem and successfully retrieve the dielectric profile of the target under test.

The Chapter is organized as follows. In Section 2, a detailed mathematical description of the reconstruction method and the relative solving procedure are provided. The imaging setup configuration and the performed scattering measurements are discussed in Section 3. Some preliminary results concerning the inversion of measured data are reported in Section 4. Finally, conclusions are outlined in Section 5.

2. Mathematical formulation

The considered approach assumes tomographic imaging conditions and it aims at reconstructing the distributions of the dielectric properties of a slice of the target (Fig. 1). A transmitting (TX) antenna is successively positioned in Sdifferent locationsrsTX, s=1,,S, and generates a set of known z-polarized incident waves, whose electric field vectors can be expressed as:

Figure 1.

Schematic representation of the considered tomographic configuration

Eincsr,ω=eincsx,y,ωz^E1

where ωis the angular working frequency. It should be noted that the proposed approach does not require plane wave illumination of the target.

The object is assumed to have cylindrical geometry, with the cylindrical axis directed along a direction parallel to the electric field vectors of the incident electric field (i.e., the z axis). Moreover, the dielectric properties are assumed to be independent from the z coordinate, i.e.,ϵr=ϵx,yandσr=σx,y, being ϵand σthe dielectric permittivity and the electric conductivity, respectively.

The object interacts with the impinging electric field. From the above hypotheses it results that the resulting total electric field is z-polarized, too, and it can be written as:

Etotsr,ω=etotsx,y,ωz^=escattsx,y,ωz^+eincsx,y,ωz^E2

where escattsx,y,ωis the scattered electric field (due to the sth illumination), which is a mathematical quantity taking into account for the interaction effect between the incident electric field and the target. The total electric field is measured, for any location of the transmitting antenna, by a receiving (RX) antenna successively positioned inM pointsrs,mRX, s=1,,S,m=1,,M.

From a mathematical point of view, the relationship between the measured total electric field and the dielectric properties of the target can be modeled by using a Lippmann-Schwinger equation(Pastorino, 2010), i.e.,

etotsx,y,ω=eincsx,y,ω+jk024Dτx',y'etotsx',y',ωH02k0ρdx'dy',
s=1,,SE3

where k0=ωϵ0μ0is the free-space wavenumber (being ϵ0and μ0the dielectric permittivity and the magnetic permeability of the vacuum, respectively). In equation (3), H02is the Hankel function of zero-th order and second kind, ρ=x-x'2+y-y'2, and τdenotes the contrast function, which is defined as:

τx,y=ϵx,y-jσx,yωϵ0-1E4

The imaging procedure, starting from the measured samples of the total electric field in the S×Mmeasurement locations, is aimed at retrieving the contrast function, which contains all information about the unknown distributions of the dielectric properties. Consequently, equation (3) is discretized (by using pulse basis function to represent the unknowns(Richmond, 1965)) and computed at the measurement positions, leading to the following set of discrete equations:

escatts=Hsdiagτetots,
s=1,,SE5

where:

escatts=escattsxs,1RX,ys,1RX,ωescattsxs,MRX,ys,MRX,ω,
s=1,,SE6

is an array containing the samples of the scattered electric field at the Mmeasurement points (beingxs,mRX,ys,mRX, m=1,,M, their positions in the transverse plane) for the sth illumination, the vector:

etots=etotsx1D,y1D,ωetotsxND,yND,ωE7

is an array containing the values of total electric field in the centersxiD,yiD, i=1,,N, of the Nsubdomains used to discretize the investigation area, and:

τ=τx1D,y1DτxND,yNDE8

is an array containing the coefficients of the discretized contrast function. Finally, the matrix Hsis given by:

Hs=jk024hs,1,1hs,1,Nhs,M,1hs,M,N,
s=1,,SE9

whose elements are provided by the following relation:

hs,m,i=DiH02k0ρs,mdx'dy'E10

being Dithe ith subdomain andρs,m=xs,mRX-x'2+ys,mRX-y'2.

It is worth noting that the array etotsis unknown into equation (5). Consequently, a second equation is needed to solve the inverse problem. Such relation is found by applying (3) to points inside the investigation area. By using the same discretization, we obtain the following matrix relation:

etots=eincs-Gdiagτetots,
s=1,,SE11

where:

eincs=eincsx1D,y1D,ωeincsxND,yND,ω,
s=1,,SE12

is an array containing the values of the incident electric field in the centers of the NsubdomainsDi. Moreover, the term:

G=jk024g1,1g1,NgN,1gN,NE13

is a matrix whose elements are given by:

gi,k=DkH02k0ρidx'dy'E14
with
ρi=xiD-x'2+yiD-y'2E15
.

Equations (5) and (11) are combined together in order to obtain the following set of nonlinear equations:

escatts=HdiagτI-Gsdiagτ-1eincs=Asτ,
s=1,,SE16

which can be written as:

escatt=escatt1escattS=A1τASτ=AτE17

Equation (16) needs to be solved in order to retrieve the contrast function τ.Once this term is obtained, the distributions of the dielectric parameters can be calculated using equation (4). The solution of equation (16) represents however a highly ill-posed problem. Consequently, a regularized inversion algorithm must be used (Autieri, Ferraiuolo, &Pascazio, 2011; Lobel, Blanc-Féraud, Pichot, &Barlaud, 1997). The solving procedure is based on a two-step iterative strategy (Bozza et al., 2006; Estatico et al., 2005), in which an outer linearization is performed by means of an Inexact-Newton scheme and a regularized solution to the obtained linear system is calculated by means of a truncated Landweber algorithm (Landweber, 1951).

The developed iterative procedure works as follows.

  1. Set iteration index to n=0and initialize the unknownτn, e.g., by choosingτn=0;

  2. Linearize the equation escatt=Aτin order to obtain a linear equationJnh=en, where escatt-Aτnand Jnis the Jacobian matrix (i.e., the discrete counterpart of the Frechét derivative) of Aat pointτn, which is given by(Remis& van den Berg, 2000):

Jn=Jn1JnS=Hn1diagetotn1HnSdiagetotnSE18

being etotns=I-Gdiagτn-1eincsthe total electric field (for the sth illumination) inside the investigation area due tothe current estimate ofthe solution τnand Hnsa inhomogeneous Green matrix given byHns=I-Gdiagτn-1Hs;

  1. Find a regularized solution hof the linearized equation by using a truncated Landweber algorithm;

  2. Update current solution withτn+1=τn+h;

  3. Check if a convergence criteria (e.g., a maximum number of iterations nmaxor a threshold on the residual) is fulfilled. Otherwise go back to step 2.

The truncated Landweber algorithm can be summarized as follows.

  1. Set iteration index to l=0and initialize the unknownhl=0;

  2. Update the current solution with

hl+1=hl-βJn*Jnhl-enE19

where β=0.5/Jn*Jn2and Jn*is the adjoint ofJn;

  1. Check if a convergence criteria (e.g., a maximum number of iterations, lmax, or a threshold on the residual) is fulfilled. Otherwise go back to step 2.

3. Imaging setup configuration and measurements

In order to prove the effectiveness of the approach, experimental validations are performed on a reference target by adopting the setup configuration in Fig. 2.A square investigation domain Daround the scatterer is assumed, with a transmitting antenna giving an incident field oriented along the longitudinal axis of the target and assuming the positions indicated in Fig. 2as black circles. A receiving antenna is adopted to collect the complex (amplitude and phase) field scattered by the test object on the void circles, at spacingsΔx, Δysatisfying the Shannon’s sampling theorem(Bucci&Franceschetti, 1989).

3.1. X-band measurements

The imaging configuration described in Fig. 2is realized into the anechoic chamber of Microwave Laboratory at University of Calabria. Two standard X-band horn antennas are adopted as transmitting and receiving antennas, as providing a sufficiently large pattern impinging on the reference target. This is given by a wood cylinder of length equal about to 60 cm (~ 20 λ @ 10 GHz) and square cross section of side equal about to 4 cm. A photograph showing the imaging setup is reported inFig. 3.

Scattered field measurements are performed at a frequency equal to 10 GHz, on a square measurement domain of side equal to39 cm, which has been discretized into4Mpoints (M=53) with spacingsΔx=Δy=λ/4. As highlighted above, both transmitting and receiving horns are oriented with the field parallel to the cylinder axis. The amplitude and phase behavior of the measured scattered field is reported in Figs. 4-7 for two different positions of the illuminating horn, along two different sides of the acquisition domain scanned by the receiving probe. The positions of transmitting and receiving antennas are visible in the picture within the same figures.Similar behaviors are obtained for the other positions of the transmitting and the receiving antennas.

Figure 2.

Imaging setup configuration

Figure 3.

Photograph of measurement setup into the Microwave Laboratory at University of Calabria

Figure 4.

Amplitude (a) and phase (b) of measured scattered field: configuration in the picture

Figure 5.

Amplitude (a) and phase (b) of measured scattered field: configuration in the picture

Figure 6.

Amplitude (a) and phase (b) of measured scattered field: configuration in the picture

Figure 7.

Amplitude (a) and phase (b) of measured scattered field: configuration in the picture

4. Preliminary reconstruction results

In this section, preliminary reconstruction results are reported. Figure 8 provides the reconstructed image of the object described in Section 3. It is obtained by inverting the real measured data described in Section 3 by using the procedure described in Section 2.

Figure 8.

Reconstructed distribution of the relative dielectric permittivity inside the square investigation domain

For every side of the measurement domain, three positions of the TX antenna are used. In particular, for the first side (i.e., the one characterized by coordinatex=-19.5 cm), the y-positions of TX antenna are equal to-3.625, 0, and3.625 cm. For these three source positions, only the Mmeasurement points located on the opposite side of the measurement domain (i.e., for the first side, those characterized by coordinatex=19.5 cm) are used. The remaining views are constructed in a similar way. The total number of illuminations is S=12and the total number of measured samples isS×M=636. The investigation area is assumed of square shape and side equal to 0.06 m. It is partitioned into 20×20 square subdomains. The algorithm is initialized by using a rough estimate obtained by means of a back-propagation algorithm (Lobel, Kleinman, Pichot, Blanc-Feraud, &Barlaud, 1996). The inner loop (Landweber algorithm) is stopped after a fixed number of iterationslmax=5. The outer loop (Newton linearization) is stopped according to the L-curve criteria (Vogel, 2002), leading to an estimated optimal number of outer iterations ofnlc=3. Maximum, minimum and mean values of the retrieved dielectric permittivity distribution are reported in Table 1. Finally, Fig. 9 reports the profiles obtained by cutting the 2D distribution along two horizontal and vertical axes passing from the center of the investigation domain. As can be seen, the presence of the target can be suitably retrieved.

Figure 9.

Horizontal and vertical profiles of the reconstructed distribution of the relative dielectric permittivity along lines passing from the center of the investigation domain

MinMaxMeanVariance
Whole domain1.002.901.670.48
Object1.102.902.300.20
Background1.001.411.080.01

Table 1.

Values of the retrieved relative dielectric permittivity

5. Conclusion

The non-invasive inspection of dielectric objects has been considered in this Chapter to provide an accurate characterization of the permittivity profile at microwave frequencies. A mathematical formulation in terms of a Fredholm integral equation of the first order has been assumed, and a suitable discretization has been performed in order to numerically solve the resulting inverse problem, with a regularization approach adopted to overcome the intrinsic ill-posedness. The proposed imaging technique has been experimentally assessed by performing scattered field measurements on a square investigation domain surrounding a cylindrical dielectric target of known properties. Measured X-band data acquired by a standard horn antenna have been collected for different positions of the illuminating horn, and a successful reconstruction of the expected dielectric profile has been obtained from the application of the proposed technique.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sandra Costanzo, Giuseppe Di Massa, Matteo Pastorino, Andrea Randazzo and Antonio Borgia (November 14th 2012). Non-Invasive Microwave Characterization of Dielectric Scatterers, Microwave Materials Characterization, Sandra Costanzo, IntechOpen, DOI: 10.5772/50842. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="http://www.intechopen.com/embed/microwave-materials-characterization/non-invasive-microwave-characterization-of-dielectric-scatterers" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

1267total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Nondestructive Evaluations by Using a Prototype of a Microwave Tomograph

By R. Monleone, M. Pastorino, S. Poretti, A. Randazzo, A. Massimini and A. Salvade

Related Book

First chapter

Localized Bessel Beams: Basic Properties and Emerging Communication Applications

By Sandra Costanzo

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us