CN104914422A - Adaptive TBD radar weak target detection method - Google Patents

Adaptive TBD radar weak target detection method Download PDF

Info

Publication number
CN104914422A
CN104914422A CN201510359853.9A CN201510359853A CN104914422A CN 104914422 A CN104914422 A CN 104914422A CN 201510359853 A CN201510359853 A CN 201510359853A CN 104914422 A CN104914422 A CN 104914422A
Authority
CN
China
Prior art keywords
noise ratio
signal
ratio
threshold
accumulation
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.)
Pending
Application number
CN201510359853.9A
Other languages
Chinese (zh)
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.)
724th Research Institute of CSIC
Original Assignee
724th Research Institute of CSIC
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 724th Research Institute of CSIC filed Critical 724th Research Institute of CSIC
Priority to CN201510359853.9A priority Critical patent/CN104914422A/en
Publication of CN104914422A publication Critical patent/CN104914422A/en
Pending legal-status Critical Current

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
    • 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/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter

Abstract

The present invention relates to an adaptive TBD radar weak target detection method which is mainly applied to the detection of a heavy clutter zone weak group target in a sea guard radar. For an original data, according to the needed detection probability and false alarm rate, a first threshold is set, points which are lower than the threshold are abandoned, and the coordinates and amplitudes of points which are higher than the threshold are recorded. After 6 frames of data are accumulated, all recorded points are subjected to Hough transform, in the process of carrying out target detection in a transform domain, the accumulated value is determined according to the signal to noise ratio of an echo, energy accumulation is employed for a low signal to noise ratio, binary accumulation is employed for a high signal to noise ratio, amplitude accumulation is employed for a middle signal to noise ratio, a second threshold is set, the target is announced when the threshold is exceeded, a related point is found through inverse transformation, at the same time the track of the target is announced, the binary accumulation is employed for the high signal to noise ratio, the energy accumulation is employed for the low signal to noise ratio, and the amplitude accumulation is employed for the middle signal to noise ratio.

Description

A kind of self-adaptation TBD radar Weak target detecting method
Technical field
The invention belongs to sea guard radar object detection field, particularly heavy clutter district multiple targets based on tracking technique before the detection of Hough transform.
Background technology
Dim targets detection under clutter environment is the significant challenge that modern sensor faces, traditional detection technique detects every frame signal, and by not having the data associated to lose, owing to not accumulating the echo information of interframe, be difficult to improve the detection perform to Small object.And based on following the tracks of (TBD) technology before the detection of Hough transform, achieving the interframe non-inherent accumulation of target, can detection perform be improved.
1994, Hough transform is applied in search radar target detection by the people such as the Carlson of the U.S. first, propose the TBD algorithm based on Hough transform, and analyze the detection and tracking performance of this algorithm, give the computing formula of detection probability and false-alarm probability.Compared with classic method, the method has lot of advantages, as improve detection probability, solving range walk problem, being easier to realize; By emulation, for dissimilar target, signal to noise ratio (S/N ratio) can improve 4dB, but there is the problem that Small object covered by general objective, improved afterwards based in the binary integration device of HT, although can address this problem, performance declines to some extent, have the SNR of 3dB to improve for non-fluctuating target, strong fluctuating target has the SNR of 1dB to improve.
On the basis of Carlson work, the scholars such as Kabakchiev discuss in different clutter background situation, the detection and tracking performance of the combined detector that CFAR detects and Hough transform is formed.But above-mentioned algorithm all hypothetical target does linear uniform motion along radar radial direction, and within multiple antenna cycle, target echo does not cross over localizer unit.For the situation that target leap localizer unit is walked about, Garvanov proposes a kind of polar coordinates Hough transform TBD algorithm, directly processes the radar measurement data under azimuth-range dimension polar coordinates, and this algorithm is more suitable for process radar measured data.2011, the Moyer of the U.S. proposed a kind of multidimensional HT-TBD algorithm for weak target under strong clutter background, and this algorithm can overcome the defect of traditional HT-TBD algorithm, strengthened the target detection probability under random noise background.
However, when adopting the TBD algorithm based on Hough transform to detect multiple targets, if the RCS difference between target is larger, large target can produce Small object covers, if reduction detection threshold, improve the detection probability to Small object, then can comprise more false-alarm in testing result.
Summary of the invention
On the basis of the present invention's track algorithm before heavy clutter district is based on the detection of the multiple targets of Hough transform, propose a kind of adaptive detection method.The amplitude range of target echo is larger, under different signal to noise ratio (S/N ratio) conditions, adopt different detection methods to improve detection probability, reduces false alarm rate.Traditional based on track algorithm before the detection of Hough transform, the echoed signal of time domain is transformed to parameter field, accumulates at parameter field, adopt the weights of single accumulation, be difficult in the clutter situation of change, obtain best detection perform.Herein by the adaptive mode selecting accumulation, effectively solve target occlusion and the low problem of detection probability.
Accompanying drawing explanation
Accompanying drawing 1 is trace flow before the self-adapting detecting based on Hough transform.
Embodiment
Tracking before the adaptive detection of clutter district multiple targets, concrete implementation step following (see accompanying drawing 1):
1. the echo data of continuous 6 frames of sequential storage surveyed area.
Hough transform will carry out multiframe data processing, and then target echo accumulation is inadequate very little for frame number, and target detection probability is low.If frame number is too many, then cause the problems such as memory space is large, calculated amount is large, false-alarm is high.Choosing comprehensively selects 6 frame data Combined Treatment.
2. calculate the signal to noise ratio (S/N ratio)/signal to noise ratio of echo.
The computing formula of the signal to noise ratio (S/N ratio) of echo is SNR=E/N 0, wherein E is the power of target echo, N 0it is Carrier To Noise Power Density.
3. filter raw radar data.
If all carry out Hough transform to every frame data of luv space, calculated amount will be very huge.In fact raw radar data contains a large amount of noises, adopts one-level thresholding to filter luv space data, only carries out Hough transform to the data crossing thresholding.Calculated amount increases with one-level thresholding and declines.Lower one-level thresholding, owing to having processed more data, improve detection perform, but calculated amount is large.
In distance to a declared goal-time domain, the initial false-alarm probability of each resolution element is p f, then:
p f = ∫ η ∞ p ( x ) d x = ∫ η ∞ e - x d x = e - η
Thus can to obtain one-level thresholding be η=-ln (p f), each resolution element detection probability p dfor:
p d = ∫ η ∞ p ( x | S ) d x = ∫ η ∞ e - x / ( 1 + S ) 1 + S d x = e - η / ( 1 + S )
In formula, S represents signal to noise ratio (S/N ratio), as S=0, and can from p dobtain p f.
4. adaptive H ough converts
By the data Hough transform of η of one-level thresholding excessively all in distance verses time territory to parameter field, when parameter field is weighted accumulation, according to the signal to noise ratio (S/N ratio) that step 2 calculates, according to the selection weighting accumulative means in the interval at signal to noise ratio (S/N ratio) place.By the detection perform curve of Monte-Carlo Simulation three kinds of accumulative means along with signal to noise ratio (S/N ratio), signal to noise ratio, between the application area determining three kinds of accumulative means.Signal amplitude accumulates: using the echo amplitude of each point as accumulation item; Binary integration: having crossed one-level thresholding is then 1, but is then 0; Energy accumulation: accumulate again after the amplitude square of signal.
5. detection probability and false alarm rate calculate
Suppose that the data in distance verses time territory have N number of unit to parameter field contributive, wherein have the data (x of m unit i, i=1 ..., probability m) having exceeded initial threshold η is and this m unit sum must be over the second thresholding ξ, namely in parameter field, a certain unit ζ more than the probability of the second thresholding ξ is:
Pr ( ζ > ξ ) = Σ m = 1 N C N m ( p d ) m ( 1 - p d ) N - m × Pr ( y = Σ i = 1 m x i > ξ )
Carry out approximate processing to the Section 2 of above formula, finally obtaining detection probability is:
P D = 1 - ( 1 - p d ) N - Σ m = 1 m * - 1 C N m ( p d ) m ( 1 - p d ) N - m × [ 1 - e - ξ - m η 1 + S Σ k = 0 m - 1 ( ξ - m η 1 + S ) k / k ! ]
M in formula *=ceil (ξ/η) is for being not less than the smallest positive integral of ξ/η. and when S is 0, can obtain false-alarm probability is:
P F = 1 - ( 1 - p f ) N - Σ m = 1 m * - 1 C N m ( p f ) m ( 1 - p f ) N - m × [ 1 - e - ( ξ - m η ) Σ k = 0 m - 1 ( ξ - m η ) k / k ! ]
6. the parameter corresponding according to check point, the track of backtracking target, announces testing result.
For the ease of data processing, by parameter space discretize, form grid (unit) one by one, owing to there being the coordinate points of Cartesian coordinates that the Curves of common intersection is corresponding must be point-blank in parameter space, second thresholding ξ is set in parameter field again, as final detection threshold.By inverse Hough transform, check point mapped time data space of parameter space, the targetpath of formation time-distance plane, can realize the detection to target.

Claims (3)

1. a self-adaptation TBD radar Weak target detecting method, it is characterized in that: testing result is not announced to single frames target, by Hough transform Combined Treatment multiframe data, in the mode of transform domain according to the adaptive selection weighting accumulation of the signal to noise ratio (S/N ratio) of echo: high s/n ratio adopts binary weighting accumulation, middle signal to noise ratio (S/N ratio) adopts amplitude weighting accumulation, and low signal-to-noise ratio adopts the weighting accumulation of energy.
2. self-adaptation TBD radar Weak target detecting method according to claim 1, is characterized in that: the determination methods in the height of signal to noise ratio (S/N ratio), for determining two threshold values, is greater than high threshold, is judged to be high s/n ratio situation; Be less than threshold ones and be judged to be low signal-to-noise ratio; Be in the middle of two threshold values is middle signal to noise ratio (S/N ratio).
3. self-adaptation TBD radar Weak target detecting method according to claim 1 and 2, it is characterized in that the defining method of described threshold value is: in the clutter meeting K distribution, add two batches of varying strength targets, detection probability under three kinds of weighting accumulative means is drawn by Monte-Carlo Simulation, false alarm rate, with the change curve of signal to noise ratio (S/N ratio), determines the threshold value of two signal to noise ratio (S/N ratio)s by the distribution of curve.
CN201510359853.9A 2015-06-25 2015-06-25 Adaptive TBD radar weak target detection method Pending CN104914422A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510359853.9A CN104914422A (en) 2015-06-25 2015-06-25 Adaptive TBD radar weak target detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510359853.9A CN104914422A (en) 2015-06-25 2015-06-25 Adaptive TBD radar weak target detection method

Publications (1)

Publication Number Publication Date
CN104914422A true CN104914422A (en) 2015-09-16

Family

ID=54083654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510359853.9A Pending CN104914422A (en) 2015-06-25 2015-06-25 Adaptive TBD radar weak target detection method

Country Status (1)

Country Link
CN (1) CN104914422A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911542A (en) * 2016-07-04 2016-08-31 中国人民解放军海军航空工程学院 Hypersonic velocity target TBD detection method for polynomial Hough conversion
CN106093904A (en) * 2016-06-17 2016-11-09 电子科技大学 Clutter map CFAR Methods based on multiframe double threshold hierarchical detection mechanism
CN106651878A (en) * 2016-12-21 2017-05-10 福建师范大学 Method for extracting straight line from local invariant feature points
CN109143184A (en) * 2018-10-29 2019-01-04 北京理工大学 A kind of double threshold detection method of scanning radar
CN109669180A (en) * 2018-12-14 2019-04-23 河海大学 A kind of continuous wave radar unmanned plane detection method
CN109839621A (en) * 2017-11-24 2019-06-04 西安艾索信息技术有限公司 A kind of improved TBD algorithm
CN109917375A (en) * 2019-03-11 2019-06-21 西安电子工程研究所 Low repetition it is short it is resident under the conditions of Hovering Helicopter detection method
CN112327297A (en) * 2020-10-23 2021-02-05 北京理工大学 Radar detection correlation method for fluctuation characteristics of insects
CN112526501A (en) * 2019-09-19 2021-03-19 苏州豪米波技术有限公司 Radar system for detecting life breath
CN113075636A (en) * 2021-04-02 2021-07-06 中国人民解放军海军航空大学 Parallel line coordinate transformation and weak target detection method for measuring points
CN113702965A (en) * 2021-08-31 2021-11-26 中国人民解放军海军航空大学 Improved accumulation method based on peak value convergence and simultaneous detection method of strong and weak targets

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412295A (en) * 2013-08-30 2013-11-27 西安电子科技大学 High-speed maneuvering weak target detection method based on echo precise model
CN103558597A (en) * 2013-11-15 2014-02-05 武汉大学 Method for detecting dim target in sea clutter on basis of spectrum kurtosis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412295A (en) * 2013-08-30 2013-11-27 西安电子科技大学 High-speed maneuvering weak target detection method based on echo precise model
CN103558597A (en) * 2013-11-15 2014-02-05 武汉大学 Method for detecting dim target in sea clutter on basis of spectrum kurtosis

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
邓晓波 等: "基于TBD方法的弱小目标预警研究", 《现代雷达》 *
郝有东 等: "一种与杂波类型无关的延展目标CFAR检测器", 《空军雷达学院学报》 *
黄勇 等: "基于Hough变换的检测前跟踪算法的性能分析", 《现代雷达》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093904A (en) * 2016-06-17 2016-11-09 电子科技大学 Clutter map CFAR Methods based on multiframe double threshold hierarchical detection mechanism
CN105911542A (en) * 2016-07-04 2016-08-31 中国人民解放军海军航空工程学院 Hypersonic velocity target TBD detection method for polynomial Hough conversion
CN106651878B (en) * 2016-12-21 2019-06-11 福建师范大学 A method of for extracting straight line from local invariant feature point
CN106651878A (en) * 2016-12-21 2017-05-10 福建师范大学 Method for extracting straight line from local invariant feature points
CN109839621A (en) * 2017-11-24 2019-06-04 西安艾索信息技术有限公司 A kind of improved TBD algorithm
CN109143184A (en) * 2018-10-29 2019-01-04 北京理工大学 A kind of double threshold detection method of scanning radar
CN109669180B (en) * 2018-12-14 2020-09-15 河海大学 Continuous wave radar unmanned aerial vehicle detection method
CN109669180A (en) * 2018-12-14 2019-04-23 河海大学 A kind of continuous wave radar unmanned plane detection method
CN109917375A (en) * 2019-03-11 2019-06-21 西安电子工程研究所 Low repetition it is short it is resident under the conditions of Hovering Helicopter detection method
CN109917375B (en) * 2019-03-11 2023-01-03 西安电子工程研究所 Method for detecting hovering helicopter under condition of low repetition frequency and short residence
CN112526501A (en) * 2019-09-19 2021-03-19 苏州豪米波技术有限公司 Radar system for detecting life breath
CN112327297A (en) * 2020-10-23 2021-02-05 北京理工大学 Radar detection correlation method for fluctuation characteristics of insects
CN112327297B (en) * 2020-10-23 2024-04-12 北京理工大学 Radar detection correlation method for insect gamma fluctuation characteristics
CN113075636A (en) * 2021-04-02 2021-07-06 中国人民解放军海军航空大学 Parallel line coordinate transformation and weak target detection method for measuring points
CN113075636B (en) * 2021-04-02 2022-06-24 中国人民解放军海军航空大学 Parallel line coordinate transformation and weak target detection method for measuring points
CN113702965A (en) * 2021-08-31 2021-11-26 中国人民解放军海军航空大学 Improved accumulation method based on peak value convergence and simultaneous detection method of strong and weak targets

Similar Documents

Publication Publication Date Title
CN104914422A (en) Adaptive TBD radar weak target detection method
CN104215951B (en) System and method for detecting low-speed small target under sea cluster background
CN102879766B (en) Method and apparatus for detecting and tracking faint target of high frequency ground wave radar
CN103217673B (en) CFAR detecting method under inhomogeneous Weibull clutter background
CN104931934B (en) A kind of radar plot condensing method based on PAM cluster analyses
CN101452075B (en) At-sea small target detecting method based on average period
CN104316914B (en) Radar target self-adaptation detection method depending on shape parameters
CN107015221B (en) A kind of low false alarm rate fast target detection method for ground surveillance radar
CN104914433A (en) Linked list sorting-based OS-CFAR multi-target extraction realization method
CN104881561A (en) Hough transform-based track-before-detect method of multidimensional parameters
CN105319537A (en) Navigation radar co-frequency interference inhibition method based on spatial correlation
CN106468771A (en) A kind of multi-target detection and tracking method under the high clutter conditions of low Observable
CN106054150A (en) First-initiation second-confirmation radar track initiation method
CN104714215A (en) Clutter intensity zoning based constant-false-alarm detection method
CN106569193A (en) Sea-surface small target detection method based on front-back revenue reference particle filter
CN109342828A (en) A kind of lightening pulse signal detecting method based on frequency domain constant false alarm
CN102546052B (en) Self-adaptive threshold controlling and parameter calculating and measuring technology and method
CN109839621A (en) A kind of improved TBD algorithm
CN106353743A (en) Asymptotically optimal radar target detection method matched to equivalent shape parameter
CN102621535B (en) High-efficiency method for estimating covariance matrix structures
CN104931951A (en) Detection method based on Hough transform domain clutter map for small target in heavy clutter region
CN103605120A (en) Range-extended target intrinsic mode local singular value orderly fusion detection method
CN104898097A (en) FPGA-based phase demodulation constant false alarm rate (CFAR) radar signal detection method
CN107122764A (en) ShipTargets detection method based on KpN models
CN105242272A (en) Vehicle-mounted millimeter wave crashproof radar constant false alarm detection method based on autoregression time sequence model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150916

WD01 Invention patent application deemed withdrawn after publication