GRFT-Based Moving Ship Target Detection and Imaging in Geosynchronous SAR

Geosynchronous synthetic aperture radar (GEO SAR) has great potentials in ship surveillance due to its high time resolution and wide swath coverage. However, the remote slant range will result in a very low signal-to-noise ratio (SNR) of echoes that need to be enhanced by long-time coherent integration. The generalized Radon-Fourier transform (GRFT) can realize the coherent integration of moving target under long integration time by jointly parameter searching along range and velocity directions. Unfortunately, in GEO SAR, the very large slant range and long synthetic aperture will cause the curved synthetic aperture trajectory and non-negligible signal round-trip delay, leading to the failure of the traditional slant range and GRFT signal model for moving targets. This paper proposes an improved GRFT-based approach to realize the detection and imaging of moving ship targets in GEO SAR. Firstly, the accurate slant range for moving ship targets is constructed and the GRFT signal is redefined considering the curved trajectory and signal round-trip delay in GEO SAR. Then, GRFT responses to different motion parameters are analyzed. The procedures of moving ship targets detection and imaging in GEO SAR are presented through the detection with coarse-searched motion parameters in GRFT and the following imaging with fine-searched motion parameters based on minimum entropy. Finally, computer simulations verify the proposed GRFT-based method.

As is well known, ship surveillance needs high-frequent and large-scale monitoring of vast sea areas.Unfortunately, low earth orbit SAR (LEO SAR) has coarse time resolution (days to weeks) and limited swath coverage (tens to hundreds of kilometers).Although spaceborne multifunctional SAR has been applied in the detection of moving target (such as TerraSAR-X [10] and RADARSAT-2 [11]), achieving the monitoring of ships [12], the further application is limited by the revisiting time and swath width.
Geosynchronous SAR (GEO SAR) runs on the orbit with the orbit height of around 36,500 km [13], which can provide daily coverage for approximately one third of Earth surface.Areas of interest can be observed within the range of every 1.5-2.5 h per day [14].Thus, GEO SAR has great advantages in ship surveillance over vast ocean area.
Obviously, due to high time resolution and wide swath coverage, GEO SAR has great potential in monitoring moving ships.However, the remote slant range of GEO SAR will result in a very low signal-to-noise ratio (SNR) of echoes.The defect of low SNR can be compensated for the static target through long time integration [32].Unfortunately, the parameters of moving targets are unknown, leading to the failure of coherent integration and the resultantly SNR loss.Therefore, traditional detection and imaging methods for moving targets used in LEO SAR and airborne SAR are not suitable for GEO SAR.
Moving targets detection and imaging can be divided into single-channel and multi-channel SAR schemes according to SAR system configuration.Because of the limitation of low SNR, multichannel configuration cannot be realized currently, as it is mainly designed for clutter suppression where ships have higher radar cross section (RCS) against the surrounding sea background (in general, ships has higher backscattering energy as the metallic materials and corner reflection structures [33,34]).Therefore, moving ship detection and estimation can be realized without sea clutter suppression by single-channel SAR scheme, which can be divided into short-time and long-time integration.Short-time integration utilizes the difference of Doppler frequency distribution between moving targets and clutter to achieve moving targets detection [35], such as moving target indication (MTI) and moving target detection (MTD).This kind of methods only integrates signal energy in the same range unit, so it cannot be realized in moving targets detection in GEO SAR where the range migration is much severer even for static target and the energy is dispersed in tens of range units [36,37].Thus, long-time integration approaches should be employed.
Focus-before-detection (FBD) technique [38] is a kind of novel long-time integration method.It has been applied to maneuvering targets detection [39,40], targets detection under strong clutter background [41,42], estimation of ground target motion parameters [43], etc. FBD can greatly improve SNR by long-time coherent integration, thereby realizing the detection and estimation of weak targets.GRFT [44] is a typical FBD method, which can effectively overcome the coupling between the range walk and phase modulations by jointly searching along range and velocity directions.
However, in GEO SAR, the integration time is at 100 s level, so the trajectory cannot be approximated as a straight line anymore; meantime, its signal round-trip delay is in 0.1 s order, so the platform motion during the signal propagation cannot be ignored.Thus, the curved synthetic aperture trajectory and non-negligible signal round-trip delay will lead to the failure of the traditional slant range model for moving targets, which has been commonly used in LEO and Airborne SAR.The GRFT-based moving ship target detection and imaging method cannot be employed directly.
This paper proposes an improved GRFT-based approach to realize moving ship targets detection and imaging in GEO SAR.The main contributions are as follows: the accurate slant range model for moving ship targets in GEO SAR is established; then, based on the slant range model, the GRFT-based signal integration method is modified and the processing for moving ship targets detection and imaging is given.In the paper, we assume that the moving ship can be focused when there exists the rotation of the moving ship caused by ocean waves.This is under the assumption that the ship has huge inertia and the motion of hull centroid is not affected by ocean waves.The fine analysis of the ocean waves' influences will be considered in future work.
The paper is organized as follows.In Section 2, the accurate slant range model for moving ship targets in GEO SAR is established and then simplified by analyzing influences of target motion.Subsequently, the GRFT signal model in GEO SAR is modified.In Section 3, responses of GRFT for different motion parameters are analyzed, and then two kinds of motion parameters' searching intervals (including coarse and fine motion parameters' searching) are defined.Then, the minimum entropy is used as the moving ship target imaging criterion.The processing of moving ship targets detection and imaging in GEO SAR, which can be summarized as detection with motion parameters' coarse searching in GRFT and the following imaging with fine-searched motion parameters, is given.In Section 4, simulations are carried out to verify the proposed method.Section 5 concludes the paper.

GRFT Model for Moving Ship Targets in GEO SAR
In order to make equations more clearly in this section, the symbol list is given as Table 1.

Mathematical Symbols
The Mean of Symbol

Subscript n nth pulse transmitting moment Subscript MT
Moving ship target Subscript rm The echo after the range compression Superscript ˆSearching value Superscript T The matrix transposition c The velocity of light R MT,n The slant range from the satellite to moving ship targets The satellite velocity vector t The slow time

∆
The relative motion between GEO SAR and the target within signal round-trip delay r MT (t) The Taylor series expansion of slant range r 0 , k 1MT , k 2MT , k 3MT and k 4MT Zero-order to fourth-order coefficients of the Taylor series expansion to moving ship target ∆r 0MT , ∆k 1MT , ∆k 2MT and ∆k 3MT Zero-order to third-order coefficients of Taylor series expansion to moving ship target in the relative motion The unit vector of radial velocity r mt

Range position of the target at ACM a mt
Azimuth position of the target at ACM v r The radial velocity of target motion v a The azimuth velocity of target motion f drm The Doppler FM rate of moving ship targets f dr The Doppler FM rate of static targets T I The integration time T p The pulse duration of the transmitting signal B s The pulse bandwidth rmt,n r Range position searching grid ât,m a Azimuth position searching grid n r The number of range searching grid m a The number of azimuth searching grid ρ r

Range resolution ρ a
Azimuth resolution λ Wavelength

Accurate Slant Range Model
In GEO SAR, the movement of static target with the Earth rotation must be considered.Therefore, the moving target velocity can be expressed as the resultant velocity, which combines the Earth rotation velocity and the velocity of target motion.The moving ship target have a huge inertia, so here we assume its motion is uniform within the integration time [45,46].→ v MT0 is defined as the velocity of ship targets' motion.The slant range model of moving ship targets considering the curved synthetic aperture trajectory and non-negligible signal round-trip delay, which can be obtained by modifying the slant range model of static targets [47], can be written as Figure 1 shows the geometry of moving ship targets in GEO SAR.This model takes into account the Earth rotation.Thus, it can be considered that the resultant velocity of moving ship targets can be divided into

Accurate Slant Range Model
In GEO SAR, the movement of static target with the Earth rotation must be considered.Therefore, the moving target velocity can be expressed as the resultant velocity, which combines the Earth rotation velocity and the velocity of target motion.The moving ship target have a huge inertia, so here we assume its motion is uniform within the integration time [45,46].
is defined as the velocity of ship targets' motion.The slant range model of moving ship targets considering the curved synthetic aperture trajectory and non-negligible signal round-trip delay, which can be obtained by modifying the slant range model of static targets [47], can be written as Figure 1 shows the geometry of moving ship targets in GEO SAR.This model takes into account the Earth rotation.Thus, it can be considered that the resultant velocity of moving ship targets can be divided into and gn v  .The slant range of moving targets in GEO SAR can be approximated by the Taylor series expansion of the two part in right of (1) on the slow time [47] Figure 1.Geometry of moving ship targets and geosynchronous synthetic aperture radar (GEO SAR).
The slant range of moving targets in GEO SAR can be approximated by the Taylor series expansion of the two part in right of (1) on the slow time [47] r MT (t) ≈ (r 0 The influence of moving targets velocity on different coefficients is not the same.Accordingly, the influence is analyzed to simplify the coefficients of the slant range model in (2) (detailed deduction can be found in the Appendix B.1).The influence of the moving ship targets velocity on k 3MT and k 4MT , which are far less than π/8 and π/16, respectively, can be negligible.Therefore, k 3MT and k 4MT can be expressed in the same form as the static target.The Taylor series expansion of first part in (2) can be expressed as When the synthetic aperture time for GEO SAR imaging is at 100 s level, the influence of the moving ship targets velocity on the relative motion can be approximated as constant and it does not need compensation (see Appendix ?? in detail).So ∆r 0MT , ∆k 1MT , ∆k 2MT , and ∆k 3MT in (2) can be rewritten as ) Therefore, (2) can be rewritten as It is well known that the target motion influences the SAR azimuth focusing, due to the influence of the radial velocity and the azimuth velocity on the Doppler centroid and the Doppler FM rate, respectively.Therefore, two-dimensional (2D) velocity, Doppler centroid, and the Doppler FM rate of moving ship targets in GEO SAR are defined as following.
When GEO SAR in side-looking mode, the resultant velocity of moving ship targets can be expressed as is the radar platform velocity along the azimuth direction.Accordingly, the Doppler FM rate of moving ship targets can be expressed as Remote Sens. 2018, 10, 2002 The Doppler FM rate of static targets can be expressed as So, k 1MT and k 2MT can be rewritten as It can be known from ( 15) that k 1MT is determined by the radial velocity v r .At the same time, from (16), the Doppler FM rate f drm has an influence on k 2MT , which is determined by azimuth velocity → v a .That is, k 2MT is determined by v a .Therefore, it should be noted that only the influence of ship targets motion on k 1MT and k 2MT cannot be ignored.

GRFT Model for Moving Ship Targets
GRFT performs better than other signal integration methods in detection [40] and it can be approximately equal to the maximum likelihood estimation in parameter estimation [44].Although the echo SNR in GEO SAR is very low, GRFT can achieve the requirements of the weak target's signal integration.However, due to the characteristics of the slant range model, the existing GRFT cannot be applied to GEO SAR directly, resulting in the GRFT needing to be redefined to make it suitable for the signal integration of moving ship targets.
Based on the analysis in Section 2.1, the difference of the motion parameters will cause the slant range model to be changed.Therefore, it is necessary to define parametric slant range model for moving ship targets as where r 0 = r 2 mt + a 2 mt .Note that r mt , a mt , v r , and v a are unknown parameters that are caused by the motion of target, resulting in they are used as parameters for the parametric slant range model.
After the range compression, the received signal can be given as where ρ r = c/2B s , B s = γT p , γ is the frequency modulation rate of linear frequency modulated (LFM) signal, A rm = T p B s A r is the amplitude after the range compression, and A r is the amplitude before the range compression.Subsequently, the GRFT can be defined as s rm (t, r(t, rmt , âmt , vr , va )) exp j4πr(t, rmt , âmt , vr , va ) According to the previous analysis, the parameters in (17) have influence on r 0 , k 1MT , and k 2MT , and the coefficients, such as ∆r 0 , ∆k 1 , ∆k 2 , k 3 , ∆k 3 , and k 4 can be obtained from the platform directly and accurately.In this case, substitute (17), (18), into (19) GRFT s rm (r mt , âmt , vr , va ) = where A C = A rm exp(j4π(r 0 − r 0 )/λ).Based on ( 13), (15), and ( 16), (20) can be rewritten as GRFT s rm (r mt , âmt , vr , va ) ≈ (21), the difference of vr and va will cause the GRFT responses changed.

Procedure of the Moving Ship Target Detection and Imaging
The moving ship target detection and imaging procedure can be summarized, as in Figure 2.
Step 1: Coarse searching motion parameter for GRFT-based detection After the range compression, the target's 2D echoes are processing by (19) with the motion parameters' coarse searching interval.Motion parameters' searching interval is obtained by ( 23) and (24) in Section 3.3.2.The radial and azimuth velocity searching scope are vr ∈ v r,min v r,max and va ∈ v a,min v a,max , respectively.When the searching velocities have satisfied the focusing requirement (both | vr − v r | ≤ ρ r /T I and |v a − va | ≤ ρ a /2T I are shown in Section 3.3.1),the amplitude of GRFT response is greater than the CFAR detection threshold.From [48], the CFAR detection threshold can be obtained as When the GRFT response is greater than threshold, the radial and the azimuth searched velocity are vr_r and va , respectively.Unfortunately, due to the accuracy of the radial searched velocity is not enough, the azimuth shift occurs in GRFT response.So, it is necessary to relocate the moving ship target in imaging.
Step 2: Fine searching motion parameter for GRFT-based imaging Using GRFT to realize the moving ship target imaging with fine searching motion parameter, the radial velocity's searching interval is given by ( 25) in Section 3.3.2;the radial velocity's searching scope is vr_ f ∈ vr_r − ρ r /T I vr_r + ρ r /T I ; and, the azimuth searching velocity is set as va .The minimum entropy (26), which has obtained by Section 3.4, is used as criterion to determine whether the relocation requirement (| vr − v r | ≤ λ/2T I is shown in Section 3.3.1) is satisfied.When the entropy of GRFT response is minimum, the radial searched velocity is vr_ f .The GRFT response to vr_ f and va can be used as the moving ship target's image.
and, the azimuth searching velocity is set as ˆa v .The minimum entropy (26), which has obtained by Section 3.4, is used as criterion to determine whether the relocation requirement ( ˆ/ 2

GRFT-Based Signal Integration
The schematic diagram of GRFT-based signal integration is given in Figure 3. Generally, the echo of moving target has cross range cell migration.When the signal has been submerged in the noise, the traditional integration method, such as MTD, which integrates just in a range unit, cannot be used.From Figure 3, the echoes after range compression are distributed in multiple range units that can be described by ( ) , , , ,

GRFT-Based Signal Integration
The schematic diagram of GRFT-based signal integration is given in Figure 3. Generally, the echo of moving target has cross range cell migration.When the signal has been submerged in the noise, the traditional integration method, such as MTD, which integrates just in a range unit, cannot be used.From Figure 3, the echoes after range compression are distributed in multiple range units that can be described by r(t, r mt , a mt , v r , v a ).GRFT use (19) to integrate signal by searching rmt , âmt , vr , va .That is, s rm (t, r(t, rmt , âmt , vr , va )) represents the signal in r(t, rmt , âmt , vr , va ) at slow time t and exp(j4πr(t, rmt , âmt , vr , va )/λ) is used to compensate phase fluctuation.At the same time, the integral realize the signal integration in whole T I .Since r(t, rmt , âmt , vr , va ) is close enough to r(t, r mt , a mt , v r , v a ), the signal will be integrated in the corresponding range and azimuth position searching grid, such as r mt and a mt , where the moving target located in ACM.Therefore, GRFT can realize coherent integration in T I , which makes it superior than the traditional integration method.
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 27 searching grid, such as mt r and mt a , where the moving target located in ACM.Therefore, GRFT can realize coherent integration in I T , which makes it superior than the traditional integration method.The GRFT-based signal integration processing for moving ship targets is given, as follows: (1) Determine range-azimuth 2D position and radial-azimuth 2D velocity searching grid according to requirements of the resolution and GRFT responses to different motion parameters.(2) Based on 2D position and 2D velocity searching grid, signal integration can be realized by (19).( 3) Each of 2D velocity searching grid produces a GRFT response, which is the integration gain in 2D position searching grid.
There are two key issues in signal integration processing, as follows.Firstly, the 2D position searching grid is divided based on resolution.According to the SAR principle [1], the range and the azimuth resolution are / 2 , respectively.The position of moving targets at ACM just needs to be located in the resolution cell.Therefore, the range position searching grid is and the azimuth position searching grid is ˆt a a mρ ∈ ± .
Secondly, 2D velocity searching grid needs to be divided.The radial and azimuth velocity The GRFT-based signal integration processing for moving ship targets is given, as follows: (1) Determine range-azimuth 2D position and radial-azimuth 2D velocity searching grid according to requirements of the resolution and GRFT responses to different motion parameters.(2) Based on 2D position and 2D velocity searching grid, signal integration can be realized by (19).(3) Each of 2D velocity searching grid produces a GRFT response, which is the integration gain in 2D position searching grid.
There are two key issues in signal integration processing, as follows.Firstly, the 2D position searching grid is divided based on resolution.According to the SAR principle [1], the range and the azimuth resolution are ρ r = c/2B s and ρ a = λ → r s0 − → r MT0 /2v s T I , respectively.The position of moving targets at ACM just needs to be located in the resolution cell.Therefore, the range position searching grid is rt ∈ r 0c ± nρ r and the azimuth position searching grid is ât ∈ ±mρ a .
Secondly, 2D velocity searching grid needs to be divided.The radial and azimuth velocity searching scope are vr ∈ v r,min v r,max and va ∈ v a,min v a,max , respectively.Unfortunately, the 2D velocity searching grid dividing requires the searching interval.Therefore, the GRFT response to different motion parameters must be studied to obtain 2D velocity searching interval.In GRFT, the parameters' searching interval is very important.When the parameters' searching interval is too large, GRFT cannot be focused, leading to detection failure [49].On the contrary, the calculation is too huge.Therefore, the GRFT response to different searching parameters needs to be studied.
Since the searching parameters, such as rmt , âmt , vr , and va are different with the motion parameters, there are several kinds of GRFT response, such as moving targets completely focused, azimuth shift, completely unfocused, and azimuth defocused.According to (21), three different cases that are based on different radial and azimuth searching velocity are analyzed in Appendix C.
The analysis results can be given, as following.Since | vr − v r | > ρ r /T I , the GRFT response will be approximately zero.Since both radial and azimuth velocity are searched inaccurately, the azimuth shift and the azimuth defocus in GRFT response equal to 2( vr − v r )T I /λ and 2(v a − va )T I /ρ a , respectively.In this case, since | vr − v r | ≤ λ/2T I and |v a − va | > ρ a /2T I , the GRFT response of moving target is located accurately and the azimuth defocus is equal to 2(v a − va )T I /ρ a .When λ/2T < | vr − v r | ≤ ρ r /T I and |v a − va | ≤ ρ a /2T I , the GRFT response of the moving target is focused and the azimuth shift equals to 2( vr − v r )T I /λ.Since the 2D velocity is searched accurately, the moving target can be completely focused where it is located at ACM.The GRFT response to different searching motion parameter can be summarized in Table 2.

Motion Parameters' Searching Interval
Unfortunately, only the completely focused requirements are used as the searching interval, resulting in the searching computation is too large.From Table 2, when the 2D searching velocity does not satisfy | vr − v r | ≤ ρ r /T I or |v a − va | ≤ ρ a /2T I , the GRFT response will become completely unfocused or azimuth defocused, which leads to its amplitude decrease.Due to the constant false-alarm rate (CFAR) detection being sensitive to signal SNR, it can be used to determine whether the 2D searching velocity satisfies above requirements.
The moving ship target relocation is also an important task of ship surveillance.Therefore, the accuracy of motion parameters needs to satisfy not only the focusing requirements (| vr − v r | ≤ ρ r /T I and |v a − va | ≤ ρ a /2T I ), but also the relocation requirements (| vr − v r | ≤ λ/2T I ).According to ρ r λ, the coarse and the fine searching interval are proposed to satisfy the focusing and the relocation requirements, respectively.In this case, the computation of the GRFT processing can be greatly reduced.Two kinds of motion parameters' searching interval are presented, as follows.
1. Motion parameters' coarse searching interval (1) Radial velocity's coarse searching interval Form Table 2, since | vr − v r | ≤ ρ r /T I is not satisfied, the GRFT response will be approximately zero.The radial velocity's coarse searching interval can be expressed as (2) Azimuth velocity searching interval When |v a − va | > ρ a /2T I , the target is defocused along azimuth direction and the azimuth defocus equals to 2(v a − va )T I /ρ a .Therefore, the azimuth velocity searching interval can be expressed as The 2D velocity searching grid is divided by ( 23) and ( 24), and then the moving target signal is integrated using (19).When the 2D searching velocity satisfies both | vr − v r | ≤ ρ r /T I and |v a − va | ≤ ρ a /2T I , the moving targets are focused without azimuth defocus, which can realize target detection.At this time, the radial and azimuth searching velocity of GRFT are the motion parameters searched value in coarse searching.
However, since | vr − v r | ≤ λ/2T I is not satisfied, the azimuth shift appears in the GRFT response.Note that the relative position between the static targets is not changed, that is the static targets have the same azimuth shift.Unfortunately, the azimuth shift, which is related to the motion between targets and platform, is not the same between the static and the moving targets.Therefore, the moving targets imaging needs relocation and a fine searching of the radial velocity is required.

Motion parameters' fine searching interval
When | vr − v r | ≤ λ/2T I is not satisfied, the GRFT response is shifted along the azimuth.Consequently, the radial velocity's fine searching interval can be obtained as GRFT is a kind of FBD-based methods, that is the GRFT response can be used as the imaging of moving ship targets, since is satisfied by coarse-searched motion parameter.Unfortunately, the radial velocity's fine searching will obtain multiple GRFT responses for a moving target, resulting in CFAR criteria that cannot be used to determine whether the searching velocity satisfies | vr − v r | ≤ λ/2T I .At the same time, the ship target is a kind of extended target that cannot be estimated by finding the maximum value within a resolution cell.Therefore, in Section 3.4, the minimum entropy will be used as the criterion to realize moving targets imaging.

Imaging Based on the Minimum Entropy
The imaging of moving ship targets need to satisfy both focusing and relocation requirements.However, CFAR criterion cannot determine whether the relocation requirement is satisfied, that is, other kind of criterion is required to the motion parameters' fine searching.
The entropy is a measure of energy divergence [50].It is well known that when the motion searching parameter is closer to the true value, the signal integration is focused better and the entropy is smaller.Thus, the minimum entropy can be the criterion to determine whether the radial searching velocity is closest to the true velocity.From (21), the minimum entropy for GRFT can be obtained as |GRFT srm (r mt,nr , âmt,ma , vr , va )| 2 ln|GRFT srm (r mt,nr , âmt,ma , vr , va )| 2 (26) where |GRFT s rm (r t,n r , ât,m a , vr , va )| 2 , va is the azimuth searched velocity obtained by the coarse-searched motion parameters.Because the entropy calculates the energy integration in the region of interest, the motion parameters estimation of the extended targets can be realized and the moving ship target can be relocated in imaging.

Numerical Experiments and Performance Analysis
In this section, numerical experiments are carried out to verify the effectiveness of the proposed GRFT-based detection and imaging methods for moving ship targets in GEO SAR.In Section 4.1, the GRFT responses is verified, which proposed in Section 3.3.Subsequently, the 2D raw echoes of a moving ship target are generated in apogee of GEO SAR to verify the detection and imaging methods that are proposed by Section 3.1.In Section 4.2, the GRFT-based moving ship targets detection is obtained.In Section 4.3, the GRFT-based imaging with the minimum entropy and evaluation of imaging are given.The system parameters of GEO SAR for the simulation are given as Table 3.

GRFT Responses
In Section 3.3, GRFT responses to different searching motion parameters are analyzed.In this section, the simulation is obtained to verify the validity of GRFT responses.
Suppose that GEO SAR is near apogee; the moving ship target is located in the center of the scene at ACM; the radial velocity is v r = 10 m/s; and, the azimuth velocity is v a = 10 m/s.In order to analyze the experimental results easily, the ship target is represented by the point target.The GRFT responses to three kinds of motion parameters' searching interval are simulated, as follows.
(1) GRFT response to radial velocity's coarse searching According to (23), the radial velocity's coarse searching grid can be divided as vr = 9 + nρ v_r (m/s), n = 0, ±1, ±2; the azimuth searching velocity is set as va = 10 m/s.The GRFT response can be obtained by (19).The simulation results are shown in Figure 4.The GRFT responses to the radial velocity's fine searching are shown in Figure 5.The azimuth profile of GRFT responses with radial velocity's coarse searching is shown in Figure 4.The azimuth profile of GRFT responses with all coarse-searched radial velocities are shown in Figure 4a.When | vr − v r | ≤ ρ r /T I is satisfied, the GRFT response is much larger than other search velocities.When the radial searching velocities do not satisfy | vr − v r | ≤ ρ r /T I , the enlarged view of GRFT responses are shown in Figure 4b.
From the simulation results, it can be seen that the moving ship target can be focused in an azimuth resolution cell, since | vr − v r | ≤ ρ r /T I .Otherwise, the GRFT response is very weak.
(2) GRFT response to radial velocity's fine searching Since | vr − v r | ≤ λ/2T I is not satisfied, the GRFT response is shifted along azimuth and the azimuth shift equal to 2( vr − v r )T I /λ.Based on (25), the radial velocity's fine searching grid can be set as vr = 10 + nρ v_ f (m/s), n = 0, ±1, ±2, ±3; the azimuth searching velocity is set as va = 10 m/s.The GRFT responses to the radial velocity's fine searching are shown in Figure 5.The GRFT responses to the radial velocity's fine searching are shown in Figure 5.  (3) GRFT response to azimuth velocity's searching Since |v a − va | ≤ ρ a /2T I , the GRFT response of moving ship targets do not defocus along azimuth.Otherwise, the azimuth defocus is equal to 2(v a − va )T I /ρ a .Based on (24), the azimuth velocity searching grid can be divided as va = nρ v_a (m/s), n = 0, ±1, ±2, ±3; the radial searching velocity is set as vr = 10 m/s.The simulation results are shown in Figure 6.The above three cases are verified the GRFT response to radial velocity's coarse searching, fine searching, and azimuth velocity's searching, respectively.Accordingly, the simulation results are consistent with the theoretical analysis results that are shown in Table 2.

Simulation Results of the GRFT-Based Detection
The moving ship target is represented by five point targets and is located in the center of the scene at ACM; the radial velocity is 10 m/s due to the maximum velocity of the ship is 33 knots (about 16.98 m/s); the 2D searching velocity grid is divided according to (23) and (24).The GRFT response with coarse-searched motion parameters is simulated by (19).When the 2D searching parameters are equal to _ ˆ9.9228 r r v = m/s and ˆ9.8911 a v = m/s, the corresponding GRFT response is greater than (22).The corresponding GRFT response is shown as Figure 8.The 'red cross' in the center of scene represents the location of target at ACM.From simulation results, although the accuracy of 2D searched velocity satisfy the focusing requirement, the GRFT response cannot be used as the imaging of moving ship targets due to the azimuth shift.Therefore, it is necessary to relocate the target.From Figure 6, since |v a − va | ≤ ρ a /2T I ( va = 11.3162m/s), the target is focused in a resolution cell and QPE ≤ π/4.When |v a − va | ≤ ρ a /2T I is not satisfied, the amplitude of the GRFT response is decreased due to being partially focused.The azimuth defocus is equal to 2(v a − va )T I /ρ a .
The above three cases are verified the GRFT response to radial velocity's coarse searching, fine searching, and azimuth velocity's searching, respectively.Accordingly, the simulation results are consistent with the theoretical analysis results that are shown in Table 2.

Simulation Results of the GRFT-Based Detection
The moving ship target is represented by five point targets and is located in the center of the scene at ACM; the radial velocity is v r = 10 m/s; and, the azimuth velocity is v a = 10 m/s.We generate the echoes of the moving ship target (a target with five points) into the background raw data of GEO SAR.After the range compression, the received signal of moving ship target is completely submerged in the noise.The simulation results are shown in Figure 7.
The range of the radial and the azimuth searching velocity is −v max v max ; where v max = 17 m/s due to the maximum velocity of the ship is 33 knots (about 16.98 m/s); the 2D searching velocity grid is divided according to (23) and (24).The GRFT response with coarse-searched motion parameters is simulated by (19).When the 2D searching parameters are equal to vr_r = 9.9228 m/s and va = 9.8911 m/s, the corresponding GRFT response is greater than (22).The corresponding GRFT response is shown as Figure 8.The 'red cross' in the center of scene represents the location of target at ACM.From simulation results, although the accuracy of 2D searched velocity satisfy the focusing requirement, the GRFT response cannot be used as the imaging of moving ship targets due to the azimuth shift.Therefore, it is necessary to relocate the target.Due to the orbital height of GEO SAR varies with the orbital position, we set the scope of slant range ( MT R ) in all orbit positions as 33,000 km-40,000 km.From [35], the signal SNR after range compression can be given as ( ) 10 log 4 where k is Boltzmann constant and s T is receiver temperature.Therefore, the corresponding scope of rm SNR is 21.3 − dB-24.6 − dB.After signal integration by GRFT, the SNR of GRFT response will become 10.12 dB-13.46dB and 20.12 dB-23.45dB, since the integration times are equal to 10 s and 100 s, respectively.In order to prove the superiority of proposed integration method in GEO SAR, MTD, and Radon transform (RT) are used as comparisons that can be seem as the typical method of short-time integration and long-time incoherent integration, respectively.After signal integration by RT, the SNR of output will become 4.64 dB-8.08 dB and 9.72 dB-13.06dB, since the integration time equal to 10 s and 100 s, respectively.After processing by MTD with any integration time, the MTD output has been submerged in noise.
Furthermore, the detection probability of MTD, RT and GRFT can be calculated with a given constant false alarm ratio , which are shown as Figure 9. From the simulation results, when the integration time equal to 10 s, the GRFT can obtain detection probability more than 80% around 33, 000 km and the other two methods cannot be used in target detection.When the integration time equal to 100 s, the GRFT can realize high performance detection in any orbital position.At the same time, RT can realize detection around 33, 000 km and MTD cannot realize detection.Because MTD can integrate signal in just in a single range unit, it cannot realize detection in GEO SAR.Although Due to the orbital height of GEO SAR varies with the orbital position, we set the scope of slant range (R MT ) in all orbit positions as 33,000 km-40,000 km.From [35], the signal SNR after range compression can be given as where k is Boltzmann constant and T s is receiver temperature.Therefore, the corresponding scope of SNR rm is −21.3 dB-−24.6 dB.After signal integration by GRFT, the SNR of GRFT response will become 10.12 dB-13.46dB and 20.12 dB-23.45dB, since the integration times are equal to 10 s and 100 s, respectively.
In order to prove the superiority of proposed integration method in GEO SAR, MTD, and Radon transform (RT) are used as comparisons that can be seem as the typical method of short-time integration and long-time incoherent integration, respectively.After signal integration by RT, the SNR of output will become 4.64 dB-8.08 dB and 9.72 dB-13.06dB, since the integration time equal to 10 s and 100 s, respectively.After processing by MTD with any integration time, the MTD output has been submerged in noise.
Furthermore, the detection probability of MTD, RT and GRFT can be calculated with a given constant false alarm ratio P f a = 10 −6 , which are shown as Figure 9. From the simulation results, when the integration time equal to 10 s, the GRFT can obtain detection probability more than 80% around 33, 000 km and the other two methods cannot be used in target detection.When the integration time equal to 100 s, the GRFT can realize high performance detection in any orbital position.At the same time, RT can realize detection around 33, 000 km and MTD cannot realize detection.Because MTD can integrate signal in just in a single range unit, it cannot realize detection in GEO SAR.Although RT realizes signal integration across range unit, it does not compensate phase fluctuation.Therefore, GRFT is superior to RT and MTD in GEO SAR.
RT realizes signal integration across range unit, it does not compensate phase fluctuation.Therefore, GRFT is superior to RT and MTD in GEO SAR.

Simulation Results of the GRFT-Based Imaging
The GRFT-based imaging with fine-searched motion parameters is processed after the detection.
Where the range of fine-searched radial velocity is  4.Meanwhile, the evaluation of each point targets in corresponding GRFT response are given in Table 5.  12. From the simulation results, since the radial searching velocity do not satisfy the relocation requirement, the azimuth shift is obtained and the PSLR of GRFT response is not much different.Therefore, due to the entropy simulation results in this section, the minimum entropy can realize the relocation of the moving ship target, which can be used as the criterion of the GRFT-based imaging.

Simulation Results of the GRFT-Based Imaging
The GRFT-based imaging with fine-searched motion parameters is processed after the detection.
Where the range of fine-searched radial velocity is vr_ f ∈ vr_r − ρ r /T I vr_r + ρ r /T I = 9.8394 10.0051 (m/s); the radial searching velocity grid is divided by (25); and, the azimuth searching velocity is va = 9.8911 m/s.The entropy of the GRFT response to fine-searched motion parameters can be simulated by (26) and the simulation results are shown as Figure 10a.Since vr_r = 10.0002m/s, the entropy of GRFT response is minimum and the corresponding GRFT response is shown as Figure 10b.It can be seen from Figure 10b that, since the entropy of GRFT response is minimum, the accuracy of radial searched velocity can realize the moving ship target relocation, which satisfy | vr − v r | ≤ λ/2T I .The parameters of the moving ship target are summarized in Table 4.Meanwhile, the evaluation of each point targets in corresponding GRFT response are given in Table 5. Figure 11 is obtained PSLR of point target imaging in range and azimuth directions.
To verify the necessity of minimum entropy, we set the radial searching velocity as vr_r = 10.0002+ n•ρ rv_ f (m/s), n = ±5, ±10, ±15.The corresponding GRFT responses and azimuth profiles are shown in Figure 12.From the simulation results, since the radial searching velocity do not satisfy the relocation requirement, the azimuth shift is obtained and the PSLR of GRFT response is not much different.Therefore, due to the entropy simulation results in this section, the minimum entropy can realize the relocation of the moving ship target, which can be used as the criterion of the GRFT-based imaging.

Conclusions
Due to short revisiting time and large wide coverage, GEO SAR has a greater researching value than other platforms in the moving ship target surveillance.However, some problems have not been resolved, such as invalid traditional range model and weak signal echo.In this paper, the proposed accurate slant range has considered the characteristics of curved trajectory and non-negligible signal round-trip delay; the redefined GRFT model has resolved the problem of moving ship targets' signal integration in GEO SAR; the proposed processing realize detection, relocation, and imaging of the moving ship target.
In order to solve the problems of the moving ship target detection and imaging in GEO SAR.Firstly, from the analysis results, only the influence of approximately uniform moving ship target velocity on first-order and second-order coefficients cannot be negligible, and then the accurate slant range with four-order coefficients is given.Secondly, the GRFT model is redefined for moving ship targets.For the sake of realizing the signal integration, GRFT responses, such as completely unfocused, azimuth defocused, azimuth shift, and completely focused are analyzed.Resulting in two kinds of motion parameters' searching interval, such as coarse and fine searching interval, are given to realize the focusing (overcoming completely unfocused and azimuth defocused) and relocation (overcoming azimuth shift and realizing completely focused) of moving ship targets in the GRFT

Conclusions
Due to short revisiting time and large wide coverage, GEO SAR has a greater researching value than other platforms in the moving ship target surveillance.However, some problems have not been resolved, such as invalid traditional range model and weak signal echo.In this paper, the proposed accurate slant range has considered the characteristics of curved trajectory and non-negligible signal round-trip delay; the redefined GRFT model has resolved the problem of moving ship targets' signal integration in GEO SAR; the proposed processing realize detection, relocation, and imaging of the moving ship target.
In order to solve the problems of the moving ship target detection and imaging in GEO SAR.Firstly, from the analysis results, only the influence of approximately uniform moving ship target velocity on first-order and second-order coefficients cannot be negligible, and then the accurate slant range with four-order coefficients is given.Secondly, the GRFT model is redefined for moving ship targets.For the sake of realizing the signal integration, GRFT responses, such as completely unfocused, azimuth defocused, azimuth shift, and completely focused are analyzed.Resulting in two kinds of motion parameters' searching interval, such as coarse and fine searching interval, are given to realize the focusing (overcoming completely unfocused and azimuth defocused) and relocation (overcoming azimuth shift and realizing completely focused) of moving ship targets in the GRFT response, respectively.Because GRFT is a kind of the FBD method, when both focusing and relocation requirement are satisfied, the corresponding GRFT response can be used for the imaging of moving targets.Thirdly, the processing of detection and imaging is proposed.It can be summarized as the detection with coarse-searched parameters and imaging with fine-searched parameters which satisfy requirements of focusing and relocation, respectively.
Finally, numerical experiments and scene simulations are provided to demonstrate the effectiveness of the proposed methods in GEO SAR with 53 • orbit inclination.From the simulation results, the proposed method can realize the moving ship targets imaging with azimuth PSLR and ISLR of about −13.27 dB and −9.77 dB, respectively.At the same time, the proposed imaging method can satisfy the relocation requirement, resulting in the moving ship target can be accurately located in an azimuth resolution cell of SAR imaging.From Figure A1, it can be seen that 3mt k is far less than / 8 π , resulting in it can be negligible.
Then The influence on 4 MT k in (A5) can be defined as ( )      3 is 300 s.
The simulation results of (A16) are shown as Figure A3.From the analysis and experimental results, it can be seen that the influence of moving ship target motion on the relative motion can be approximated as a constant and do not need compensation.The following conclusions can be obtained through analysis and simulation.
1.The influence of the moving ship target velocity on 3MT k and 4 MT k can be negligible.
2. The influence of the moving ship target velocity on the relative motion caused by signal roundtrip delay can be approximated as a constant and do not need compensation.

Appendix C. GRFT Response to Different Searching Motion Parameters
According to (21), three different cases based on different radial and azimuth searching velocity are discussed as follows Case 1: when ˆr The following conclusions can be obtained through analysis and simulation. 1.
The influence of the moving ship target velocity on k 3MT and k 4MT can be negligible.

2.
The influence of the moving ship target velocity on the relative motion caused by signal round-trip delay can be approximated as a constant and do not need compensation.Therefore, when |v a − va | > ρ a /2T I , the GRFT response of moving target is defocused along azimuth and the azimuth defocus equal to 2(v a − va )T I /ρ a .
velocity vector of target motion (uniform) → v gn The velocity vector of target with Earth rotation → v sn

k 3 and k 4
First-order to fourth-order coefficients of the Taylor series expansion to static target ∆r 0 , ∆k 1 , ∆k 2 and ∆k 3Zero-order to third-order coefficients of Taylor series expansion to static target in the relative motion → r MT0 Moving target position vector at the aperture center moment (ACM) to fourth-order velocity vectors of the satellite at ACM → to fourth-order velocity vectors of the static target with Earth rotation at ACM r(t, r mt , a mt , v r , v a )Parametric slant range model→ u T sMT

Figure 1 .
Figure 1.Geometry of moving ship targets and geosynchronous synthetic aperture radar (GEO SAR).
shown in Section 3.3.1) is satisfied.When the entropy of GRFT response is minimum, the radial searched velocity is _ ˆr f v .The GRFT response to _ ˆr f v and ˆa v can be used as the moving ship target's image.

Figure 2 .
Figure 2. Flowchart of the moving ship target detection and imaging.
a v v .GRFT use(19) to integrate signal by searching ˆmt r , ˆmt a , ˆr v , ˆa v .r a v v π λ is used to compensate phase fluctuation.At the same time, the integral realize the signal integration in whole a v v , the signal will be integrated in the corresponding range and azimuth position

Figure 2 .
Figure 2. Flowchart of the moving ship target detection and imaging.

3. 3 .
Influence of Different Motion Parameters on Signal Integration 3.3.1.GRFT Response to Different Searching Motion Parameters

Figure 4 .Figure
Figure 4. Azimuth profile of GRFT responses to radial velocity's coarse searching: (a) all coarsesearched radial velocities (the pink line represents ˆ/ r r r I v v T ρ − ≤

( 2 )
GRFT response to radial velocity'not satisfied, the GRFT response is shifted along azimuth and the azimuth shift equal to ( ) on (25), the radial velocity's fine searching grid can be n = ± ± ± ; the azimuth searching velocity is set as ˆ10 a v = m/s.

Figure 5 .
Figure 5. Azimuth profile of GRFT responses to radial velocity's fine searching.

Figure 4 .
Figure 4. Azimuth profile of GRFT responses to radial velocity's coarse searching: (a) all coarse-searched radial velocities (the pink line represents | vr − v r | ≤ ρ r /T I ); (b) enlarged view of GRFT response when | vr − v r | > ρ r /T I .

Figure 4 .
Figure 4. Azimuth profile of GRFT responses to radial velocity's coarse searching: (a) all coarsesearched radial velocities (the pink line represents ˆ/ r r r I v v T ρ − ≤

( 2 )
GRFT response to radial velocity's fine searchingSince not satisfied, the GRFT response is shifted along azimuth and the azimuth shift equal to () n = ± ± ± ; the azimuth searching velocity is set as ˆ10 a v = m/s.

Figure 5 .
Figure 5. Azimuth profile of GRFT responses to radial velocity's fine searching.

Figure 5 .
Figure 5. Azimuth profile of GRFT responses to radial velocity's fine searching.It can be seen from simulation results that since | vr − v r | ≤ λ/2T I ( vr = 10 m/s) the moving ship targets is located accurately.On the contrary, the azimuth shift equal to 2( vr − v r )T I /λ.

( 3 )
GRFT response to azimuth velocity's response of moving ship targets do not defocus along azimuth.Otherwise, the azimuth defocus is equal to () n = ± ± ± ; the radial searching velocity is set as ˆ10 r v = m/s.The simulation results are shown in Figure6.

Figure 6 .From Figure 6
Figure 6.Azimuth profile of GRFT responses to azimuth velocity searching.
generate the echoes of the moving ship target (a target with five points) into the background raw data of GEO SAR.After the range compression, the received signal of moving ship target is completely submerged in the noise.The simulation results are shown in Figure7.The range of the radial and the azimuth searching velocity is max max

Figure 6 .
Figure 6.Azimuth profile of GRFT responses to azimuth velocity searching.

27 Figure 7 .
Figure 7. Outputs of signal data after the range compression.

Figure 7 .
Figure 7. Outputs of signal data after the range compression.

Figure 7 .
Figure 7. Outputs of signal data after the range compression.

Figure 8 .
Figure 8.The corresponding GRFT response for target detection.

Figure 8 .
Figure 8.The corresponding GRFT response for target detection.

Figure 9 .
Figure 9.The detection probability for different integration methods: (a) when integration time equals to 10 s; and, (b) when integration time equals to 100 s.
/s); the radial searching velocity grid is divided by(25); and, the azimuth searching velocity is ˆ9.8911 a v = m/s.The entropy of the GRFT response to fine-searched motion parameters can be simulated by(26) and the simulation results are shown as Figure10a.Since _ ˆ10.0002 r r v = m/s, the entropy of GRFT response is minimum and the corresponding GRFT response is shown as Figure 10b.It can be seen from Figure 10b that, since the entropy of GRFT response is minimum, the accuracy of radial searched velocity can realize the moving ship target relocation, which satisfy ˆ/ 2 of the moving ship target are summarized in Table Figure 11  is obtained PSLR of point target imaging in range and azimuth directions.To verify the necessity of minimum entropy, we set the radial searching velocity as n = ± ± ± .The corresponding GRFT responses and azimuth profiles are shown in Figure

Figure 9 .
Figure 9.The detection probability for different integration methods: (a) when integration time equals to 10 s; and, (b) when integration time equals to 100 s.

Figure 12 .
Figure 12.GRFT responses and azimuth profiles to different fine-searched radial velocities.

Figure A1 .
Figure A1.The influence of Target Motion on 3MT k

Figure A1 .
Figure A1.The influence of Target Motion on k 3MT of slant range.

k and 2 k
are first-order and second-order coefficients in slant range of the static target, ) can be negligible.The simulation results is shown as FigureA2.4mt k is far less than / 16 π , that is the influence can be negligible.

Figure A2 .
Figure A2.The influence of target motion on 4 MT k of slant range.Figure A2.The influence of target motion on k 4MT of slant range.

Figure A2 .
The influence of target motion on k 4MT of slant range.
maximum integration time of GEO SAR with parameters shown in Table

Figure A3 .
Figure A3.The influence of target motion on the relative motion caused by signal round-trip delay.

Figure A3 .
Figure A3.The influence of target motion on the relative motion caused by signal round-trip delay.

Table 4 .
The parameters of the moving ship target.

Table 5 .
Evaluation of GRFT-based point target imaging.

Table 4 .
The parameters of the moving ship target.

Table 5 .
Evaluation of GRFT-based point target imaging.