CN106918807B - A kind of Targets Dots condensing method of radar return data - Google Patents
A kind of Targets Dots condensing method of radar return data Download PDFInfo
- Publication number
- CN106918807B CN106918807B CN201710111371.0A CN201710111371A CN106918807B CN 106918807 B CN106918807 B CN 106918807B CN 201710111371 A CN201710111371 A CN 201710111371A CN 106918807 B CN106918807 B CN 106918807B
- Authority
- CN
- China
- Prior art keywords
- echo data
- target
- binary
- sliding window
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000001514 detection method Methods 0.000 claims abstract description 96
- 238000012545 processing Methods 0.000 claims abstract description 43
- 238000009833 condensation Methods 0.000 claims abstract description 37
- 230000005494 condensation Effects 0.000 claims abstract description 37
- 238000009825 accumulation Methods 0.000 claims abstract description 31
- 230000001427 coherent effect Effects 0.000 claims abstract description 12
- 230000000916 dilatatory effect Effects 0.000 claims abstract description 3
- 230000008569 process Effects 0.000 claims description 13
- 238000013139 quantization Methods 0.000 claims description 11
- 230000002776 aggregation Effects 0.000 claims description 10
- 238000005054 agglomeration Methods 0.000 claims description 8
- 230000003628 erosive effect Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000000717 retained effect Effects 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 19
- 238000001914 filtration Methods 0.000 description 10
- 238000005260 corrosion Methods 0.000 description 9
- 230000007797 corrosion Effects 0.000 description 9
- 230000010339 dilation Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006740 morphological transformation Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达信号处理技术领域,公开了一种雷达回波数据的目标点迹凝聚方法,包括:获取雷达多个扫描周期的原始回波数据;对原始回波数据依次进行包络检波和多脉冲非相参积累,得到脉冲积累后的回波数据;对脉冲积累后的回波数据进行恒虚警检测,进而对恒虚警检测后的回波数据依次进行二值量化处理以及沿方位维的二值滑窗检测;将二值滑窗检测后的回波数据作为二值图像数据,对二值图像数据进行先膨胀后腐蚀的操作,从而得到二值图像数据中的所有连通域;滤除第一连通域和第二连通域,得到包含目标点迹的剩余连通域;根据包含目标点迹的剩余连通域得到目标信息,实现对目标点迹的凝聚;能够提高目标参数估计精度和目标分辨能力。
The invention belongs to the technical field of radar signal processing, and discloses a method for condensing target point traces of radar echo data. The pulse non-coherent accumulation is used to obtain the echo data after pulse accumulation; the echo data after pulse accumulation is subjected to constant false alarm detection; binary sliding window detection; take the echo data detected by the binary sliding window as the binary image data, and perform the operation of dilating and then corroding the binary image data, so as to obtain all connected domains in the binary image data; filter In addition to the first connected domain and the second connected domain, the remaining connected domain containing the target point trace is obtained; the target information is obtained according to the remaining connected domain containing the target point trace, so as to realize the condensation of the target point trace; it can improve the estimation accuracy of the target parameters and the target Distinguishing ability.
Description
技术领域technical field
本发明属于雷达信号处理技术领域,尤其涉及一种雷达回波数据的目标点迹凝聚方法,用于雷达回波数据的预处理和目标参数提取。The invention belongs to the technical field of radar signal processing, and in particular relates to a method for condensing target point traces of radar echo data, which is used for preprocessing of radar echo data and extraction of target parameters.
背景技术Background technique
随着雷达技术的发展,人们要求雷达具有更多的功能和作用。传统意义上的雷达主要用于测量目标的距离、方位、俯仰等基本信息,现代雷达要求在更为复杂的环境(干扰与杂波)下实现对目标的检测、点迹录取和目标跟踪。With the development of radar technology, people require radar to have more functions and roles. The radar in the traditional sense is mainly used to measure the basic information such as distance, azimuth, and pitch of the target. Modern radar requires detection, point trace recording and target tracking to be realized in a more complex environment (interference and clutter).
现代雷达技术向着智能化、小型化、信息化和精密化方向发展。人们要求雷达能够获取更为精确的目标信息,并且实现更为稳定的目标跟踪。点迹凝聚技术就是通过对雷达回波数据进行处理实现目标参数估计和点迹录取的过程,雷达数据处理系统对录取的点迹进行航迹处理实现目标的跟踪显示。国内对雷达数据处理技术的研究相比国外起步较晚,点迹凝聚处理技术的研究主要针对MTD体制的雷达系统。现阶段点迹凝聚主要采用的方法有基于滑窗法的点迹凝聚、基于图像处理的点迹凝聚等。Modern radar technology is developing in the direction of intelligence, miniaturization, informatization and precision. People require radar to be able to obtain more precise target information and achieve more stable target tracking. Point trace condensation technology is the process of estimating target parameters and recording point traces by processing radar echo data. The radar data processing system processes the recorded point traces to achieve target tracking display. Compared with foreign countries, the domestic research on radar data processing technology started later, and the research on point trace condensation processing technology is mainly aimed at the radar system of the MTD system. At this stage, the main methods of point trace condensation are point trace condensation based on sliding window method and point trace condensation based on image processing.
点迹凝聚现阶段主要存在的问题是复杂环境下如何对目标点迹的准确判别,以及如何改进不同雷达体制下点迹凝聚处理的方法。基于滑窗法的点迹凝聚方法是较为常用的一种点迹凝聚方法,主要存在的问题是目标终了门限值较大时会出现目标点迹的分裂,目标在方位上的分裂常常通过降低门限值的方法进行解决,而门限值选取又不能过低,过低的门限值造成虚假点迹过多,影响点迹质量。经典的基于图像处理的方法是采用图像轮廓查找求取质心的方法实现点迹凝聚,无法改善点迹在距离上和方位上的分裂问题。The main problems of point trace condensation at this stage are how to accurately discriminate target point traces in complex environments, and how to improve the processing methods of point trace condensation under different radar systems. The point trace condensation method based on the sliding window method is a commonly used point trace condensation method. The main problem is that the target point trace will be split when the target end threshold is large. The method of threshold value is used to solve the problem, and the threshold value should not be too low. Too low threshold value will cause too many false traces and affect the quality of traces. The classical method based on image processing is to use the method of image contour search to obtain the centroid to achieve point trace condensation, which cannot improve the splitting problem of point traces in distance and azimuth.
发明内容SUMMARY OF THE INVENTION
针对上述现有技术的缺点,本发明的目的在于提供一种雷达回波数据的目标点迹凝聚方法,不仅采用滑窗法中的二值滑窗检测,同时采用了数学中的形态滤波算法对连通域进行膨胀腐蚀操作,通过选取合适的结构元素不仅可以改善滑窗法出现的目标点迹在方位上的分裂现象,同时可以实现改善目标点迹在距离上的分裂现象。In view of the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide a target point trace condensation method for radar echo data, which not only adopts the binary sliding window detection in the sliding window method, but also adopts the morphological filtering algorithm in mathematics. The connected domain is expanded and corroded. By selecting appropriate structural elements, the splitting phenomenon of the target point trace in the sliding window method can be improved not only in azimuth, but also in distance.
为达到上述目的,本发明采用如下技术方案予以实现。In order to achieve the above object, the present invention adopts the following technical solutions to achieve.
一种雷达回波数据的目标点迹凝聚方法,所述方法包括如下步骤:A method for condensing target point traces of radar echo data, the method comprises the following steps:
步骤1,获取雷达经过脉冲压缩处理的多个扫描周期的原始回波数据,所述原始回波数据为包含距离维和方位维的二维数据,且所述多个扫描周期的原始回波数据沿方位维依次排列;Step 1: Obtain the raw echo data of multiple scanning cycles of the radar after pulse compression processing, the raw echo data is two-dimensional data including the distance dimension and the azimuth dimension, and the raw echo data of the multiple scanning cycles are The orientation dimension is arranged in order;
步骤2,对所述原始回波数据依次进行包络检波和非相参积累,得到脉冲积累后的回波数据;Step 2, sequentially performing envelope detection and non-coherent accumulation on the original echo data to obtain echo data after pulse accumulation;
步骤3,对所述脉冲积累后的回波数据进行恒虚警检测,得到恒虚警检测后的回波数据,进而对所述恒虚警检测后的回波数据依次进行二值量化处理以及沿方位维的二值滑窗检测,得到二值滑窗检测后的回波数据。Step 3, performing constant false alarm detection on the echo data after the pulse accumulation, to obtain echo data after constant false alarm detection, and then performing binary quantization processing on the echo data after constant false alarm detection in turn; The binary sliding window detection along the azimuth dimension is used to obtain the echo data after the binary sliding window detection.
步骤4,将所述二值滑窗检测后的回波数据转化为对应的二值图像数据,且二值图像的长度对应所述二值滑窗检测后的回波数据的距离单元,二值图像的宽度对应所述二值滑窗检测后的回波数据的脉冲数,对所述二值图像数据进行先膨胀后腐蚀的操作,从而得到所述二值图像数据中的所有连通域;Step 4: Convert the echo data detected by the binary sliding window into corresponding binary image data, and the length of the binary image corresponds to the distance unit of the echo data detected by the binary sliding window. The width of the image corresponds to the pulse number of the echo data detected by the binary sliding window, and the binary image data is first expanded and then eroded to obtain all connected domains in the binary image data;
步骤5,在所述所有连通域中,滤除第一连通域和第二连通域,所述第一连通域为只包含一个孤立点迹的连通域,所述第二连通域为在距离维上扩展超过门限值的点迹组成的连通域,从而得到包含目标点迹的剩余连通域;根据所述包含目标点迹的剩余连通域得到目标信息,实现对目标点迹的凝聚。Step 5: In all the connected domains, filter out the first connected domain and the second connected domain, the first connected domain is a connected domain that contains only one isolated point trace, and the second connected domain is a connected domain in the distance dimension. The connected domain composed of the point traces exceeding the threshold value is upwardly expanded, thereby obtaining the remaining connected domain containing the target point trace; the target information is obtained according to the remaining connected domain containing the target point trace, so as to realize the agglomeration of the target point trace.
本发明技术方案的特点和进一步的改进为:The characteristics and further improvement of the technical solution of the present invention are:
(1)在所述步骤1之后,且所述步骤2之前,所述方法还包括:对所述原始回波数据进行二次对消,抑制原始回波数据中的杂波。(1) After the step 1 and before the step 2, the method further includes: performing secondary cancellation on the original echo data to suppress clutter in the original echo data.
(2)所述步骤2具体为:(2) Described step 2 is specifically:
对所述原始回波数据进行包络检波,从而得到原始回波数据在对应距离-方位单元的幅度信息;所述距离单元为所述原始回波数据在距离维上的采样点,所述方位单元为所述原始回波数据在方位维上的采样点。Envelope detection is performed on the original echo data, so as to obtain the amplitude information of the original echo data in the corresponding distance-azimuth unit; the distance unit is the sampling point of the original echo data in the distance dimension, and the azimuth The unit is the sampling point of the original echo data in the azimuth dimension.
对包络检波后的回波数据在相邻脉冲上进行滑窗积累,得到脉冲积累后的回波数据。The echo data after envelope detection is accumulated on adjacent pulses by sliding window to obtain echo data after pulse accumulation.
(3)所述步骤3具体为:(3) the step 3 is specifically:
对所述脉冲积累后的回波数据沿距离单元进行滑窗检测,将窗口内所有参考单元的均值乘以门限因子得到的值设置为滑窗检测的门限值;将待检测单元的幅值与所述滑窗检测的门限值进行比较,若待检测单元的幅值大于所述滑窗检测的门限值,则将待检测单元的数据保留,否则将所述待检测单元的数据设置为0;Sliding window detection is performed on the echo data after the pulse accumulation along the distance unit, and the value obtained by multiplying the mean value of all reference units in the window by the threshold factor is set as the threshold value of sliding window detection; the amplitude value of the unit to be detected is set. Compared with the threshold value of the sliding window detection, if the amplitude of the unit to be detected is greater than the threshold value of the sliding window detection, the data of the unit to be detected is retained, otherwise the data of the unit to be detected is set is 0;
进而将幅值大于所述滑窗检测的门限值的待检测单元的数据设置为1,完成对恒虚警检测后的回波数据进行二值量化处理的过程;Then, the data of the unit to be detected whose amplitude value is greater than the threshold value of the sliding window detection is set to 1, and the process of performing binary quantization processing on the echo data after the constant false alarm detection is completed;
对进行二值量化处理后的回波数据沿方位维进行二值滑窗检测时,采用M/N准则,其中,M表示二值滑窗检测的检测门限,N表示二值滑窗检测时所用窗口的长度。When binary sliding window detection is performed on the echo data after binary quantization processing along the azimuth dimension, the M/N criterion is adopted, where M represents the detection threshold of binary sliding window detection, and N represents the detection threshold of binary sliding window detection. The length of the window.
(4)步骤4中,对所述二值图像数据进行先膨胀后腐蚀的操作,从而得到所述二值图像数据中的所有连通域,具体包括:(4) In step 4, the binary image data is first expanded and then eroded, so as to obtain all connected domains in the binary image data, specifically including:
选取结构元素,用选取的结构元素对二值图像数据进行先膨胀后腐蚀的操作,然后对先膨胀后腐蚀的二值图像数据进行连通域的查找。Select the structuring elements, use the selected structuring elements to dilate and then corrode the binary image data, and then search the connected domain of the binary image data that is dilated and then corroded.
(5)步骤5中,所述根据所述包含目标点迹的剩余连通域得到目标信息,实现对目标点迹的凝聚,具体包括:(5) In step 5, the target information is obtained according to the remaining connected domain including the target point trace, so as to realize the agglomeration of the target point trace, which specifically includes:
获取所述包含目标点迹的剩余连通域的质心位置,将所述质心位置在原始回波数据中对应的距离单元作为目标的距离信息,将所述质心位置在原始回波数据中对应的方位单元作为目标的方位信息;Obtain the centroid position of the remaining connected domain including the target trace, take the distance unit corresponding to the centroid position in the original echo data as the distance information of the target, and take the centroid position in the original echo data The corresponding azimuth The orientation information of the unit as the target;
获取所述包含目标点迹的剩余连通域对应的原始回波数据,根据所述包含目标点迹的剩余连通域对应的原始回波数据所在的距离-方位单元以及对应的包络检波结果,求得目标的距离信息、方位信息以及幅度信息。Obtain the original echo data corresponding to the remaining connected domain containing the target point trace, according to the distance-azimuth unit where the original echo data corresponding to the remaining connected domain containing the target point trace is located and the corresponding envelope detection result, find Get the distance information, bearing information and amplitude information of the target.
本发明与现有技术相比,具有以下优点:Compared with the prior art, the present invention has the following advantages:
(1)由于本发明将回波数据转换为二值图像进行处理,可以采用软件实现和离线处理;具有基于图像处理和基于滑窗法点迹凝聚方法各自的优点,并且在经典点迹凝聚方法的基础上采用形态滤波方法进一步对点迹进行处理,弥补了两种经典方法的不足,进一步改善了目标点迹的分裂现象,同时具有较高的点迹质量;(1) Since the present invention converts echo data into binary images for processing, software implementation and offline processing can be used; it has respective advantages based on image processing and point trace condensation methods based on the sliding window method, and is in the classical point trace condensation method. On the basis of the morphological filtering method, the point trace is further processed, which makes up for the shortcomings of the two classical methods, further improves the splitting phenomenon of the target point trace, and at the same time has a higher point trace quality;
(2)本发明提供的目标点迹凝聚方法采用二值图像连通域查找方法,能精确地确定目标点迹分布区域,对目标点迹有较高的分辨率,查找后采用质心法确定区域质心或者采用回波幅度加权方法对目标参数进行估计,目标参数估计比滑窗法更准确;(2) The target point trace agglomeration method provided by the present invention adopts a binary image connected domain search method, which can accurately determine the distribution area of the target point trace, and has a higher resolution for the target point trace. After searching, the centroid method is used to determine the regional centroid Or use the echo amplitude weighting method to estimate the target parameters, and the target parameter estimation is more accurate than the sliding window method;
(3)对于目标点迹的判别和虚假点迹的滤除更为有效,首先结合基于滑窗法目标点迹的判别和点迹凝聚准则以及目标分布特点,可以制定这样的目标判别准则:对距离上延伸较大的目标点迹进行滤除和孤立的点迹进行滤除,其次对连通域面积较小的进行滤除(具体参数要根据情况进行设定,没有统一的标准),这些目标判断和点迹滤除准则能有效识别目标点迹,同时能灵活控制点迹数量;(3) It is more effective for the discrimination of target traces and the filtering of false traces. First, combining the discrimination of target traces based on the sliding window method, the point trace aggregation criterion and the characteristics of target distribution, the following target discrimination criteria can be formulated: The target point traces with a larger extension in the distance are filtered out and the isolated point traces are filtered out, and then the small connected domain area is filtered out (the specific parameters should be set according to the situation, there is no unified standard). Judgment and dot trace filtering criteria can effectively identify target dot traces, and at the same time can flexibly control the number of dot traces;
(4)通过图像膨胀腐蚀操作进行连通域合并能有效地减少目标分裂造成的影响,进一步改善基于滑窗法点迹凝聚的效果。(4) Connected domain merging through image expansion and erosion operation can effectively reduce the impact of target splitting, and further improve the effect of point trace condensation based on the sliding window method.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,,还可以根据这些附图获得其他的附图。在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention, and for those of ordinary skill in the art, other drawings can also be obtained according to these drawings. On the premise of no creative work, other drawings can also be obtained from these drawings.
图1为本发明实施例提供的一种雷达回波数据的目标点迹凝聚方法的流程示意图;1 is a schematic flowchart of a method for condensing target point traces of radar echo data according to an embodiment of the present invention;
图2为本发明实施例提供的二值滑窗检测过程示意图;2 is a schematic diagram of a binary sliding window detection process provided by an embodiment of the present invention;
图3为本发明实施例提供的不同结构元素对同一目标的膨胀结果示意图;3 is a schematic diagram of the expansion results of different structural elements on the same target provided by an embodiment of the present invention;
图4为本发明实施例提供的不同结构元素对同一目标的腐蚀结果示意图;4 is a schematic diagram of the corrosion results of different structural elements on the same target according to an embodiment of the present invention;
图5为本发明实施例提供的像素邻接关系示意图;5 is a schematic diagram of a pixel adjacency relationship provided by an embodiment of the present invention;
图6为本发明实施例提供的P、Q、S点的连通关系示意图;6 is a schematic diagram of the connectivity relationship of points P, Q, and S provided by an embodiment of the present invention;
图7为本发明实施例提供的连通域查找过程示意图;7 is a schematic diagram of a connected domain search process provided by an embodiment of the present invention;
图8为本发明实施例提供的原始回波数据示意图;FIG. 8 is a schematic diagram of original echo data provided by an embodiment of the present invention;
图9为本发明实施例提供的非相参积累后的回波数据示意图;9 is a schematic diagram of echo data after non-coherent accumulation provided by an embodiment of the present invention;
图10为本发明实施例提供的恒虚警检测结果后的回波数据示意图;10 is a schematic diagram of echo data after a constant false alarm detection result provided by an embodiment of the present invention;
图11为本发明实施例提供的二值滑窗检测结果示意图;11 is a schematic diagram of a binary sliding window detection result provided by an embodiment of the present invention;
图12为本发明实施例提供的二值图像膨胀腐蚀结果示意图;12 is a schematic diagram of a binary image dilation corrosion result provided by an embodiment of the present invention;
图13为本发明实施例提供的目标点迹凝聚结果示意图。FIG. 13 is a schematic diagram of a result of target point trace condensation according to an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
雷达技术的发展要求雷达系统提供更高的目标参数,天线波束对目标进行扫描后的回波数据经信号处理和目标检测后会产生分裂现象,分为距离分裂和方位分裂。对于两坐标雷达来说点迹凝聚就是对目标一次信息(距离方位)的提取。通常对目标信息的提取采用二元滑窗检测器,滑窗检测器减少虚警的程度取决于检测门限和窗口长度。点迹凝聚一方面要求对虚假点迹进行滤除,另一方面要求有较高的分辨能力。能否减少目标分裂带来的影响同时提高目标参数质量对于点迹凝聚来说是一个重要的指标。本发明实施例提供的基于图像处理的点迹凝聚方法一方面算法在二值图像中实现,不受其他条件影响,较为容易实现;另一方面图像处理与基于滑窗法点迹凝聚相结合,增加了对点迹的判别准则,有利于提高目标点迹质量。本发明实施例提供的基于图像处理的点迹凝聚方法是基于滑窗法点迹凝聚的思想和原理与图像处理相关算法相结合,并在此基础上进行改进的结果。The development of radar technology requires the radar system to provide higher target parameters. After the antenna beam scans the target, the echo data will be split after signal processing and target detection, which can be divided into distance splitting and azimuth splitting. For two-coordinate radar, point trace condensation is the extraction of primary information (range and azimuth) of the target. Usually, a binary sliding window detector is used to extract the target information, and the degree of false alarm reduction of the sliding window detector depends on the detection threshold and the window length. On the one hand, dot trace condensation requires filtering out false dot traces, and on the other hand, it requires high resolution. Whether the impact of target splitting can be reduced and the quality of target parameters can be improved is an important indicator for point trace condensation. The image processing-based point trace condensation method provided by the embodiments of the present invention is, on the one hand, implemented in a binary image, and is not affected by other conditions, so it is relatively easy to implement; The criterion for the discrimination of the point trace is added, which is beneficial to improve the quality of the target point trace. The image processing-based point trace condensation method provided by the embodiment of the present invention is the result of combining the idea and principle of point trace condensation based on the sliding window method with an image processing related algorithm, and improving on this basis.
本发明实施例提供一种雷达回波数据的目标点迹凝聚方法,如图1所示,所述方法包括如下步骤:An embodiment of the present invention provides a method for condensing target point traces of radar echo data. As shown in FIG. 1 , the method includes the following steps:
步骤1,获取雷达经过脉冲压缩处理的多个扫描周期的原始回波数据,所述原始回波数据为包含距离维和方位维的二维数据,且所述多个扫描周期的原始回波数据沿方位维依次排列。Step 1: Obtain the raw echo data of multiple scanning cycles of the radar after pulse compression processing, the raw echo data is two-dimensional data including the distance dimension and the azimuth dimension, and the raw echo data of the multiple scanning cycles are The orientation dimensions are arranged in order.
需要补充的是,在所述步骤1之后,且所述步骤2之前,所述方法还包括:对所述原始回波数据进行二次对消,抑制原始回波数据中的杂波。It should be supplemented that, after the step 1 and before the step 2, the method further includes: performing secondary cancellation on the original echo data to suppress clutter in the original echo data.
步骤2,对所述原始回波数据依次进行包络检波和多脉冲非相参积累,得到脉冲积累后的回波数据。Step 2: Envelope detection and multi-pulse non-coherent accumulation are sequentially performed on the original echo data to obtain echo data after pulse accumulation.
脉冲积累分为相参积累和非相参积累。相参积累能较好地改善目标信噪比,但由于运算复杂度较高,而且相参积累后的点迹凝聚估计误差较大,尤其是当目标多普勒频率出现在多个通道中时或者目标沿切向运动时凝聚的效果不好,对于两坐标雷达来说,目标主要参数为距离、方位和速度,根据信号目标检测后的数据对目标参数进行提取。采用非相参积累运算复杂度低,在很多场景下适应性很强,因此本发明技术方案采用多脉冲非相参积累。Pulse accumulation is divided into coherent accumulation and non-coherent accumulation. Coherent accumulation can improve the target signal-to-noise ratio well, but due to the high computational complexity and the large error in the estimation of point trace aggregation after coherent accumulation, especially when the target Doppler frequency appears in multiple channels Or when the target moves tangentially, the condensing effect is not good. For the two-coordinate radar, the main parameters of the target are distance, azimuth and speed, and the target parameters are extracted according to the data after the signal target is detected. Adopting non-coherent accumulation has low computational complexity and strong adaptability in many scenarios, so the technical solution of the present invention adopts multi-pulse non-coherent accumulation.
所述步骤2具体为:The step 2 is specifically:
对所述原始回波数据进行包络检波,从而得到原始回波数据在对应距离-方位单元的幅度信息;所述距离单元为所述原始回波数据在距离维上的采样点,所述方位单元为所述原始回波数据在方位维上的采样点。Envelope detection is performed on the original echo data, so as to obtain the amplitude information of the original echo data in the corresponding distance-azimuth unit; the distance unit is the sampling point of the original echo data in the distance dimension, and the azimuth The unit is the sampling point of the original echo data in the azimuth dimension.
对包络检波后的回波数据在相邻脉冲上进行滑窗积累,得到脉冲积累后的回波数据。The echo data after envelope detection is accumulated on adjacent pulses by sliding window to obtain echo data after pulse accumulation.
步骤3,对所述脉冲积累后的回波数据进行恒虚警检测,得到恒虚警检测后的回波数据,进而对所述恒虚警检测后的回波数据依次进行二值量化处理以及沿方位维的二值滑窗检测,得到二值滑窗检测后的回波数据,进一步降低虚假点迹带来的影响。所述二值滑窗检测后的回波数据为二维数据,其中一维表示距离单元,另一维表示方位维的脉冲数;Step 3, performing constant false alarm detection on the echo data after the pulse accumulation, to obtain echo data after constant false alarm detection, and then performing binary quantization processing on the echo data after constant false alarm detection in turn; The binary sliding window detection along the azimuth dimension is used to obtain echo data after binary sliding window detection, which further reduces the influence of false point traces. The echo data detected by the binary sliding window is two-dimensional data, wherein one dimension represents the distance unit, and the other dimension represents the number of pulses in the azimuth dimension;
所述恒虚警检测采用单元平均恒虚警检测(CA-CFAR)。The constant false alarm detection adopts cell averaged constant false alarm detection (CA-CFAR).
所述步骤3中,进而对所述恒虚警检测后的回波数据依次进行二值量化处理以及沿方位维的二值滑窗检测,得到二值滑窗检测后的回波数据,具体为:In the step 3, further perform binary quantization processing and binary sliding window detection along the azimuth dimension on the echo data after the constant false alarm detection, to obtain the echo data after the binary sliding window detection, specifically: :
首先对经过步骤2处理得到的脉冲积累后的回波数据进行采用单元平均恒虚警检测处理,单元平均恒虚警检测是沿距离单元进行滑窗检测,将窗口内所有参考单元的均值乘以门限因子得到的值设置为滑窗检测的门限值,将待检测单元的幅值与所述滑窗检测的门限值进行比较,对过门限数据保留,未过门限数据设置为0;而后将过门限值的数据设置为1,从而完成对所述脉冲积累后的回波数据进行恒虚警检测以及二值量化处理的过程;First, the echo data after the pulse accumulation obtained by the processing in step 2 is processed by using the unit average constant false alarm detection. The value obtained by the threshold factor is set as the threshold value of the sliding window detection, the amplitude of the unit to be detected is compared with the threshold value of the sliding window detection, the data that has passed the threshold is retained, and the data that has not passed the threshold is set to 0; then The data that has passed the threshold value is set to 1, so as to complete the process of performing constant false alarm detection and binary quantization processing on the echo data after the pulse accumulation;
对进行二值量化处理后的回波数据沿方位维进行二值滑窗检测时,采用M/N准则,其中,M表示二值滑窗检测的检测门限,N表示二值滑窗检测时所用窗口的长度。When binary sliding window detection is performed on the echo data after binary quantization processing along the azimuth dimension, the M/N criterion is adopted, where M represents the detection threshold of binary sliding window detection, and N represents the detection threshold of binary sliding window detection. The length of the window.
二值滑窗检测的目标起始门限对目标点迹的分辨影响不大,所以这里设置目标终了门限与起始门限相同,M与N的关系按照经验公式M=1.5*sqrt(N)进行选取,但是为了减少目标点迹分裂,二值滑窗检测的检测门限应适当的降低。The target start threshold of binary sliding window detection has little effect on the resolution of target traces, so the target end threshold is set to be the same as the start threshold, and the relationship between M and N is selected according to the empirical formula M=1.5*sqrt(N) , but in order to reduce the target point trace splitting, the detection threshold of the binary sliding window detection should be appropriately reduced.
需要说明的是,目标终了门限值是指滑窗检测过程中对目标的点迹进行分辨时,二值滑窗积累值是一个变化过程,通过选取起始门限和终了门限判断目标点迹在方位上的起始点和终了点,从而完成目标点迹在方位上的分辨。It should be noted that the target end threshold value refers to the process of distinguishing the target traces in the sliding window detection process, and the accumulated value of the binary sliding window is a process of change. The starting point and ending point on the azimuth, so as to complete the resolution of the target trace in the azimuth.
具体的,二值滑窗检测过程如图2所示,其中,图2(a)为二值量化处理后的输入信号,图2(b)为二值滑窗检测积累过程,图2(c)为二值滑窗检测结果示意图;M/N准则是指窗口长度N=5,检测门限M=3,当窗口内二值累加结果大于等于门限值时就认为检测到目标记为1,否则记为0。Specifically, the binary sliding window detection process is shown in Figure 2, in which Figure 2(a) is the input signal after binary quantization processing, Figure 2(b) is the binary sliding window detection accumulation process, Figure 2(c) ) is a schematic diagram of the binary sliding window detection result; the M/N criterion refers to the window length N=5, the detection threshold M=3, when the binary accumulation result in the window is greater than or equal to the threshold value, it is considered that the detected target is marked as 1, Otherwise, it is recorded as 0.
需要补充的是,所述原始回波数据指的是天线扫描后经过脉冲压缩处理后的数据,对回波数据的处理包括杂波抑制、信号积累和恒虚警检测。It should be added that the original echo data refers to the data after the antenna is scanned and processed by pulse compression, and the processing of the echo data includes clutter suppression, signal accumulation and constant false alarm detection.
对于两坐标雷达来说,目标参数指的是目标的距离和方位信息,这里对回波数据进行非相参积累,杂波抑制采用对消器实现。For the two-coordinate radar, the target parameter refers to the distance and azimuth information of the target. Here, the echo data is non-coherently accumulated, and the clutter suppression is realized by a canceller.
但是根据具体情况,比如为了不影响对切向飞行的目标检测,可以对回波数据只进行简单的信号积累而省略二次对消过程,本发明对回波数据处理采用滑窗积累、恒虚警检测和二值滑窗检测的方法。However, according to specific circumstances, for example, in order not to affect the target detection of tangential flight, the echo data can be simply accumulated and the secondary cancellation process can be omitted. The present invention uses sliding window accumulation, constant virtual methods for alarm detection and binary sliding window detection.
步骤4,将所述二值滑窗检测后的回波数据作为二值图像数据,对所述二值图像数据进行先膨胀后腐蚀的操作,从而得到所述二值图像数据中的所有连通域。Step 4: Take the echo data detected by the binary sliding window as binary image data, and perform an operation of dilating and then corroding the binary image data, so as to obtain all connected domains in the binary image data. .
具体的,对所述二值图像数据采用形态滤波算法进行膨胀腐蚀操作,目的是实现连通域的合并。Specifically, a morphological filtering algorithm is used to perform a dilation and erosion operation on the binary image data, in order to realize the merging of connected domains.
图像的膨胀与腐蚀操作,从图像处理角度来看,二值图像处理就是将一个小的结构元素在大的二值图像上逐点进行移动并比较,根据比较的结果做出相应的处理。结构元素是指具有一定尺寸的背景图像,结构元素没有固定的形状,在设计形态变换算法的同时需要根据输入图像和所需的信息进行结构元素的设计。通常结构元素形状有正方形、矩形、圆形和线性等。简单地说,膨胀就是将物体周围的背景点合并到物体中,物体通过膨胀向外扩展,这样对于两个相邻的物体就有可能被连在一起,在二值图像中就表现为连通域的合并。相反的,图像的腐蚀操作就是消除物体的边界点,另外消除较小的物体。Image expansion and erosion operations, from the perspective of image processing, binary image processing is to move and compare a small structural element on a large binary image point by point, and make corresponding processing according to the comparison results. Structural elements refer to background images with a certain size. Structural elements have no fixed shape. When designing morphological transformation algorithms, structural elements need to be designed according to the input image and the required information. Usually the shape of the structuring element is square, rectangle, circle and linear etc. Simply put, expansion is to merge the background points around the object into the object, and the object expands outward through expansion, so that it is possible for two adjacent objects to be connected together, which is represented as a connected domain in a binary image. 's merger. In contrast, the erosion operation of the image is to remove the boundary points of objects, and also to remove smaller objects.
(1)图像的膨胀处理可以简单地描述为结构元素在二值图像上进行遍历,以二值图像黑色点为例,当结构元素原点与二值图像中某点像素相同,那么结构元素与二值图像对应的所有点都变为黑色点。如图3所示为不同结构元素对同一目标的膨胀结果示意图,图3(a)为结构元素1,图3(b)为结构元素2,图3(c)为原二值图像,图3(d)为采用结构元素1对原二值图像进行膨胀的结果,图3(e)为采用结构元素2对原二值图像进行膨胀的结果,其中,结构元素中符号“+”表示坐标原点所在位置。(1) The expansion processing of the image can be simply described as the traversal of the structuring element on the binary image. Taking the black point of the binary image as an example, when the origin of the structuring element is the same as the pixel of a certain point in the binary image, then the structuring element is the same as the binary image. All points corresponding to the value image become black points. Figure 3 is a schematic diagram of the expansion results of different structural elements for the same target, Figure 3 (a) is the structural element 1, Figure 3 (b) is the structural element 2, Figure 3 (c) is the original binary image, Figure 3 (d) is the result of expanding the original binary image with structural element 1, and Figure 3(e) is the result of expanding the original binary image with structural element 2, where the symbol "+" in the structural element represents the origin of coordinates location.
(2)图像的腐蚀操作可以简单描述为结构元素在二值图像上进行遍历时,以二值图像黑色点为例,当结构元素原点与二值图像中某点像素相同,如果结构元素中其余像素有一个及以上的像素与二值图像对应像素点不同,那么原点所对应的二值图像的点变为白色点。如图4所示为不同结构元素对同一目标的腐蚀结果示意图,图4(a)为结构元素1,图4(b)为结构元素2,图4(c)为原二值图像,图4(d)为采用结构元素1对原二值图像进行腐蚀的结果,图4(e)为采用结构元素2对原二值图像进行腐蚀的结果。(2) The corrosion operation of the image can be simply described as the traversal of the structuring element on the binary image, taking the black point of the binary image as an example, when the origin of the structuring element is the same as the pixel of a certain point in the binary image, if the rest of the structuring elements If one or more pixels of the pixel are different from the corresponding pixels of the binary image, the point of the binary image corresponding to the origin becomes a white point. Figure 4 shows the schematic diagram of the corrosion results of different structural elements on the same target, Figure 4(a) is the structural element 1, Figure 4(b) is the structural element 2, Figure 4(c) is the original binary image, Figure 4 (d) is the result of corroding the original binary image by using structural element 1, and FIG. 4(e) is the result of corroding the original binary image by using structural element 2.
本发明实施例采用的结构元素为线段形结构元素,中心点在线段中心位置,长度采用相参处理脉冲数的1/4,对二值图像先膨胀后腐蚀,膨胀腐蚀各进行一次。二值图像膨胀实现连通域的合并,腐蚀操作对膨胀后的二值图像进行恢复,但恢复后的图像与原来的二值图像是不一样的。The structural elements used in the embodiments of the present invention are line-segment structural elements, the center point is at the center of the line segment, and the length is 1/4 of the number of coherent processing pulses. The binary image dilation realizes the merging of connected domains, and the erosion operation restores the dilated binary image, but the restored image is different from the original binary image.
具体的,步骤4中,得到所述二值图像数据中的所有连通域时需要对膨胀腐蚀后的连通域进行查找。Specifically, in step 4, when all the connected domains in the binary image data are obtained, it is necessary to search the connected domains after expansion and corrosion.
邻接与连通是像素间的基本关系,二值图像中除边缘像素以外,每个像素周围都有8个邻接点,根据邻接点的位置进行区分又有4邻接点。如图5所示,图5(a)黑色点周围为8邻接示意图,图5(b)黑色点周围为4邻接示意图。Adjacency and connectivity are the basic relationships between pixels. Except for edge pixels in a binary image, there are 8 adjacent points around each pixel, and 4 adjacent points are distinguished according to the position of the adjacent points. As shown in Fig. 5, Fig. 5(a) is a schematic diagram of 8 adjacencies around the black dots, and Fig. 5(b) is a schematic diagram of 4 adjacencies around the black dots.
在二值图像中如果两个点通过一系列灰度相同的点序列根据4邻接关系或者8邻接关系将连在一起那么就称这两个点是连通域的,这样所有与这两个点连通的点的集合就构成了连通域,如图6所示,像素值为1的点P、Q、S,其中P与Q是8连通的,S与Q是4连通的,S与P是8连通的。In a binary image, if two points are connected together by a series of point sequences with the same gray level according to the 4-adjacency relationship or the 8-adjacency relationship, then the two points are said to be connected domains, so that all the points are connected to these two points. The set of points constitutes a connected domain. As shown in Figure 6, the points P, Q, and S whose pixel value is 1, where P and Q are 8 connected, S and Q are 4 connected, and S and P are 8 connected.
以8连通为例,对像素值为1的像素进行连通域查找,假设A为二值图像中的一个连通域,A中一个点已知为P,那么对于连通域A的查找可以通过下述迭代式进行:Taking 8 connectivity as an example, the connected domain is searched for the pixel whose pixel value is 1. Assuming that A is a connected domain in the binary image, and a point in A is known as P, then the search for the connected domain A can be done through the following Iteratively:
X0=PX 0 =P
当Xk=Xk-1时,算法收敛,并且A=Xk。其中,B表示结构元素,表示图像膨胀操作,Y为原始二值图像。The algorithm converges when Xk = Xk-1 , and A= Xk . where B represents a structural element, Indicates the image dilation operation, and Y is the original binary image.
如图7所示为连通域查找示意图,采用8连通关系,结构元素如图7(a),图7(b)为连通域查找的起始点P,黑色元素为已提取的元素,灰色为尚未提取的元素,图7(c)为查找连通域第一次迭代的结果,图7(d)为查找连通域第二次迭代的结果,图7(e)为连通域查找最终结果示意图。根据连通域查找的方法,查找后的连通域需要进行标记,直到所有的像素为1的点全部都标记完为止,这样就完成了连通域的查找了。Figure 7 is a schematic diagram of the connected domain search, using 8 connectivity relationships, the structural elements are shown in Figure 7(a), Figure 7(b) is the starting point P of the connected domain search, the black elements are the extracted elements, and the gray is not yet. The extracted elements, Figure 7(c) is the result of the first iteration of finding the connected domain, Figure 7(d) is the result of the second iteration of finding the connected domain, and Figure 7(e) is a schematic diagram of the final result of finding the connected domain. According to the connected domain search method, the searched connected domain needs to be marked until all the points whose pixels are 1 are marked, thus completing the connected domain search.
步骤5,在所述所有连通域中,滤除第一连通域和第二连通域,所述第一连通域为只包含一个孤立点迹的连通域,所述第二连通域为在距离维上扩展超过门限值的点迹组成的连通域,从而得到包含目标点迹的剩余连通域;根据所述包含目标点迹的剩余连通域得到目标信息,实现对目标点迹的凝聚。Step 5: In all the connected domains, filter out the first connected domain and the second connected domain, the first connected domain is a connected domain that contains only one isolated point trace, and the second connected domain is a connected domain in the distance dimension. The connected domain composed of the point traces exceeding the threshold value is upwardly expanded, thereby obtaining the remaining connected domain containing the target point trace; the target information is obtained according to the remaining connected domain containing the target point trace, so as to realize the agglomeration of the target point trace.
步骤5中,所述根据所述包含目标点迹的剩余连通域得到目标信息,实现对目标点迹的凝聚,具体包括:In step 5, the target information is obtained according to the remaining connected domain including the target point trace, and the condensation of the target point trace is realized, which specifically includes:
获取所述包含目标点迹的剩余连通域的质心位置,将所述质心位置在原始回波数据中对应的距离单元作为目标的距离信息,将所述质心位置在原始回波数据中对应的方位单元作为目标的方位信息;Obtain the centroid position of the remaining connected domain including the target trace, take the distance unit corresponding to the centroid position in the original echo data as the distance information of the target, and take the centroid position in the original echo data The corresponding azimuth The orientation information of the unit as the target;
获取所述包含目标点迹的剩余连通域对应的原始回波数据,根据所述包含目标点迹的剩余连通域对应的原始回波数据所在的距离单元、方位单元以及对应的包络检波结果,求得目标的距离信息、方位信息以及幅度信息。Obtain the original echo data corresponding to the remaining connected domain including the target point trace, according to the distance unit, the azimuth unit and the corresponding envelope detection result where the original echo data corresponding to the remaining connected domain including the target point trace is located, Obtain the distance information, bearing information and amplitude information of the target.
具体的,本发明实施例中,当实测数据目标回波较弱,为了尽可能保留目标点迹,从而采用的点迹滤除准则为:Specifically, in the embodiment of the present invention, when the target echo of the measured data is weak, in order to retain the target spot trace as much as possible, the adopted spot trace filtering criterion is:
(1)对独立的点迹进行滤除,根据目标回波分布特点,通常目标在方位上或者距离上会发生扩展,因此可以根据目标点迹在距离或方位扩展的大小扩展滤除一些虚假点迹;(1) Filter out independent point traces. According to the distribution characteristics of target echoes, the target usually expands in azimuth or distance. Therefore, some false points can be filtered out according to the size of the target point trace expansion in distance or azimuth. trace;
(2)对回波数据对应的距离上扩展较多的点迹进行滤除,距离单元的最大值根据目标和系统的估计参数来确定;(2) Filter out the point traces with more expansion on the distance corresponding to the echo data, and the maximum value of the distance unit is determined according to the estimated parameters of the target and the system;
(3)在系统允许的情况下,应尽可能保留较多的点迹,所以点迹分辨与滤除准则要根据实际情况增加或者减少。(3) When the system allows, as many dot traces should be retained as possible, so the dot trace discrimination and filtering criteria should be increased or decreased according to the actual situation.
实验内容与结果Experiment content and results
实验1,回波数据为某对空雷达录取的回波数据,如图8所示为某对空数据的截取部分,数据由雷达9个扫描周期内截取的相同方位和距离段上的数据拼接而成。Experiment 1. The echo data is the echo data recorded by a certain air-to-air radar. Figure 8 shows the intercepted part of a certain air-to-air data. The data is spliced from the data on the same azimuth and distance segment intercepted by the radar within 9 scanning periods. made.
目标回波数据显示如图8所示,滑窗积累后的回波数据如图9所示。如图10所示为恒虚警处理后的结果,如图11所示为二值滑窗检测后的结果,图11(a)为二值滑窗检测M/N准则中,M=2,N=5时的检测结果,图11(b)为二值滑窗检测M/N准则中,M=3,N=5时的检测结果。如图12所示为二值图像膨胀腐蚀后的结果,如图13所示为本发明实施例提供的目标点迹凝聚方法的点迹凝聚结果。The target echo data display is shown in Figure 8, and the echo data accumulated by the sliding window is shown in Figure 9. Figure 10 shows the result after constant false alarm processing, Figure 11 shows the result after binary sliding window detection, Figure 11(a) is the binary sliding window detection M/N criterion, M=2, The detection result when N=5, Fig. 11(b) is the detection result when M=3 and N=5 in the binary sliding window detection M/N criterion. FIG. 12 shows the result of the binary image after dilation and corrosion, and FIG. 13 shows the point trace condensation result of the target point trace condensation method provided by the embodiment of the present invention.
由图8可以看出,目标周围存在较强的杂波分布,单次截取的数据很难判断目标点迹的位置,因此通过拼接后的数据可以大致看到目标的位置,可以通过对目标位置的估计对点迹凝聚效果进行验证。As can be seen from Figure 8, there is a strong clutter distribution around the target, and it is difficult to determine the position of the target trace from a single intercepted data. Therefore, the position of the target can be roughly seen through the spliced data. The estimation of point trace condensation effect is verified.
由图9可以看出,目标经过非相参积累后,对于弱目标的改善情况并不明显,但是还是有一定好处的,主要是对沿径向分布的气象杂波有一定平滑作用,对接下来进行恒虚警检测是有利的。It can be seen from Figure 9 that after the target is accumulated non-coherently, the improvement of the weak target is not obvious, but it still has certain benefits, mainly because it has a certain smoothing effect on the meteorological clutter distributed along the radial direction, which is very important for the next step. It is advantageous to perform constant false alarm detection.
由图10可以看出,剩余点迹中目标点迹和虚假点迹受信号处理和恒虚警检测影响比较大,恒虚警检测后与原始回波数据进行对比,通过对目标位置的大致判断可以看到同一目标的点迹分裂现象较为明显,图中标记的地方即为目标分裂的点迹。It can be seen from Figure 10 that the target traces and false traces in the remaining traces are greatly affected by signal processing and constant false alarm detection. After the constant false alarm detection, it is compared with the original echo data. It can be seen that the dot trace splitting phenomenon of the same target is more obvious, and the marked place in the figure is the dot trace of the target splitting.
由图11可以看出,二值滑窗检测能有效改善目标分裂带来的影响,同时降低虚警,但是采用M/N准则时,门限值过大会损失较弱目标的点迹,滑窗检测时会出现目标点迹分裂现象,通常通过降低门限值的方法解决,门限值过低虚假点迹会增多,对点迹质量影响较大。所以门限值需要折中考虑。由图12可以看出,通过对滑窗检测后的二值图像进行膨胀腐蚀操作后一方面保证了目标分辨能力另一方面对目标分裂的改善也是相当有效的。As can be seen from Figure 11, the binary sliding window detection can effectively improve the impact of target splitting and reduce false alarms, but when the M/N criterion is used, if the threshold value is too large, the point traces of weaker targets will be lost. The target point trace splitting phenomenon will occur during detection, which is usually solved by reducing the threshold value. If the threshold value is too low, false traces will increase, which has a great impact on the quality of the point traces. Therefore, the threshold value needs to be considered as a compromise. It can be seen from Figure 12 that, by performing the dilation and corrosion operation on the binary image detected by the sliding window, on the one hand, the target resolution is guaranteed, and on the other hand, the improvement of target splitting is also quite effective.
由图13可以看出,改进后的方法成功实现了点迹凝聚,目标点迹凝聚取得了较好的效果。It can be seen from Figure 13 that the improved method successfully achieves point trace condensation, and the target point trace condensation achieves good results.
综上,由实验结果可知,本发明方法能够有效的实现点迹凝聚,同时减少目标分裂对点迹凝聚的影响。To sum up, it can be seen from the experimental results that the method of the present invention can effectively realize the point trace agglomeration, and at the same time reduce the influence of target splitting on the point trace agglomeration.
本发明实施例提供的技术方案相比滑窗法的点迹凝聚,基于图像处理的点迹凝聚技术大大提高了目标参数精度和目标分辨能力,而且图像中相关算法比如本发明采用的二值图像膨胀腐蚀方法可以进一步改善滑窗法点迹凝聚中可能存在目标分裂的问题。滑窗法的好处是能有效克服目标分裂和虚警带来的影响,但是滑窗法门限值的选取对点迹凝聚影响较大,这一点可以通过图像处理中形态滤波进行改善。所以结合滑窗法和图像处理的优点,将点迹变换到图像域进行处理,采用连通域查找算法和图像膨胀腐蚀处理,一方面采用二值滑窗检测对目标回波数据进行处理和点迹判决,另一方面可以利用图像处理方法克服点迹分裂,从而获得较高质量的目标估计参数。本专利所提的方法结合了滑窗法和图像处理方法,改进方法主要体现在一是采用滑窗检测器和恒虚警检测对回波数据的处理,第二是采用图像处理中形态滤波方法改善目标点迹分裂的影响。改进的点迹凝聚方法对点迹凝聚技术的研究有着重要的参考价值和研究意义,实测数据验证了所提方法的有效性。Compared with the point trace condensation of the sliding window method, the technical solutions provided by the embodiments of the present invention greatly improve the target parameter accuracy and target resolution ability based on the image processing-based point trace condensation technology, and the related algorithms in the image, such as the binary image used in the present invention The dilation corrosion method can further improve the problem of target splitting in the point trace condensation of the sliding window method. The advantage of the sliding window method is that it can effectively overcome the influence of target splitting and false alarms, but the selection of the threshold value of the sliding window method has a great influence on the point trace condensation, which can be improved by morphological filtering in image processing. Therefore, combining the advantages of the sliding window method and image processing, the point trace is transformed into the image domain for processing, and the connected domain search algorithm and image expansion and corrosion processing are used. On the one hand, binary sliding window detection is used to process the target echo data and point traces. On the other hand, image processing methods can be used to overcome dot trace splitting, so as to obtain higher quality target estimation parameters. The method proposed in this patent combines the sliding window method and the image processing method. The improved method is mainly reflected in the first is to use the sliding window detector and the constant false alarm detection to process the echo data, and the second is to use the morphological filtering method in the image processing. Improves the effect of target trail splitting. The improved point trace condensation method has important reference value and research significance for the research of point trace condensation technology, and the measured data verify the effectiveness of the proposed method.
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。The above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited to this. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed by the present invention. should be included within the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111371.0A CN106918807B (en) | 2017-02-28 | 2017-02-28 | A kind of Targets Dots condensing method of radar return data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111371.0A CN106918807B (en) | 2017-02-28 | 2017-02-28 | A kind of Targets Dots condensing method of radar return data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106918807A CN106918807A (en) | 2017-07-04 |
CN106918807B true CN106918807B (en) | 2019-07-09 |
Family
ID=59454618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710111371.0A Active CN106918807B (en) | 2017-02-28 | 2017-02-28 | A kind of Targets Dots condensing method of radar return data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106918807B (en) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107561518B (en) * | 2017-07-27 | 2020-10-09 | 中国船舶重工集团公司第七二四研究所 | Three-coordinate radar point trace condensation method based on two-dimensional sliding window local extreme value |
CN107861107B (en) * | 2017-10-23 | 2021-02-05 | 电子科技大学 | Double-threshold CFAR (computational fluid dynamics) and trace point agglomeration method suitable for continuous wave radar |
CN108254727A (en) * | 2017-12-08 | 2018-07-06 | 西安电子科技大学 | A Contour Tracking-Based Coalescing Method of Radar Spots |
CN108490410B (en) * | 2018-03-06 | 2021-08-17 | 中国船舶重工集团公司第七二四研究所 | Two-coordinate radar sea target joint detection and tracking method |
CN108120961B (en) * | 2018-03-08 | 2023-07-14 | 合肥乾升电子科技有限公司 | Method and device for realizing PPI display enhancement of detection video based on target azimuth expansion |
CN108680904A (en) * | 2018-05-02 | 2018-10-19 | 深圳市华讯方舟雷达技术装备有限公司 | A kind of suppressing method of marine navigation radar co-channel interference |
CN108614264B (en) * | 2018-05-04 | 2023-08-04 | 武汉华讯国蓉科技有限公司 | Radar target point trace condensation method based on connected marking rule |
CN108663667A (en) * | 2018-05-11 | 2018-10-16 | 清华大学 | Range extension target detection method and system under the incomplete observation of radar |
CN109061593B (en) * | 2018-07-31 | 2020-12-01 | 安徽四创电子股份有限公司 | Target extraction method based on ground monitoring radar spectrum data |
CN109100697B (en) * | 2018-07-31 | 2020-12-22 | 安徽四创电子股份有限公司 | Target condensation method based on ground monitoring radar system |
CN109375177B (en) * | 2018-08-30 | 2020-12-08 | 安徽四创电子股份有限公司 | Moving target detection method for airport scene surveillance radar system |
CN109254271B (en) * | 2018-08-30 | 2021-05-04 | 安徽四创电子股份有限公司 | Static target suppression method for ground monitoring radar system |
CN109375185A (en) * | 2018-10-12 | 2019-02-22 | 中国电子科技集团公司第十四研究所 | A kind of point mark filtering technique based on bianry image connected domain area measurement |
CN111175709B (en) * | 2019-12-27 | 2023-02-24 | 中国电子科技集团公司第十四研究所 | Error suppression-based large-range meteorological radar-oriented jigsaw method |
CN111366924B (en) * | 2020-05-09 | 2022-11-25 | 苏州雷森电子科技有限公司 | Radar target point trace aggregation method based on multi-dimensional connected domain processing |
CN112051568B (en) * | 2020-09-01 | 2023-06-30 | 西安电子工程研究所 | Pitching angle measurement method of two-coordinate radar |
CN112541872B (en) * | 2020-12-08 | 2024-02-27 | 陕西长岭电子科技有限责任公司 | Meteorological radar display effect optimization method |
CN112305508B (en) * | 2020-12-23 | 2021-04-30 | 北京海兰信数据科技股份有限公司 | Echo splitting processing method and system based on marine radar |
CN112835001B (en) * | 2020-12-29 | 2023-08-08 | 西安电子科技大学 | A sea surface target radar spot condensing method, electronic equipment, and storage medium |
CN113868496A (en) * | 2021-09-14 | 2021-12-31 | 中国人民解放军海军航空大学航空作战勤务学院 | Data association method and device, electronic equipment and computer readable storage medium |
CN113820685B (en) * | 2021-10-19 | 2024-07-19 | 中安锐达(北京)电子科技有限公司 | Primary trace point condensation method based on fan block structure |
CN114415123B (en) * | 2022-04-01 | 2022-07-19 | 北京海兰信数据科技股份有限公司 | Non-coherent neighborhood based weighting pulse accumulation processing method and system |
CN118884356B (en) * | 2024-09-23 | 2025-01-03 | 安徽隼波科技有限公司 | Point trace condensation method of millimeter wave radar extended target based on morphology |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288945A (en) * | 2011-05-11 | 2011-12-21 | 成都成电电子信息技术工程有限公司 | Image enhancing method for ship radar |
CN104036146A (en) * | 2014-06-26 | 2014-09-10 | 中国电子科技集团公司第二十八研究所 | Trace point clustering method for clustering trace points of radar targets |
CN104391294A (en) * | 2014-11-27 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | Connection domain characteristic and template matching based radar plot correlation method |
CN105044699A (en) * | 2015-07-08 | 2015-11-11 | 中国电子科技集团公司第二十八研究所 | Radar plot centroid method based on Radon-Fourier transformation |
CN105259553A (en) * | 2015-11-11 | 2016-01-20 | 西安电子科技大学 | Micro-motion target scattering point track association method based on distance-instant Doppler image |
RU2596610C1 (en) * | 2015-06-16 | 2016-09-10 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" | Method of search and detection of object |
CN106093946A (en) * | 2016-05-27 | 2016-11-09 | 四川九洲空管科技有限责任公司 | A kind of target condensing method being applicable to scene surveillance radar and device |
CN106093890A (en) * | 2016-08-26 | 2016-11-09 | 零八电子集团有限公司 | The method setting up residual clutter suppression radar residual spur |
-
2017
- 2017-02-28 CN CN201710111371.0A patent/CN106918807B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288945A (en) * | 2011-05-11 | 2011-12-21 | 成都成电电子信息技术工程有限公司 | Image enhancing method for ship radar |
CN104036146A (en) * | 2014-06-26 | 2014-09-10 | 中国电子科技集团公司第二十八研究所 | Trace point clustering method for clustering trace points of radar targets |
CN104391294A (en) * | 2014-11-27 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | Connection domain characteristic and template matching based radar plot correlation method |
RU2596610C1 (en) * | 2015-06-16 | 2016-09-10 | Федеральное государственное унитарное предприятие "Государственный научно-исследовательский институт авиационных систем" | Method of search and detection of object |
CN105044699A (en) * | 2015-07-08 | 2015-11-11 | 中国电子科技集团公司第二十八研究所 | Radar plot centroid method based on Radon-Fourier transformation |
CN105259553A (en) * | 2015-11-11 | 2016-01-20 | 西安电子科技大学 | Micro-motion target scattering point track association method based on distance-instant Doppler image |
CN106093946A (en) * | 2016-05-27 | 2016-11-09 | 四川九洲空管科技有限责任公司 | A kind of target condensing method being applicable to scene surveillance radar and device |
CN106093890A (en) * | 2016-08-26 | 2016-11-09 | 零八电子集团有限公司 | The method setting up residual clutter suppression radar residual spur |
Also Published As
Publication number | Publication date |
---|---|
CN106918807A (en) | 2017-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106918807B (en) | A kind of Targets Dots condensing method of radar return data | |
CN107145874B (en) | Ship target detection and identification method in complex background SAR image | |
CN110031843B (en) | ROI (region of interest) -based SAR (synthetic Aperture Radar) image target positioning method, system and device | |
CN112130142B (en) | Method and system for extracting micro Doppler features of complex moving target | |
CN110852992B (en) | Ship trail detection method based on top hat transformation and Radon transformation | |
CN109597065B (en) | False alarm suppression method and device for through-wall radar detection | |
CN112731307B (en) | RATM-CFAR detector based on distance-angle joint estimation and detection method | |
CN107479037B (en) | PD radar clutter area distinguishing method | |
CN113050071A (en) | Laser radar data processing method, device, equipment and storage medium | |
CN109597045B (en) | Static target steady recognition method based on two times of clutter suppression | |
CN110208767A (en) | A kind of radar target rapid detection method based on fitting correlation coefficient | |
CN108305265B (en) | Real-time processing method and system for weak and small target images | |
CN112630743A (en) | Sea surface target point trace condensation method and system based on contour tracking area growth | |
CN113534096B (en) | LiDAR signal geometric feature extraction method and system based on spline function | |
CN110133680A (en) | A Method and System for Determining Initial Value of Effective Signal of Airborne Laser Sounding Received Waveform | |
CN115407282B (en) | A SAR active deception jamming detection method based on interferometric phase under short baseline | |
CN107367730B (en) | Self-focusing method suitable for imaging scene targets by strip synthetic aperture sonar | |
CN113589283B (en) | A method for extracting ship kelvin wake elevation based on space-borne interferometric imaging altimeter | |
CN113625266A (en) | Method, device, storage medium and equipment for detecting low-speed target by using radar | |
CN111624605B (en) | A Navigational Radar Target Detection Method Based on Echo Feature in Angle Dimension | |
CN108985288A (en) | A kind of SAR image oil spilling detection method based on TGMSERs | |
CN111369507B (en) | Wake Detection Method Based on Normalized Gray Hough Transform and Local CFAR | |
Wu | Two-parameter CFAR ship detection algorithm based on Rayleigh distribution in SAR images | |
CN115220002A (en) | A fixed single-station multi-target data association tracking method and related device | |
CN114820384A (en) | Infrared weak and small target detection method based on local contrast and gradient |
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 |