Open access

Non-Invasive Microwave Characterization of Dielectric Scatterers

Written By

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

Submitted: 07 June 2012 Published: 14 November 2012

DOI: 10.5772/50842

From the Edited Volume

Microwave Materials Characterization

Edited by Sandra Costanzo

Chapter metrics overview

2,127 Chapter Downloads

View Full Metrics

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.

Advertisement

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,y andσ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μ0 is the free-space wavenumber (being ϵ0 and μ0 the 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 Hs is 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 Di the ith subdomain andρs,m=xs,mRX-x'2+ys,mRX-y'2.

It is worth noting that the array etots is 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τn and Jn is 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-1eincs the total electric field (for the sth illumination) inside the investigation area due tothe current estimate ofthe solution τn and Hns a 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 nmax or 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.

Advertisement

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 into4M points (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

Advertisement

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

Advertisement

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.

References

  1. 1. AliM. AMoghaddamM20103D Nonlinear Super-Resolution Microwave Inversion Technique Using Time-Domain DataIEEE Trans. Antennas Propag., 58page numbers (2327 EOF2336 EOF
  2. 2. AutieriRFerraiuoloGPascazioV2011Bayesian Regularization in Nonlinear Imaging: Reconstructions From Experimental Data in Nonlinearized Microwave TomographyIEEE Trans. Geosci. Remote Sens., 49page numbers (801 EOF813 EOF
  3. 3. BellizziGBucciO. MCatapanoI2011Microwave Cancer Imaging Exploiting Magnetic Nanoparticles as Contrast AgentIEEE Trans. Biomed. Eng., 58page numbers (2528 EOF2536 EOF
  4. 4. BozzaGEstaticoCPastorinoMRandazzoA2006An Inexact Newton Method for Microwave Reconstruction of Strong ScatterersIEEE Antennas Wirel. Propag. Lett., 5page numbers (61 EOF64 EOF
  5. 5. BucciCardace, N., Crocco, L., & Isernia, T. (2001Degree of Nonlinearity and A New Solution Procedure in Scalar Two-Dimensional Inverse Scattering Problems.J. Opt. Soc. Am. A, 18page numbers (1832 EOF43 EOF
  6. 6. BucciFranceschetti, G. (1989On the Degrees of Freedom of Scattered FieldsIEEE Trans. Antennas Propag., 37page numbers (918 EOF926 EOF
  7. 7. CatapanoICroccoLIserniaT2007On Simple Methods for Shape Reconstruction of Unknown ScatterersIEEE Trans. Antennas Propag., 55page numbers (1431 EOF1436 EOF
  8. 8. ChenX2008Signal-Subspace Method Approach to the Intensity-Only Electromagnetic Inverse Scattering ProblemJ. Opt. Soc. Am. A, 25page numbers (2018 EOF2024 EOF
  9. 9. CostanzoSDi Massa, G. (2011Advanced numerical techniques for near-field antenna measurements, In: Numerical simulations of physical and engineering processes, J. Awrejcewicz Ed., 321338INTECH, 978-9-53307-620-1Croatia
  10. 10. El-ShenaweeMDornOMoscosoM2009An Adjoint-Field Technique for Shape Reconstruction of 3-D Penetrable Object Immersed in Lossy MediumIEEE Trans. Antennas Propag., 57page numbers (520 EOF534 EOF
  11. 11. EstaticoCBozzaGMassaAPastorinoMRandazzoA2005ATwo-stepIterative Inexact-Newton Method for Electromagnetic Imaging of Dielectric Structures from Real Data. Inv. Probl., 21page numbers (S81 EOFS94 EOF
  12. 12. FerrayeRDauvignacJYPichotC. (2003An Inverse Scattering Method Based on Contour Deformations by Means of a Level Set Method Using Frequency Hopping TechniqueIEEE Trans. Antennas Propag., 51page numbers (1100 EOF
  13. 13. FranchoisAPichotC1997Microwave imaging-complex permittivity reconstruction with a Levenberg-Marquardt method. IEEE Trans. Antennas Propag., 45page numbers (203 EOF
  14. 14. GilmoreCMojabiPLoVetri, J. (2009Comparison of an Enhanced Distorted Born Iterative Method and the Multiplicative-Regularized Contrast Source Inversion methodIEEE Trans. Antennas Propag., 57page numbers (2341 EOF2351 EOF
  15. 15. HabashyT. MAbubakarA2004A General Framework for Constraint Minimization for the Inversion of Electromagnetic MeasurementsProg. Electromag. Res., 46page numbers (265 EOF312 EOF
  16. 16. IserniaTPascazioVPierriR2001On the local minima in a tomographic imaging techniqueIEEE Trans. Geosci. Remote Sens., 39page numbers (1596 EOF1607 EOF
  17. 17. KharkovskySZoughiR2007Microwave and millimeter wave nondestructive testing and evaluation- Overview and recent advances. IEEE Instrum. Meas. Mag., 10page numbers (26-38).
  18. 18. Landweber, L. (1951). An Iteration Formula for Fredholm Integral Equations of the First Kind. Am. J. Math., 73 (1951) page numbers (615-624).
  19. 19. LesselierDBowlerJ2002Special section on electromagnetic and ultrasonic nondestructive evaluationInv. Probl., 18
  20. 20. LitmanALesselierDSantosaF1998Reconstruction of a two-dimensional binary obstacle by controlled evolution of a level-setInv. Probl., 14page numbers (685 EOF
  21. 21. LobelPBlanc-féraudLPichotCBarlaudM1997A New Regularization Scheme for Inverse ScatteringInverse ProblemsInv. Probl., 13page numbers (403 EOF
  22. 22. LobelPKleinmanR. EPichotCBlanc-feraudLBarlaudM1996Conjugate-Gradient Method for Soliving Inverse Scattering with Experimental DataIEEE Antennas Propag. Mag., 38
  23. 23. OliveriGLizziLPastorinoMMassaA2012A Nested Multi-Scaling Inexact-Newton Iterative Approach for Microwave ImagingIEEE Trans. Antennas Propag., 60page numbers (971 EOF983 EOF
  24. 24. PastorinoM2010Microwave imagingJohn Wiley, Hoboken N.J.
  25. 25. PaulsenK. DPoplackS. PLiDFanningM. WMeaneyP. M2000A clinical prototype for active microwave imaging of the breastIEEE Trans. Microw. Theory Tech., 48page numbers (1841 EOF1853 EOF
  26. 26. RekanosI. T2008Shape Reconstruction of a Perfectly Conducting Scatterer Using Differential Evolution and Particle Swarm OptimizationIEEE Trans. Geosci. Remote Sens., 46page numbers (1967 EOF1974 EOF
  27. 27. RemisR. Fvan den Berg, P. M. (2000On the equivalence of the Newton-Kantorovich and distorted Born methodsInv. Probl., 16page numbers (L1 EOFL4 EOF
  28. 28. RichmondJ1965Scattering by a Dielectric Cylinder of Arbitrary Cross Section Shape. IEEE Trans. Antennas Propag., 13page numbers (334-341).
  29. 29. RoccaPBenedettiMDonelliMFranceschiniDMassaA2009Evolutionary Optimization as Applied to Inverse Scattering ProblemsInv. Probl., 25
  30. 30. SchilzWSchiekB1981Microwave Systems for Industrial Measurements. Advances in Electronics and Electron Physics, 309381Elsevier.
  31. 31. SheaJ. DKosmasPVan VeenB. DHagnessS. C2010Contrast-Enhanced Microwave Imaging of Breast Tumors: A Computational Study Using 3D Realistic Numerical PhantomsInv. Probl., 26
  32. 32. van den BergP. M., & Abubakar, A. (2001Contrast Source Inversion Method: State of Art. Prog. Electromag. Res., 34page numbers (189-218).
  33. 33. VogelC. R2002Computational methods for inverse problemsSociety for Industrial and Applied Mathematics, Philadelphia.
  34. 34. ZhangZ. QLiuQ. H2004Three-Dimensional Nonlinear Image Reconstruction for Microwave Biomedical Imaging.IEEE Trans. Biomed. Eng., 51page numbers (544 EOF8 EOF
  35. 35. ZhouHTakenakaTJohnsonJTanakaT2009A Breast Imaging Model Using Microwaves and a Time Domain Three Dimensional Reconstruction MethodProg. Electromag. Res., 93page numbers (57 EOF70 EOF
  36. 36. ZoughiR2000Microwave Non-Destructive Testing and EvaluationKluwer Academic Publishers, Dordrecht.

Written By

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

Submitted: 07 June 2012 Published: 14 November 2012