Open AccessArticle

Research on Configuration Constraints of Airborne Bistatic SARs

Yidi Chen

Renwen Chen

Hao Liu

Jiapeng Guo

Yujie Wang

and

Junyi Zhang

College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Author to whom correspondence should be addressed.

Sensors 2022, 22(17), 6534; https://doi.org/10.3390/s22176534

Submission received: 17 July 2022 / Revised: 18 August 2022 / Accepted: 24 August 2022 / Published: 30 August 2022

(This article belongs to the Section Radar Sensors)

Download

Browse Figures

Figure 1
Imaging geometry model of bistatic forward-looking SAR. "> Figure 2
Description of RCMC. (a) Bistatic imaging geometry model, (b) before LRCMC, (c) after LRCMC. "> Figure 3
Flowchart of bistatic SAR imaging algorithm. "> Figure 4
Layout of scenic spots. "> Figure 5
Simulation of scenic spots. (a)Traditional azimuth pulse compression, (b) pulse compression after space variation correction. "> Figure 6
Schematic diagram of azimuth and range pulse at point B. (a) Schematic diagram of azimuth and range pulse at point C. (b) Schematic diagram of azimuth and range pulse at point B. "> Figure 7
Geometric model of bistatic forward-looking SAR imaging. "> Figure 8
Schematic diagram of area of resolution element. "> Figure 9
Range resolution and azimuth resolution. (a) Range resolution, (b) azimuth resolution. "> Figure 10
Angle of resolution. "> Figure 11
Range resolution and Azimuth resolution. (a) Range resolution, (b) azimuth resolution. "> Figure 12
Angle of resolution. ">

Versions Notes

Abstract

Based on the analysis of the airborne bistatic synthetic aperture radar (SAR) imaging geometric mode, an extended nonlinear chirp scaling algorithm is employed to simulate and verify the imaging effect of the bistatic SARs. A gradient theory-based two-dimensional resolution bistatic SAR model is proposed, and the constraints of the multi-platform flight trajectory parameters meeting the imaging accuracy of the bistatic SAR are analyzed. Finally, through the bistatic SAR imaging simulation of cooperative flight trajectories under various situations, the spatial configuration constraint envelope between the flight vehicles to achieve the optimal resolution is revealed. The results of this paper will provide a theoretical reference for the SAR application in formation flight control.

Keywords:

bistatic SAR; airborne; spatial configuration; radar imaging

1. Introduction

Bistatic synthetic aperture radar (SAR) puts the receiver and transmitter on different platforms. The transmitter actively transmits radar signals to the target area and the receiver passively receives the echo signals to achieve high-resolution two-dimensional imaging of the target area [1,2]. Bistatic SAR has outstanding advantages such as good concealment, strong target detection ability, strong anti-interference ability, forward-looking imaging [3], etc., and thus has a good application prospect.

Many researchers have made great efforts to promote bistatic SAR imaging accuracy and range. Martino studied an imaging geometry model and the processing algorithm of a spaceborne bistatic SAR system for formation flight in large-baseline mode [4], and Servidia proposed a control rate method of configuration reconstruction of a spaceborne bistatic SAR system for formation flight in different Earth observation modes [5]. Lu established an imaging model of the GEO-LEO bistatic interferometric SAR system, designed its formation configuration, and analyzed its interferometric measurement accuracy [6]. The above methods are all based on the satellite platform with an accurate and controllable trajectory. However, it will not be suitable in some cases when the flight trajectory is unknown in advance. To overcome this problem, Li established an imaging geometric model of a missile-borne forward bistatic SAR system and analyzed the influence relationship between the two-dimensional resolution and system parameters [7]. Meng established a general bistatic forward SAR signal model and analyzed the azimuth space variation effect and its spectral characteristics [8]. Ding analyzed a bistatic SAR signal model in the special mode, when the transmitting platform horizontally flies while the receiving platform pitches down, and proposed an extended chirp scaling (CS) imaging algorithm [9]. S. Li analyzed a signal model under the non-parallel flight trajectory of the receiving platform, and a nonlinear CS (NLCS) imaging algorithm based on a frequency extension was proposed [10].

The high mobility of the moving platform makes the bistatic SAR imaging resolution change significantly. Therefore, it is more important than ever to obtain the balance between the high mobility of the moving platform and the optimal resolution of bistatic SAR imaging by reasonably planning the motion trajectory of the transceiver platform. However, because the resolution of the bistatic SAR forward-looking imaging is significantly influenced by the geometric model parameters of the platform and the target, for practical applications, it is necessary to analyze the influence of the parameters of the receiving platform and the target geometric model and find the configuration boundary of the moving platform under the constraints of the optimal or necessary resolution.

This paper investigates the airborne bistatic SAR imaging geometry model and the extended NLCS algorithm and verifies the imaging effects by simulation. Under the premise of the imaging effect, the constraint relationship between the aircraft spatial configuration and the imaging resolution is analyzed. Finally, the configuration parameter envelope that meets with the predefined resolution is obtained by a bistatic SAR imaging simulation through all the cooperative flight routes.

2. Imaging Geometry Model of Bistatic SAR

A general bistatic geometric model [11] is shown in Figure 1. It is assumed that the transmitter and the receiver fly along an unparallel straight line at a constant speed. The transmitter’s flight speed is V_T and the receiver’s is V_R, and the two speeds are different. A right-hand coordinate system is adopted, and the x-y plane is the ground plane, at which the target point is located. The flight direction of the transmitter is parallel to the y axis.

In the bistatic model, both the radar wave beams of the transmitter and the receiver are at the same direction and cover the reference targets. The instantaneous slope range between the transmitter and the target is

R_{T} (η)

, the corresponding receiving instantaneous slope range is

R_{R} (η)

, the time at the range direction is

τ

, and the time at azimuth direction is

η

. When the radar beam center irradiates on the reference target point, the azimuth time is

η = 0

, at which the transmitting angle and receiving angle are

θ_{s q T}

and

θ_{s q R}

, respectively.

Assuming that the transmitted broadband signal is

s (τ)

, up-conversion of the signal at frequency

f_{o}

will result in

s_{t} (τ) = Re {s (τ) \exp (j 2 π f_{0} τ)}

(1)

and the two-way range

R (η)

of the echo signal from the target will be

s_{r} (τ, η) = w_{a z} s_{t} (τ - \frac{R (η)}{c})

(2)

where

w_{a z}

is the azimuth envelope, and

c

is the speed of light. Here, the amplitude of the signal is ignored.

R (η)

[12] may be described as

\begin{array}{l} R (η) & = R_{T} (η) + R_{R} (η) \\ = \sqrt{V_{T}^{2} η^{2} + R_{T c e n}^{2} - 2 V_{T} η R_{T c e n} \sin (θ_{s q T})} \\ + \sqrt{V_{R}^{2} η^{2} + R_{R c e n}^{2} - 2 V_{R} η R_{R c e n} \sin (θ_{s q R})} \end{array}

(3)

In Equation (3),

R_{T c e n}

is the range between the transmitting platform and the target at

η = 0

moment, and

R_{R c e n}

is the range between the receiving platform and the target at

η = 0

moment.

Through down conversion, one can obtain the equation below

s (τ, η) = w_{a z} (η) w_{r} (τ - \frac{R (η)}{c}) \exp {- j 2 π \frac{R (η)}{λ}}

(4)

where

w_{r} (*)

is the range envelope and

w_{a z} (*)

is the azimuth envelope, which are determined by the composed antenna azimuth diagram.

3. Imaging Principle Model of Bistatic SAR

3.1. Analysis of Bistatic SAR Echo Signal

When the bistatic SAR echo signal model is analyzed, the time domain expression of target signal may be written as [11]

s_{A} (τ, η) = w_{a z} (η) s_{r} (τ - \frac{R (η)}{c}) \exp {- j 2 π \frac{R_{1} (η)}{λ}}

(5)

where

R_{1} (η) = R_{c e n} + k_{2} η^{2} + k_{3} η^{3} + k_{4} η^{4} + \dots

(6)

R_{1} (η)

is the range after the linear components’ compensation, and Fourier series expansion is performed at azimuth zero,

R_{c e n}

is the sum of

R_{T c e n}

and

R_{R c e n}

Using series inversion method [11], the two-dimensional spectrum expression of echo signal can be derived as [11]

S (f_{τ}, f_{η}) = W_{r} (f_{τ}) W_{a z} (f_{η} + (f_{o} + f_{τ}) \frac{k_{1}}{c}) \exp {j ϕ (f_{τ}, f_{η})}

(7)

where

ϕ (f_{τ}, f_{η})

is the phase [11]

ϕ (f_{τ}, f_{η}) = - 2 π (\frac{f_{o} + f_{τ}}{c}) R_{c e n} - \frac{π f_{τ}^{2}}{K_{r}} + 2 π \frac{c}{4 k_{2} (f_{o} + f_{τ})} {(f_{η} + (f_{o} + f_{τ}) \frac{k_{1}}{c})}^{2} + 2 π \frac{c^{2} k_{3}}{8 k_{2}^{3} {(f_{o} + f_{τ})}^{2}} {(f_{η} + (f_{o} + f_{τ}) \frac{k_{1}}{c})}^{3} + 2 π \frac{c^{3} (9 k_{3}^{2} - 4 k_{2} k_{4})}{64 k_{2}^{5} {(f_{o} + f_{τ})}^{3}} {(f_{η} + (f_{o} + f_{τ}) \frac{k_{1}}{c})}^{4}

(8)

3.2. Bistatic SAR Imaging Algorithm

For the bistatic SAR imaging processing algorithm, Wang et al. derived an accurate bistatic point target reference spectrum and presented a frequency domain-based focusing algorithm [13]; some algorithms such as fast filtering back-projection and a nonlinear CS (NLCS) algorithm were proposed in papers [14,15,16,17,18,19,20,21] in which the bistatic SAR imaging of aircraft platform is processed based on the extended NLCS algorithm.

Suppose that four point targets are placed in the bistatic SAR imaging model. When the compound beam center passes through target B and target C, the minimum slope ranges of the two targets are the same, and the compound beam center passes through target B and target A at the same time. The slope of linear range cell migration (LRCM) of targets A, C, and D is consistent, as shown in Figure 2.

The signal expression of target point A after linear range cell migration correction (LRCMC) and linear phase removal is as follows [12]

s_{A} (τ, η) \approx ρ_{r} {τ - \frac{R_{l r c m A} (η)}{C}} w_{a z} (η) \exp {- j 2 π \frac{R_{l r c m A} (η)}{λ}}

(9)

where

ρ_{r} (*)

is the range pulse compression envelope.

Similarly, the signal expression of target point C after LRCMC may be deduced as follows

s_{c} (τ, η) \approx ρ_{r} {τ - \frac{R_{l r c m C} (η)}{c}} w_{a z} (η - η_{C}) \exp {- j 2 π \frac{R_{l r c m C} (η)}{λ}}

(10)

where

R_{l r c m C} (η) = R_{c e n A} + k_{C 2} {(η - η_{C})}^{2} + k_{C 3} {(η - η_{C})}^{3} + k_{C 4} {(η - η_{C})}^{4} + \dots

(11)

Because

k_{A 2} \neq k_{C 2}

, the Doppler frequency modulation slope of target point A and target point C is inconsistent, so a single azimuth matched filter cannot be applied for azimuth compression. Before azimuth matched filtering, the quadratic compression term of range and azimuth coupling is removed.

NLCS algorithm’s disturbance function is used to equalize the quadratic phase term of the target point, and the disturbance function of target point A is introduced as follows [12]

s_{A p e r t} (τ, η) = s_{A} (τ, η) \exp {j π α η^{3}}

(12)

where

α = \frac{1}{3} [\frac{V_{T}^{2} \cos^{2} θ_{s q T}}{λ R_{T c e n A} R_{c e n A}} + \frac{V_{R}^{2} \cos^{2} θ_{s q R}}{λ R_{R c e n A} R_{c e n A}}] \times (V_{T} \sin θ_{s q T} + V_{R} \sin θ_{s q R})

(13)

Then, after the phase term of target point C is multiplied by the disturbance function

s_{C p e r t}^{'} (τ, η_{1}) \approx ρ_{r} (τ - \frac{R_{c e n A}}{c}) w_{a z} (η_{1}) \times \exp {j π α η^{3}} \exp {j π α 3 η_{C}^{2} η_{1}} \exp {j π α η_{C}^{3}} \times \exp {j π (\frac{V_{T}^{2} \cos^{2} θ_{s q T}}{R_{T c e n A}} + \frac{V_{R}^{2} \cos^{2} θ_{s q R}}{R_{R c e n A}}) η_{1}^{2}}

(14)

Among them, the first term is the cubic phase modulation term, and all target points are consistent.

The second term is the Doppler offset, and all target points have the same offset. Usually, the azimuth sampling frequency is 20% greater than the Doppler bandwidth, which will not cause aliasing [22]. Therefore, it is necessary to control the offset within a certain range to avoid aliasing.

The third term is a constant, which depends on the position of the target point and has no effect on the focus of the other point.

The fourth term is azimuth linear frequency modulation, and all targets are tuned at the same frequency.

Because the disturbance function only considers the quadratic phase term and does not consider the higher-order term, the higher-order phase higher than the quadratic is not equalized. Therefore, the influence of the higher-order term should be considered in azimuth matched filtering. Meanwhile, in the shift-invariant region of imaging, the influence of higher-order terms can be ignored. Therefore, only the target point A can be used to deduce the azimuth matched filter function.

The azimuth Fourier transform of the target point A is carried out as [12]

S_{a z A} (τ, f_{η}) \approx \int w_{r} (τ - \frac{R_{c e n A}}{c}) w_{a z} (η) \exp {- j 2 π f_{η} η} \times \exp {- j \frac{2 π}{λ} (k_{A 2} η^{2} + k_{A 3} η^{3} + k_{A 4} η^{4}) + j π α η^{3}} d η

(15)

Using the stationary phase principle and series inversion, the azimuth matched filter is obtained as the conjugate term, as is shown below

h_{a m f} (f_{η}) = \exp {- j ϕ_{a m f} (η (f_{η}))}

(16)

where

ϕ_{a m f} (f_{η}) = - \frac{2 π}{λ} (k_{A 2} η^{2} + k_{A 3} η^{3} + k_{A 4} η^{4}) + π α η^{3} - 2 π f_{η} η

(17)

and the relationship with frequency is shown in the following formula

η = A_{1} f_{η} + A_{2} f_{η}^{2} + A_{3} f_{η}^{3} + \dots

(18)

a_{1} = - \frac{2}{λ} k_{A 2} a_{2} = (\frac{3}{2} α - \frac{3}{λ} k_{A 3}) a_{3} = - \frac{4}{λ} k_{A 4}

(19)

A_{1} = \frac{1}{a_{1}} A_{2} = - \frac{a_{2}}{a_{1}^{3}} A_{3} = \frac{2 a_{2}^{2} - a_{1} a_{3}}{a_{1}^{5}}

(20)

The azimuth matched filter under each range element can be calculated numerically (k value depends on the range element). Finally, the focused images of all range units of the whole shift-invariant scene are obtained.

The flowchart of the proposed bistatic SAR imaging algorithm is shown in Figure 3. In the distance direction, the fast Fourier transform (FFT) algorithm is firstly performed on the original radar echo signal, and the echo signal is transformed into the frequency domain for processing. The received pulse is then compressed in the distance direction, so that the main energy of the signal is concentrated into a narrower duration. In order to eliminate the coupling between range and azimuth direction, the trajectory of the range migration curve is corrected to a straight line parallel to the azimuth, and LRCMC is used for correction. Then, perform inverse fast Fourier transform (IFFT) transformation on it, remove the linear phase, then FFT transformation. In the azimuth direction, first perform FFT transformation on the signal of the above operation and transform it into the frequency domain for processing. In order to eliminate the phase coupling distortion under strabismus or large aperture, the Second Range Compression (SRC) algorithm is used. The signal is then transformed into the time domain for residual calibration. In order to solve the nonlinear spatial variation of range migration (RCM) with range and the linear spatial variation of SRC with range of motion-invariant bistatic forward-looking SAR, the NLCS imaging algorithm is adopted. Then, perform FFT, azimuth compression, and IFFT transformation on it, and finally obtain focused image data.

3.3. Simulation Verification of Bistatic SAR Imaging Algorithm

The transmitter position error, receiver position error, and transmission line length measurement error affected the quality of imaging only slightly. All of these measurement errors caused an image shift along the range direction [23]. Therefore, the error analysis of bistatic SAR imaging is not discussed in the simulation. The configuration and system parameters are initialized as Table 1. The transverse range between the transmitting and the receiving platforms is 5 km, the transmitting platform is 1 km behind the receiving platform in longitude direction, and the heights of the two platforms are the same. Based on this configuration, the bistatic SAR imaging simulation verification is carried out. On the ground, nine points are arranged according to the receiver coordinate system, as shown in Figure 4.

Points B and C are the edge points of the field-focusing diagram, as is shown in Figure 5. Figure 5a is the imaging effect diagram without the azimuth variation correction, and Figure 5b is the imaging effect diagram with the two-dimensional variation correction. As can be seen from the focusing effect, the azimuth focusing effect is good after the azimuth space variation correction, while the edge points without the two-dimensional space variation correction show the defocus phenomenon.

Figure 6 shows the pulse profiles of points B and C along the azimuth and range directions. The blue solid line and the orange dotted line represent the result after and before the correction, respectively. As can be seen from the figure, the two-dimensional focusing effect of the target point is better after the space variation correction, but without the space variation correction, the azimuth is raised to the first sidelobe, the peak sidelobe ratio is lost, and the azimuth resolution is decreased. In addition, it can be seen from the PSLR and ISLR calculation results in Table 2 that the focusing effect of the edge points of the dual or multiple base SAR imaging algorithm in this configuration is consistent with that of the scene center point, which fully demonstrates the feasibility of this configuration.

In addition, the peak sidelobe ratio and integral sidelobe ratio index parameters of point A, point B, and point C of the scenic spots in the central field were calculated, as shown in the following table. It can be seen that the focusing effect of points A, B, and C is good.

4. Analysis of Cooperative Configuration Constraints for Aircrafts

4.1. Analysis of Bistatic SAR Imaging Architecture

The bistatic SAR imaging effect is mutually restricted with the spatial configuration of aircrafts, so it is necessary to select appropriate matching parameters between the imaging quality and the platform flight control to ensure a good configuration of the navigation control and to create conditions for high-precision imaging. Therefore, it is necessary to study the oblique angle of view, flight track, irradiation band width, and beam angle planning of the signals transmitted by the aircraft platform in cooperative flight to the target, as well as the accurate alignment and reception of the reflected signals of the target. The geometric model of the bistatic SAR forward-looking imaging of aircraft is shown in Figure 7.

Based on the simulation analysis of the imaging characteristics under dynamic conditions, the design method of the cooperative flight trajectory is proposed as follows:

First, determine the standard flight track of the radar signal receiving platform according to the flight destination.
According to the preset flight trajectory of the receiving platform, design the flight airspace range of the transmitting platform that meets the requirements of the imaging resolution.
Further, determine the actual flight envelope of the transmitting platform based on the flight airspace range of the transmitting platform in (2), considering the deviation of the actual navigation and control system of the transmitting platform, as well as the engineering constraints, such as the change of the range between the flight platform and the target, and the communication between the platforms.
Finally, plan the cooperative flight scheme and flight parameters of the receiving platform and the transmitting platform.

4.2. Resolution Analysis of Multi-Base SAR Imaging Based on Gradient Principle

Due to the influence of spatial relationship among target, transmitting platform, and receiving platform, the time resolution of aircraft bistatic SAR is firstly analyzed under given flight trajectory conditions in order to achieve two-dimensional high-resolution imaging.

According to the gradient principle, the forward-looking range resolution of the double base can be expressed as [24]

ρ_{r g} = \frac{c}{B (i_{R g} + i_{T g})}

(21)

i_{R_{g}}

and

i_{T_{g}}

, respectively, represent the ground projection of the range gradient of the transceiver platform,

c

represents the propagation speed of radar wave,

B

the signal bandwidth.

ρ_{r g}

the direction of range resolution, and its modulus

| ρ_{r g} |

is the size of range resolution, which reflects the ability of this configuration to distinguish targets in the horizontal plane. Then, on the premise of a certain signal bandwidth,

| ρ_{r g} |

is related to the bistatic projection angle

β_{g}

formed by the receiving platform and the target, as well as the projection component of the range gradient on the ground. All these factors are implicit and cannot be directly applied to the subsequent configuration design. Therefore, we convert the above formula to

| ρ_{r g} | = \frac{c}{B \sqrt{\sin^{2} φ_{r} + \sin^{2} φ_{t} + 2 \sin φ_{r} \sin φ_{t} \cos β_{g}}}

(22)

According to the above formula, under the condition of a single variable, the larger

β_{g}

is, the worse the ground range resolution will be. The smaller

φ_{t}

is, or the higher the flight height

H_{t}

of the aircraft is, the worse the ground range resolution will be.

According to the gradient principle, at the azimuth moment

t_{n}

, the azimuth resolution of the scene center point is [24]

ρ_{a g} = 1 / (T_{a} \nabla f_{g})

(23)

where

\nabla f_{g} = \frac{1}{λ} {\frac{1}{R_{t}} [V_{t g} - (V_{t g} \cdot i_{T g}) i_{T g}] + \frac{1}{R_{r}} [V_{r g} - (V_{r g} \cdot i_{R g}) i_{R g}]}

(24)

\nabla f_{g}

is the project of the Doppler gradient vector of the system on the ground.

V_{r g}

and

V_{t g}

represent the projection of the velocity vector of the receiver and the transmitter platform on the ground, respectively.

R_{r}

and

R_{t}

are the ranges between the transceiver platform and the target at the center of aperture, respectively.

T_{a}

the synthetic aperture time. The direction of

ρ_{a g}

is the direction of azimuth resolution, and its modulus

| ρ_{a g} |

is the magnitude of azimuth resolution.

According to Equation (24),

\nabla f_{g}

is composed of the system wavelength

λ

and the angular velocity vector of the transceiver platform, which is not only related to their modulus values but also related to their respective directions. For the convenience of analysis,

\nabla f_{g}

can be updated as

\nabla f_{g} = \nabla f_{t g} + \nabla f_{r g} = \frac{1}{λ R_{t}} (V_{t g} \cos θ_{t g}) i_{T w g} + \frac{1}{λ R_{r}} (V_{r g} \cos φ_{r g}) i_{R w g}

(25)

where

i_{R w g}

and

i_{T w g}

is the projection on the ground of the angular velocity unit vector of the transceiver platform, respectively.

V_{t g} \cos θ_{t g} ∕ λ R_{t}

represents the Doppler component mode

\nabla f_{t g}

provided by the transmitter, and

V_{r g} \cos φ_{r g} ∕ λ R_{r}

represents the Doppler component mode

\nabla f_{r g}

provided by the receiver. When the magnitude difference of the Doppler components provided by the transceiver platform is large, the azimuth resolution of the system can be approximately considered to be determined by the Doppler components of a single platform.

In order to further constrain the double-base resolution, a concept of resolving unit area is proposed, as shown in Figure 8. Because when the angle between range resolution direction and azimuth resolution direction is too small or too large, two-dimensional high-resolution images cannot be formed, the resolution unit area

S

can be defined by

S = \frac{| ρ_{r g} | \cdot | ρ_{a g} |}{\sin ϑ}

(26)

where

ϑ

is the included angle of range and azimuth resolution, that is, the included angle of direction

ρ_{r g}

and

ρ_{a g}

. It can be found that when the range and azimuth resolutions are orthogonal, the two-dimensional resolution formed by the same resolution becomes the best.

In the bistatic forward configuration, the direction of the two-dimensional resolution cannot remain orthogonal all the time. Therefore,

α

needs to be rationally planned. When the target is far away, the included angle

α

of the 2D resolution of all points in the scene is evenly distributed, which can be designed according to the center point of the scene. However, when the target distance is relatively close, the distribution of

α

is not uniform due to the large Doppler gradient vector direction span of the receiver. In order to improve this situation, it is necessary to increase the Doppler component of the transmitter as much as possible by reducing the angle of view

θ_{t}

of the transmitter when the target is close to the target.

5. Bistatic Imaging Simulation Based on Configuration Constraints for Aircraft Platform

5.1. Simulation Verification of Bistatic SAR Imaging in Different Configurations

In this section, the influence of the spatial distribution of a multi-base SAR platform on imaging resolution is analyzed through the bistatic SAR imaging simulation verification of a cooperative flight trajectory, and then the constraints of the platform configuration are obtained.

Bistatic SAR imaging simulation verification is carried out for different flight paths that meet the configuration requirements. The simulation parameters and index requirements of the system are shown in Table 1.

The resolution simulation in the case of a fixed longitude distance and variable transverse distance between two platforms.
The transverse distance of the two vehicles ranges from 5 km to 100 km, and a simulation is conducted at every 5 km step, as shown in Figure 9. Figure 9a shows how the range resolution varies with the distance between the receiving platform and the target, while Figure 9b shows the performance of the azimuth resolution. Figure 10 shows the relationship between the resolution and the distance between the receiving platform and the target.
It can be seen from the simulation that the larger the transverse distance between two vehicles, the better the final range resolution will be, because the larger the interval between two vehicles, the larger the ground projection component provided by the transmitting platform will be. The larger the transverse distance between the two vehicles, the worse the azimuth resolution of the final stage will be, because the transmitting platform with a short interval can provide a larger angular velocity component under the same velocity size. When the transverse distance between two platforms is small, the resolution angle is less than 30° in most cases, while when the transverse distance between the two platforms is large, the resolution angle remains 30° to 150° most of the time (except for the rear period), which meets the requirements of high-resolution imaging.
Resolution simulation in the case of a fixed transverse distance and variable longitude distance between two platforms.
In the collaborative flight of two platforms, the data transmission distance between the platforms is limited, and the platform spacing is generally within 20–40 km. Here, the transverse distance between the two platforms is fixed at 40 km for the simulation. The position of the receiving platform is fixed, and the distance between the transmitting platform and the receiving platform changes from −50 km to 50 km (negative indicates the backward direction). The simulation analysis is conducted at every 10 km step, and the results are shown in Figure 11. Among them, Figure 11a shows how the range resolution varies with the distance between the receiving platform and the target, while Figure 11b indicates the results of the azimuth resolution. Figure 12 shows the situation that the resolution-included angle in the imaging field varies with the distance between the receiving platform and the target.

Figure 9. Range resolution and azimuth resolution. (a) Range resolution, (b) azimuth resolution.

Figure 10. Angle of resolution.

Figure 11. Range resolution and Azimuth resolution. (a) Range resolution, (b) azimuth resolution.

Figure 12. Angle of resolution.

As can be seen from the above simulation, the resolution of the whole range becomes better when the transmitting platform is pulled back at a large distance, and the greater the distance is, the higher the resolution is. However, when the transmitting platform is pulled back, the azimuth resolution deteriorates, and the deterioration becomes more obvious at a larger distance. This is mainly because when the transmitting platform is pulled back at a larger distance, the included angle between the ground projection component provided by the transmitting platform and the ground projection component provided by the receiving platform decreases, while the angular velocity component provided by the transmitting platform decreases. In addition, after the transmitting platform is pulled back, the resolution angle of the scene center will decrease correspondingly. The larger the distance is, the smaller the resolution angle of the scene center will be.

On the contrary, when the transmitting platform is pulled forward, the azimuth resolution will improve, while the range resolution will deteriorate (especially at the ending stage). The greater the distance is pulled forward, the azimuth resolution will improve more obviously, while the final range resolution will deteriorate more obviously. This is because, when the transmitting platform moves forward, the azimuth velocity provided by the transmitting platform increases, while the included angle between the ground shadow component provided by the transmitting platform and the ground projection component provided by the receiving platform increases, and it will be even larger than 90° at the end of the stage, which will largely decrease the range resolution.

5.2. Configuration Constraints Optimization

By the simulation results that go through all the cooperative flight tracks in this paper, the changing tendencies of the range resolution, azimuth resolution, and resolution angle with transverse and longitude distances between platforms can be obtained.

Regarding the transverse distance, it is recommended to be as large as possible on the basis of ensuring the communication network. When the longitude distance of the platforms keeps constant and the transverse distance changes, it can be seen from the simulation that the larger the transverse distance is, the better the range and the azimuth resolution will be. When the transverse distance between the two platforms is small, the resolution angle is less than 30° most of the time, while when the transverse distance between the two platforms is large, the resolution angle will go between 30° and 150° in the most time (except at the ending period). This will satisfy the imaging requirements.

As to the longitude distance, it is recommended that the transmitting station be separated as large as −10~30 km. When the longitude distance of the platforms changes and the transverse distance is kept constant, the simulation shows that the range and the azimuth resolution have conflicting requirements for the distance and direction of the front and the rear platforms. With a preset resolution constraint, the distance between the platform and the target is analyzed from the flight segment of 80–40 km. For the range resolution, when the platform and the target are 80 km far away, the preset resolution can be met. When the target distance of the platform is 40 km, the distance between the front and the rear transmitting platforms is −50~30 km, which can meet the requirements. For the azimuth resolution, when the platform target distance is 80 km, the longitude distance of the transmitting platforms is −10~50 km, which can meet the resolution requirements; when the platform target distance is 40 km, it can meet the preset resolution requirements. In summary, it is recommended that the two transmitting stations be separated from −10 km to 30 km away.

6. Conclusions

Based on the construction of the geometric model of bistatic SAR imaging for an aircraft platform, this paper uses the extended NLCS algorithm to simulate and verify the effect of aircraft bistatic SAR imaging. A geometric model of an aircraft and target is established based on forward-looking bistatic SAR imaging, and a trajectory design method for platform cooperative flight is proposed. Based on the resolution gradient model, the relationship between the precision of the bistatic SAR imaging and the cooperative flight trajectory is studied, and the envelope range of the platform flight position for imaging the preset target point is given under a desirable resolution. The work of this paper will provide a theoretical basis for the implementation of the bistatic SAR imaging algorithm in engineering applications.

Author Contributions

Conceptualization, Y.C. and R.C.; methodology, Y.C.; software, Y.C.; validation, H.L., J.G. and Y.W.; formal analysis, H.L.; investigation, J.G.; resources, R.C.; data curation, Y.W.; writing—original draft preparation, Y.C.; writing—review and editing, J.Z.; visualization, Y.C.; supervision, R.C.; project administration, R.C.; funding acquisition, R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

Conflicts of Interest

The authors declare no conflict of interest.

References

Fan, C.Y.; He, S.S.; Guo, Q. Research status of radar forward-looking imaging technology. Electron. Opt. Control 2021, 28, 59–64. (In Chinese) [Google Scholar]
Huber, S.; de Almeida, F.Q.; Villano, M.; Younis, M.; Krieger, G.; Moreira, A. Tandem-L: A technical perspective on future spaceborne SAR sensors for Earth observation. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4792–4807. [Google Scholar]
Zhang, H.; Deng, Y.; Wang, R.; Li, N.; Zhao, S.; Hong, F.; Wu, L.; Loffeld, O. Spaceborne/stationary bistatic SAR imaging with TerraSAR-X as an illuminator in staring-spotlight mode. IEEE Trans. Geosci. Remote Sens. 2016, 54, 5203–5216. [Google Scholar]
Di Martino, G.; Di Simone, A.; Grassi, M.; Grasso, M.; Graziano, M.D.; Iodice, A.; Moccia, A.; Renga, A.; Riccio, D.; Ruello, G. Formation-Flying SAR Receivers in Far-from-Transmitter Geometry: Signal Model and Processing Scheme. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 2711–2714. [Google Scholar]
Servidia, P.A.; Espana, M. On Autonomous Reconfiguration of SAR Satellite Formation Flight With Continuous Control. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 3861–3873. [Google Scholar]
Lu, Z.; Wang, Y.; Kuang, H.; Huang, Y.; Wang, Z.; Zhu, Y. Formation Design and Accuracy Analysis of GEO-LEO Bistatic InSAR System. In Proceedings of the 2019 IEEE International Conference on Signal, Information and Data Processing (ICSIDP), Chongqing, China, 11–13 December 2019; pp. 1–4. [Google Scholar]
Li, Y.; Meng, Z.; Xing, M.; Bao, Z. Configuration study of missile-borne bistatic forward-looking SAR. In Proceedings of the 2014 IEEE China Summit & International Conference on Signal and Information Processing (ChinaSIP), Xi’an, China, 9–13 July 2014; pp. 184–188. [Google Scholar]
Meng, Z.; Li, Y.; Xing, M.; Bao, Z. Property analysis of bistatic forward-looking SAR with arbitrary geometry. J. Syst. Eng. Electron. 2016, 27, 111–127. [Google Scholar]
Ding, J.; Li, Y.; Quan, Y.; Wu, C.; Wang, Z.; Mei, H. Analysis of Diving Configuration Of Bistatic Forward-Looking SAR Based on Nonlinear Chirp Scaling Algorithm. In Proceedings of the 2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Xiamen, China, 26–29 November 2019; pp. 1–5. [Google Scholar]
Li, S.; Zhong, H.; Yang, C.; Song, H.; Zhao, R.; Cao, J.; Xu, X. Focusing Nonparallel-Track Bistatic SAR Data Using Modified Frequency Extended Nonlinear Chirp Scaling. IEEE Geosci. Remote Sens. Lett. 2020, 19, 1–5. [Google Scholar] [CrossRef]
Neo, Y.L.; Wong, F.; Cumming, I.G. A two-dimensional spectrum for bistatic SAR processing using series reversion. IEEE Geosci. Remote Sens. Lett. 2007, 4, 93–96. [Google Scholar] [CrossRef]
Wong, F.H.; Cumming, I.G.; Neo, Y.L. Focusing bistatic SAR data using the nonlinear chirp scaling algorithm. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2493–2505. [Google Scholar] [CrossRef]
Wang, R.; Loffeld, O.; Neo, Y.L.; Nies, H.; Walterscheid, I.; Espeter, T.; Klare, J.; Ender, J.H.G. Focusing bistatic SAR data in airborne/stationary configuration. IEEE Trans. Geosci. Remote Sens. 2009, 48, 452–465. [Google Scholar]
Eldhuset, K. A new fourth-order processing algorithm for spaceborne SAR. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 824–835. [Google Scholar]
Neo, Y.L.; Wong, F.H.; Cumming, I.G. Processing of azimuth-invariant bistatic SAR data using the range Doppler algorithm. IEEE Trans. Geosci. Remote Sens. 2007, 46, 14–21. [Google Scholar]
Vu, V.T.; Pettersson, M.I. Fast backprojection algorithms based on subapertures and local polar coordinates for general bistatic airborne SAR systems. IEEE Trans. Geosci. Remote Sens. 2015, 54, 2706–2712. [Google Scholar]
Li, Y.; Duan, C.; Wang, W.; Yang, X.; Yang, Z. Improvement of non-linear chirp-scaling algorithm for highly squint bistatic SAR data focusing. J. Eng. 2019, 2019, 6712–6717. [Google Scholar]
Liang, M.; Su, W.; Gu, H. Focusing high-resolution high forward-looking bistatic SAR with nonequal platform velocities based on keystone transform and modified nonlinear chirp scaling algorithm. IEEE Sens. J. 2018, 19, 901–908. [Google Scholar] [CrossRef]
Deng, H. Study on Imaging and Positioning of Missile-Borne Monostatic/Bistatic SAR with Curve Tracks. Ph.D. Thesis, Xidian University, Xi’an, China, 2018. (In Chinese). [Google Scholar]
Mei, H.W. Imaging Algorithm and Position Method Study on Bistatic/Multistatic SAR. Ph.D. Thesis, Xidian University, Xi’an, China, 2019. (In Chinese). [Google Scholar]
Meng, Z.Q.; Li, Y.C.; Xing, M.D.; Bao, Z. Imaging method for the extended scene of missile-borne bistatic forward-looking SAR. J. Xidian Univ. Electron. Sci. 2016, 43, 31–37. [Google Scholar]
Sun, G.; Xing, M.; Wang, Y.; Wu, Y.; Wu, Y.; Bao, Z. Sliding spotlight and TOPS SAR data processing without subapertur. IEEE Geosci. Remote Sens. Lett. 2011, 8, 1036–1040. [Google Scholar] [CrossRef]
Shao, Y.; Wang, R.; Deng, Y.K.; Liu, Y.; Chen, R.; Liu, G.; Balz, T.; Loffeld, O. Error analysis of bistatic SAR imaging and stereoscopy bistatic SAR. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4518–4543. [Google Scholar] [CrossRef]
Moccia, A.; Renga, A. Spatial resolution of bistatic synthetic aperture radar: Impact of acquisition geometry on imaging performance. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3487–3503. [Google Scholar] [CrossRef]

Figure 1. Imaging geometry model of bistatic forward-looking SAR.

Figure 2. Description of RCMC. (a) Bistatic imaging geometry model, (b) before LRCMC, (c) after LRCMC.

Figure 3. Flowchart of bistatic SAR imaging algorithm.

Figure 4. Layout of scenic spots.

Figure 5. Simulation of scenic spots. (a)Traditional azimuth pulse compression, (b) pulse compression after space variation correction.

Figure 6. Schematic diagram of azimuth and range pulse at point B. (a) Schematic diagram of azimuth and range pulse at point C. (b) Schematic diagram of azimuth and range pulse at point B.

Figure 7. Geometric model of bistatic forward-looking SAR imaging.

Figure 8. Schematic diagram of area of resolution element.

Table 1. Bistatic SAR forward-looking simulation parameters.

Parameter	Transmitter Down Angle	Receiver Down Angle	Transmitter Squint Angle	Receiver Squint Angle	Transmitter Initial Position	Receiver Initial Position
Value	15.23°	65.49°	−61.88°	89.85°	(1, 0.5, 15) km	(0, 0, 10) km

Table 2. Bistatic SAR forward-looking simulation parameters.

Parameter	Central Point (A)		Edge Point (B)		Edge Point (C)
Parameter	Azimuth Direction	Range Direction	Azimuth Direction	Range Direction	Azimuth Direction	Range Direction
PSLR (dB)	−13.20	−13.26	−13.07	−13.31	−13.15	−13.33
ISLR (dB)	−10.07	−9.87	−9.77	−10.21	−9.78	−10.08

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Chen, Y.; Chen, R.; Liu, H.; Guo, J.; Wang, Y.; Zhang, J. Research on Configuration Constraints of Airborne Bistatic SARs. Sensors 2022, 22, 6534. https://doi.org/10.3390/s22176534

AMA Style

Chen Y, Chen R, Liu H, Guo J, Wang Y, Zhang J. Research on Configuration Constraints of Airborne Bistatic SARs. Sensors. 2022; 22(17):6534. https://doi.org/10.3390/s22176534

Chicago/Turabian Style

Chen, Yidi, Renwen Chen, Hao Liu, Jiapeng Guo, Yujie Wang, and Junyi Zhang. 2022. "Research on Configuration Constraints of Airborne Bistatic SARs" Sensors 22, no. 17: 6534. https://doi.org/10.3390/s22176534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Research on Configuration Constraints of Airborne Bistatic SARs

Abstract

1. Introduction

2. Imaging Geometry Model of Bistatic SAR

3. Imaging Principle Model of Bistatic SAR

3.1. Analysis of Bistatic SAR Echo Signal

3.2. Bistatic SAR Imaging Algorithm

3.3. Simulation Verification of Bistatic SAR Imaging Algorithm

4. Analysis of Cooperative Configuration Constraints for Aircrafts

4.1. Analysis of Bistatic SAR Imaging Architecture

4.2. Resolution Analysis of Multi-Base SAR Imaging Based on Gradient Principle

5. Bistatic Imaging Simulation Based on Configuration Constraints for Aircraft Platform

5.1. Simulation Verification of Bistatic SAR Imaging in Different Configurations

5.2. Configuration Constraints Optimization

6. Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI