Compressed Sensing and Defect-Based Dictionaries for ...

0 downloads 0 Views 375KB Size Report
achieved by applying compressed sensing techniques. In this work, using effective modelling of possible defects, such as air gaps between layers, we constructΒ ...
Compressed Sensing and Defect-Based Dictionaries for Characteristics Extraction in mm-Wave Non-Destructive Testing Edison Cristofani1,2, Mathias Becquaert1,2, Gokarna Pandey2, Marijke Vandewal1, Nikos Deligiannis2,4 and Johan Stiens2,3 1

CISS Department, Royal Military Academy, 30 Av. de la Renaissance, 1000 Brussels, Belgium 2 ETRO Department, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium 3 SSET-department, Imec, Kapeldreef 75, 3001 Leuven, Belgium 4 iMinds, 9050 Ghent, Belgium

Abstractβ€”In ultra-wideband non-destructive testing of large multilayered polymers, data collection and reduction can be achieved by applying compressed sensing techniques. In this work, using effective modelling of possible defects, such as air gaps between layers, we construct defect dictionaries and use them as support data for a signal similarity-based classifier, which will automatically extract the main characteristics of the inspected defect.

I. INTRODUCTION

N

ON-Destructive Testing (NDT) of industrial materials or parts is a fundamental aspect in maintenance to guarantee and maximize life-span of larger structures. Teflon, ABS or fiberglass can be found in critical parts which undergo strong mechanical stress such as wind turbine blades, aircraft wings, or certain gas and oil pipelines. These materials present a very high transparency at mm-wave frequencies enabling see-through testing capabilities [1]. Although testing at higher mm waves is an attractive option, parts of interest can be relatively large, thus performing any type of measurement becomes cumbersome if not impractical. This situation can be even more accentuated when inspection systems produce resolutions down to the cm or mm. A reduction in the measurement time and data collection burden is therefore a desirable solution, provided that inspection performance remains unaffected and the material characteristics are properly extracted. II. COMPRESSED SENSING AND ULTRA-WIDEBAND NDT Sensors providing millimeter scanning accuracies generate enormous amounts of data due to the high number of measurement positions required. Achieving full signal reconstruction using lower sampling rates can be decisive when it comes to guaranteeing NDT using Ultra-Wideband (UWB) sensors. The sparse sampling theorem that forms the crux of Compressed Sensing (CS) proves to be an effective technique for reducing the amount of collected data, assuming that the signal of interest has a sparse representation in a given basis or learned dictionary [2]. Let π‘₯ ∈ ℂ𝑁 be a vector containing a complete measurement and 𝑦 ∈ ℂ𝑃 be a partial or incomplete observation satisfying 𝑃 < 𝑁. Let 𝑨 be a 𝑃 Γ— 𝑁 matrix defining the relation 𝑦 = 𝑨π‘₯, which corresponds to an underdetermined linear system with infinite solutions. This system may only be solved by imposing a priori knowledge of unknown vector π‘₯ being sparse, that is, most of its elements are zero. Additionally, matrix 𝑨 must satisfy one of the following conditions: the Restricted Isometry Property (RIP), the null space property, and the mutual coherence property [4]. In order to solve this problem, a minimization problem can be formulated as

π‘₯Μ‚ = arg minβ€–π‘₯β€–0 𝑠. 𝑑. 𝑨π‘₯ = 𝑦, π‘₯

(1)

where β€–π‘₯β€–0 is the β„“0 -pseudo-norm of vector π‘₯, which imposes sparsity. The two described sets of constraints become impractical but can be relaxed using the following considerations [5]: i. ii.

The β„“0 -pseudo-norm makes the minimization in (1) an NP-hard problem. The β„“1 -norm is used instead. Satisfying the constraints for 𝑨 is unrealistic for large datasets. Populating matrix 𝑨 with randomly generated elements obeying a Gaussian distribution allows vector recovery with high probability.

In the application hereby proposed, we consider the vector π‘₯ to be the fully sampled vector according to the ShannonNyquist sampling theorem [6]. In a conventional, noiseless measurement, matrix 𝑨 would be the identity matrix of size 𝑁 Γ— 𝑁 with vectors π‘₯ and 𝑦 of length 𝑁. For practical reasons, 𝑨 cannot be a random matrix as stated in consideration ii) and random sampling is implemented by randomly selecting subsets of the identity matrix. The sparsity condition may not always be satisfied in the measurement domain and, thus, a 𝑁 Γ— 𝑁 transformation matrix 𝚿 is added to the problem statement in (1) resulting in 𝑦 = π‘¨πšΏπœŽ, where 𝜎 = 𝚿 T π‘₯ is defined as the transformed and sparsified representation of π‘₯. We solve the proposed problem by applying CS algorithms based on orthogonal matching pursuit [7] or convex relaxation [8]. UWB sensors can benefit from CS principles in two ways: Firstly, lower sampling rates produce less data to be stored or transferred, and thus a simpler, less hardware-demanding detection chain is required [3], [9]. Secondly, CS is of special interest when using UWB stepped-frequency systems for NDT, for which discrete frequencies within the operation bandwidth are measured multiple times to stabilize the signal and improve Signal-to-Noise Ratio (SNR). Therefore, reducing measured frequencies thanks to CS implies reducing measurement times as well. III. SCENE DESCRIPTION AND EXPERIMENTAL SETUP Air gaps produced by defective gluing during the production process are common artifacts in multilayered materials and can affect structural integrity at a larger scale. In order to understand the possibilities of applying mm-wave inspection and CS, two 20-mm (𝑀1 = 𝑀2 = 20 π‘šπ‘š) ABS pieces (refractive indices 𝑛1 = 𝑛2 β‰ˆ 1.7 [10], see Fig. 1) are produced to simulate controlled air gaps 𝑔 ranging from 1 mm to 15 mm and including the no-defect case, or zero-thickness air gap. These configurations are measured with a mm-wave vector network analyzer operating a in the 40 βˆ’ 75 𝐺𝐻𝑧 range and a frequency step of 5.83 𝑀𝐻𝑧.

Random DS factor

0.9

sensor

0.2

0.7 0.4

0.6 0.5

0.6

0.4 0.3

0.8 1 -10

In-depth dimension

Fig. 1: Schema of the proposed measurement setup. Two materials 𝑀1 and 𝑀2 are placed in parallel at a certain distance (air gap).

0.8

0.2 0.1 0

10 SNR (dB)

20

30

Fig. 2: Successful reconstruction rates of two-layer measurements over 50 iterations. 𝐷𝑆 varies between 0.05 and 1. Operation point for typical SNR yields a minimum DS factor of 0.2.

The in-depth resolution in free space is defined as π‘π‘œ ⁄2𝐡 β‰ˆ 4 π‘šπ‘š, where π‘π‘œ is the speed of light in free space and 𝐡 is the system’s bandwidth.

Tests on real measurements show that sparse data recovery using CS theory allows extraction of the thickness and location of air gaps within the explored materials by reconstructing the media interfaces. To achieve such results, we use a Fourier sparsifying matrix 𝚿 with the purpose to increase sparsity since the data recorded directly from the sensor is not sparse enough. Fig. 2 shows actual results for successful air gap reconstruction (1 perfect match; 0 no match) for a range of random downsampling factors and SNR. The downsampling factor is defined as 𝐷𝑆 = π‘˜/𝑁, where π‘˜ is the number of randomly measured samples. The no operation region can be clearly identified in the left-hand-side of the figure, due to low SNR and/or low 𝐷𝑆. For a typical measurement under controlled conditions, SNR may vary between 15 and 25 dB. Thanks to applying CS principles, our approach enables a downsampling factor as low as 0.2 without lessening the overall sensor performance. Fig. 3 depicts a successful reconstruction of a scene presenting a 9-mm air gap under typical SNR (20 𝑑𝐡) and keeping only 20 % (𝐷𝑆 = 0.2) of the original samples. Although the amplitudes of the media interfaces show a slight mismatch, the locations are perfectly determined. Thus, by choosing the adequate parameters, our CS-based approach enables faster exploitation of the see-through capabilities of UWB sensors in the mm-wave range while preserving system performances. Regarding future work, automatic air gap determination thanks to robust classification via the β„“1 minimization [11] will be possible by using either real observations or a simulator to model and generate ideal observations of air gaps that will populate a dictionary of possible defects in (1). This technique, which operates like a similarity-based classifier, can be then used during the sparse data reconstruction process to recover not only the original sensed signal data but also to determine which of the elements in the dictionary provides the highest likelihood. Combining this approach with different sparsifying bases and dictionary learning techniques could reduce even more the amount of required data for perfect measurement reconstructions.

Arbitrary amplitude

IV. RESULTS AND FUTURE WORK

1

Original signal Reconstructed signal

0.9 0.8 0.7 0.6 0.5

In-depth dimension

Fig. 3: Reconstruction of a 9-mm air gap using a two-layer configuration in ABS and typical SNR and 𝐷𝑆 = 0.2. Interfaces are perfectly located.

REFERENCES [1]

E. Cristofani, A. Brook, M. Vandewal, C. Matheis, and J. Jonuscheit. Assessment of 3D signal and image processing using FMCW THz signals. In Proc. of the OPTRO 2012, 5th International Symposium on Optronics in Defense and Security, Paris, France, February 2012. [2] D. L. Donoho, β€œCompressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006. [3] E. Cristofani, M. Becquaert, and M Vandewal. Performance of 2D compressive sensing on wide-beam through-the-wall imaging. Journal of Electrical and Computer Engineering, 2013(636972), 2013., p1-11, 2013. [4] E. J. CandΓ¨s, J. K. Romberg, and T. Tao, β€œStable signal recovery from incomplete and inaccurate measurements,” Communications on Pure and Applied Mathematics, vol. 59, no. 8, pp. 1207–1223, 2006. [5] H. Rauhut, β€œCompressive sensing and structured random matrices,” in Theoretical Foundations and Numerical Methods for Sparse Recovery, M. Fornasier, Ed., vol. 9 of Radon Series on Computational and Applied Mathematics, pp. 1–92, deGruyter, Berlin, Germany, 2010. [6] C. E. Shannon, β€œA mathematical theory of communication,” The Bell Systems Technical Journal, vol. 27, pp. 379–423, 623–656, 1948. [7] M. A. Davenport and M. B. Wakin, "Analysis of Orthogonal Matching Pursuit Using the Restricted Isometry Property," in IEEE Transactions on Information Theory, vol. 56, no. 9, pp. 4395-4401, Sept. 2010. [8] Y. C. Eldar and H. Rauhut, "Average Case Analysis of Multichannel Sparse Recovery Using Convex Relaxation," in IEEE Transactions on Information Theory, vol. 56, no. 1, pp. 505-519, Jan. 2010. [9] M. Becquaert, E. Cristofani, G. Pandey, M. Vandewal, J. Stiens, N. Deligiannis, β€œCompressed sensing mm-wave SAR for non-destructive testing applications using side information”, in IEEE Radar Conference, PA, May 2016. [10] Berins, M.L., SPI Plastics Engineering Handbook of the Society of the Plastics Industry, Inc., New York: Van Nostrand Reinhold, 1991 [11] J. Wright, A. Y. Yang, A. Ganesh, S. Sastry, Y. Ma, β€œRobust face recognition via sparse representation.” IEEE transactions on pattern analysis and machine intelligence, vol. 31, no. 2, pp. 210-227, 2009.