Next Article in Journal
A Practical Guide to Source and Receiver Locations for Surface Wave Transmission Measurements across a Surface-Breaking Crack in Plate Structures
Next Article in Special Issue
Noise Suppression for GPR Data Based on SVD of Window-Length-Optimized Hankel Matrix
Previous Article in Journal
Neural Activities Classification of Human Inhibitory Control Using Hierarchical Model
Previous Article in Special Issue
Target Doppler Rate Estimation Based on the Complex Phase of STFT in Passive Forward Scattering Radar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Focusing Bistatic Forward-Looking Synthetic Aperture Radar Based on an Improved Hyperbolic Range Model and a Modified Omega-K Algorithm

School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(17), 3792; https://doi.org/10.3390/s19173792
Submission received: 15 July 2019 / Revised: 26 August 2019 / Accepted: 30 August 2019 / Published: 1 September 2019
(This article belongs to the Special Issue Recent Advancements in Radar Imaging and Sensing Technology)

Abstract

:
For parallel bistatic forward-looking synthetic aperture radar (SAR) imaging, the instantaneous slant range is a double-square-root expression due to the separate transmitter-receiver system form. The hyperbolic approximation provides a feasible solution to convert the dual square-root expression into a single-square-root expression. However, some high-order terms of the range Taylor expansion have not been considered during the slant range approximation procedure in existing methods, and therefore, inaccurate phase compensation occurs. To obtain a more accurate compensation result, an improved hyperbolic approximation range form with high-order terms is proposed. Then, a modified omega-K algorithm based on the new slant range form is adopted for parallel bistatic forward-looking SAR imaging. Several simulation results validate the effectiveness of the proposed imaging algorithm.

1. Introduction

Synthetic aperture radar (SAR) attracts massive research enthusiasm among researchers due to its excellent ability to detect targets without the limitation of the external environment [1]. The penetration ability of SAR makes it irreplaceable compared with optical imaging, while it is challenging in traditional monostatic SAR to obtain excellent imaging performance in forward-looking imaging mode, which limits the application of SAR technology. To solve the problem, bistatic SAR has been widely used for forward-looking imaging due to its particular system configuration. The separate transmitter and receiver configuration provides extra advantages like reliable hiding power and system flexibility [2].
One-stationary bistatic SAR, as a special form of general bistatic SAR, was first studied for forward-looking imaging. Several methods have been proposed, such as the squint minimization [3,4], the keystone transform [5], and the ellipse model [2,4]. The Doppler frequency is decided by the moving transmitter or the moving receiver, which is similar to monostatic SAR. Then, the bistatic SAR was proposed where both the transmitter and the receiver are moving. The azimuth resolution is determined by both platforms. For bistatic forward-looking SAR, the difficulty of imaging algorithms lies in the solution of the two-dimensional spectrum because of its unique double-square-root form of echo signal expression [6,7]. Some basic studies of bistatic SAR were proposed to illustrate the advantages [6,8]. Compared with the monostatic situation, the principle of stationary phase (POSP) cannot be applied to solve the derivative zero point when performing azimuth Fourier transform. Several methods have been proposed to solve the problem. Loffeld’s bistatic formula (LBF) was proposed to solve the double-square-root expression [6]. Respective stationary points of the transmitter and receiver are obtained first to transform the double-square-root expression into Taylor expansion form. Then, the ultimate spectrum is solved based on the joint stationary point of the Taylor expression. The contributions of the transmitter and receiver are assumed to be the same, which leads to approximation errors. The extended Loffeld’s bistatic formula (ELBF) [9] and the modified Loffeld’s bistatic formula (MLBF) [10] were proposed later to improve the solution process of stationary points. These two methods assign different weights on the transmitter and receiver. However, all three LBF methods need to solve the stationary points three times, which leads to deduction complexity. The method of series reversion (MSR) [11] is a widely-used method for precisely solving those equations with series terms. In SAR imaging algorithms, Taylor expansion is regarded as a common operation, and MSR can be applied to solve the Fourier transform composed of Taylor expansion. However, it is still challenging to conduct imaging algorithm deduction due to the series form.
To simplify the solution of the spectrum, the hyperbolic approximation was utilized to transform the echo expression with the double-square-root form into the expression with the single-square-root form. In the first version, a parameter named the equivalent speed was defined [12]. In the improved version, two more parameters (the equivalent slant range and the equivalent squint angle) [13] were added in the hyperbolic function to approximate the range more accurately. Moreover, an improved hyperbolic approximation model with additional parameters was proposed for residual compensation [14]. However, considering the solution process in the methods mentioned above, the defined parameters are solved by setting the constant term, the linear term, and the quadratic term of the Taylor expansion of echo equal, which means the influence of the cubic term, the quartic term, and the remaining terms is ignored. In this article, we propose a new model to finish the hyperbolic approximation.
As for imaging algorithms, range Doppler (RD) imaging algorithms, chirp scaling (CS) imaging algorithms, back-projection (BP) imaging algorithms, and omega-K imaging algorithms based on the LBF spectrum, the MSR spectrum, and the hyperbolic approximation spectrum have been proposed in the past few years [9,14,15,16]. For RD imaging algorithms, it is too fundamental to handle the complex situation of bistatic SAR system. The calculation time consumption is a severe problem for real-time processing when applying BP imaging algorithms. For CS imaging algorithms, it is difficult for researchers to conduct formula derivation. Thus, the omega-K imaging algorithm is selected in this article to finish imaging.
To approximate the slant range more accurately, the cubic term and the quartic term are taken into account in this article. An equivalent hyperbolic range model is introduced first to lay the foundation of the imaging algorithm. The range error analysis is provided to demonstrate the approximation ability of the proposed range model immediately. Then, the modified omega-K imaging algorithm including the signal model and detailed processing steps are presented. Finally, some experimental simulations are given to prove the efficiency of the proposed algorithm.
This article is organized as follows. Section 2 gives the geometry of the bistatic forward-looking SAR and the equivalent hyperbolic range model corresponding to the bistatic system. Section 3 gives the detailed modified omega-K imaging algorithm. Simulation results are given in Section 4 to validate the proposed algorithm. Section 5 provides the conclusion.

2. Geometry and Equivalent Slant Range Model

The parallel bistatic forward-looking SAR system diagram in the Cartesian coordinate system and the derived equivalent slant range model are established first. Then, the analysis of range error based on the equivalent slant range model is provided.

2.1. Equivalent Slant Range Model

Figure 1 shows the geometry of parallel bistatic forward-looking SAR. The transmitter T and the receiver Rmove along the parallel red lines parallel to the x-axis. η p c is the synthetic aperture center time of the imaging scene. x c , y c , 0 is the location coordinate of the imaging center, and x p , y p , 0 is the location of an arbitrary target P in the imaging scene. R T c is the slant range between the transmitter and the target P at the phase center crossing time η p c , and R R c is the range between the receiver and the target P at η p c . The approximate forward-looking angle of the receiver is θ R , and the approximate squint angle of the transmitter is θ T . V T and V R represent the speed of the transmitter and the receiver, respectively. It is assumed that both the transmitter and the receiver can cover the imaging scene during the aperture synthesis.
The instantaneous slant ranges from the transmitter and the receiver to the target P are:
R T η = R T c 2 + V T 2 η η p c 2 2 R T c V T η η p c sin θ T , R R η = R R c 2 + V R 2 η η p c 2 2 R R c V R η η p c sin θ R ,
where η is the slow time.
Thus, the total range is:
R η = R T η + R R η .
It is challenging to solve the two-dimensional spectrum due to the double-square-root expression form of R η . The hyperbolic approximation [12] can be used to convert the double-square-root form to a single-square-root form by defining the equivalent speed and equivalent angle. Traditional hyperbolic approximation [12,13,14] ignored the high-order terms of the Taylor expansion of R η . To realize a more accurate compensation, an improved equivalent slant range with high-order terms is proposed. The range model is expressed as:
R e η = R e 2 + V e 2 η η p c 2 2 R e V e η η p c sin θ e + E ( η η p c ) 3 + F ( η η p c ) 4 ,
R η = 2 R e η ,
where R e , V e , and θ e are the new equivalent slant range at phase crossing time, the new equivalent speed, and the new equivalent squint angle. Compared with existing hyperbolic approximation algorithms, the proposed range model adds two additional high-order terms for range error compensation. To solve the unknown variables, we first expand Equations (1) and (3) into a fourth-order Taylor series at η = η p c . Then, we get:
R T η = R T c V T sin θ T η η p c + V T 2 cos θ T 2 2 R T c η η p c 2 + V T 3 sin θ T cos 2 θ T 2 R T c 2 η η p c 3 + V T 4 cos 2 θ T 5 sin 2 θ T 1 8 R T c 3 η η p c 4 ,
R R η = R R c V R sin θ R η η p c + V R 2 cos θ R 2 2 R R c η η p c 2 + V R 3 sin θ R cos 2 θ R 2 R R c 2 η η p c 3 + V R 4 cos 2 θ R 5 sin 2 θ R 1 8 R R c 3 η η p c 4 ,
R e η = R e V e sin θ e η η p c + V e 2 cos θ e 2 2 R e η η p c 2 + V e 3 sin θ e cos 2 θ e 2 R e 2 η η p c 3 + V e 4 cos 2 θ e 5 sin 2 θ e 1 8 R e 3 η η p c 4 + E ( η η p c ) 3 + F ( η η p c ) 4 .
Substituting Equations (5)–(7) into Equations (2) and (4) and letting the first five terms of Taylor expansion be equal, then we get:
R T c + R R c = 2 R e V T sin θ T + V R sin θ R = 2 V e sin θ e V T 2 cos θ T 2 2 R T c + V R 2 cos θ R 2 2 R R c = 2 V e 2 cos θ e 2 2 R e V T 3 sin θ T cos 2 θ T 2 R T c 2 + V R 3 sin θ R cos 2 θ R 2 R R c 2 = 2 V e 3 sin θ e cos 2 θ e 2 R e 2 + E V T 4 cos 2 θ T 5 sin 2 θ T 1 8 R T c 3 + V R 4 cos 2 θ R 5 sin 2 θ R 1 8 R R c 3 = 2 V e 4 cos 2 θ e 5 sin 2 θ e 1 8 R e 3 + F .
Solving the five equations in Equation (8), then we get:
R e = 1 2 R T c + R R c V e = A 2 + B θ e = arcsin A / V e E = C V e 3 sin θ e cos 2 θ e 2 R e 2 F = D V e 4 cos 2 θ e 5 sin 2 θ e 1 8 R e 3 ,
where:
A = V T sin θ T + V R sin θ R / 2 B = V T 2 cos 2 θ T R T c + V R 2 cos 2 θ R R R c R e / 2 C = V T 3 sin θ T cos 2 θ T 4 R T c 2 + V R 3 sin θ R cos 2 θ R 4 R R c 2 D = V T 4 cos 2 θ T 5 sin 2 θ T 1 16 R T c 2 + V R 4 cos 2 θ R 5 sin 2 θ R 1 16 R R c 2 .
At this point, all defined variables are solved. The range error analysis based on the new equivalent range model is presented next.

2.2. Range Error Analysis

To evaluate the proposed equivalent range model, an analysis of the range error based on an X-band bistatic SAR system is given. The simulated parameters are listed in Table 1. The results of the equivalent hyperbolic slant range error are shown in Figure 2.
Figure 2a is the approximation error of the traditional hyperbolic approximation range model [13], where the high-order terms are ignored. Figure 2b is the approximation error of the proposed hyperbolic range model. The constant term, the linear term, and the quadratic term in Equations (5)–(7) are used to solve the defined variables. Thus, the residual terms lead to the approximation slant range error. To prove that the proposed model can reduce the range error compared with the traditional model, we first give the expression of the traditional model and its corresponding Taylor expansion, which are:
R t η = R t 2 + V t 2 η η p c 2 2 R t V e η η p c sin θ t ,
R t η = R t V t sin θ t η η p c + V t 2 cos θ t 2 2 R t η η p c 2 + V t 3 sin θ t cos 2 θ t 2 R t 2 η η p c 3 + V t 4 cos 2 θ t 5 sin 2 θ t 1 8 R t 3 η η p c 4
where R t η , R t , V t , and θ t are the variables in traditional range model. The error in Figure 2a is the difference between the sum of the cubic terms, the quartic terms, and the residual terms in Equations (5) and (6) and the sum of the cubic term, the quartic term, and the residual term in Equation (12), while the error in Figure 2b is the difference between the sum of the residual terms in Equations (5) and (6) and the residual term in Equation (7). The error caused by the cubic and quartic terms is eliminated. From Figure 2, it can be found that the error in Figure 2a is up to 1.9 m, while the error in Figure 2b is less than 0.15 m. According to the parameters listed in Table 1, the approximation slant range error of the proposed model is much less than a range solution cell. Therefore, the proposed equivalent slant range model is more accurate than the traditional range model. The following imaging algorithm is derived based on the proposed range model.

3. Imaging Algorithm

According to the previous analysis, the improved hyperbolic approximation model can equal the true slant range better than traditional hyperbolic approximate models. In this section, a modified omega-K algorithm based on the improved equivalent range model is proposed for the parallel bistatic forward-looking SAR imaging.

3.1. Signal Model

Assume that a linear frequency-modulated signal is transmitted from the transmitter to the receiver. Then, the base-band echo signal of an arbitrary target P is given as:
S 1 t r , η = exp j π γ t r 2 R e η c exp j 4 π R e η λ
where γ is the range chirp rate, c is the light speed, λ is the wavelength, t r is the fast time, and η is the slow time. To simplify the expression and further derivation, the envelopes of the range and azimuth are ignored.
Transforming Equation (13) into the range-frequency azimuth-time domain yields:
S 2 ( f r , η ) = exp j π f r 2 γ exp j 4 π f r + f c c R e η
where f r is the frequency domain variable corresponding to t r and f c is the carrier frequency. From Equation (14), it can be easily found that the first exponential term is the range frequency modulation term. This term can be compensated by multiplying its complex conjugate in the range frequency domain. Thus, the first frequency modulation compensation function is:
H 1 F M f r , η = exp j π f r 2 γ .
Multiplying Equation (14) by Equation (15) yields:
S 3 ( f r , η ) = exp j 4 π f r + f c c R e η .
The exponential term in Equation (16) indicates the severe coupling between range and azimuth. To finish the phase focusing, a modified omega-K algorithm based on the signal model is presented.

3.2. Modified Omega-K Imaging Algorithm

To analyze the exponential term in Equation (16), Equation (3) is substituted into Equation (16) firstly. Then, we get:
S 4 f r , η = exp j 4 π f r + f c c [ R e 2 + V e 2 η η p c 2 2 R e V e η η p c sin θ e = + E η η p c 3 + F η η p c 4 ] .
Equation (17) shows that the signal consists of the traditional hyperbolic term and high-order terms. The traditional omega-K can handle the hyperbolic term well, but cannot handle the high-order terms. The first step of the omega-K algorithm is the compensation of the cubic term and the quartic term. Variable substitution is performed on Equation (17), and then, we get:
S 5 k r , X = exp j k r R e 2 + X X p c 2 2 R e X X p c sin θ e = + E V e 3 X X p c 3 + F V e 4 X X p c 4 ]
where k r = 4 π f r + f c c is the wavenumber, X = V e η , and X p c = V e η p c . Then, we get R e η = R e X .
Transforming Equation (18) into two-dimensional wavenumber domain yields:
S 6 k r , k x = S 5 k r , X exp j k x X d X = exp j k r R e X exp j k x X d X = exp j ϕ k r , k x , X d X
where k x = 2 π f a V e , f a is the azimuth frequency, and:
ϕ k r , k x , X = k r R e 2 + X X p c 2 2 R e X X p c sin θ e + E V e 3 X X p c 3 + F V e 4 X X p c 4 + k x X .
To solve Equation (19), the stationary phase point of ϕ k r , k x , X should be obtained firstly. However, the existence of high-order terms complicates the solution process. For further analysis, the phase is first rewritten as:
ϕ k r , k x , X = ϕ t k r , k x , X + k r E V e 3 X X p c 3 + F V e 4 X X p c 4
where ϕ t k r , k x , X is the traditional phase term. It is widely accepted that if the phase error is smaller than π / 4 [1], the imaging performance will not be affected much by the approximation. The phase error simulation is given in Figure 3.
From Figure 3, it can been seen that all absolute phase errors are less than π / 4 . Thus, the stationary phase point of ϕ t k r , k x , X is regarded as the approximate stationary phase point of ϕ k r , k x , X . The approximate stationary phase point of ϕ k r , k x , X is:
X * = k x R e sin θ e k r 2 k x 2 + R e sin θ e + X p c ,
where X * is only a designation of the solution and ( * ) is not an operator.
Substituting Equation (22) in Equation (19) and applying POSP yield the two-dimensional wavenumber domain signal as:
S 7 k r , k x = exp j k r 2 k x 2 R e cos θ e j k x R e sin θ e + X p c = j k r E V e 3 X * X p c 3 + F V e 4 X * X p c 4 .
The cubic term and quartic term in Equation (23) can be easily compensated by multiplying its conjugate form. Therefore, the high-order filter is:
H 2 k r , k x = exp j k r E V e 3 X * * X p c 3 + F V e 4 X * * X p c 4
where X * * is the value of X * at the reference range and ( * * ) is not an operator.
Multiplying Equations (23) and (24), we get the compensated signal for the further omega-K imaging algorithm. The signal is:
S 8 k r , k x = exp j k r 2 k x 2 R e cos θ e j k x R e sin θ e + X p c .
A two-step omega-K is performed on Equation (25) to finish the imaging focusing.
The first step is the bulk focusing. A reference function is designed based on the reference range to finish coarse focusing. This filter can compensate the phase of signals of those points at the reference range. The reference function is:
H r f k r , k x = exp j k r 2 k x 2 R r e f cos θ e + j k x R r e f sin θ e + X p c .
Multiplying Equations (25) and (26) gets:
S 9 k r , k x = exp j k r 2 k x 2 cos θ e R e R r e f j k x sin θ e R e R r e f .
After bulk focusing, the residual phase at the reference range is removed. However, the residual phase of points not at the reference range remains. Moreover, the phase contains coupling terms between range and azimuth. For precise focusing of all points, the Stolt interpolation function is given as:
k y = k r 2 k x 2 cos θ e + k x sin θ e .
After Stolt interpolation, the resampled signal becomes:
S 10 k r , k x = exp j k y R e R r e f .
From Equation (29), it is evident that the coupling between range and azimuth has been removed. The phase is a linear function of k y . Then, the inverse fast Fourier transform is implemented on Equation (29) to complete imaging.
According to the analysis mentioned above, the whole imaging process is shown in Figure 4.
The specific steps are as follows:
(1)
Performing range fast Fourier transform (FFT) on SAR data gets S 2 ( f r , η ) .
(2)
Multiplying Equation (15) and S 2 ( f r , η ) gets S 3 ( f r , η ) .
(3)
Performing azimuth fast Fourier transform (FFT) on S 3 ( f r , η ) gets S 7 ( k r , k x ) .
(4)
Multiplying Equation (24) and S 7 ( k r , k x ) gets S 8 ( k r , k x ) .
(5)
Multiplying Equation (26) and S 8 ( k r , k x ) gets S 9 ( k r , k x ) .
(6)
Performing Stolt interpolation on S 9 ( k r , k x ) gets S 10 ( k r , k x ) .
(7)
Performing 2D-IFFT on S 10 ( k r , k x ) gets output SAR focusing results.

4. Simulation Results

In this section, to demonstrate the effectiveness of the proposed imaging algorithm, experimental simulations of parallel bistatic forward-looking SAR are carried out. The system parameters are listed in Table 1. Four points at different locations were chosen to compare the imaging performance. They were P 0 ( 0 , 0 ) , P 1 ( 0 , 500 ) , P 2 ( 200 , 0 ) , and P 3 ( 200 , 500 ) . The unit of the coordinates is meters. The omega-K imaging algorithm based on the traditional three-parameters hyperbolic range model [13] was selected as the reference.
Figure 5 is the comparison of the overall imaging performance before geometric correction. Figure 5a is the result of the traditional imaging algorithm, and Figure 5b is the result of the proposed imaging algorithm. In Figure 5a, although the four points can be successfully focused, the quality of the right two points has distortion. In contrast, Figure 5b shows that the proposed algorithm achieves a better focus quality on the right two points than the traditional algorithm.
To observe the imaging performance more intuitively, the sub-images of the four points extracted from Figure 5 are given by Figure 6. Figure 6a–c presents the imaging results of P 0 , P 2 , and P 3 achieved by the traditional hyperbolic range model given in [13], respectively. Figure 6e,f shows the imaging quality of the three targets obtained by the proposed modified omega-K imaging algorithm. From Figure 6a,d, both algorithms can obtain an excellent focusing quality of the scene center P 0 . For the omega-K algorithm, the scene center is always chosen as the reference point to perform bulk focusing. For the points away from the center ( P 2 and P 3 ), it is evident that the proposed algorithm performs much better than the traditional algorithm. For further analysis, the azimuth impulse response of the farthest point P 3 is given in Figure 7. Table 2 gives out the peak sidelobe ratio (PSLR) and the integrated sidelobe ratio (ISLR) of targets P 3 .
Figure 7a is achieved by the traditional hyperbolic omega-K algorithm. Figure 7b is achieved by the proposed hyperbolic omega-K algorithm. Compared with the traditional omega-K algorithm, the proposed omega-K algorithm can improve the performance of the azimuth impulse response. The objective image quality values demonstrated the effectiveness of the proposed omega-K algorithm.

5. Conclusions

In this article, an improved hyperbolic range model was proposed to deal with the particular form of the echo of bistatic forward-looking SAR. The modified omega-K imaging algorithm based on the hyperbolic range model was used to finish focusing. The high-order terms were taken into account to reduce the range approximation error. Extra phase compensation benefited the focusing of the omega-K algorithm. Compared with the range model without high-order compensation terms, the proposed method showed the effectiveness of imaging quality by simulation results.

Author Contributions

Conceptualization, W.S. and H.G.; methodology, C.W.; software, C.W.; validation, all authors; formal analysis, C.W.; investigation, J.Y.; writing–original draft preparation, C.W.; writing—review and editing, W.S. and J.Y.; supervision, H.G.; project administration, W.S.; funding acquisition, W.S., H.G. and J.Y.

Funding

This research was funded by the National Natural Science Foundation of China under Grant Numbers 61471198 and 61671246 and the Natural Science Foundation of Jiangsu Province under Grant Numbers BK20160847 and BK20170855.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  2. Zhong, H.; Yan, A.; Sun, M.; Zhang, X.; Zhang, J. An improved ENLCS algorithm based on a quadratic ellipse model for high-resolution one-stationary bistatic SAR. Remote. Sens. Lett. 2018, 9, 867–876. [Google Scholar] [CrossRef]
  3. Ma, C.; Gu, H.; Su, W.; Zhang, X.; Li, C. Focusing one-stationary bistatic forward-looking synthetic aperture radar with squint minimisation method. IET Radar Sonar Navig. 2015, 9, 927–932. [Google Scholar] [CrossRef]
  4. Liang, M.; Su, W.; Ma, C.; Gu, H.; Yang, J. Focusing one-stationary bistatic forward-looking synthetic aperture radar based on squint minimization and modified nonlinear chirp scaling algorithm. Int. J. Remote Sens. 2018, 39, 7830–7845. [Google Scholar] [CrossRef]
  5. Wu, J.; Li, Z.; Huang, Y.; Yang, J.; Yang, H.; Liu, Q.H. Focusing Bistatic Forward-Looking SAR with Stationary Transmitter Based on Keystone Transform and Nonlinear Chirp Scaling. IEEE Geosci. Remote Sens. Lett. 2014, 11, 148–152. [Google Scholar] [CrossRef]
  6. Loffeld, O.; Nies, H.; Peters, V.; Knedlik, S. Models and useful relations for bistatic SAR processing. IEEE Trans. Geosci. Remote Sens. 2004, 42, 2031–2038. [Google Scholar] [CrossRef]
  7. Walterscheid, I.; Ender, J.H.G.; Brenner, A.R.; Loffeld, O. Bistatic SAR Processing and Experiments. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2710–2717. [Google Scholar] [CrossRef]
  8. Wang, R.; Loffeld, O.; Ul-Ann, Q.; Nies, H.; Ortiz, A.M.; Samarah, A. A Bistatic Point Target Reference Spectrum for General Bistatic SAR Processing. IEEE Geosci. Remote Sens. Lett. 2008, 5, 517–521. [Google Scholar] [CrossRef]
  9. Wang, R.; Loffeld, O.; Neo, Y.L.; Nies, H.; Dai, Z. Extending Loffeld’s bistatic formula for the general bistatic SAR configuration. IET Radar Sonar Navig. 2010, 4, 74–84. [Google Scholar] [CrossRef]
  10. Ma, C.; Gu, H.; Su, W.; Li, C.; Chen, J. Focusing bistatic forward-looking synthetic aperture radar based on modified Loffeld’s bistatic formula and chirp scaling algorithm. J. Appl. Remote Sens. 2014, 8, 083586. [Google Scholar] [CrossRef]
  11. 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]
  12. Bamler, R.; Boerner, E. On the use of numerically computed transfer functions for processing of data from bistatic SARs and high squint orbital SARs. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Seoul, Korea, 29 July 2005. [Google Scholar]
  13. Qiu, X.; Hu, D.; Ding, C. Focusing Bistaitc Images use RDA based on Hyperbolic Approximating. In Proceedings of the CIE International Conference on Radar, Shanghai, China, 16–19 October 2006. [Google Scholar]
  14. Ma, C.; Gu, H.; Su, W.; Li, C. Focusing bistatic forward-looking synthetic aperture radar based on modified hyperbolic approximating. Acta Phys. Sin. 2014, 63, 028403. [Google Scholar]
  15. Chen, S.; Yuan, Y.; Zhang, S.; Zhao, H.; Chen, Y. A New Imaging Algorithm for Forward-Looking Missile-Borne Bistatic SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 543–1552. [Google Scholar] [CrossRef]
  16. Feng, D.; An, D.; Huang, X. An Extended Fast Factorized Back Projection Algorithm for Missile-Borne Bistatic Forward-Looking SAR Imaging. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2724–2734. [Google Scholar] [CrossRef]
Figure 1. Geometry of forward-looking bistatic SAR.
Figure 1. Geometry of forward-looking bistatic SAR.
Sensors 19 03792 g001
Figure 2. Approximation error of the bistatic slant range. (a) Approximation error of the traditional range model. (b) Approximation error of the proposed range model.
Figure 2. Approximation error of the bistatic slant range. (a) Approximation error of the traditional range model. (b) Approximation error of the proposed range model.
Sensors 19 03792 g002
Figure 3. Phase error simulation.
Figure 3. Phase error simulation.
Sensors 19 03792 g003
Figure 4. Flowchart of modified omega-K.
Figure 4. Flowchart of modified omega-K.
Sensors 19 03792 g004
Figure 5. Imaging results. (a) Imaging results of the traditional hyperbolic omega-K algorithm. (b) Imaging results of the proposed hyperbolic omega-K algorithm.
Figure 5. Imaging results. (a) Imaging results of the traditional hyperbolic omega-K algorithm. (b) Imaging results of the proposed hyperbolic omega-K algorithm.
Sensors 19 03792 g005
Figure 6. Imaging results. (a) Imaging result of P 0 by the traditional algorithm. (b) Imaging result of P 2 by the traditional algorithm. (c) Imaging result of P 3 by the traditional algorithm. (d) Imaging result of P 0 by the proposed algorithm. (e) Imaging result of P 2 by the proposed algorithm. (f) Imaging result of P 3 by the proposed algorithm.
Figure 6. Imaging results. (a) Imaging result of P 0 by the traditional algorithm. (b) Imaging result of P 2 by the traditional algorithm. (c) Imaging result of P 3 by the traditional algorithm. (d) Imaging result of P 0 by the proposed algorithm. (e) Imaging result of P 2 by the proposed algorithm. (f) Imaging result of P 3 by the proposed algorithm.
Sensors 19 03792 g006
Figure 7. Azimuth impulse response of P 3 . (a) Traditional hyperbolic omega-K algorithm. (b) Proposed hyperbolic omega-K algorithm.
Figure 7. Azimuth impulse response of P 3 . (a) Traditional hyperbolic omega-K algorithm. (b) Proposed hyperbolic omega-K algorithm.
Sensors 19 03792 g007
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersValuesParametersValues
Carrier frequency9 GHzTransmitter center slant range4300 m
Pulse duration μ s Transmitter squint angle 7
Bandwidth200 MHzReceiver center slant range3600 m
Sampling frequency300 MHzReceiver forward-looking angle 33
Pulse repetition frequency1 kHzSensor speed200 m/s
Table 2. Image quality parameters of P 3 . PSLR, peak sidelobe ratio; ISLR, integrated sidelobe ratio.
Table 2. Image quality parameters of P 3 . PSLR, peak sidelobe ratio; ISLR, integrated sidelobe ratio.
TargetsPSLR (dB)ISLR (dB)
AzimuthRangeAzimuthRange
Traditional omega-K algorithm−1.705---
Proposed omega-K algorithm−12.87−13.33−8.86−9.9558

Share and Cite

MDPI and ACS Style

Wang, C.; Su, W.; Gu, H.; Yang, J. Focusing Bistatic Forward-Looking Synthetic Aperture Radar Based on an Improved Hyperbolic Range Model and a Modified Omega-K Algorithm. Sensors 2019, 19, 3792. https://doi.org/10.3390/s19173792

AMA Style

Wang C, Su W, Gu H, Yang J. Focusing Bistatic Forward-Looking Synthetic Aperture Radar Based on an Improved Hyperbolic Range Model and a Modified Omega-K Algorithm. Sensors. 2019; 19(17):3792. https://doi.org/10.3390/s19173792

Chicago/Turabian Style

Wang, Chenchen, Weimin Su, Hong Gu, and Jianchao Yang. 2019. "Focusing Bistatic Forward-Looking Synthetic Aperture Radar Based on an Improved Hyperbolic Range Model and a Modified Omega-K Algorithm" Sensors 19, no. 17: 3792. https://doi.org/10.3390/s19173792

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

Article Metrics

Back to TopTop