CN104914422A - Adaptive TBD radar weak target detection method - Google Patents
Adaptive TBD radar weak target detection method Download PDFInfo
- 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
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/41—Details 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/414—Discriminating 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
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:
Thus can to obtain one-level thresholding be η=-ln (p
f), each resolution element detection probability p
dfor:
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:
Carry out approximate processing to the Section 2 of above formula, finally obtaining detection probability is:
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:
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.
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)
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)
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 |
-
2015
- 2015-06-25 CN CN201510359853.9A patent/CN104914422A/en active Pending
Patent Citations (2)
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)
Title |
---|
邓晓波 等: "基于TBD方法的弱小目标预警研究", 《现代雷达》 * |
郝有东 等: "一种与杂波类型无关的延展目标CFAR检测器", 《空军雷达学院学报》 * |
黄勇 等: "基于Hough变换的检测前跟踪算法的性能分析", 《现代雷达》 * |
Cited By (16)
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 |