Next Article in Journal
Optical Temperature Sensor Based on Polysilicon Waveguides
Next Article in Special Issue
Low-Altitude Windshear Wind Speed Estimation Method Based on KASPICE-STAP
Previous Article in Journal
Salient Targets and Fear of Falling Changed the Gait Pattern and Joint Kinematic of Older Adults
Previous Article in Special Issue
Airborne Radar Anti-Jamming Waveform Design Based on Deep Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method for Interferometric Phase Estimation in Dual-Channel Cancellation

Nanjing Research Institute of Electronic Technology, Nanjing 210039, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(23), 9356; https://doi.org/10.3390/s22239356
Submission received: 26 September 2022 / Revised: 9 November 2022 / Accepted: 28 November 2022 / Published: 1 December 2022

Abstract

:
Multichannel SAR systems have grown rapidly over the past decade due to their powerful high-resolution and wide-swath (HRWS) capabilities. Because spatially separated channels also have the potential to suppress jamming, dual-channel cancellation is a general method that is effective regardless of the type of jamming signal. In this paper, the principle of dual-channel cancellation (DCC) is introduced, and several practical problems using DCC are also discussed. Moreover, this paper emphasizes interferometric phase estimation, which is the key to DCC. If the jamming-to-signal ratio (JSR) is high, the interferometric phase can be estimated accurately from the interferometry of two channel signals, but estimation becomes rather difficult when the JSR decreases. To solve the problem of interferometric phase estimation under a low JSR, a novel interferometric phase estimation method using cosine similarity is proposed in this paper. L-band airborne dual-channel SAR is performed to investigate the applicability of the method. The results not only prove that cosine similarity is an effective method for interferometric phase estimation, but also demonstrate the potential of DCC in the SAR anti-jamming processing.

1. Introduction

Multichannel synthetic aperture radar (SAR) allows continuous observation of the Earth’s surface regardless of weather and daylight, which has motivated the rapid development of SAR technology. Among the different methods, ground moving target indication is used to monitor the velocities of moving targets [1,2,3,4,5,6], while interferometric SAR (InSAR) generates digital elevation maps of the ground surface [7,8,9].
SAR enjoys great benefits in the face of electronic counter measures due to its coherent processing and wide bandwidth. However, if the jamming power is high enough, SAR performance can still be seriously affected. In most instances, the jamming signal is similar and may even be identical to the radar signal. It is difficult to distinguish the interfering signal from the radar echo because they overlap in both the time and frequency domains [1,10,11].
Several anti-jamming methods have been proposed, such as the RFI suppression method [12], WBI suppression method [13], STAP method [14], spatial location feature recognition method [15], and other methods [16,17,18,19]. Evaluation methods for the anti-jamming performance of SAR have also been studied [20,21]. Taking full advantage of multichannel approaches, several researchers have proposed multichannel SAR methods [14,22,23,24,25,26,27,28,29], but they are usually limited to specific jamming types. To deal with this problem, Ma [30] proposed a general method, dual-channel cancellation, which is effective regardless of the jamming type when dealing with an active jammer. Theoretical analysis and simulations have demonstrated the effectiveness of DCC, but problems have arisen in practice. Ma pointed out that without an accurate estimate of the interferometric phase of the radiation source (IPSR), the jamming signal will not be completely cancelled [30].
In this paper, other practical problems that may be faced using DCC are discussed, such as the estimation of the IPSR, the effect of co-registration of two images on DCC, and resolution degradation after DCC. Among these problems, the most intractable one is the estimation of IPSR.
Although DCC has been proven to be effective in simulations [28,30,31], no flight experiment has ever been reported. To validate the DCC method and demonstrate the above issues, an experiment was conducted using an L-band airborne dual-channel SAR sensor. We used the DCC method and compared the SAR images before and after DCC, and the results of the experiment were quite satisfactory.
The remainder of this paper is organized as follows. In Section 2, the signal model of the dual channel cancellation is introduced. Combined with simulation, the estimation of IPSR based on cosine similarity is then studied. After that, the effects of co-registration, resolution degradation, and dark straps are analyzed. In Section 3, the airborne experimental results at the L-band are shown to prove the validity and practicality of our method. Finally, we discuss our findings and draw conclusions in Section 4 and Section 5.

2. Materials and Methods

For dual-channel SAR, the complex images that are processed are highly correlated, regardless of whether the two channels are along-track or cross-track. Interferometric images of the two channels contain certain information about the height and velocity of the imaging area, and the two intensity images are almost identical if the amplitude-phase errors are neglected or compensated for. On the other hand, the jamming signals that the two channels receive have a small time lag and a corresponding phase shift, which change slowly according to the geometrical relationship between the SAR channels and the source of radiation. If the IPSR is estimated precisely, then the jamming signals can be eliminated by phase compensation and cancellation.

2.1. Signal Model of the Dual Channel Cancellation

From the study in [32], for an arbitrary point target, we assume that the echoes of channel 1 and channel 2 received during the slow time η are written as follows:
s pt 1 ( t ) = A pt ( t ) exp [ j ϕ pt ( t ) ]
s pt 2 ( t ) = A pt ( t Δ t pt ) exp [ j ϕ pt ( t ) + j Δ ϕ pt ( η ) ]
where t is quick time, A pt ( t ) and ϕ pt ( t ) are the amplitude and phase of the echo that channel 1 received, Δ t pt is the time lag of the echoes between the two channels, and Δ ϕ pt ( η ) is the phase shift of the echoes between the two channels.
Let a source of radiation be deployed in the imaging area, transmitting signals of the same frequency band toward the dual-channel SAR. The jamming signals received by channel 1 and channel 2 can be written as
s j 1 ( t ) = A j ( t ) exp [ j ϕ j ( t ) ]
s j 2 ( t ) = A j ( t Δ t j ) exp [ j ϕ j ( t ) + j Δ ϕ j ( η ) ]
where A j ( t ) and ϕ j ( t ) are the amplitude and phase of the jamming signals, respectively, and Δ ϕ j ( t ) is the IPSR. Note that Δ ϕ j ( η ) is the slow time variant. Usually, if the two channels are not too far away, the time lag Δ t j is neglected, and A j ( t Δ t j ) in (4) becomes A j ( t ) . The signals that the two channels receive are the summation of the ground echoes and the jamming signals, which are
s 1 ( t ) = A pt ( t ) exp [ j ϕ pt ( t ) ] + A j ( t ) exp [ j ϕ j ( t ) ]
s 2 ( t ) = A pt ( t ) exp [ j ϕ pt ( t ) + j Δ ϕ pt ( η ) ] + A j ( t ) exp [ j ϕ j ( t ) + j Δ ϕ j ( η ) ]
Compensate (5) with exp [ Δ ϕ j ( η ) ] , and form the difference with (6). This results in the following difference signal:
s d ( t ) = s 2 ( t ) s 1 ( t ) exp [ j Δ ϕ j ( η ) ]   = A pt ( t ) exp [ j ϕ pt ( t ) + j Δ ϕ pt ( η ) ] A pt ( t ) exp [ j ϕ pt ( t ) + j Δ ϕ j ( η ) ]
From (7), we can see that the jamming signals are eliminated. Meanwhile, the phase of the echoes is affected by Δ ϕ j ( η ) , wherein if Δ ϕ pt ( η ) is close to Δ ϕ j ( η ) , and s d ( t ) is close to 0, it means that some of the echoes would also be cancelled.

2.2. Estimation of the Interferometric Phase of the Source of Radiation

The most fundamental step of DCC is the estimation of the IPSR, which is denoted as Δ ϕ j ( η ) in the equations above.
The interferometric phase of the dual channels (IPDC) is the phase difference of the vector addition of the ground echoes and the jamming signals.
On the one hand, if the jamming-to-signal ratio (JSR) is high enough, the jamming signals dominate the echoes, and IPSR dominates the interferometric phase, which means that IPSR is very close to IPDC. Since IPDC is easy to compute, in this case, IPSR is considered equal to IPDC [30]. The cost of this approximation is that IPSR is not precise and will cause residual interference after DCC.
The jamming residuals after DCC with different IPSR errors are obtained by simulation, as shown in Figure 1. From Figure 1, we can see that the IPSR error decreases with the reduction of the jamming residual after DCC. Generally, cancellation with residual interference less than −10 dB can be considered effective. In such a case, the IPSR error should not be greater than 18°, which corresponds to the blue curve in Figure 1.
On the other hand, if the JSR is not high enough, the estimation of IPSR will become rather intractable. The difference between IPSR and IPDC shows relatively strong ground echoes.
Figure 2 and Figure 3 show the IPDC under different JSRs, in which the color map represents the interferometric phase ranging from 0 to 2π. It is clear that a higher JSR can result in a smoother IPDC. If one examines the arbitrary azimuth bin in Figure 2, it can be found that the IPDCs along the range bins are nearly the same, which is consistent with the previous analysis. In addition, the distribution of IPDCs along the azimuth direction can be easily identified as being nearly linear.
However, the performance of IPDC will be seriously affected by the value of the JSR. As is clearly shown in Figure 3, when the JSR is reduced to 0 dB, the IPDC image will behave more noisily than in Figure 2, in which the IPDC distribution pattern along the azimuth direction can barely be identified.
Figure 4 and Figure 5 show the IPSR and IPDC distribution comparison along the azimuth direction under different JSRs.
As is shown in Figure 4, when the JSR is 10 dB, there is a good match between the IPSR and IPDC, in which we can note that the IPDC is wrapped from 0° to 360°. However, in Figure 5, where the JSR is 0 dB, the IPSRs are quite different from the IPDCs and will be clearly not effective in DCC.
To address this problem, a traversal algorithm using cosine similarity is presented. For an arbitrary pulse repeat time (PRT), we traverse the IPSR from 0 to 2π and perform DCC. As the interference signals are different from the SAR echoes in certain aspects (for example, in the frequency domain), then we can use the similarity of the spectrum of pure SAR echoes and that of DCC signals to determine the optimal IPSR estimation.
In this method, cosine similarity is used to describe the similarity in the spectrum shapes, which is calculated as below:
γ ( I P S R ) = k = 1 n F 1 k F 2 k ( I P S R ) k = 1 n F 1 k 2 k = 1 n F 2 k 2 ( I P S R )
where F 1 k demotes the spectrum of the pure SAR echoes, F 2 k ( I P S R ) denotes the spectrum of the DCC signal corresponding to IPSR, and γ ( I P S R ) is the cosine similarity corresponding to IPSR. When the IPSR meets the optimal value, the cosine similarity will reach its peak.
Since the IPSR changes slowly along the azimuth direction, one does not need to perform the traversal algorithm for each PRT. To simplify the calculation, one can repeat this process along the azimuth direction at a suitable PRT interval and use interpolation to achieve IPSR of each PRT.
To verify this method, a preliminary simulation based on real data is conducted. The data were obtained using an L-band airborne SAR, whose main parameters are listed in Table 1. The original data were achieved without jamming, and the spectrum is shown in the upper side of Figure 6. In the simulation, a jamming signal with 30 MHz bandwidth is added to the original data, in which the JSR is 0 dB. The shape of the spectrum is then distorted, as shown on the lower side of Figure 6. After the jamming operation, we traverse the IPSR from 0 to 2π, perform DCC, and examine the shape of the spectrum until a close shape is matched.
In this simulation, the ideal IPSR is set to 90°, and the IPSR interval is set to 1°. Under this situation, the cosine similarity of different IPSR values could be obtained, as clearly shown in Figure 7. As we can see from the Figure 7, the cosine similarity will reach the peak value of 0.8414 when the IPSR is 90°, which is consistent with the simulation condition.
To verify the validity of the cosine similarity method in different situations, we perform IPSR estimation simulations under different JSR, as shown in Figure 8. In Figure 8, it is demonstrated that the IPSR estimation results are quite accurate under different IPSR angles. Moreover, we can find that the cosine similarity method can also achieve good performance, even at −10 dB and −20 dB for different JSRs.

2.3. Effects of Co-Registration

In Section 2.1, the time lag between two SAR channels is neglected, which is true for the along-track DCC, since the distance between the two channels is usually less than 1 m, and such an along-track distance would cause a time lag that is much less than that of a sample interval. However, for cross-track DCC, the channel distance is usually in meters for airborne platforms and can even reach hundreds of meters for spaceborne platforms. The two channels are distributed in the cross-track direction, which makes it possible to cause a time lag of tens to thousands of sample intervals. Thus, range co-registration for interference signals will be needed for cross-track DCC.
Assuming that IPSR estimation is precise, the jamming residual after DCC suffers from a co-registration error. The residual is dependent on the jamming bandwidth to signal bandwidth ratio since the narrowband signal phase changes slower than that of the wideband signal.
Figure 9 shows the jamming residual with the co-registration error of different jamming bandwidth to signal bandwidth ratios. The curves show that the jamming residual suffers from the jamming bandwidth to signal bandwidth ratio and co-registration error, which will grow with either the jamming bandwidth to signal bandwidth ratio or the co-registration error. As we can see from Figure 9, when the jamming bandwidth is equal to the signal bandwidth, a 0.2 resolution cell co-registration error will approximately cause −12 dB jamming residual, which is commonly considered to be acceptable.

2.4. Ghost Images

Along-track DCC usually leads to ghost images. There is a co-registration step in along-track multichannel SAR processing. However, since the jamming signals must be processed in the same PRT, there will be no azimuth co-registration step, which causes ghost images. The degree of the ghost image is proportional to the channel distance since the channel distance of a multichannel SAR is typically close to its azimuth resolution. Figure 10 gives a simulation demonstration that shows how, if the channel distance is much too large, then clear ghost images would appear.

2.5. Dark Strap

There is usually more than one dark strap on DCC images. The direction of the dark straps mainly depends on the channel configuration. In cross-track DCC, the interferometric phase map is periodic along the range direction. Those pixels whose interferometric phase is close to IPSR are cancelled, forming dark straps along the azimuth direction and cycle along the range direction. In along-track DCC, however, the dark straps are located along the range direction and cycle along the azimuth direction.
The dark strap is the negative effect of DCC, which will slightly affect the interpretation of SAR images. The period of the dark straps is inversely proportional to the channel distance [22]. Ma proposed a three-channel cancellation method that can reduce the dark straps. Moreover, Ma also generalizes the method to N-channel cancellation [31]. Another method is inverse power compensation [31], in which the amplitudes are compensated inverse to the powers, which will also mitigate the effect.

2.6. Complexity Analysis

To demonstrate the innovation more comprehensively, we discuss the complexities of our method and the classical Chirp Scaling (CS) algorithm and assume that the input data are D R N × K . For fast Fourier transform (FFT) computing, we always extend the size to the nearest power of two: D R N 2 × K 2 , N 2 = 2 m N , K 2 = 2 n K . The time cost of 128 FFT points and 128 dot product points on DSP are 4 μs and 2 μs, respectively. The major time cost of the CS algorithm is summarized in Table 1. From the table, we can obtain the total time cost of CS algorithm as
T C S = 3 K N 2 + 8 K 2 N 2 256
Similarly, we can summarize the major time cost of our method in Table 2. Additionally, the total time cost of our method is
T o u r = 3 N K 2 256
Comparing T C S and T o u r , it is well-known that the complexity of our method is less than that of the CS algorithm. Our method just determines the FFT and the dot products once, and CS algorithm conducts them multiple times.

3. Experimental Studies

Nriet-SAR (SAR of Nanjing Research Institute of Electronic Technology), an airborne dual-band (L, X) multichannel polarimetric SAR sensor, was applied to conduct the experiment [33]. The principal parameters are shown in Table 3.
The geometric relationship between the SAR sensor and the jammer is briefly illustrated in Figure 11, in which the radar moves along a linear flight path with a typical flight height of 7000 m, and the jammer is in the SAR imaging area.
Figure 12 shows the optical image (from Google Earth) of the experimental area in Weinan City, Shanxi, China. The blue icon in the middle of Figure 12 indicates the location of the jammer (source of radiation) used in the experiment, and the blue arrow shows the antenna orientation of the transmitting and receiving antennas. Additionally, the flight direction of the Nriet-SAR platform is from south to north (from right to left in Figure 12).
Figure 13 shows the experimental placement of the jammer, the transmitting antenna, and the receiving antenna. To clearly explain the working process of the jammer in detail, the schematic diagram and processing chain are shown in Figure 14. After the receiving antenna of the jammer receives the chirp signal transmitted from the radar, the jammer is triggered to generate the jamming signals. The jamming signals are transmitted by the transmitting antenna to the radar after being amplified by the power amplifier.
To verify the effectiveness of our proposed method, two linear SAR trajectories have been designed to observe the imaging area shown in Figure 12 in this experiment, in which the parameter configuration of SAR is shown in Table 1. In the first linear SAR trajectory, the jammer is turned off, which can be seen as a set of comparative experiments, so we can achieve a perfect SAR image of the imaging area with commonly used SAR imaging algorithms. In the other linear SAR trajectory, the jammer is turned on and transmits jamming signals to the receiving antenna of SAR according to the processing chain shown in Figure 14.
In this airborne experiment, the original SAR data obtained by Nriet-SAR are processed using along-track DCC. Figure 15 shows the SAR imaging result without jamming, and Figure 16 shows the SAR image with jamming. As we can see from Figure 15 and Figure 16, when the jammer is turned off, SAR can perfectly reconstruct the SAR image of the ground features in the imaging area. However, it can be clearly found that parts of the ground features are not clearly visible and that they are seriously influenced by the jamming signals in Figure 16.
To eliminate the influence of the jamming signals and to reconstruct the SAR image of the imaging area, a cosine similarity-based DCC method has been applied in the airborne experimental data. Figure 17 and Figure 18 show the DCC images, while the IPSR error of Figure 18 is set to 18° on purpose, and the IPSR of Figure 17 is relatively precise.
As we can see from Figure 17, most of the interference images have been eliminated, and the ground features have become visible after accurate cosine similarity-based along-track DCC. Due to the set of the deliberate IPSR errors of 18°, which can be considered as the wrong estimations of the IPSR, Figure 18 suffers from the slight influence of the jamming signals. However, a dark strap still appears in the middle of Figure 17 and Figure 18 along the range direction of the jammer (from top to bottom), which is consistent with our analysis in Section 2.5.
Compared with the SAR image without jamming in Figure 15, the ground features in the imaging area can still be well reconstructed with our proposed method, while the SAR is severely influenced by the jammer.
To demonstrate the effectiveness of our method, the information entropy of image was used to reflect the potential of suppressing SAR jamming. Image entropy is a statistic based on information theory. The smaller the image entropy is, the more orderly the image is and the more information it contains. First, for a SAR image, we can calculate its entropy by
H = i = 0 255 p i log 2 ( p i )
where p i is the probability of different grey values in the image. It was easy to calculate the entropy of the clear image, the jamming image, the image after DCC, and the image after DCC with an IPSR error of 18°: 4.8648, 5.1113, 4.9347, and 5.0248, respectively. The entropy of the image DCC is less than that of the jamming image, and it is greater than the entropy of the original image. Jamming would cause an increase in image entropy, and DCC could decrease the entropy of the jamming image. Overall, the results indicate that our method can suppress SAR jamming.

4. Discussion

In the practical application of along-track multichannel SAR, the channel error between channels mainly arises from the SAR system error, the baseline error, and the squint angle error. Related research illustrates that the channel error can influence the performance of channel cancellation operations. Specifically, multichannel SAR working at a higher frequency could be more sensitive to channel errors. In this situation, channel errors should be compensated for precisely before the multichannel SAR channel cancellation operation is performed.
This paper mainly focuses on the research on dual-channel SAR systems. As for N-channel SAR systems, the problem of dark straps can be perfectly solved after channel cancellation operations. In the N-channel SAR system, one channel can be chosen as the counter reference. Therefore, the channel cancellation operation of the N-channel SAR system can be regarded as a concatenation of multiple dual-channel SAR systems in which the proposed estimation method of IPSR based on cosine similarity can still be applied effectively.
In our future work, we will focus on the combination of channel error compensation and DCC operation for dual-channel SAR system. Moreover, the channel cancellation of the N-channel SAR system with much lower JSR will also represent the main direction of our future research.

5. Conclusions

In this paper, we study a dual-channel cancellation method for SAR jamming suppression that uses the interferometric phase of the jammer to eliminate the jamming signals via phase compensation and cancellation. To solve the problem of interferometric phase estimation under a low JSR, a novel method for IPSR estimation using cosine similarity is proposed, and several practical problems are discussed. An airborne experiment is conducted, and the preliminary results verify the effectiveness and practicality of the proposed method and show the great potential of the DCC method in SAR interference suppression.

Author Contributions

All the authors have made contributions to the article in different areas. Conceptualization, L.H. and A.L.; methodology, L.H.; software, L.H.; validation, L.H. and H.X.; investigation, L.H.; data curation, L.H. and Z.H.; writing—original draft preparation, L.H. and D.H.; writing—review and editing, L.H., A.L. and D.H.; funding acquisition, Z.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Jiangsu Province under grants BK20190143 and JSGG20190823094603691, and by the National Natural Science Foundation of China under grant 62001436.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available upon request from the corresponding author. The data are not publicly available due to privacy restrictions.

Acknowledgments

The authors would like to thank the staff of the Nanjing Research Institute of Electronic Technology for their valuable conversations and comments.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

AbbreviationFull Name
HRWSHigh-resolution and wide-swath
DCCDual-channel cancellation
JSRJamming-to-signal-ratio
SARSynthetic aperture radar
InSARInterferometric SAR
IPSRInterferometric phase of the source of radiation
IPDCInterferometric phase of the dual channels
PRTPulse repeat time
CSChirp Scaling
FFTFast Fourier transform
IFFTInverse fast Fourier transform
PRFPulse repeat frequency
NESZNoise equivalent sigma zero

References

  1. Huang, L.; Dong, C.; Shen, Z.; Zhao, G. The Influence of Rebound Jamming on SAR GMTI. IEEE Geosci. Remote Sens. Lett. 2015, 12, 399–403. [Google Scholar] [CrossRef]
  2. Huang, Z.; Liu, A.; Xia, X.; Xu, H.; Wu, G. Azimuth Relocation for Multichannel SAR Ground Moving Targets via Noncoregistrated Inteferometry. IEEE Geosci. Remote Sens. Lett. 2020, 18, 652–656. [Google Scholar] [CrossRef]
  3. Tan, W.; Xu, W.; Huang, P.; Huang, Z.; Qi, Y.; Han, K. Investigation of Azimuth Multichannel Reconstruction for Moving Targets in High Resolution Wide Swath SAR. Sensors 2017, 17, 1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wang, C.; Liao, G.; Zhang, Q. First Spaceborne SAR-GMTI Experimental Results for the Chinese Gaofen-3 Dual-Channel SAR Sensor. Sensors 2017, 17, 2683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wang, L.; Li, Y.; Wang, W.; An, D. Moving Target Indication for Dual-Channel Circular SAR/GMTI Systems. Sensors 2020, 20, 158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Xu, J.; Huang, Z.; Yan, L.; Zhou, X.; Zhang, F.; Long, T. SAR Ground Moving Target Indication Based on Relative Residue of DPCA Processing. Sensors 2016, 16, 1676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Giese, C.; Ulrich, D. The TanDEM-X space segment. In Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, 24–29 July 2011. [Google Scholar] [CrossRef]
  8. Lachaise, M.; Bachmann, M.; Schweisshelm, B.; Fritz, T. The Tandem-X Change Dem: Status of the Change Raw Dems Production. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021. [Google Scholar] [CrossRef]
  9. Xu, H.; Li, S.; You, Y.; Liu, A.; Liu, W. Unwrapped Phase Estimation via Normalized Probability Density Function for Multibaseline InSAR. IEEE Access 2018, 7, 4979–4988. [Google Scholar] [CrossRef]
  10. Tian, T.; Zhou, F.; Zhao, B. Multi-receiver deception jamming against synthetic aperture radar. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016. [Google Scholar] [CrossRef]
  11. Liu, Y.; Bu, Y.; Liu, H.; Wang, R. SAR Seeker Deception Jamming and Anti-jamming Measures. In Proceedings of the 2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Xiamen, China, 26–29 November 2019. [Google Scholar] [CrossRef]
  12. Liu, H.; Li, D. RFI Suppression Based on Sparse Frequency Estimation for SAR Imaging. IEEE Geosci. Remote Sens. Lett. 2016, 13, 63–67. [Google Scholar] [CrossRef]
  13. Yang, Z.; Du, W.; Liu, Z.; Liao, G. WBI Suppression for SAR Using Iterative Adaptive Method. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 1008–1014. [Google Scholar] [CrossRef]
  14. Yu, C.; Ma, X.; Zhang, Y.; Dong, Z.; Liang, D. Multichannel SAR ECCM based on Fast-time STAP and Pulse diversity. In Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, 24–29 July 2011. [Google Scholar]
  15. Liu, Z.; Zhang, Q.; Li, K. An Anti-Jamming Method against Two-Dimensional Deception Jamming by Spatial Location Feature Recognition. Sensors 2021, 21, 7702. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, Y.; Zhang, Q.; Liu, Z.; Li, G.; Xiong, S.; Luo, Y. An Anti-Jamming Method against Interrupted Sampling Repeater Jamming Based on Compressed Sensing. Sensors 2022, 22, 2239. [Google Scholar] [CrossRef] [PubMed]
  17. Qiu, X.; Wang, P.; Jiang, J.; Li, X.; Tang, Z.; Wang, S. Research on SAR Anti-jamming Technique Based on Orthogonal LFM-PC Signals with Adaptive Initial Phase. In Proceedings of the 2021 2nd China International SAR Symposium (CISS), Shanghai, China, 3–5 November 2021. [Google Scholar] [CrossRef]
  18. Qiu, X.; Zhang, T.; Li, S.; Yuan, T.; Wang, F. SAR Anti-Jamming Technique Using Orthogonal LFM-PC Hybrid Modulated Signal. In Proceedings of the 2018 China International SAR Symposium (CISS), Shanghai, China, 10–12 October 2018. [Google Scholar] [CrossRef]
  19. Xin, J.; Bin, Y.C.; Dong, Q.W.; Hai, G.Z. Anti-jamming property of converse beam cross sliding spotlight SAR. In Proceedings of the 2009 2nd Asian-Pacific Conference on Synthetic Aperture Radar, Xi’an, China, 26–30 October 2009. [Google Scholar] [CrossRef]
  20. Zhao, J.; Zhang, Q.; Zhang, X.; He, H.; Zhou, T.; Zheng, J. Evaluation of the anti-jamming performance of SAR based on ELECTRE-III. In Proceedings of the 2021 IEEE International Conference on Advances in Electrical Engineering and Computer Applications (AEECA), Dalian, China, 27–28 August 2021. [Google Scholar] [CrossRef]
  21. Zhou, M.; Wang, Q.; He, F.; Meng, J. Impacts of Phase Noise on the Anti-Jamming Performance of Power Inversion Algorithm. Sensors 2022, 22, 2362. [Google Scholar] [CrossRef] [PubMed]
  22. Qin, J.; Yang, J.; He, Z.; Cai, C. Analysis of Target Loss Due to Suppressing SAR Jamming Using Dual-channel Cancellation. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006. [Google Scholar] [CrossRef]
  23. Wang, W.; Wu, J.; Pei, J.; Sun, Z.; Yang, J. Anti-Deceptive Jamming of Jammer on the Coast for Multistatic SAR. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021. [Google Scholar] [CrossRef]
  24. Ji, P.; Xing, S. An Anti-Jamming Method Against SAR Stationary Deceptive Targets Based on DPCA Processing. In Proceedings of the 2019 IEEE 4th International Conference on Signal and Image Processing (ICSIP), Wuxi, China, 19–21 July 2019. [Google Scholar] [CrossRef]
  25. Kong, Z.; Li, P.; Yan, X.; Hao, X. Anti-Sweep Jamming Design and Implementation Using Multi-Channel Harmonic Timing Sequence Detection for Short-Range FMCW Proximity Sensors. Sensors 2017, 17, 2042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wang, R.; Chen, J.; Wang, X.; Sun, B. High-Performance Anti-Retransmission Deception Jamming Utilizing Range Direction Multiple Input and Multiple Output (MIMO) Synthetic Aperture Radar (SAR). Sensors 2017, 17, 123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Wang, R.; Sun, B.; Yi, C.; Chen, J.; Zhou, Y. Multi-Channnel and Mimo Sar Anti-Jamming Analysis. In Proceedings of the 2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018. [Google Scholar] [CrossRef]
  28. Sun, B.; Li, J. A New Interference Elimination Method for Multi-Satellite SAR System. In Proceedings of the 2008 IEEE International Geoscience and Remote Sensing Symposium, Boston, MA, USA, 7–11 July 2008. [Google Scholar] [CrossRef]
  29. Wang, W.; Wu, J.; Li, Y.; Pei, J.; Yang, J.; Liang, C. A Novel Anti-Deceptive Jamming Method for Multistatic SAR. In Proceedings of the 2019 IEEE International Geoscience and Remote Sensing Symposium, Yokohama, Japan, 28 July–2 August 2019. [Google Scholar] [CrossRef]
  30. Ma, Y. Research on SAR Countermeasure Techniques; UEST of China: Chengdu, China, 2004. [Google Scholar]
  31. Ma, X.; Qin, J.; He, Z.; Yang, J.; Lu, Q. Three-Channel Cancellation of SAR Blanketing Jamming Suppression. Acta Electron. Sin. 2007, 35, 1015–1020. [Google Scholar] [CrossRef]
  32. Huang, L.; Dong, C.; Shen, Z.; Zhao, G. Investigation on Countermeasure against InSAR Dual-Channel Cancellation Technique with Multi-antenna Jamming. J. Electron. Inf. Technol. 2015, 37, 913–918. [Google Scholar]
  33. Liu, A.; Wang, F.; Xu, H.; Li, L. N-SAR: A New Multichannel Multimode Polarimetric Airborne SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3155–3166. [Google Scholar] [CrossRef]
Figure 1. Jamming residual with IPSR error.
Figure 1. Jamming residual with IPSR error.
Sensors 22 09356 g001
Figure 2. IPDC image of 10 dB JSR.
Figure 2. IPDC image of 10 dB JSR.
Sensors 22 09356 g002
Figure 3. IPDC image of 0 dB JSR.
Figure 3. IPDC image of 0 dB JSR.
Sensors 22 09356 g003
Figure 4. IPSRs and IPDCs along azimuth (10 dB JSR).
Figure 4. IPSRs and IPDCs along azimuth (10 dB JSR).
Sensors 22 09356 g004
Figure 5. IPSRs and IPDCs along azimuth (0 dB JSR).
Figure 5. IPSRs and IPDCs along azimuth (0 dB JSR).
Sensors 22 09356 g005
Figure 6. Spectrum without jamming and with 0 dB jamming.
Figure 6. Spectrum without jamming and with 0 dB jamming.
Sensors 22 09356 g006
Figure 7. Cosine similarity with IPSR.
Figure 7. Cosine similarity with IPSR.
Sensors 22 09356 g007
Figure 8. IPSR estimation under different JSRs.
Figure 8. IPSR estimation under different JSRs.
Sensors 22 09356 g008
Figure 9. Jamming residual with co-registration error of different jamming bandwidth to signal bandwidth ratio.
Figure 9. Jamming residual with co-registration error of different jamming bandwidth to signal bandwidth ratio.
Sensors 22 09356 g009
Figure 10. Demonstration of azimuth resolution degradation. (a) SAR image before DCC. (b) SAR image after DCC.
Figure 10. Demonstration of azimuth resolution degradation. (a) SAR image before DCC. (b) SAR image after DCC.
Sensors 22 09356 g010
Figure 11. Geometric relationship of the SAR and the jammer.
Figure 11. Geometric relationship of the SAR and the jammer.
Sensors 22 09356 g011
Figure 12. Optical image of the imaging area and the location and orientation of the jammer.
Figure 12. Optical image of the imaging area and the location and orientation of the jammer.
Sensors 22 09356 g012
Figure 13. The jammer, the transmitting antenna, and the receiving antenna.
Figure 13. The jammer, the transmitting antenna, and the receiving antenna.
Sensors 22 09356 g013
Figure 14. The schematic diagram and processing chain. (a) Schematic diagram. (b) Processing chain.
Figure 14. The schematic diagram and processing chain. (a) Schematic diagram. (b) Processing chain.
Sensors 22 09356 g014
Figure 15. SAR image without jamming.
Figure 15. SAR image without jamming.
Sensors 22 09356 g015
Figure 16. SAR image with jamming.
Figure 16. SAR image with jamming.
Sensors 22 09356 g016
Figure 17. SAR image after DCC with relatively precise IPSR.
Figure 17. SAR image after DCC with relatively precise IPSR.
Sensors 22 09356 g017
Figure 18. SAR image after DCC with IPSR error of 18°.
Figure 18. SAR image after DCC with IPSR error of 18°.
Sensors 22 09356 g018
Table 1. The major time cost of the CS algorithm.
Table 1. The major time cost of the CS algorithm.
OperationInput Data SizeOutput Data Size Time   Cost   ( μ s )
Azimuth FFT N × K N 2 × K K × N 2 / 128 × 4
Range K points dot product N 2 × K N 2 × K N 2 × K / 128 × 2
Range FFT N 2 × K N 2 × K 2 N 2 × K 2 / 128 × 4
Range K 2 points dot product N 2 × K 2 N 2 × K 2 N 2 × K 2 / 128 × 2
Range IFFT N 2 × K 2 N 2 × K 2 N 2 × K 2 / 128 × 4
Azimuth N 2 points dot product N 2 × K 2 N 2 × K 2 N 2 × K 2 / 128 × 2
Azimuth IFFT N 2 × K 2 N 2 × K 2 N 2 × K 2 / 128 × 4
Table 2. The major time cost of our method.
Table 2. The major time cost of our method.
OperationInput Data SizeOutput Data SizeTime Cost (μs)
N times FFT N × K N × K 2 N × K 2 / 128 × 4
N times K 2 points dot product N × K 2 N × K 2 N × K 2 / 128 × 2
Table 3. Main parameters of Nriet-SAR L-band sensor.
Table 3. Main parameters of Nriet-SAR L-band sensor.
Parameter Value
Frequency bandL
Bandwidth100 MHz
Pulse repeat time (PRF)800 Hz
Noise equivalent sigma zero (NESZ)<−24 dB
Range resolution1.5 m
Azimuth resolution0.6 m
Channel number2
Channel distance0.6 m
Typical flight height7000 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, L.; Liu, A.; Huang, Z.; Xu, H.; Han, D. A Novel Method for Interferometric Phase Estimation in Dual-Channel Cancellation. Sensors 2022, 22, 9356. https://doi.org/10.3390/s22239356

AMA Style

Huang L, Liu A, Huang Z, Xu H, Han D. A Novel Method for Interferometric Phase Estimation in Dual-Channel Cancellation. Sensors. 2022; 22(23):9356. https://doi.org/10.3390/s22239356

Chicago/Turabian Style

Huang, Long, Aifang Liu, Zuzhen Huang, Hui Xu, and Dong Han. 2022. "A Novel Method for Interferometric Phase Estimation in Dual-Channel Cancellation" Sensors 22, no. 23: 9356. https://doi.org/10.3390/s22239356

APA Style

Huang, L., Liu, A., Huang, Z., Xu, H., & Han, D. (2022). A Novel Method for Interferometric Phase Estimation in Dual-Channel Cancellation. Sensors, 22(23), 9356. https://doi.org/10.3390/s22239356

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