CN103675768B - A kind of suppression method for satellite-borne SAR radio-frequency interference - Google Patents
A kind of suppression method for satellite-borne SAR radio-frequency interference Download PDFInfo
- Publication number
- CN103675768B CN103675768B CN201310693687.7A CN201310693687A CN103675768B CN 103675768 B CN103675768 B CN 103675768B CN 201310693687 A CN201310693687 A CN 201310693687A CN 103675768 B CN103675768 B CN 103675768B
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- radio
- decomposition
- data matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a kind of suppression method for satellite-borne SAR radio-frequency interference.Technical scheme comprises: the first step: build data matrix, and the range line echo according to gained builds data matrix; Second step: estimate covariance correlation matrix, the data matrix constructed by utilization estimates its covariance correlation matrix; 3rd step: feature decomposition, obtain eigenwert and characteristic of correspondence vector and sort; 4th step: wavelet decomposition estimates dominant eigenvalue number, the proper vector structure interference space corresponding to the dominant eigenvalue estimated by utilization; 5th step: AF panel and data reconstruction, completes wide band radio-frequency AF panel.The present invention all can obtain good inhibition to narrow radio frequency AF panel and wide band radio-frequency AF panel, and Robust Performance.
Description
Technical field
The invention belongs to the interleaving techniques field of space flight and microwave remote sensing, particularly Radio frequency interference (RFI) suppressing method in a kind of satellite-borne SAR (Synthetic Aperture Radar, synthetic-aperture radar) data.
Background technology
Spaceborne SAR System can complete wide swath, high resolving power earth observation, is the effective means of carrying out military mapping and battlefield surveillance.In order to meet growing military and civilian demand, the particularly active demand of more high resolving power and concealed target detection, the working frequency range of satellite-borne SAR is expanded gradually toward higher and lower both direction, the electromagnetic environment faced will become increasingly complex, particularly low band Spaceborne SAR System, a large amount of radio and television and communication system is there is within the scope of frequency spectrums of operation, they are in continuous duty, and emissive power is higher, even the bandwidth of operation of some system is wider, such as DVB-T (Digital VideoBroadcasting--Terrestrial, DTTB) system, bandwidth has 8MHz, these undesired signals enter into Spaceborne SAR System receiver, Space-borne SAR Imaging quality will be had a strong impact on, cause image interpretation and interpretation difficulty.
The SAR Radio frequency interference (RFI) suppressing method of the feature based Subspace Decomposition that Feng Zhou proposes is relatively good to narrow radio frequency interference suppressioning effect.But wide band radio-frequency is disturbed, due to be with wide, power is high, the large eigenwert that number is more can be produced, the error that the dominant eigenvalue number decision method that existing eigendecomposition method adopts can cause dominant eigenvalue number to judge when carrying out Eigenvalue Decomposition.Based on the advantage of wavelet decomposition in detection of Singular Point, the present invention proposes the suppression method for satellite-borne SAR radio-frequency interference of a kind of wavelet decomposition and eigendecomposition associating.This method all has good inhibition to narrow radio frequency interference and wide band radio-frequency interference, and Robust Performance.
Summary of the invention
The object of the invention is: the satellite-borne SAR wide band radio-frequency disturbance restraining method proposing a kind of wavelet decomposition and eigendecomposition associating, the method takes full advantage of the advantage of wavelet decomposition in detection of Singular Point, good inhibition is all had to narrow radio frequency interference and wide band radio-frequency interference, and Robust Performance.
The thinking of technical solution of the present invention is: the covariance correlation matrix first constructing raw data, feature decomposition is carried out to it, the method of wavelet decomposition is utilized to carry out dominant eigenvalue judgement, what dominant eigenvalue was corresponding is interference components, utilize dominant eigenvalue characteristic of correspondence vector structure interference space, then the mapping component of original signal on feature interference subspace is deducted, undesired signal can be curbed.
Technical solution of the present invention is:
Suppose that x represents the data that a certain orientation moment satellite-borne SAR receives, x=[x
1..., x
m..., x
m]
t, wherein M represents that distance is to sampling number, and subscript T represents matrix transpose operation.
The first step, structure data matrix
Data x is divided into the subvector of K L dimension, L determines according to actual needs, usual 32≤L≤128, wherein K=M-L+1.A definition kth sub-vector x
kfor:
x
k=[x
k,…x
k+L-1]
T,k=1,2,…,K
Then, the data matrix X of L × K is built:
X=[x
1,…,x
k]
Second step, estimate covariance correlation matrix
Data matrix X is utilized to estimate its covariance correlation matrix
3rd step, feature decomposition
Right
carry out Eigenvalues Decomposition, eigenwert is arranged according to descending order, obtain characteristic value sequence λ
1, λ
2..., λ
l, its corresponding proper vector is u
1, u
2..., u
l.
4th step, estimation dominant eigenvalue number, structure interference space
To characteristic value sequence λ
1, λ
2..., λ
lcarry out multiple dimensioned one-dimensional discrete wavelet decomposition, the wavelet coefficient c of i-th layer that obtains
ifor
c
i=[cA
i,cD
i,cD
i-1,…,cD
1],i=1,2,…,N,
Wherein cA
irepresent the low frequency coefficient that i-th layer of decomposition obtains; CD
irepresent the high frequency coefficient that i-th layer of decomposition obtains.N gets at most
,
represent and round downwards.
Utilize high frequency coefficient to be reconstructed signal and obtain detail signal d, obtain the sequence number that in detail signal d, maximum amplitude is corresponding, this sequence number is exactly the number r of dominant eigenvalue, then construct interference space F to be:
F=span{u
1,u
2,…,u
r},1≤r<L
5th step, AF panel and data reconstruction
Following formula is utilized to obtain interfering data I
k:
Following formula is utilized to calculate the signal after suppressing Radio frequency interference (RFI)
Thus, the data matrix after the Radio frequency interference (RFI) that is inhibited
Rearrange data matrix
namely the inverse process that the first step builds data matrix is implemented, the signal vector x ' after the Radio frequency interference (RFI) that is inhibited.
Employing the present invention is desirable obtains following technique effect:
The suppression method for satellite-borne SAR radio-frequency interference that the small echo that the present invention proposes is combined with eigendecomposition, can suppress, in interfering process, to utilize wavelet decomposition accurately to estimate dominant eigenvalue number in employing eigendecomposition method.Overcome existing Subspace Decomposition method and inaccurate defect is judged to broadband interference dominant eigenvalue number, all can obtain good inhibition to narrow radio frequency AF panel and wide band radio-frequency AF panel, and Robust Performance.
Accompanying drawing explanation
Fig. 1 is the principle flow chart of suppression method for satellite-borne SAR radio-frequency interference provided by the invention;
When Fig. 2 is for utilizing emulated data, the dominant eigenvalue number that the inventive method and existing method generate is estimated to contrast schematic diagram;
Fig. 3 is the imaging results of the raw data used in emulation experiment;
Fig. 4 is the imaging results of Fig. 3 being added to the rear data of wide band radio-frequency interference;
Fig. 5 utilizes existing method to carry out the imaging results after AF panel to Fig. 4;
Fig. 6 utilizes method provided by the invention to carry out the imaging results after AF panel to Fig. 4.
Embodiment
Fig. 1 is the principle flow chart of suppression method for satellite-borne SAR radio-frequency interference provided by the invention.Whole flow process is divided into five steps: with SAR data range line echoed signal for input, the first step: build data matrix, and the range line echo according to gained builds data matrix; Second step: estimate covariance correlation matrix, the data matrix constructed by utilization estimates its covariance correlation matrix; To second step, 3rd step: feature decomposition, estimates that the covariance correlation matrix obtained carries out feature decomposition, obtain eigenwert and characteristic of correspondence vector, carry out descending order arrangement to eigenwert, corresponding proper vector order also corresponds to the order of eigenwert.4th step: wavelet decomposition estimates dominant eigenvalue number, build interference space, multiple dimensioned one-dimensional discrete wavelet decomposition is carried out to the characteristic value sequence after sequence, obtain corresponding low frequency coefficient and high frequency coefficient, utilize high frequency coefficient to be reconstructed characteristic value sequence, obtain the detailed information of characteristic value sequence, obtain the sequence number that maximum amplitude in detail signal is corresponding, just can determine the number of dominant eigenvalue, the proper vector structure interference space corresponding to the dominant eigenvalue estimated by utilization; 5th step: AF panel and data reconstruction, data matrix is projected to the interfering data estimated by interference space acquisition, the interfering data deducting estimation from data matrix can obtain the data matrix after AF panel, data matrix after AF panel is rearranged and can obtain the vector of the echo after AF panel, complete the wide band radio-frequency AF panel of this range line.
When Fig. 2 is for utilizing emulated data, the dominant eigenvalue number that the inventive method and existing method generate is estimated to contrast schematic diagram.Emulated data utilizes one-dimensional linear FM signal to increase the DVB-T interfere generation of 8MHz, is with rhombus solid line in figure
represent the characteristic value sequence (according to arranging from big to small) of covariance correlation matrix gained after feature decomposition; Band square frame dotted line
represent and adopt the dominant eigenvalue estimated sequence that obtains of existing ratio calculation method, namely what wherein maximal value was corresponding 5 is judged to dominant eigenvalue number; The line of band Diamond spot
represent and adopt the dominant eigenvalue estimated sequence that obtains of method provided by the invention, namely what wherein maximal value was corresponding 7 is judged to dominant eigenvalue number; Dotted line " ... " represents that dominant eigenvalue number that existing ratio calculation method estimates corresponds to the position of primitive character value sequence; Dot-and-dash line
represent that dominant eigenvalue number that method provided by the invention obtains corresponds to the position of primitive character value sequence.As seen from the figure, compared with the conventional method, the inventive method can accurately estimate dominant eigenvalue number under wide band radio-frequency interference exists situation, and estimated result is better than existing method.
Fig. 3, Fig. 4, Fig. 5 and Fig. 6 are the results of carrying out emulation experiment.Emulation experiment have chosen a glitch-free RadarSat-1 echo measured data, and in these data, with the addition of the DVB-T undesired signal that bandwidth is 8MHz, jamming-to-signal ratio is 20dB.
Fig. 3 is the imaging results of the raw data used in emulation experiment.As can be seen from the figure, clear picture, contrast is high.Fig. 4 is the imaging results of Fig. 3 being added to the rear data of wide band radio-frequency interference.As can be seen from the figure, due to the existence of interference, image object information is capped, and picture contrast declines.Fig. 5 utilizes existing method to carry out the imaging results after AF panel to Fig. 4, and what wherein dominant eigenvalue number estimation method adopted is existing radiometer algorithm, and after AF panel, picture quality recovers less, and contrast is not high.Fig. 6 utilizes method provided by the invention to carry out the imaging results after AF panel to Fig. 4, and after AF panel, picture quality obviously promotes, and detailed information recovers comparatively obvious, and picture contrast is significantly improved.
Above emulation experiment describes the inventive method still estimating dominant eigenvalue number exactly by under wide band radio-frequency disturbing effect, ensure that the accuracy that interference space constructs, and wide band radio-frequency is disturbed can be effectively suppressed.Meanwhile, for narrow radio frequency interference, the present invention also can be suppressed result preferably.
Claims (2)
1. a suppression method for satellite-borne SAR radio-frequency interference, is characterized in that, comprises the steps,
Suppose that x represents the data that a certain orientation moment satellite-borne synthetic aperture radar receives, x=[x
1..., x
m..., x
m]
t, wherein M represents that distance is to sampling number, and subscript T represents matrix transpose operation;
The first step, structure data matrix:
Data x is divided into the subvector of K L dimension, wherein K=M-L+1; A definition kth sub-vector x '
kfor:
x′
k=[x
k,…x
k+L-1]
T,k=1,2,…,K
Then, the data matrix X of L × K is built:
X=[x′
1,…,x′
K]
Second step, estimate covariance correlation matrix:
Data matrix X is utilized to estimate its covariance correlation matrix
3rd step, feature decomposition:
To covariance correlation matrix
carry out Eigenvalues Decomposition, eigenwert is arranged according to descending order, obtain characteristic value sequence λ
1, λ
2..., λ
l, its corresponding proper vector is u
1, u
2..., u
l;
4th step, estimation dominant eigenvalue number, structure interference space:
To characteristic value sequence λ
1, λ
2..., λ
lcarry out multiple dimensioned one-dimensional discrete wavelet decomposition, the wavelet coefficient c of i-th layer that obtains
ifor
c
i=[cA
i,cD
i,cD
i-1,…,cD
1],i=1,2,…,N,
Wherein cA
irepresent the low frequency coefficient that i-th layer of decomposition obtains; CD
irepresent the high frequency coefficient that i-th layer of decomposition obtains; N maximal value is
represent and round downwards;
Utilize high frequency coefficient to be reconstructed signal and obtain detail signal d, obtain the sequence number that in detail signal d, maximum amplitude is corresponding, this sequence number is exactly the number r of dominant eigenvalue, then construct interference space F to be:
F=span{u
1,u
2,…,u
r},1≤r<L
5th step, AF panel and data reconstruction:
Following formula is utilized to obtain interfering data I
k:
Following formula is utilized to calculate the signal after suppressing Radio frequency interference (RFI)
Thus, the data matrix after the Radio frequency interference (RFI) that is inhibited
Rearrange data matrix
namely the inverse process that the first step builds data matrix is implemented, the signal vector x ' after the Radio frequency interference (RFI) that is inhibited.
2. suppression method for satellite-borne SAR radio-frequency interference according to claim 1, is characterized in that, 32≤L≤128.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310693687.7A CN103675768B (en) | 2013-12-17 | 2013-12-17 | A kind of suppression method for satellite-borne SAR radio-frequency interference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310693687.7A CN103675768B (en) | 2013-12-17 | 2013-12-17 | A kind of suppression method for satellite-borne SAR radio-frequency interference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675768A CN103675768A (en) | 2014-03-26 |
CN103675768B true CN103675768B (en) | 2015-10-21 |
Family
ID=50313909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310693687.7A Active CN103675768B (en) | 2013-12-17 | 2013-12-17 | A kind of suppression method for satellite-borne SAR radio-frequency interference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675768B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104101868B (en) * | 2014-06-30 | 2016-08-24 | 西安电子科技大学 | Radar multi-false-target jamming suppressing method based on interference space reconstruct |
US9537530B1 (en) * | 2015-06-26 | 2017-01-03 | Intel Corporation | Transceiver device and method of processing signals |
CN105974376B (en) * | 2016-05-09 | 2018-05-25 | 中国人民解放军国防科学技术大学 | A kind of SAR radio frequency interferences suppressing method |
CN106646472B (en) * | 2016-12-30 | 2019-01-18 | 中国人民解放军国防科学技术大学 | SAR radio frequency interference suppressing method based on null tone cascade filtering |
CN109672487B (en) * | 2018-12-13 | 2021-04-16 | 宁波连鸿电子科技有限公司 | Interference adjusting method of robustness self-adaptive variable load filter |
CN113655443A (en) * | 2021-08-26 | 2021-11-16 | 北京理工大学 | Low-frequency band SAR radio frequency interference suppression method for broadband digital television signal |
CN116299205B (en) * | 2023-05-17 | 2023-09-01 | 西安电子科技大学 | Time domain sliding window subspace projection SAR broadband interference suppression method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6292592B1 (en) * | 1998-10-19 | 2001-09-18 | Raytheon Company | Efficient multi-resolution space-time adaptive processor |
CN102243300B (en) * | 2011-04-15 | 2012-10-10 | 中国人民解放军国防科学技术大学 | Method for radio frequency interference suppression and error correction of low-frequency synthetic aperture radar |
-
2013
- 2013-12-17 CN CN201310693687.7A patent/CN103675768B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN103675768A (en) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675768B (en) | A kind of suppression method for satellite-borne SAR radio-frequency interference | |
CN105137410B (en) | The waveform optimization method of high-resolution radar communicating integral based on OFDM | |
CN104698459B (en) | Stripe SAR (specific absorption resolution) compressed sensing and imaging method for missing data | |
CN107271993B (en) | Scanning radar angle super-resolution imaging method based on maximum posterior | |
CN105974376B (en) | A kind of SAR radio frequency interferences suppressing method | |
CN103472449B (en) | BP wideband synthesis method based on MIMO image domain | |
CN103235295B (en) | Method for estimating small-scene radar target range images on basis of compression Kalman filtering | |
CN103630888A (en) | High-precision real-time microwave velocity and distance measuring device based on symmetrical triangle LFMCW (Linear Frequency Modulation Continuous Wave) radar | |
CN103364783B (en) | Moving target radial velocity non-fuzzy estimation method based on single-channel SAR (synthetic aperture radar) | |
CN107621635B (en) | Forward-looking sea surface target angle super-resolution method | |
CN104698431B (en) | Based on the multichannel SAR orientation ambiguity solution method that obscuring component DOA estimates | |
CN110275158A (en) | Wideband radar echo-signal method for parameter estimation based on Bayes's compressed sensing | |
CN102879783A (en) | Sparse detection frequency signal-based inverse synthetic aperture radar (ISAR) imaging method | |
CN103344948B (en) | Method for computing external illuminator radar cross-ambiguity function utilizing sparse Fourier transform | |
CN106338731A (en) | Multi-sub band concurrent MIMO-SAR radar imaging method | |
CN104237886B (en) | High-precision synthetic aperture radar imaging method | |
CN105301589A (en) | High-resolution wide-swath SAR (synthetic aperture radar) ground motion object imaging method | |
CN103983959A (en) | SAR system movement target radial speed estimation method based on data reconstruction | |
Li et al. | High speed maneuvering target detection based on joint keystone transform and CP function | |
CN107102328B (en) | FPGA-based real-time imaging signal processing method and FPGA | |
Lu et al. | Ionospheric decontamination and sea clutter suppression for HF skywave radars | |
CN115877382A (en) | Motion error estimation method based on adjacent pulse transformation difference of frequency modulated continuous wave | |
CN108919263B (en) | ISAR high-resolution imaging method based on maximum mutual information criterion | |
CN114545351A (en) | Maneuvering target coherent detection method and system based on range frequency axis inversion transformation and second-order WVD (WVD) | |
CN102508265A (en) | Signal separation estimation theory-based satellite navigation signal multipath interference suppression method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |