Next Article in Journal
Automatic User Preferences Selection of Smart Hearing Aid Using BioAid
Previous Article in Journal
Enteral Nutrition and Hydration in Patients with Acute Stroke: Efficacy of an Automatic Pump System for Water Administration and Flushes—A Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on Millimeter Wave SAR Imaging for Non-Destructive Testing of Rebar in Reinforced Concrete

1
Department of Smart Information Technology Engineering, Kongju National University, Cheonan 31080, Korea
2
Department of Green Smart Architect Engineering, Kongju National University, Cheonan 31080, Korea
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(20), 8030; https://doi.org/10.3390/s22208030
Submission received: 13 September 2022 / Revised: 10 October 2022 / Accepted: 19 October 2022 / Published: 20 October 2022
(This article belongs to the Section Electronic Sensors)

Abstract

:
In this study, we investigate a millimeter wave (mmWave) synthetic aperture radar (SAR) imaging scheme utilizing a low-cost frequency modulated continuous wave (FMCW) radar to take part in non-destructive testing which could be a useful tool for both civilian and military demands. The FMCW radar working in the frequency range from 76 GHz to 81 GHz is equipped with a 2-D moving platform aiming to reconstruct the 2-D image of the shape of the target object. Due to the lab environment containing several devices and furniture, various noise and interference signals from the floor are not avoidable. Therefore, the digital signal processing algorithms are joined to remove the undesired signals as well as improve the target recognition. This study adopts the range migration algorithms (RMAs) on the processed reflected signal data to form the image of the target because of its verified ability in this type of mission. On the other hand, the integration of compressed sensing (CS) algorithms into the SAR imaging system is also researched which helps to improve the performance of the system by reducing the measurement duration while still maintaining the image quality. Three minimization algorithms are used involving the imaging system as the CS solvers reconstruct the radar data before being processed by RMA to form the image. The proposed imaging scheme demonstrates its good ability with high azimuth resolution in the mission of detecting tiny cracks in the rebar of reinforced concrete. In addition, the participation of CS algorithms improves the performance of the scheme as the cracks on the rebar can be located on the images, which are reconstructed from only 30% of the dataset. The comparison of CS solvers shows that ADMM outperforms the other candidates in the reconstruction task.

1. Introduction

According to international regulations, to ensure the serviceability of any building materials, products, and elements over their entire useful lives, their specific physical and strength properties should be defined [1]. Hence, the constant monitoring of the structures such as lifetime assessment of industrial construction plays a crucial role in the guarantee of safety [2]. Although it has the disadvantage of being an indirect method, non-destructive testing (NDT) is widely applied in civil engineering because of its beneficial characteristics, including its non-destructive nature, repeatability, and economical aspect. NDT methods can conduct a condition assessment to provide information for the structural performance of the material, e.g., equipment integrity analysis, corrosion monitoring of structures and equipment, corrosion damage evaluation, fatigue and creep damage prediction, and fitness-for-service evaluations [2]. There are various types of NDT methods employed in routine services for the integrity of structural component detection, such as radiography testing, eddy current testing, ultrasonic, and resonant near-field probes. However, these methods not only require expertise but have their own limitations. For example, the operator safety factor must be taken into account as they must pay careful attention while using the radiography method [3]. On the other hand, eddy current is known as one of the effective NDT techniques for testing the steel members [4,5], which is based on the change of voltage or current in the sensing coil. Nevertheless, a highly skillful and experienced inspector is one of the main requirements. Additionally, some special signals occur as the defect at the root of threads and the signal lags severely as the increment of the detection depth [6,7] makes them hard to process. The ultrasonic method is based on measuring the transit time that is required for the ultrasonic pulse to go through the material several times. Hence, the knowledge about the precise physical characteristics of the material, such as the ultrasonic velocity, plays an important role [8] which can lead to inaccurate inspection results [9].
Recently, SAR imaging has become one of the most effective NDT methods and appealed to more and more attention not only in ground or sea detection but also short distance applications including indoor detection or structure health monitoring (SHM) since it takes advantage of the natural characteristics of the electromagnetic (EM) wave with good penetration performance. In general, radars are utilized to emit the EM wave for both civilian and military applications in the form of impulse radars and continuous wave (CW) radars. In terms of CW radars, the FMCW radar is preferable as this type of radar is not restricted to military applications as well as it can be adapted for commercial tasks such as measuring material thickness [10], monitoring terrain displacement [11,12], see-through-wall applications [13], life activity monitoring [14], as well as automotive applications in collision avoidance [15], controlling adaptive cruise [16], and all-weather cruise [17]. Moreover, FMCW radars are extensively employed in imaging applications. The combination of FMCW radar and SAR imaging techniques results in a low-cost, effective imaging sensor with high resolution. The research in [18] also studies the application of SAR imaging in detecting the rebar in the reinforced concrete but due to the low resolution of the imaging scheme, it can only identify the presence of the rebar without recognizing any defects. The authors in [19] proposed a 2-D array antenna system to reconstruct the SAR images of small structures, which may be a highly cost-effective method. In [20], a scheme employing FMCW SAR imaging for the NDT application to detect a defect on the 3-D printed object was presented. However, a huge number of measurements in the experiment can cause a time consumption problem.
The appearance of CS in the previous two decades [21,22] has brought many CS-based potential applications not only for image processing but also for signal processing, and SAR imaging is also one of the candidates [23,24]. The author in [25] proposed an intensity-hue-saturation-based (IHS-based) pan-sharpening method using ripplet transform and CS to solve the evident spectral distortion problem caused by the IHS transform. In this research, the CS is adopted to reconstruct the intensity component to integrate the local information from both the intensity component and the panchromatic (PAN) image. The research in [26] studies the integration of CS in spare reconstruction on the randomly sub-sampled raw data of automobile FMCW SAR for far-field imaging. Their CS scheme is based on the basis pursuit denoising (BPDN) solver to recover the Ku-band FMCW signal data. CS was also utilized by the study in [27] aiming to archive the high azimuth resolution for the low-THz band using the convex solver. By applying the CS method, the SAR image can be reconstructed even though the radar data are acquired in an irregular spatial interval. The CS-SAR imaging scheme can form the images of sparse or compressible targets from a randomly under-sampled data collection and thus it can help to reduce the number of measurements efficiently. The study in [28] applies the CS to recover the FMCW signals before SAR processing from two of the in-phase and quadrature (I/Q) components using the l 1 - n o r m minimization [29], which may take a considerable duration of the recovery process. In this study, we adopted the RMA algorithm as the SAR imaging construction technique which has been proven as an effective approach in near-field SAR imaging [30] to identify the location of a crack in the steel rebar which is used in reinforced concrete. The 2-D shape of the object is formed in the image result, which causes numerous measurements to accumulate adequate signal data. Therefore, by directly applying CS to the complex signal data gathering from the measurement platform, this study investigates a method to reduce the great number of measurements and the side-lobe effect caused by the non-uniform spacing of measurement points [31], which result in saving the measurement and the processing period while maintaining the quality of the images. Additionally, we investigated three different CS-based solvers, namely the alternating direction method of multipliers (ADMM), iterative reweighted least squares (IRLS), and coordinate descent (CD) to recover the signal data from small random measurement collections. To maintain the image quality and reduce environmental interference, the finite impulse response (FIR) filter along with the smoothing window and signal-to-noise improvement techniques are employed. The quality of the images formed from the data recovered by the CS-based methods was compared to that of the original image which was reconstructed by utilizing the entire data. In addition, the effectiveness of different CS-based algorithms is also revealed.
This paper is organized as follows: in Section 2, the methodology of the proposed imaging scheme is presented. Section 3 describes the experiment setup along with the results. Finally, Section 4 gives the summarization and conclusion of this paper as well as the future works to enhance the system performance.

2. Theory and Methodology

In this section, we describe the theoretical approaches which are employed in our study. We first present the FMCW echo containing the reflectivity information of the target scatterer. Then, the method of adopting the signal data to reconstruct the SAR image by RMA is depicted followed by the compressed sensing algorithms to improve the performance of our imaging scheme. The overall process of the proposed scheme is illustrated in Figure 1.

2.1. FMCW Signal and RMA Approach

The characteristics of the FMCW signal that can be utilized by SAR imaging are presented in this section. As shown in Figure 2, an ( x , y , z ) Cartesian coordinate system is established consisting of the x-axis, y-axis, and z-axis, which denote horizontal, vertical, and range directions, respectively. We assume that a single scatterer located at position ( x , y , z 0 ) with a complex reflectivity ρ ( x , y ) while the location of the FMCW radar is ( x R , y R , 0 ) . The radar sensor moving on the Oxy plane results in the rectangle grid of measurement points. The radar is considered to carry a single full-duplex antenna at the mid-point between the transmitting and receiving antennas since they are so close to each other [32]. In the FMCW SAR imaging configuration, by performing the data collection at all measurement points, a 3-D (three-dimensional) data cube can be accumulated as the backscattered signal data. The intermediate frequency (IF) signal at each measurement point is denoted as:
s I F ( x , y , t ) = ρ ( x , y ) R 2 exp j 2 π ( f 0 τ + α t τ 1 2 α τ 2 )
where R is the distance between the radar sensor and the scatterer, f 0 is the carrier frequency at 0, τ is the time delay between the transmitted and the received signal with the speed of light c, and α is the slope of the frequency in the chirp duration T with the bandwidth B. It is clearly seen that the IF signal here is a sinusoidal signal with the beat frequency f b = α τ which is proportional to the time delay ( τ ) so that it can provide the range information. The residual video phase (RVP), which is the last term in (1), is negligible and can be ignored [33]. Therefore, the IF signal can be redefined in the wavenumber domain as:
s I F ( x , y , k ) = ρ ( x , y ) exp j 2 k R R 2 , 2 π f 0 c k 2 π f T c
where f T = f 0 + α T is the maximum swept frequency of FMCW radar and is the corresponding wave number at each specific frequency f.
As a 2-D target can be characterized by utilizing the reflectivity function ρ ( x , y ) , the RMA employs the ρ ( x , y ) at all measurement points to form the image of the target object. Here we assume that the measurement system continuously captures data over the O x y plane which covers the target dimension. Hence, the backscattered signal data from the target plane at z 0 distance can be derived by taking the integral of (2)
s ( x , y , k ) = ρ ( x , y ) exp ( j 2 k R ) R 2 d x d y
where the distance between the transceiver antenna and the general point on the target at distance z 0 from the aperture plane is R = ( x x ) ( 2 y y ) 2 + z 0 2 . As the considered single full-duplex antenna, R 2 in (3) is approximately ( z 0 R 1 ) . Due to the near-field application adopted in this study, the error in this approximation can be negligible [34]. In addition, z 0 is the constant and as the target is stationary, it can also be ignored. As a result, the signal data for image construction will be processed by a combination of phase terms and the remaining R 1 dependence. Now, the backscattered signal is
s ( x , y , k ) = ρ ( x , y ) exp ( j 2 k R ) R d x d y
According to Weyl’s study in [35], there is a representation of a spherical wave as the superposition of plane waves. In the case of the homogeneous plane waves, the backscattered signal can become
s ( x , y , k ) = j 2 π ( k x , k y ) × exp ( j k z z 0 ) k z exp ( j ( k x x + k y y ) ) d k x d k y
where
( k x , k y ) = ρ ( x , y ) exp ( j ( k x x + k y y ) ) d x d y = FT 2 D [ ρ ( x , y ) ]
The backscattered data can be inferred by using 2-D IFT as
s ( x , y , k ) = IFT 2 D ( k x , k y ) exp ( j k z z 0 ) k z
After dropping the constant components, it yields that
( k x , k y ) = S ( k x , k y , k ) exp ( j k z z 0 ) ,
where S ( k x , k y , k ) = FT 2 D [ s ( s , y , k ) ] . To extract the complex reflectivity at a point on the O x y plane, the 2-D IFT is applied to (8). By coherently summating all complex reflectivity at all points within the target plane, the 2-D image is constructed
ρ ( x , y ) = IFT 2 D S ( k x , k y , k ) exp ( j k z z 0 ) d k

2.2. Data Collecting from 2-D Scanning System and RMA Involvement

As mentioned in Section 2.1, the radar scans on the 2-D aperture plane. Let us denote the positions of the platform as ( x 0 , y 0 ) and ( x H N , y V M ) at the initial and final points on the aperture plane, respectively, where H N and V M denote the number of measurement points on each direction. The FMCW radar moves along the horizontal and vertical directions on the x-axis and y-axis in the zigzag trajectory, as shown in Figure 3. In this figure, the radar is moved with the horizontal and vertical steps of Δ x and Δ y , respectively. The measurement step must satisfy the Nyquist sampling which is dependent on the aperture size ( D x , D y ) , target dimension ( D t h and D t W are the height and the width of the target, respectively), and the distance between the antenna and the target ( z 0 )
Δ x λ ( D x + D t w ) 2 / 4 + z 0 2 2 ( D x + D t w ) , Δ y λ ( D y + D t h ) 2 / 4 + z 0 2 2 ( D y + D t h ) ,
where λ is the wavelength. To maintain the quality of the image, the spatial sampling intervals should be based on the worst cases of the Nyquist sampling where the z 0 0 , the Δ x = Δ y λ / 4 . However, to achieve high-quality imaging but with fewer measurements, these sampling intervals can be larger than the worst case of the Nyquist criterion or we can apply the CS-based method, which is presented in the next section. Since we are considering the 2-D imaging on the aperture plane, the range resolution ( δ z ) is not so important. Instead, we take the cross-range resolutions (or azimuth resolutions) into account, which can be expressed as:
δ x λ 2 D x , δ y λ 2 D y ,
where δ x and δ y correspond to the horizontal and vertical resolutions, respectively.
The 3-D data cube of the FMCW echo is accumulated after measuring all measurement points. Let us denote the received signal at the measurement point ( x i , y j ) as s ( i , j ) , and each chirp contains N s samples, the collection of chirp signal data at all point having size N s × V M × H N can be formed as
S = s ( H N , 1 ) s ( H N , 2 ) s ( H N , V M ) s ( 1 , 2 ) s ( 2 , 2 ) s ( 2 , V M ) s ( 1 , 1 ) s ( 1 , 2 ) s ( 1 , V M )
In this study, we applied the FIR filter and window smoothing method to reduce interference from the environment. Besides that, we also apply an integral technique to improve the signal-to-noise ratio [36] as
s ( x , y , k ) = 0 T s ( x , y , t ) exp ( j 2 π f b t ) d t
where f b is the beat frequency of the IF signal at the corresponding distance z 0 . Now, the RMA is applied by firstly transforming the data using range-FFT (fast Fourier transform) which can be expressed as F ( S ) . After that, we perform the range focusing on the target to get the planar data at distance z 0 , denoted as M = F ( S , k z 0 ) with the dimension of V M × H N . Next, we define the matched filter as
MF = exp ( j 2 k x 2 + y 2 + ( z 0 ) 2 ) .
The image can be constructed as
M ^ = IFT 2 D S t o l t FT 2 D ( M MF )
where ⊗ denotes the element-wise multiplication operator for matrices. We adopt STOLT interpolation here into the image reconstruction to take the uniformly distributed data at the k z domain.

2.3. Compressed Sensing

In this section, we describe the participation of CS in the image scheme to improve the performance of the system. The involvement of CS contributes two main benefits to the imaging scheme. First, during the measurement process, the radar collects the echo signal by scanning on the 2-D aperture plane consisting of numerous measurement points which leads to a significantly long duration to accumulate all the signal data. With CS, the number of measurement points can be considerably decreased as the radar only emits and receives the signal at certain points as on the random measurement matrix. Secondly, CS helps to diminish the side-lobe effect caused by the non-uniform spacing of measurement points [31] which can bring about the ghost effect on the image result. The details of CS algorithms are presented in the next subsections.

2.3.1. CS Theory

Signal acquisition is the main topic in signal processing. Sampling theorems provide the bridge between the continuous and the discrete-time worlds. The most famous theorem is often attributed to Shannon (but is usually called the Nyquist theorem) and states that the sampling rate must be twice the maximum frequency present in the signal to perfectly recover the signal. In 2004, David Dohono [21] proved that a signal can be reconstructed with fewer samples than the sampling theorem, which introduced the concept of compressive sampling. CS allows compressing signals with sparse or compressible representation while they are sampled. It initiates from the idea that it is not necessary to invest a lot of power into observing the entries of a sparse signal because most of them will be zero. The fundamental principle of the CS method can be described in Figure 4, which contains two phases, namely compress and recovery operations.
(1) A signal x with length N ( x R N × 1 ) needs to be reconstructed from M measurements such that M N .
(2) The signal is then made sparse through some transform matrix Ψ as:
x = Ψ s
where s is the k sparse coefficient vector of size N × 1 ( s R N × 1 ) , and Ψ is the transform matrix whose size of N × N .
(3) With knowledge of the sparse vectors, it is possible to reconstruct signal x from (16). Thus, the goal of compressed sensing is to find the sparsest vector s that is consistent with the measurements vector. The measured signal vector is calculated as:
y = C x = C Ψ s = θ s
where y is the measurement vector or the compressed signal of size M × 1 , θ = C Ψ is the random sensing matrix of size M t i m e s N with C R M × N being the measurement matrix. It is important to note that the measurement matrix C is not selected randomly, it must possess the RIP (restricted isometry property) condition, which means that it is guaranteed to hardly change the length of vector x as long as the vector x is at least k-sparse [37]. In [38], researchers show that the measurement matrix C adopting a random Gaussian matrix can probably satisfy the RIP condition. Another study [39] employs the random Bernoulli matrix as it extends the symmetric signs ensemble which also satisfies RIP with high probability.
(4) The reconstruction of the original signal can be seen as the process of solving y = θ s . As only the sparsity k is given, this equation is an under-determined equation system. Hence, s then be recovered by finding the optimal l 0 - n o r m to minimize the number of non-zero items in s which is a difficult problem. However, if we turn this optimization problem into an l 1 - n o r m problem, it can be solved by using convex optimization [29]. The recovering process of s is as follows:
s * = arg min θ s y 2 < ϵ
where . p is the p norm and ϵ is the residual.

2.3.2. SAR Imaging with Compressed Sensing

For the application of CS, the sparse data are formed by reducing the dataset of reflected signals, which means that a smaller number of measurement points should be randomly taken in the aperture plane. In this paper, we selected the measurement points by applying the random Bernoulli matrix as the measurement matrix C. Let us express the data matrix with irregular data acquisition as S ˜ , it can be seen as a matrix in which random positions in S are set to zero. As most of the elements in S ˜ are zeros, we can recover the data by using the l 1 m i n i m i z a t i o n method. Assuming the size of S ˜ is N s × V ˜ M × H ˜ N , where V ˜ M V M and H ˜ N H N . We also assume that at the ith sample of all received signals at all measurement points, the received signal vector is expressed as S ˜ ( u , v , w = i ) . The signal vector s * at the ith sample can be restored by solving the minimization problem. In this study, three methods, namely ADMM, IRLS, and CD, are investigated to solve our problem. These methods can be expressed as below:
  • The ADMM algorithm has been shown as a simple, effective, and fast convergent approach [40,41]. The optimization problem applied by the ADMM algorithm is given as
    s * = arg min 1 2 θ s y 2 + λ z 1 subject to s = z
    The regularization parameter λ (> 0) controls the balance between the data fidelity term 1 2 y θ s 2 and the regularization term z 1 . The update equations of ADMM over iteration for (19) are given by
    s [ k + 1 ] = arg min 1 2 θ s y 2 + η 2 s z [ k ] + w [ k ] 2
    = θ T θ + η I 1 θ T + η ( z [ k ] + w [ k ] ) z [ k + 1 ] = arg min λ z 1 + η 2 s [ k + 1 ] z + w [ k ] 2
    = prox ( λ η f ) s [ k + 1 ] + w [ k ]
    w [ k + 1 ] = w [ k ] + s [ k + 1 ] z [ k + 1 ]
    where k ( = 0 , 1 , 2 , ) is the iteration index in the algorithm, ( . ) T is the transpose of the matrix, I is the identify matrix, and η ( > 0 ) is the parameter. The initial values of w [ 0 ] and z [ 0 ] are 0.
  • In the IRLS method, the weight least square in each iteration is used to infer the next estimate with the weights derived from the last iteration. The optimization problem in each iteration is as follows: s * n + 1 = arg min s T W n s subject to y = θ s , n 1 , where W is the diagonal weight matrix. At the kth iteration, this matrix is computed from the solution of the current iteration s k , W i k = | s i k | 1 , 1 i N . The closed form solution for s k + 1 can be inferred as s k + 1 = W k 1 θ T θ W k 1 θ T y .
  • The coordinate descent method can be applied in multi-variable minimization by solving a sequence of scalar minimization sub-problems. By minimizing each sub-problem along a selected coordinate while all other coordinates are fixed, the estimate of the solution is improved [42].

3. Experiment Setup and Results

In this chapter, the measurement system for FMCW SAR imaging is presented including the descriptions of all the hardware components and their integrations into the system to collect echo data from the target. After that, the system configuration of the radar and the experiment setups for collecting data are also introduced, followed by their experiment results.
The mmWave radar in our system is the Texas Instrument AWR1443Boost, which is an evaluation module based on a single-chip AWR1443 mmWave sensor. There are three transmit antennas and four receive ones integrated into the radar. It can emit an FMCW signal with a 4 GHz available bandwidth in the frequency range from 76 GHz to 81 GHz [43]. Besides that, it is also configured to be interfaced with the DCA1000 evaluation module, as shown in Figure 5a. This DCA1000 module is a real-time data capture board for interfacing with Texas Instrument’s mmWave sensors. The raw ADC data recorded by AWR1443Boost are then captured by the DCA1000 module and streamed to the host PC over Ethernet as packetized data. The designed moving platform in our system has two orthogonal linear actuators, as shown in Figure 5b.

3.1. Experiment Setup

In this study, we take into consideration the feasibility of integrating the system into the NDT applications. To archive that, we conduct an NDT experiment by taking the involvement of the reinforced concrete as the MUT (material under test). This experiment strives to locate the cracks in the steel rebar inside the concrete samples. We conduct the measurements on three reinforced concrete samples, one sample has two gaps (cracks) on the rebar with dimensions of 1 mm and 3 mm, and another sample contains two cracks of 2 mm and 4 mm on its rebar (we named them sample A and sample B). To make the concrete sample easier to produce, every crack is filled with a silicon layer whose low reflection ability is compared to the steel rebars. Each rebar is 150 mm long and its diameter is 13 mm. The distance between the surface of the concrete to its rebars is around 2 mm. Figure 6 shows the structure of the reinforced concrete samples with a crack on their rebars.
Based on these characteristics of the concrete, we consider D t w = 150 mm and D t h = 13 mm. Therefore, to make sure the coverage of the scanning region on the target, the dimension of the scanning plane is defined as D x = 150 mm and D y = 40 mm. The movement steps in both horizontal and vertical directions are Δ x = 1 mm and Δ y = 2 mm, respectively, which results in 3000 measurement points including 150 points in the horizontal direction and 20 points in the vertical one. Because of the concrete layer covering the rebar, we decide to utilize the maximum antenna gain, which is set to 48 dB. To validate the different azimuth resolutions on the image which play a crucial role in detecting hidden cracks, we carry on the experiments on two different values of z 0 , 100 mm and 50 mm. The experiment setup for detecting cracks in the rebar of concrete is shown in Figure 7.

3.2. Experiment Results

According to the theory of the azimuth resolution in (11), as the range between the radar to the target ( z 0 ) increases, the resolution is decreased correspondingly. Our image results in Figure 8 confirm the negative correlation between these two factors in the SAR algorithm. These range values of z 0 at 50 mm and 100 mm and experiments configuration in Section 2.2 theoretically lead to the theoretical resolutions of ( δ x 1 mm, δ y 2.4 mm) and ( δ x 1.3 mm, δ y 4.9 mm), respectively. For all the experiments, the same dimensions of the 3-D data cubes are captured, 512 × 20 × 150 corresponding to N s × V M × H N . Figure 8a,b illustrates the image results of the two first experiments on concrete samples where the distance from the FMCW radar to the steel rebar ( z 0 ) is 50 mm. It can be seen that the proposed imaging scheme can determine the positions of the rebars hidden within the concretes. In terms of crack detection application, these experiments show inefficient abilities since it is only able to determine the position of 3 mm and 4 mm gaps, while the small cracks whose sizes of 1 mm and 2 mm do not appear in these figures. Even though the theoretical horizontal resolution for this experiment is around 1 mm, the ghost effect on the image that comes from multi-static measurements makes it so the tiny gaps are not located. In addition, another reason for this result is that the reflection of the silicon layers filled in the gaps can be shown as parts of the rebars on the SAR images. Although the weak echo is reflected from those layers, the gaps cannot be distinguished. The SAR image results of the experiments carried out on the same concrete samples but different from the z 0 = 100 mm are illustrated in Figure 8c,d. As shown in these images, although there are weaker reflections at the positions of the cracks than in the other parts of the rebars, it is difficult to conclude whether these positions contain the cracks or not without reference to the real structure of the concrete. There are two main factors that can explain this result. The first one is the lower theoretical resolution from the experiment setup, while the second one is the echo from the silicon layer within the gaps as in the previous experiment on concrete samples.
Next, we validate the integration of CS algorithms into the imaging scheme. We first take smaller parts of the dataset from the original one and then try to reconstruct the SAR image from the recovered data by adopting CS. We apply the Bernoulli random matrix ( C ) to take the amount of 10%, 20%, 30%, 40%, and 50% from the original data. We use the discrete cosine transform (DCT) as the transform matrix ( Ψ ) . After that, the three aforementioned CS-based solvers are used, which involved restoring the data from these partitions. The performance of these three solvers is depicted in Figure 9, including the evaluation of their reconstruction errors and their convergence. To validate the reconstruction error, for each CS solver method, we infer the error norm between the reconstructed data from them and the one from the convex solver at all five amounts of the dataset, which is expressed as E r r = x k x C V X 2 , where k (= 1, 2, 3) represents for the order of three solvers (ADMM, IRLS, and CD). It can be seen in Figure 9a that, for all solvers, as we increment the amount of data in of dataset, the reconstructed error is reduced. More specifically, ADMM brings about the lowest error from the recovered data at all parts which means that it performs better than its counterparts in our problem. On the other hand, Figure 9b reveals the convergence of this solver after 100 iterations. It again confirms the negative correlation between the amount of data and the objective function value in terms of finding the minimum solution as, the higher percentage of the dataset, the lower value of the objective function. Overall, CD takes more iterations than other methods to converge which leads to a longer duration to recover the data, while IRLS takes the least time in the same task.
Finally, the formed SAR images from the recovered data by adopting ADMM solvers on four parts of data including 10%, 30%, 40%, and 50% from the original one are shown in Figure 10. Here, we reveal the images of the rebar consisting of cracks of 1 mm and 3 mm when the distance z 0 = 50 mm. In Figure 10a, due to the high reconstructed error from the small amount of data, it illustrates the distortion on the shape of the rebar resulting in the undetectable crack positions on it. On the contrary, the SAR images formed by using 30% to 50% of the data (Figure 10b,d) present higher quality with less distorted segments, which can help to determine the location of the defect on the MUT. Nevertheless, the low contrast in Figure 10b, where only 30% of data is employed, still makes the defect slightly difficult to be recognized compared to Figure 10c,d, where the 3 mm crack can be shown as the same as the original in Figure 8a.

4. Conclusions

In this paper, we describe the design, development, and evaluation of the SAR imaging system adopting the mmWave signal from FMCW radar. Our scheme demonstrates its ability to detect the defect in the object by locating its position and constructing its 2-D surface shape, which can support finding and determining the failures when that object is invisible to human eyes. Specifically, we adopted the imaging scheme to the NDT application for the reinforced concrete and although it failed to detect the small cracks in its rebar, it works well in cases of larger cracks. Additionally, this proposed system can take part in various see-through-wall applications from detecting concealed weapons, structure health monitoring, and material assessment during manufacturing to finding and saving living objects under the wreckage after a disaster. To archive the 2-D imaging of the object, the RMA method was investigated along with the digital signal processing techniques applied to remove the undesired signal and improve the image quality. Besides that, the involvement of the CS algorithm in the imaging system helps to enhance the performance of the system as it significantly reduces the requirements for spatial data acquisition and the side-lobe effect from the non-uniform spacing of measurement points. Three CS-based reconstruction algorithms, namely ADMM, IRLS, and CD, are presented to reconstruct the data with spare samples. The results confirm that the scheme can form the image even with a sparsity rate of 30%, while the ADMM algorithm shows better performance. In the future, the methods to improve the azimuth resolution of the system need to be researched for identifying tiny defects by joining the AI algorithms and placing virtual antenna elements on wider non-uniform samples.

Author Contributions

Conceptualization, I.-P.H. and K.-H.K.; methodology, software, T.-H.P.; validation, T.-H.P.; formal analysis, investigation, T.-H.P.; writing—original draft preparation, T.-H.P.; writing—review and editing, I.-P.H.; visualization, T.-H.P.; supervision, I.-P.H.; project administration, I.-P.H.; funding acquisition, I.-P.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Basic Science Research Program under Grant 2020R1I1A3057142, in part by the Priority Research Center Program through the National Research Foundation under Grant 2019R1A6A1A03032988, and in part by the Underground City of the Future Program funded by the Ministry of Science and ICT.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Runkiewicz, L. Application of non-destructive testing methods to assess properties of construction materials in building diagnostics. Archit. Civ. Eng. Environ. 2009, 2, 79–86. [Google Scholar]
  2. Einav, I. Non-Destructive Testing for Plant Life Assessment; International Atomic Energy Agency: Vienna, Austria, 2005; pp. 1–61. [Google Scholar]
  3. Ekinci, S.; Aksu, M.; Bingoeldag, M.; Dogruoez, M.; Kurtcebe, T.; Yildirim, A.; Saricam, S.; Yilmaz, N. Development of Protocols for Corrosion and Deposits Evaluation in Pipes by Radiography; IAEA-TECDOC-1445; International Atomic Energy Agency: Vienna, Austria, 2006. [Google Scholar]
  4. Haibo, L.; Yongqing, W.; Meng, L.; Tongyu, Z.; Baoliang, L. Thickness Measurement Using Ultrasonic Scanning Method for Large Aerospace Thin-Walled Parts. In Proceedings of the 2019 IEEE 5th International Workshop on Metrology for AeroSpace (MetroAeroSpace), Turin, Italy, 19–21 June 2019; pp. 243–247. [Google Scholar]
  5. Zhongzhu, L.; Chunguang, X.; Bolong, M. Application of SAGE Algorithm to Estimate Thin Layer Material’s Thickness in Ultrasonic NDE. In Proceedings of the 2011 International Conference on Electronic & Mechanical Engineering and Information Technology, Harbin, China, 12–14 August 2011; pp. 2869–2872. [Google Scholar]
  6. Yao, Y.; Tung, S.-T.; Glisic, B. Crack Detection and Characterization Techniques - An Overview. Struct. Control Health Monit. 2014, 21, 1387–1413. [Google Scholar] [CrossRef]
  7. Klepka, A.; Staszewski, W.J.; Jenal, R.B.; Szwedo, M.; Iwaniec, J.; Uhl, T. Nonlinear Acoustics for Fatigue Crack Detection Experimental Investigations of Vibro-acoustic Wave Modulations. Struct. Health Monit. 2012, 11, 197–211. [Google Scholar] [CrossRef]
  8. Bustillo, J.; Achdjian, H.; Arciniegas, A.; Blanc, L. Simultaneous determination of wave velocity and thickness on overlapped signals using Forward Backward algorithm. NDT&E Int. 2017, 86, 100–105. [Google Scholar]
  9. Bertovic, M.; Gaal, M.; Müller, C.; Fahlbruch, B. Investigating human factors in manual ultrasonic testing: Testing the human factors model. Insight Non-Destr. Test. Cond. Monit. 2011, 53, 673–676. [Google Scholar] [CrossRef]
  10. Galin, N.; Worby, A.; Markus, T.; Leuschen, C.; Gogineni, P. Validation of Airborne FMCW Radar Measurements of Snow Thickness Over Sea Ice in Antarctica. IEEE Trans. Geosci. Remote Sens. 2012, 50, 3–12. [Google Scholar]
  11. Iglesias, R.; Aguasca, A.; Fabregas, X.; Mallorqui, J.J.; Monells, D.; Lpez-Martnez, C.; Pipia, L. Ground-Based Polarimetric SAR Interferometry for the Monitoring of Terrain Displacement Phenomena; Part II: Applications. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 994–1007. [Google Scholar] [CrossRef] [Green Version]
  12. Luo, Y.; Song, H.; Wang, R.; Deng, Y.; Zhao, F.; Xu, Z. Arc FMCW SAR and Applications in Ground Monitoring. IEEE Trans. Geosci. Remote Sens. 2014, 52, 5989–5998. [Google Scholar] [CrossRef]
  13. Charvat, G.L.; Kempell, L.C.; Coleman, C. A Low-Power High-Sensitivity X-Band Rail SAR Imaging System. IEEE Antennas Propag. Mag. 2008, 50, 108–115. [Google Scholar] [CrossRef]
  14. Wang, G.; Gu, C.; Inoue, T.; Li, C. A Hybrid FMCW-Interferometry Radar for Indoor Precise Positioning and Versatile Life Activity Monitoring. IEEE Trans. Microw. Theory Tech. 2014, 62, 2812–2822. [Google Scholar] [CrossRef]
  15. Boukari, B.; Moldovan, E.; Affes, S.; Wu, K.; Bosisio, R.G.; Tatu, S.O. Six-Port FMCW Collision Avoidance Radar Sensor Configurations. In Proceedings of the 2008 Canadian Conference on Electrical and Computer Engineering, Niagara Falls, ON, Canada, 4–7 May 2008; pp. 305–308. [Google Scholar]
  16. Polychronopoulos, A.; Amditis, A.; Floudas, N.; Lind, H. Integrated object and road borders tracking using 77GHz automotive radars. IEE Proc. Radar Sonar Navig. 2004, 151, 375–381. [Google Scholar] [CrossRef]
  17. Russell, M.E.; Crain, A.; Curran, A.; Campbell, R.A.; Drubin, C.A.; Miccioli, W.F. Millimeter-Wave Radar Sensor for Automotive Intelligent Cruise Control (ICC). IEEE Trans. Microw. Theory Tech. 1997, 45, 2444–2453. [Google Scholar] [CrossRef]
  18. Ann, F.; Peter, L.; Luc, T. 2-D Near-field SAR Non-destructive Testing of Rebars in a Concrete Wall. Int. J. Appl. Electromagn. Mech. 2004, 19, 333–338. [Google Scholar]
  19. Sato, M.; Takahashi, K.; Iitsuka, Y.; Koyama, C.N. Near range polarimetric SAR for non-destructive inspection of wooden buildings. In Proceedings of the IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Singapore, 1–4 September 2015; pp. 306–309. [Google Scholar]
  20. Lsaadouny, M.E.; Barowski, J.; Jebramcik, J.; Rolfes, I. Millimeter Wave SAR Imaging for the Non-Destructive Testing of 3D-printed Samples. In Proceedings of the 2019 International Conference on Electromagnetics in Advanced Applications (ICEAA), Granada, Spain, 24 October 2019; pp. 1283–1285. [Google Scholar]
  21. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  22. Candès, E.J.; Romberg, J.K.; Tao, T. Stable signal recovery from incomplete and inaccurate measurements. Commun. Pure Appl. Math. 2006, 59, 1207–1223. [Google Scholar] [CrossRef] [Green Version]
  23. Zhang, X.; Xi, X.; Pu, Y.; Zheng, C.; Wu, D. Compressive sensing for through-wall imaging with spread spectrum radar. In Proceedings of the 2016 IEEE MTT-S International Wireless Symposium (IWS), Shanghai, China, 14–16 March 2016; pp. 1–4. [Google Scholar]
  24. Kazemi, P.; Modarres-Hashemi, M.; Naghsh, M.M. Block compressive sensing in Synthetic Aperture Radar (SAR). In Proceedings of the 2016 24th Iranian Conference on Electrical Engineering (ICEE), Shiraz, Iran, 10–12 May 2016; pp. 1324–1329. [Google Scholar]
  25. Yang, C.; Zhan, Q.; Liu, H.; Ma, R. An IHS-Based Pan-Sharpening Method for Spectral Fidelity Improvement Using Ripplet Transform and Compressed Sensing. Sensors 2018, 18, 3624. [Google Scholar] [CrossRef] [Green Version]
  26. Jung, D.-H.; Kang, H.-S.; Kim, C.-K.; Park, J.; Park, S.-O. Sparse Scene Recovery for High-Resolution Automobile FMCW SAR via Scaled Compressed Sensing. IEEE Trans. Geosci. Remote Sens. 2019, 57, 10136–10146. [Google Scholar] [CrossRef]
  27. Phippen, D.; Daniel, L.; Hoare, E.; Cherniakov, M.; Gashinova, M. Compressive sensing for automotive 300GHz 3D imaging radar. In Proceedings of the 2020 IEEE Radar Conference (RadarConf20), Florence, Italy, 21–25 September 2020; pp. 1–6. [Google Scholar]
  28. Lee, S.; Jung, Y.; Lee, M.; Lee, W. Compressive Sensing-Based SAR Image Reconstruction from Sparse Radar Sensor Data Acquisition in Automotive FMCW Radar System. Sensors 2021, 21, 7283. [Google Scholar] [CrossRef]
  29. Candés, E.; Romberg, J. L1-Magic: Recovery of Sparse Signals via Convex Programming; Caltech: Pasadena, CA, USA, 2005; pp. 1–19. [Google Scholar]
  30. Yanik, M.E.; Torlak, M. Near-field MIMO-SAR millimeter-wave imaging with sparsely sampled aperture data. IEEE Access. 2019, 7, 31801–31819. [Google Scholar] [CrossRef]
  31. Kim, H.; You, S.; Jeong, B.J.; Byun, W. Azimuth Angle Resolution Improvement Technique with Neural Network. In Proceedings of the 2020 International Conference on Information and Communication Technology Convergence (ICTC), Jeju, Korea, 21–23 October 2020; pp. 1384–1387. [Google Scholar]
  32. Li, J.; Stoica, P. MIMO Radar Signal Processing, 1st ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2008. [Google Scholar]
  33. Wang, G.; Munoz-Ferreras, J.-M.; Gu, C.; Li, C.; Gomez-Garcia, R. Application of linear-frequency-modulated continuous-wave (LFMCW) radars for tracking of vital signs. IEEE Trans. Microw. Theory Tech. 2014, 62, 1387–1399. [Google Scholar] [CrossRef]
  34. Silverstein, S.D.; Zheng, Y. Near-field inverse coherent imaging problems: Solutions, simulations, and applications. In Proceedings of the Thrity-Seventh Asilomar Conference on Signals, Systems & Computers, Pacific Grove, CA, USA, 9–12 November 2003; pp. 1193–1197. [Google Scholar]
  35. Brekhovskikh, L.M.; Godin, O. Acoustics of Layered Media II: Point Sources and Bounded Beams; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  36. Yanik, M.E.; Torlak, M. Near-Field 2-D SAR Imaging by Millimeter-Wave Radar for Concealed Item Detection. In Proceedings of the 2019 IEEE Radio and Wireless Symposium (RWS), Orlando, FL, USA, 20–23 January 2019; pp. 1–4. [Google Scholar]
  37. Stephen, L.; Bruton, J.; Nathan, K. Data Driven Science & Engineering; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  38. Candes, E.J.; Romberg, J.; Tao, T. Exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 2004, 52, 489–509. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, G.; Jiao, S.; Xu, X.; Wang, L. Compressed sensing and reconstruction with Bernoulli matrices. In Proceedings of the 2010 IEEE International Conference on Information and Automation, Harbin, China, 20–23 June 2010; pp. 455–460. [Google Scholar]
  40. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers; Now Foundations and Trends: Norwell, MA, USA, 2011; pp. 1–122. [Google Scholar]
  41. Liu, H.; Song, B.; Qin, H.; Qiu, Z. An adaptive-ADMM algorithm with support and signal value detection for compressed sensing. IEEE Signal Process. Lett. 2013, 20, 315–318. [Google Scholar] [CrossRef]
  42. Li, Y.; Osher, S. Coordinate descent optimization for l 1 minimization with application to compressed sensing; a greedy algorithm. Inverse Probl. Imaging 2009, 3, 487–503. [Google Scholar] [CrossRef]
  43. mmWave Radar Sensors|TI.com. 2022. Available online: http://www.ti.com/sensors/mmwave/overview.html (accessed on 1 May 2022).
Figure 1. Block diagram of the imaging system.
Figure 1. Block diagram of the imaging system.
Sensors 22 08030 g001
Figure 2. Target measurement using mmWave radar system.
Figure 2. Target measurement using mmWave radar system.
Sensors 22 08030 g002
Figure 3. The measurement trajectory.
Figure 3. The measurement trajectory.
Sensors 22 08030 g003
Figure 4. Compressed sensing processes: (a) signal compression and (b) signal reconstruction.
Figure 4. Compressed sensing processes: (a) signal compression and (b) signal reconstruction.
Sensors 22 08030 g004
Figure 5. Measurement system (a) the TI AWR1440Boost with DCA1000 data card and (b) the 2-D motion platform.
Figure 5. Measurement system (a) the TI AWR1440Boost with DCA1000 data card and (b) the 2-D motion platform.
Sensors 22 08030 g005
Figure 6. Reinforced concrete sample structures consisting of: (a) 1 mm and 3 mm cracks (sample A) and (b) 2 mm and 4 mm cracks (sample B).
Figure 6. Reinforced concrete sample structures consisting of: (a) 1 mm and 3 mm cracks (sample A) and (b) 2 mm and 4 mm cracks (sample B).
Sensors 22 08030 g006
Figure 7. Non-destructive testing measurement setup: (a) the top view and (b) the real setup.
Figure 7. Non-destructive testing measurement setup: (a) the top view and (b) the real setup.
Sensors 22 08030 g007
Figure 8. SAR image results at distance z 0 = 50 mm of (a) sample A, (b) sample B, and at z 0 = 100 mm of (c) sample A and (d) sample B.
Figure 8. SAR image results at distance z 0 = 50 mm of (a) sample A, (b) sample B, and at z 0 = 100 mm of (c) sample A and (d) sample B.
Sensors 22 08030 g008
Figure 9. (a) The reconstruction error and (b) the convergence of three CS solvers.
Figure 9. (a) The reconstruction error and (b) the convergence of three CS solvers.
Sensors 22 08030 g009
Figure 10. SAR images from the data recovered by the ADMM method employing (a) 10% data, (b) 30% data, (c) 40% data, and (d) 50% data.
Figure 10. SAR images from the data recovered by the ADMM method employing (a) 10% data, (b) 30% data, (c) 40% data, and (d) 50% data.
Sensors 22 08030 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pham, T.-H.; Kim, K.-H.; Hong, I.-P. A Study on Millimeter Wave SAR Imaging for Non-Destructive Testing of Rebar in Reinforced Concrete. Sensors 2022, 22, 8030. https://doi.org/10.3390/s22208030

AMA Style

Pham T-H, Kim K-H, Hong I-P. A Study on Millimeter Wave SAR Imaging for Non-Destructive Testing of Rebar in Reinforced Concrete. Sensors. 2022; 22(20):8030. https://doi.org/10.3390/s22208030

Chicago/Turabian Style

Pham, The-Hien, Kil-Hee Kim, and Ic-Pyo Hong. 2022. "A Study on Millimeter Wave SAR Imaging for Non-Destructive Testing of Rebar in Reinforced Concrete" Sensors 22, no. 20: 8030. https://doi.org/10.3390/s22208030

APA Style

Pham, T. -H., Kim, K. -H., & Hong, I. -P. (2022). A Study on Millimeter Wave SAR Imaging for Non-Destructive Testing of Rebar in Reinforced Concrete. Sensors, 22(20), 8030. https://doi.org/10.3390/s22208030

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