CN110221265B - Distance extension target detection method based on strong scattering point self-adaptive estimation - Google Patents

Distance extension target detection method based on strong scattering point self-adaptive estimation Download PDF

Info

Publication number
CN110221265B
CN110221265B CN201910476624.3A CN201910476624A CN110221265B CN 110221265 B CN110221265 B CN 110221265B CN 201910476624 A CN201910476624 A CN 201910476624A CN 110221265 B CN110221265 B CN 110221265B
Authority
CN
China
Prior art keywords
threshold
gamma
value
eta
strong scattering
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
CN201910476624.3A
Other languages
Chinese (zh)
Other versions
CN110221265A (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.)
Xi'an Changyuan Electron Engineering Co ltd
Xian Electronic Engineering Research Institute
Original Assignee
Xi'an Changyuan Electron Engineering Co ltd
Xian Electronic Engineering Research Institute
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 Xi'an Changyuan Electron Engineering Co ltd, Xian Electronic Engineering Research Institute filed Critical Xi'an Changyuan Electron Engineering Co ltd
Priority to CN201910476624.3A priority Critical patent/CN110221265B/en
Publication of CN110221265A publication Critical patent/CN110221265A/en
Application granted granted Critical
Publication of CN110221265B publication Critical patent/CN110221265B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a distance extension target detection method based on self-adaptive estimation of strong scattering points, aiming at the problem that the detection performance of a traditional extension target detector is not stable when the information of the scattering points of a target is unknown. Compared with the traditional algorithm, the algorithm provided by the invention has higher robustness, and does not need any prior information.

Description

Distance extension target detection method based on strong scattering point self-adaptive estimation
Technical Field
The invention relates to the field of radar target detection, in particular to range expansion target detection, and provides a novel range expansion target detection method based on strong scattering point self-adaptive estimation.
Background
The range high-resolution radar can reduce the clutter power of each range unit and more accurately characterize the structural characteristics of the target, thereby beneficially affecting target detection, classification, identification and the like, and therefore, the range high-resolution radar is widely applied. When the distance resolution is far smaller than the target, the target energy is dispersed to the plurality of distance resolution units to form a distance extended target, and the detection performance is greatly lost by adopting the traditional point target constant false alarm detection algorithm, so that the constant false alarm detection algorithm for researching the distance extended target has great significance.
Two classic extended target detection algorithms proposed in the eighties of the last century are an energy accumulation detector (integrator) and a binary detector (M/N), respectively, wherein the energy accumulation detector accumulates all scattering points in a window to be detected, and has better detection performance when signal capacity is uniformly distributed in the detection window, and has larger detection loss when the scattering points of the target are sparse; the binary detector is an extension of a narrow-band constant false alarm detector, is easy to understand and realize, and has two disadvantages, namely, the number of strong scattering centers needing a priori target is large in limitation in practical application, and the detection performance is reduced by only using the information of the number of the scattering centers without using energy information during the second detection.
There are several researchers studying the optimal detector when the scattering point density is a priori. There is a document that proposes a Scattering point Density prior maximum likelihood estimation detector (SDD-GLRT), which can effectively detect extended targets with different sparsity/Scattering point densities, and is a statistically optimal detection, and does not fully utilize Scattering point information. In the literature, a double-Threshold Constant False Alarm detector (DT-CFAR) is proposed by using density and amplitude information of scattering points, and double-Threshold detection is performed by using double thresholds, and a Constant False Alarm is used for each detection, so that the detection performance of a sparse target is greatly improved, but the density information of prior scattering points is also required.
In recent years, many scholars have attempted to design extended target detectors that do not rely on a priori information on scattering points. Some researchers have proposed a distance extended Target detector (OS-RSTD) based on Order statistics, which first performs descending Order arrangement on the scattering point capability in the detection window, and then performs accumulation detection on all possible scattering point numbers exhaustively based on the sorted scattering points, until making a decision, and solves the problem that the traditional method depends on the Target scattering point information. However, in order to maintain constant false alarm, the false alarm rate of each detection channel is higher than that of the detector, and when the number of scattering points is small, the loss is large compared with that of a single-point detector. An improved dual-threshold maximum likelihood estimation detector (DT-GLRT) is proposed by the scholars, which does not use prior information of scattering points to estimate a first threshold, but uses maximum likelihood estimation and AIC criteria, and uses noise power as the first threshold, and the second threshold is calculated according to a false alarm rate and the number of strong scattering points. Compared with the traditional method, the performance of the method in different scattering point distribution environments is greatly improved, but the target scattering point information is not utilized in the first threshold calculation, so that the false alarm rate of the first detection is high, and the detection performance of the system is reduced because the false alarm point participates in the second detection.
In summary, the main problem of the conventional extended target detection method is that firstly, prior scattering point information is required, which cannot be satisfied in many application occasions; secondly, when the prior information is unknown, the detection performance of the traditional method is not stable when different scattering points are distributed.
Disclosure of Invention
Technical problem to be solved
The invention provides a distance extension target detection method based on strong scattering point self-adaptive estimation, aiming at the problems that the traditional extension target detection method needs scattering point prior information or the detection performance is not stable enough when different scattering points are distributed.
Technical scheme
A distance extension target detection method based on strong scattering point self-adaptive estimation is characterized by comprising the following steps:
step 1:the echo signal after square rate detection is input to the detector, and is recorded as Y = { Y = 1 ,y 2 ,...y J };
Step 2: performing k-means clustering with the cluster of 2 on Y to obtain two cluster sets, and marking the cluster with the large mean value as C 1 And another cluster is marked as C 2
And step 3: estimating the value of the number of strong scattering centers K: k = card (C) 1 ) Where card represents the number of elements in the matrix;
and 4, step 4: determining the value of a first threshold gamma as set C 2 The maximum value of the medium element, i.e.;
Figure BDA0002082478120000031
where max (, is taken to be the maximum value;
and 5: determining the value of the second threshold η:
the relationship among the number K of strong scattering points, the primary threshold gamma and the second threshold eta:
Figure BDA0002082478120000032
wherein, P fa Is the false alarm probability, σ, of the system 2 The power is noise power, and J is the number of input detection points;
will P fa 、σ 2 J and K are used as known quantities, a gamma value can be obtained by calculation for each eta, and the relation between gamma and eta is stored into a table for table lookup; inquiring corresponding eta according to the gamma obtained in the step 4 by adopting a lookup table method, and if the gamma in the table is not accurate, inquiring the eta corresponding to the closest gamma in the table to be used as a second threshold;
and 6: carrying out non-coherent accumulation on the strong scattering points and obtaining the quantity to be detected
Figure BDA0002082478120000033
And 7: and when D is larger than eta, judging that the target is present, otherwise, judging that the target is absent.
Advantageous effects
The invention provides a distance extended target detection method based on strong scattering point self-adaptive estimation, which can self-adaptively estimate the number of strong scattering points of a target without any prior information, and solves the problem that the traditional extended target method needs the prior scattering point information; the method has high accuracy in estimating the strong scattering points, so that the detection performance is stable when different scattering points are distributed, and the problem that the detection performance is not stable when the traditional method is distributed at different scattering points is solved.
Compared with the traditional group target identification method, the method has the advantages that:
1) In the threshold design, not only noise information but also target information are utilized, and the strong scattering point can be estimated in a self-adaptive manner, so that the detection performance under different scattering point distribution environments is improved;
2) The prior information of scattering points is not relied on, so that the application range of the scattering point is greatly expanded;
3) The number of the threshold and the strong scattering points is estimated by adopting a kmeans clustering algorithm, the algorithm operation amount is small, and the engineering realization is easy.
Drawings
FIG. 1 is a flow chart of an extended target detection method proposed by the present invention
Detailed Description
The invention will now be further described with reference to the following examples and drawings:
the method adopts a double-threshold idea, and firstly self-adaptively estimates the number of strong scattering points and a first threshold by using a kmeans clustering algorithm; then, a second threshold is determined according to the first threshold and the number of scattering points, so that the constant false alarm performance is ensured; and finally, completing detection through two judgments. The method comprises the following specific steps:
1) The input to the detector is an echo signal after square rate detection, and is denoted as Y = { Y = 1 ,y 2 ,...y J };
2) Performing k-means clustering (cluster is 2) on Y to obtain two cluster sets, and marking the cluster with large mean value as C 1 And another cluster is marked as C 2
3) Estimate of the number of strong scattering centers K, K = card (C) 1 ) Card stands for momentThe number of elements in the array;
4) The value of the first threshold gamma is set C 2 The maximum value of the medium element, i.e.;
Figure BDA0002082478120000041
where max (, is taken to be the maximum value;
5) The formula (1) is a relation among the number K of strong scattering points, a primary threshold gamma and a second threshold eta:
Figure BDA0002082478120000042
in the formula, P fa The false alarm probability of the system is preset. Sigma 2 The noise power can be estimated in real time, J is the number of input detection points, K is the number of strong scattering points estimated in step 3), and γ is the value of the first threshold estimated in step 4), so that the variable in the formula (1) only has a quadratic threshold η. But since eta cannot contain known quantity P fa 、σ 2 J, K and gamma, and using the method of lookup table to express P fa 、σ 2 J and K are known quantities, a gamma value can be calculated for each eta, and the relation between gamma and eta is stored into a table for table lookup. In the table lookup process, for a given P fa 、σ 2 And J, manufacturing a plurality of tables according to different K values for different K values estimated in the step 3). When the gamma value corresponding to which eta cannot be accurately corresponded, the estimated first threshold is used
Figure BDA0002082478120000051
Closest to the gamma value in the table
Figure BDA0002082478120000053
As a second threshold η;
6) Carrying out non-coherent accumulation on the strong scattering points and obtaining the quantity to be detected
Figure BDA0002082478120000052
7) And finishing detection, judging that the target is present if D is larger than eta, and otherwise, judging that the target is absent.

Claims (1)

1. A distance extension target detection method based on strong scattering point self-adaptive estimation is characterized by comprising the following steps:
step 1: the echo signal after square rate detection is input to the detector, and is recorded as Y = { Y = 1 ,y 2 ,...y J };
And 2, step: performing k-means clustering with the cluster of 2 on Y to obtain two cluster sets, and recording the cluster with large mean value as C 1 And another cluster is marked as C 2
And 3, step 3: estimating the value of the number of strong scattering centers K: k = card (C) 1 ) Where card represents the number of elements in the matrix;
and 4, step 4: determining the value of a first threshold gamma as set C 2 The maximum value of the medium element, i.e.;
Figure FDA0003877473720000011
where max (, is taken to be the maximum value;
and 5: determining a value of a second threshold η:
the relationship among the number K of strong scattering points, the first threshold gamma and the second threshold eta:
Figure FDA0003877473720000012
wherein, P fa Is the false alarm probability, σ, of the system 2 The power is noise power, and J is the number of input detection points;
will P fa 、σ 2 J and K are used as known quantities, a gamma value can be obtained by calculation for each eta, and the relation between gamma and eta is stored into a table for table lookup; when the gamma value corresponding to which eta can not be accurately corresponded, the estimated first threshold is used
Figure FDA0003877473720000013
With gamma in the tableThe values being closest
Figure FDA0003877473720000014
As a second threshold η;
step 6: carrying out non-coherent accumulation on the strong scattering points and obtaining the quantity to be detected
Figure FDA0003877473720000015
And 7: and if D is larger than eta, judging that the target is present, otherwise, judging that the target is absent.
CN201910476624.3A 2019-06-03 2019-06-03 Distance extension target detection method based on strong scattering point self-adaptive estimation Active CN110221265B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910476624.3A CN110221265B (en) 2019-06-03 2019-06-03 Distance extension target detection method based on strong scattering point self-adaptive estimation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910476624.3A CN110221265B (en) 2019-06-03 2019-06-03 Distance extension target detection method based on strong scattering point self-adaptive estimation

Publications (2)

Publication Number Publication Date
CN110221265A CN110221265A (en) 2019-09-10
CN110221265B true CN110221265B (en) 2022-12-06

Family

ID=67819165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910476624.3A Active CN110221265B (en) 2019-06-03 2019-06-03 Distance extension target detection method based on strong scattering point self-adaptive estimation

Country Status (1)

Country Link
CN (1) CN110221265B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115453484B (en) * 2022-08-29 2023-08-25 电子科技大学 Self-adaptive scattering center estimated distance expansion target detection method
CN116299401B (en) * 2023-05-19 2023-10-17 成都航空职业技术学院 Constant false alarm method and device based on target scattering point position and storage medium thereof
CN116736256B (en) * 2023-08-11 2023-10-24 南京隼眼电子科技有限公司 Radar identification method and device and electronic equipment

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5703592A (en) * 1995-01-17 1997-12-30 Racal Radar Defence Systems Limited Method and apparatus for estimating the detection range of a radar
JP2002236172A (en) * 2000-12-06 2002-08-23 Mitsubishi Electric Corp Radar device
CN102230961A (en) * 2011-03-31 2011-11-02 北京航空航天大学 Method for detecting spread target based on phase compensation processing
CN102323577A (en) * 2011-09-08 2012-01-18 北京理工雷科电子信息技术有限公司 High-resolution radar dual-threshold detector based on energy accumulation
CN104680183A (en) * 2015-03-14 2015-06-03 西安电子科技大学 SAR target identification method based on scattering point and K-center one-class classifier
CN107390196A (en) * 2017-07-28 2017-11-24 西安电子科技大学 Fast fluctuating target double threshold CFAR detection method based on multistatic radar
CN107861107A (en) * 2017-10-23 2018-03-30 电子科技大学 A kind of double threshold CFAR suitable for continuous wave radar and Plot coherence method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5703592A (en) * 1995-01-17 1997-12-30 Racal Radar Defence Systems Limited Method and apparatus for estimating the detection range of a radar
JP2002236172A (en) * 2000-12-06 2002-08-23 Mitsubishi Electric Corp Radar device
CN102230961A (en) * 2011-03-31 2011-11-02 北京航空航天大学 Method for detecting spread target based on phase compensation processing
CN102323577A (en) * 2011-09-08 2012-01-18 北京理工雷科电子信息技术有限公司 High-resolution radar dual-threshold detector based on energy accumulation
CN104680183A (en) * 2015-03-14 2015-06-03 西安电子科技大学 SAR target identification method based on scattering point and K-center one-class classifier
CN107390196A (en) * 2017-07-28 2017-11-24 西安电子科技大学 Fast fluctuating target double threshold CFAR detection method based on multistatic radar
CN107861107A (en) * 2017-10-23 2018-03-30 电子科技大学 A kind of double threshold CFAR suitable for continuous wave radar and Plot coherence method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Distributed Target Detection Exploiting Persymmetry in Gaussian Clutter;Jun Liu Et al.;《IEEE Transactions on Signal Processing》;20190215;第67卷(第4期);第1022-1033页 *
基于强散射点径向积累的高分辨极化目标检测研究;李永祯等;《电子学报》;20010331;第29卷(第3期);第307-310页 *
距离扩展目标的双门限恒虚警检测器及性能分析;顾新锋等;《电子与信息学报》;20120630;第34卷(第6期);第1318-1323页 *

Also Published As

Publication number Publication date
CN110221265A (en) 2019-09-10

Similar Documents

Publication Publication Date Title
CN110221265B (en) Distance extension target detection method based on strong scattering point self-adaptive estimation
CN103995259B (en) Radar target adaptive-filtering fusion detection method under intensive interference environment
CN113721219B (en) Radar signal sorting method and system based on multi-parameter clustering
CN116166960B (en) Big data characteristic cleaning method and system for neural network training
CN104462819A (en) Local outlier detection method based on density clustering
CN111275132A (en) Target clustering method based on SA-PFCM + + algorithm
CN109145114B (en) Social network event detection method based on Kleinberg online state machine
CN102243641A (en) Method for efficiently clustering massive data
CN114036610A (en) Penetration depth prediction method based on data enhancement
CN111611293B (en) Outlier data mining method based on feature weighting and MapReduce
CN104751171A (en) Method of classifying Naive Bayes scanned certificate images based on feature weighting
CN105095908A (en) Video image group behavior characteristic processing method and apparatus
CN112799028B (en) False target identification method based on RCS fluctuation statistical characteristic difference
CN110988856B (en) Target detection trace agglomeration algorithm based on density clustering
CN114355298A (en) Radar composite modulation pulse signal identification method
US20220114460A1 (en) Apparatus of Identifying Heterogeneous Time-Series Data Expression with High Efficiency
CN114519372B (en) One-dimensional range profile target recognition method based on support vector machine
Dai et al. Clustering of DOA data in radar pulse based on SOFM and CDbw
CN112014821B (en) Unknown vehicle target identification method based on radar broadband characteristics
Zhang et al. Maximum focal inter-class angular loss with norm constraint for automatic modulation classification
CN110488279B (en) Pedestrian and tree distinguishing method and device based on radar
CN115453484B (en) Self-adaptive scattering center estimated distance expansion target detection method
CN113466798B (en) Sea surface small target detection method based on polarized signal multi-scale entropy characteristics
CN110784887A (en) Method for detecting number of abnormal signal sources in gridding radio signal monitoring system
Zhi et al. A Self-Adaptive OPTICS Clustering Algorithm Based on the Lightning Distribution

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
TA01 Transfer of patent application right

Effective date of registration: 20221026

Address after: 710199 Weiqu Fengqi East Street, Chang'an District, Xi'an City, Shaanxi Province

Applicant after: Xi'an Electronic Engineering Research Institute

Applicant after: XI'AN CHANGYUAN ELECTRON ENGINEERING Co.,Ltd.

Address before: Changan District Fengqi road 710100 Shaanxi city of Xi'an Province

Applicant before: Xi'an Electronic Engineering Research Institute

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant