Next Article in Journal
Multiple Electric Energy Consumption Forecasting Using a Cluster-Based Strategy for Transfer Learning in Smart Building
Previous Article in Journal
A Simple High-Resolution Near-Field Probe for Microwave Non-Destructive Test and Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel High-Frequency Vibration Error Estimation and Compensation Algorithm for THz-SAR Imaging Based on Local FrFT

1
Terahertz Technology Innovation Research Institute, University of Shanghai for Science and Technology, Shanghai 200093, China
2
Shanghai Key Laboratory of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China
3
Shanghai Institute of Intelligent Science and Technology, Tongji University, Shanghai 200092, China
4
Shanghai Radio Equipment Research Institute, Shanghai 201109, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(9), 2669; https://doi.org/10.3390/s20092669
Submission received: 17 March 2020 / Revised: 2 May 2020 / Accepted: 4 May 2020 / Published: 7 May 2020
(This article belongs to the Section Fault Diagnosis & Sensors)

Abstract

:
Compared with microwave synthetic aperture radar (SAR), terahertz SAR (THz-SAR) is easier to achieve ultrahigh-resolution image due to its higher frequency and shorter wavelength. However, higher carrier frequency makes THz-SAR image quality very sensitive to high-frequency vibration error of motion platform. Therefore, this paper proposes a novel high-frequency vibration error estimation and compensation algorithm for THz-SAR imaging based on local fractional Fourier transform (LFrFT). Firstly, the high-frequency vibration error of the motion platform is modeled as a simple harmonic motion and THz-SAR echo signal received in each range pixel can be considered as a sinusoidal frequency modulation (SFM) signal. A novel algorithm for the parameter estimation of the SFM signal based on LFrFT is proposed. The instantaneous chirp rate of the SFM signal is estimated by determining the matched order of LFrFT in a sliding small-time window and the vibration acceleration is obtained. Hence, the vibration frequency can be estimated by the spectrum analysis of estimated vibration acceleration. With the estimated vibration acceleration and vibration frequency, the SFM signal is reconstructed. Then, the corresponding THz-SAR imaging algorithm is proposed to estimate and compensate the phase error caused by the high-frequency vibration error of the motion platform and realize high-frequency vibration error estimation and compensation for THz-SAR imaging. Finally, the effectiveness of the novel algorithm proposed in this paper is demonstrated by simulation results.

1. Introduction

Terahertz waves usually refers to electromagnetic waves from 100 GHz to 10 THz, which lie between the microwave and infrared, known as the THz gap. However, recent research results show that the THz gap has great development potential and application value due to its unique electromagnetic wave band [1,2]. In recent years, THz waves have been widely used in explosives detection, drug detection, imaging, radar and wireless broadband communication [3,4,5,6,7].
Compared with conventional microwave SAR imaging, THz-SAR imaging has considerable advantages in two-dimensional or three-dimensional ultrahigh-resolution imaging and detection and the recognition of weak-small-slow maneuvering targets due to its higher frequency and shorter wavelength. Ultrahigh range resolution can be obtained with a large bandwidth signal, and ultrahigh azimuth resolution can be obtained with a large Doppler bandwidth. In practical engineering application, airborne SAR systems are usually influenced by atmospheric disturbance, which introduces motion errors in radar echo signals. Low-frequency motion errors can be measured by motion sensors and compensated for in imaging processing. Due to the small amplitude, the influence of high-frequency vibration error on microwave SAR imaging is generally negligible. However, a higher frequency also means that the imaging quality of THz-SAR is more sensitive to the vibration error of the motion platform. Although the high-frequency vibration error will not affect the track of THz-SAR, it will seriously influence the phase of the radar echo signal, thus worsening the image quality [8,9]. Therefore, high-frequency vibration error estimation and compensation is the key for THz-SAR ultrahigh resolution imaging.
Since the high-frequency vibration error of a motion platform can be modeled as simple harmonic motion, the phase of THz-SAR-received echo signal in each range gate is modulated periodically. Meanwhile, the received echo signal of THz-SAR in each range gate is considered as sinusoid frequency modulation (SFM). Different from low-frequency motion error, high-frequency vibration error may produce a paired echo due to periodic modulation of the echo signal phase, which seriously affects the image quality. When the vibration frequency is very high, the traditional motion compensation algorithms will fail [10,11]. In [12], a THz-SAR imaging method under helicopter platform vibration was proposed. The compensation function is constructed by modern signal processing methods such as short-time Fourier transform and parametric spatial projection after pair echo is focused by adopting Keystone transform. Due to the non-stationary characteristics of SFM signals, the commonly used time-frequency analysis method cannot effectively analyze SFM signal. In [13], an SFM signal parameter estimation algorithm based on discrete sinusoidal frequency modulation transformation was proposed. It needs to search in three-dimensional space, resulting in huge computational burden. In [14], the simulated annealing algorithm was used to reduce the computational load, but converging to a local minimum occasionally is inevitable when it is used in practice, which leads to a large error in the parameter estimation results. In [15,16], sinusoidal frequency modulation Fourier transform and sinusoidal frequency modulation Fourier-Bessel transform were proposed to estimate SFM signal parameters, respectively. However, both these two algorithms directly make use of the SFM signal phase, which causes the algorithm performance to be seriously affected by signal-to-noise ratio (SNR) due to the phase unwrapping step.
Fractional Fourier transform (FrFT) is considered as a transform of signal from time domain into fractional Fourier domain (FrFD) between time and frequency [17]. As the FrFT of a signal is regarded as the decomposition in terms of an orthonormal basis set formed by chirp signals in FrFD, FrFT is a more suitable analysis tool for non-stationary signal processing. However, FrFT could not be directly used for the detection and parameter estimation of SFM signals. This is due to the fact that the phase order of SFM signals is more than two over the entire sample period and FrFT of SFM signals cannot have a good energy concentration property. Therefore, as the short-time Fourier transform (STFT) can determine the time–frequency relation of a time-varying signal, performing FrFT on SFM signals in a small-time window can determine the time–chirp rate relation. Through searching matched order in each sliding window, the relation of time–chirp rate of SFM signals can be obtained. Therefore, a novel algorithm for the parameter estimation of an SFM signal based on local FrFT (LFrFT) is proposed. The instantaneous chirp rate of SFM signal is determined by the matched order of LFrFT in sliding small-time window and the vibration acceleration is obtained. Hence, the vibration frequency is estimated by the spectrum analysis of vibration acceleration. With the estimated vibration acceleration and vibration frequency, the SFM signal is reconstructed and the high-frequency vibration error is obtained. Then, the corresponding THz-SAR imaging algorithm is proposed to estimate and compensate the phase error caused by high-frequency vibration error of motion platform and realize THz-SAR high-frequency vibration error estimation and compensation imaging. Finally, simulation experiments are made to verify the effectiveness of high-frequency vibration error estimation method and THz-SAR imaging algorithm with high-frequency vibration error.
This paper is organized as follows. THz-SAR imaging model with a high-frequency vibration error is established in Section 2. In Section 3 and Section 4, the LFrFT-based high-frequency vibration estimation method and the THz-SAR high-frequency vibration error compensation imaging algorithm are proposed, respectively. In Section 5, simulation experimental results are given to demonstrate the effectiveness of the proposed algorithm. Finally, a conclusion to this paper is drawn in Section 6.

2. THz-SAR Imaging Model with High-Frequency Vibration Error

In airborne SAR systems, the radar platform is usually influenced by atmospheric disturbance, which introduces motion errors in the radar echo signal. Among them, high-frequency vibration error is generally referred to meet the following conditions:
| f v T s | 1
where f v is the frequency of high-frequency vibration error and T s is the synthetic aperture time of the imaging target.
The high-frequency vibration error of aircraft platform is close to a superposition of multiple harmonic motions with an octave relationship. The platform vibration amplitude of a helicopter is the composition of each order’s harmonic. However, the vibration amplitude of the first harmonic plays a leading role in the platform vibration [18]. Thus, only the first harmonic vibration is considered here. The high-frequency vibration error can be expressed as
d ( t ) = A v s i n ( 2 π f v t + φ v )
where A v is the amplitude of high-frequency vibration error and φ v is the initial phase of high-frequency vibration error.
Considering that the influence of high-frequency vibration error perpendicular to the line of sight on the imaging quality of THz-SAR can be ignored [19], only high-frequency vibration error in the plane along the line of sight is considered here, as shown in Figure 1. The high-frequency vibration error is in the YOZ plane, and the angle between the vibration direction and the Z axis is ϕ . The period of vibration is T v = 1 / f v . The radar platform flies horizontally along the X direction with a speed v . The height of the radar platform is H . The coordinates of the radar platform at the azimuth moment t can be expressed as ( v t , d s i n ϕ , H + d c o s ϕ ) . For arbitrary point target P ( v t 0 , y 0 , 0 ) in the imaging region, whose zero doppler time is t 0 , the instantaneous slant range R ( t ) can be expressed as
R ( t ) = v 2 ( t t 0 ) 2 + ( y 0 + d s i n ϕ ) 2 + ( H + d c o s ϕ ) 2 = r s 2 ( t ) + d 2 + 2 ( H c o s ϕ + y 0 s i n ϕ ) d
where r s ( t ) is the instantaneous slant range of the point target with the ideal motion state of radar platform and can be expressed as
r s ( t ) = v 2 ( t t 0 ) 2 + y 0 2 + H 2 = v 2 ( t t 0 ) 2 + r 0 2
where r 0 is the nearest slant range of the point target.
Since the synthetic aperture length of THz-SAR is short and the amplitude of high-frequency vibration error is far less than the nearest slant range r 0 , (3) can be approximated as
R ( t ) r 0 + v 2 ( t t 0 ) 2 2 r 0 + ( H c o s ϕ + y 0 s i n ϕ ) r s ( t ) d r 0 + v 2 ( t t 0 ) 2 2 r 0 + C d
where C = c o s ( θ 0 ϕ ) < 1 and is constant, and θ 0 = a r c t a n ( y 0 H ) is the view angle of the point target.
Suppose that radar antenna transmits a linear frequency modulation (LFM) signal, whose center frequency is f c , pulse duration time is T , and chirp rate is K r . After mixing demodulation, the echo signal of point target received by radar antenna can be expressed as
s ( τ , t ; r 0 , t 0 ) = w r ( τ 2 R ( t ) c ) w a ( t t 0 ) e x p { j 4 π λ R ( t ) } e x p { j π K r ( τ 2 R ( t ) c ) 2 }
where c is the speed of light, λ is the carrier wavelength, and w r and w a are the antenna pattern of range direction and azimuth direction, respectively.
After range compression, it can be obtained from (5) and (6)
s ( τ , t ; r 0 , t 0 ) = s i n c ( π B ( τ 2 r 0 c v 2 ( t t 0 ) 2 c r 0 2 C d c ) ) w a ( t t 0 ) e x p { j 4 π λ [ r 0 + v 2 ( t t 0 ) 2 2 r 0 + C d ] }
where B is the signal bandwidth.
The amplitude of high-frequency vibration error is very small, usually in the millimeter or sub-millimeter level, and the range resolution of THz-SAR is above the centimeter level. As the parameters listed in the simulation experiment, the vibration error amplitude is 0.5 mm and less than 1 / 8 of the range resolution, 6.64 cm. Therefore, the range cell migration caused by high-frequency vibration error can be ignored. After performing the range cell migration correction through interpolation in the range-doppler domain, (7) can be expressed as
s ( τ , t ; r 0 , t 0 ) = s i n c ( π B ( τ 2 r 0 c ) ) w a ( t t 0 ) e x p { j 4 π λ [ r 0 + v 2 ( t t 0 ) 2 2 r 0 + C d ] }
After the dechirp operation and residual video phase (RVP) compensation [20], (8) can be expressed as
s ( τ , t ; r 0 , t 0 ) = s i n c ( π B ( τ 2 r 0 c ) ) w a ( t t 0 ) e x p { j [ 4 π λ r 0 2 π 2 v 2 t 0 2 λ r 0 + 4 π λ C d ] }
For (9), the spectral analysis algorithm using the fast Fourier transform (FFT) operation can focus the azimuth signal in the doppler domain. However, due to the existence of high-frequency vibration error, the focused image will generate paired echoes, leading to a decrease in the image quality. If the parameters of the high-frequency vibration error are known, the effect can be eliminated by the phase compensation algorithm, thus ensuring the imaging quality of THz-SAR.

3. LFrFT-Based High-Frequency Vibration Error Estimation Method

According to (9), the phase error of an echo signal caused by high-frequency vibration error in each range gate can be regarded as the SFM signal, which can be expressed as
g ( t ) = A e x p { j 2 π 2 λ C 1 s i n ( 2 π f v t + φ v ) }
where C 1 = C A v is constant.
Suppose that d L O S ( t ) is the high-frequency vibration error along the line of sight, we can obtain that
d L O S ( t ) = C d = C 1 s i n ( 2 π f v t + φ v )
On the one hand, if we take the second derivative of d L O S ( t ) directly, the vibration acceleration along the line-of-sight direction is
a ( t ) = 2 d L O S ( t ) t 2 = 4 π 2 f v 2 C 1 s i n ( 2 π f v t + φ v ) = 4 π 2 f v 2 d L O S ( t )
Therefore, if the vibration acceleration a ^ ( t ) and vibration frequency f ^ v along the line-of-sight direction are both obtained, the high-frequency vibration d ^ L O S ( t ) along the line-of-sight direction can be estimated by
d ^ L O S ( t ) = a ^ ( t ) 4 π 2 f ^ v 2
On the other hand, in a short time [ t 0 , t 0 + Δ t ] , d L O S ( t ) can be approximated by quadratic polynomials, namely
d L O S ( t ) = x ( t 0 ) + v ( t 0 ) t + 1 2 a ( t 0 ) t 2
where x ( t 0 ) is the initial vibration displacement, v ( t 0 ) is the initial vibration velocity and a ( t 0 ) is the vibration acceleration.
Therefore, in the short time [ t 0 , t 0 + Δ t ] , (10) can be expressed as
g ( t ) = A e x p { j 2 π [ 2 λ x ( t 0 ) + 2 λ v ( t 0 ) t + 1 λ a ( t 0 ) t 2 ] }
From (15), we see that it is an LFM signal. If the chirp rate c ^ ( t 0 ) has been obtained, the vibration acceleration a ( t 0 ) can be estimated through the following equation.
a ^ ( t 0 ) = λ 2 c ^ ( t 0 )
According to the definition of FrFT, the FrFT of the LFM signal y ( t ) with matched order p 0 is an impulse function at the fractional domain frequency u 0 . At the same time, the matched order p 0 and the fractional domain frequency u 0 can be estimated by solving the following optimal equation [21], i.e.,
{ p ^ 0 , u ^ 0 } = arg m a x p , u | F p { y ( t ) } |
where F p { y ( t ) } is the pth FrFT of the signal y ( t ) . Therefore, the chirp rate c ( t 0 ) can be determined by the matched order p 0 , i.e.,
c ^ ( t 0 ) = c o t ( π 2 p ^ 0 )
By selecting the starting time t 0 continuously and performing FrFT on the azimuth signal in the corresponding time interval [ t 0 , t 0 + Δ t ] to estimate the local vibration acceleration a ( t 0 ) , the vibration acceleration a ^ ( t ) of the whole azimuth signal can be obtained. At the same time, it can be seen from (12) that the frequency of vibration acceleration along the line-of-sight direction is also f v , and the vibration frequency f ^ v can be estimated through spectral analysis of the estimated vibration acceleration a ^ ( t ) . Since the FrFT is carried out in the local region of azimuth signal, we call it the local FrFT, denoted as LFrFT.
To sum up, the procedure of the LFrFT-based high-frequency vibration error estimation method is as follows:
(a) select the appropriate time interval Δ t and divide the signal g ( t ) into N segments;
(b) for each segment signal, the grading iterative search method [22] is used to estimate the matched order p ^ 0 through (17);
(c) local vibration acceleration a ^ ( t 0 ) can be estimated by matched order p ^ 0 through (16) and (18);
(d) repeat steps (b) and (c) until the vibration acceleration of all segments’ signal is estimated;
(e) adopting the moving average method to filter the vibration acceleration to improve the estimation accuracy;
(f) perform FFT on the vibration acceleration and estimate the vibration frequency f ^ v from the spectrum;
(g) the estimation of high-frequency vibration error d ^ ( t ) can be obtained through (11) and (12).

4. High-Frequency Vibration Error Compensation Imaging Algorithm

It can be seen from (9) that conventional imaging algorithms such as range Doppler (RD) and chirp scaling (CS) will generate paired echoes in SAR images and reduce the image quality when high-frequency vibration error exists in THz-SAR platform. Therefore, combined with the high-frequency vibration estimation method based on LFrFT, the high-frequency vibration compensation imaging algorithm of THz-SAR is proposed as follows:
Step 1. After range compression and range cell migration correction, the echo signal can be expressed as
s ( τ , t ; r 0 , t 0 ) = s i n c ( π B ( τ 2 r 0 c ) ) w a ( t t 0 ) e x p { j 4 π λ [ r 0 + v 2 ( t t 0 ) 2 2 r 0 + C A v s i n ( 2 π f v t + φ v ) ] }
Step 2. Select the range gate with the dominant scatter, multiply it by the reference function s r e f shown in (19) to complete the dechirp operation, and carry out RVP compensation, we can obtain the SFM signal shown in (10);
s r e f = e x p { j 4 π λ [ r 0 + v 2 t 2 2 r 0 ] }
Step 3. Perform the LFrFT-based estimation method proposed in Section 3 and obtain the high-frequency vibration error;
Step 4. Adopt the estimated high-frequency vibration error to reconstruct the compensation signal s c o m and multiply it by (18) to complete the high-frequency vibration error compensation;
s c o m = e x p { j 4 π λ C A s i n ( 2 π f v t + φ v ) }
Step 5. Azimuth compression is done through matched filtering or dechirp operation, and the well-focused image is obtained.
Therefore, the flowchart of THz-SAR high-frequency vibration error estimation and compensation imaging algorithm is shown in Figure 2.

5. Simulation Experiments

In the section, the high-frequency vibration error estimation algorithm based on LFrFT is used to estimate the high-frequency vibration error, and then the new proposed THz-SAR high-frequency vibration compensation imaging algorithm is adopted for focusing on the imaging of a point target and an area target. The THz-SAR system’s simulation parameters are given in Table 1.

5.1. Simulation Analysis of LFrFT-Based High-Frequency Vibration Error Estimation

Considering the high-frequency vibration case with the form of (10) and supposing that, f v = 20   Hz ,   C 1 = 5 × 10 4   m and φ v = 0 . It can be seen from Table 1 that the synthetic aperture time of the THz-SAR system is 0.4 s. The high-frequency vibration acceleration estimated by LFrFT and its spectrum are shown in Figure 3 and Figure 4, respectively. From the comparison between the estimated vibration acceleration and the true vibration acceleration, it can be found that the estimation error of vibration acceleration is larger at the extreme value points. This is due to the fact that the approximation error is larger when the second order approximation is adopted. At the same time, the estimated values of vibration acceleration after moving the average filter is smoother and closer to the true value of vibration acceleration than that of vibration acceleration before filtering. It is obvious from the spectrum diagram that the frequency of vibration acceleration is 20.04 Hz, which is almost consistent with the actual vibration frequency. Therefore, combining the estimated vibration acceleration and vibration frequency, the estimated vibration displacement can be obtained, as shown in Figure 5. According to the comparison between the estimated value and the true value of displacement, the LFrFT-based high-frequency vibration estimation method can estimate the SFM signal well.
In the following, we compare the performance of the proposed method with that of discrete sinusoid frequency modulation transform (DSFMT)-based method in [14]. One metric to evaluate the estimation performance is the normalized root-mean-square error (NRMSE), as in [15]. It is defined as
NRMSE = d e s t i m a t e ( t ) d t r u e ( t ) d t r u e ( n )
where denotes the 2 -norm operator. The NRMSE values the normalized Euclidean distance between the estimated high-frequency vibration error and the true one.
In the simulation experiment, the SNR of the signal is varied in 5 dB increments between 0 and 15 dB. For each SNR value, we apply the LFrFT-based method and DSFMT-based method to obtain the high-frequency vibration error, respectively, and 100 simulations are performed. Table 2 shows the average NRMSE of estimated high-frequency vibration error for LFrFT-based method and DSFMT-based method. From Table 2, we can see that the estimation accuracy degrades slowly with the decrease in the SNR, and the estimate performance of our proposed method is better than that of DSFMT-based method.

5.2. Simulation Analysis of High-Frequency Vibration Error Estimation and Compensation Imaging Algorithm

The THz-SAR system parameters in Table 1 and the high-frequency vibration error in simulation 5.1 are used to generate the echo signal for a single point target. Firstly, the RD algorithm is used for focusing on imaging, as shown in Figure 6a. It can be seen that the single point target has many paired echoes in the case of high-frequency vibration error. Then, the proposed high-frequency vibration error estimation and compensation imaging algorithm is adopted for focusing on imaging, as shown in Figure 6b. It is obvious that the point target is well focused and paired echoes is eliminated after the high-frequency vibration error is estimated and compensated.
In order to test the focus quality of the point target more intuitively, a rectangular window is adopted to intercept the point target in Figure 6b. After up-sampling eight times, the resulting two-dimensional sections (amplitude and phase) are shown in Figure 7, respectively. As can be seen from the two profiles, the range resolution and the azimuth resolution are 0.0664 and 0.0696 m, the peak sidelobe ratios are −13.2643 and −11.8592 dB, and the integral sidelobe ratios are −30.2416 dB and −25.6668 dB, respectively. The theoretical range resolution and azimuth resolution are known to be 0.0664 and 0.0692 m, respectively. Hence, the point target is well focused. Since the estimation of high-frequency vibration error is not ideal, the residual error after compensation makes the azimuth peak sidelobe ratio slightly higher. Therefore, it is necessary to further improve the estimation accuracy of the high-frequency vibration error or adopt the window function for suppressing the sidelobe.
Then, the multi-point targets are used for the simulation experiment. The THz-SAR imaging results of the RD algorithm and the new proposed imaging algorithm in this paper are showed in Figure 8a,b, respectively. It can be seen from Figure 8 that the new proposed imaging algorithm in this paper is still applicable to the case of multiple point targets and can effectively suppress paired echoes.
Finally, since there are no measured data of THz-SAR, the area targets are simulated to generate echo data in the experiment. The THz-SAR imaging results of the RD algorithm and the new proposed imaging algorithm in this paper are showed in Figure 9b,c, respectively. It can be seen that the new proposed imaging algorithm in this paper can effectively suppress paired echoes. However, the image quality after high-frequency vibration error compensation is slightly worse than the raw image shown in Figure 9a due to the residual error.
Remarks: The estimation accuracy of vibration displacement directly affects the focusing quality of THz-SAR image and is determined by the estimation accuracy of vibration acceleration with the LFrFT-based method. However, due to the existence of many factors, such as the signal-to-clutter-noise ratio of selected dominant scatter signal and window size of LFrFT, the estimation of high-frequency vibration displacement is not ideal, as shown in Table 2. After compensation with estimated high-frequency vibration displacement, the paired echoes in azimuth are suppressed in the THz-SAR image, but the image quality may still be a little defocused due to the residual error, as shown in Figure 9. In actual THz-SAR image processing, the residual phase error can be further extracted and compensated by using existing autofocus algorithms such as PGA.

6. Conclusions

Considering that THz-SAR image quality is sensitive to high-frequency vibration error of the motion platform, a novel THz-SAR high-frequency vibration error estimation and compensation imaging algorithm based on LFrFT is proposed in this paper. The THz-SAR echo signal received in each range pixel is modeled as SFM signal. The instantaneous chirp rate of SFM signal is estimated by the matched order of LFrFT in sliding small-time window and the vibration acceleration is obtained. Meanwhile, the vibration frequency is estimated by the spectrum analysis of vibration acceleration. With the estimated vibration acceleration and vibration frequency, high-frequency vibration error of motion platform is reconstructed. Then, the corresponding THz-SAR imaging algorithm is proposed to compensate the phase error caused by high-frequency vibration error of motion platform and realize THz-SAR high-frequency vibration error estimation and compensation imaging. Finally, the effectiveness of LFrFT-based high-frequency vibration error estimation and compensation imaging algorithm is verified by simulation experiments.
Although this paper proposes a corresponding solution to the problem of high-frequency vibration error for THz-SAR, the solution is also applicable to other wavebands, such as the millimeter band. In fact, the high-frequency vibration error of THz-SAR platform may contain a variety of vibration frequency components. Therefore, it is necessary to further study the high-frequency vibration compensation imaging algorithm of THz-SAR in the case of multi-frequency vibration error. The premise of the proposed method is that there are dominant scatter points in the imaging scene. The case where there is no dominant scatter point is also worth studying in future work.

Author Contributions

Conceptualization, Y.L.; methodology, Y.L.; funding acquisition, Y.Z. and J.S.; software, Y.L.; validation, Y.L. and Q.Z.; resources, L.D. and Y.Z.; writing—original draft preparation, Y.L.; writing—review and editing, Y.L., Q.Z. and L.D.; project administration, Y.Z. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National key R&D Program of China (2018YFF01013003), in part by the National Natural Science Foundation of China (41601508, 61722111, 61731020, 61605113, 61601291), the 111 Project (D18014), the International Joint Lab Program supported by Science and Technology Commission Shanghai Municipality (17590750300), Natural Science Foundation of Shanghai (17ZR1428700).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zandonella, C. Terahertz imaging: T-ray specs. Nature 2003, 424, 721–722. [Google Scholar] [CrossRef] [PubMed]
  2. Redo, S.A.; Zhang, X.C. Terahertz science and technology trends. IEEE J. Sel. Quantum Elect 2008, 14, 260–269. [Google Scholar] [CrossRef]
  3. Balacey, H.; Recur, B.; Perraud, J.-B.; Sleiman, J.B.; Guillet, J.-P.; Mounaix, P. Advanced processing sequence for 3-D THz imaging. IEEE Trans. THz Sci. Technol. 2016, 6, 191–198. [Google Scholar] [CrossRef]
  4. Gao, J.; Cui, Z.; Cheng, B.; Qin, Y.; Deng, X.; Deng, B.; Li, X.; Wang, H. Fast three-dimensional image reconstruction of a standoff screening system in the terahertz regime. IEEE Trans. THz Sci. Technol. 2018, 8, 38–51. [Google Scholar] [CrossRef]
  5. Rappaport, T.S.; Xing, Y.; Kanhere, O.; Ju, S.; Madanayake, A.; Mandal, S.; Alkhateeb, A.; Trichopoulos, G.C. Wireless communications and applications above 100 GHz: Opportunities and challenges for 6G and beyond. IEEE Access 2019, 7, 78729–78757. [Google Scholar] [CrossRef]
  6. Gao, H.; Li, C.; Zheng, S.; Wu, S.; Fang, G. Implementation of the phase shift migration in MIMO-sidelooking imaging at terahertz band. IEEE Sens. J. 2019, 19, 9384–9393. [Google Scholar] [CrossRef]
  7. 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]
  8. 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. Proceedings of 2018 China International SAR Symposium (CISS), Shanghai, China, 10–12 October 2018; pp. 1–5. [Google Scholar]
  9. Batra, A.; Wiemeler, M.; Kreul, T.; Goehringer, D.; Kaiser, T. SAR signal processing architecture and effects of motion errors for mmWave and THz frequencies. Proceedings of 2019 Second International Workshop on Mobile Terahertz Systems (IWMTS), Bad Neuenahr, Germany, 1–3 July 2019; pp. 1–6. [Google Scholar]
  10. Xing, M.; Jiang, X.; Wu, R.; Zhou, F.; Bao, Z. Motion compensation for UAV SAR based on raw radar data. IEEE Trans. Geosci. Remote Sens. 2008, 47, 2870–2883. [Google Scholar] [CrossRef]
  11. Chen, J.; Liang, B.; Yang, D.; Zhao, D.; Xing, M.; Sun, G. Two-step accuracy improvement of motion compensation for airborne SAR with ultrahigh resolution and wide swath. IEEE Trans. Geosci. Remote Sens. 2019, 57, 7148–7160. [Google Scholar] [CrossRef]
  12. Wang, H.; Zhang, Y.; Wang, B.B.; Sun, J.P. A novel helicopter-borne Terahertz SAR imaging algorithm based on keystone transform. Proceedings of 2014 12th International Conference on Signal Processing (ICSP), Hangzhou, China, 19–23 October 2014; pp. 1958–1962. [Google Scholar]
  13. Sun, Z.G.; Chen, J.; Cao, X.; Guo, L.L.; Shen, L.R. Parameter estimation of multicomponent SFM signals based on discrete sinusoidal frequency modulation transform. System Eng. Electron. 2012, 34, 1973–1979. [Google Scholar]
  14. Wang, Y.; Wang, Z.F.; Zhao, B.; Xu, L. Enhancement of azimuth focus performance in terahertz SAR imaging based on the compensation for high frequency vibration of platform. IEEE Sens. J. 2016, 16, 6333–6345. [Google Scholar] [CrossRef]
  15. Peng, B.; Wei, X.Z.; Deng, B.; Chen, H.W.; 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]
  16. He, Q.F.; Zhang, Q.; Luo, Y.; Sun, L. Sinusoidal frequency modulation Fourier-Bessel series for multicomponent SFM signal estimation and separation. Mathematical Problems in Eng. 2017, 2017, 1–14. [Google Scholar] [CrossRef] [Green Version]
  17. Almeida, L.B. The fractional Fourier transform and time-frequency representations. IEEE Trans. Signal Process. 1994, 42, 3084–3091. [Google Scholar] [CrossRef]
  18. Ceruti, A.; Liverani, A.; Recanatesi, L. Improving helicopter flight simulation with rotor vibrations. Proceedings of International Conference on Innovative Methods in Product Design, Venice, Italy, 15–17 June 2011; pp. 636–645. [Google Scholar]
  19. Zhao, Y.L.; Zhang, Q.Y.; Li, C. Vibration error analysis and motion compensation of video synthetic aperture radar. J. Radars 2015, 4, 230–239. [Google Scholar]
  20. Zhang, Y.; Sun, J.P.; Lei, P.; Hong, W. SAR-based paired echo focusing and suppression of vibrating targets. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7593–7605. [Google Scholar] [CrossRef]
  21. Qi, L.; Tao, R.; Zhou, S.Y. Detection and parameter estimation of multi-component LFM signals based on the fractional Fourier transform. Sci. China Series F Inf. Sci. 2004, 47, 184–198. [Google Scholar] [CrossRef]
  22. Chen, X.L.; Guan, J. A fast FRFT based detection algorithm of multiple moving targets in sea clutter. Proceedings of International Conference on Radar, Washington, DC, USA, 10–14 May 2010; pp. 402–406. [Google Scholar]
Figure 1. Terahertz synthetic aperture radar (THz-SAR) imaging geometry with high-frequency vibration error.
Figure 1. Terahertz synthetic aperture radar (THz-SAR) imaging geometry with high-frequency vibration error.
Sensors 20 02669 g001
Figure 2. THz-SAR high-frequency vibration compensation imaging algorithm based on local fractional Fourier transform (LFrFT).
Figure 2. THz-SAR high-frequency vibration compensation imaging algorithm based on local fractional Fourier transform (LFrFT).
Sensors 20 02669 g002
Figure 3. Real and estimated values of high-frequency vibration acceleration.
Figure 3. Real and estimated values of high-frequency vibration acceleration.
Sensors 20 02669 g003
Figure 4. Frequency spectrum of high-frequency vibration acceleration.
Figure 4. Frequency spectrum of high-frequency vibration acceleration.
Sensors 20 02669 g004
Figure 5. Real and estimated values of high-frequency vibration displacement.
Figure 5. Real and estimated values of high-frequency vibration displacement.
Sensors 20 02669 g005
Figure 6. Point target image focused by (a) RD algorithm; (b) the proposed algorithm.
Figure 6. Point target image focused by (a) RD algorithm; (b) the proposed algorithm.
Sensors 20 02669 g006
Figure 7. Focusing quality analysis of point target. (a) Range; (b) azimuth.
Figure 7. Focusing quality analysis of point target. (a) Range; (b) azimuth.
Sensors 20 02669 g007
Figure 8. Imaging results of multi-point targets. (a) RD algorithm; (b) the propose algorithm.
Figure 8. Imaging results of multi-point targets. (a) RD algorithm; (b) the propose algorithm.
Sensors 20 02669 g008
Figure 9. Imaging results of area targets. (a) Raw Image; (b) RD algorithm; (c) the proposed algorithm.
Figure 9. Imaging results of area targets. (a) Raw Image; (b) RD algorithm; (c) the proposed algorithm.
Sensors 20 02669 g009
Table 1. System parameter table.
Table 1. System parameter table.
ParameterValue
Center frequency200 GHz
Platform height2000 m
Platform speed50 m/s
View angle30°
Pulse width1.5 µs
Signal bandwidth2.0 GHz
Sample frequency2.5 GHz
Pulse repetition frequency1000 Hz
Table 2. Average normalized root-mean-square error (NRMSE) of estimated high-frequency vibration error under typical signal-to-noise ratios (SNRs).
Table 2. Average normalized root-mean-square error (NRMSE) of estimated high-frequency vibration error under typical signal-to-noise ratios (SNRs).
SNRAverage NRMSE of DSFMTAverage NRMSE of LFrFT
0 dB0.22850.1973
5 dB0.17710.1234
10 dB0.08290.0678
15 dB0.05780.0352

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/s20092669

AMA Style

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(9):2669. https://doi.org/10.3390/s20092669

Chicago/Turabian Style

Li, Yinwei, Li Ding, Qibin Zheng, Yiming Zhu, and Jialian Sheng. 2020. "A Novel High-Frequency Vibration Error Estimation and Compensation Algorithm for THz-SAR Imaging Based on Local FrFT" Sensors 20, no. 9: 2669. https://doi.org/10.3390/s20092669

APA Style

Li, Y., Ding, L., Zheng, Q., Zhu, Y., & Sheng, J. (2020). A Novel High-Frequency Vibration Error Estimation and Compensation Algorithm for THz-SAR Imaging Based on Local FrFT. Sensors, 20(9), 2669. https://doi.org/10.3390/s20092669

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