Joint Two-Dimensional Ambiguity Resolving Based on Space-Time Filtering for MIMO-SAR

In order to resolve the range and azimuth ambiguity in MIMO-SAR imaging, this paper proposed a joint two-dimensional ambiguity resolving method based on space-time 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 space-time 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 MIMO-SAR system.


Introduction
Multiple-input and multiple-output spaceborne synthetic aperture radar (MIMO-SAR), as a high resolution and wide swath microwave remote sense system [1][2][3][4][5], is proved to be an extremely useful surveillance tool for moving target indication, wide swath imaging, and three-dimensional 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 two-dimensional 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 two-dimensional ambiguity resolving method for MIMO-SAR 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 MIMO-SAR.The two-dimensional ambiguity resolving algorithm is analyzed in Section 3. Section 4 gives the simulations and validates the proposed processing approach for the MIMO-SAR system.Finally, Section 5 concludes this paper. ×  subapertures, where  is the number of azimuth subapertures and its length is  1 . is the number of azimuth subapertures and its height is  2 .We denoted the pth column and the qth row as   =  0 + ( − 1) 1 ,  = 1, . . ., , with  0 the azimuth beginning coordinate (  ,   , ), and   =  0 + ( − 1) 2 ,  = 1, . . ., , with  0 the range beginning coordinate.

MIMO-SAR System Mode and Signal Model
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    = (  +   0 )/2,   0 =  0 +  2 ( − 1)/2, 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    = ( 0 +   )/2 and    = (  +   0 )/2.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 Lth pulse period, the first range dimension subapertures transmit signal in the same time and make the wave beam light the Lth subswath through controlling phase weighting.The transmitting signal can be expressed as follows: where   is the downlooking angle of the center line of lth subswath,   is the carrier frequency,  = /  is the wavelength, Δ  =   1 + ⋅ ⋅ ⋅ +   −1 is the delay time of the lth 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 (  ) = sin( 2 (sin   − sin )/2)/ sin( 2 (sin   − sin )/2),   =   is slow time,  is a integer, and the value scope of  is from 1 to , corresponding to the swath from far to near. = V  , V is the flat velocity, and   is the length of synthetic aperture.
is the coordinates of one scatter point in the Lth scene, and  =  0 −   .Assume that scene is plane,   = 0.
Performing the range match filtering to (2), we can gain where is the instantaneous distance from the scatter point on the lth subswath to ( 0 ,  0 ,  0 ).As we know, the MIMO-SAR 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.
After azimuth fast Fourier transforms (FFT), we can gain the Doppler ambiguity signal   (  ,   ): where 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.
According to the forgoing analysis, if the scatter instantaneous slant distance satisfies the preceding equation, (7) can be shown as where exp(−(2  /V)√ 2  −  2  ) and exp((4   cos   (  )/)) 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.

Two-Dimensional Ambiguity Resolving
3.1.Resolving Azimuth Ambiguity.Because ( 8) is range and Doppler ambiguity signal, we express the matrix vector z  (  ) as In this paper, we use joint 2-dimensional method to resolve ambiguity based on static weighting vector.It is assumed that the vector matrix is W × 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 where () + denotes the pseudoinverse operation,  +  = pinv(  ).Consider S = [ 11 , . . .,  1 ,  12 , . . .,  2 , . . .,  1 , . . .,   ] 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  − 1 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.

Resolving Range Ambiguity.
According to (8), the range steering vector and the azimuth matching function are concerned.We expressed matrix vector z  (  ) as follows: It is assumed that the vector matrix is W × and the th column weighting vector is shown as: where () + denotes the pseudoinverse operation,  +  = pinv(  ).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 Δ  =   1 + ⋅ ⋅ ⋅ +   −1 to the range envelope, which make the real position of scatter point (  ,   ,   ) was get.Then the whole wide swath signal was obtained.

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  1 and azimuth angle  1 ; we assume there is another point   with a height difference ℎ between , so the phase difference of the weighting vector pointing to  and   is where  2 ,  2 are the downlooking angle and azimuth angle in the position of   .From the geometry configuration of Figure 4, we can get that cos  1 = sin  1 sin  and cos  2 = sin  2 sin , so the phase difference can be written as The clutter which remained caused by the phase difference is With the simulating parameter of  = 700 km,  = 880 km,  = 0.03 m,  1 = 1 m, and  2 = 0.8 m, 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.

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 3 × 3 subapertures.It is supposed that these scatter points have the same backscattering coefficient.Through calculation one can get  1 =  Δ /2 +  2 =  Δ +  3 , so the three scatter points are range-ambiguous with each other.The Doppler band width is   = 2V/ = 3600 Hz and the azimuth sample frequency is 1200 Hz, 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 two-dimensional 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 two-dimensional 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.

Conclusion
In this paper, a novel joint two-dimensional ambiguity resolving method based on space-time filtering was proposed, which can resolve the range and azimuth ambiguity of MIMO-SAR for high-resolution 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.

Figure 1 2 InternationalFigure 1 :
Figure 1 presents the geometry of a multichannel SAR.It supposed that the whole antenna aperture was divided into

2 −PRF/ 2
The Doppler is not ambiguity PRF/Band width

Figure 2 :
Figure 2: The relationship between cosine of cone angle and Doppler.

Figure 3 :
Figure 3: The relationship between downlooking angle and instantaneous slant distance.

Figure 4 :
Figure 4: The geometry of terrain height variance.

Figure 5 :
Figure 5: The image of joint two-dimensional processing.

Table 1 :
Phase difference in the terrain with different height.