Research Article  Open Access
Joint TwoDimensional Ambiguity Resolving Based on SpaceTime Filtering for MIMOSAR
Abstract
In order to resolve the range and azimuth ambiguity in MIMOSAR imaging, this paper proposed a joint twodimensional ambiguity resolving method based on spacetime filtering, which according to the ambiguity in different subscene is relative to different nearest echo range, and different azimuth time is relative to different instantaneous incidence angle, to make twodimensional spacetime steering vector for resolving ambiguity of the echo, which also can finish the imaging processing at the same time. Simulation results validate the proposed processing approach for the MIMOSAR system.
1. Introduction
Multipleinput and multipleoutput spaceborne synthetic aperture radar (MIMOSAR), as a high resolution and wide swath microwave remote sense system [1–5], is proved to be an extremely useful surveillance tool for moving target indication, wide swath imaging, and threedimensional imaging.
ScanSAR also can obtain wide unambiguous swath coverage at the cost of a degraded azimuth resolution [6]. Multichannel line arrays can be arranged by range or azimuth [7], for the range multichannel can obtain wide swath coverage, utilizing digital beamforming to resolve range ambiguity, but there is still blind zone [8, 9]. With the emergence of techniques such as displaced phase center [10] and multiple aperture reconstruction [11, 12], a number of displaced azimuth subapertures are used to receive echoes in one pulse repetition interval (PRI), thereby increasing the effective azimuth sampling rate without changing the pulse repetition frequency (PRF). In this way, one can lower the requirement on PRF to achieve wider swath imaging without impairing azimuth resolution; in other words, one can broaden the Doppler spectrum to obtain finer azimuth resolution without reducing the swath coverage.
To obtain wide swath coverage and high resolution, one can utilize twodimensional array to transmit and receive SAR signals and then resolve the range and azimuth ambiguity [13]. However, this method makes range point target downlooking angle as range steer vector in approximate treatment and does not consider the changing of the downlooking angle with the azimuth time.
In this paper, a joint twodimensional ambiguity resolving method for MIMOSAR is proposed, which uses deferent targets relative to deferent range as steer vector which is independent of the downlooking angle. Firstly, the azimuth ambiguity is resolved. Then, the range ambiguity resolving and azimuth match filtering are performed instantaneously to get target image, which makes the imaging algorithms simpler.
This paper is organized as follows. Section 2 reviews the signal model of MIMOSAR. The twodimensional ambiguity resolving algorithm is analyzed in Section 3. Section 4 gives the simulations and validates the proposed processing approach for the MIMOSAR system. Finally, Section 5 concludes this paper.
2. MIMOSAR System Mode and Signal Model
Figure 1 presents the geometry of a multichannel SAR. It supposed that the whole antenna aperture was divided into subapertures, where is the number of azimuth subapertures and its length is . is the number of azimuth subapertures and its height is . We denoted the th column and the th row as , , with the azimuth beginning coordinate , and , , with the range beginning coordinate.
Use the subcaptures on the same column range dimension transmit signals at the same time; through adjusting phase weighting to control the wave beam, the center of transmit wave beam is , , and the whole antenna receives the echo signal. According to the principle of equivalent phase center, this mode can be known as each equivalent phase center transmitting and receiving itself. The equivalent phase centers located at , where and . The deferent range subapertures light the deferent subswaths from far to nearby adjusting phase weightings. The ellipses in Figure 1 are the subswathes.
Assume that, during the th pulse period, the first range dimension subapertures transmit signal in the same time and make the wave beam light the th subswath through controlling phase weighting. The transmitting signal can be expressed as follows: where is the downlooking angle of the center line of th subswath, is the carrier frequency, is the wavelength, is the delay time of the th subpulse, is the whole time, is the fast time, is a integer, is pulse repetition interval, and is the chirp rate.
The baseband signal received by the pth row and the qth column subaperture was where is slow time, is a integer, and the value scope of is from 1 to , corresponding to the swath from far to near. is the flat velocity, and is the length of synthetic aperture. is the instantaneous slant range, is the coordinates of one scatter point in the th scene, and . Assume that scene is plane, .
Performing the range match filtering to (2), we can gain where is the instantaneous distance from the scatter point on the th subswath to .
As we know, the MIMOSAR can be considered as a space signal sampling which should satisfy the sampling equation as follows: where denotes the position change of antenna element and denotes the wavelength of the signal.
So the time delay should be where denotes the carrier frequency.
And the time resolution can be calculated by where denotes the bandwidth.
In a narrow band system is much larger than , so we can get much smaller than , which means that the time delay made by the position change of antenna element is very small compared with the time resolution.
The time delay made by the position change of antenna element can be ignored, which is reasonable because the distance between the elements is very small.
This system uses low PRF to get wide swath, but the low PRF can cause azimuth Doppler ambiguity; Figure 2 illuminates this.
(a) The Doppler is not ambiguity
(b) The Doppler is ambiguity
After azimuth fast Fourier transforms (FFT), we can gain the Doppler ambiguity signal : where , is the vertical distance form scatter point to flight track, is the instantaneous distance of scatter point, is the angle of instantaneous slant distance and flight track, and , .
When is satisfied, different scatter points on subswathes can make aliasing and then cause range ambiguity. But the different scatter points are corresponding to different downlooking angles, as shown in Figure 3.
(a) The range is not ambiguity
(b) The range is ambiguity
According to the forgoing analysis, if the scatter instantaneous slant distance satisfies the preceding equation, (7) can be shown as where and are caused by the position of azimuth and range subapertures, and they are the steer vector for resolving ambiguity. One can know that there is the azimuth frequency in the range steer vector, so one needs to resolve the azimuth ambiguity firstly and then resolve the range ambiguity.
3. TwoDimensional Ambiguity Resolving
3.1. Resolving Azimuth Ambiguity
Because (8) is range and Doppler ambiguity signal, we express the matrix vector as
In this paper, we use joint 2dimensional method to resolve ambiguity based on static weighting vector. It is assumed that the vector matrix is and the th column weighting vector is shown as
It means that the corresponding downlooking angle and azimuth angle position is 1 and other ambiguities are zero. Consider where expresses matrix conjugate transpose operator and , , when , . Consider where denotes the pseudoinverse operation, . Consider
Using the weighting vector to resolve ambiguity, the following is gained: After this, the resolving azimuth ambiguity operation is finished.
According to array signal theory, the Doppler ambiguity time must satisfy the following formula: . To the member of azimuth subaperture , it can make one restriction direction at most and the other is zero, so only time’s ambiguity can be resolved on azimuth Doppler dimension. With the same principle, only time ambiguity can be resolved on azimuth Doppler dimension.
After resolving the azimuth ambiguity, the signal is combined according to the Doppler spectrum. The range resolving ambiguity is followed.
3.2. Resolving Range Ambiguity
According to (8), the range steering vector and the azimuth matching function are concerned. We expressed matrix vector as follows: It is assumed that the vector matrix is and the th column weighting vector is shown as: It means that the corresponding downlooking angle and azimuth angle position is 1 and other ambiguities are zero. Consider where expresses matrix conjugate transpose operator and , , when . Consider where denotes the pseudoinverse operation, . Consider
Utilizing the weighting vector to resolve the ambiguity, we can gain At this time, the range ambiguity was resolved completely. The azimuth was finished when range was resolved. Then after the azimuth inverse FFT (IFFT), we can finish the azimuth pulse compression.
For different point target slant distance , the subpulse time delay should be added. For the th subswath, added the time delay to the range envelope, which make the real position of scatter point was get. Then the whole wide swath signal was obtained.
3.3. Impact of Terrain Height Variance
Forenamed static weighting vector is based on the assumption that the terrain is flat. Here we can analyze the impact of terrain height variance on the weighting vector.
As shown in Figure 4, there is a point with the downlooking angle and azimuth angle ; we assume there is another point with a height difference between , so the phase difference of the weighting vector pointing to and is where , are the downlooking angle and azimuth angle in the position of . From the geometry configuration of Figure 4, we can get that and , so the phase difference can be written as
The clutter which remained caused by the phase difference is
With the simulating parameter of , and , we get that the simulation result of the clutter which remained varied with the height difference shown in Table 1.

From Table 1 we can see that if the change of terrain height is below 1 km, our proposed method can still be effective.
4. Simulation Experiment
To validate the proposed approach, a simulation experiment on point targets is carried out. The main system parameters are listed in Table 2. In this simulation, the whole antenna aperture was divided into subapertures.

It is supposed that these scatter points have the same backscattering coefficient. Through calculation one can get , so the three scatter points are rangeambiguous with each other. The Doppler band width is and the azimuth sample frequency is , which makes the azimuth Doppler spectrum ambiguous with three times. Figure 5(a) shows the image of an original single channel. The azimuth and the range both have ambiguity. Range ambiguity makes three scatter points with different range positions mixed together, so only one point can be seen in Figure 5(a); because the system Doppler ambiguity three times, three focused points can be seen along azimuth. Both sides were ambiguity points and their energy was smaller than the center point obviously. Through the joint twodimensional ambiguity resolving method proposed in this paper, the unambiguity image can be obtained as shown in Figure 5(b). After twodimensional ambiguity resolving, the scatter points were located in their real positions. Figure 5(c) is azimuth profile of the far scatter point. In the figure, the standard sinc function was obtained after the azimuth pulse compressed, which shows the unambiguous Doppler spectrum was recovered by the joint twodimensional processing. Figure 5(d) is range profile of the three range points. From the figure, one can note that the range ambiguity is also resolved.
(a) The image with 2D ambiguity
(b) The image after twodimensional ambiguity resolving
(c) Azimuth profile of far scatter point
(d) Range profile of middle azimuth cell
5. Conclusion
In this paper, a novel joint twodimensional ambiguity resolving method based on spacetime filtering was proposed, which can resolve the range and azimuth ambiguity of MIMOSAR for highresolution and wide swath imaging. This method does not use any approximate operation, which can finish resolving ambiguity and two dimensions focused instantaneously. Imaging results of the simulated point target validate the proposed approach. Because this method does not depend on SAR system work mode, it has certain general use.
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 Grant 61201433 and supported by Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT14B09).
References
 G. Krieger, N. Gebert, M. Younis, and A. Moreira, “Advanced synthetic aperture radar based on digital beamforming and waveform diversity,” in Proceedings of the IEEE Radar Conference (RADAR '08), pp. 767–772, Rome, Italy, May 2008. View at: Publisher Site  Google Scholar
 A. Ossowska, J. H. Kim, and W. Wiesbeck, “Modeling of nonidealities in receiver frontend for a simulation of multistatic SAR system,” in Proceedings of the 4th European Radar Conference (EURAD '07), pp. 13–16, Munich, Germany, October 2007. View at: Publisher Site  Google Scholar
 G. Krieger, N. Gebert, and A. Moreira, “Multidimensional waveform encoding: a new digital beamforming technique for synthetic aperture radar remote sensing,” IEEE Transactions on Geoscience and Remote Sensing, vol. 46, no. 1, pp. 31–46, 2008. View at: Publisher Site  Google Scholar
 W. Q. Wang, “Spacetime coding MIMOOFDM SAR for highresolution imaging,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 8, pp. 3094–3104, 2011. View at: Publisher Site  Google Scholar
 W. Wang, “Mitigating range ambiguities in high PRF SAR with OFDM waveform diversity,” IEEE Geoscience and Remote Sensing Letters, vol. 10, no. 1, pp. 101–105, 2013. View at: Publisher Site  Google Scholar
 P. Huang and W. Xu, “A new spaceborne burst synthetic aperture radar imaging mode for wide swath coverage,” Remote Sensing, vol. 6, no. 1, pp. 801–814, 2014. View at: Google Scholar
 P. Huang, S. Li, and W. Xu, “Investigation on fullaperture on multichannel azimuth data processing in TOPS,” IEEE Geoscience and Remote Sensing Letters, vol. 4, no. 11, pp. 728–732, 2014. View at: Google Scholar
 G. D. Callaghan and I. D. Longstaff, “Wideswath spaceborne SAR and range ambiguity,” in Proceedings of the IEEE Radar Conference, pp. 248–252, Edinburgh, Scotland, October 1997. View at: Publisher Site  Google Scholar
 G. D. Callaghan and I. D. Longstaff, “Wideswath spaceborne SAR using a quadelement array,” IEE Proceedings: Radar, Sonar and Navigation, vol. 146, no. 3, pp. 159–165, 1999. View at: Publisher Site  Google Scholar
 G. Krieger, N. Gebert, and A. Moreira, “Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling,” IEEE Geoscience and Remote Sensing Letters, vol. 1, no. 4, pp. 260–264, 2004. View at: Publisher Site  Google Scholar
 G. Krieger, N. Gebert, and A. Moreira, “SAR signal reconstruction from nonuniform displaced phase centre sampling,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium Proceedings: Science for Society: Exploring and Managing a Changing Planet (IGARSS '04), vol. 3, pp. 1763–1766, Anchorage, AK, USA, September 2004. View at: Publisher Site  Google Scholar
 N. Gebert, G. Krieger, and A. Moreira, “SAR signal reconstruction from nonuniform displaced phase centre sampling in the presence of perturbations,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '05), pp. 1034–1037, Seoul, Korea, July 2005. View at: Publisher Site  Google Scholar
 P. Huang, W. Xu, and W. Qi, “Two dimension digital beamforming pre processing in multibeam Scansar,” Progress in Electromagnetics Research, vol. 136, pp. 495–508, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Pingping Huang 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.