Abstract

Synthetic aperture radar (SAR) was originally exploited to image stationary scenes. However, it is important to derive target information of velocity for many applications. The fractional Fourier transform (FrFT) is a generalization of the classical Fourier transform and is well-known as a useful tool to estimate the chirp rate of linear frequency-modulated (LFM) signals. Motion compensation is critical to moving target imaging. It is difficult for us to obtain the actual motion parameters in real scenarios. Based on the moving target echo model in airborne along-track interferometric SAR (ATI-SAR) and expression of the ATI phase, a method is proposed to estimate the ship velocity by combining the ATI phase with FrFT. First, we use the FrFT to evaluate the chirp rate of the moving target echo. Then, we construct an equation to estimate the ship velocity using the chirp rate estimation, peak response time, and ATI phase. Finally, the simulation experiments are used to validate the effectiveness of the proposed method.

1. Introduction

Sea traffic monitoring is one of the most important applications of synthetic aperture radar (SAR) imagery. Moving ship detection is an important issue in sea traffic monitoring, and imagery quality has a great impact on detection performance. Unfortunately, conventional imaging algorithms, such as Range-Doppler (RD), Chirp Scaling (CS), and SPECAN [1], have focused on static scenes. The performances of these imaging algorithms for moving targets worsen due to the motion of the target. SAR signatures of moving targets suffer from azimuthal displacement, smearing, defocusing caused by motion in the range direction, and loss of focus due to azimuthal velocity [2]. One solution that can address this problem for moving target imaging is motion compensation. By combining motion compensation with conventional imaging algorithms, it is possible to focus on moving targets. Motion compensation relies on the acquiring of motion parameters, i.e., velocities, for a target. Generally, utilizing the Automatic Identification System (AIS) conveniently provides the motion parameters of a moving target. Unfortunately, in many cases, no AIS data are available for moving targets. Therefore, an alternative method to realize velocity estimations should be explored.

Research on velocity estimation was first applied using ground moving target indication (GMTI). The challenge for GMTI is separating the returned signal of the moving target from the stationary background [3]. However, compared with a target moving across the ground, the stationary hypothesis is invalid in the wide ocean area due to different phenomena, such as ocean waves or moving ships [4]. Accordingly, it is not an easy task to estimate velocity for moving ships in the ocean.

Estimation methods for ship velocity have been extensively investigated in recent years. Methods can be roughly divided into two categories. The first category is based on ship wakes. Ship wakes have an important function in ship detection and can be used to measure ship velocity. In SAR images, a ship wake appears like a fan of bright line-shaped features around a central dark line [5]. Tunaley [6] exploited the cross-range separation between a ship’s location and its wake to estimate ship speed in images. Zilman et al. [7] observed Kelvin and turbulent wakes, and analysed the spectrum peak of ship-generated waves along the Kelvin cusp-lines to estimate ship speed. Copeland et al. [8] applied localized Radon transform to detect ship wakes in SAR. Courmontagne [9] exploited the stochastic match filtering technique to improve the detection performance of Radon transform. Carona and Marques [10] presented a strategy using the Radon transform to detect ship wakes and then estimate the range velocity component. Nevertheless, ship wakes are not always visible in SAR images. For example, with specific reference to RADARSAT-2 imagery, ship wakes associated with detected ships have been detected in less than 5% of cases [11].

The second category is based on Doppler parameter estimation because Doppler parameters are directly related to the motion parameters of a moving target. Compared with methods based on ship wakes, the methods in the second category have wider applications to ship velocity estimation. Renga and Moccia [12] proposed a ship velocity method for SAR images. The method measures the Doppler centroid of a ship candidate when the main feature of the ship is clearly imaged and the relevant complex pixels can be isolated from the background. Radius and Marques [13] proposed estimating the velocity vector of a moving ship in single-look complex (SLC) SAR data from antenna pattern effects induced by target motion. Dragosevic and Vachon [14] described a novel method for the radial component of a ship’s velocity using single-channel synthetic aperture radar acquired by RADARSAT-1, which is applicable in the absence of visible wakes. Kirscht [15] assessed the temporal correlation between sublook images, which can be used to track vessel motion. Time-frequency processing is a powerful tool for analysing the properties of a linear frequency-modulated (LFM) signal, and the fractional Fourier transform (FrFT) is the representative in the processing. Pelich et al. [16] estimated the azimuth velocity component in SAR images based on the FrFT. However, they neglected the effect of range velocity and range acceleration components. Huang et al. [17] utilized Radon-high-order time-chirp rate transforms to refocus on fast moving ground targets and estimate high-order motion parameters. The method concurrently addresses the range walk and Doppler frequency shift in the presence of Doppler ambiguity, which could be generalized to estimate ship velocity. However, in the case of low signal-to-clutter ratio (SCR), the estimation performance of low-order motion parameters degraded because of the error propagation effect. In addition to these two categories of methods, there are others. For example, Soccorsi and Lehner [18] applied change detection techniques to a time series to estimate ship speed, and Ao et al. [19] proposed a moving ship velocity estimation method based on subaperture decomposition. Filippo [20] applied pixel tracking and convex optimization to achieve micromotion and inclination angle estimation. Back et al. [21] used the Doppler parameter estimation to estimate two-dimensional vessel velocity.

To summarize, these methods are usually based on single-channel SAR. However, the performance of these algorithms is limited by the signal-to-clutter-plus-noise ratio (SCNR) [22]; it is difficult to separate the target from the clutter background. These methods can barely estimate multiple motion parameters simultaneously with high accuracy; often, they achieve only one parameter estimation, such as range velocity or azimuth velocity. With the advent of multichannel SAR systems, additional echo information can be obtained and clutter can be suppressed more; therefore, more attention has been attracted for moving target indication (MTI) in multichannel than in single-channel applications. Three techniques have dominated the literature over the last few decades, displaced phase centre antenna (DPCA) [22], along-track interferometry (ATI) [23], and space-time adaptive processing (STAP) [24]. This has led to a breakthrough in applications to ship detection and velocity estimation. For example, in [25], Dragosevic et al. combined adaptive beamforming and frequency tracking, with ATI to estimate ship speed. Jiang et al. [26] applied FrFT to processing and used DPCA ship signatures in the range-compressed domain to achieve their goals. Chiu [27] proposed a settlement based on FrFT and ATI in the range-compressed domain. Cerutti-Maori et al. generalized the DPCA and STAP processing from a dual-channel SAR system to any multichannel SAR system, termed extended DPCA (EDPCA) and imaging STAP (ISTAP) in [27, 28], respectively. These modified techniques can also be used to estimate ship speed.

The SAR moving target velocity estimation can be achieved using two means, for both single-channel and multichannel. The first class of algorithm is implemented on the raw data or range-compressed domain, and the second class is implemented on the focused SAR image. For estimating velocity in multichannel SAR, which is the objective in this study, the algorithm proposed in [27] for data or the range-compressed domain has excellent proven performance. However, operator-friendly image-type products are sometimes preferred [29]. This paper focuses on developing a ship velocity estimation algorithm for the image domain in the dual-channel ATI-SAR mode. The proposed algorithm works on processed images for dual-channel SAR systems rather than range-compressed domain, which could estimate multiple motion parameters simultaneously.

The remainder of this paper is organized as follows. In Section 2, we provide an overview on the moving target echo model and ATI phase. Section 3 introduces the definition and properties of FrFT and the proposed method that combines the ATI phase with FrFT to estimate the motion parameters of a moving ship. In Section 4, we compare our method with the method provided in [16] in the image domain. Finally, summaries and conclusions are provided in Section 5.

2. Overview on the Moving Target Echo Model and ATI Phase

ATI-SAR is a technique initially applied to measuring ocean currents [30] and has since been used to detect slow-moving objects, such as vehicles and ships [31], and estimate the velocity of moving targets and river [3235]. Obtaining velocity estimates of moving targets is an important goal of SAR-based moving target indication (MTI). Estimated outputs are often a combination of range, bearing, and Doppler frequency. The moving target echo model has a huge impact on velocity estimation; that is, the model accuracy has a great influence on algorithm accuracy, and the ATI phase has a proven ability to estimate radial velocity [35]. Therefore, in this section, we provide a brief overview of the echo model for moving targets and expression of the ATI phase used in our proposed method.

2.1. Moving Target Echo Model

The target and radar geometry for airborne SAR is illustrated in Figure 1, modified from [36]. In this figure, denotes the along-track direction, represents the across track direction, the initial coordinate of platform is , and the target is located at when , assuming the target moves with velocity components and acceleration components . The symbol indicates the platform velocity, and is the distance between the target and radar at , i.e., the nearest distance.

According to the geometry shown in Figure 1, we can obtain the range equation between the target and radar as [36]

In general, we apply the approximate expression of (1) using the Taylor series expression [37]; applying the second-order Taylor series expression, equation (1) can be rewritten as [36]

The approximation shown in the above expression is based on the following fact, i.e., the influence of the third and higher order terms of the Taylor series expression is relatively small because the values of acceleration and the rate of change of acceleration are too small for the distance between the radar platform and the target. The distance between the phase centres of two antennas in a dual-channel system is , and distance from the radar platform to the target is assumed to be large enough so that the far-field approximation can be applied [36].

The moving target echo model in terms of range history can be expressed as follows [36]:where is the magnitude of the fore channel and aft channel, respectively, is the range from the transmitting antenna to the moving target and back to the two receive antennas, is the wave number, and is a rectangle window centred at of length .The difference in the received echo between the fore and aft channel results from . For the fore channel, the range can be denoted by

On the basis of (4), the range of the aft channel can be expressed using a function of as

Furthermore, the expression of was derived in [38] as

In [39], the first-order Taylor series expression was used to approximate (6), that is,

For the any given sample time for the fore channel, the range to the target at time should be determined for the aft data to line up the channels. The registered aft channel can be expressed as

Substituting (2) and (8) into (4) and (5), the following equations can be obtained:

Substituting (9) and (10) into (3), we obtain detailed expressions of the moving target echo for the fore and aft channels; then, the established moving target echo model is used to calculate the ATI phase. The model characterized with formula (24) is an approximation of the real scenario. The more precise model should consist of motion in three directions, i.e., roll, yaw, and heave. Roll and yaw are considered in formula (24), but heave motion is neglected, for more complex models making it more difficult to retrieve parameters from the model. At present, it is not realistic to solve the task including six unknown parameters, i.e., velocity and acceleration in three directions. To the best of our knowledge, the model in [40, 41] has not been used in the published literature to estimate motion parameters in three directions simultaneously. Moreover, the velocity information for range and azimuth contained in the model shown in formula (24) might be sufficient for practical applications, such as ship location and refocusing.

2.2. ATI Phase

From these described derivations, we can obtain the ATI phase using conjugate multiplication. However, the ATI operator is often achieved in the azimuth-compressed domain [36]. The relationship between the original echo and azimuth-compressed echo is as follows:where is the compressed azimuth signal for the channel and is the reference function.

According to (11), we can derive the expression or peak response time for the fore channel and ATI phase as

The detailed derivation can be found in [36]. It should be noted that the fitting condition of the equation (12) is target’s azimuth coordinate is 0 when . The rationality of this fitting condition should be further assessed. However, considering that equation (12) has been used to the practical applications of velocity estimation in [36], we follow the same fitting condition in this paper.

Evaluating the expression of the ATI phase, the ATI phase is determined from three motion parameters of the moving target. That is, azimuth velocity , range velocity , and acceleration . However, in the other references, such as [38, 42, 43], equation (13) is modified as

Comparing (13) with (14), we find that in (14), the effect of and on the ATI phase is ignored, which is because one equation cannot solve three unknown parameters simultaneously. The error due to this approximation is acceptable in many cases. As discussed in [36], when and are nonzero, the overall shift in the ATI phase increases or decreases depending on their sign. In practice, we cannot guarantee the sign of the motion parameters, so replacing (13) with (14) is a reasonable choice. The above conclusions are derived for airborne SAR only, but it is easy to be extended to spaceborne SAR. We apply (13) for further discussion in the following sections.

3. The FrFT and the Proposed Method

In this section, we briefly introduce the definition and applications of FrFT and propose a method for estimating ship motion parameters.

3.1. FrFT

The FrFT has many applications in solving differential equations, quantum mechanics, and quantum optics, pattern recognition, and studying time-frequency distribution [44]. Chirp rate estimates from the LFM signal is one of the most significant applications for FrFT.

The FrFT represents a generalized form of the classical Fourier transform, and the FrFT definition is given as follows:where is the order of FrFT, corresponding to a rotation angle , and is the kernel function of the FrFT. The intact definition is

Applied to LFM signals, compared with the energy spread in the Fourier transform domain, the energy in the FrFT domain is more concentrated. The order of the FrFT facilitates finding the rotation that optimizes the representation in the time-frequency domain. Herein, optimization means that the target energy is the most concentrated.

3.2. Proposed Method

In the general linear chirp formula , the chirp rate is , and the optimum transform angle that best focuses on the fractional Fourier spectrum of the chirp signal is defined as [44]where is the sampling frequency and is the number of samples.

Adapting the equation for the azimuth LFM signal, (17) is replaced [44] aswhere is the pulse repetition frequency, is the number of samples in the azimuth, and represents the azimuth rate.

There is a relationship between phase and chirp rate:where is chirp rate and is echo phase.

Similar to the derivation in [16], on the basis of (9) and (10), we can obtain the chirp rate of the moving target echo after range compression for the fore channels aswhere is the chirp rate for the stationary scene.

After azimuth compression, the chirp rate expression for the fore channel can be obtained as follows:

Equation (21) is the theoretical chirp rate. Clearly, from (18) and (21), we can obtain the following equation:where is the optimum transform angle corresponding to the fore channel.

Combining (12), (13), and (22), we obtain an equation set relating , , and ay0. From this equation set, we can solve the three unknown motion parameters. Generally speaking, it is difficult to obtain analytic solutions from a nonlinear equation set. Fortunately, the equation set derived in our paper can be obtained using mathematical manipulation.

We use the following variable substitution:

Then, equations (12), (13), and (22) can be rewritten as follows:

First, from (26), we can easily obtain the value of according to equation (27), i.e.,

We can estimate from (24):

Next, we notice that can be resolved from (25) as (29), because the estimation of has been already calculated:

Finally, we obtain the estimation of from (23) and (26) as

The overall motion parameter estimation procedure is summarized in Figure 2, which involves four steps detailed as follows:Step 1. Generate the ATI image of the moving target: conjugate multiply the fore channel image with aft channel image to obtain the ATI image.Step 2. Obtain peak time and ATI phase: determine the peak time from the fore channel image on the basis of pixel magnitude. Then, retrieve the ATI phase from the maximum response in the ATI image.Step 3. Calculate the chirp rate: utilize global FrFT map to obtain the optimum rotation angle from the fore channel image (the strategy in [16] can be exploited to obtain the global FrFT map, see [16] for details). Then, calculate the corresponding chirp rate of the moving target echo.Step 4. Estimate the motion parameters: exploit formulae (28)–(30) to estimate the azimuth velocity, range velocity, and acceleration for the moving target.

4. Simulation Results

In this section, we compare our proposed method with the method presented in [16] and illustrate our method’s superiority. The method in [16] is implemented in image domain and the methods based on multichannel SAR are often achieved in RD domain [17]. In the general case, the SAR systems such as TanDEM-X could provide high-resolution complex image data [45], so the method implemented in image domain might have more practicability than the method achieved in RD domain. Therefore, in the scope of this study, we consider the velocity estimation realized in image domain rather than RD domain.

4.1. Method Description in [16]

The moving target echo model in the contrast method is very similar to that proposed in this study. For the single-channel SAR, the model must be simplified. Compared with (20), the influences of and on the chirp rate are ignored in [16]. Hence, the chirp rate of the moving target echo after the range compression for the fore channels is

Therefore, the chirp rate after azimuth compression can be expressed as

According to the above derivation, azimuth velocity can be estimated from (33) aswhere is the optimum transform angle corresponding to the single channel.

4.2. Simulation Experiments

We have described the principle of the method in [16], which generalized the FrFT from the signal domain to the image domain. As a result, we can compare the two methods in the image domain. Table 1 provides the simulation parameters for the ATI imaging; the simulation parameters are identical to the DC-8 platform in the [46]. The data collected from the system have been widely applied to ship detection [47, 48]. Figures 3(a) and 3(b) correspond to the simulated images of the fore and aft channels, respectively. X-axis of images represents the range, and y-axis represents azimuth.

From Figure 3, we can obtain the complex image data after conjugate multiplying. Then, the peak response time and ATI phase can be resolved from the ATI image. For the fore channel image, we extract the azimuth line containing the pixel with the greatest magnitude to calculate the azimuth chirp rate based on the FrFT.

Assuming the target moves along the azimuth line, our method can be reduced to the contrast method because the values of range velocity and acceleration are zero. Table 2 provides the estimation and Figure 4 presents the absolute error of the estimation, x-axis denotes the true azimuth velocity, and y-axis denotes the absolute errors between the true and estimated azimuth velocity. In Table 2, the percent error is calculated as follows:where is the percent error, is the estimation of azimuth velocity, and is the true value of azimuth velocity.

From Table 2 and Figure 4, the contrast method has a comparable estimation accuracy in the absence of range velocity and acceleration. However, in practice, the ship usually moves with an angle from the azimuth line. Thus, the range velocity and acceleration are nonzero, and the approximation in [16] may result in an additional error. Next, we consider the estimation results from our method and the contrast method in the presence of range velocity and acceleration. Table 3 provides the estimation results from our method and the contrast method; the range velocity and acceleration are focused on 4 m/s and 0.1 m/s2, respectively, and the azimuth velocity varies from 1 m/s to 15 m/s. Analogously, in Table 4, the azimuth velocity and range acceleration are focused on 4 m/s and 0.1 m/s2, respectively, and the range velocity varies from 1 m/s to 10 m/s. The azimuth and range velocity in Table 5 are simultaneously focused on 4 m/s, and the range acceleration varies from 0.01 m/s2 to 0.1 m/s2. The estimated results of different simulated scenario are shown in Figures 5 and 6. Figure 5(a) represents the absolute errors between the true and estimated azimuth velocity, Figure 5(b) represents the contrast trial between the methods in [16] and our proposed method, Figure 5(c) represents the absolute errors between the true and estimated range velocity, and Figure 5(d) represents the absolute errors between the true and estimated range acceleration.

As shown in Table 3 and Figure 5, our proposed method provides a more accurate estimation than the contrast method for azimuth velocity. The value estimated using the contrast method is biased from the truth above 2 m/s; as the azimuth velocity increases, the bias becomes larger. A similar observation can be made for the proposed method, but the bias is only 0.21 m/s when the azimuth velocity is 15 m/s. The difference is due to the moving target echo model. In comparison, the contrast method derived the azimuth velocity estimation expression by neglecting the range velocity and acceleration.

Comparing Table 2 with Table 3, the performance of the contrast method significantly degrades with range velocity and acceleration. In addition to the azimuth velocity, the proposed method can guarantee estimation accuracy for range velocity and acceleration.

From Table 4 and Figure 6, the influence of range velocity on azimuth velocity estimation can be analysed more clearly. When the range velocity is 1 m/s, the absolute error in the azimuth velocity estimation is about 0.4 m/s; as the range velocity increases to 10 m/s, the absolute error increases to 0.8 m/s. This observation further verifies the effect resulting from the range velocity.

Table 5 and Figure 7 indicate the influence of range acceleration on the azimuth velocity estimation. When the range acceleration is 0.01 m/s2, the absolute error is about 0.1 m/s, and as the range acceleration increases to 0.1 m/s2, the absolute error increases to 0.6 m/s. In contrast to range velocity, the range acceleration has a smaller impact on the azimuth velocity estimation. The chirp rate is likely dominated by the range velocity, except for the range acceleration has parameters shown in Table 1. Ignoring the constant, the chirp rate is dependent on the three motion parameters, as shown in Section 3:

Equation (35) indicates that the influence of azimuth velocity on chirp rate is multiplied by the platform velocity. Therefore, it may be appropriate in [16] to neglect the influence of range velocity and acceleration for the data being collected from Radarsat-2, whose SAR sensor velocity is 7545 m/s. However, for the airborne platform, the influence of range velocity and acceleration cannot be ignored in some conditions, especially when the value of acceleration is large, such as 0.1 m/s2 in Table 3.

All the tables show evident trends in the estimations calculated using the contrast method. It is clear that the chirp rate is a monotonic function of azimuth velocity (for small values), range velocity, and acceleration. In Tables 4 and 5, as the range velocity and acceleration, respectively, increase, the absolute errors in the azimuth velocity estimation also increase due to a lower estimation in the chirp rate. However, from Table 3, the absolute errors in the azimuth velocity are almost 0.6 m/s. Small fluctuations are probably caused by different FrFT estimation precisions. In contrast to the trends in the results from percentage of the azimuth velocity, the proposed method is not evident. This is because the proposed method requires measuring more quantities from the images, which leads to inevitable measurement errors, so the estimation fluctuates more and does not have evident trends.

5. Conclusion

In this study, we introduced a novel velocity estimation algorithm for moving ships based on the FrFT and ATI. Compared with the method in [16], the proposed method fully exploits the information advantages of multichannel SAR, combining the FrFT and ATI to construct equation sets that can achieve multiple motion parameter estimations. From the experiments in Section 4, the proposed method obtained a more accurate estimation for azimuth velocity compared to the contrast method, and the estimated range velocity and acceleration are relatively accurate.

Moreover, the chirp rates measured at different parts of a moving ship are not constant because they are modulated by the ship’s three-axis motion.

The proposed solution to this problem is similar to the contrast method. In [16], the authors applied the FrFT locally to segments extracted around maximum values from the input signal. The estimation accuracy is determined based on the constructed model. In this study, the model is the same as the commonly used model in [36] and a good performance is obtained. However, for a moving ship (or even a “stationary” ship), the heave motion must be included or modelled to be of practical use. In addition, even without including the heave motion components in the ship range history model, the acceleration component in the along-track direction introduces nonlinearity in the signal chirp. We have no reason to doubt that a more accurate estimation could be theoretically acquired using the model in [36]; however, this is challenging at present and should be the subject of future research. Therefore, in this study, we only consider azimuth velocity, range velocity, and acceleration. In this case, the received echo can be approximated well using the LFM signal. Moreover, our proposed method is not suitable for spaceborne SAR because of earth surface curvature and topography. There is a large deviation if the flat earth approximation is applied to the spaceborne SAR data [49].The multichannel SAR, which contains three or more channels, and nonlinearity and nonlinearity in the signal chirp, caused by azimuth acceleration, will be considered in future research to achieve accurate motion parameter estimations.

Data Availability

The SAR simulation data used to support the findings of this study are currently under embargo.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Conceptualization and funding acquisition were carried out by K.H. and C.C.; methodology was developed by C.C. and Y.L.; validation was done by Y.L., C.C., and L.Z.; formal analysis was conducted by Y.L. and K.O.; investigation was carried out by K.O., K.H., and C.C.; writing (original draft preparation) was performed by Y.L. and K.H.; writing (review and editing) was performed by C.C., K.H., and L.Z; supervision was provided by K.H; project administration was carried out by C.C.

Acknowledgments

This research was funded by the National Science Foundation of China under grant no. 71471174. This work was supported in part by the National Natural Science Foundation of China under Project 41822105, in part by Key Research Plan of Hunan Province under Project 2019SK2173, in part by the Field Foundation under Project 61404160109, and in part by the Fundamental Research Funds for the Central Universities under Project 2682020ZT34.