CN106772273B - A kind of SAR false target disturbance restraining method and system based on dynamic aperture - Google Patents

A kind of SAR false target disturbance restraining method and system based on dynamic aperture Download PDF

Info

Publication number
CN106772273B
CN106772273B CN201611040047.6A CN201611040047A CN106772273B CN 106772273 B CN106772273 B CN 106772273B CN 201611040047 A CN201611040047 A CN 201611040047A CN 106772273 B CN106772273 B CN 106772273B
Authority
CN
China
Prior art keywords
dynamic
target point
sar
matrix
false
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
Application number
CN201611040047.6A
Other languages
Chinese (zh)
Other versions
CN106772273A (en
Inventor
赵博
黄磊
张亮
廖斌
何春龙
张沛昌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen University
Original Assignee
Shenzhen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen University filed Critical Shenzhen University
Priority to CN201611040047.6A priority Critical patent/CN106772273B/en
Publication of CN106772273A publication Critical patent/CN106772273A/en
Application granted granted Critical
Publication of CN106772273B publication Critical patent/CN106772273B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9027Pattern recognition for feature extraction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Abstract

The invention discloses a kind of SAR false target disturbance restraining method and system based on dynamic aperture, method include: that the blended data of basis generates data matrix and reference signal;Data matrix and reference signal carry out imaging, generate the one-dimensional picture in orientation, one-dimensional as generating aiming spot set to be estimated according to orientation;Dynamic observing matrix is generated according to preset dynamic aperture window and aiming spot set;Minimum solution, the isolated one-dimensional picture of the corresponding true bearing of real target point vector and the corresponding one-dimensional picture of false bearing of false target point vector are carried out using the observation characterization matrix and dynamic observing matrix of construction;It repeats the above process, the real goal two dimensional image after being reconstructed and false target two dimensional image.The present invention, by constructing the objective function of cheating interference separation problem and optimizing, effectively inhibits the interference of false target, reduces the high expense of system upgrade under the premise of not improving existing radar system complexity.

Description

A kind of SAR false target disturbance restraining method and system based on dynamic aperture
Technical field
The present invention relates to signal processing technology field more particularly to a kind of SAR false target interference based on dynamic aperture Suppressing method and system.
Background technique
SAR deception jammer carries out delay and phase tune according to SAR signal of the pre-designed false scene to intercepting and capturing System, it is lower than power needed for traditional compacting interference, and endanger bigger.SAR, that is, English " Synthetic Aperture The abbreviation of Radar ", means synthetic aperture radar.Jammer generates false target true to nature in real scene, to upset The acquisition of information and decision of SAR.Therefore the extensive concern of scholars is caused for the research of SAR cheating interference, also correspondingly More significant achievement is achieved, such that SAR deception jammer is more life-like under the operating mode of Various Complex, more increases Effect ground generates more false targets, to increase the menace of cheating interference.Therefore, in order to promote the target of SAR system Identification capability needs to study corresponding false target suppressing method and real scene is reconstructed.
Zhao Shanshan etc. is in document " Discrimination of Deception Targets in Multistatic Radar Based on Clustering Analysis”(IEEE Sensors Journal,Vol.16, No.8, Apr.2016:2500-2508) in the system architecture based on more base radars, in amplitude proportional feature space using layering The method of classification analysis analyzes the dispersion characteristic of cheating interference target and real goal, and by optimizing classification number The classification of false target and real goal is realized with the method that designs different classes of minimum cost.Although this method can be one It is completed in a pulse recurrence interval, but needs multistatic radar system as support, system complexity and higher cost.
Xiong Wei etc. is in document " Trilinear decomposition-based spatial- polarisational filter method for deception jamming suppression of radar”(IET Radar, Sonar& Navigation, Vol.10, No.4, Mar.2016:765-773) utilize three linear unmixed models acquisition mesh The Combined estimator of the pairs of two-dimentional arrival direction and polarization parameter of mark and jammer, and based on this construction spatial polarizations oblique projection filter Blended data is projected to each sublinear space of target to inhibit to interfere by wave device.Although this method can effectively inhibit deception dry It disturbs, but system structure is complicated, higher cost.Therefore SAR in the prior art is in reconstruction process, since cheating interference result presses down Effect processed is poor, poor so as to cause real goal imaging results, and real scene quality reconstruction is poor, to influence the acquisition of information of radar With decision, make troubles later period acquisition of information.
Therefore, the existing technology needs to be improved and developed.
Summary of the invention
In view of the deficiencies in the prior art, it is an object of that present invention to provide a kind of, and the SAR false target based on dynamic aperture is dry Disturb suppressing method and system, it is intended to SAR in the prior art is solved in reconstruction process, due to cheating interference result inhibitory effect Difference, poor so as to cause real goal imaging results, real scene quality reconstruction is poor, to influence the acquisition of information of radar and determine Plan, the technical issues of making troubles later period acquisition of information.
Technical scheme is as follows:
A kind of SAR false target disturbance restraining method based on dynamic aperture, wherein method comprising steps of
A, the blended data for obtaining true SAR scene echoes and cheating interference is corrected blended data raw after handling At the data matrix of image scene, the reference signal of frequency modulation removal is generated according to SAR system parameter;
B, the signal that frequency modulation removal is generated after a distance unit in data matrix is multiplied point by point with reference signal, to going to adjust The signal of frequency carries out imaging, generates the one-dimensional picture in orientation and is generated to be estimated according to the local peaking position of the one-dimensional picture in orientation Aiming spot set;
C, the signal of frequency modulation removal is extracted according to pre-set dynamic aperture window, and carried out at frequency modulation removal imaging Reason generates dynamic aperture imaging results;Extract the data width of dynamic aperture imaging results and aiming spot set corresponding position Degree generates dynamic observing matrix;
D, the amplitude loss coefficient of real target point Yu false target point is calculated separately, construction observes characterization matrix, according to After characterization matrix and dynamic observing matrix carry out minimization problem solution, isolated real target point vector is corresponding true The one-dimensional picture in orientation and the corresponding one-dimensional picture of false bearing of false target point vector;
E, according to aiming spot set that the element in vector is one-dimensional to the one-dimensional picture of true bearing and false bearing respectively Corresponding point carries out assignment as in;
F, it repeats the above steps B, step C, step D, step to the one-dimensional picture in orientation in data matrix different distance unit E, real goal two dimensional image and false target two dimensional image after respectively obtaining reconstruct.
The SAR false target disturbance restraining method based on dynamic aperture, wherein the step A specifically includes step It is rapid:
A1, the blended data for obtaining true SAR scene echoes and SAR cheating interference, carry out carrier wave, Range compress, away from After migration correction process, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), column direction table in data matrix Show distance dimension, line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;
A2, the reference signal that azimuth dimension deramp processing is generated according to SAR system parameter, are denoted as s0(ta);Specifically:
Wherein fdcFor doppler centroid, γ is doppler frequency rate.
The SAR false target disturbance restraining method based on dynamic aperture, wherein the step B is specially to walk It is rapid:
B1, data matrix s (t is obtainedr,ta) in k-th of distance unit data sk(ta), with reference signal s0(ta) by Point generates the signal s of frequency modulation removal after being multipliedDk(ta);
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates target point Scattering coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is radar signal Wavelength;B2, to sDk(ta) Fourier transformation is carried out, it is one-dimensional as I to generate the orientation comprising real target point and false target pointk (fa);
WhereinIndicate Fourier transformation;
B3, I is obtainedk(fa) local peaking position, generate aiming spot set D to be estimatedf
The SAR false target disturbance restraining method based on dynamic aperture, wherein the step C specifically includes step It is rapid:
C1, acquisition center are ηp, width is the dynamic aperture window w of Tp(ta)
C2, according to wp(ta) to sDk(ta) extract, after extraction carry out frequency modulation removal imaging, generate dynamic aperture at As result Ikp(fa)
C3, to dynamic aperture imaging results Ikp(fa) in aiming spot set DfThe data amplitude of corresponding position carries out It extracts, is stored as the pth row that dynamic observes data matrix U;
C4, step C1, C2, C3 are repeated, until obtaining complete dynamic observing matrix U.
The SAR false target disturbance restraining method based on dynamic aperture, wherein the step D specifically includes step It is rapid:
D1, the amplitude loss factor alpha (η for calculating separately real target pointpf) with the amplitude loss factor beta of false target point (ηp)
D2, the center η according to dynamic aperture windowpAnd aiming spot set DfCorresponding τf, calculate separately α (ηpf) With β (ηp) observation characterization matrix A and B are constructed afterwards;
D3, minimization problem is solvedIsolated real target point vector p is corresponding The one-dimensional picture of the true bearing one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q= diag(q)。
A kind of SAR false target Interference Suppression System based on dynamic aperture, wherein system includes:
Reference signal generation module, for obtaining the blended data of true SAR scene echoes and cheating interference, to mixed number According to the data matrix for generating image scene after processing is corrected, the reference signal of frequency modulation removal is generated according to SAR system parameter;
Aiming spot set generation module generates the letter of frequency modulation removal after being multiplied point by point for data matrix with reference signal Number, imaging is carried out to the signal of frequency modulation removal, generates the one-dimensional picture in orientation, it is raw according to the local peaking position of the one-dimensional picture in orientation At the aiming spot set comprising aiming spot;
Dynamic observing matrix generation module, for being mentioned according to pre-set dynamic aperture window to the signal of frequency modulation removal It takes, and carries out frequency modulation removal imaging, generate dynamic aperture imaging results;Extract dynamic aperture imaging results and aiming spot Gather the data amplitude of corresponding position, generates dynamic observing matrix;
False target point separation module, for calculating separately the amplitude loss coefficient of real target point Yu false target point, Construction observation characterization matrix, after carrying out minimization problem solution according to characterization matrix and dynamic observing matrix, what is separated is true The one-dimensional picture of the corresponding true bearing of real target point vector and the corresponding one-dimensional picture of false bearing of false target point vector;
Assignment module, for according to aiming spot set by the element in vector respectively to the one-dimensional picture of true bearing and void Corresponding point carries out assignment in the one-dimensional picture of false orientation;
Image reconstruction module obtains after handling the one-dimensional picture in orientation in data matrix different distance unit Real goal two dimensional image and false target two dimensional image after reconstruct.
The SAR false target Interference Suppression System based on dynamic aperture, wherein the reference signal generation module It specifically includes:
Data correction unit is unloaded for obtaining the blended data of true SAR scene echoes and SAR cheating interference After wave, Range compress, range migration correction processing, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), data Column direction indicates distance dimension in matrix, and line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;
Reference signal generation unit, for generating the reference signal of azimuth dimension deramp processing, note according to SAR system parameter For s0(ta);Specifically:
Wherein fdcFor doppler centroid, γ is doppler frequency rate.
The SAR false target Interference Suppression System based on dynamic aperture, wherein the aiming spot collection symphysis It is specifically included at module:
Frequency modulation removal signal generation unit, for obtaining data matrix s (tr,ta) in k-th of distance unit data sk (ta), with reference signal s0(ta) it is point-by-point be multiplied after generate the signal s of frequency modulation removalDk(ta);
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates target point Scattering coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is radar signal Wavelength;
Orientation is one-dimensional as generation unit, for sDk(ta) Fourier transformation is carried out, it generates comprising real target point and void The orientation of decoy point is one-dimensional as Ik(fa);
WhereinIndicate Fourier transformation;
Aiming spot set generation unit, for obtaining Ik(fa) local peaking position, generate target point to be estimated Location sets Df
The SAR false target Interference Suppression System based on dynamic aperture, wherein the dynamic observing matrix generates Module specifically includes:
Dynamic aperture window acquiring unit is η for obtaining centerp, width is the dynamic aperture window w of Tp(ta)
Dynamic aperture imaging results generation unit, for according to wp(ta) to sDk(ta) extract, it carries out adjusting after extraction Frequency imaging generates dynamic aperture imaging results Ikp(fa)
Dynamic observing matrix generation unit, for dynamic aperture imaging results Ikp(fa) in aiming spot set Df The data amplitude of corresponding position extracts, and is stored as the pth row that dynamic observes data matrix U, is obtained according to the data of each row Complete dynamic observing matrix U.
The SAR false target Interference Suppression System based on dynamic aperture, wherein described image reconstructed module is specific Include:
Loss coefficient computing unit, for calculating separately the amplitude loss factor alpha (η of real target pointpf) and false mesh Amplitude loss factor beta (the η of punctuatep)
Characterization matrix structural unit, for the center η according to dynamic aperture windowpAnd aiming spot set DfIt is corresponding τf, calculate separately α (ηpf) and β (ηp) observation characterization matrix A and B are constructed afterwards;
Minimization problem solves unit, for solving minimization problemWhat is separated is true The corresponding one-dimensional picture of true bearing of the real target point vector p one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q= diag(q)。
The SAR false target disturbance restraining method and system, the present invention that the present invention provides a kind of based on dynamic aperture can Under the premise of not improving existing single channel radar system complexity, according to the real target point and falseness at different direction position The amplitude loss variation characteristic of target point constructs the objective function of cheating interference separation problem and optimizes, reduces The high expense of system upgrade.
Detailed description of the invention
Fig. 1 is a kind of preferred embodiment of SAR false target disturbance restraining method based on dynamic aperture of the invention Flow chart.
Fig. 2 is that a kind of concrete application of SAR false target disturbance restraining method based on dynamic aperture of the invention is implemented Example there are the real goal imaging results schematic diagrames of cheating interference.
Fig. 3 is that a kind of concrete application of SAR false target disturbance restraining method based on dynamic aperture of the invention is implemented Real goal imaging results schematic diagram after the cheating interference separation of example.
Fig. 4 is that a kind of concrete application of SAR false target disturbance restraining method based on dynamic aperture of the invention is implemented Example there are the false target imaging results schematic diagrames of cheating interference.
Fig. 5 is a kind of preferred embodiment of SAR false target Interference Suppression System based on dynamic aperture of the invention Functional schematic block diagram.
Specific embodiment
To make the purpose of the present invention, technical solution and effect clearer, clear and definite, below to the present invention further specifically It is bright.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not intended to limit the present invention.
The preferred embodiment of the present invention also provides a kind of SAR false target Interference Suppression System based on dynamic aperture Flow chart, as shown in Figure 1, method includes:
Step S100, the blended data for obtaining true SAR scene echoes and cheating interference, is corrected place to blended data The data matrix that image scene is generated after reason, the reference signal of frequency modulation removal is generated according to SAR system parameter.
When it is implemented, step S100 is specifically included:
Step S101, the blended data for obtaining true SAR scene echoes and SAR cheating interference carries out carrier wave, distance pressure After contracting, range migration correction processing, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), column side in data matrix To distance dimension is indicated, line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;
Step S102, the reference signal that azimuth dimension deramp processing is generated according to SAR system parameter, is denoted as s0(ta);Specifically Are as follows:
Wherein fdcFor doppler centroid, γ is doppler frequency rate.
When it is implemented, trIt is distance to the fast time, i.e., after imaging radar emits signal to ground scene, to scene reflectivity The signal returned carries out digital sample, it is assumed that sample rate Fs, sampling number Nr, then the fast time can be expressed as [- Nr/2:Nr/ 2-1]/Fs, wherein [- Nr/2:Nr/ 2-1] it indicates by-Nr/ 2 are incremented to NrThe vector of/2-1, incremental spacing 1;
taFor the orientation slow time, i.e., imaging radar, which is repeated as many times, emits signal to ground scene and receives, duplicate frequency Rate is known as pulse recurrence frequency (PRF, pulse repetition frequency), number of repetition Na, then the slow time indicates For [- Na/2:Na/2-1]/PRF。
The wideband-radar signal of imaging radar transmitting are as follows:
Wherein Tp refers to the pulse width of radar emission signal, and fc refers to that the carrier frequency of signal, K refer to the frequency modulation rate of transmitting signal;
Carrier processing is gone to refer to signal s (tr) and radiofrequency signal exp (- j2 π fctr) be multiplied point by point
Range compress refers to that the broadband imaging signal to imaging radar transmitting carries out matched filtering or frequency modulation removal, frequency modulation removal It can be regarded as solution line frequency modulation, English is that Dechirp is handled, and hereinafter replaces frequency modulation removal with Dechirp.Concrete operations are as follows:
Matched filtering refers to signal s1(tr) and reference signalIt carries out as follows Operation
WhereinIndicate Fourier transformation,Indicate inverse Fourier transform, i.e., by s1(tr) and sref1(tr) respectively Be multiplied point by point after Fourier transformation, then the result of multiplication is subjected to inverse Fourier transform.
Dechirp refers to signal s1(tr) and reference signalIt proceeds as follows
I.e. by s1(tr) and sref2(tr) it is point-by-point be multiplied after carry out Fourier transformation.
Range migration correction refer to according to radar to observation scene distance change rule adjust the distance compressed signal into Row correction, specific bearing calibration is the prior art, and details are not described herein again.
Step S200, the letter of frequency modulation removal is generated after a distance unit in data matrix is multiplied point by point with reference signal Number, imaging is carried out to the signal of frequency modulation removal, generates the one-dimensional picture in orientation, it is raw according to the local peaking position of the one-dimensional picture in orientation At aiming spot set to be estimated.
When it is implemented, wherein to the signal of frequency modulation removal carry out imaging can be Fourier transformation or using CS or CLEAN technology carries out imaging.Wherein CS refers to that compressed sensing (compressive sensing) is imaged, specifically, CS imaging technique and CLEAN technology are also the prior art in SAR related fields, and details are not described herein again.
In further embodiment, step S200 is specifically included:
Step S201, data matrix s (t is obtainedr,ta) in k-th of distance unit sk(ta), with reference signal s0(ta) by Point generates the signal s of frequency modulation removal after being multipliedDk(ta);Assuming that including M real target point and N number of false target point in the row:
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates target point Scattering coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is radar signal Wavelength;
Step S202, to sDk(ta) Fourier transformation is carried out, generate the orientation comprising real target point and false target point One-dimensional picture Ik(fa);
WhereinIndicate Fourier transformation.
Step S203, I is obtainedk(fa) local peaking position, generate aiming spot set D to be estimatedf
Step S300, the signal of frequency modulation removal is extracted according to pre-set dynamic aperture window, and carries out frequency modulation removal Imaging generates dynamic aperture imaging results;Extract dynamic aperture imaging results and aiming spot set corresponding position Data amplitude generates dynamic observing matrix.
Specifically, DfIt is a set, the position of the element representation target point of the inside, in Ik(fa) in find local peaking, By the Position Number of local peaking (or corresponding delay τf) store into set DfLocal peaking is target point, but this local peaks Value is likely to be real target point and false target point is superimposed being formed as a result, so back separates it.
When it is implemented, step S300 is specifically included:
Step S301, obtaining center is ηp, width is the dynamic aperture window w of Tp(ta)
Step S302, according to wp(ta) to sDk(ta) extract, frequency modulation removal imaging is carried out after extraction, generates dynamic Aperture imaging result Ikp(fa)
Step S303, to dynamic aperture imaging results Ikp(fa) in aiming spot set DfThe data width of corresponding position Degree extracts, and is stored as the pth row that dynamic observes data matrix U;
Step S304, step step S301, step S302, step S303 are repeated, until obtaining complete dynamic observation square Battle array U.
Step S400, the amplitude loss coefficient of real target point Yu false target point, construction observation characterization square are calculated separately Battle array, after carrying out minimization problem solution according to characterization matrix and dynamic observing matrix, isolated real target point vector pair The one-dimensional picture of the true bearing answered and the corresponding one-dimensional picture of false bearing of false target point vector.
When it is implemented, step S400 is specifically included:
Step S401, the amplitude loss factor alpha (η of real target point is calculated separatelypf) damaged with the amplitude of false target point Lose factor beta (ηp)
Step S402, according to the center η of dynamic aperture windowpAnd aiming spot set DfCorresponding τf, calculate separately α (ηpf) and β (ηp) observation characterization matrix A and B are constructed afterwards;
Assuming that set DfLength be N, the dynamic aperture center η of selectionpNumber be that M (has the η that M is differentp's Value), utilize M different ηpValue according to α (ηpf)、β(ηp) can calculate separately to obtain the column vector a and b of 2 M × 1. Then in different τfValue under repeat the step, N number of column vector can be respectively obtained, these column vectors are arranged successively to obtain Matrix A and B.So matrix A and B are the matrixes of M × N.
Step S403, minimization problem is solvedIt is separated The corresponding one-dimensional picture of true bearing of the real target point vector p one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q= diag(q)。
Specifically, minimization problem is solved according to the CVX kit in MATLAB, isolated real goal The corresponding one-dimensional picture of true bearing of the point vector p one-dimensional picture of false bearing corresponding with false target point vector q.
Step S500, according to aiming spot set by the element in vector respectively to the one-dimensional picture of true bearing and falseness side Corresponding point carries out assignment in the one-dimensional picture in position.
When it is implemented, according to aiming spot set DfBy the element in each vector respectively to the one-dimensional picture of true bearing and Corresponding point carries out assignment in the one-dimensional picture of false bearing, and wherein the one-dimensional picture of true bearing and the one-dimensional picture of false bearing can also be weighed referred to as The one-dimensional picture of structure.The one-dimensional of reconstruct seems the real target point vector sum false target point vector referred to after separation, can most started Zero vector is helped in first definition, and the value of target point out to be estimated is later to corresponding element assignment.
Step S600, it repeats the above steps S200, step to the one-dimensional picture in orientation in data matrix different distance unit S300, step S400, step S500, real goal two dimensional image and false target two dimensional image after respectively obtaining reconstruct.
When it is implemented, repeating the above steps S200, step to the one-dimensional picture in orientation in data matrix different distance unit S300, step S400, step S500, until reconstruct is respectively completed to distance unit all in data matrix, after obtaining reconstruct Real goal two dimensional image and false target two dimensional image.
The concrete application of SAR false target disturbance restraining method of the present invention also by a kind of based on dynamic aperture is implemented Example, in particular emulation experiment, when emulation, are emulated using MATLAB software.
The parameter for emulating data is as follows: simulation objectives have 5 × 5 scattering dot matrix to constitute, and wherein real goal is arranged in Z Shape, false target are arranged in S-shaped.Fig. 2 show the real goal imaging results there are cheating interference, and it is dry that Fig. 3 show deception Real goal imaging results after disturbing separation, Fig. 4 show the false target imaging results after cheating interference separation.
By above method embodiment it is found that dynamic aperture of the present invention observation divides the SAR data there are cheating interference Analysis, can be under the premise of not improving existing single channel radar system complexity, according to the real goal at different direction position The amplitude loss variation characteristic of point and false target point constructs the objective function of cheating interference separation problem and optimizes and asks Solution improves conventional one-channel broadband imaging radar electronic warfare and takes advantage of to realize the cheating interference separation in single-channel SAR imaging Deceive the ability of interference.Different from the cheating interference suppressing method based on multi channel imaging radar, method proposed by the present invention avoids The high expense that single channel upgrades to multi-channel system is a kind of more economical effective anti-deceptive interference method.
A kind of preferred embodiment of SAR false target Interference Suppression System based on dynamic aperture is additionally provided in the present invention Functional schematic block diagram, as shown in figure 5, system includes:
Reference signal generation module 100, for obtaining the blended data of true SAR scene echoes and cheating interference, to mixed The data matrix that data are corrected generation image scene after processing is closed, is believed according to the reference that SAR system parameter generates frequency modulation removal Number;Specifically as described in embodiment of the method.
Aiming spot set generation module 200, generates frequency modulation removal after being multiplied point by point for data matrix with reference signal Signal, imaging is carried out to the signal of frequency modulation removal, generates the one-dimensional picture in orientation, according to the local peaking of the one-dimensional picture in orientation position It sets, generates the aiming spot set comprising aiming spot;Specifically as described in embodiment of the method.
Dynamic observing matrix generation module 300, for according to pre-set dynamic aperture window to the signal of frequency modulation removal into Row extracts, and carries out frequency modulation removal imaging, generates dynamic aperture imaging results;Extract dynamic aperture imaging results and target point The data amplitude of location sets corresponding position generates dynamic observing matrix;Specifically as described in embodiment of the method.
False target point separation module 400, for calculating separately the amplitude loss system of real target point Yu false target point Number, construction observation characterization matrix are separated after carrying out minimization problem solution according to characterization matrix and dynamic observing matrix The one-dimensional picture of the corresponding true bearing of real target point vector and the corresponding one-dimensional picture of false bearing of false target point vector;Specifically such as Described in embodiment of the method.
Assignment module 500, for according to aiming spot set by the element in vector respectively to the one-dimensional picture of true bearing Assignment is carried out with point corresponding in the one-dimensional picture of false bearing;Specifically as described in embodiment of the method.
Image reconstruction module 600 obtains after handling the one-dimensional picture in orientation in data matrix different distance unit Real goal two dimensional image and false target two dimensional image after to reconstruct;Specifically as described in embodiment of the method.
The SAR false target Interference Suppression System based on dynamic aperture, wherein the reference signal generation module It specifically includes:
Data correction unit is unloaded for obtaining the blended data of true SAR scene echoes and SAR cheating interference After wave, Range compress, range migration correction processing, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), data Column direction indicates distance dimension in matrix, and line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;Specifically As described in embodiment of the method.
Reference signal generation unit, for generating the reference signal of azimuth dimension deramp processing, note according to SAR system parameter For s0(ta);Specifically:
Wherein fdcFor doppler centroid, γ is doppler frequency rate;Specifically as described in embodiment of the method.
The SAR false target Interference Suppression System based on dynamic aperture, wherein the aiming spot collection symphysis It is specifically included at module:
Frequency modulation removal signal generation unit, for obtaining data matrix s (tr,ta) in k-th of distance unit data sk (ta), with reference signal s0(ta) it is point-by-point be multiplied after generate the signal s of frequency modulation removalDk(ta);
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates target point Scattering coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is radar signal Wavelength;Specifically as described in embodiment of the method.
Orientation is one-dimensional as generation unit, for sDk(ta) Fourier transformation is carried out, it generates comprising real target point and void The orientation of decoy point is one-dimensional as Ik(fa);
WhereinIndicate Fourier transformation;Specifically as described in embodiment of the method.
Aiming spot set generation unit, for obtaining Ik(fa) local peaking position, generate target point to be estimated Location sets Df;Specifically as described in embodiment of the method.
The SAR false target Interference Suppression System based on dynamic aperture, wherein the dynamic observing matrix generates Module specifically includes:
Dynamic aperture window acquiring unit is η for obtaining centerp, width is the dynamic aperture window w of Tp(ta)
Specifically as described in embodiment of the method.
Dynamic aperture imaging results generation unit, for according to wp(ta) to sDk(ta) extract, it carries out adjusting after extraction Frequency imaging generates dynamic aperture imaging results Ikp(fa)
Specifically as described in embodiment of the method.
Dynamic observing matrix generation unit, for dynamic aperture imaging results Ikp(fa) in aiming spot set Df The data amplitude of corresponding position extracts, and is stored as the pth row that dynamic observes data matrix U, is obtained according to the data of each row Complete dynamic observing matrix U;Specifically as described in embodiment of the method.
The SAR false target Interference Suppression System based on dynamic aperture, wherein described image reconstructed module is specific Include:
Loss coefficient computing unit, for calculating separately the amplitude loss factor alpha (η of real target pointpf) and false mesh Amplitude loss factor beta (the η of punctuatep)
Specifically as described in embodiment of the method.
Characterization matrix structural unit, for the center η according to dynamic aperture windowpAnd aiming spot set DfIt is corresponding τf, calculate separately α (ηpf) and β (ηp) observation characterization matrix A and B are constructed afterwards;Specifically as described in embodiment of the method.
Minimization problem solves unit, for solving minimization problemWhat is separated is true The corresponding one-dimensional picture of true bearing of the real target point vector p one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q= diag(q);Specifically as described in embodiment of the method.
In conclusion the invention proposes a kind of SAR false target disturbance restraining method and system based on dynamic aperture, Method includes: to generate data matrix and reference signal according to the blended data of true SAR scene echoes and cheating interference;Data Matrix carries out imaging after being multiplied point by point with reference signal, generate the one-dimensional picture in orientation, one-dimensional to be estimated as generating according to orientation Aiming spot set;It is mentioned according to the data amplitude of preset dynamic aperture window and aiming spot set corresponding position It takes, generates dynamic observing matrix;Construction observation characterization matrix, is minimized using observation characterization matrix and dynamic observing matrix It solves, the isolated one-dimensional picture of the corresponding true bearing of real target point vector and the corresponding falseness side of false target point vector The one-dimensional picture in position;Respectively in different distance unit the one-dimensional picture of true bearing and the one-dimensional picture of false bearing repeat the above process, obtain Real goal two dimensional image and false target two dimensional image after to reconstruct.The present invention can not improve existing single channel radar system Under the premise of complexity of uniting, changed according to the amplitude loss of real target point and false target point at different direction position special Sign constructs the objective function of cheating interference separation problem and optimizes, reduces the high expense of system upgrade.
It should be understood that the application of the present invention is not limited to the above for those of ordinary skills can With improvement or transformation based on the above description, all these modifications and variations all should belong to the guarantor of appended claims of the present invention Protect range.

Claims (10)

1. a kind of SAR false target disturbance restraining method based on dynamic aperture, which is characterized in that the method includes the steps:
A, the blended data for obtaining true SAR scene echoes and cheating interference, blended data is corrected after processing generate at The data matrix of image field scape generates the reference signal of frequency modulation removal according to SAR system parameter;
B, the signal that frequency modulation removal is generated after a distance unit in data matrix is multiplied point by point with reference signal, to frequency modulation removal Signal carries out imaging, generates the one-dimensional picture in orientation according to the local peaking position of the one-dimensional picture in orientation and generates target to be estimated Point location sets;
C, the signal of frequency modulation removal is extracted according to pre-set dynamic aperture window, and carries out frequency modulation removal imaging, it is raw At dynamic aperture imaging results;The data amplitude of dynamic aperture imaging results and aiming spot set corresponding position is extracted, it is raw At dynamic observing matrix;
D, the amplitude loss coefficient of real target point Yu false target point, construction observation characterization matrix, according to characterization are calculated separately After matrix and dynamic observing matrix carry out minimization problem solution, the isolated corresponding true bearing of real target point vector One-dimensional picture and the corresponding one-dimensional picture of false bearing of false target point vector;
E, according to aiming spot set by the element in vector respectively in the one-dimensional picture of true bearing and the one-dimensional picture of false bearing Corresponding point carries out assignment;
F, it repeats the above steps B, step C, step D, step E to the one-dimensional picture in orientation in data matrix different distance unit, point Real goal two dimensional image and false target two dimensional image after not reconstructed.
2. the SAR false target disturbance restraining method according to claim 1 based on dynamic aperture, which is characterized in that institute It states step A and specifically includes step:
A1, the blended data for obtaining true SAR scene echoes and SAR cheating interference, carry out carrier wave, Range compress, distance and move After dynamic correction process, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), in data matrix column direction indicate away from From dimension, line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;
A2, the reference signal that azimuth dimension deramp processing is generated according to SAR system parameter, are denoted as s0(ta);Specifically:
Wherein fdcFor doppler centroid, γ is doppler frequency rate.
3. the SAR false target disturbance restraining method according to claim 2 based on dynamic aperture, which is characterized in that institute Stating step B is specially step:
B1, data matrix s (t is obtainedr,ta) in k-th of distance unit data sk(ta), with reference signal s0(ta) be multiplied point by point The signal s of frequency modulation removal is generated afterwardsDk(ta);
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates the scattering of target point Coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is the wave of radar signal It is long;B2, to sDk(ta) Fourier transformation is carried out, it is one-dimensional as I to generate the orientation comprising real target point and false target pointk(fa);
WhereinIndicate Fourier transformation;
B3, I is obtainedk(fa) local peaking position, generate aiming spot set D to be estimatedf
4. the SAR false target disturbance restraining method according to claim 3 based on dynamic aperture, which is characterized in that institute It states step C and specifically includes step:
C1, acquisition center are ηp, width is the dynamic aperture window w of Tp(ta)
C2, according to wp(ta) to sDk(ta) extract, frequency modulation removal imaging is carried out after extraction, generates dynamic aperture imaging knot Fruit Ikp(fa)
C3, to dynamic aperture imaging results Ikp(fa) in aiming spot set DfThe data amplitude of corresponding position extracts, It is stored as the pth row that dynamic observes data matrix U;
C4, step C1, C2, C3 are repeated, until obtaining complete dynamic observing matrix U.
5. the SAR false target disturbance restraining method according to claim 4 based on dynamic aperture, which is characterized in that institute It states step D and specifically includes step:
D1, the amplitude loss factor alpha (η for calculating separately real target pointpf) with the amplitude loss factor beta (η of false target pointp)
D2, the center η according to dynamic aperture windowpAnd aiming spot set DfCorresponding τf, calculate separately α (ηpf) and β (ηp) observation characterization matrix A and B are constructed afterwards;
D3, minimization problem is solvedThe isolated corresponding true side of real target point vector p The one-dimensional picture in the position one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q=diag (q)。
6. a kind of SAR false target Interference Suppression System based on dynamic aperture, which is characterized in that system includes:
Reference signal generation module, for obtaining the blended data of true SAR scene echoes and cheating interference, to blended data into The data matrix that image scene is generated after row correction process, the reference signal of frequency modulation removal is generated according to SAR system parameter;
Aiming spot set generation module, it is raw after being multiplied point by point for a distance unit in data matrix with reference signal At the signal of frequency modulation removal, imaging is carried out to the signal of frequency modulation removal, the one-dimensional picture in orientation is generated, according to the part of the one-dimensional picture in orientation Peak position generates aiming spot set to be estimated;
Dynamic observing matrix generation module, for being extracted according to pre-set dynamic aperture window to the signal of frequency modulation removal, And frequency modulation removal imaging is carried out, generate dynamic aperture imaging results;Extract dynamic aperture imaging results and aiming spot collection The data amplitude of corresponding position is closed, dynamic observing matrix is generated;
False target point separation module is constructed for calculating separately the amplitude loss coefficient of real target point Yu false target point Characterization matrix is observed, after carrying out minimization problem solution according to characterization matrix and dynamic observing matrix, the true mesh that is separated The one-dimensional picture of the corresponding true bearing of punctuate vector and the corresponding one-dimensional picture of false bearing of false target point vector;
Assignment module, for according to aiming spot set by the element in vector respectively to the one-dimensional picture of true bearing and falseness side Corresponding point carries out assignment in the one-dimensional picture in position;
Image reconstruction module is reconstructed after handling the one-dimensional picture in orientation in data matrix different distance unit Real goal two dimensional image and false target two dimensional image afterwards.
7. the SAR false target Interference Suppression System according to claim 6 based on dynamic aperture, which is characterized in that institute Reference signal generation module is stated to specifically include:
Data correction unit, for obtaining the blended data of true SAR scene echoes and SAR cheating interference, carry out carrier wave, away from After tripping contracting, range migration correction processing, the data matrix of K × L dimension image scene is generated, s (t is denoted asr,ta), in data matrix Column direction indicates distance dimension, and line direction indicates azimuth dimension, trIt is distance to fast time, taFor the orientation slow time;
Reference signal generation unit is denoted as s for generating the reference signal of azimuth dimension deramp processing according to SAR system parameter0 (ta);Specifically:
Wherein fdcFor doppler centroid, γ is doppler frequency rate.
8. the SAR false target Interference Suppression System according to claim 7 based on dynamic aperture, which is characterized in that institute Aiming spot set generation module is stated to specifically include:
Frequency modulation removal signal generation unit, for obtaining data matrix s (tr,ta) in k-th of distance unit data sk(ta), with Reference signal s0(ta) it is point-by-point be multiplied after generate the signal s of frequency modulation removalDk(ta);
Wherein aa() refers to the orientation envelope of radar signal, is determined by the antenna radiation pattern of radar, and σ indicates the scattering of target point Coefficient, τ indicate the orientation delay of target point, RkK-th of distance unit is indicated to the distance of radar, λ is the wave of radar signal It is long;Orientation is one-dimensional as generation unit, for sDk(ta) Fourier transformation is carried out, generating includes real target point and false target The orientation of point is one-dimensional as Ik(fa);
WhereinIndicate Fourier transformation;
Aiming spot set generation unit, for obtaining Ik(fa) local peaking position, generate aiming spot to be estimated Set Df
9. the SAR false target Interference Suppression System according to claim 8 based on dynamic aperture, which is characterized in that institute Dynamic observing matrix generation module is stated to specifically include:
Dynamic aperture window acquiring unit is η for obtaining centerp, width is the dynamic aperture window w of Tp(ta)
Dynamic aperture imaging results generation unit, for according to wp(ta) to sDk(ta) extract, after extraction carry out frequency modulation removal at As processing, dynamic aperture imaging results I is generatedkp(fa)
Dynamic observing matrix generation unit, for dynamic aperture imaging results Ikp(fa) in aiming spot set DfIt is corresponding The data amplitude of position extracts, and is stored as the pth row that dynamic observes data matrix U, is obtained completely according to the data of each row Dynamic observing matrix U.
10. the SAR false target Interference Suppression System according to claim 9 based on dynamic aperture, which is characterized in that institute Image reconstruction module is stated to specifically include:
Loss coefficient computing unit, for calculating separately the amplitude loss factor alpha (η of real target pointpf) and false target point Amplitude loss factor beta (ηp)
Characterization matrix structural unit, for the center η according to dynamic aperture windowpAnd aiming spot set DfCorresponding τf, point It Ji Suan not α (ηpf) and β (ηp) observation characterization matrix A and B are constructed afterwards;
Minimization problem solves unit, for solving minimization problemThe true mesh separated The one-dimensional picture of the corresponding true bearing of the punctuate vector p one-dimensional picture of false bearing corresponding with false target point vector q;
Wherein p, q respectively indicate real target point vector and false target point vector to be solved, P=diag (p), Q=diag (q)。
CN201611040047.6A 2016-11-11 2016-11-11 A kind of SAR false target disturbance restraining method and system based on dynamic aperture Active CN106772273B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611040047.6A CN106772273B (en) 2016-11-11 2016-11-11 A kind of SAR false target disturbance restraining method and system based on dynamic aperture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611040047.6A CN106772273B (en) 2016-11-11 2016-11-11 A kind of SAR false target disturbance restraining method and system based on dynamic aperture

Publications (2)

Publication Number Publication Date
CN106772273A CN106772273A (en) 2017-05-31
CN106772273B true CN106772273B (en) 2019-03-05

Family

ID=58975013

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611040047.6A Active CN106772273B (en) 2016-11-11 2016-11-11 A kind of SAR false target disturbance restraining method and system based on dynamic aperture

Country Status (1)

Country Link
CN (1) CN106772273B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108318865B (en) * 2018-01-22 2022-05-03 西安电子科技大学 Multichannel SAR deception jamming identification and self-adaptive suppression method
CN109001732A (en) * 2018-06-07 2018-12-14 西北工业大学 A kind of compressed sensing Step Frequency SAR imaging restoration and reconstruction method of optimization
CN110850401B (en) * 2019-08-27 2022-06-28 天津大学 RFID label positioning method based on motion model and synthetic aperture
CN110852521B (en) * 2019-11-18 2022-06-07 中国人民解放军国防科技大学 Method and device for generating deception path
CN113203991B (en) * 2021-04-29 2022-05-31 电子科技大学 Anti-deception jamming method of multi-base SAR (synthetic aperture radar) in multi-jammer environment
CN113311431B (en) * 2021-05-27 2023-09-05 深圳大学 Deception jamming inhibition method and system based on single-channel SAR single imaging
CN114895264A (en) * 2022-06-23 2022-08-12 南京慧尔视智能科技有限公司 Method and device for suppressing false radar target, electronic equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102353945A (en) * 2011-03-31 2012-02-15 北京航空航天大学 Three-dimensional position reconstructing method based on ISAR (inverse synthetic aperture radar) image sequence for scattering point
CN105467368A (en) * 2015-11-25 2016-04-06 深圳大学 Multi-receiver equidistant rectangular distribution-based SAR deception jamming method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102353945A (en) * 2011-03-31 2012-02-15 北京航空航天大学 Three-dimensional position reconstructing method based on ISAR (inverse synthetic aperture radar) image sequence for scattering point
CN105467368A (en) * 2015-11-25 2016-04-06 深圳大学 Multi-receiver equidistant rectangular distribution-based SAR deception jamming method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
双通道极化SAR干扰抑制;郭睿 等;《电子与信息学报》;20100615;第32卷(第6期);第1343-1349页

Also Published As

Publication number Publication date
CN106772273A (en) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106772273B (en) A kind of SAR false target disturbance restraining method and system based on dynamic aperture
Ponce et al. Fully polarimetric high-resolution 3-D imaging with circular SAR at L-band
US9075129B2 (en) Method and system for forming images by comparing subsets of image data
McCorkle et al. Order N^ 2 log (N) backprojector algorithm for focusing wide-angle wide-bandwidth arbitrary-motion synthetic aperture radar
CN108051809B (en) Moving target imaging method and device based on Radon transformation and electronic equipment
US8125370B1 (en) Polarimetric synthetic aperture radar signature detector
Ash et al. Wide-angle synthetic aperture radar imaging: Models and algorithms for anisotropic scattering
US9250323B2 (en) Target detection utilizing image array comparison
Tang et al. Multipolarization through-wall radar imaging using low-rank and jointly-sparse representations
CN106680815B (en) MIMO radar imaging method based on tensor sparse representation
Zhou et al. Three-dimensional scattering center extraction based on wide aperture data at a single elevation
d'Alessandro et al. Phenomenology of ground scattering in a tropical forest through polarimetric synthetic aperture radar tomography
CN107576961A (en) A kind of relatively prime down-sampled sparse imaging method of interval synthetic aperture radar
CN112415515A (en) Method for separating targets with different heights by airborne circular track SAR
Ma et al. A sequential approach for Sentinel-1 TOPS time-series co-registration over low coherence scenarios
Karakuş et al. Generalized Bayesian model selection for speckle on remote sensing images
JP5134878B2 (en) Radar equipment
Garry et al. Passive ISAR part I: framework and considerations
Li et al. Azimuth super-resolution of forward-looking imaging based on Bayesian learning in complex scene
CN108845318B (en) Satellite-borne high-resolution wide-range imaging method based on Relax algorithm
Parker et al. Preliminary experimental results for RF tomography using distributed sensing
Bouzerdoum et al. A low-rank and jointly-sparse approach for multipolarization through-wall radar imaging
Huyghebaert et al. Determination of the azimuthal extent of coherent E-region scatter using the ICEBEAR linear receiver array
Pavlikov et al. Algorithm of Formation Radio Images from Aerospace Carriers
Sanjuan-Ferrer Detection of coherent scatterers in SAR data: Algorithms and applications

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant