New Approach for Unambiguous High-Resolution Wide-Swath SAR Imaging

The high-resolution wide-swath (HRWS) SAR system uses a small antenna for transmitting waveform and multiple antennas both in elevation and azimuth for receiving echoes. It has the potential to achieve wide spatial coverage and fine azimuth resolution, while it suffers from elevation pattern loss caused by the presence of topographic height and impaired azimuth resolution caused by nonuniform sampling. A new approach for HRWS SAR imaging based on compressed sensing (CS) is introduced. The data after range compression of multiple elevation apertures are used to estimate direction of arrival (DOA) of targets via CS, and the adaptive digital beamforming in elevation is achieved accordingly, which avoids the pattern loss of scan-on-receive (SCORE) algorithm when topographic height exists. The effective phase centers of the system are nonuniformly distributed when displaced phase center antenna (DPCA) technology is adopted, which causes Doppler ambiguities under traditional SAR imaging algorithms. Azimuth reconstruction based on CS can resolve this problem via precisely modeling the nonuniform sampling. Validation with simulations and experiment in an anechoic chamber are presented.


Introduction
As an active microwave remote sensing imaging tool, synthetic aperture radar (SAR) has the unique capability of obtaining abundant electromagnetic information from ground objects throughout the day and night and in all weather conditions.As such, it has been widely used in military reconnaissance, mapping, resource exploration, environmental protection, disaster management, and so forth.Two most important characteristics of modern SAR systemsazimuth resolution and swath coverage-are contradictive in system design [1].To coordinate this contradiction, highresolution wide-swath (HRWS) SAR is proposed [2][3][4].It uses separate small antenna transmitting radar signals, multiple elevation, and azimuth subantennas receiving echoes.Unlike traditional monostatic SAR, HRWS SAR adopts adaptive beamforming technology and displaced phase center antenna (DPCA) technology in elevation and azimuth, respectively.
HRWS SAR system uses multiple elevation antennas and DBF to form a sharp and high gain pattern which follows the direction of echoes in real time.The weighting coefficients for each channel are set as a function of time.This method is called SCORE (scan-on-receive).The direction of arrival (DOA) is calculated based on the vertical slant-range plane acquisition geometry, under the hypothesis of a strict spherical earth model; that is, no topographic height is considered.While, in real acquisition scenarios, characterized by mountains and valleys, displacements exist between the actual DOA and the preset DOA, which will cause pattern loss with respect to the ideal operational conditions [5,6].HRWS SAR system adopts DPCA technique which uses multiple apertures in azimuth and acquires, for each transmitted pulse, additional samples along the synthetic aperture.As a result, the pulse repetition frequency (PRF) enables the unambiguous mapping of a wider image swath.To ensure that these azimuth samples are uniformly distributed, there is a rigid restriction posed on its selection of PRF: SAR platform 2 International Journal of Antennas and Propagation moves just one half of its total antenna length between subsequent radar pulses.Such a rigid selection of the PRF is often in conflict with the optimal design of other system parameters.If this rigid condition is not satisfied, there will be periodically nonuniform sampling in azimuth and azimuth ambiguities will appear under traditional imaging algorithms based on matched filter [7,8].
In this paper, a new approach for unambiguous HRWS SAR imaging is proposed based on compressed sensing (CS).Multiple elevation spatial sampling points are used to estimate the real-time DOA of echoes, so that the displacement caused by acquisition scenarios is eliminated.As the echoes in elevation are limited by the numbers of subantennas, DOA estimation can be achieved effectively based on CS methods.Real-time adaptive DBF is achieved according to the estimated DOA. CS is introduced into azimuth DPCA SAR imaging.Although nonuniform sampling exists in azimuth, the observation model can be precisely constructed according to system and platform parameters.DPCA imaging algorithm based on CS can well recover the target scene and effectively suppress the azimuth ambiguities.
The rest of the paper is organized as follows.In Section 2 overviews of CS are presented.In Section 3, the model of elevation DBF and azimuth DPCA reconstruction are described in detail.Algorithms are validated via simulations and experiments in Section 4. Finally, Section 5 sums up the paper and presents conclusions.

Compressed Sensing
The theory of compressed sensing (CS) mainly focuses on the perfect reconstruction of the signal which has certain degree of sparsity.Assume that signal x (x ∈ C  ) is  sparse under the basis Ψ × , which indicates that it can be represented by Ψ with only  ( ≪ ) nonzero coefficients where  is an ×1 column vector and there are only  nonzero elements in .The observation model of CS is expressed as where Φ is the measurement matrix which contains  ×  elements.It is an ill-posed problem to recover x when  < .However, according to the principle of CS, x can be exactly recovered by solving an ℓ 1 optimization problem if the sensing matrix Θ = ΦΨ satisfies the restricted isometry property (RIP) condition.The number of measurements (i.e., the number of elements of y) should be  =  for the assurance of perfect reconstruction, in which  = (log(/)) is the oversampling factor determined by the sparsity of target signal (the sparsity is defined by /).If signal x is nonideally sparse or the measurement y contains noise, the problem turns into [9,10]  where  describes the error level introduced by the noise or the approximation of x which is nonideally sparse.
If the sensing matrix Θ satisfies  order RIP, then, for each  sparse signal , there exists the restricted isometry constant (RIC)   ∈ (0, 1) as It is a combinational problem to determine whether a specified matrix satisfies certain order RIP and the RIC cannot be obtained easily.The noncoherence feature is a more widely used criterion which is numerically defined by the maximum cross-correlation coefficient among the columns of Θ where Θ() and Θ(V) represent the th and the Vth column of Θ, respectively, and ⟨ ⋅, ⋅ ⟩ denotes the inner-product of two vectors.In fact, (Θ) equals the maximum value of the nondiagonal elements in the matrix G = Θ  Θ after normalized processing.A larger value of (Θ) indicates the more similar two columns in the sensing matrix, which may cause errors when the reconstruction algorithm tries to locate the support of the signal.The target signal can be recovered exactly if  satisfies the following condition [9]: ) .
It has been proved that CS has the potential of reducing the radar raw data volume and simplifying the system complexity.The application of CS combined with imaging radar technology has become the research focus in recent years [11].

HRWS SAR Signal Model Based on CS
The antenna structure of HRWS system is shown in Figure 1.
There are  subantennas in elevation direction and  subantennas in azimuth direction (for simplicity, the values of  and  are set to be the odd numbers).A small antenna in the antenna system center transmits the signal and all the subantennas receive the echoes.The length of transmitting antenna in the elevation direction is designed to be small enough to cover the wide swath.

Elevation Signal Model.
To increase the antenna gain and reduce the ambiguity, we adopt the digital beam forming (DBF) for the multiple receiving channels to obtain a narrow beam and change the weighting factors of the antenna array instantly to make the beam scan from the near end to the far end within the swath to track and receive the echo of observation terrain.The geometric model of DBF in the elevation direction is shown in Figure 2.After the pulse compression for each echo signal from the multiple channels, the wave direction in each range cell is sparse.According to the CS theory, the direction of the wave can be exactly estimated with few data.In the following, we establish the CS-based elevation signal model.In order to simplify the description, the signal model in a single range cell is considered.The transmitted signal is expressed as where  is the range time,  is the pulse duration time,   is the frequency modulation rate, and   is the carrier frequency.Taking the  subantennas (which locate in the center of azimuth direction) for example, and assuming a zero slant look angle, the distance between the target point and the th ( = 1, . . ., ) subantenna in the elevation direction is where  (+1)/2 is the distance between the ( + 1)/2 subantenna and the target point.  is the distance of the adjacent phase centers in the elevation direction. is the angle between the antenna point direction and the vertical direction to earth center. is the angle between the wave direction from target point and its vertical direction to earth center.In case of far-field,   is written as where the third item can be ignored since it is relatively small compared with the former two items, For the point target in the nearest slant range, the two-way delay in time is where  0 is the light speed and  0 is the two-way delay in time for the ( + 1)/2th subantenna.As a conclusion, regardless of the electromagnetic loss in the air, the system loss, and the envelope change caused by the platform position, the baseband form of the received signal from the th subantenna is where  is the backscattering coefficients of the point target.
After demodulation and range compression, the signal is written as where  is the carrier wavelength.By registration processing, Thus, for the samples at  =  0 , the compressed signal is expressed as In the aforementioned formula, the exponential item exp (2   0 ) is included into .Considering the additive noise, the vector form of the multichannel data is expressed as where y, a(), and e r are  × 1 vectors.e r is the noise vector.a() is the point vector which includes the phase information among different channels and has a direct connection with the wave direction.The th element of a() is written as The system may receive echoes from many directions at a single time, which is caused by the range ambiguity International Journal of Antennas and Propagation and folding effect.In order to obtain the correct wave direction, the wave direction range is discretely divided into several angle grids where  0 is the wave direction of the near range of the swath and Δ is the minimum discrete angle step.The expression which contains echo from all the wave directions can be written as In matrix-vector form, where A is an  ×  matrix, The number of the wave directions is finite and sparse compared with all the directions in real application, which means that the number of nonzero elements in  is far less than .By applying CS, they can be estimated and then the wave direction   can be calculated accordingly.The estimation of  is modeled as min (|α|) , st.

Azimuth Signal
Model.The HRWS SAR system adopts displaced phase center antenna (DPCA) technology to retrieve the azimuth signal by rearranging the multichannel echoes in the order of equivalent phase center and avoids the Doppler ambiguity caused by limited sampling.However, DPCA requires that SAR platform moves just one half of its total antenna length between subsequent radar pulses.This condition is often difficult to meet in real applications, which makes the arrangement of the equivalent phase center nonuniform.And the azimuth ambiguity emerges when directly processing with the nonuniform azimuth signal.Currently, there are two methods to solve the problem, the phase compensation method and the spectrum reconstruction method.The phase compensation method corrects the phase of the multiaperture data so that it corresponds to single uniform SAR signal.The main disadvantage of this method is the severe impact on the resolution [12].The spectrum reconstruction method retrieves the unambiguous signal from aliasing Doppler spectrum by solving a serial of filter banks.It can be done in time domain or frequency domain, but the calculation complexity of this algorithm is time-consuming [13].
The CS-based azimuth nonuniform sampling reconstruction method builds the exact observation model by using the system and platform parameters and then achieves the recovery through solving an ℓ 1 norm minimizing problem.Taking the  azimuth subantennas locating in the elevation center of the antenna array, for example, with the ignorance of the influence on the signal amplitude caused by range compression, the echo model of the th ( = 1, . . ., ) subantenna is where  0 is the nearest slant range, V is the platform velocity,  is the azimuth time, and   is the length of the receiving subantenna.Let Δ 0 = 1/; then the distance that the platform moves within a pulse repetition interval is  0 = Δ 0 V.For the system that has one transmitting channel and  receiving channels, the discrete grid of the observation scene can be set to   =  0 / to achieve the high-resolution reconstruction.As a conclusion, assume Δ = Δ 0 /; then discrete form of the azimuth samples is The azimuth observation model of the th subantenna is expressed with (26), in which the time interval of the row vector is Δ, while the time interval of the column vector is Δ.Consider Accordingly, the observation matrix Φ of the whole azimuth observation model should be reorganized from observation matrices of subapertures.After reorganization, the (, ) elements of Φ are where ⌊/⌋ is to get the largest integer less than / and After azimuth scene discretization, the reflectivities of scene can be expressed as an  × 1 column vector If we rearrange the data of each subaperture according to the subsequence of effective phase centers, the measurements can be expressed as follows: Thus, the observation model can be expressed as where e a is the corresponding additive observation noise vector.Finally, the observed target scene can be reconstructed through solving a constrained ℓ 1 minimization problem min (|x|) , st. ‖w − Φx‖ <     e a    2 . (32)

HRWS SAR Imaging Flowchart.
In order to reduce the amount of data satellite transmission, real-time estimation of DOA and digital beamforming are carried out with range compressed echoes from different channels.Only 1/ of data needs to be transmitted to the ground.Azimuth unambiguous reconstruction is achieved on the ground.So, HRWS SAR imaging processing is divided into two parts: onboard realtime processing and ground processing, as shown in Figure 3.

HRWS SAR Simulation
4.1.Elevation DOA and Adaptive Beamforming.According to the simulation parameters in Table 1, DOA of the echoes from targets at different heights can be calculated.Angular displacements in the presence of topographic height at different slants are illustrated in Figure 4.It can be found that angular displacement reaches 0.1 ∘ to 0.2 ∘ at height of 1 km, and it reaches 0.3 ∘ to 0.4 ∘ at height of 2 km.Angular displacement increases with increasing height of target.SCORE pattern loss in the presence of topographic height can also be calculated as shown in Figure 5.We find that when the target height reaches 2 km which is a common situation, the pattern loss is close to 2 dB.Therefore, it is necessary to use real-time digital beam forming techniques to estimate the DOA and avoid the pattern loss due to target height.We consider three targets with the same slant range distance but with different incidence angles, which are 30 ∘ ,  32 ∘ , and 33 ∘ , respectively, and the incidence angle space [29.6, 34.9] ∘ according to the whole swath is divided into 530 subangels evenly with each angle interval being 0.01 ∘ .The signals are pure and without noise pollution.Laplace optimization algorithm is adopted for CS reconstruction in our simulation [14].Results of DOA estimation based on CS are shown in Figure 6.We find that the three targets with different DOA are correctly recovered.
In real acquisition scenarios, we need to estimate DOA at each range sampling time and the number of used samples is decided by the number of subantennas in elevation direction influence the estimation result.The performance of DOA estimation is calculated empirically via simulation.We randomly set three targets with different topographic height and incident angle for each trail.The trail is declared a success if the sum of the differences between estimated DOA and preset DOA is within the range of 0.1 ∘ .The empirical probability of success is determined by repeating the CS reconstruction trials 1000 times and calculating the percentage of success.This percentage of correct DOA estimation is plotted as a function of SNR in Figure 7 (labeled as monoazimuth observation).The results show that the method with monoazimuth sampling can by 100% correctly estimate DOA only when SNR is higher than 15 dB.This level of SNR for range compressed signal is hard to achieve.If more independent observations are used, the DOA estimation performance will be much improved.As the satellite is moving, observation geometry for the same target at adjacent azimuth time changes.And the satellite echo delay varies because of the existence of range cell migration (RCM), which affects the correct estimation of DOA.In fact, if we adopt nearby azimuth samplings to estimate DOA jointly, it is possible that the error is in the tolerable range.Meanwhile, if target position and amplitude variations caused by RCM exceed the tolerance range, they can be compensated by registration.These variations caused by platform motion can be illustrated by the following simulation in which 100 adjacent azimuth sampling points are adapted.DOA change and RCM against number of azimuth samples are shown in Figures 8 and 9, respectively.We can find that DOA   changes and RCM is very small and can be ignored in actual processing.
The performance of DOA estimation from multiple azimuth samples is calculated in the same way as monoazimuth sampling.The number of adjacent azimuth samplings is set as 16, so the number of effective samples for each estimation reach 16 ×  = 300.From results illustrated in Figure 7, we find that joint use of multiple azimuth samples can by 100% correctly estimate DOA even when SNR is higher than −3 dB.range Doppler (RD) algorithm [15] is shown in Figure 10(a), where azimuth ambiguity reaches to −19 dB.The CS reconstruction results in azimuth are shown in Figure 10(b).As seen in Figure 10, azimuth ambiguity caused by nonuniform sampling is effectively suppressed by CS methods to a level below −38 dB.

Experiment in an Anechoic Chamber
Experiments for azimuth ambiguity suppression using CS is also carried out in an anechoic chamber.Antenna configuration and trihedral corner reflector of this ground-based DPCA SAR system are shown in Figure 11.It has one aperture transmitting stepped frequency and three apertures receiving echoes in azimuth.These four apertures are 3.2 cm long and tightly arranged.The system parameters are listed in Table 2.
International Journal of Antennas and Propagation    Effective azimuth centers of this DPCA SAR system are shown in Figure 12 and they are nonuniformly distributed.The observed trihedral corner reflector located at the slant range of 2.5 meters.The experimental results with traditional RDA and CS algorithm are shown in Figure 13.We use impulse response width (IRW), peak sidelobe ratio (PSLR), and integrated sidelobe ratio (ISLR) as the performance indexes, and the results are shown in Table 3.It becomes clear that DPCA SAR imaging result with traditional RDA is strongly blurred because the azimuth effective phase centers are nonuniformly spaced.Our proposed CS imaging algorithm can effectively suppress the azimuth ambiguity caused by nonuniform sampling while keeping high azimuth resolution.

Conclusion
In this paper, DOA estimation from multiple-elevationchannel echoes and adaptive beamforming in real time are proposed, which can effectively solve elevation pattern loss problems in HRWS SAR system.In low SNR, the joint use of multiple azimuth observation data can ensure the correct estimation of DOA.HRWS SAR's one aperture transmitting and several apertures receiving mode caused nonuniform sampling in azimuth and further brought azimuth ambiguity.
Our proposed method based on CS can effectively suppress this ambiguity when precise observation models are established according to systems and platforms parameters by solving an optimization problem.

Figure 4 :Figure 5 :
Figure 4: Angular displacement in the presence of topographic height.

Figure 6 :Figure 7 :
Figure 6: Estimation results of DOA based on CS.

Figure 8 :
Figure 8: DOA variation versus number of azimuth samples.

Figure 10 :
Figure 10: Azimuth slice of imaging results when nonuniform sampling exits in azimuth.

Figure 11 :
Figure 11: Antenna configuration and trihedral corner reflector in ground-based DPCA SAR system.

Figure 12 :
Figure 12: Azimuth effective phase centers of ground-based SAR system.

Figure 13 :
Figure 13: Experimental results of azimuth profile.

Table 3 :
Comparison of performance for imaging results.