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:
where
is the frequency of high-frequency vibration error and
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
where
is the amplitude of high-frequency vibration error and
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
axis is
. The period of vibration is
. The radar platform flies horizontally along the
direction with a speed
. The height of the radar platform is
. The coordinates of the radar platform at the azimuth moment
can be expressed as
. For arbitrary point target
in the imaging region, whose zero doppler time is
, the instantaneous slant range
can be expressed as
where
is the instantaneous slant range of the point target with the ideal motion state of radar platform and can be expressed as
where
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
, (3) can be approximated as
where
and is constant, and
is the view angle of the point target.
Suppose that radar antenna transmits a linear frequency modulation (LFM) signal, whose center frequency is
, pulse duration time is
, and chirp rate is
. After mixing demodulation, the echo signal of point target received by radar antenna can be expressed as
where
is the speed of light,
is the carrier wavelength, and
and
are the antenna pattern of range direction and azimuth direction, respectively.
After range compression, it can be obtained from (5) and (6)
where
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
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
After the dechirp operation and residual video phase (RVP) compensation [
20], (8) can be expressed as
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
where
is constant.
Suppose that
is the high-frequency vibration error along the line of sight, we can obtain that
On the one hand, if we take the second derivative of
directly, the vibration acceleration along the line-of-sight direction is
Therefore, if the vibration acceleration
and vibration frequency
along the line-of-sight direction are both obtained, the high-frequency vibration
along the line-of-sight direction can be estimated by
On the other hand, in a short time
,
can be approximated by quadratic polynomials, namely
where
is the initial vibration displacement,
is the initial vibration velocity and
is the vibration acceleration.
Therefore, in the short time
, (10) can be expressed as
From (15), we see that it is an LFM signal. If the chirp rate
has been obtained, the vibration acceleration
can be estimated through the following equation.
According to the definition of FrFT, the FrFT of the LFM signal
with matched order
is an impulse function at the fractional domain frequency
. At the same time, the matched order
and the fractional domain frequency
can be estimated by solving the following optimal equation [
21], i.e.,
where
is the
pth FrFT of the signal
. Therefore, the chirp rate
can be determined by the matched order
, i.e.,
By selecting the starting time continuously and performing FrFT on the azimuth signal in the corresponding time interval to estimate the local vibration acceleration , the vibration acceleration 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 , and the vibration frequency can be estimated through spectral analysis of the estimated vibration acceleration . 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 and divide the signal into segments;
(b) for each segment signal, the grading iterative search method [
22] is used to estimate the matched order
through (17);
(c) local vibration acceleration can be estimated by matched order 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 from the spectrum;
(g) the estimation of high-frequency vibration error 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
Step 2. Select the range gate with the dominant scatter, multiply it by the reference function
shown in (19) to complete the dechirp operation, and carry out RVP compensation, we can obtain the SFM signal shown in (10);
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
and multiply it by (18) to complete the high-frequency vibration error compensation;
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.
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.