Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Comparative Study on Predicting Topography from Gravity Anomaly and Gravity Gradient Anomaly
Next Article in Special Issue
SAR Multi-Angle Observation Method for Multipath Suppression in Enclosed Spaces
Previous Article in Journal
A Lightweight Man-Overboard Detection and Tracking Model Using Aerial Images for Maritime Search and Rescue
Previous Article in Special Issue
A Modified Iteration-Free SPGA Based on Removing the Linear Phase
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Studies on High-Resolution Airborne Synthetic Aperture Radar Image Formation with Pseudo-Random Agility of Interpulse Waveform Parameters

College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(1), 164; https://doi.org/10.3390/rs16010164
Submission received: 2 November 2023 / Revised: 16 December 2023 / Accepted: 25 December 2023 / Published: 30 December 2023
(This article belongs to the Special Issue Advances in Synthetic Aperture Radar Data Processing and Application)
Figure 1
<p>Spotlight SAR data collection geometry.</p> ">
Figure 2
<p>Flowchart of airborne SAR image formation with agile interpulse waveform parameters.</p> ">
Figure 3
<p>4.5 kg microminiature SAR: (<b>a</b>) physical photo; (<b>b</b>) front view; (<b>c</b>) rear view; (<b>d</b>) assembly view.</p> ">
Figure 4
<p>Dual-polarization stripmap SAR image, 0.5 m resolution: (<b>a</b>) VV polarization; (<b>b</b>) VH polarization.</p> ">
Figure 5
<p>Real-time stripmap image of rural area at 0.3 m resolution.</p> ">
Figure 6
<p>Real-time spotlight images of urban area at 0.15 m resolution.</p> ">
Figure 7
<p>Agile interpulse waveform parameters at resolution 0.5 m: (<b>a</b>) pulse width; (<b>b</b>) PRI; (<b>c</b>) initial phase; (<b>d</b>) chirp rate polarity.</p> ">
Figure 8
<p>Amplitude equalization factor for compensating for the fluctuations in processing gain and average power resulting from the random agility of pulse width and PRI at a resolution of 0.5 m.</p> ">
Figure 9
<p>Rural area SAR images at 0.5 m resolution with agile interpulse waveform parameters: (<b>a</b>) no compensation for agility of pulse width or chirp rate; (<b>b</b>) no compensation for agility of chirp rate polarity; (<b>c</b>) no compensation for agility of initial phase; (<b>d</b>) no compensation for agility of PRI; (<b>e</b>) no compensation for all agile interpulse waveform parameters; (<b>f</b>) compensation for all agile interpulse waveform parameters according to the flowchart shown in <a href="#remotesensing-16-00164-f002" class="html-fig">Figure 2</a>.</p> ">
Figure 10
<p>The comparison of the chirp rate polarity of the jamming signal matches and mismatches the transmitted signal: (<b>a</b>) result of pulse compression and rotating-platform processing with no compensation of chirp rate polarity agility; (<b>b</b>) local enlarged result of the area marked by yellow dotted line box of (<b>a</b>); (<b>c</b>) result of pulse compression and rotating-platform processing with compensation of chirp rate polarity agility; (<b>d</b>) local enlarged result of the area marked by yellow dotted line box of (<b>c</b>).</p> ">
Figure 11
<p>Agile interpulse waveform parameters at resolution 0.15 m: (<b>a</b>) pulse width; (<b>b</b>) PRI; (<b>c</b>) initial phase; (<b>d</b>) chirp rate polarity.</p> ">
Figure 12
<p>Amplitude equalization factor for compensating the fluctuations of processing gain and average power resulted from the random agility of pulse width and PRI at a resolution of 0.15 m.</p> ">
Figure 13
<p>Spotlight SAR image of a power transmission line scene at 0.15 m resolution with agile interpulse waveform parameters.</p> ">
Review Reports Versions Notes

Abstract

:
By means of alteration of the transmitted linear frequency modulation (LFM) signal waveform parameters, such as pulse width or chirp rate, initial phase, pulse repetition interval (PRI), and chirp rate polarity at every position of synthetic apertures, the pseudo-random agility technology of interpulse waveform parameters in airborne Synthetic Aperture Radar (SAR) actively increases the complexity and uncertainty of radar waveforms. This technology confuses jamming interception receivers, thus improving its anti-interference ability for active coherent jamming, which is one of the main research interests of airborne SAR technology. But the pseudo-random agility technology for interpulse waveform parameters faces certain challenges of large computation and complex system design, which need to be further studied and solved. To address these issues, a processing scheme of high-resolution SAR image formation which is appropriate for agile interpulse waveform parameters is proposed in this paper. This method can deal with multiple agile parameters, not only single ones as in most existing literature. Its computation load is nearly comparable to that of traditional SAR image formation with constant waveform parameters. The high-resolution SAR imaging results obtained by processing SAR raw data with agile interpulse waveform parameters demonstrate the effectiveness of the proposed method. In addition, real SAR images with resolutions of 0.5 m and 0.15 m, which are rarely found in the public literature, are shown under the circumstance of randomly changing the transmitted wideband LFM signal pulse parameters one by one.

1. Introduction

Conventional airborne SAR transmits wideband LFM signals with constant waveform parameters during the synthetic aperture data collection period, which are easy to be intercepted, sorted and identified by the jammer. Thus, it has the risk of interference, which reduces the reconnaissance efficiency of airborne SAR in electronic countermeasure environments [1]. Therefore, electronic countermeasures (ECMs) [1,2,3,4] have become the main challenge faced by airborne SAR sensors. To address this issue, electronic counter-countermeasures (ECCMs) [5,6,7] have also become one of the main research interests in the field of airborne SAR. Among various SAR ECM technologies, the most threatening one is active coherent interference technology [4,8]. It usually obtains accurate radar waveform parameters by interception and modulates the phase to repeat the radar signal or independently generates false target jamming signals which are coherent with the real radar echo signals, so that after this kind of jamming signals enter the radar, they can obtain high-enough coherent processing gain in range and azimuth dimensions. In this way, a good interference effect can be achieved with low jamming power, so it is an ideal SAR jamming method at present. The literature [8] points out that multi-channel SAR can provide more spatial freedom, and it is a common active coherent interference suppression method [9]. But a multi-channel system will inevitably need more complex hardware and a heavier computation load, which is inconvenient for engineering implementation. The active waveform countermeasure method of single-channel SAR, which suppresses the active coherent interference by random agility of the interpulse waveform parameters, is a desirable, simple and effective anti-active coherent interference method. Generally, SAR transmits wideband LFM signals, and its waveform parameters include carrier frequency, bandwidth, pulse width, chirp rate, chirp rate polarity, initial phase and pulse repetition interval (PRI). While the bandwidth should be kept constant to ensure the equal range resolution of each pulse during the synthetic aperture data collection period, other parameters, such as pulse width or chirp rate [6], chirp rate polarity [10], initial phase [11] and PRI [12], can be set to random agility in principle. The random agility increases the complexity and uncertainty of the transmitted signal, confuses the jamming interception, disturbs the timing of coherent jamming repeater and destroys the coherence of the jamming signal in range and azimuth dimensions; finally, the purpose of suppressing coherent interference is achieved. The literature [13] points out that the pseudo-random agility technology for interpulse waveform parameters is faced with certain challenges of large computation and complex system design, which needs to be further studied and solved. Most of the existing works I the literature only study SAR anti-coherent interference technology [6,9,10,11,12] itself with random agility of interpulse waveform parameters. However, the influence of random agility of interpulse waveform parameters on airborne SAR image formation is ignored, and no definite conclusion, i.e., whether airborne SAR with random agility of interpulse waveform parameters can acquire high-quality high-resolution SAR images just like conventional airborne SAR with constant waveform parameters, was given. Although there are some works in the literature [12] that study SAR imaging based on non-uniform PRI, they mostly focus on a frequency domain reconstruction [14] for non-uniformly sampled signals in situations where the time offsets of each sampling instance are known to restore the uniform sampling, and then they carry out the normal imaging process. The validity of this reconstruction algorithm for non-uniformly sampled signals resulting from the random agility of interpulse PRI was not mentioned yet [12]. In addition, the process of transforming non-uniform azimuth sampling into the desired uniform azimuth sampling requires more computation and memory resources, which is also inconvenient for engineering implementation.
The purpose of this paper is to study a high-resolution airborne SAR image formation technique under the condition of pseudo-random agility of interpulse waveform parameters, and the main innovations are as follows.
(1)
Considering the pseudo-random agility technology of interpulse waveform parameters in SAR, only a single parameter or two parameters, such as pulse width and/or initial phase, was changed in most previous works. The compensation for agile interpulse waveform parameters in SAR imaging processing is relatively easy, but its anti-interference ability for active coherent jamming is limited. The influence of agile interpulse waveform parameters, such as pulse width, chirp rate polarity, initial phase and PRI, on SAR imaging processing and their compensation methods are analyzed in detail. The computation load of the method proposed in this paper was comparable to that of traditional SAR image formation with constant waveform parameters. Obviously, the anti-interference ability for active coherent jamming should be improved through simultaneous random agility of pulse width, chirp rate polarity, initial phase and PRI.
(2)
An innovatively designed and highly integrated Ku band 4.5 kg microminiature SAR sensor is briefly introduced, which consists of a co-aperture dual-polarized patch array antenna, two-axis gimbal, 80 W solid state power amplifier, dual-channel radio frequency (RF) receiver including a built in limiter, on-board processing board performing real-time image formation, 2 TB data storage board to store dual-channel raw data and real-time SAR images, necessary power conversion module and a commercial off-the-shelf micro electro mechanical systems (MEMS)-based inertial measurement unit (IMU)/Global Position System (GPS) for motion measurement. Its performance specifications and several example SAR images with constant interpulse waveform parameters are shown to illustrate the proper function mode and excellent behavior of this 4.5 kg microminiature SAR sensor. Considerations to upgrade this SAR with the capability of the random agility of interpulse waveform parameters are given.
The rest of this paper is structured as follows. The influences of agile interpulse waveform parameters on SAR image formation are discussed in detail in Section 2; then, a flowchart of airborne SAR image formation is proposed. An innovative 4.5 kg microminiature SAR and the upgrading considerations on it with the capability of random agility of interpulse waveform parameters are introduced briefly in Section 3. Experimental verification results of SAR imaging, which were obtained through processing this upgraded SAR raw data according to the proposed flowchart of airborne SAR image formation under the circumstance of randomly changing the transmitted wideband LFM signal waveform pulse parameters one by one, are shown in Section 4. Section 5 gives some discussion. Section 6 summarizes the whole paper.

2. Influences of Agile Interpulse Waveform Parameters on SAR Image Formation

At every pulse repetition frequency (PRF) moment, SAR transmits wideband LFM signals, whose waveform parameters, such as carrier frequency, pulse width or chirp rate, polarity, initial phase and PRI, can be randomly altered within a certain range, except that the bandwidth should be kept constant. The impacts of these variations in waveform parameters on the imaging process are analyzed as follows.

2.1. Agility of Carrier Frequency

If we change carrier frequency randomly within the tolerable range of radar hardware at every PRF moment, the anti-jamming ability of the radar can be improved in principle. Assuming that the nominal carrier frequency of the radar is f c , f c i denotes the carrier frequency of the ith pulse and can be expressed as
f c i = f c + Δ f i
where Δ f i is the random frequency offset of the ith pulse from the nominal carrier frequency.
The spotlight SAR data collection geometry is illustrated in Figure 1. Without the loss of generality, assume that the radar undergoes an arbitrary flight path. The instantaneous position of the radar antenna phase center (APC) is denoted by the vector r o t = x a t , y a t , z a t at slow time t. The variables θ = θ t and φ = φ t represent the instantaneous azimuth angle and grazing angle, respectively. Denote θ r e f as the reference azimuth angle and φ r e f as the corresponding grazing angle at slow time t = 0 . Consider a ground point target scatterer, whose coordinates are defined by the vector r c = x c , y c , 0 . The vector between the APC and this target scatterer is r t = r o t r c .
In the case of constant carrier frequency f c , and ignoring a minor phase term often referred to as the residual video phase error (RVPE), the 2-D phase history [15,16] delivered to the Polar Format Algorithm (PFA) after motion compensation with respect to the scene center is expressed as
S 1 t , f τ = C exp j 4 π f c + f τ c r o t r t
where t represents slow time, f τ is range frequency, c is the speed of propagation, and C includes the nonessential factors of the transmitted pulse envelop and azimuth antenna pattern.
Similarly, the 2-D phase history when the radar carrier frequency is agile can be expressed as
S 2 t , f τ = C exp j 4 π f c i + f τ c r o t r t     = C exp j 4 π f c + f τ c r o t r t exp j 4 π Δ f i c r o t r t
Comparing Equation (2) with Equation (3), we can see that the first phase term in the second line of Equation (3) is the same as Equation (2), which is the basic imaging term containing the useful information about locating and focusing the scene target. But the second phase term resulting from the random frequency offset Δ f i is undesirable and needs to be eliminated to restore the correct 2-D phase history.
At every synthetic aperture time instant of t, the random frequency offset Δ f i and the instantaneous distance r o t are fixed values, and the instantaneous distance r t is a variable related to the position of every target scatterer in the imaging scene, so the second phase term is space-variant. Obviously, the compensation of the second phase term to restore the correct 2-D phase history requires a heavy computation load, which makes it impractical for engineering implementation. So, we do not recommend changing carrier frequency randomly.

2.2. Agility of Pulse Width and Chirp Rate Polarity

Assume that the bandwidth B of the transmitted LFM signal is kept constant. a i denotes the chirp rate polarity, which is randomly +1 or −1, and the superscript i represents the ith pulse. Bandwidth is expressed as
B = K r i τ p i
where K r i and τ p i denote the chirp rate and pulse width, respectively. So, the range match filter H is
H r i f τ = exp j a i π K r i f τ 2
where f τ is the range frequency. As can be seen from Equation (5), the range match filter needs to be changed pulse by pulse according to the current chirp rate and chirp rate polarity. Obviously, this is different from the traditional range match filter with constant waveform parameters.
The random agility of the pulse width and chirp rate polarity will not affect the maximum range processing gain obtained by pulse compression of the echo signal, while the jamming signal cannot obtain sufficient range processing gain due to chirp rate or chirp rate polarity mismatch; active coherent interference can be suppressed effectively in this way [6].
Because the processing gain [17] of the range match filter is τ p i B , the agility of the pulse width will lead to fluctuation in the processing gain, resulting in non-uniform image intensity from pulse to pulse. Moreover, the non-constant pulse width τ p i and P R I i will lead to a variable duty factor, which is represented as
d i = τ p i PRI i
and it means the interpulse fluctuation of average power P a v i transmitted during the synthetic aperture data collection period. P a v i is expressed as
P a v i = P t d i
where P t represents the peak transmitter power and is usually a constant value, which will also lead to non-uniform image intensity from pulse to pulse. Therefore, the random agility of the pulse width or PRI will eventually affect the uniformity of the SAR image, which needs to be compensated for.
Assume that the pulse width τ p 0 and pulse repetition interval P R I 0 were chosen according to the requirements of the operating range and Signal to Noise Ratio (SNR) of the airborne SAR with constant waveform parameters. The processing gain equalization factor and the power equalization factor with agile waveform parameters are, respectively, expressed as
A 1 i = τ p i B τ p 0 B = τ p i τ p 0
A 2 i = d i d 0 = τ p i PRI 0 τ p 0 PRI i
Therefore, the equalization of the processing gain and average power fluctuations resulting from the random agility of pulse width and PRI can be achieved by multiplying the amplitude of the ith echo signal by a factor 1 A 1 i A 2 i .

2.3. Agility of Initial Phase

Airborne SAR can randomly change the initial phase ϕ n i of the transmitted wideband LFM signal on purpose at every PRF moment. ϕ n i represents the initial phase of the ith pulse, which is a uniformly distributed random number between 0 , 2 π and is statistically independent from pulse to pulse. After receiving the echo signal, the random initial phase can be compensated for by multiplying the echo signal by a phase exp j ϕ n i according to the prior information of the initial phase of the ith pulse, so as to restore the coherence of the azimuth signal during the synthetic aperture data collection period. The subsequent azimuth processing is the same as that of conventional SAR. This random initial phase will cause the azimuth phase of the jamming signal to be not perfectly coherent, and the azimuth processing gain cannot be obtained sufficiently, so that the jamming signal can be suppressed to a certain extent [11].

2.4. Agility of PRI

The agility of PRI makes it difficult for the jammer to apply fully coherent jamming signal with airborne SAR, resulting in a certain degree of ineffectiveness of the coherent jamming signal. However, for airborne SAR itself, the random agility of PRI causes random fluctuation in the duty factor, resulting in variation in average power, which leads to non-uniform image intensity from pulse to pulse. A compensation technique was given in Section 2.2. Moreover, non-uniform azimuth sampling caused by the random agility of PRI also leads to a poor azimuth sidelobe if the image formation process is not appropriate for non-uniform azimuth sampling, and neglecting this non-uniformity results in the deterioration of SAR image quality [12].
According to Ref. [15], the 2-D phase history in Equation (2) after range resampling in PFA is
S 3 t , f τ = C exp j 4 π f c + f τ c cos φ r e f x c cot θ + y c
where cot θ is a function of t. As analyzed in Section 2.3 of Ref. [15], even for constant PRI, which means uniformly spaced time instants of t, the corresponding values of cot θ are, in general, non-uniformly spaced for the realistic motion of the APC. However, the range-resampled phase history in Equation (10) can be interpolated on specific time instants to achieve uniformly spaced values for the term cot θ , which is a linear function of t.
Finally, after azimuth resampling [15], the 2-D phase history becomes
S 4 t , f τ = C exp j 4 π c cos φ r e f Ω x c f c t + y c f τ exp j 4 π f c c cos φ r e f y c
where Ω is a constant value [15] determined by the azimuth interpolation process.
In Equation (11), there has been no cross-product quadratic phase term of t and f τ . The 2-D sinusoid in Equation (11) can be converted into a focused target response via the 2-D Fourier transform with respect to t and f τ .
Obviously, the azimuth resampling procedure in PFA leads to azimuth sample spacing equalization in SAR data processing. So, the PFA is naturally suitable for non-uniform azimuth sampling caused by the random agility of PRI. The smooth embedment of the equalization into the image formation algorithm is a unique advantage that the PFA has over the other precise SAR image formation algorithms, such as the Range Migration Algorithm (RMA) [18,19,20].

2.5. Airborne SAR Image Formation

There are several typical and commonly used SAR image formation algorithms, such as the Range Doppler Algorithm (RDA) [18,19], Chirp Scaling Algorithm (CSA) [18,19,20,21], RMA [22] and PFA [23,24,25,26,27,28].
Considering that PFA performs phase adjustment and compensation in the time domain, it is not only simple, efficient and easy to implement in engineering, but it can also effectively compensate for the non-ideal trajectory of a radar platform and is naturally suitable for non-uniform azimuth sampling caused by the random agility of PRI [15]. Thus, the time-consuming processing step for transforming non-uniform azimuth sampling into the desired uniform azimuth sampling [14] can be omitted. However, PFA has the space-variant phase error caused by planar wavefront assumption [23,24,25], which limits the scene size of SAR images without defocus effects and causes the obvious geometric distortion of SAR images. Space-Variant Post-Filtering (SVPF) [16,25,26,27,28] is adopted to correct wavefront curvature and thus to enlarge the scene size limits and eliminate geometric distortion in polar-formatted SAR imagery. Therefore, we propose to use PFA + SVPF + PGA [29,30,31] as a processing scheme of high-resolution SAR image formation, which is appropriate for agile interpulse waveform parameters.
The proposed imaging process flowchart is shown in Figure 2, in which random initial phase and pulse-to-pulse amplitude fluctuation that are caused by variable pulse width and PRI are compensated for during preprocessing, as shown in the blue dotted line box. Pulse compression and rotating platform compensation are shown in the green box, which completes the range match filter, as shown in Equation (5), and stripmap-to-spotlight conversion. It should be noted that the range match filter needs to be changed pulse by pulse according to the interpulse agility of pulse width and chirp rate polarity, which is different from traditional pulse compression with constant waveform parameters. The standard polar format processing flowchart [32] is shown in the blue box, and SVPF plus PGA are shown in the yellow box to accomplish wavefront curvature correction [16,28] and phase gradient autofocus [33]. It can be seen that compared with conventional airborne SAR processing, only the random initial phase and amplitude compensation in the blue dotted line box are added as a preprocessing step, and the range match filter in the green box needs to be calculated pulse by pulse, so its resulting computation can be neglected.

3. Considerations on Upgrading A SAR with Agile Interpulse Waveform Parameters

3.1. Brief Introduction of 4.5 kg Microminiature SAR

A Ku band microminiature SAR is chosen to improve the capability of pseudo-random agile interpulse waveform parameters. This SAR was innovatively designed and integrated all of the following function modules into a single unit, as shown in Figure 3: an antenna/gimbal, transmitter/receiver, processing/data storage board, power conversion module and MEMS based IMU/GPS. Figure 3a is a physical photo, Figure 3b is the front view, Figure 3c is the rear view and Figure 3d is the assembly view of this microminiature SAR. In addition, the inclusion of a host computer and display (either onboard or via a remote link) completes this SAR system. Table 1 outlines the top-level performance specifications for this SAR system.
This Ku band 4.5 kg microminiature SAR has been flight-demonstrated on a variety of manned aircraft, such as the Y12, Cessna 208 B and Robinson R44, and small unmanned aerial vehicles (SUAVs), such as fixed-wing or multi-rotor UAVs. A large amount of high-resolution stripmap and spotlight dual-polarization SAR raw data and images were acquired. Here, we show just a few examples of the data product, illustrating the proper function mode and excellent performance of this SAR sensor with constant waveform parameters.
The 0.5 m × 0.5 m resolution dual-polarization stripmap SAR images are shown in Figure 4, which were obtained in the first flight demonstration of this 4.5 kg microminiature SAR in September 2020 on a Y12. From the rural area, we can see the obvious RCS difference between the VV-polarized SAR image shown in Figure 4a and the VH-polarized SAR image shown in Figure 4b.
Next, the real-time stripmap SAR image of a rural area at 0.3 m × 0.3 m resolution generated by this radar sensor is shown in Figure 5, which was flight tested on a small fixed-wing UAV.
Finally, four real-time spotlight SAR images of an urban area at 0.15 m × 0.15 m resolution acquired from different view angles of the same track are shown in Figure 6, whose scene center is the ground parking lot of a hotel. Their corresponding azimuth angle is 50°, 53.71°, 57.81° and 62.31°, respectively.
It can be seen that the SAR images of different resolutions shown above from Figure 4, Figure 5 and Figure 6 are well focused, and they have clear and rich texture and uniform intensity under the condition of constant waveform parameters.

3.2. Upgrading the 4.5 kg Microminiature SAR with Agile Interpulse Waveform Parameters

Considering its wideband LFM signal is generated by a digital waveform synthesizer, it is convenient to upgrade this radar with the capability of pseudo-random agility of interpulse waveform parameters such as pulse width, chirp rate polarity, initial phase and PRI at every PRF moment. A 0.5 m and 0.15 m resolution mode was chosen to be upgraded. Table 2 is the original waveform parameters of 0.5 m and 0.15 m resolution mode.
When the radar waveform parameters are agile, in order to ensure that the duty factor change caused by variable pulse width or PRI do not cause too much change in radar range distance or SNR, it is necessary to limit the variation range of the pulse width and PRI. According to the radar equation [17,22] of SAR, R 3 P a v d , if the variation of radar range is required to be less than 10%, then the estimated allowable variation range of pulse width and PRI is about 9~11 μ s and 900~1100 μ s , respectively. Therefore, we have determined the random agility mode for pulse width, PRI, initial phase and chirp rate polarity. Table 3 is a description of four agile parameters in our experiment.

4. Experimental Verification of SAR Image Formation with Agile Interpulse Waveform Parameters

The 4.5 kg microminiature SAR was upgraded according to the considerations introduced in Section 3.2. Owing to no suitable jammer at hand, and given that the main purpose of this work is to study and verify whether airborne SAR with agile interpulse waveform parameters can acquire high-quality high-resolution SAR images just like conventional airborne SAR with constant waveform parameters, this upgraded 4.5 kg SAR was flight tested independently, and SAR raw data with resolutions of 0.5 m and 0.15 m were collected with agile interpulse waveform parameters.
The variations in the agile pulse width, PRI, initial phase and chirp rate polarity of the first 1000 pulses in one frame of raw data collected by the upgraded 4.5 kg SAR with a resolution of 0.5 m are shown in Figure 7. It can be shown that four waveform parameters, including pulse width, PRI, initial phase and chirp rate polarity, are agile simultaneously at every PRF moment; therefore, the complexity and uncertainty of radar waveform can be greatly increased, and the anti-interference ability for active coherent jamming can be significantly improved in principle.
The amplitude equalization factor required for compensating the fluctuations in processing gain and average power resulting from the random agility of pulse width and PRI is shown in Figure 8. It can be found that the range of the amplitude equalization factor is about 0.87 to 1.17.
The 0.5 m resolution rural area SAR images are shown in Figure 9. Figure 9a is a small area from the imaging result obtained according to the processing flowchart shown in Figure 2, except that the agility of the pulse width or chirp rate is not compensated for, which corresponds to the white dotted line box in the SAR image given in Figure 9f. Similarly, Figure 9b–d are the same area when the agility of chirp rate polarity, initial phase and PRI are also not compensated for, respectively. It can be seen that Figure 9a is severely defocused, which is due to the mismatch of chirp rate resulting from pulse width agility, and Figure 9b–d have various degrees of SNR deterioration. The SNR deterioration in Figure 9c is due to the incoherence in the azimuth caused by initial phase agility, and the image deterioration in Figure 9d is due to poor sidelobe in the azimuth caused by PRI agility. The reason for the image deterioration shown in Figure 9b will be explained in the following section. Figure 9e shows the rural area SAR image obtained by processing the SAR raw data with no compensation for agile interpulse waveform parameters; obviously, it is completely defocused. Figure 9f shows the rural area SAR image obtained by processing the same SAR raw data with compensation for agile interpulse waveform parameters according to the proposed processing flowchart shown in Figure 2. It can be seen that the SAR image shown in Figure 9f is well focused and has clear and rich texture and uniform intensity.
Due to the random agility of chirp rate polarity, when the chirp rate polarity of the jamming signal matches the transmitted signal, the jamming signal can normally obtain the range processing gain, but when it is mismatched, the jamming signal cannot obtain the processing gain. The result of pulse compression and rotating-platform processing with no compensation for chirp rate polarity agility is shown in Figure 10a, and the result of pulse compression and rotating-platform processing with compensation for chirp rate polarity agility is shown in Figure 10c. Figure 10b,d are the local enlarged results of their corresponding area marked by yellow dotted line boxes. Comparing the vertical bright line in Figure 10b with that in Figure 10d, we can see that the vertical lines are intermittent and continuous, respectively. It can be found that the processing result as indicated by Figure 10a resulted from a chirp rate polarity mismatch and is equivalent to the random loss of samples or unequally spaced sampling in the azimuth dimension, so after normal azimuth processing, the image will be deteriorated by azimuth mainlobe broadening and a poor sidelobe, just like the effect shown in Figure 9b.
The contrast value adopted in the objective evaluation is as follows:
C = σ I m , n E I m , n = E I m , n E I m , n 2 E I m , n
where σ ( ) represents the standard deviation of the image intensity, E ( ) represents the mean of the image intensity and I m , n represents the intensity of the pixel at (m, n) in the image. The higher the contrast value, the higher the image quality.
The contrast values corresponding to Figure 9a–d and the same area indicated by the white dotted line box in Figure 9f are presented in Table 4. It can be seen that the area indicated by the white dotted line box in Figure 9f has the highest contrast value, which is at least twice higher than that of Figure 9a–d.
The variations in the agile pulse width, PRI, initial phase and chirp rate polarity of the first 1000 pulses in one frame of raw data collected by the upgraded 4.5 kg SAR with resolution 0.15 m are shown in Figure 11.
The amplitude equalization factor required for compensating the fluctuations in processing gain and average power resulting from the random agility of pulse width and PRI is shown in Figure 12. As can be seen from Figure 12, the range of the amplitude equalization factor is also about 0.87 to 1.17.
According to the same processing flow shown in Figure 2, not surprisingly, the same high-quality SAR image with 0.15 m resolution can be obtained. Figure 13 shows a rural area scenario, with the power transmission line, the pylon reflection, its shadow and power transmission line shadow clearly visible in the spotlight SAR image at a resolution 0.15 m with agile interpulse waveform parameters.
The high-resolution SAR imaging results shown in Figure 9f and Figure 13, and the contrast values shown in Table 4, indicate that airborne SAR can acquire the same satisfactory high-quality SAR images definitely under the condition of random agility of interpulse waveform parameters, just as the images shown in Figure 4, Figure 5 and Figure 6 with constant waveform parameters.

5. Discussion

In most previous works, the active waveform countermeasure method of SAR generally uses only single agile parameter or two agile parameters such as pulse width and/or initial phase, which is relatively easy for engineering implementation, but with limited anti-active coherent interference ability. The simultaneous random agility of four waveform parameters of SAR, namely pulse width, chirp rate polarity, initial phase and PRI, was proposed and verified in this paper; therefore, the complexity and uncertainty of the radar waveform were greatly increased, and the anti-interference ability can certainly be improved in principle.
Based on the detailed analysis of the 2-D phase history of SAR as illustrated in Equation (2), we find that random agility of carrier frequency is not recommended due to its compensation requiring a heavy computation load. And this finding has rarely been seen in previous studies in the literature.
Comparing the results of experimental verification of SAR image formation with agile interpulse waveform parameters, as shown in Figure 9f and Figure 13, with the example images of Figure 4, Figure 5 and Figure 6, acquired by the 4.5 kg microminiature SAR sensor with constant waveform parameters, it can be seen that airborne SAR can acquire satisfactory high-quality SAR images definitely under the condition of random agility of interpulse waveform parameters through our proposed processing flowchart shown in Figure 2.
We had focused on the study of high-resolution airborne SAR image formation under the condition of pseudo-random agility of interpulse waveform parameters in this paper. In the future, we will try to carry out practical countermeasures against active coherent jammer-type systems to verify the actual suppression effect of our proposed method.

6. Conclusions

The pseudo-random agility technology of interpulse waveform parameters can certainly improve the anti-active coherent interference ability of airborne SAR. However, this pseudo-random agility technology itself will bring some difficulties and challenges to accurate SAR imaging. The main purpose of this work is to study and verify whether airborne SAR with agile interpulse waveform parameters can acquire high-quality high-resolution SAR images just like conventional airborne SAR with constant waveform parameters. The main findings of this paper are as follows.
(1)
The influences of agile interpulse waveform parameters on SAR imaging were analyzed. It can be seen that pulse width or chirp rate, polarity, initial phase and PRI can be randomly altered within a certain range and also can be easily compensated during the SAR imaging process. But, changing carrier frequency randomly is not recommended due to its compensation requiring a heavy computation load. From the above analysis, an efficient processing flowchart of high-resolution SAR image formation which was appropriate for agile interpulse waveform parameters was proposed. Its computation load was comparable to that of traditional SAR image formation with constant waveform parameters.
(2)
A Ku band 4.5 kg microminiature SAR, which was innovatively designed and highly integrated, was introduced briefly. Several example SAR images with constant interpulse waveform parameters were shown. Considerations to upgrade this SAR with the capability of random agility of interpulse waveform parameters were given.
(3)
SAR raw data with agile interpulse waveform parameters were acquired using this upgraded SAR sensor. Then, real high-resolution SAR images with resolutions of 0.5 m and 0.15 m were processed successfully from the acquired SAR raw data and were shown under the circumstance of randomly changing the transmitted wideband LFM signal waveform pulse parameters one by one.
The above research results fully indicate that the proposed processing scheme of high-resolution SAR image formation for agile interpulse waveform parameters studied in this paper is simple and effective, which lays a solid foundation for future fielding airborne SAR to counter active coherent interference in complex electromagnetic environments by means of active waveform countermeasure, and it has been of great practical value in engineering.

Author Contributions

Conceptualization, Z.Y. and D.Z.; methodology, Z.Y.; software, Z.Y.; validation, Z.Y.; formal analysis, S.N.; investigation, S.N.; data curation, Z.Y.; writing—original draft preparation, Z.Y.; writing—review and editing, Z.Y.; supervision, J.L.; project administration, D.Z.; funding acquisition, D.Z. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by Guangdong Basic and Applied Basic Research Foundation under Grant 2020B1515120060 and in part by Funding for Outstanding Doctoral Dissertation in NUAA under Grant BCXJ23-07.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goj, W.W. Synthetic Aperture Radar and Electronic Warfare; Artech House: Norwood, MA, USA, 1993. [Google Scholar]
  2. Schroer, R. Electronic warfare. IEEE Aerosp. Electron. Syst. Mag. 2003, 18, 49–54. [Google Scholar]
  3. Li, N.J.; Zhang, Y.T. A survey of radar ECM and ECCM. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 1110–1120. [Google Scholar]
  4. Roome, S. Digital radio frequency memory. Electron. Commun. Eng. J. 1990, 2, 147–153. [Google Scholar] [CrossRef]
  5. Liu, G.S.; Gu, H.; Zhu, X.H.; Sun, Y.B.; Sun, G.M.; Su, W.M. A contrast experiment of ECCM of random to pseudo-random binary phase coded CW radar. In Proceedings of the 1996 International Conference on Signal Processing, Beijing, China, 14–18 October 1996; pp. 1621–1624. [Google Scholar]
  6. Soumekh, M. SAR-ECCM using Phase-Perturbed LFM Chirp Signals and DRFM Repeat Jammer Penalization. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 191–205. [Google Scholar] [CrossRef]
  7. Garmatyuk, D.S.; Narayanan, R.M. ECCM Capabilities of an Ultrawideband Bandlimited Random Noise Imaging Radar. IEEE Trans. Aerosp. Electron. Syst. 2002, 38, 1243–1255. [Google Scholar] [CrossRef]
  8. Huang, Y.; Zhao, B.; Tao, M.L.; Chen, Z.Y.; Hong, W. Review of Synthetic Aperture Radar Interference Suppression. J. Radars 2020, 9, 86–106. (In Chinese) [Google Scholar]
  9. Yu, C.R.; Ma, X.L.; Zhang, Y.S.; Dong, Z.; Liang, D.N. Multichannel SAR ECCM based on Fast-time STAP and Pulse diversity. In Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Vancouver, BC, Canada, 24–29 July 2011; pp. 921–924. [Google Scholar]
  10. Li, W.; Liang, D.N.; Dong, Z. SAR Anti-jamming Method Based on Jittering Slope Polarity of LFM and Thresholding. J. Remote Sens. 2007, 11, 171–176. (In Chinese) [Google Scholar]
  11. Feng, Q.Q.; Xu, H.P.; Wu, Z.F.; Sun, B. Deceptive jamming suppression for SAR on time-varying initial phase. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Beijing, China, 10–15 July 2016; pp. 4996–4999. [Google Scholar]
  12. Zhang, Y.; Guan, Z.Q.; Li, W.Y. SAR ECCM and Imaging Based on Non-uniform PRI. In Proceedings of the 2016 International Conference on Wireless Communications WCNM, Chennai, India, 23–25 March 2016; pp. 2292–2294. [Google Scholar]
  13. Cui, G.L.; Fan, T.; Kong, Y.K.; Yu, X.X.; Sha, M.H.; Kong, L.J. Pseudo-random Agility Technology for Interpulse Waveform Parameters in Airborne Radar. J. Radars 2022, 11, 213–226. (In Chinese) [Google Scholar]
  14. Jenq, Y.C. Perfect Reconstruction of Digital Spectrum from Nonuniformly Sampled Signals. IEEE Trans. Instrum. Meas. 1997, 46, 649–652. [Google Scholar] [CrossRef]
  15. Zhu, D.Y.; Mao, X.H.; Li, Y.; Zhu, Z.D. Far-Field Limit of PFA for SAR Moving Target Imaging. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 917–929. [Google Scholar] [CrossRef]
  16. Mao, X.H.; Zhu, D.Y.; Zhu, Z.D. Polar format algorithm wavefront curvature compensation under arbitrary radar flight path. IEEE Geosci. Remote Sens. Lett. 2012, 9, 526–530. [Google Scholar] [CrossRef]
  17. Doerry, A.W. Performance Limits for Synthetic Aperture Radar; Sandia National Laboratories: Albuquerque, NM, USA, 2006. [Google Scholar]
  18. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2004. [Google Scholar]
  19. Franceschetti, G.; Lanari, R. Synthetic Aperture Radar Processing; CRC Press: Carabas, FL, USA, 1999. [Google Scholar]
  20. Carrara, W.G.; Goodman, R.S.; Majewski, R.M. Spotlight Synthetic Aperture Radar Signal Processing Algorithms; Artech House: Norwood, MA, USA, 1995. [Google Scholar]
  21. Raney, R.K.; Runge, H.; Bamler, R.; Cumming, I.G.; Wong, F.H. Precision SAR Processing Using Chirp Scaling. IEEE Trans. Geosci. Remote Sens. 1994, 32, 786–799. [Google Scholar] [CrossRef]
  22. Soumekh, M. Synthetic Aperture Radar Signal Processing with MATLAB Algorithms; Wiley-Interscience: Hoboken, NJ, USA, 1999. [Google Scholar]
  23. Jakowatz, C.V.; Wahl, D.E.; Eichel, P.H.; Ghiglia, D.C.; Thompson, P.A. Spotlight-Mode Synthetic Aperture Radar: A Signal Processing Approach; Kluwer Academic Publishers: Norwell, MA, USA, 1996. [Google Scholar]
  24. Gorham, L.; Rigling, B. Scene Size Limits for Polar Format Algorithm. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 73–84. [Google Scholar] [CrossRef]
  25. Horvath, M.S.; Gorham, L.A.; Rigling, B.D. Scene Size Bounds for PFA Imaging with Postfiltering. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1402–1406. [Google Scholar] [CrossRef]
  26. Gorham, L.A.; Rigling, B.D. Fast Corrections for Polar Format Algorithm With a Curved Flight Path. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 2815–2824. [Google Scholar] [CrossRef]
  27. Garber, W.L.; Hawley, R.W. Extensions to Polar Formatting with Spatially Variant Post-Filtering. In Proceedings of the 2011 Conference on Algorithms for Synthetic Aperture Radar Imagery XVIII SPIE, Orlando, FL, USA, 12–17 February 2011; pp. 1–12. [Google Scholar]
  28. Doren, N.E.; Jakowatz, C.V.; Wahl, D.E.; Thompson, P.A. General formulation for wavefront curvature correction in polar formatted spotlight-mode SAR images using space-variant post-filtering. In Proceedings of the 1997 International Conference on Image Processing, Santa Barbara, CA, USA, 26–29 October 1997; pp. 861–864. [Google Scholar]
  29. Mao, X.H.; Zhu, D.Y.; Zhang, Y.D. Knowledge-Aided Two-Dimensional Autofocus for Synthetic Aperture Radar. In Proceedings of the 2013 IEEE Radar Conference, Ottawa, ON, Canada, 29 April–3 May 2013; pp. 1–6. [Google Scholar]
  30. Mao, X.H.; Zhu, D.Y.; Zhu, Z.D. Autofocus correction of APE and residual RCM in Spotlight SAR polar format imagery. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 2693–2706. [Google Scholar] [CrossRef]
  31. Mao, X.H.; Zhu, D.Y. Two-dimensional Autofocus for Spotlight SAR Polar Format Imagery. IEEE Trans. Comput. Imag. 2016, 2, 524–539. [Google Scholar] [CrossRef]
  32. Zhu, D.Y.; Ye, S.H.; Zhu, Z.D. Polar Format Algorithm using Chirp Scaling for Spotlight SAR Image Formation. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 1433–1448. [Google Scholar] [CrossRef]
  33. Wahl, D.E.; Eichel, P.H.; Ghiglia, D.C.; Jakowatz, C.V. Phase Gradient Autofocus—A Robust Tool for High Resolution SAR Phase Correction. IEEE Trans. Aerosp. Electron. Syst. 1994, 30, 827–834. [Google Scholar] [CrossRef]
Figure 1. Spotlight SAR data collection geometry.
Figure 1. Spotlight SAR data collection geometry.
Remotesensing 16 00164 g001
Figure 2. Flowchart of airborne SAR image formation with agile interpulse waveform parameters.
Figure 2. Flowchart of airborne SAR image formation with agile interpulse waveform parameters.
Remotesensing 16 00164 g002
Figure 3. 4.5 kg microminiature SAR: (a) physical photo; (b) front view; (c) rear view; (d) assembly view.
Figure 3. 4.5 kg microminiature SAR: (a) physical photo; (b) front view; (c) rear view; (d) assembly view.
Remotesensing 16 00164 g003
Figure 4. Dual-polarization stripmap SAR image, 0.5 m resolution: (a) VV polarization; (b) VH polarization.
Figure 4. Dual-polarization stripmap SAR image, 0.5 m resolution: (a) VV polarization; (b) VH polarization.
Remotesensing 16 00164 g004
Figure 5. Real-time stripmap image of rural area at 0.3 m resolution.
Figure 5. Real-time stripmap image of rural area at 0.3 m resolution.
Remotesensing 16 00164 g005
Figure 6. Real-time spotlight images of urban area at 0.15 m resolution.
Figure 6. Real-time spotlight images of urban area at 0.15 m resolution.
Remotesensing 16 00164 g006
Figure 7. Agile interpulse waveform parameters at resolution 0.5 m: (a) pulse width; (b) PRI; (c) initial phase; (d) chirp rate polarity.
Figure 7. Agile interpulse waveform parameters at resolution 0.5 m: (a) pulse width; (b) PRI; (c) initial phase; (d) chirp rate polarity.
Remotesensing 16 00164 g007
Figure 8. Amplitude equalization factor for compensating for the fluctuations in processing gain and average power resulting from the random agility of pulse width and PRI at a resolution of 0.5 m.
Figure 8. Amplitude equalization factor for compensating for the fluctuations in processing gain and average power resulting from the random agility of pulse width and PRI at a resolution of 0.5 m.
Remotesensing 16 00164 g008
Figure 9. Rural area SAR images at 0.5 m resolution with agile interpulse waveform parameters: (a) no compensation for agility of pulse width or chirp rate; (b) no compensation for agility of chirp rate polarity; (c) no compensation for agility of initial phase; (d) no compensation for agility of PRI; (e) no compensation for all agile interpulse waveform parameters; (f) compensation for all agile interpulse waveform parameters according to the flowchart shown in Figure 2.
Figure 9. Rural area SAR images at 0.5 m resolution with agile interpulse waveform parameters: (a) no compensation for agility of pulse width or chirp rate; (b) no compensation for agility of chirp rate polarity; (c) no compensation for agility of initial phase; (d) no compensation for agility of PRI; (e) no compensation for all agile interpulse waveform parameters; (f) compensation for all agile interpulse waveform parameters according to the flowchart shown in Figure 2.
Remotesensing 16 00164 g009
Figure 10. The comparison of the chirp rate polarity of the jamming signal matches and mismatches the transmitted signal: (a) result of pulse compression and rotating-platform processing with no compensation of chirp rate polarity agility; (b) local enlarged result of the area marked by yellow dotted line box of (a); (c) result of pulse compression and rotating-platform processing with compensation of chirp rate polarity agility; (d) local enlarged result of the area marked by yellow dotted line box of (c).
Figure 10. The comparison of the chirp rate polarity of the jamming signal matches and mismatches the transmitted signal: (a) result of pulse compression and rotating-platform processing with no compensation of chirp rate polarity agility; (b) local enlarged result of the area marked by yellow dotted line box of (a); (c) result of pulse compression and rotating-platform processing with compensation of chirp rate polarity agility; (d) local enlarged result of the area marked by yellow dotted line box of (c).
Remotesensing 16 00164 g010
Figure 11. Agile interpulse waveform parameters at resolution 0.15 m: (a) pulse width; (b) PRI; (c) initial phase; (d) chirp rate polarity.
Figure 11. Agile interpulse waveform parameters at resolution 0.15 m: (a) pulse width; (b) PRI; (c) initial phase; (d) chirp rate polarity.
Remotesensing 16 00164 g011
Figure 12. Amplitude equalization factor for compensating the fluctuations of processing gain and average power resulted from the random agility of pulse width and PRI at a resolution of 0.15 m.
Figure 12. Amplitude equalization factor for compensating the fluctuations of processing gain and average power resulted from the random agility of pulse width and PRI at a resolution of 0.15 m.
Remotesensing 16 00164 g012
Figure 13. Spotlight SAR image of a power transmission line scene at 0.15 m resolution with agile interpulse waveform parameters.
Figure 13. Spotlight SAR image of a power transmission line scene at 0.15 m resolution with agile interpulse waveform parameters.
Remotesensing 16 00164 g013
Table 1. 4.5 kg microminiature SAR specifications.
Table 1. 4.5 kg microminiature SAR specifications.
ParameterValueComments
Weight4.3 kgmeasured value
Size (length × width × height)324.7 mm × 98 mm × 164.5 mmmeasured value
Power+28 VDC at 125 W maxmeasured value
Frequency BandKu
RF Bandwidth1500 MHz max
Co-aperture dual-polarized patch array AntennaAperture size250 mm × 90 mmMeasured value at center frequency.
Gain21.5 dB
Beamwidth8.1° AZ/19.6° EL
Sidelobe−20.5 dB
Isolation−35 dB
Gimbal TypeAzimuth and Elevation Dual AxisAzimuth coverage: −90° ± 30° or +90° ± 30°, Elevation coverage: 0~180°.
Transmitter/Receiver with Exciter moduleTransmitter Power80 W peakSolid state power amplifier
Noise Figure of RF Receiver4 dBIncludes built in limiter and dual channel receiver.
RF losses~3 dBTransmitter output to LNA/Limiter input.
Processing/data storage moduleProcessing board1 FPGA + 2 multi-core DSPOn-board real-time image formation
Data storage board2 TB SSDStore dual-channel raw data and real-time SAR images
Power conversion moduleOutput all kinds of low-voltage power supply required by each module and cooling fanThe input is 28 VDC.
Motion sensorMEMS based IMU/GPSCOTS (commercial off-the-shelf)
SAR Imaging modesSpotlight and stripmap
ResolutionSpotlight0.15 m, 0.3 m, 0.5 m optional
Stripmap0.3 m, 0.5 m, 1 m, 3 m optional
Maximum Range≮15 kmAssume 3 m resolution, −25 dB noise equivalent reflectivity
Image Size8 K × 8 K pixels maxBoth along-track and cross-track for spotlight mode
SAR Imaging formation on-boardCS-PFA + SVPF + PGASpace-Variant Post-Filtering [16,28] based Polar Format Algorithm using Chirp Scaling [32] with Phase Gradient Autofocus [33]
Table 2. Radar original waveform parameters.
Table 2. Radar original waveform parameters.
Waveform ParametersValue
Resolution mode/m0.5/0.15
Bandwidth/MHz400/1200
Chirp rate polarity+
Pulse width/ μ s 10
PRI/ μ s 1000
Table 3. The description of four agile parameters.
Table 3. The description of four agile parameters.
Agile ParametersVariation RangeStep Width
Pulse width/ μ s 9~1180 ns
PRI/ μ s 900~110080 ns
Initial phase/°0~360°360°/214
Chirp rate polarity−1 or +1
Table 4. Contrast values corresponding to Figure 9a–d and the same area indicated by the white dotted line box in Figure 9f.
Table 4. Contrast values corresponding to Figure 9a–d and the same area indicated by the white dotted line box in Figure 9f.
Figure 9aFigure 9bFigure 9cFigure 9dFigure 9f
Contrast value1.00510.6025.04514.94630.859
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ye, Z.; Zhu, D.; Niu, S.; Lv, J. Studies on High-Resolution Airborne Synthetic Aperture Radar Image Formation with Pseudo-Random Agility of Interpulse Waveform Parameters. Remote Sens. 2024, 16, 164. https://doi.org/10.3390/rs16010164

AMA Style

Ye Z, Zhu D, Niu S, Lv J. Studies on High-Resolution Airborne Synthetic Aperture Radar Image Formation with Pseudo-Random Agility of Interpulse Waveform Parameters. Remote Sensing. 2024; 16(1):164. https://doi.org/10.3390/rs16010164

Chicago/Turabian Style

Ye, Zheng, Daiyin Zhu, Shilin Niu, and Jiming Lv. 2024. "Studies on High-Resolution Airborne Synthetic Aperture Radar Image Formation with Pseudo-Random Agility of Interpulse Waveform Parameters" Remote Sensing 16, no. 1: 164. https://doi.org/10.3390/rs16010164

APA Style

Ye, Z., Zhu, D., Niu, S., & Lv, J. (2024). Studies on High-Resolution Airborne Synthetic Aperture Radar Image Formation with Pseudo-Random Agility of Interpulse Waveform Parameters. Remote Sensing, 16(1), 164. https://doi.org/10.3390/rs16010164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop