Azimuth Multichannel Reconstruction for Moving Targets in Geosynchronous Spaceborne–Airborne Bistatic SAR

: In a multichannel geosynchronous spaceborne–airborne bistatic synthetic aperture radar (GEO ‐ SA ‐ BiSAR) system, the airborne receiver can obtain high ‐ resolution microwave images with good signal ‐ to ‐ noise ratios (SNRs) by passively receiving echoes from the desired area. Since the Doppler modulation and range history of a moving target are obviously different from a stationary target, a signal geometry model for moving targets in multichannel GEO ‐ SA ‐ BiSAR is established in this paper. According to simulation results, the along track velocity introduces target defocusing in azimuth, and the slant range velocity mainly causes multiple false targets. To resolve these problems, a modified multichannel reconstruction method in azimuth channel GEO ‐ SA ‐ BiSAR is proposed according to the azimuth multichannel impulse response of the imaged moving target. Before azimuth multichannel raw data combination, both spatial ‐ variant range cell migration correction (RCMC) and azimuth nonlinear chirp scaling (ANLCS) should be performed to reduce the influence of the range offset and lower the Doppler bandwidth of the whole raw data, respectively. Afterward, a novel azimuth multichannel reconstruction algorithm is carried out via the modified reconstruction matrix based on the estimated target velocity. The target slant range velocity estimation is implemented by introducing the signal intensity ratio (SIR). Compared with the conventional method for the stationary target to handle the raw data of the moving target, the false targets could be obviously suppressed by using the proposed approach. Imaging results on both simulated point and distributed scene targets validate the proposed multichannel reconstruction approach.


Introduction
Synthetic aperture radar (SAR) can get microwave images of the desired scene under all weather conditions, day and night [1]. For spaceborne SAR, most SAR systems are in a low earth orbit (LEO) at heights ranging from 300 km to 1200 km. With the development of the spaceborne SAR system, the disadvantages of LEO SAR, such as small instantaneous beam coverage and long revisit time, are becoming more and more obvious. To resolve these problems, the concept of geostationary SAR (GEO SAR) was proposed [2][3][4][5], which has an orbit height of nearly 36,500 km. Consequently, compared with LEO SAR, a larger wide-swath could be illuminated with a shorter revisit time by GEO SAR [6,7]. Therefore, GEO SAR has many advantages in military reconnaissance, ocean applications, and earth dynamics [8].
Although GEO SAR could overcome major shortcomings of LEO SAR, it is very difficult to get microwave images with high azimuth resolution and signal-to-noise ratio (SNR) [9]. Fortunately, compared with monostatic GEO SAR with the same transmit peak power, the bistatic SAR imaging scheme consisting of a GEO transmitter and an airborne receiver, named as geosynchronous spaceborne-airborne bistatic synthetic aperture radar (GEO-SA-BiSAR), could obtain SAR images with finer azimuth resolution and higher SNR [10,11]. Furthermore, compared with traditional monostatic SAR systems, the GEO-SA-BiSAR system has enormous advantages, such as well concealment, strong anti-attack ability, real-time SAR imaging capacity, and well flexibility [12,13]. Therefore, the unique characteristics of GEO-SA-BiSAR would make it be widely adopted in future. Compared with conventional LEO SAR, azimuth precisely data focusing in GEO-SA-BiSAR is more difficult. This is because the slant range history is more complex than other SAR systems, and the low operating pulse repetition frequency (PRF) in GEO SAR results in Doppler spectrum aliasing [14]. In GEO-SA-BiSAR systems, phase and time synchronization are very important issues, especially for the imaged scene with the low SNR level. In recent years, multiple phase/time synchronization processing approaches have been proposed [15][16][17]. In the GEO-SA-BiSAR system, the direct-path signal is received as the reference signal by an additional receive channel for phase and time synchronization. For the ocean scene with the low SNR, phase and time synchronization are usually implemented according to strong scatters such as small islands and stationary ships.
The operating PRF of the GEO-SA-BiSAR system is determined by the GEO transmitter, and a low PRF is chosen to obtain wide swath coverage and to lessen the burden of raw data storage onboard and down link transmission. The wide Doppler bandwidth for a fine azimuth resolution depends on the airborne receiver, because the Doppler modulation rate corresponding to the airborne receiver is obviously larger than the GEO transmitter. Therefore, the GEO-SA-BiSAR system brings a contradiction between the low PRF of the transmitter and the wide Doppler bandwidth determined by the receiver. Similar to the high resolution wide swath (HRWS) SAR system [18], the azimuth multichannel imaging scheme could resolve this contradiction by introducing additional azimuth spatial samples. The azimuth multichannel receiving technology is composed of multiple receiving sub-apertures on the receiving airborne platform, which can improve the target Doppler bandwidth for a better azimuth resolution but with a low PRF sampling rate. To combine azimuth multichannel SAR together before the single channel GEO-SA-BiSAR processor, azimuth multichannel reconstruction usually should be taken because of non-uniform sampling. In recent years, multiple reconstruction and imaging approaches for the stationary target have been proposed [19][20][21], but few studies have focused on the moving target in GEO-SA-BiSAR.
Compared with the stationary target in multichannel GEO-SA-BiSAR [22,23], the moving target has a different range history and Doppler modulation [24,25]. If azimuth multichannel reconstruction and imaging approaches for stationary scenes were adopted to handle the echoes of moving targets, imaging quality would be seriously reduced. To resolve this issue, the geometry model of the moving target in GEO-SA-BiSAR is analyzed in this paper, and the effects of the target velocity including the along track and slant range components on the focusing quality are analyzed. According to the azimuth multichannel impulse function, a novel azimuth multichannel reconstruction method for moving targets in GEO-SA-BiSAR is proposed. Before combining azimuth multichannel raw data, spatial-variant range cell migration correction (RCMC) is performed to reduce the influence of the range offset, and the Doppler bandwidth of the whole raw data is decreased via azimuth nonlinear chirp scaling (ANLCS). The key fact of the novel reconstruction approach is the modified reconstruction matrix based on the estimated moving target velocity. The along track velocity component is easily obtained by azimuth data autofocusing similarly to the single-channel SAR system. For the estimation of the slant range velocity component, a signal intensity ratio (SIR) is introduced to describe the ratio between the in-band and out-of-band signal after azimuth multichannel reconstruction with the estimated target velocity [26,27]. After azimuth multichannel reconstruction with the accurate target velocity, the moving target is well focused via the traditional image processing technique for moving targets, and its corresponding false targets due to the target velocity are obviously suppressed. This paper is organized as follows. The signal model of the moving target in azimuth multichannel GEO-SA-BiSAR is established in Section 2. The effects of the moving target velocity on the target focusing quality are analyzed in Section 3. The proposed reconstruction approach is presented in detail in Section 4. Simulation results for both point and distributed targets are analyzed in Section 5 to validate the proposed reconstruction approach. Conclusions are reported in Section 6.

Geometry and Azimuth Signal Model for Moving Targets in GEO-SA-BiSAR
Because of the low operating PRF in GEO SAR, the azimuth raw data would be under-sampled and its spectrum would be aliased, if an azimuth signal with a large azimuth bandwidth were received by the airborne platform. There are two ways to resolve this contradiction. The first way is to improve the operating PRF value of the GEO SAR system, although this raises system complexity. The other way is to use azimuth multichannel reception technology on the airborne receiving platform to increase azimuth spatial samplings. The azimuth multichannel reception technique can improve the effective azimuth sampling rate without increasing PRF, and it has been implemented by multiple airborne/spaceborne SAR systems.

Geometry and Slant Range
The moving target imaging geometry in GEO-SA-BiSAR is illustrated in Figure 1. GEO SAR transmits radar signals, and the airborne receiver receives echoes reflected by the moving target. The GEO SAR moves with the velocity T v , and the velocity of the airborne is R v in the y-direction. The direct-path signal from the GEO transmitter is received by the airborne platform for time/phase synchronization. The M receiving sub-apertures are aligned in azimuth, and the physical interval between the m-th receiving sub-aperture and the reference is Rxm Figure 1. The observation geometry in azimuth multichannel geosynchronous spaceborne-airborne bistatic synthetic aperture radar (GEO-SA-BiSAR).
where 0 R denotes the shortest slant range from the receiver to the imaged moving target, and t is the azimuth time. The slant range history of (1) can be expanded into a Taylor series as follows: The slant range history of (4) expanded into a Taylor series can be expressed as Combining (3) with (5), the following expression can be obtained: For a GEO SAR with a high orbit, a fourth-order approximation of   T R t is usually modeled and established to accurately estimate the actual slant range between the GEO SAR system and the target with a long synthetic aperture time [28], and it is written as where 1T R denotes the shortest distance between the moving target and the GEO transmitter, and   1, 2, 3, 4 nT k n  represents the polynomial coefficients of the n-th order. However, the synthetic aperture time of the GEO-SA-BiSAR system is much smaller than monostatic GEO-SAR, and the phase errors caused by the last two terms are less than 4 1 10    [24], which could be neglected. Consequently, where t A , t v , and t r indicate the acceleration, velocity, and range vectors of the illuminated moving target in the earth-centered inertial coordinate system, respectively, while s A , s v , and s r represent the acceleration, velocity, and range vectors of the GEO SAR platform, respectively. The distance obtained from (6) can be substituted into (9) to obtain: With the parameters in Table 1, the maximum value of 2T k is about 0.03, while the maximum slant range errors of the second and third terms in the whole synthetic aperture time are 2.1 10 m  , respectively. These range errors are much less than a quarter of the wavelength and can be ignored. Consequently, combining (6) and (12) together, the total slant range can be rewritten as follows:

Azimuth Multichannel Response Model
For an illuminated moving target, the baseband raw data of the m-th receiving sub-aperture can be written as where A is a complex constant,   r w  and   a w  indicate the transmit and receive pulse envelopes, respectively,  is the fast time, c is the light speed,  represents the carrier wavelength, r K denotes the chirp rate, and p  indicates the duration of the transmitted pulse.
Since this paper is mainly focused on multichannel reconstruction in azimuth, the azimuth signal component of (15) is given as Substituting (14) into (16), the following expression is obtained: For the stationary target, the slant range history   fixed,m R t for the m-th receiving channel of the stationary target can be written as stationary target can be written as Consequently, the azimuth impulse response of the m-th receiving sub-aperture for the stationary target can be rewritten as [24]   fixed, fixed, fixed, According to (17) and (22), the multichannel impulse response for an imaged moving target can also be obtained from the single channel GEO-BiSAR azimuth response as follows:

Doppler Bandwidth and PRF Analysis
The azimuth bandwidth of GEO-SA-BiSAR can be written as where T is the illumination time for a point target in GEO-SA-BiSAR. With the parameters listed in Table 1, the azimuth bandwidth corresponding to the airborne multichannel receiver is 337.5 Hz, the Doppler bandwidth determined by the GEO transmitter is only 1.1 Hz (inclination is 60°), and the adopted synthetic aperture time is assumed to be 4.5 s.
The timing diagram selection result of the GEO SAR system for the swath coverage of 500 km is shown in Figure 2. The blue and red areas indicate the interference areas caused by nadir echoes and transmitted pulses, respectively, while these areas must be avoided during the timing diagram design. Generally speaking, because of the narrow instantaneous processed azimuth beam bandwidth of GEO monostatic SAR, the operating PRF in GEO monostatic SAR systems is usually less than 200 Hz in order to obtain the wide swath GEO SAR raw data with azimuth unambiguous sampling [29], [30]. To obtain the desired area with 500 km ground range swath width, selected PRFs are indicated by the regions with black lines as shown in Figure 2, and the values of the selected PRF are in the range of 160-200 Hz. Therefore, the selected PRF value for the following SAR data simulation is set to 180 Hz. According to the above analysis, the exploited azimuth bandwidth of the GEO-SA-BiSAR system is 338.6 Hz and larger than the adopted system PRF. Consequently, Doppler spectrum aliasing would occur as a result of the contradiction between the wider Doppler bandwidth because of the receiver velocity and the lower operating PRF in GEO SAR. In order to overcome the contradiction and resolve the Doppler spectrum aliasing, the azimuth multichannel receiving technology can be applied on the receiver platform. According to the principle of azimuth multichannel SAR [31], [32], the total effective azimuth sampling rate with M azimuth channels can be increased to . Although just two azimuth receiving channels are needed in the airborne platform, three azimuth receiving channels are designed for a higher azimuth over-sampling rate in order to obtain a better Doppler spectrum ambiguity suppression effect.

Effects of the Target Velocity on Imaging Results
In the azimuth multichannel GEO-SA-BiSAR system, the different range history of the moving target causes the different multichannel impulse response in (23) with azimuth time delay and phase shift. To analyze the effects of the target velocity on these parameters in detail, designed simulation experiments are implemented, and Table 1 lists the corresponding simulation parameters.  Figure 3 shows the influences of the along track velocity component on the three parameters aforementioned. The relative sensor velocity and azimuth time delay are sensitive to the along track velocity, but the phase shift is independent of the along track velocity. In Figure 3(a), the relative sensor velocity varies from approximately 300 m/s to 280 m/s as the along track velocity changes. The influence of the along track velocity on the azimuth time delay is shown in Figure 3(b), and the variation of the azimuth time delay with the change of the along track velocity is responsible for the changed phase function for azimuth multichannel reconstruction. In Figure 3(c), the phase shift is independent of the along track velocity. In the previous simulation in Figure 3 to present the azimuth time delay and the phase shift varying with the along track velocity, parameters of the third channel are calculated, while the middle azimuth channel is considered as the reference channel.

Effects on Imaging Results
To validate the above analysis, imaging simulation experiments on point targets are designed according to the parameters in Table 1. As shown in Figures 5-7, both slant range and along track velocities lead to azimuth defocusing and false targets. However, the azimuth defocusing phenomenon due to the azimuth velocity component is much more serious, as shown in Figure 6(c), while the level of false targets introduced by the slant range velocity is obviously higher, as shown in Figure 7(d). With the same value, the Azimuth time delay(ms) along track velocity introduces a larger relative sensor velocity change, which results in the larger additional phase error responsible for azimuth defocusing as shown in Figure 8(a). Meanwhile, Figure 8(b) shows the maximum azimuth mismatch phase error in the synthetic aperture time caused by the slant range velocity, which leads to azimuth defocusing. However, the formation mechanisms of false targets due to the azimuth and slant range velocity component are quite different. The reason for false targets caused by the azimuth velocity component is azimuth non-uniform sampling because of the changed azimuth relative velocity, whereas the reason for false targets introduced by the slant range velocity is azimuth multichannel phase imbalance as shown in Figure 8(c), since the phase shift is sensitive to the along track velocity. Taking the azimuth middle channel as the reference, Figure  8 Besides azimuth defocusing and false targets, the slant range also leads to the moving target location shift as shown in Figure 7(b) and (d), and the azimuth location shift is caused by the additional Doppler shift due to the slant range and the changed azimuth modulation rate. Similar to the conventional SAR, the azimuth location shift is computed as As simulation results shown in Figure 5-8, the azimuth velocity component mainly causes target defocusing in azimuth because of the azimuth Doppler modulation rate mismatch, whereas the slant range velocity component mainly causes azimuth false targets because of azimuth multichannel phase imbalance.

Azimuth Doppler Spectrum Reconstruction
As the above analysis indicates, the total effective azimuth Doppler sampling rate could be multiplied by arranging multiple receiving sub-apertures along the azimuth direction in a GEO-SA-BiSAR system. However, if the optimum PRF is not selected, additional azimuth multichannel reconstruction would be required to resolve the non-uniform sampling problem before the imaging processor. If the traditional azimuth multichannel reconstruction and imaging approaches were still adopted to process the echoes of moving targets, as shown in Figures 6 and 7, imaging quality would be seriously reduced. Therefore, an innovative reconstruction approach based on the moving target azimuth multichannel response in (23) is proposed. Before combining azimuth multichannel raw data together, spatial-variant RCMC is performed to reduce the effects of the range cell migration on multichannel reconstruction and the Doppler bandwidth of the whole raw data is reduced by ANLCS. The key point of this approach is the azimuth multichannel reconstruction filters with the estimated velocity of the moving target. After the modified azimuth multichannel reconstruction, the moving target could be well focused, and its corresponding false targets could be obviously suppressed. The flowchart of the presented multichannel reconstruction approach is shown in Figure  9, and it includes three major processing steps: RCMC, azimuth de-ramping, and azimuth multichannel raw data combination. Compared with the large stationary scene, the moving target is always of a small size, and its corresponding echoes are submerged in reflected echoes of the large stationary scene. After range compression, spatial-variant RCMC should be applied to reduce the effects of range cell migration (RCM) before azimuth multichannel reconstruction. The spatial-variant RCMC in GEO-SA-BiSAR could be achieved via keystone transform, and its detailed processing steps can be found in [25,33].
Furthermore, the instantaneous Doppler centroid varies with both the slant range and the azimuth location of the GEO transmitting platform, which extends the Doppler bandwidth of the raw data. Therefore, the varied instantaneous Doppler centroid should be removed before azimuth multichannel reconstruction by ANLCS, and ANLCS is implemented by multiplying the raw data by where   2,3, 4 i q i  are the coefficients of the ANLCS operation [24].
Transforming the multichannel impulse response of the moving target in (23) where f represents the Doppler frequency. Consequently, the azimuth multichannel system matrix for the moving target in GEO-SA-BiSAR is given as follows: Afterward, the azimuth multichannel reconstruction function ( ) f P is expressed as The reconstruction function   Afterward, the unambiguous spectrum of the imaged moving targets could be obtained by combining the echoes from all azimuth receiving channels. Finally, the equivalent raw data received by a single channel are obtained after azimuth inverse fast Fourier transform (IFFT).

Moving Target Velocity Estimation
According to (34), the azimuth multichannel reconstruction matrix in (36) for the imaged moving target is obtained from the given target velocity components in both directions. Therefore, the moving target velocity estimation is very important for azimuth multichannel reconstruction. As the relative sensor velocity change is mainly related to the along track velocity component, this component can be predicted by azimuth auto-focusing just like in the traditional case. However, for each azimuth channel, the azimuth signal is under-sampled and its corresponding spectrum is aliased, and the slant velocity cannot be obtained via the conventional Doppler centroid estimation approach.
With the system parameters listed in Table 1, both Figures 10 and 11 show the azimuth multichannel reconstructed simulation results with different estimated slant range velocities when the designed slant range velocity is 10 m/s. Figure 10 shows the azimuth band-limited case, in which only the exploited azimuth beam width is considered, while Figure 11 shows the non-band-limited case, in which both in-band and out-of-band azimuth signals are considered during echo simulation of moving targets. In both cases, for the well reconstructed Doppler spectra with the accurate slant velocity, most signal power would gather around the estimated Doppler centroid as shown in Figures  10 and 11. (a) (b) (c) According to the signal power distribution phenomenon shown in Figures 10 and 11, the SIR could be introduced to judge whether the multichannel spectrum is well reconstructed. Furthermore, the SIR could be expressed as two types as follows: with where a B is the processed Doppler bandwidth,

 
; r W f u indicates the reconstructed Doppler spectrum, and ˆr u is the estimated slant range velocity. The SIR 1  represents the power ratio between in-band and out-of-band data, while the SIR 2  indicates the average amplitude ratio between in-band and out-of-band data. To improve the slant range velocity estimation accuracy, the product 1 2     of these two ratio parameters is used to estimate.
The introduced SIR  varying with the estimated slant range velocity is shown in Figure 12, and all SIRs reach the maximum value only when the accurate slant range velocity is adopted for azimuth multichannel reconstruction. Furthermore, the curve of  is the steepest among the three curves, which would result in the most accurate velocity estimation result. From Figure 12, it can be found that SIR  can well estimate the velocity of the moving target in both cases of band-limited and non-band-limited. As we know, the moving target velocity estimation accuracy is mainly affected by the signal-toclutter-and-noise ratio (SCNR) and ambiguity-to-signal ratio (ASR) levels, since echoes of the moving target are mixed with its corresponding ambiguity signal and echoes of the stationary scene in the receiving window. For SAR systems with airborne receivers, only azimuth ambiguity energy should be considered, and a higher operating PRF will obtain a lower ASR level. Consequently, to evaluate the accuracy of the presented velocity estimation approach based on the introduced SIR  , two sets of experiments on point targets are carried out. One is carried out to estimate the moving target slant range velocity under different SCNRs and selected PRF values, and the other is designed for different velocity estimations with the same SCNR and PRF. An L-band SAR image of the ocean is used to simulate the clutter signal, which is added to the simulated raw data of the moving target to obtain different SCNRs. In the first group of simulation experiments, the designed slant range velocity is set to 10 m/s. The velocity estimation results under different SCNRs and selected PRF values are obtained and summarized in Table 2. As shown in Table 2, the higher SCNR and PRF are with, the more accurate velocity estimation result would be achieved. Compared with the cases with different SCNRs, the effects of different PRF values are very small and can even be neglected, since all selected PRF values in Table 2 have a high azimuth over-sampling rate, which results in a low ASR level. To obtain the relative estimation error below 5%, the SCNR level should be greater than -10 dB. Therefore, clutter suppression should be carried out before the slant range velocity estimation, and the SCNR usually would be more than 0 dB after clutter suppression. As several clutter suppression methods in HRWS SAR have been proposed in recent years [34][35][36], clutter suppression will not be discussed in this paper. For simulation experiments in Table 3, the moving target is designed with the same SCNR but different slant range velocities. The simulated SCNR is set to -10 dB, and the estimated velocities are summerized in Table 3. As shown in Table 3, all deviations between the designed and estimated slant range velocities are below 0.15 m/s, and all relative errors are below 3%. The largest relative error of the estimation result is the sixth target, because the designed slant range velocity of the sixth moving target is too small at 1 m/s.

Moving Target Imaging in GEO-SA-BiSAR
In the real SAR raw data, echoes of the moving target are mixed with echoes of the stationary scene. Since the stationary target and the moving target have different Doppler characteristics and range cell migrations (RCMs) in azimuth multichannel GEO-SA-BiSAR, azimuth multichannel reconstruction and two-dimension (2D) focusing for the stationary scene and the moving target should be carried out separately. Therefore, the proposed imaging approach for moving targets in azimuth multichannel GEO-SA-BiSAR is shown in Figure 13, and the proposed imaging approach consists of three major parts: stationary scene processing, moving target processing, and image data fusion. Azimuth multichannel imbalance, including amplitude, phase, and time delay, would seriously reduce the imaging quality, so multichannel imbalance correction should be carried out before azimuth multichannel reconstruction. Afterward, azimuth multichannel reconstruction and 2D focusing approaches suitable for azimuth multichannel GEO-SA-BiSAR in [24] should be performed. According to the aforementioned analysis and simulation results, the along track velocity component mainly leads to defocusing in azimuth, while the slant range velocity mainly results in false targets. Besides recently proposed moving target detection approaches [37], we could use this phenomenon to detect the moving target, especially for ocean scenes. After azimuth multichannel reconstruction and clutter suppression via Doppler filtering, moving targets can be detected in the obtain SAR image from the conventional imaging processor, and their corresponding positions can be approximately obtained. Since echoes of these moving targets are compressed in very limited range samples after range compression and RCMC, echoes of moving targets with very limited samples in both azimuth and range directions around their approximate position are extracted in the range azimuth-time domain before azimuth focusing. After extracting the echoes of these detected moving targets, the data of the corresponding areas is set to zero, and then echoes of the stationary scene is obtained after inverse focusing operation. Consequently, echoes of the stationary scene and moving targets could be processed separately. If different moving targets are far enough to each other, the echoes of these moving targets can be extracted and processed separately. However, the echoes of moving targets cannot be extracted individually, if they are too close to each other. In this case, the speed of one moving target or the average speed of multiple closely distributed moving targets can be used for imaging, at the cost of the imaging quality. But in practice, moving targets in the ocean scene usually will not be too close, unless they move together with almost the same speed.
In terms of processing extracted echoes of moving targets, both velocity components of moving targets should be alternately estimated at first. The along track velocity is estimated by azimuth autofocusing, while the SIR  is introduced for the slant range velocity estimation. With the accurately estimated target velocity, azimuth multichannel echoes of the moving target could be well reconstruction. Finally, the reconstructed azimuth multichannel raw data of the moving target could be easily handled by the single channel imaging processor of GEO-SA-BiSAR [38].
After removing echoes of moving targets from the whole raw data of the imaged scene, the stationary imaged scene should be focused again. With the estimated target velocity including both velocity components in the azimuth and range directions, the moving target location correction can be easily implemented according to (30). Finally, the superposition result of the focused stationary scene and moving target focusing result after location correction can be obtained. It is noticed that the echo location of the moving target in the whole SAR raw data is known during extracting the echo signal of the moving target. Therefore, it is only need to superimpose the focused moving target on the stationary scene according to its corresponding location.

Simulation Experiment
In this section, simulation experiments on both point and distributed targets are designed to evaluate the proposed processing approach, and simulation parameters are summarized in Table 1. The height of the GEO transmit platform is 36,500 km, and its corresponding velocity along the azimuth direction is 2200 m/s. The height of the airborne receiver platform with three azimuth receiving channels is 3 km, and its velocity is 300 m/s. The synthetic aperture time for this GEO-SA-BiSAR is about 4.5 s, and its corresponding Doppler bandwidth is about 337.5 Hz. Because of the transmitted PRF of 180 Hz, the Doppler spectrum in each azimuth channel is aliased.

Simulation on Point Targets
First, two simulation experiments on a single moving target are carried out to demonstrate the proposed reconstruction approach. One single-point target is designed with only an along track velocity component of 10 m/s, and the other single-point target is designed with only a slant range velocity of 10 m/s. Their corresponding multichannel reconstruction and imaging results are shown in Figures 14 and 15, respectively. Compared with the imaging results in Figures 6 and 7, the imaging qualities are obviously improved as shown in Figures 14 and 15   Furthermore, a simulation experiment on a designed scene with four point targets is carried out, and the relative geometrical relationship of the four point targets is shown in Figure 16(a). All point targets have different moving velocities: target P1 has an along track velocity of 10 m/s, target P2 has a slant range velocity of 10 m/s, target P3 has an along track velocity of 5 m/s and a slant range velocity of 10 m/s, and target P4 has an along track velocity of 10 m/s and a slant range velocity of 5 m/s. With the conventional azimuth multichannel reconstruction approach, all point targets are defocused with obviously false targets as shown in Figure 16(b). After adopting the proposed azimuth multichannel reconstruction and imaging approaches, the four targets are well focused and their false targets are obviously suppressed as shown in Figure 16(c). Interpolated contour plots of the four points are shown in Figures 16(d)-(g), it can be seen that each target is well focused, and their corresponding imaging quality parameters including resolution (res), integrated-side-lobe ratio (ISLR), peak-sidelobe ratio (PSLR) and maximum false target amplitude (MFTA) are measured and summarized in Table 4. Compared with the results obtained from the conventional imaging processor for the stationary scene and the theoretical values, the proposed approach shows better imaging capacity and false target suppression ability for moving targets.   Figure 17 shows the maximum false target amplitude by using the conventional and proposed approaches in both band-limited and non-band-limited cases. In the band-limited case, false targets are suppressed to below -60 dB and can be ignored. In the non-band-limited condition, false targets are suppressed to below -25 dB, the reduced amplitude level of false targets is less than the maximum azimuth ambiguity power in conventional SAR, and they almost would not affect following target detection and recognition.

Simulation of Distributed Targets
To further evaluate the proposed multichannel reconstruction approach, a distributed ocean imaging scene simulation including three ships (denoted as ship 1, ship 2, and ship 3) is designed. The gray amplitude of a spaceborne ocean SAR image is set to be the input radar cross section (RCS) of the designed ocean, while the input RCSs of the three ships are also obtained from the real focused SAR image. The first ship is stationary and designed for comparison, the second ship moves along the slant range direction, and the third ship has both along track and slant range velocity components. The slant range velocity of ship 2 is 10 m/s, while both the along track and slant range velocity components of ship 3 are set to 10 m/s. The parameters of the multichannel GEO-SA-BiSAR system for this simulation experiment are listed in Table 1. Figure 18(a) shows the imaging result after conventional azimuth multichannel reconstruction. Since ship 2 and ship 3 have slant range velocity components, their imaging result are with a pair of serious false targets and azimuth location shift. Furthermore, ship 3 is defocused in azimuth because of the along track velocity component. In Figure 18(a), ship 2 and ship 3 have the same slant range velocity, but their corresponding distances between the real target and the adjacent false target are different. The reason for this phenomenon is that the distance is determined by both the slant range velocity and the slant range from the moving target to the airborne receiver. In Figure 18, the slant range of ship 2 and ship 3 are 6490 m and 4860 m, respectively, and their corresponding distances between the real target and the false target in theory are 934.58 m and 723.99 m, respectively. After measurement in Figure 18(a), measured distances agree with the theoretical values. At the same time, according to (30), the resulting azimuth location shifts of ship 2 and ship 3 are also different, and their azimuth location shifts are 216.3 m and 167.6 m, respectively. Figure 18(b) shows the imaging result of the proposed azimuth multichannel reconstruction approach, the false targets of moving targets are well suppressed, while the defocused ship 3 in Figure 18(a) is well focused in Figure 18(b). According to the location shift in (30), the moving target location can be corrected as shown in Figure  18(c).

Conclusions
GEO-SA-BiSAR systems can obtain SAR images with finer azimuth resolution and higher SNR than monostatic GEO SAR systems with the same power consumption. However, the operating PRF of a GEO-SA-BiSAR system is less than the wide Doppler bandwidth determined by the airborne receiver, which leads to azimuth spectrum aliasing. In order to resolve this problem, azimuth multichannel technology is adopted in the airborne receiver platform in this study. Moreover, azimuth multichannel reconstruction of the echo signal is usually required, when the optimum PRF is not selected. In recent years, many researchers have studied conventional azimuth multichannel reconstruction approaches to address this problem.
However, because of the different Doppler modulation and range history between the stationary target and the moving target, the imaging geometry model of the moving target in GEO-SA-BiSAR is established in this paper. According to signal analysis and simulation results, the along track velocity mainly leads to target defocusing in azimuth, and the slant range velocity mainly results in false targets. To solve these problems, the novel multichannel reconstruction approach based on the azimuth multichannel response of the moving target is proposed, and the key point of the proposed approach is the modified multichannel reconstruction matrix corresponding to the estimated target velocity.
As the slant range velocity component of the imaged moving target should be known before azimuth multichannel reconstruction, a signal identity ratio between the in-band and out-of-band spectra is introduced for the slant range velocity estimation, and the velocity estimation accuracy is mainly determined by the SCNR. For ocean scenes with a high SCNR level, clutter suppression may not be required, but for most land scenes, clutter suppression must be performed before the slant range velocity estimation.