Next Article in Journal
Numerical Investigation of Heat Transfer on Unsteady Hiemenz Cu-Water and Ag-Water Nanofluid Flow over a Porous Wedge Due to Solar Radiation
Next Article in Special Issue
Multiscale Compact Modelling of UTC-Photodiodes Enabling Monolithic Terahertz Communication Systems Design
Previous Article in Journal
GI Guidelines for the Metropolitan City of Cagliari (Italy): A Method for Implementing Green Areas
Previous Article in Special Issue
Flexible Ultra-Wide Electro-Optic Frequency Combs for a High-Capacity Tunable 5G+ Millimeter-Wave Frequency Synthesizer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A High-Frequency Vibration Error Compensation Method for Terahertz SAR Imaging Based on Short-Time Fourier Transform

1
Terahertz Technology Innovation Research Institute, Terahertz Spectrum and Imaging Technology Cooperative Innovation Center, Shanghai Key Lab of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China
2
School of Optoelectronic Information and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
3
Shanghai Radio Equipment Research Institute, Shanghai 201109, China
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2021 / Revised: 14 November 2021 / Accepted: 15 November 2021 / Published: 17 November 2021
(This article belongs to the Special Issue Applications of Millimeter-Wave and Terahertz Technologies)

Abstract

:
High-frequency vibration error of a moving radar platform easily introduces a non-negligible phase of periodic modulation in radar echoes and greatly degrades terahertz synthetic aperture radar (THz-SAR) image quality. For solving the problem of THz-SAR image-quality degradation, the paper proposes a multi-component high-frequency vibration error estimation and compensation approach based on the short-time Fourier transform (STFT). To improve the robustness of the method against noise effects, STFT is used to extract the instantaneous frequency (IF) of a high-frequency vibration error signal, and the vibration parameters are coarsely obtained by the least square (LS) method. To reduce the influence of the STFT window widths, a method based on the maximum likelihood function (MLF) is developed for determining the optimal window width by a one-dimensional search of the window widths. In the case of high noise, many IF estimation values seriously deviate from the true ones. To avoid the singular values of IF estimation in the LS regression, the random sample consensus (RANSAC) is introduced to improve estimation accuracy. Then, performing the STFT with the optimal window width, the accurate vibration parameters are estimated by LS regression, where the singular values of IF estimation are excluded. Finally, the vibration error is reconstructed to compensate for the non-negligible phase of the platform-induced periodic modulation. The simulation results prove that the error compensation method can meet THz-SAR imaging requirements, even at a low signal-to-noise ratio (SNR).

1. Introduction

Terahertz (THz) waves are the electromagnetic waves between 0.1 THz and 10 THz; they are in the gap between microwave and infrared. In terms of their energy, they are between electrons and photons [1,2]. Therefore, THz waves have many unique advantages that microwaves and infrared do not have. For example, THz waves can generate a higher carrier frequency and larger absolute bandwidth than microwaves. At the same time, THz radiation can penetrate materials such as ceramics, fats, fabrics, plastics, etc., with very little attenuation. In addition, when it penetrates the human body, THz radiation is harmless due to its lower photon energy. These unique characteristics allow THz imaging to be used in many fields [3,4,5,6,7].
The concept of a synthetic aperture imaging technique in the microwave band has been introduced in the THz band [8,9]. Compared to microwave synthetic aperture radar (SAR), THz-SAR can achieve higher two-dimensional image resolution. Additionally, with its shorter synthetic aperture time, THz-SAR imaging has a higher frame rate [10,11]. Theoretically, the SAR principle requires uniform linear movement of the antenna phase center (APC). In practice, due to the influence of atmospheric turbulence and artificial control, APC always deviates from the ideal motion state, resulting in motion error. According to the variations of the motion error in the whole synthetic aperture, it can be divided into low-frequency motion error and high-frequency vibration error. The low-frequency motion error ranges from the centimeter level to the meter level, but now there are many compensation methods to eliminate its influence on SAR imaging. Meanwhile, the vibration amplitude of high-frequency vibration error is very small, usually at the millimeter level or even the sub-millimeter level, such that its influence can be ignored in microwave SAR imaging. However, since the wavelengths of THz waves are relatively short and comparable to the vibration amplitude of high-frequency vibration error, its influence has to be considered in THz-SAR imaging.
Currently, there exist two kinds of high-frequency vibration-error compensation methods. Similar to motion-error compensation in microwave SAR, phase errors caused by high-frequency vibration error is estimate directly through extended phase gradient autofocus (PGA) [12,13], where the specific vibration information is still ambiguous after phase-error compensation. Therefore, some researchers have focused on another class of methods based on parameter estimation, whose main idea is to estimate the vibration-signal parameters. Sinusoidal frequency modulation Fourier transform (SFMFT) and wavelet multi-resolution analysis are directly utilized to operate on the phase term of the vibration signal through phase unwrapping in [14,15,16]. However, these methods will fail when the phase ambiguity is serious.
In [17], the discrete fractional Fourier transform (DFrFT) is used to reconstruct the vibration signal with the estimated instantaneous vibration accelerations and vibration frequency. In [18], the vibration parameters are estimated based on the wavelet transform. The adaptive chirplet decomposition and the short-time Fourier transform (STFT) are used in [19] and [20], respectively, and an estimation method considering the variation of vibration amplitude is proposed in [21]. However, the above methods may fail in dealing with multiple frequency component vibration.
The discrete sinusoidal frequency modulation transforms (DSFMT) and the optimization method are combined to estimate the multi-component vibration parameters in [22,23]. However, it inevitably results in large estimation errors because of the convergence to a local minimum. In [24], an improved DFrFT-based estimation method is proposed to estimate multi-component vibration parameters. However, there exists a great correlation between the method accuracy and the window width.
In reality, all the above estimation methods must first extract the vibration signal of the strong point target in the coarsely focused image. The extracted vibration signal will be disturbed by the surrounding clutter signal, and the signal-to-clutter ratio may be fairly low. Although it is biased in estimating instantaneous frequency (IF), STFT is comparatively robust against noise effects [25]. Thus, a novel parameter estimation method of multi-component high-frequency vibration error based on the STFT is presented in this paper. After IF estimates using STFT, least square (LS) regression is performed on the IF estimates to estimate the high-frequency vibration parameters. The estimated parameters are insensitive to the noise effects. However, the method accuracy is influenced by the window width. To overcome the problem, the parameters are estimated for various window widths in the STFT. Then, a method based on the maximum likelihood function (MLF) is developed for determining the optimal window width by a one-dimensional search of these window widths. Here, the multi-dimensional search for a set of parameters is avoided, thus reducing the computing time. In addition, the method accuracy is also influenced by the accuracy of the IF-estimate samples in the LS regression. If the noise is significant, the estimation error is aggravated by the IF estimates that deviate significantly from the true value. To avoid singular IF estimates in the LS regression, a random sample consensus (RANSAC) [26,27] is introduced to randomly select the IF-estimate samples for multiple times. The previous MLF can also be used as a measure of the selection results. Finally, to show the superiority of the proposed method, it is compared with the original STFT-based method.
The rest of this paper is organized as follows. Section 2 establishes the high-frequency vibration error model and the THz-SAR echo-signal model. The proposed estimation and compensation method is described in detail in Section 3. Section 4 gives the simulation results, and conclusions are drawn in Section 5.

2. THz-SAR Echoes Signal Model with Vibration Error

2.1. High-Frequency Vibration Error Model of Platform

In SAR imaging, besides the translational motion error, the vibration error of the motion platform is also inevitable because of the effects of air flow and engine vibration. Compared with traditional translational motion error, the motion platform’s vibration is of the higher vibration frequency that usually meets the conditions.
| f · T s | 1
where T s is the synthetic aperture time of SAR.
So far, the vibration of motion platform is usually modeled as a simple harmonic vibration, which is expressed as
Δ R ( t ) = a sin ( 2 π f t + ϕ )
where t represents the azimuth time, a denotes the vibration amplitude, ϕ represents the initial phase, and f is vibration frequency.
However, in reality, the high-frequency vibration of the motion platform is complex and modeled as a superposition of multiple harmonic vibration. Therefore, to characterize the motion platform vibration more accurately in this paper, we model it with the form of
Δ R ( t ) = i = 1 I a i sin ( 2 π f i t + ϕ i )
where I represents the number of vibration components, a i , f i , and ϕ i denote the vibration amplitude, vibration frequency and initial phase of the i’th vibration component.

2.2. THz-SAR Echo Signal Model

Figure 1 gives the THz-SAR imaging geometric model when there is high-frequency vibration error. The radar operates in right side-looking strip-map mode. Ideally, the radar platform moves horizontally in the X direction at a reference altitude H and a constant velocity v .
The coordinate of the ideal position of the APC is ( v t ,   0 , H ) . Take a point scatterer, P, with the coordinate ( v t 0 ,   y 0 , 0 ) as an example, where t 0 is the zero doppler time and y 0 is the Y-axis position. The ideal instantaneous slant range r ( t ) from the APC to the point is
r ( t ) = r 0 2 + v 2 ( t t 0 ) 2 = r 0 + v 2 ( t t 0 ) 2 2 r 0
where r 0 = y 0 2 + H 2 is the nearest slant range of point scatterer.
However, in practical applications, the actual trajectory of the APC deviates from the ideal motion state because of the vibration error and is indicated by the curve line. Thus, the instantaneous slant range is affected by high-frequency vibration error. The actual instantaneous slant range R ( t ) is
R ( t ) = r ( t ) + Δ R ( t )
Assume that the linear frequency modulation (LFM) pulse signal emitted by THz-SAR is
s ( τ ) = rect ( τ T ) exp ( j 2 π f 0 τ + j π γ τ 2 )
where rect ( · ) denotes the rectangular window function, τ represents the fast time, T denotes the signal pulse width, f 0 represents the carrier frequency, and γ denotes the chirp rate.
In THz-SAR imaging, the radar system usually uses the dechirp heterodyne receiver [28]. After compensating the residual video phase, the recorded baseband echo signal is
s ( τ , t ) = σ · rect [ τ 2 R ( t ) / c T ] exp { j 4 π f 0 c R ( t ) + j 4 π γ R ( t ) c τ }
where c is the speed of light.
Then, through performing the fast Fourier transform (FFT) upon τ and a systematic range cell-migration correction (RCMC), the range-compressed signal is
s ( f , t ) = δ ( f γ 2 r 0 c ) exp { j 4 π λ r } exp { j 4 π λ Δ R }
where δ ( · ) denotes the impulse function and λ represents the wavelength.
Substituting (3) and (4) into (8), we can obtain
s ( f , t ) = δ ( f γ 2 r 0 c ) exp { j 4 π λ ( r 0 + v 2 ( t t 0 ) 2 2 r 0 ) } exp { j 4 π λ i = 1 I a i sin ( 2 π f i t + ϕ i ) }
Then, after performing the dechirp operation and the residual video phase compensation in azimuth, (9) can be rewritten as
s ( f , t ) = δ ( f γ 2 r 0 c ) exp { j 4 π λ ( r 0 + v 2 t 0 r 0 t ) } exp { j 4 π λ i = 1 I a i sin ( 2 π f i t + ϕ i ) }
Finally, the azimuth compression in the Doppler domain can be accomplished by performing the FFT upon t in (10). However, the high-frequency vibration error brings the additional period-modulated phase into the ideal echo, as shown in the last term of (10). The phase error results in the paired echo in the THz-SAR images.

3. Compensation Method

3.1. STFT-Based Rough Estimation

First, we select the dominant point target from the coarsely focused image. The signal of interest (SoI) is obtained through performing inverse FFT to convert the coarsely focused image to the range-compressed phase history domain. To eliminate the linear phase in SoI, similar to PGA [29], the dominant point target is circularly shifted to the image center before performing inverse FFT. Thus, after accomplishing the above operations, we obtain the SoI, which is expressed as
s ( t ) = exp { j 4 π λ i = 1 I a i sin ( 2 π f i t + ϕ i ) }
The signal is sampled with the sampling time Δ t = 1 / f P R F , where f P R F is the pulse repetition frequency. The discretized SoI signal s ( n ) is represented as
s ( n ) = s ( n · Δ t ) ,   n [ N / 2 , N / 2 ]
where N represents the sample point number.
On the one hand, the IF can be estimated by using the STFT. That is to say
f ^ h ( n ) = arg   max f | STFT h ( n , f ) |
where STFT h ( n , f ) is the STFT of the signal with window width h .
On the other hand, the IF of the signal can be acquired by taking the derivative of the phase in (11):
f ( t ) = 2 λ i = 1 I 2 π f i a i cos ( 2 π f i t + ϕ i )
After discretization, the IF in (14) can be rewritten as
f ( n ) = i = 1 I [ C 0 , i cos ( 2 π f i n Δ t ) + S 0 , i sin ( 2 π f i n Δ t ) ]
where
{ C 0 , i = 2 λ 2 π f i a i cos φ i S 0 , i = 2 λ 2 π f i a i sin φ i
After estimating the IF f ^ h ( n ) through (13), the high-frequency vibration parameters are obtained by solving the following linear equations, namely:
f ^ = M x
where
{ f ^ = [ f ^ ( Q / 2 ) ,   f ^ ( Q / 2 + 1 ) , , f ^ ( Q / 2 1 ) , f ^ ( Q / 2 ) ] T x = [ C 0 , 1 , S 0 , 1 , C 0 , 2 , S 0 , 2 , C 0 , I , S 0 , I ]   T M = [ a 0 , 1 , b 0 , 1 , a 0 , 2 , b 0 , 2 , a 0 , I , b 0 , I ] a 0 , i = [ c o s ( 2 π f i ( Q / 2 ) Δ t ) , c o s ( 2 π f i ( Q / 2 ) Δ t ) ] T b 0 , i = [ s i n ( 2 π f i ( Q / 2 ) Δ t ) , s i n ( 2 π f i ( Q / 2 ) Δ t ) ] T
and Q + 1 is the number of azimuth instants.
According to (14), the frequencies of the IF are the same as those of the vibration. Thus, through performing the spectral analysis on the estimated IF, the vibration frequencies are obtained. Based on the estimated vibration frequency, the LS solution of (17) is
x ^ = ( M T M ) 1 M T   f ^
Thus, the other high-frequency vibration parameters are obtained by
a i = x ^ 2 ( 2 i 1 ) + x ^ 2 ( 2 i ) 2 λ 2 π f ^ i
ϕ i = arctan ( x ^ ( 2 i ) x ^ ( 2 i 1 ) )
Although the obtained parameters are insensitive to noise, the method’s accuracy is affected by the window width. It is due to that the STFT has the limitation that it cannot satisfy the optimal frequency resolution and time resolution at the same time [30,31]. The window width h of STFT will bring errors to the estimation results. Meanwhile, the accuracy of the IF estimation samples in the LS regression also influences the method’s accuracy, especially in weak targets surrounded by clutter. Owing to these two factors, there may be a large deviation between the estimated values and the actual values, which cannot meet the requirement of THz-SAR vibration error compensation. Thus, THz-SAR images are still influenced by residual phase errors.

3.2. Improved STFT-Based Fine Estimation

Owing to the window width, h , and the singular values of IF estimation, the estimated vibration parameters values may have a large deviation from the true ones. Here, MLF and RANSAC are introduced to the STFT-based estimation, which can improve the estimation accuracy.
Firstly, the IF is estimated by the STFT for various window widths. The set of various window widths in the STFT is given as
H = { h | h m i n h h m a x }
where h m i n and h m a x are the minimum and maximum window width, which can be set according to the size of the sampling data. For each window width h H , vibration parameters are estimated through the above STFT-based estimation method, and a set of vibration parameters { a ^ i , f ^ i , ϕ ^ i } can be obtained. Then, to get the optimal window width, one-dimensional search is performed on the window width of STFT. In this paper, the MLF is adopted as the evaluation function.
G ( h ) = | n s ( n ) · exp ( j 4 π λ i I a ^ i sin ( 2 π f ^ i n Δ t + ϕ ^ i ) ) |
When the value of MLF is maximal, the corresponding window width is the optimal window, and the vibration parameters under the optimal window are to be, themselves, optimal, obtained through
{ a ^ i , f ^ i , ϕ ^ i } = arg   max h G ( h )
In this way, the estimation error caused by window width of STFT is greatly reduced and the parameter estimation accuracy is improved. What is more, the 3I-dimensional search over the parameter space is replaced by a one-dimensional search of window width. Compared with the STFT-based estimation method, the determination of the optimal window only adds a small amount of computation.
Although the refinement of the optimal window width determination can improve the estimation accuracy, the singular values of IF estimation still have a negative influence on the method performance. As the noise level increases, the percentage of the singular values of IF estimation rapidly increases and, with it, the greater the influence on the estimation accuracy. Thus, on the basis of determining optimal window width, RANSAC is used to reduce the estimation error caused by the singular values of IF estimation. A summary of the proposed STFT-based fine estimation method is described next.
Step 1. Initialize a set H of window width, and set the current value of MLF to G = 0 .
Step 2. For window width h j H , the STFT of the signal s ( n ) is calculated and IF estimates are obtained through (13).
Step 3. Vibration parameters { a ^ i , f ^ i , ϕ ^ i } are estimated by the LS method.
Step 4. On the basis of the estimated vibration parameters, use (23) to calculate the MLF G ( h ) . If G ( h ) > G , set G = G ( h ) , h = h j and record the corresponding estimated parameters values.
Step 5. Determine whether to traverse the window width completely. If not, repeat Step 2 to Step 4; otherwise, the current window width is optimal.
Step 6. Set the iteration number M, set the parameter values and the criterion function under the optimal window width as the initialization vibration parameters and the initialization G , separately.
Step 7. For the optimal window width h , use the STFT to estimate the IF through (13).
Step 8. When the vibration parameters are estimated by LS for the m’th time, the 2I IF estimation values are randomly selected from the Q + 1 IF values estimated in Step 7. The corresponding discrete time of these 2I IF values is n m ( q ) , q = 1 , 2 , , 2 I and n m ( q ) < n m ( q + 1 ) < Q . These selected IF estimation values form the set f ^ ( m ) .
f ^ ( m ) = { f ^ ( n m ( q ) · Δ t ) , q = 1 , 2 , , 2 I }
Step 9. LS method is used to estimate the vibration parameters, namely
x ^ ( m ) = ( K T K ) 1 K T f ^ ( m )
where
{ K = [ c 0 , 1 , d 0 , 1 , c 0 , 2 , d 0 , 2 , c 0 , I , d 0 , I ] c 0 , i = [ cos ( 2 π f i n m ( 1 ) Δ t ) , cos ( 2 π f i n m ( 2 I ) Δ t ) ] T d 0 , i = [ sin ( 2 π f i n m ( 1 ) Δ t ) , sin ( 2 π f i n m ( 2 I ) Δ t ) ] T
Then, use (20) and (21) to calculate the parameters { a ^ i , m , φ ^ i , m } .
Step 10. On the basis of the estimated vibration parameters, use (23) to calculate the MLF G ( m ) . If G ( m ) > G , set G = G ( m ) and record the corresponding estimated parameters values.
Step 11. If m < M , repeat Step 8 to Step 10. Otherwise, the current estimation values are the final estimation results.
After the above improved STFT-based fine estimation, the effects of the STFT window width and the singular values of IF estimation have been reduced as much as possible. The estimate accuracy of vibration parameters has been greatly improved.

3.3. Vibration Compensation Method

In THz-SAR imaging, the high-frequency vibration error leads to the paired echo and reduces the image quality when the imaging algorithms are directly performed on the echo data. Here, a multi-component high-frequency vibration error compensation method based on STFT is proposed for solving the above problem.
First, the dominant point target is selected from the coarsely focused image and the SoI is extracted for the vibration-parameters estimation. Then, the accurate vibration error parameters are estimated by the improved STFT-based fine estimation method. Finally, on the basis of the optimal estimated parameters values, the vibration error is compensated, and the THz-SAR focused image is obtained. Therefore, the flow chart of the proposed vibration error compensation method is given in Figure 2.

4. Simulation Analysis

In this section, the STFT-based rough estimation, the STFT-based rough estimation with the optimal window width and the improved STFT-based fine estimation are firstly used to estimate the vibration parameters. In a different SNR, the estimation performances of these three methods are compared. Then, through using these methods to estimate and compensate for vibration error in THz-SAR imaging, the effectiveness of the proposed method is further verified.

4.1. Estimation Performance Comparison

In this sub-section, simulations are given to evaluate the performance of the proposed parameter estimation method. Firstly, a single-component high-frequency vibration signal is considered and takes the form of (11). Signal parameters are set as a 1 = 2.00   mm , f 1 = 10   Hz , ϕ 1 = π / 3 ; and, in the noisy scenario, the signal-to-noise ratio (SNR) is defined as SNR = 1 / σ 2 , where σ 2 is the variance of white Gaussian noise.
The SNR of the high-frequency vibration signal changes by 1-dB increments between −2 dB and 6 dB. For each SNR value, the above parameters estimation methods are performed to obtain the signal parameters, and one hundred simulations are performed. The normalized root–mean–square error (NRMSE), as shown in [32], is used to evaluate the estimation performance.
NRMSE = 1 N i = 1 N ( x ^ i x ) 2 | x |
where N denotes the number of simulations, x ^ i is the estimated value for the i’th simulation, and x is the true value.
For different parameters estimation methods, Figure 3 gives the measured NRMSEs of the high-frequency vibration-signal parameters within different SNR circumstances. Here, in the STFT-based rough estimation method, we chose 20 and 30 as the window width. In the latter two methods, the optimal window width,   h o p , is 26, which is determined by maximizing MLF in (23). The STFT-based rough estimation method produces different estimation accuracies for different window widths. It can be obviously seen from Figure 3 that the improved STFT-based fine estimation has the best performance in rather low SNR and the performance difference decreases as SNR increases. When SNR is higher than 2 dB, all three methods can effectively estimate the SFM signal parameters.
The high-frequency vibration signal in the form of (11) and I = 2 is considered. Besides the first vibration component mentioned above, the second vibration parameters are set as a 2 = 0.3   mm ,   f 2 = 20.0   Hz , and φ 2 = π / 6 . In addition, the signal is also corrupted with Gaussian white noise. The SNR of the high-frequency vibration signal changes in 1-dB increments between −2 dB and 6 dB. For each SNR value, the above parameters estimation methods were performed to obtain the signal parameters, and one hundred simulations were performed.
Regarding different parameter-estimation methods, Figure 4 gives the measured NRMSEs of the high-frequency vibration signal parameters in different SNR circumstances. Here, in the STFT-based rough estimation method, we also chose 20 and 30 as window widths. In the latter two methods, the optimal window width,   h o p , is 26, which is determined by maximizing MLF in (23). Different window widths make in the STFT-based rough estimation method produce different estimation accuracies. It can be seen from the Figure 4 that the improved STFT-based fine estimation had the best performance, except for the initial phase, ϕ 2 . The proposed method is based on MLF to estimate the optimal vibration parameters. Since the amplitude of the second vibration error is relatively small, the MLF may be insensitive to the initial phase ϕ 2 of the second vibration error. Therefore, when the optimal MLF value is obtained, there is still a great difference between the initial phase and the true value. In any case, even in rather low SNR, the proposed estimation method has better estimation accuracy and stability.
To sum up, from the perspective of parameter estimation, the proposed improved STFT-based fine estimation method showed better performance than that of the rough estimation method for both single-component and multi-component high-frequency vibration signals.

4.2. Single-Component Vibration Error Compensation

In this sub-section, simulations are given to evaluate the compensation performance of the proposed vibration error-compensation method for THz-SAR imaging. The THz-SAR system parameters are shown in Table 1. When the echo signal of multi-point targets is generated, single-component vibration error is added to it. To estimate and compensate the vibration error, the dominant scatter point is firstly selected according to the energies of range bins from the unfocused image. Then, after the target is moved to the image center by a circular shift, the SoI in range-compressed phase history domain is achieved through inverse FFT. To verify the outstanding performance of the proposed vibration error compensation method, the results obtained by the STFT-based rough estimation method and the STFT-based rough estimation method with the optimal window width are also given here, respectively.
First, the IF of the SoI is estimated by the STFT and the vibration frequency is obtained from the amplitude spectrum of the IF. On the basis of the estimated vibration frequency, we can obtain the vibration amplitude and initial phase. Table 2 lists the estimated parameters and MLF values. Figure 5 shows the true vibration error and the estimated vibration error. We can see from the Table 2 that the estimated values by the STFT with h = 20 have a large deviation from the true values. Hence, the red line in Figure 5 shows that the estimated vibration displacement also deviates from the true one, especially at extreme points of vibration. Accordingly, the MLF value is much smaller than the true MLF value, also indicating that the performance of the method is poor. Although the performance of the STFT with the optimal window width h o p = 26 is improved, the estimates result still diverges from the true values, due to the singular IF estimates. Fortunately, however, Table 2 shows that the parameter-estimate values using the improved STFT-based fine estimation method are consistent with the true ones. It means that the improved method greatly suppresses both the effects of window width and the singular values of IF estimation. It is also verified by the MLF value in Table 2 and the estimated vibration displacement in the green line of Figure 5.
The vibration error is compensated for in the imaging process, and the performance of different estimation methods are then analyzed and compared. Figure 6a shows the imaging results of the RD algorithm without considering the vibration error. It is obvious that the image quality is seriously degraded due to the paired echoes caused by uncompensated vibration error. After the vibration error is estimated and compensated based on the above methods, Figure 6b–d gives the imaging results. The phase difference between the true vibration phase and the estimated vibration phase is given in Figure 7. We can see from Figure 6a,b that the image quality, after vibration error compensation, has improved, but it is still affected by paired echoes. This is caused by the residual period-modulated phase greater than π / 4 , which is shown in the red line of Figure 7. Meanwhile, the green line of Figure 7 shows that the residual phase is no greater than π / 4 , but is still periodically modulated. Therefore, the imaging quality of Figure 6c is further improved, but its side lobes are still relatively high. Fortunately, the remainder phase error is smaller than π / 4 and no longer periodically modulated after compensation by the proposed method. Therefore, the paired echoes are perfectly suppressed, as shown in Figure 6d.
Finally, we performed a quantitative analysis to assess image quality. The center point target of Figure 6b–d is intercepted using a rectangular window. Figure 8 shows the resulting azimuth-dimensional amplitude profiles. Table 3 lists the impulse response width (IRW), peak sidelobe ratio (PSLR) and integral sidelobe ratios (ISLR) for all vibration compensation methods. It is clear from Figure 8 and Table 3 that the STFT-based rough estimation method with h = 20 is too poor to improve image quality. Correspondingly, after vibration error compensation by the STFT with h o p = 26 , the image quality was improved, but its PSLR and ISLR were still high, as shown in red line of Figure 8 and Table 3. The great thing is that there are no longer paired echoes in the SAR images and the imaging quality indices are consistent with the theoretical ones after vibration error compensation by the STFT-based fine estimation method, which validate that the improved method can solve the problem of the former two methods and has the best performance.

4.3. Multi-Component Vibration Error Compensation

In this sub-section, the simulation results of multi-component vibration error compensation are provided to assess the performance of the proposed vibration error compensation method for THz-SAR imaging. Similarly, when the echo signal of multi-point targets is generated, the two components high-frequency vibration error is added. The SoI containing vibration error is extracted firstly according to the above steps. Similarly, the results obtained by the rough estimation method based on the STFT-based rough estimation method and the STFT-based rough estimation method with the optimal window width are also given here.
Using the STFT, we can obtain the IF of the SoI, and the vibration frequency is obtained from the amplitude spectrum of the IF. On the basis of the estimated vibration frequency, the vibration amplitude and initial phase are obtained through the LS regression. Table 4 lists all the estimated parameters and MLF values. Figure 9 shows the true vibration error and the estimated vibration error. It can be seen from the Table 4 that the estimated values by the STFT with h = 20 have a large deviation from the true ones. Hence, we can see from the red line of Figure 9 that the estimated vibration displacement also deviates from the true ones, especially at extreme point of vibration. Accordingly, the MLF value is much smaller than the true MLF value, implying that the performance of the STFT-based rough estimation method is poor. Although the performance of the STFT with the optimal window width h o p = 26 is improved, the estimates still diverge from the true values due to the singular IF estimates. However, of value is that the estimated values of parameters by the proposed improved STFT-based method are consistent with the true ones, as shown in Table 4. It implies that the improved method greatly suppresses both the effects of window width and the singular values of IF estimation, which is also verified by the MLF value in Table 4 and the estimated vibration displacement in the green line of Figure 9.
The vibration error is compensated for in the THz-SAR imaging process, and the performance of different estimation methods are then analyzed and compared. Figure 10a shows the focusing results of the RD algorithm. Due to vibration error, the SAR image are completely defocused. Figure 10b–d shows the imaging results after the vibration error is compensated for based on the above methods, respectively. Figure 11 also gives the phase difference between the true vibration phase and the estimated vibration phase. We can see from Figure 10a,b that the image quality after vibration error compensation is improved, but it is still affected by paired echoes. This is caused by the residual period-modulated phase error greater than π / 4 , which is shown in the red line of Figure 11. Meanwhile, the green line of Figure 11 shows that the residual phase is no greater than π / 4 , but is still periodically modulated. Therefore, the imaging quality of Figure 10c is further improved, but its side lobes are still relatively high. Fortunately, the residual phase error is smaller than π / 4 and no longer periodically modulated after vibration error compensation by the proposed improved method. Therefore, the SAR image is fully focused, which is shown in Figure 10d.
Finally, we will make a quantitative analysis to assess the image quality. The center point target of Figure 10b–d is intercepted using a rectangular window. Figure 12 gives the resulting azimuth-dimensional amplitude profiles. Table 5 lists the imaging quality indices for all vibration compensation methods. We can see from Figure 12 and Table 5 that the STFT-based rough estimation method with h = 20 is too poor to improve the image quality. Correspondingly, after vibration error compensation by the STFT with h o p = 26 , the image quality is improved, but its PSLR and ISLR are still high, as shown by red line of Figure 12 and Table 5. The value of the proposed method is demonstrated in that there are no longer paired echoes in the SAR images and the imaging quality indices are consistent with the theoretical ones after vibration error compensation by the STFT-based fine estimation method, validating that the improved method can solve the problem of the former two methods and has the best performance.

5. Conclusions

To compensate for the non-negligible phase of periodic modulation induced by high-frequency vibration error of a radar platform, a novel vibration error compensation method based on STFT with the combination of the optimal window width and RANSAC was proposed in this paper. The method not only reduced the influence of the STFT window widths but also avoided the singular values of IF estimation in the LS regression. Therefore, regardless of whether the SNR is high or low, the proposed method can fulfill the focusing requirements of the image. Finally, our simulation results have verified its validity. However, at present, the low power of terahertz radiation sources means that such a radar must operate at a close range, and there are few THz-SAR systems installed on aircraft for the time being. Therefore, we will use the measured SAR data to further test the effectiveness of the algorithm in the future. Meanwhile, all the existing high-frequency vibration error estimation methods ignore the effect of low-frequency motion errors. Thus, research on the synchronous estimation of low-frequency motion error and high-frequency vibration error in THZ-SAR imaging is the key to work advancing this technology in the future.

Author Contributions

Conceptualization, Y.L.; Funding acquisition, X.D. and Y.Z.; Methodology, Y.L.; Project administration, X.D. and Y.Z.; Resources, Q.Z.; Software, Q.W.; Validation, J.J. and Q.Z.; Writing—original draft, Q.W. and J.J.; Writing—review and editing, Y.L., Q.Z. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Natural Science Foundation of Shanghai (Grant No. 21ZR1444300), in part by the National key R and D Project of China (Grant No. 2018YFF01013003), in part by the National Natural Science Foundation of China (Grant No. 61722111, 61731020, 42005110), the 111 Project (Grant No. D18014), the International Joint Lab Project of Science and Technology Commission Shanghai Municipality (Grant No. 17590750300), the Key Project of Science and Technology Commission Shanghai Municipality (Grant No. YDZX20193100004960).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, P.H. Terahertz technology. IEEE Trans. Microw. Theory Tech. 2002, 50, 910–928. [Google Scholar] [CrossRef]
  2. Redo-Sanchez, A.; Zhang, X. Terahertz science and technology trends. IEEE J. Sel. Top. Quantum Electron. 2008, 14, 260–269. [Google Scholar] [CrossRef]
  3. Chan, W.; Deibel, J.; Mittleman, D.M. Imaging with terahertz radiation. Rep. Progr. Phys. 2007, 70, 1325–1379. [Google Scholar] [CrossRef]
  4. Balacey, H.; Recur, B.; Perraud, J.; Sleiman, J.B.; Guillet, J.; Mounaix, P. Advanced processing sequence for 3-D THz imaging. IEEE Trans. Terahertz Sci. Technol. 2016, 6, 191–198. [Google Scholar] [CrossRef]
  5. Gao, H.; Li, C.; Zheng, S.; Wu, S.; Fang, G. Implementation of the phase shift migration in MIMO-side looking imaging at terahertz band. IEEE Sens. J. 2019, 19, 9384–9393. [Google Scholar] [CrossRef]
  6. Yang, Q.; Zeng, Y.; Zhang, Y.; Wang, H.; Deng, B.; Qin, Y. Envelope correction of micro-motion targets based on multi-layer perceptron during THz-ISAR sensing. IEEE Access 2019, 7, 183596–183603. [Google Scholar] [CrossRef]
  7. Gao, H.; Li, C.; Wu, S.; Geng, H.; Zheng, S.; Qu, X.; Fang, G. Study of the extended phase shift migration for three-dimensional MIMO-SAR imaging in terahertz band. IEEE Access 2020, 8, 24773–24783. [Google Scholar] [CrossRef]
  8. Ding, J.; Kahl, M.; Loffeld, O.; Bolívar, P.H. THz 3-D image formation using SAR techniques: Simulation, processing and experimental results. IEEE Trans. Terahertz Sci. Technol. 2013, 3, 606–616. [Google Scholar] [CrossRef]
  9. Zhang, B.; Pi, Y.; Li, J. Terahertz imaging radar with inverse aperture synthesis techniques: System structure, signal processing, and experiment results. IEEE Sens. J. 2015, 15, 290–299. [Google Scholar] [CrossRef]
  10. Kim, S.; Fan, R.; Dominski, F. ViSAR: A 235 GHz radar for airborne applications. In Proceedings of the IEEE Radar Conference (RadarConf18), Oklahoma City, OK, USA, 23–27 April 2018; pp. 1549–1554. [Google Scholar]
  11. Zuo, F.; Min, R.; Pi, Y.; Li, J.; Hu, R. Improved method of video synthetic aperture radar imaging algorithm. IEEE Geosci. Remote Sens. Lett. 2019, 16, 897–901. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Sun, J.; Lei, P.; Wang, H. High-frequency vibration compensation of helicopter-borne THz-SAR. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 1460–1466. [Google Scholar] [CrossRef]
  13. Liang, Y.; Li, G.; Zhang, G.; Xiang, C.; Wu, J. A nonparametric paired echo suppression method for helicopter-borne SAR imaging. IEEE Geosci. Remote Sens. Lett. 2020, 17, 2080–2084. [Google Scholar] [CrossRef]
  14. Peng, B.; Wei, X.; Deng, B.; Chen, H.; Liu, Z.; Li, X. A sinusoidal frequency modulation Fourier transform for radar-based vehicle vibration estimation. IEEE Trans. Instrum. Meas. 2014, 63, 2188–2199. [Google Scholar] [CrossRef]
  15. Shi, S.; Li, C.; Hu, J.; Zhang, X.; Fang, G. A high frequency vibration compensation approach for terahertz SAR based on sinusoidal frequency modulation Fourier transform. IEEE Sens. J. 2021, 21, 10796–10803. [Google Scholar] [CrossRef]
  16. Xia, H.; Chen, Q.; Li, Y.; Fu, C.; Wang, H. A high frequency vibration compensation approach in Terahertz SAR based on wavelet multi-resolution analysis. In Proceedings of the 2018 China International SAR Symposium (CISS), Shanghai, China, 10–12 October 2018; pp. 1–5. [Google Scholar]
  17. Li, Y.; Ding, L.; Zheng, Q.; Zhu, Y.; Sheng, J. A novel high frequency vibration error estimation and compensation algorithm for THz-SAR imaging based on local FrFT. Sensors 2020, 20, 2669. [Google Scholar] [CrossRef] [PubMed]
  18. Xia, W.; Huang, L. Target vibration estimation in SAR based on phase-analysis method. EURASIP J. Adv. Signal Process 2016, 94, 1–11. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Y.; Wang, Z.; Zhao, B.; Xu, L. Compensation for high-frequency vibration of platform in SAR imaging based on adaptive chirplet decomposition. IEEE Geosci. Remote Sens. Lett. 2016, 13, 792–795. [Google Scholar] [CrossRef]
  20. Djurović, I.; Popović-Bugarin, V.; Simeunović, M. The STFT-based estimator of micro-doppler parameters. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1273–1283. [Google Scholar] [CrossRef]
  21. Gao, Y.; Zhang, Z.; Xing, M.; Zhang, Y.; Li, Z. Paired echo suppression algorithm in helicopter-borne SAR imaging. IET Radar Sonar Navig. 2017, 11, 1605–1612. [Google Scholar] [CrossRef]
  22. Wang, Y.; Wang, Z.; Zhao, B.; Xu, L. Enhancement of azimuth focus performance in high-resolution SAR imaging based on the compensation for sensors platform vibration. IEEE Sens. J. 2016, 16, 6333–6345. [Google Scholar] [CrossRef]
  23. Wang, Y.; Wang, Z.; Zhao, B.; Xu, L. Parameters estimation of sinusoidal frequency modulation signal with application in synthetic aperture radar imaging. J. Appl. Remote Sens. 2016, 10, 020502. [Google Scholar] [CrossRef]
  24. Li, Y.; Wu, Q.; Wu, J.; Li, P.; Zheng, Q.; Ding, L. Estimation of high-frequency vibration parameters for Terahertz SAR imaging based on FrFT with combination of QML and RANSAC. IEEE Access 2021, 9, 5485–5496. [Google Scholar] [CrossRef]
  25. Djurović, I.; Stanković, L.J. STFT-based estimator of polynomial phase signals. Signal Process 2012, 92, 2769–2774. [Google Scholar] [CrossRef]
  26. Chum, O.; Matas, J. Optimal randomized RANSAC. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 30, 1472–1482. [Google Scholar] [CrossRef] [Green Version]
  27. Rahman, M.; Li, X.; Yin, X. DL-RANSAC: An improved RANSAC with modified sampling strategy based on the likelihood. In Proceedings of the 2019 IEEE 4th International Conference on Image, Vision and Computing (ICIVC), Xiamen, China, 5–7 July 2019; pp. 463–468. [Google Scholar]
  28. Caris, M.; Stanko, S.; Palm, S.; Sommer, R.; Wahlen, A.; Pohl, N. 300 GHz radar for high resolution SAR and ISAR applications. In Proceedings of the 2015 16th International Radar Symposium (IRS), Dresden, Germany, 24–26 June 2015; pp. 577–580. [Google Scholar]
  29. Wahl, D.E.; Eichel, P.H.; Ghiglia, D.C.; Jakowatz, C.V., Jr. Phase gradient autofocus-A robust tool for high resolution SAR phase correction. IEEE Trans. Aerosp. Electron. Syst. 1994, 30, 827–834. [Google Scholar] [CrossRef] [Green Version]
  30. Pei, S.; Huang, S. STFT with adaptive window width based on the chirp rate. IEEE Trans. Signal Process. 2012, 60, 4065–4080. [Google Scholar] [CrossRef] [Green Version]
  31. Kwok, H.K.C.; Jones, D.L. Instantaneous frequency estimation using an adaptive short-time Fourier transform. In Proceedings of the Conference Record of The Twenty-Ninth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 30 October–1 November 1995; pp. 543–547. [Google Scholar]
  32. Kajbaf, H.; Case, J.T.; Zheng, Y.R.; Kharkovsky, S.; Zoughi, R. Quantitative and qualitative comparison of SAR images from incomplete measurements using compressed sensing and nonuniform FFT. In Proceedings of the 2011 IEEE Radar Con (RADAR), Kansas City, MO, USA, 23–27 May 2011; pp. 592–596. [Google Scholar]
Figure 1. THz-SAR imaging geometric model.
Figure 1. THz-SAR imaging geometric model.
Applsci 11 10862 g001
Figure 2. Flow chart of the proposed vibration error compensation estimate method.
Figure 2. Flow chart of the proposed vibration error compensation estimate method.
Applsci 11 10862 g002
Figure 3. The estimation accuracy of signal parameters for different estimation methods. (a) vibration amplitude a 1 . (b) vibration frequency f 1 . (c) initial phase ϕ 1 .
Figure 3. The estimation accuracy of signal parameters for different estimation methods. (a) vibration amplitude a 1 . (b) vibration frequency f 1 . (c) initial phase ϕ 1 .
Applsci 11 10862 g003
Figure 4. The estimation accuracy of signal parameters for different estimation methods. (a) vibration amplitude a 1 ; (b) vibration frequency f 1 ; (c) initial phase ϕ 1 ; (d) vibration amplitude a 2 ; (e) vibration frequency f 2 ; and (f) initial phase ϕ 2 .
Figure 4. The estimation accuracy of signal parameters for different estimation methods. (a) vibration amplitude a 1 ; (b) vibration frequency f 1 ; (c) initial phase ϕ 1 ; (d) vibration amplitude a 2 ; (e) vibration frequency f 2 ; and (f) initial phase ϕ 2 .
Applsci 11 10862 g004
Figure 5. True and estimated single-component vibration error.
Figure 5. True and estimated single-component vibration error.
Applsci 11 10862 g005
Figure 6. Comparison of compensation effect. (a) No compensation; (b) compensation based on STFT with h = 20 ; (c) compensation based on STFT with h o p = 26 ; (d) compensation based on STFT with h o p = 26 and RANSAC.
Figure 6. Comparison of compensation effect. (a) No compensation; (b) compensation based on STFT with h = 20 ; (c) compensation based on STFT with h o p = 26 ; (d) compensation based on STFT with h o p = 26 and RANSAC.
Applsci 11 10862 g006
Figure 7. Phase estimate error.
Figure 7. Phase estimate error.
Applsci 11 10862 g007
Figure 8. Azimuth-dimensional amplitude profiles of center-point target.
Figure 8. Azimuth-dimensional amplitude profiles of center-point target.
Applsci 11 10862 g008
Figure 9. True and estimated multi-component vibration error.
Figure 9. True and estimated multi-component vibration error.
Applsci 11 10862 g009
Figure 10. Comparison of compensation effect. (a) No compensation; (b) compensation based on STFT with h = 20 ; (c) compensation based on STFT with h o p = 26 ; (d) compensation based on STFT with h o p = 26 and RANSAC.
Figure 10. Comparison of compensation effect. (a) No compensation; (b) compensation based on STFT with h = 20 ; (c) compensation based on STFT with h o p = 26 ; (d) compensation based on STFT with h o p = 26 and RANSAC.
Applsci 11 10862 g010
Figure 11. Phase estimate error.
Figure 11. Phase estimate error.
Applsci 11 10862 g011
Figure 12. Quality evaluation diagram of multi-component compensation.
Figure 12. Quality evaluation diagram of multi-component compensation.
Applsci 11 10862 g012
Table 1. System parameters.
Table 1. System parameters.
ParameterValue
carrier frequency220 GHz
signal bandwidth3 GHz
pulse width1.5 μs
platform velocity50 m/s
slant range2296 m
sample frequency2.8 GHz
pulse repetition frequency1050 Hz
Table 2. Estimated parameter values.
Table 2. Estimated parameter values.
Parameter a 1   ( mm ) f 1   ( Hz ) ϕ 1   ( rad ) G
true values2.001.001.05511.39
h = 202.309.960.85119.31
hop = 262.129.981.07430.24
hop = 26 with RANSAC2.029.981.05495.60
Table 3. Imaging quality evaluation after vibration error compensation.
Table 3. Imaging quality evaluation after vibration error compensation.
ParameterIRW (m)PSLR (dB)ISLR (dB)
theoretical values0.050−11.95−13.73
h = 200.064−0.21−5.38
hop = 260.053−7.10−8.57
hop = 26 with RANSAC0.051−12.37−11.48
Table 4. Estimated parameter values.
Table 4. Estimated parameter values.
Parametera1|a2 (mm)f1|f2 (Hz) ϕ 1 | ϕ 2   ( rad ) G
true values2.00|0.6010.00|20.001.05|0.52511.48
h = 201.85|0.669.96|20.230.93|0.62172.14
hop = 261.90|0.649.97|19.941.09|0.62386.56
hop = 26 with RANSAC2.01|0.599.97|19.941.06|0.56464.39
Table 5. Quality evaluation table of multi-component compensation.
Table 5. Quality evaluation table of multi-component compensation.
ParameterIRW (m)PSLR (dB)ISLR (dB)
theoretical values0.050−11.95−13.73
h = 200.074−2.29 −5.76
hop = 260.071−11.20−10.03
hop = 26 with RANSAC0.062−12.72−10.09
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Wu, Q.; Jiang, J.; Ding, X.; Zheng, Q.; Zhu, Y. A High-Frequency Vibration Error Compensation Method for Terahertz SAR Imaging Based on Short-Time Fourier Transform. Appl. Sci. 2021, 11, 10862. https://0-doi-org.brum.beds.ac.uk/10.3390/app112210862

AMA Style

Li Y, Wu Q, Jiang J, Ding X, Zheng Q, Zhu Y. A High-Frequency Vibration Error Compensation Method for Terahertz SAR Imaging Based on Short-Time Fourier Transform. Applied Sciences. 2021; 11(22):10862. https://0-doi-org.brum.beds.ac.uk/10.3390/app112210862

Chicago/Turabian Style

Li, Yinwei, Qi Wu, Jiawei Jiang, Xia Ding, Qibin Zheng, and Yiming Zhu. 2021. "A High-Frequency Vibration Error Compensation Method for Terahertz SAR Imaging Based on Short-Time Fourier Transform" Applied Sciences 11, no. 22: 10862. https://0-doi-org.brum.beds.ac.uk/10.3390/app112210862

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