Abstract

In the squinted synthetic aperture radar (SAR) imaging of the near-field environment, range-dependent characteristic of squint angle cannot be ignored, which causes azimuth-dependent range cell migration (RCM) after linear range walk correction (LRWC). In this study, an efficient SAR imaging algorithm applied in the near-field environment is proposed. In the processing of the range focusing, LRWC is firstly used to remove the linear RCM. Then, the residual LRCM is expanded into azimuth-invariant and azimuth-variant terms in consideration of the residual LRCM of azimuth-dependent. Range cell migration azimuth scaling (RCMAS) is designed to remove the azimuth-variant term before secondary range compression (SRC) and range compression (RC). In the azimuth focusing, azimuth distortion compensation (ADC) is performed to compensate the azimuth distortion, following which azimuth nonlinear chirp scaling (ANCS) is applied to equalize the frequency modulation (FM) rate for azimuth compression (AC). The simulated results show that more accurate and improved imaging result can be obtained with the proposed algorithm.

1. Introduction

Due to the ability of high-resolution imaging for the observed scene with any time and any weather, synthetic aperture radar (SAR) has become one of the most attractive radar techniques [1, 2]. There are many imaging modes for SAR, such as stripmap [3, 4], spotlight [5, 6], and sliding spotlight [7, 8]. Particularly, the squinted SAR in the stripmap is widely applied in military and civilian applications. By adjusting the squint angle, SAR can flexibility choose the imaging scene to meet diverse needs. However, the squinted SAR also brings new difficulties and challenges. As the squinted angle increases, the signal property varies a lot and the accurate focusing becomes difficult [1, 9, 10]. Therefore, it is necessary to study imaging algorithm for high-resolution squinted SAR.

A lot of algorithms have been proposed in the past years, for example, range-Doppler (RD) algorithm [1113], chirp scaling (CS) algorithm [1419], omega-K algorithm [20, 21], and backprojection (BP) algorithm [22, 23]. Theoretically, omega-K algorithm and BP algorithm can be applied in any imaging geometry and any squint angle without any approximation. However, the low efficiency and high computational load limit its practical application. As the squint angle increases, the range cell migration (RCM) and range-azimuth coupling become severe. This degrades the performance of RD algorithm. CS algorithm and its modifications have shown significant potential to process high-resolution squinted SAR data. Sun et al. [15] analyzed the spectrum of squinted data and used the squint minimization processing proposed in [24] to reduce the range-azimuth coupling. Then, an azimuth scaling function was created to equalize the Doppler rate of the targets in the same range cell. Liu et al. [16] first used linear range walk correction (LRWC), secondary range compression (SRC), and range compression (RC) to accomplish the range focusing. Then, in order to solve the problem of the azimuth-dependent frequency modulation (FM) rate caused by LRWC, a filter function and a scaling function, called the nonlinear chirp scaling (NCS) algorithm, were constructed. Chen et al. [17] built the accurate signal model of the missile-borne SAR which took into account not only the missile velocity but also the acceleration. In addition, the azimuth distortion of image was considered in the modified azimuth NCS algorithm. Li et al. [18] and An et al. [19] analyzed the characteristics of the azimuth phase terms in the case of bistatic and monostatic SAR, respectively. The modified scaling coefficients were derived by eliminating the higher-order phase approximation and azimuth-independent cubic phase term.

All of the above proposed algorithms processed the squinted SAR data acquired in a far-field environment. In this case, the wavefront of the electromagnetic wave transmitted by a radar was usually considered as a plane because of the very large distance between the radar and target. However, in a near-field environment, the spherical wavefront of the electromagnetic wave cannot be ignored. This will result in the squint angle of target varying with the range. Therefore, the imaging algorithm of high-resolution squinted SAR in the near-field environment needs to be redesigned.

The remainder of this paper is organized as follows. The radar signal model is described in Section 2. Section 3 gives the proposed imaging algorithm including seven procedures. Simulation results are provided to verify the effectiveness of the proposed algorithm in Section 4. Finally, this study is concluded in Section 5.

2. Signal Model

The general geometry of squinted SAR imaging in the near-field environment is shown in Figure 1. The radar moves along the direction parallel to the X-axis at a constant velocity . Or is the synthetic aperture center where the azimuth time ta is equal to zero. θ is the radar squint angle. The imaging scene is located at OXZ plane. L is the distance between the radar trajectory and OXZ plane, and H is the height of the imaged scene center. r0 is the slant range from the target to the radar trajectory along the radar line-of-sight direction. R is the instantaneous slant range between the radar and the target T when the radar is at azimuth time ta. OrF is the half of the synthetic aperture length, whose value is vta. The target T is located at the azimuth time tn, which satisfies HT = vtn. The instantaneous squint angle of target T is denoted by θq, which is expressed as . The targets A∼E are used in the following simulation.

In Figure 1, the instantaneous range R and its Taylor expansion at ta = tn are given bywhere k1, k2, k3, and k4 are specifically expressed as

Assuming a Gaussian pulse is transmitted by the radar, the demodulated radar echo signal received from target is given bywhere tr is the fast time, ta is the slow time, is the range time envelope, is the azimuth time envelope, c is the speed of light, is the Gaussian pulse, and fc is the carrier frequency.

As shown in Figure 1, the instantaneous squint angle θq of target T satisfies . It is obvious that θq is dependent on slant range r0 when L and θ are invariant. In order to clearly explain the influence of range-dependent squint angle on the RCM, simulation with the parameters listed in Table 1 was carried out. Supposed that there were two target points A (1.15 m, 0 m, 0.10 m) and C (1.15 m, 0 m, 1.50 m). As shown in Figure 1, the targets A and C were in the lowest and highest position of the imaging scene, respectively. They were located in the same azimuth cell. The squint angles of A and C were 29.97° and 24.79°, respectively. Figure 2 shows linear RCM (LRCM) and LRCM error between targets A and C. Figure 3 shows quadratic RCM (QRCM) and QRCM error between targets A and C. In Figure 2(b), LRCM error varied linearly with the azimuth time. The maximum LRCM error was approximately equal to the range resolution cell (0.0882 m). Therefore, in order to acquire the high-resolution image, the impact of range-dependent squint angle on the LRCM must be taken into account. In Figure 3(b), QRCM error was much smaller than the range resolution cell. Therefore, the variance of the QRCM and higher-order RCM caused by squint angle can be ignored.

3. Imaging Algorithm

The proposed algorithm includes the range focusing and the azimuth focusing, as shown in Figure 4. The range focusing consists of four steps. First, LRWC is used to remove most of the LRCM. Second, azimuth-dependent residual LRCM is removed through range cell migration azimuth scaling (RCMAS) processing. Third, SRC is applied to remove the coupling of range and azimuth. Finally, the range focusing is accomplished by RC processing. The azimuth focusing is resolved by three steps. The azimuth distortion compensation (ADC) is firstly performed to compensate the azimuth distortion. Then, the azimuth-dependent FM rate is equalized by azimuth NCS (ANCS). Finally, AC is operated to finish the azimuth focusing. The detailed derivations and explanations are introduced in the following.

3.1. Linear Range Walk Correction (LRWC)

Substituting (1) into (3) and applying the range fast Fourier transform (FFT) yieldswhere fr is the range frequency, is the range frequency envelope, and is the frequency response of the Gaussian pulse. The second exponential term in (4) represents the LRCM term, which is the major component of RCM. The third and fourth exponential terms are the QRCM and higher-order RCMs. According to (4), LRWC function is given as follows:where θs is the squint angle of the reference point. Multiplying (4) with (5) yieldswhere rLRCM = r0 + sinθstn is the new slant range corresponding to r0 after LRWC processing. k1res = (−sinθq) − (−sinθs) is the coefficient of residual LRCM.

In (6), the existence of coefficient k1res indicates that LRCM has not been completely removed after LRWC. Figure 5 shows the schematic diagram of LRWC processing. There are five targets P0∼P4, and target P0 is regarded as the reference point. Targets P1, P0, and P2 have the same beam center crossing time ta = 0, and targets P1 and P3 and P2 and P4 are located in the same range cell, respectively. In Figure 5(b), the black dotted lines represent the target LRCMs. LRCM trajectories of P1 and P3 have the same slope because these two targets have the same slant angle. This can also be explained for P2 and P4. Since the target squint angle varies with the slant range, targets P1, P0, and P2 have different LRCMs. As shown in Figure 5(c), the red dotted line denotes LRWC function. After LRWC processing, targets P3, P0, and P4 are located in the same range cell. However, targets P3 and P4 still have the residual LRCMs, and the slopes of their residual LRCM trajectories are different. This problem will be discussed in the next section.

3.2. Range Cell Migration Azimuth Scaling (RCMAS)

From previous analysis, after LRWC processing, the residual LRCMs of targets P3, P0, and P4 located in the same range cell are azimuth-dependent. Not only the residual LRCM but also the QRCM and higher-order RCM vary with the target azimuth position. As introduced in Section 2, azimuth-dependent residual LRCM needs to be considered, and azimuth-dependent QRCM and higher-order RCM are too smaller which are ignored. Therefore, in order to equalize the azimuth-dependent residual LRCM, k1res is expanded in terms of tn as follows:where k1c is the azimuth-invariant term and k1a is the coefficient of the azimuth-variant term. Based on (7), a perturbation function is designed to remove the azimuth-variant term of residual LRCM as follows:where p2 is the coefficient in the perturbation function and will be determined later. Multiplying (6) with (8) yields

Then, making and substituting (7) into (9) yieldswhere is expressed as

In order to remove the azimuth-dependent residual LRCM, the azimuth-variant coefficients k1a should be eliminated. In other words, the azimuth-variant terms should be zero in (11), and it satisfies the following equation:

Through solving (12), the coefficient p2 in the perturbation function is obtained as

Substitute (13) into (10) and rewrite (10) aswhere the coefficients k1n, k2n, k3n, and k4n are expressed as

In (14), the first exponential term is a constant which depends on target position. This term has no effect on the focusing process. It can be ignored if a magnitude image is the final product [14, 25]. Therefore, this term is ignored in the following algorithm.

3.3. Secondary Range Compression (SRC)

With the principle of stationary phase (POSP) and the method of series reversion (MSR) [26, 27], (14) is transformed into the 2D frequency domain as follows:where is expressed as (17), and the second exponential term represents the target azimuth position.

Expanding (17) into a power series of fr yieldswhere , , , and can be found in Appendix.

In (18), is the azimuth modulated term, contains the information of residual RCM, represents the coupling of range and azimuth, also called secondary SRC term, and (n ≥ 3) shows the higher-order cross coupling terms. In this paper, the range of reference point was used to realize the range compensation, and SRC function is given as follows:

Multiplying (19) with (16) yields

3.4. Range Compression (RC)

RC is a procedure in which the echo signal is multiplied by the conjugate of the radar transmitting signal to improve the range resolution. In this study, RC function is given as follows:where is the conjugate of the Gaussian pulse. Multiplying (21) with (20) and transforming the result into the range time domain by range inverse FFT (IFFT) yieldswhere Br is the signal bandwidth.

3.5. Azimuth Distortion Compensation (ADC)

Rewriting as and expanding into Taylor’s series of fa at fa = 0 yieldswhere ψ0, ψ1, ψ2, ψ3, and ψ4 are expressed in Appendix. The first term is the inherent phase information, the second term represents the azimuth distortion, the third one is the azimuth modulation term, and the fourth and fifth show the cubic and quartic phase terms. The azimuth distortion term may cause the deviation of the target azimuth position. In order to eliminate the azimuth distortion term, the azimuth distortion compensation function is described as follows:

Substituting (23) into (22) and then multiplying with (24) yieldswhere Ka is the azimuth FM rate expressed as follows:

3.6. Azimuth Nonlinear Chirp Scaling (ANCS)

It can be seen that the azimuth FM rate Ka contains k1n, k2n, k3n, and k4n, which are expressed in (15). Except k1n, which is a constant, the other three are azimuth-dependent. Therefore, Ka also varies with the target azimuth position. It is necessary to model the azimuth-dependent Ka. First, k2n, k3n, and k4n are expanded in terms of tn as follows:where k2nc, k3nc, and k4nc are the azimuth-invariant terms and k2na, k3na, and k4na are the coefficients of the azimuth-variant terms. Then, Ka is expanded in terms of tn and kept up to the second-order as follows:

It can be seen from (28) that the azimuth-dependent characteristic of the azimuth FM rate is very similar to the range-dependent characteristic of range FM rate. Thus, the range nonlinear chirp scaling algorithm can also be applied to azimuth signal for equalizing the azimuth FM rate, which has been described in [1518]. The derivation of ANCS is expressed in the following.

First, the cubic phase term ψ3 is evaluated along the azimuth position as follows:where ψ3c is the azimuth-invariant term and ψ3a is the coefficient of the azimuth-variant term. In order to eliminate the azimuth-independent cubic phase term and quartic phase term, the phase compensation factor is given as follows:

Multiplying (30) with (25) yields

Then, a fourth-order azimuth filter function is designed to filter the phase term, and this function is shown as follows:where Y3 and Y4 are the coefficients which will be determined later. Multiplying (32) with (31) yields

Transforming (33) into the azimuth time domain yields

Then, a fourth-order chirp scaling factor is constructed to eliminate the azimuth-dependent azimuth FM rate, and this factor is given by

Multiplying (35) with (34) yields

The POSP is applied to transform (36) into the azimuth frequency domain:where

Substitute (28) into (38); is approximately expressed as a power series of tn and fa as follows:

In (39), the first term is the azimuth-independent phase modulation. The second term is the target real azimuth position. The third can cause geometric distortion along target azimuth direction. The fourth, fifth, and sixth are the azimuth-dependent phase modulation terms. The seventh is independent of the azimuth frequency. And the last one, which represents the summation of other expansion series terms, is so small that it can be neglected.

To eliminate the azimuth-dependent modulation and azimuth distortion, the coefficients C1, C2, C3, C4, and C5 are set as follows:where α is the scaling factor ().

Through solving (39) and (40), the coefficients q2, q3, q4, Y3, and Y4 are obtained:

Substituting (28) and (41) into (37) yields

It can be seen from (42) that the azimuth-dependent modulation and azimuth distortion have been completely eliminated. In (42), the first exponential term contains the azimuth position of target. The residual exponential term is the azimuth compression factor.

3.7. Azimuth Compression (AC)

AC is performed after ANCS processing. AC function is given by

Multiplying (43) with (42) and transforming the result into the azimuth time domain yieldswhere Ba is the Doppler bandwidth. From (44), it is seen that the target is focused at (tn/2α, rLRCM), which needs to be corrected by geometry correction to obtain the real position using

4. Simulation

To verify the effectiveness of the proposed algorithm, the simulation with parameters in Table 1 was carried out. The Gaussian pulse was defined as (46), and its waveforms in the time and frequency domain are shown in Figure 6.

There were five targets in the simulation, whose 3D coordinates were A (1.15 m, 0 m, 0.10 m), B (1.15 m, 0 m, 1.04 m), C (1.15 m, 0 m, 1.50 m), D (0.67 m, 0 m, 1.50 m), and E (1.64m, 0m, 0.10m), respectively. The distribution of A∼E in the azimuth and range direction is shown in Figure 7. The targets A, B, and C are located in the same azimuth cell. The targets D, B, and E will be lay in the same range cell after LRWC. The target B was chosen as the reference point. The coordinates of these five targets are set so that the performance of the proposed algorithm can be easily shown. The scaling factor α is set to 0.505. In this study, the imaging results are obtained without any processing of the weighting function or side-lobe control.

Figure 8 shows the RCM trajectories of targets D, B, and E before and after RCMAS. In Figure 8(a), it was observed that RCM trajectory of target B, as well as targets D and E, was approximated into a straight line, which illustrated LRCM was the major component of RCM. After the reference range LRWC, targets D, B, and E were located in the same range cell. LRCM of the target B had been eliminated. However, because of the range-dependent characteristic of the squint angle, targets D and E still had different residual LRCM. As shown in Figure 8(c), after RCMAS, RCM trajectories of D and E were azimuth-invariant, which proved that the azimuth-dependent residual LRCM caused by range-dependent squint angle had been eliminated.

In order to evaluate the performance of the proposed algorithm, two methods, the extended NCS (ENCS) algorithm in [19] and the proposed algorithm, were utilized to process the simulated SAR data. Figure 9 shows the imaging results using these two algorithms. The color in Figure 9 represents the scattering intensity of target. In the simulation imaging, the deep red denoted the target point, and the dark blue showed there was no target. Because of the influence of side lobe, there also existed light blue. It can be easily seen that five targets were focused on the true position in Figure 9(b). However, in Figure 9(a), targets A, C, D, and E all had the azimuth shift compared with the true position. Targets A and E shifted along the azimuth positive direction, and C and D shifted along the azimuth negative direction. In the ENCS algorithm, because the range-dependent characteristic of squint angle was not considered, there was no the residual LRCM term after LRWC, that is, k1res in (6) was equal to zero. This could cause that ψ1 was equal to zero and azimuth distortion compensation processing was not operated. Therefore, the focused targets had the azimuth shift using ENCS algorithm. In addition, there was defocus in targets A, C, D, and E because their RCMs were not fully removed. The proposed algorithm improved the image quality as shown in Figure 9(b). The imaging result contained the range and azimuth information of targets. Table 3 shows the target position error in the simulated images using ENCS algorithm and the proposed algorithm. The range errors of D and E with the ENCS algorithm were much larger than those with the proposed algorithm because RCMs of D and E did not make a precise compensation. In the azimuth, because of the impact of the azimuth shift, the errors of targets A, C, D, and E with ENCS algorithm were much larger than those with the proposed algorithm. This proved that the proposed algorithm performed well in the near-field environment.

5. Conclusion

In the near-field environment, the range-dependent squint angle of the target needs to be considered, which causes that there exists the residual LRCM of azimuth-dependent after reference range LRWC. In order to overcome the azimuth-dependent characteristic of the residual LRCM, the residual LRCM is expanded into azimuth-invariant and azimuth-variant terms. A novel perturbation function in RCMAS is designed to remove the azimuth-variant term. In the azimuth focusing, a filter function and a scaling function are applied to eliminate the high-order Doppler phases and equalize FM rate. The whole procedures only consist of FFT, IFFT, and multiplication operations without interpolation, which means higher efficiency and easier implementation. The simulation proves that more accurate and improved imaging result can be obtained with the proposed algorithm.

Appendix

, , , and in (18) are expressed as (A.1)∼(A.5). ψ0, ψ1, ψ2, ψ3, and ψ4 in (23) are expressed as (A.5).

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was financially supported by the funds of the Natural Science Foundation of China (grant nos. 51874279 and 51974290) and by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).