Research Article  Open Access
MIMOBased ForwardLooking SAR Imaging Algorithm and Simulation
Abstract
Multipleinput multipleoutput (MIMO) radar imaging can provide higher resolution and better sensitivity and thus can be applied to targets detection, recognition, and tracking. Missileborne forwardlooking SAR (MFLSAR) is a new and special MIMO radar mode. It has advantage of twodimensional (2D) imaging capability in forward direction over monostatic missileborne SAR and airborne SAR. However, it is difficult to obtain accurate 2D frequency spectrum of the target echo signal due to the high velocity and descending height of this platform, which brings a lot of obstacles to imaging algorithm design. Therefore, a new imaging algorithm for MFLSAR configuration based on the method of series reversion is proposed in this paper. This imaging method can implement range compression, secondary range compression (SRC), and range cell migration correction (RCMC) effectively. Finally, some simulations of point targets and comparison results confirm the efficiency of our proposed algorithm.
1. Introduction
Radar imaging is an emerging technology which can provide a high resolution radar image of targets in long distance independent of weather conditions and sunlight illumination, and it is being more and more widely used in military and civil fields [1–5]. Multipleinput multipleoutput (MIMO) radar is an antenna system which transmits multiple probing signals via its antennas and scattering signals are received by distributed receiving antennas [6–9]. Application of MIMO radar to radar imaging could provide higher resolution and better sensitivity and thus can be applied to detection, recognition, and tracking of targets [10–12].
We consider a new MIMO radar system, that is, missileborne multistatic synthetic aperture radar (MFLSAR) as shown in Figure 1, in which transmitting signals are transmitted simultaneously and the reflected signals at receivers can be processed independently. This configuration can be applied to missile precision terminal guidance because of its advantages of 2D imaging ability in forward direction over monostatic SAR [13, 14]. Missile precision terminal guidance could conduct good performance in targets detection, recognition, orientation tracking, and attacking [15, 16]. During the moving period, transmitter and receiver cooperate with each other, and transmitter irradiates the imaging area with some squint angle, while receiver receives the target echo signals in forwardlooking mode. Without loss of generality and for the convenience of description, missileborne bistatic forwardlooking SAR (MBFLSAR) configuration is considered in this paper.
Effective and efficient imaging in this configuration is an important technology for missile precision terminal guidance; therefore, imaging algorithm for MBFLSAR is necessary and essential. Many investigations on SAR imaging have been published and some appealing approaches have been suggested. RangeDoppler algorithm (RDA) [17–19] and chirp scaling algorithm (CSA) [20] are conventional monostatic algorithms which are applied to those configurations that have azimuthinvariant property. For bistatic SAR imaging, some other methods have been proposed. Soumekh derived the 2D spectrum expression of bistatic SAR and proposed a new RMA for bistatic configuration through two approximations in the derivation [21–23]. But this algorithm is limited to the bistatic parallel configuration because of identical velocity vectors of transmitter and receiver. Neo et al. and Davidson et al. proposed a nonlinear CS algorithm in [24, 25], but it neglected spacevariance of Doppler frequency in range direction and influence of secondary range compression (SRC) on imaging. A simple operator named “dip move out” (DMO) is used to transform a bistatic survey into a monostatic one in [26, 27]. The initial dataset is used to convolve with the socalled smile, a short operator, and the output is the equivalent monostatic data. Unfortunately, the socalled smile is spacevariant both in range direction and in azimuth direction, and the changes of wave number in range direction are also neglected.
An imaging algorithm based on the method of series reversion [28, 29] for MBFLSAR is suggested in this paper. By unfolding the Doppler frequency, the stationary phase point expansion coefficients are obtained, which makes it easy to get the 2D frequency spectrum of the target echo signal. Then imaging method can be obtained based on the 2D frequency spectrum. Range compression is implemented in the range frequency domain, and SRC can be finished in the 2D frequency domain. The RCM can be corrected in the rangeDoppler domain and finally the imaging result can be got through azimuth compression.
The rest of this paper is organized as follows. Geometric configuration and signal model of MBFLSAR are introduced in Section 2. Also, the difficulty in directusing the principle of stationary phase (POSP) is discussed. In Section 3, the derivation of 2D frequency spectrum of the target echo signal is given based on the accurate approximation of the slant range history, in which the highorder terms cannot be ignored owing to the presence of high velocities and accelerations. The proposed imaging algorithm is described in Section 4 and numerical simulations are achieved to show the effectiveness of our proposed method, and the results are given in Section 5. Finally, conclusions are drawn in Section 6.
2. Geometric Configuration and Signal Model
Figure 2 depicts geometric configuration of MBFLSAR, in which both transmitter and receiver travel curvilinear descending motion. is the origin of coordinates, is the projection of transmitter when the slow time , and . Transmitter is moving with the original velocity vectors and acceleration vectors . And the ones of receiver are and , respectively. There exists an included angle between the two planes in which the two platforms travel. The transmitter keeps illuminating the imaging area in sidelooking mode, and the receiver looks in forward direction. The receiver antenna gets the echo signal reflected from the imaging area which the receiver moves towards. Assume that transmitter and receiver are at and when , and velocity vectors of them at the moment are and , respectively, where , , and denote the velocity components of in the directions , , and , respectively; and , , and are the ones of , respectively.
The bistatic slant range at time instant can be obtained as with where , , and represent the locations of transmitter in directions , , and , respectively; and , , and are the ones of receiver, respectively.
Suppose that the transmitted waveform is the linear frequency modulation (LFM), and scattering from to the receiver can be written as where is the range time, is the wavelength, is the speed of light, and is the chirp rate. and are the range and azimuth envelopes, respectively.
By applying a range fast Fourier transform (FFT) to (3), we can obtain where represents the range frequency.
The 2D frequency spectrum of the echo signal can be obtained by an azimuth FFT, expressed as where represents the azimuth frequency, and the phase in (5) can be written as
By applying the POSP, we have with
If , we can get the stationary point , but we can see that it is very difficult to determine the stationary point from (7), so 2D frequency spectrum of the echo signal cannot be obtained through direct POSP. Some other approaches should be taken to get 2D frequency spectrum of the echo signal.
3. Derivation of 2D Frequency Spectrum Based on Series Reversion
According to [28], it needs to take efficient approximation of the slant range before using the method of series reversion. Because of the high velocity and acceleration in the missile platform, it needs to keep the terms up to the thirdorder with
To apply the method of series reversion, the linear range cell migration (LRCM) should be removed firstly. Let , we have
Then through the method of series reversion, the expression of the stationary phase point can be expressed by unfolded Doppler frequency as where , .
Having obtained the stationary phase point, we should reintroduce the LRCM term to calculate 2D frequency spectrum of the echo signal. By (5), (6), and (12), 2D frequency spectrum of the echo signal is written as where
It can be seen that the coefficient of the thirdterm in the slant range is included in , which indicates the presence of the highorder terms introduced by high velocity and acceleration of this configuration into 2D frequency spectrum. The imaging algorithm can be designed based on 2D frequency spectrum, which will be discussed in the next section.
4. Imaging Algorithm for MBFLSAR
To design imaging algorithm efficiently, the phase term in (14) should be decomposed using Taylor series firstly because of the range/azimuth frequency coupling in 2D frequency spectrum, and the series expansions are as follows: with
The term in (16) represents the range compression term which is independent of the azimuth frequency . Thus the data can be rangecompressed in the range frequency azimuth time domain or in the 2D frequency domain alternatively. The term in (17) denotes the azimuth compression term, dependent only on and used in the azimuth compression. It should be implemented after range compression and RCMC in the range time azimuth frequency domain as coefficients are range variant. The term in (18) indicates the RCM term which is linearly dependent on the range frequency . It needs to be compensated because of the coupling between range and azimuth. Similar to the azimuth compression term, this term should be implemented in the range time azimuth frequency domain. Note that is the linear phase representing the location of target point; the energy of the target will be focused within the corresponding range cell after RCMC. It should be paid attention to in the imaging process. The term in (19) is the SRC term which indicates the coupling between and . It may cause significant degradation in the imaging resolution if uncompensated and is compensated in the 2D frequency domain. The last residual term in (20) is independent of the range frequency or the azimuth frequency. It has no effect on the imaging focus and thus can be neglected in MBFLSAR configuration.
Consequently, the main procedure of the proposed algorithm follows the following steps.(1)Use 2D FFT to transform the echo signal to 2D frequency domain as shown in (13).(2)Design the match filtering function in the 2D frequency domain to compensate the range compression term, the SRC term, and the residual term. The match filtering function can be obtained as follows: (3)Apply range inverse fast Fourier transform (IFFT) to convert the signal to the range time azimuth frequency domain and correct the RCM term.(4)Finish azimuth compression by multiplying the range time azimuth frequency domain signal with the azimuth match filtering function, which can be obtained as follows: (5)Perform azimuth IFFT to obtain a wellfocused SAR image.
5. Simulations and Results
In this section, several examples are provided to illustrate the performance of the proposed method. Because of the complex configuration, the highorder terms in the slant range cannot be ignored. If inappropriately approximated, the imaging focusing resolution may be affected significantly. Some simulations are conducted to compare the influences of keeping the terms up to the quadratic term (quadratic range model) and the higher order term (proposed range model). Parameters used in the simulations are listed in Table 1. The imaging results of center point and fringe point in the imaging area using the proposed range model are given in Figures 3 and 4. In order to confirm the necessity and validity using the higher order range model in this special configuration, some simulations are also conducted to compare the differences using the quadratic range model and the proposed range model, and the results are shown in Figures 5 and 6 (here, the center point is selected in comparison).

(a) 2D impulse response
(b) Azimuth impulse response
(c) Range impulse response
(a) 2D impulse response
(b) Azimuth impulse response
(c) Range impulse response
(a) Quadratic range model
(b) Proposed range model
(a) Quadratic range model
(b) Proposed range model
Figure 3 depicts the results of the center point using the proposed range model. Figure 3(a) is the 2D impulse response; Figures 3(b) and 3(c) are the azimuth impulse response and the range impulse response, respectively. The 2D impulse response of the fringe target point using the proposed range model is shown in Figure 4(a), and the azimuth impulse response and the range impulse response are given in Figures 4(b) and 4(c), respectively. We can see from Figures 3 and 4 that both center point and fringe point are well focused in range direction and azimuth direction. For the complexity of the MBFLSAR configuration, if the slant range is not approximated accurately enough, the phase term in the 2D frequency spectrum will not be compensated completely in the azimuth direction. It is easy to find that the image of the point target suffers from distortion and mainlobe broadening when the slant range keeps the terms up to the quadratic term in Figures 5(a) and 6(a). On the other hand, if we are using the proposed range model, the imaging focusing performance can be obviously improved, as given in Figures 5(b) and 6(b). It is obvious that the image is well focused in the azimuth direction, which confirms it necessary to approximate appropriately in the MBFLSAR configuration.
Table 2 lists the image quality parameters of the center point and fringe point using the proposed method to further illustrate the validity of the proposed algorithm. PSLR and ISLR represent peak sidelobe ratio and integrated sidelobe ratio, respectively. Besides, the theoretical PSLR and ISLR are −13.26 dB and −9.8 dB, respectively. We can see that the image quality parameters of both the center point and fringe point using our proposed method are all close to the theoretical ones, which indicate satisfactory imaging results and further validate the effectiveness and feasibility.

6. Conclusions
In this paper, we developed an imaging algorithm for MBFLSAR configuration, a special MIMO radar mode. The MBFLSAR imaging geometric model is established and the signal model is analyzed firstly. The highorder term in the slant range cannot be neglected because of the presence of high velocity and acceleration and the addition of two squareroot terms. Also, the difficulty in directusing the POSP in this configuration is discussed. Then 2D frequency spectrum is deduced using the method of series reversion based on the appropriate approximation of the slant range. Based on the derived spectrum, the corresponding imaging method is developed. The highorder term introduced from this configuration is compensated, together with range compression, at the beginning of the algorithm. RCMC is finished in the range time azimuth frequency domain as the rangevariant coefficients. The focused SAR image can be obtained through the azimuth compression. The feasibility and efficiency of the proposed approach are validated with the simulation experiments.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under Grants 61001211 and 61303035, by the Fundamental Research Funds for the Central Universities (K5051202016), and by the Science Foundation for Navigation (20110181004).
References
 I. G. Cumming and F. H. Wong, Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood, Mass, USA, 2005.
 D. L. Mensa, High Resolution Radar CrossSection Imaging, Artech House, Boston, Mass, USA, 1991.
 S. A. Hovanessian, Introduction to Synthetic Array and Imaging Radar, Artech House, Norwood, Mass, USA, 1980.
 R. Kapoor, A. Banerjee, G. A. Tsihrintzis, and N. Nandhakumar, “UWB radar detection of targets in foliage using alphastable clutter models,” IEEE Transactions on Aerospace and Electronic Systems, vol. 35, no. 3, pp. 819–834, 1999. View at: Publisher Site  Google Scholar
 G. Wang, X. Xia, and V. C. Chen, “Radar imaging of moving targets in foliage using multifrequency multiaperture polarimetric SAR,” IEEE Transactions on Geoscience and Remote Sensing, vol. 41, no. 8, pp. 1755–1764, 2003. View at: Publisher Site  Google Scholar
 E. Fishler, A. Haimovich, R. Blum, D. Chizhik, L. Cimini, and R. Valenzuela, “MIMO radar: An idea whose time has come,” in Proceedings of the IEEE Radar Conference, pp. 71–78, April 2004. View at: Google Scholar
 E. Fishler, A. Haimovich, R. Blum, L. Cimini, D. Chizhik, and R. Valenzuela, “Performance of MIMO radar systems: advantages of angular diversity,” in Proceedings of the 38th Asilomar Conference on Signals, Systems and Computers, pp. 305–309, November 2004. View at: Google Scholar
 R. Sharma, “Analysis of MIMO radar ambiguity functions and implications on clear region,” in Proceedings of the IEEE International Radar Conference (RADAR ’10), pp. 544–548, Washington, DC, USA, May 2010. View at: Publisher Site  Google Scholar
 J. Li and P. Stoica, MIMO Radar Signal Processing, Wiley, New York, NY, USA, 2009.
 W. Roberts, P. Stoica, J. Li, T. Yardibi, and F. A. Sadjadi, “Iterative adaptive approaches to MIMO radar imaging,” IEEE Journal on Selected Topics in Signal Processing, vol. 4, no. 1, pp. 5–20, 2010. View at: Publisher Site  Google Scholar
 J. Li, G. Liao, and H. Griffiths, “Bistatic MIMO radar spacetime adaptive processing,” in Proceedings of the IEEE Radar Conference: in the Eye of the Storm (RadarCon ’11), pp. 498–501, Kansas, Mo, USA, May 2011. View at: Publisher Site  Google Scholar
 D. W. Bliss and K. W. Forsythe, “Multipleinput multipleoutput (MIMO) radar and imaging: degrees of freedom and resolution,” in Proceedings of the 37th Asilomar Conference on Signals, Systems and Computers, pp. 54–59, November 2003. View at: Google Scholar
 I. Walterscheid, T. Espeter, J. Klare, A. R. Brenner, and J. H. G. Ender, “Potential and limitations of forwardlooking bistatic SAR,” in Proceedings of the 30th IEEE International Geoscience and Remote Sensing Symposium (IGARSS '10), pp. 216–219, July 2010. View at: Publisher Site  Google Scholar
 J. Wu, J. Yang, Y. Huang, H. Yang, and H. Wang, “Bistatic forwardlooking SAR: theory and challenges,” in Proceedings of the IEEE Radar Conference (RADAR '09), pp. 1–4, Pasadena, Calif, USA, May 2009. View at: Publisher Site  Google Scholar
 Z. Shi, H. Wang, P. Zhang, X. Tang, B. Wu, and C. Wang, “Study on the probability of successful handoff of missile trajectory from midcourse guidance to terminal guidance,” in Proceedings of the International Conference on Computational and Information Sciences (ICCIS '10), pp. 1005–1008, December 2010. View at: Publisher Site  Google Scholar
 Q. Zhao, X. Zhang, W. Song, and F. Lu, “Rotational analysis of certain terminalguidance projectile,” in Proceedings of the 8th International Conference on Electronic Measurement and Instruments (ICEMI '07), pp. 4415–4418, August 2007. View at: Publisher Site  Google Scholar
 C. Wu, “A digital system to produce imagery from SAR data,” in Proceedings of AIAA Conference on Systems Design Driven by Sensors, October 1976. View at: Google Scholar
 J. R. Bennett and I. G. Cumming, “A digital processor for the production of SEASAT synthetic aperture radar imagery,” in Proceedings of Symposium on Machine Processing of Remotely Sensed Data, June 1979. View at: Google Scholar
 C. Wu, K. Y. Liu, and M. Y. Jin, “Modeling and a correlation algorithm for spaceborne SAR signals,” IEEE Transactions on Aerospace and Electronic Systems, vol. 18, no. 5, pp. 563–575, 1982. View at: Google Scholar
 R. K. Raney, H. Runge, R. Bamler, I. G. Cumming, and F. H. Wong, “Precision SAR processing using chirp scaling,” IEEE Transactions on Geoscience and Remote Sensing, vol. 32, no. 4, pp. 786–799, 1994. View at: Publisher Site  Google Scholar
 M. Soumekh, “Bistatic synthetic aperture radar inversion with application in dynamic object imaging,” IEEE Transactions on Signal Processing, vol. 39, no. 9, pp. 2044–2055, 1991. View at: Publisher Site  Google Scholar
 M. Soumekh, “Widebandwidth continuouswave monostatic/bistatic synthetic aperture radar imaging,” in Proceedings of IEEE International Conference on Image Processing (ICIP '98), vol. 3, pp. 361–365, IEEE, Chicago, Ill, USA, October 1998. View at: Publisher Site  Google Scholar
 M. Soumekh, J. Choi, and W. Zwolinski, “Target imaging using a synthesized bistatic radar,” in Proceedings of the IEEE International Conference on Systems Engineering, pp. 39–42, August 1989. View at: Google Scholar
 Y. L. Neo, F. H. Wong, and I. G. Cumming, “Focusing bistatic SAR images using nonlinear chirp scaling,” in Proceedings of International Conference on Radar Systems, 2004. View at: Google Scholar
 G. W. Davidson, I. G. Cumming, and M. R. Ito, “A chirp scaling approach for processing squint mode SAR data,” IEEE Transactions on Aerospace and Electronic Systems, vol. 32, no. 1, pp. 121–133, 1996. View at: Publisher Site  Google Scholar
 D. D' Aria, A. M. Guarnieri, and F. Rocca, “Focusing bistatic synthetic aperture radar using dip move out,” IEEE Transactions on Geoscience and Remote Sensing, vol. 42, no. 7, pp. 1362–1376, 2004. View at: Publisher Site  Google Scholar
 A. M. Guarnieri and F. Rocca, “Reduction to monostatic focusing of bistatic or motion uncompensated SAR surveys,” IEE Proceedings: Radar, Sonar and Navigation, vol. 153, no. 3, pp. 199–207, 2006. View at: Publisher Site  Google Scholar
 Y. L. Neo, F. Wong, and I. G. Cumming, “A twodimensional spectrum for bistatic SAR processing using series reversion,” IEEE Geoscience and Remote Sensing Letters, vol. 4, no. 1, pp. 93–96, 2007. View at: Publisher Site  Google Scholar
 Y. L. Neo, F. H. Wong, and I. G. Cumming, “Processing of azimuthinvariant bistatic SAR data using the range doppler algorithm,” IEEE Transactions on Geoscience and Remote Sensing, vol. 46, no. 1, pp. 14–21, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ziqiang Meng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.