A Two Dimensional Overlapped Subaperture Polar Format Algorithm Based on Stepped-chirp Signal

In this work, a 2-D subaperture polar format algorithm (PFA) based on stepped-chirp signal is proposed. Instead of traditional pulse synthesis preprocessing, the presented method integrates the pulse synthesis process into the range subaperture processing. Meanwhile, due to the multi-resolution property of subaperture processing, this algorithm is able to compensate the space-variant phase error caused by the radar motion during the period of a pulse cluster. Point target simulation has validated the presented algorithm.


Introduction
Synthetic aperture radar (SAR) becomes an important tool in modern remote sensing for its allweather, day and night capability to provide high-resolution maps of scene of interest. The demand for radar images is constantly pushing for finer resolutions. This quest for the resolving power has two major consequences [1][2]: first, their useful bandwidth should increase in proportion with the resolution in range. Second, the length of the synthetic antenna should increase in proportion with the along-track resolution.
Due to technical limitations, particularly the limited sampling rate of the analog to digital converters, synthetic bandwidth technique [3][4][5][6] is proposed to solve the hardware challenges of wideband radar. As compared to the commonly used wideband linear frequency modulated (LFM) radar waveform, by adopting stepped-chirp waveform and applying synthetic bandwidth techniques, it is possible to achieve high range resolution while still retaining the advantages of lower instantaneous receiver bandwidth and lower analog-to-digital sampling rate. However, the relative motion between OPEN ACCESS the radar and the target during the period of a pulse cluster is not always negligible. If it is not taken into account, as a consequence, it may result in performance degradations, such as range error, loss in signal-to-noise ratio, and degraded range resolution. Unfortunately, the phase errors resulted from this motion are space-variant. The traditional synthetic bandwidth techniques only compensate these errors corresponding to a specified point target, for example, the scene center. Consequently, the purpose of this paper is to compensate this space-variant phase error.
With respect to the azimuth dimension, high resolution is obtained by coherent integration over a large aperture. The generally linear radar flight trajectory assumption, which is the basis of frequency domain image formation algorithm, is deviated, especially when nonplanar motion (NPM) occurs. Polar format algorithm (PFA) [7][8] is a popular high resolution spotlight SAR algorithm due to its efficient NPM compensation capability. However, due to the planar wavefront approximations made in PFA, the data exhibits space-variant phase errors and migration errors, which limit the focused scene size. Then, subaperture algorithm [9] is proposed. Due to its space-variant phase correction capability, fine resolution imaging of large scene becomes available.
In this paper, a new image formation algorithm which incorporates the synthetic bandwidth technique with subaperture processing is proposed. Instead of traditional pulse synthesis preprocessing, it integrates the pulse synthesis into range subaperture processing. Meanwhile, it is able to compensate the space-variant phase errors caused by the radar motion during the pulse cluster.

SAR Signal Modeling
Consider a spotlight SAR operating in the geometry of Figure 1, for simplicity, we only examine the case in which the sensor travels in a straight line at constant velocity. The central point of the scene is defined as the origin of the coordinate system, and radar antenna phase center (APC) is identified by coordinates ( , , ) x y z . The variable n α and n ϕ are the APC's instantaneous squint angle and grazing angle, respectively, at the nth LFM pulse. They equal 0 α and 0 ϕ at the center of the aperture. The distance from the APC to the scene center is cn r . In this paper, we account for the broadside mode only, then 0 0 α = . A target scatter is located at ( , , 0) x y s s , and the range from APC to this target is measured by sn r . After performing dechirp and residual video phase (RVP) elimination on the echo signal of the nth pulse, we can express the Doppler phase history as [7] ( ) where ( ) cn sn R r r ∆ = − is the differential range, c is the velocity of light, 0 f is the carrier frequency of the transmitted signal, γ is the chirp rate, s T is the sample interval in fast time, i is the range sample index with / 2 / 2 1 I i I − ≤ ≤ − , and n is the index value of pulse number with / 2 The differential range R ∆ can be expressed as [7] ( ) ( ) cos sin cos cos cos cos tan For the purpose of clearness and simplicity, we still use i instead of ' i in the following discussion. Then, equation (5) is expressed as If the space sampling position n satisfied tan n d n α α = , where dα is constant, the range resampled signal can be modeled as 3. Two Dimensional Overlapped ped ped ped Subaperture Polar format Algorithm (PFOSA) [9] Equation (7) is the phase history after range resampling. In full aperture PFA image formation, we get the image by performing an azimuth resampling followed by a 2-D DFT, or an azimuth chirp-z transform (CZT) followed by a range DFT. However, due to the space-variant phase error term ( ) n ξ , the focused scene size of interest is constrained to be very small in ultra-high resolution SAR. Subaperture algorithm, which can provide coarse resolution images before the final fine resolution image formation, has been proposed to overcome this constraint [9]. Due to the coarse information of the individual scatter's location extracted from the coarse resolution images, the compensation of space-variant phase error becomes applicable. In the following, we briefly review of the PFOSA proposed in [9].
First, we divide the azimuth and range aperture into subapertures, respectively, by making where 1 m is the azimuth intra-subaperture index limited within ∆ is the azimuth data decimation factor, 1 k is the range intra-subaperture index limited within µ is the range data decimation factor. Using equation (8), we rewrite the equation (7) as Next, applying the quadratic order approximation of ( ) n ξ and rearrange equation (9) where, the first term is a constant, which is neglected in the following discussion. The second and third terms are the azimuth and range intra-subaperture terms, which correspond to coarse resolution image. The fourth and fifth terms are the inter-subaperture terms, which correspond to fine resolution image. Also note that each exponential term contains some undesired error phase. In the subaperture terms, we choose subaperture size 1 M and 1 K in such a way that the phase error terms caused by wavefront curvature in the subaperture can be neglected. It is helpful to note that error terms in the last two terms can be compensated due to space position information extracted from the coarse resolution images. Now to facilitate the analysis, we rewrite equation (10) as following

Stepped-Chirp based PFOSA (SCPFOSA)
To reduce the transmission bandwidth, and meanwhile to achieve the high range resolution, it is possible to transmit series of narrow-band signals centered at different carrier frequencies. For example, an equivalent wideband LFM chirp can be assembled from lesser-bandwidth chirp segments in the data processing stage. These subchirp signals, which are referred to as a pulse cluster, are transmitted as separate pulses, each with their own carrier frequencies.
where 1 k is the range sample index in each chirp segment, 2 k is the chirp segment index, and n is the cluster index. It is important to note that the grazing angle 2 , n k ϕ and squint angle 2 , n k α vary not only with index n , which is desirable, but also with 2 k , which is undesirable. The latter change of 2 , n k ϕ and 2 , n k α is resulted from the radar motion during the pulse cluster. Neglecting this variation will introduce space-variant phase errors which limit the focused scene size. But in equation (12) it does not appear explicitly. To illustrate this effect, we develop a polynomial expression by using Taylor series expansion.
The key to analysis of the characteristics of the equation (12) is to have expression for  As before, we assume that the space sampling position n satisfies Phase term basic Φ contains the basic imaging information. The phase term err Φ , which causes distortion and defocus in range, is the error term introduced by the radar motion during pulse cluster.
Due to its dependence on target position ( ) Comparing with equation (6), it is clear to see that in equation (17) we get the range subaperture data naturally via the transmission and reception of the chirp segment signals. The difference is that is the phase error term resulted from radar motion during pulse cluster but after range resampling. Compared with phase error term err Φ before range resampling, parts of phase error are compensated, and the range position dependence is eliminated. Analogous to equation (11), for equation (18), we neglect the constant phase term and wavefront error terms in the subapertures and rearrange the rest terms , the two parts of ( ) 1 2 , err k k φ , are due to radar motion during the pulse cluster. These error terms are space-variant due to the dependence on azimuth position (range dependence is eliminated owing to range resampling). If these phase errors are not compensated, as a consequence, they result in displacement and defocus in range. Since the coarse location information can be extracted from the coarse resolution images, it is possible to compensate these errors by modifying the classical PFOSA. The new algorithm (we call it SCPFOSA) can be derived from equation (19), whose flow chart is illustrated in figure 2.
Step1: Perform a CZT across 1 m , get the azimuth coarse resolution estimate x s . Step2

Simulation Results
In this section, point target simulation is employed to validate the presented algorithm. The waveform parameters are chosen as: 250 The other parameters are listed as follow: standoff range is 10km, azimuth resolution is 0.1m, and radar forward velocity is 150m/s. Two point targets are simulated, the first one is the scene center point, and the other one is located at azimuth 150m away from the scene center. The new algorithm is evaluated with respect to the classical PFOSA which doesn't compensate the error terms resulted from the use of stepped-chirp signals. In subaperture algorithm, subapertures are overlapped to control the sidelobes; in particular, they are overlapped to control the amplitude to grating lobes due to data decimation. The degree of allowable overlap will depend on the window functions employed, and sidelobe toleration limits. In our paper, the overlap rate is not the problem we are concerned, so we do not employ window function in the simulation. As the phase error term  have almost the same response. However, for the azimuth displaced point target, the mainlobe of range profile is broadening for PFOSA, since the phase errors caused by motion of radar during pulse cluster are not compensated. While using SCPFOSA, due to the correction of these phase errors, its range profile has improved significantly (mainlobe reduce 12% and peak sidelobe ratio (PSR) reduces about 2.5dB).

Conclusion
In this paper, a 2-D subaperture algorithm based on stepped-chirp signal is presented. It integrates the pulse synthesis process into range subaperture processing without traditional pulse synthesis preprocessing. Meanwhile, due to the multi-resolution property of subaperture processing, this algorithm is able to compensate the space-variant phase error resulted from the motion of radar during a pulse cluster. SCPFOSA has almost the same processing flow chart with PFOSA, only the additional phase error term are added, it has the comparable computation complexity with PFOSA. Furthermore, due to its repetitive architecture in subaperture processing, the SCPFOSA is very suitable for parallel and pipeline hardware architectures.