Next Article in Journal
Sensing and 3D Mapping of Soil Compaction
Next Article in Special Issue
An Assessment of the Altimetric Information Derived from Spaceborne SAR (RADARSAT-1, SRTM3) and Optical (ASTER) Data for Cartographic Application in the Amazon Region
Previous Article in Journal
Motion Compensation of Moving Targets for High Range Resolution Stepped-Frequency Radar
Previous Article in Special Issue
SAR System for UAV Operation with Motion Error Compensation beyond the Resolution Cell
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Department of Electronic Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China
*
Author to whom correspondence should be addressed.
Sensors 2008, 8(5), 3438-3446; https://doi.org/10.3390/s8053438
Submission received: 3 April 2008 / Accepted: 19 May 2008 / Published: 26 May 2008
(This article belongs to the Special Issue Synthetic Aperture Radar (SAR))

Abstract

:
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.

1. Introduction

Synthetic aperture radar (SAR) becomes an important tool in modern remote sensing for its all-weather, 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-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 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.

2. 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 rcn. In this paper, we account for the broadside mode only, then α0 = 0. A target scatter is located at (sx, sy, 0), and the range from APC to this target is measured by rsn. 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]
f ( n , i ) = exp { j 4 π c ( f 0 + γ T s i ) R Δ }
where RΔ = (rcnrsn) is the differential range, c is the velocity of light, f0 is the carrier frequency of the transmitted signal, γ is the chirp rate, Ts is the sample interval in fast time, i is the range sample index with −I/2≤iI/2 − 1, and n is the index value of pulse number with −N/2 ≤ nN/2 −1.
The differential range RΔ can be expressed as [7]
R Δ = s x cos φ n sin α n s y cos φ n cos α n + r e = cos φ n cos α n ( s x tan α n s y + ξ ( n ) )
where ξ ( n ) = r e cos φ n cos α n, which is caused by the assumption of planar wavefront. By using the Taylor expansion, ξ(n) can be expressed as ξ (n) ≈ ε0 + ε1n + ε2n2 when the cubic and higher order terms are ignored. Since the coefficients of this polynomial are dependent on the target position, this error is space-variant.
Inserting Equation (2) into Equation (1), the sampled signal can be represented as
f ( n , i ) = exp { j 4 π c ( f 0 + γ T s i ) cos φ n cos α n ( s x tan α n s y + ξ ( n ) ) } = exp { j κ n ( 1 + β 0 i ) ( s x tan α n s y + ξ ( n ) ) }
where κ n = 4 π f 0 c cos φ n cos α n and β 0 = γ T s f 0.
Performing range resampling formulated by
κ n ( 1 + β 0 i ) = κ 0 ( 1 + β 0 i )
where κ 0 = 4 π f 0 c cos φ 0 cos α 0 is a constant, we can have the phase history as
f r ( n , í ) = exp { j κ 0 ( 1 + β 0 í ) ( s x tan α n s y + ξ ( n ) ) } .
For the purpose of clearness and simplicity, we still use i instead of í in the following discussion. Then, Equation (5) is expressed as
f r ( n , i ) = exp { j κ 0 ( 1 + β 0 i ) ( s x tan α n s y + ξ ( n ) ) }
If the space sampling position n satisfied tan αn = dαn, where is constant, the range resampled signal can be modeled as
f r ( n , i ) = exp { j κ 0 ( 1 + β 0 i ) ( s x d α n s y + ξ ( n ) ) } .

3. Two Dimensional Overlapped 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
n = m 1 + Δ 2 m 2 i = k 1 + μ 2 k 2
where m1 is the azimuth intra-subaperture index limited within −M1 /2≤ m1M1 /2 − 1, m2 is the azimuth inter-subaperture index limited within −M2/2≤m2M2/2 − 1, Δ2 is the azimuth data decimation factor, k1 is the range intra-subaperture index limited within −K1 /2 ≤ k1K1 / 2 − 1, k2 is the range inter-subaperture index limited within −K2/2≤k2K2/2−1, and m2 is the range data decimation factor. Using Equation (8), we rewrite the Equation (7) as
f r ( m 1 , m 2 ; k 1 , k 2 ) = exp { j κ 0 ( 1 + β 0 ( k 1 + μ 2 k 2 ) ) ( s x d α ( m 1 + Δ 2 m 2 ) s y + ξ ( m 1 + Δ 2 m 2 ) ) } .
Next, applying the quadratic order approximation of ξ(n) and rearrange Equation (9) following the index m1, k1, m2, k2 sequentially, we get
f r ( m 1 , m 2 ; k 1 , k 2 ) = exp { j κ 0 ( s y + ɛ 0 ) } exp { j κ 0 ( 1 + β 0 ( k 1 + μ 2 k 2 ) ) ( ( s x d α + ɛ 1 + 2 ɛ 2 Δ 2 m 2 ) m 1 + ɛ 2 m 1 2 ) } exp { j κ 0 β 0 ( s x d α Δ 2 m 2 s y + ɛ 0 + ɛ 1 Δ 2 m 2 + ɛ 2 ( Δ 2 m 2 ) 2 ) k 1 } exp { j κ 0 ( 1 + β 0 μ 2 k 2 ) ( ( s x d α Δ 2 + ɛ 1 Δ 2 ) m 2 + ɛ 2 ( Δ 2 m 2 ) 2 ) } exp { j κ 0 β 0 μ 2 ( s y + ɛ 0 ) k 2 }
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 M1 and K1 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
f r ( m 1 , m 2 ; k 1 , k 2 ) = exp { j κ 0 ( 1 + β 0 ( k 1 + μ 2 k 2 ) ) ( s x d α m 1 ) } exp { j κ 0 β 0 ( s y ) k 1 + j ϕ e 1 } exp { j κ 0 ( 1 + β 0 μ 2 k 2 ) ( s x d α Δ 2 ) m 2 + j ϕ e 2 } exp { j κ 0 β 0 μ 2 ( s y ) k 2 + j ϕ e 3 }
where ϕe1 = κ0β0sxΔ2m2k1, ϕe2 = κ0 (1 + β0μ2k2)(ε1Δ2m2+ε22m2)2), and ϕe3 = κ0β0m2ε0k2 are the undesired terms, which should be compensated in this algorithm.
From Equation (11), the processing strategy is clear. First, perform a CZT across m1, following by phase correction ϕe1, and then perform a DFT across k1 to get the coarse resolution images. After the second phase error compensation ϕe2, we perform a CZT across m2 to get azimuth fine resolution. Finally, after compensating the third phase error term ϕe3, a DFT across k2 results in the final fine resolution image.

4. 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. The carrier frequencies distribute sequentially to keep the spectrums covering the desired bandwidth.
Now assume that each pulse cluster has K2 chirp segments each with bandwidth Bs, the carrier frequency of the middle subchirp is f0, and the step carrier frequency is Δff < Bs). Then the center frequency of the k2 th (−K2/2≤k2K2/2−1) subchirp is f0 + k2Δf. Analogous to Equation (3), after preprocessing, the k2th backscattered echo signal in the nth pulse cluster can be represented by
f ( n , k 1 , k 2 ) = exp { j 4 π c ( f 0 + k 2 Δ f + γ T s k 1 ) cos φ n , k 2 cos α n , k 2 ( s x tan α n , k 2 s y + ξ ( n , k 2 ) ) }
where k1 is the range sample index in each chirp segment, k2 is the chirp segment index, and n is the cluster index. It is important to note that the grazing angle φn,k2 and squint angle αn,k2 vary not only with index n, which is desirable, but also with k2, which is undesirable. The latter change of φn,k2 and αn,k2 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 cosφn,k2 cosαn,k2 in terms of tan αn,k2. From the geometry in Figure.1, it is easy to get the following relationship
cos φ n , k 2 cos α n , k 2 = 1 1 + tan 2 φ 0 , 0 + tan 2 α n , k
where φ0,0 is the grazing angle at the aperture center corresponding to αn,k2 = 0.
Since the wavefront curvature error term ξ(n, k2) does not play an important role in this development, it is neglected. Then inserting Equation (13) into Equation (12) and perform a Taylor series expansion, we get the signal phase history in Equation (12) approximated as
Φ κ o [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] ( s x tan α n , k 2 s y + 1 2 s y cos 2 φ 0 , 0 tan 2 α n , k 2 )
where κ o = 4 π f 0 c cos φ 0 , 0, β 0 = γ T s f 0 and μ 2 = Δ f γ T s.
As before, we assume that the space sampling position n satisfies tanαn,k2 = (nK2+k2) then Equation (14) can be divided into two parts
Φ = Φ basic + Φ err
where
Φ basic = κ o [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] [ s x d α n K 2 s y + 1 2 s y cos 2 φ 0 , 0 ( d α ) 2 ( n K 2 ) 2 ] Φ err = κ o [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] [ s x d α k 2 + 1 2 s y cos 2 φ 0 , 0 ( d α ) 2 ( 2 n K 2 k 2 + k 2 2 ) ]
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 (sx, sy), this phase error is space-variant. Analogous to phase error resulted from wavefront curvature, if not compensated, it set the focused scene size limit, and particularly troublesome as resolution approaches the nominal wavelength of the radar.
For Equation (12), after range resampling, the phase history becomes
f r ( n , k 1 , k 2 ) = exp { j κ o [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] [ s x tan α n , k 2 s y + ξ ( n , k 2 ) ] }
where κ o = 4 π f 0 c cos φ 0 , 0, β 0 = γ T s f 0 and μ 2 = Δ f γ T s.
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 αn,k2 and ξ(n, k2) varies with index k2, while in Equation (6) they keep constant at specified n.
Insert tanαn,k2 = (nK2 + k2) into Equation (17), and then divide the azimuth aperture into subapertures
f r ( m 1 , m 2 ; k 1 , k 2 ) = exp { j κ 0 [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] [ s x d α ( K 2 ( m 1 + Δ 2 m 2 ) + k 2 ) s y + ξ ( m 1 , m 2 , k 2 ) ] } = exp { j κ 0 [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] [ s x d α K 2 ( m 1 + Δ 2 m 2 ) s y + ξ ( m 1 , m 2 , k 2 ) ] } exp { j ϕ e r r ( k 1 , k 2 ) }
where ϕerr(k1, k2) = κo [1 + β0(k1 + μ2k2)] sxdαk2 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
f r ( m 1 , m 2 ; k 1 , k 2 ) = exp { j κ 0 [ 1 + β 0 ( k 1 + μ 2 k 2 ) ] ( s x d α K 2 m 1 ) } exp { j κ 0 β 0 ( s y ) k 1 + j ϕ e 1 + j ϕ e r r ( k 1 ) } exp { j κ 0 ( 1 + β 0 μ 2 k 2 ) ( s x d α K 2 Δ 2 ) m 2 + j ϕ e 2 } exp { j κ 0 β 0 μ 2 ( s y ) k 2 + j ϕ e 3 + j ϕ err ( k 2 ) }
where ϕe1 = κ0β0sxdαK2Δ2m2k1, ϕe2 = κ0 (1+β0μ2k2) [ε1 Δ2m2 + ε22m2)2], and ϕe3 = κ0β0μ2ε0k2 are undesired terms just like those in PFOSA which are introduced by wavefront curvature. While error phase terms ϕerr(k1) = κoβ0sxdαk2k1 and ϕ e r r ( k 2 ) = κ o s x d α k 2 + κ o β 0 μ 2 s x d α k 2 2, the two parts of ϕerr(k1,k2), 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 m1, get the azimuth coarse resolution estimate sx.
  • Step2: Use the estimate of sx to compensate the error phase terms ϕe1 and ϕerr (k1), and then perform a FFT across k1 to obtain the range coarse resolution estimate sy.
  • Step3: Use the estimate of sx and sy to correct the error term ϕe2, then perform a CZT across m2 to get the azimuth fine resolution estimate sx.
  • Step4: Use the fine resolution estimation of sx and coarse resolution estimate of sy to compensate the error terms ϕe3 and ϕerr (k2), and then perform a FFT across k2 to get the range fine resolution estimate sy. The result is the fine resolution complex SAR image.

5. Simulation Results

In this section, point target simulation is employed to validate the presented algorithm. The waveform parameters are chosen as: Bs = 250MHz, Δf =100MHz, K2 =15. 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 ϕerr (k1, k2) results in distortion and defocus only in the range, we show the range profiles of impulse response function (IRF) for the two simulated targets (Figure. 3). For scene center point, since the error terms ϕerr (k1) and ϕerr (k2) are both zeros, the two algorithms 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).

6. 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.

References

  1. Cantalloube, H. Airborne X-band SAR Imaging with 10cm Resolution: Technical Challenge and Preliminary Results. IEE Proc-Radar Sonar Navig 2006, 153(2), 163–176. [Google Scholar]
  2. Ender, J.H.G. PAMIR – A Wide Phased Array SAR/MTI System. IEE Proc-Radar Sonar Navig. 2003, 150(3), 165–172. [Google Scholar]
  3. Lord, R.T.; Inggs, M.R. High Resolution SAR Processing Using Stepped-frequencies. IGARSS 1997, 1, 490–492. [Google Scholar]
  4. Inggs, M.R.; Van Zyl, M.W.; Knight, A. A Simulation of Synthetic Range Profile Radar, Proceedings of the 1992 South African Symposium on Communications and Signal Processing, Rondebosch, South Africa, September 1992; pp. 1–6.
  5. Schiimpf, H.; Wahlen, A. High Range Resolution by Means of Synthetic Bandwidth Generated by Frequency-stepped Chirps. Electronics Letters 2003, 39(18), 1346–1348. [Google Scholar]
  6. Wilkinson, A.J.; Lord, R.T.; Inggs, M.R. Stepped-frequency Processing by Reconstruction of Target Reflectivity Spectrum, Proceedings of the 1998 South African Symposium on Communications and Signal Processing, Rondebosch, South Africa, September 1998; pp. 101–104.
  7. Carrara, W.G.; Goodman, R.S.; Majewski, R.M. Spotlight Synthetic Aperture Radar: Signal Processing Algorithms; Artech House: Boston, USA, 1995. [Google Scholar]
  8. Jakowatz, C.V.; Wahl, D.E. Spotlight-mode Synthetic Aperture Radar: A Signal Processing Approach; Kluwer Academic Publishers: Boston, USA, 1996. [Google Scholar]
  9. Doerry, A.W. Synthetic Aperture Radar Processing with Polar Formatted Subapertures, Conference record of the 28th Asilomar conference on signals, systems&computers, Pacific Grove, California; 1994; pp. 1210–1215.
Figure 1. SAR geometry.
Figure 1. SAR geometry.
Sensors 08 03438f1
Figure 2. Flow chart of SCPFOSA.
Figure 2. Flow chart of SCPFOSA.
Sensors 08 03438f2
Figure 3. Comparison of SCPFOSA and PFOSA. (a) Range profile of the IRF for scene center point. (b) Range profile of the IRF for azimuth displaced target.
Figure 3. Comparison of SCPFOSA and PFOSA. (a) Range profile of the IRF for scene center point. (b) Range profile of the IRF for azimuth displaced target.
Sensors 08 03438f3

Share and Cite

MDPI and ACS Style

Mao, X.; Zhu, D.; Nie, X.; Zhu, Z. A Two Dimensional Overlapped Subaperture Polar Format Algorithm Based on Stepped-chirp Signal. Sensors 2008, 8, 3438-3446. https://doi.org/10.3390/s8053438

AMA Style

Mao X, Zhu D, Nie X, Zhu Z. A Two Dimensional Overlapped Subaperture Polar Format Algorithm Based on Stepped-chirp Signal. Sensors. 2008; 8(5):3438-3446. https://doi.org/10.3390/s8053438

Chicago/Turabian Style

Mao, Xinhua, Daiyin Zhu, Xin Nie, and Zhaoda Zhu. 2008. "A Two Dimensional Overlapped Subaperture Polar Format Algorithm Based on Stepped-chirp Signal" Sensors 8, no. 5: 3438-3446. https://doi.org/10.3390/s8053438

APA Style

Mao, X., Zhu, D., Nie, X., & Zhu, Z. (2008). A Two Dimensional Overlapped Subaperture Polar Format Algorithm Based on Stepped-chirp Signal. Sensors, 8(5), 3438-3446. https://doi.org/10.3390/s8053438

Article Metrics

Back to TopTop