CN114779225B - 基于迭代栅格的分布式雷达目标检测方法及装置 - Google Patents

基于迭代栅格的分布式雷达目标检测方法及装置 Download PDF

Info

Publication number
CN114779225B
CN114779225B CN202210712106.9A CN202210712106A CN114779225B CN 114779225 B CN114779225 B CN 114779225B CN 202210712106 A CN202210712106 A CN 202210712106A CN 114779225 B CN114779225 B CN 114779225B
Authority
CN
China
Prior art keywords
grid
current
channel
detection
grids
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
CN202210712106.9A
Other languages
English (en)
Other versions
CN114779225A (zh
Inventor
吕文超
宋扬
葛建军
刘光宏
韩阔业
武艳伟
裴晓帅
王欢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC Information Science Research Institute
Original Assignee
CETC Information Science 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 CETC Information Science Research Institute filed Critical CETC Information Science Research Institute
Priority to CN202210712106.9A priority Critical patent/CN114779225B/zh
Publication of CN114779225A publication Critical patent/CN114779225A/zh
Application granted granted Critical
Publication of CN114779225B publication Critical patent/CN114779225B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本公开涉及雷达目标检测领域，提供一种基于迭代栅格的分布式雷达目标检测方法及装置。本公开首先基于多个尺寸长方体的栅格离线构建三维探测空间对应的多个栅格组以及各栅格组与各通道距离单元对应的索引信息集，有效保证迭代栅格检测的时效性。之后分别利用栅格尺寸不同的栅格组对三维探测空间进行迭代检测，即不断基于上一次检测获得的存在目标的位置及其邻域范围，缩小检测范围，利用缩小后的检测范围对应的栅格组进行栅格检测，直至达到预设迭代次数，得到最终的目标检测结果，在保证检测精度的前提下，有效提升检测效率。本公开还可以不断利用尺寸更小的栅格在对应的检测范围内进行迭代检测，进一步提升检测精度。

Description

基于迭代栅格的分布式雷达目标检测方法及装置
技术领域
本公开涉及雷达目标检测技术领域，特别涉及一种基于迭代栅格的分布式雷达目标检测方法及装置。
背景技术
分布式雷达系统通常包含控制中心与多部雷达节点，控制中心与各雷达节点通过局域网形成探测网络，由控制中心控制各雷达节点，实现栅格联合检测，从而有效提升复杂电磁环境中目标探测性能。因此，基于分布式雷达系统的栅格联合检测技术已成为雷达探测领域的重点研究方向之一。
然而，传统的栅格检测方法只对探测区域进行一次栅格划分，并遍历各栅格区域完成目标检测。并且，为了保证较高的检测精度，传统的栅格检测方法中栅格尺寸往往需要与雷达节点的距离分辨率接近。因此，在面对大探测场景时，上述传统的栅格检测方法不仅要遍历大量栅格，而且计算量庞大，难以满足高效能检测的发展需求。
发明内容
本公开旨在至少解决现有技术中存在的问题之一，提供一种基于迭代栅格的分布式雷达目标检测方法及装置。
本公开的一个方面，提供了一种基于迭代栅格的分布式雷达目标检测方法，包括以下步骤：
离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与三维探测空间相对应的多个栅格组，其中，初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸，经度维尺寸、纬度维尺寸和高度维尺寸均远大于雷达距离分辨率；
分别计算各栅格组中各栅格与各雷达节点对应的方位、俯仰和距离；
离线构建各栅格组与各通道距离单元对应的索引信息集，索引信息集包括栅格组中各栅格在对应通道下的索引信息，索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息；
预处理各通道的回波数据，得到各通道的预处理结果序列，预处理包括匹配滤波处理和运动目标探测处理；
确定当前栅格组中存在目标的栅格：基于预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出当前栅格组中存在目标的栅格，其中，当前栅格组为与当前迭代次数相对应的栅格组；
基于M/N准则进行逻辑判决：将当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于预设单通道门限的栅格组成第一门限栅格集合；将第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计第一门限栅格集合中各栅格对应的单通道检测统计量大于预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断第一门限栅格集合中各栅格对应的M/N是否大于预设单通道门限，将M/N大于预设单通道门限的栅格组成第二门限栅格集合；
序贯凝聚并镜像消除：对第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将点迹集合中的点迹按累加结果降序排序，将排序后的点迹集合中的第一个点迹存入最终点迹集合，并将第一个点迹从点迹集合中删除，得到新的点迹集合，其中，累加结果由与点迹相对应的栅格的各接收通道的单通道检测统计量累加得到；
判断镜像消除是否完成：判断新的点迹集合是否为空集：若新的点迹集合不为空集，则将新的点迹集合作为当前栅格组，回到确定当前栅格组中存在目标的栅格的步骤；若新的点迹集合为空集，则判断当前迭代次数是否达到预设迭代次数；
判断当前迭代次数是否达到预设迭代次数：若当前迭代次数未达到预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，回到确定当前栅格组中存在目标的栅格的步骤，其中，指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到预设迭代次数，则将最终点迹集合作为目标检测结果。
可选的，离线构建三维空间栅格，具体包括：
分别利用尺寸为(q1×q2×q3)/(k3)0,……,(q1×q2×q3)/(k3) n ,……,(q1×q2×q3)/(k3)U的长方体的栅格，对三维探测空间进行均匀划分，得到对应的栅格组grid0,……,grid n ,……,gridU；
其中，q1代表初始栅格尺寸中栅格的经度维尺寸，q2代表初始栅格尺寸中栅格的纬度维尺寸，q3代表初始栅格尺寸中栅格的高度维尺寸，n=0,1,2,…,U表示栅格组的编号，U表示预设迭代次数，k表示预设迭代栅格尺寸比例系数。
可选的，索引信息集表示为：
Figure DEST_PATH_IMAGE001
其中，I c0,……,I cn ,……,I cU分别表示栅格组grid0,……,grid n ,……,gridU对应的索引信息集，
Figure 100002_DEST_PATH_IMAGE002
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射方位角度检索信息，
Figure 100002_DEST_PATH_IMAGE003
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收方位角度检索信息，
Figure 100002_DEST_PATH_IMAGE004
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射俯仰角检索信息，
Figure 100002_DEST_PATH_IMAGE005
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收俯仰角检索信息，I R0,……,I Rn ,……,I RU分别表示栅格组grid0,……,grid n ,……,gridU中栅格的距离检索信息。
可选的，指定区域通过以下步骤确定：
将最终点迹集合中各栅格的位置信息存储为位置集合，位置信息包括经度坐标、纬度坐标、高度坐标；
基于位置集合、初始栅格尺寸、预设迭代栅格尺寸比例系数以及上一个当前迭代次数，在上一个当前栅格组中确定指定区域。
可选的，位置集合表示为Location=[(Jq,Vq,Gq) | q=1,2,…,Nloc]，其中，Jq为栅格q的经度坐标，Vq为栅格q的纬度坐标，Gq为栅格q的高度坐标，q=1,2,…,Nloc为最终点迹集合中的栅格编号，Nloc为最终点迹集合中的栅格数量；
指定区域为上一个当前栅格组中经度处于[Jq-(3q1/2k u-1), Jq+(3q1/2k u-1)]范围、纬度处于[Vq-(3q2/2k u-1), Vq+(3q2/2k u-1)]范围、高度处于[Gq-(3q3/2k u-1), Gq+(3q3/2k u-1)]范围的区域，其中，u为当前迭代次数。
可选的，确定当前栅格组中存在目标的栅格，具体包括：
基于预处理结果序列以及与当前迭代次数对应的距离检索信息，分别计算当前栅格组中各栅格的各接收通道的单通道检测统计量；
将同一栅格的各接收通道的单通道检测统计量进行累加，得到当前栅格组中各栅格分别对应的累加结果；
将累加结果与预设检测器门限进行比较，选取累加结果大于预设检测器门限的栅格，作为当前栅格组中存在目标的栅格。
可选的，在将排序后的点迹集合中的第一个点迹存入最终点迹集合之后，所述方法还包括：根据当前栅格组对应的索引信息集查找第一个点迹的位置对应的各接收通道的预处理结果序列，将查找出的预处理结果序列置为0；
在新的点迹集合不为空集时，所述方法还包括：清除当前栅格组对应的索引信息集中不在新的点迹集合中的栅格的索引信息；
在当前迭代次数未达到预设迭代次数时，所述方法还包括：清除当前栅格组对应的索引信息集中不在指定区域的栅格的索引信息。
本公开的另一个方面，提供了一种基于迭代栅格的分布式雷达目标检测装置，包括：
第一构建模块，用于离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与三维探测空间相对应的多个栅格组，其中，初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸，经度维尺寸、纬度维尺寸和高度维尺寸均远大于雷达距离分辨率；
计算模块，用于分别计算各栅格组中各栅格与各雷达节点对应的方位、俯仰和距离；
第二构建模块，用于离线构建各栅格组与各通道距离单元对应的索引信息集，索引信息集包括栅格组中各栅格在对应通道下的索引信息，索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息；
预处理模块，用于预处理各通道的回波数据，得到各通道的预处理结果序列，预处理包括匹配滤波处理和运动目标探测处理；
确定模块，用于确定当前栅格组中存在目标的栅格：基于预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出当前栅格组中存在目标的栅格，其中，当前栅格组为与当前迭代次数相对应的栅格组；
判决模块，用于基于M/N准则进行逻辑判决：将当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于预设单通道门限的栅格组成第一门限栅格集合；将第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计第一门限栅格集合中各栅格对应的单通道检测统计量大于预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断第一门限栅格集合中各栅格对应的M/N是否大于预设单通道门限，将M/N大于预设单通道门限的栅格组成第二门限栅格集合；
消除模块，用于序贯凝聚并镜像消除：对第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将点迹集合中的点迹按累加结果降序排序，将排序后的点迹集合中的第一个点迹存入最终点迹集合，并将第一个点迹从点迹集合中删除，得到新的点迹集合，其中，累加结果由与点迹相对应的栅格的各接收通道的单通道检测统计量累加得到；
第一判断模块，用于判断镜像消除是否完成：判断新的点迹集合是否为空集：若新的点迹集合不为空集，则将新的点迹集合作为当前栅格组，重新触发确定模块；若新的点迹集合为空集，则触发第二判断模块；
第二判断模块，用于判断当前迭代次数是否达到预设迭代次数：若当前迭代次数未达到预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，重新触发确定模块，其中，指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到预设迭代次数，则将最终点迹集合作为目标检测结果。
本公开的另一个方面，提供了一种电子设备，包括：
至少一个处理器；以及，
与至少一个处理器通信连接的存储器；其中，
存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器能够执行前文记载的基于迭代栅格的分布式雷达目标检测方法。
本公开的另一个方面，提供了一种计算机可读存储介质，存储有计算机程序，计算机程序被处理器执行时实现前文记载的基于迭代栅格的分布式雷达目标检测方法。
本公开相对于现有技术而言，离线建立多个尺寸的栅格以及栅格与各通道距离单元的索引关系，有效保证了迭代栅格检测的时效性；利用多个尺寸的栅格对三维探测空间进行迭代检测，在迭代检测过程中不断基于上一次迭代检测获得的存在目标的位置及其邻域范围缩小迭代检测范围直至达到预设迭代次数，有效提升了检测效率；并且，在迭代检测过程中，还可以不断利用尺寸更小的栅格在对应的迭代检测范围内进行检测，从而有效保证了检测精度。
附图说明
一个或多个实施方式通过与之对应的附图中的图片进行示例性说明，这些示例性说明并不构成对实施方式的限定，附图中具有相同参考数字标号的元件表示为类似的元件，除非有特别申明，附图中的图不构成比例限制。
图1为本公开一实施方式提供的一种基于迭代栅格的分布式雷达目标检测方法的流程图；
图2为本公开另一实施方式提供的一种基于迭代栅格的分布式雷达目标检测方法的流程图；
图3为本公开另一实施方式提供的一种基于迭代栅格的分布式雷达目标检测装置的结构示意图；
图4为本公开另一实施方式提供的电子设备的结构示意图。
具体实施方式
现有技术中，基于空间栅格划分的协同探测系统检测方法包括以下步骤：
步骤S1：构建三维空间栅格Grid3d，将三维探测空间均匀划分为经度、纬度、高度分别为△Lo、△La、△H的长方体栅格。
步骤S2：计算各三维栅格与各雷达节点对应的方位角θ，俯仰角
Figure DEST_PATH_IMAGE006
和距离R。
步骤S3：预处理各通道回波数据，即对各通道回波数据进行匹配滤波处理和运动目标探测处理，输出预处理结果序列：
Figure DEST_PATH_IMAGE007
，其中，Y c (l)表示第c个通道中第l个距离单元的预处理结果序列，m为脉冲编号，Mp为脉冲个数。y c [l,m]为第c个通道中第l个距离单元的第m个脉冲的预处理结果序列。
步骤S4：建立三维栅格Grid3d与通道的预处理结果序列Y c (In R )的索引信息集In c ，索引信息集In c 表示为
Figure DEST_PATH_IMAGE008
，其中，
Figure DEST_PATH_IMAGE009
为发射方位角度检索信息，
Figure DEST_PATH_IMAGE010
为接收方位角度检索信息，
Figure DEST_PATH_IMAGE011
为发射俯仰角检索信息，
Figure DEST_PATH_IMAGE012
为接收俯仰角检索信息，In R 为距离检索信息。
步骤S5：计算栅格的各接收通道的单通道检测统计量
Figure DEST_PATH_IMAGE013
，通道编号c∈[1,N], c∈Z。用C表示所有覆盖栅格中心的接收通道集合，则有∀c∈C。其中，
Figure DEST_PATH_IMAGE014
为检测单元噪声。将栅格的各单通道检测统计量
Figure 509481DEST_PATH_IMAGE013
进行累加，得到累加结果
Figure DEST_PATH_IMAGE015
。将累加结果K与检测器门限γt进行比较，若累加结果K大于γt，则判决结果为对应栅格存在目标，若累加结果K小于γt，则决结果为对应栅格不存在目标；
步骤S6：（1）将各有目标的栅格中各接收通道的单通道统计量
Figure 258738DEST_PATH_IMAGE013
与给定的单通道门限γc进行比较，统计各通道检测统计量大于单通道门限γc的栅格，组成门限栅格集合
Figure DEST_PATH_IMAGE016
，其中，Np为集合P中栅格的数量。
（2）统计门限栅格集合P中各栅格的单通道检测统计量大于单通道门限γc的接收通道个数，记为M，依据M/N逻辑判决准则，依次判决集合P中各栅格对应的M/N是否大于门限γc。保留M/N大于门限γc的栅格形成新的集合Pmn。
步骤S7：（1）对集合Pmn进行序贯凝聚，将序贯凝聚得到的点迹集合Psort按照累加统计量即累加结果K降序排序，将集合Psort中的第1个点迹存入集合Pclean，根据预处理结果序列Y c (l)的索引信息集In c 查找该点迹位置对应的各接收通道预处理数据，将查找出的对应数据Y c (In R )置为0，并从Psort中删除该点迹。
（2）如果Psort不为空集，则令Grid3d=Psort，并回到步骤S4；如果Psort为空集,则执行步骤S8。
步骤S8：以Pclean作为最终检测点迹集合。
上述检测方法只对探测区域进行了一次栅格划分，通过遍历各栅格区域完成目标检测。并且，为了保证较高的检测精度，上述检测方法中的栅格尺寸往往需要与雷达节点的距离分辨率接近。因此，在面对大探测场景时，上述检测方法不仅要遍历大量栅格，而且计算量庞大，难以满足高效能检测的发展需求。
为使本公开实施方式的目的、技术方案和优点更加清楚，下面将结合附图对本公开的各实施方式进行详细的阐述。然而，本领域的普通技术人员可以理解，在本公开各实施方式中，为了使读者更好地理解本公开而提出了许多技术细节。但是，即使没有这些技术细节和基于以下各实施方式的种种变化和修改，也可以实现本公开所要求保护的技术方案。以下各个实施方式的划分是为了描述方便，不应对本公开的具体实现方式构成任何限定，各个实施方式在不矛盾的前提下可以相互结合相互引用。
本公开的一个实施方式涉及一种基于迭代栅格的分布式雷达目标检测方法，其流程如图1所示，包括：
步骤101，离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与三维探测空间相对应的多个栅格组，其中，初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸，经度维尺寸、纬度维尺寸和高度维尺寸均远大于雷达距离分辨率。
示例性的，步骤101具体包括：
分别利用尺寸为(q1×q2×q3)/(k3)0,……,(q1×q2×q3)/(k3) n ,……,(q1×q2×q3)/(k3)U的长方体的栅格，对三维探测空间进行均匀划分，得到对应的栅格组grid0,……,grid n ,……,gridU。其中，q1代表初始栅格尺寸中栅格的经度维尺寸，q2代表初始栅格尺寸中栅格的纬度维尺寸，q3代表初始栅格尺寸中栅格的高度维尺寸，n=0,1,2,…,U表示栅格组的编号，U表示预设迭代次数，k表示预设迭代栅格尺寸比例系数。
需要说明的是，q1、q2、q3作为初始栅格尺寸中不同维度的尺寸，均远大于雷达距离分辨率。
通过利用远大于雷达距离分辨率的初始栅格尺寸将三维探测空间划分为尺寸依次减小的栅格组成的栅格组，可以使后续的栅格检测不断利用尺寸较小的栅格进行，从而有效保证了检测精度。
示例性的，在步骤101中，可以将不同尺寸栅格的经度坐标分别存储至矩阵J0、……、J n 、……、JU，将不同尺寸栅格的纬度坐标分别存储至矩阵V0、……、V n 、……、V U，将不同尺寸栅格的高度坐标分别存储至矩阵G0、……、G n 、……、G U。矩阵J0、……、J n 、……、JU可以表示为：
Figure DEST_PATH_IMAGE017
其中，矩阵J0、……、J n 、……、JU分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的经度坐标，i0、…、in、…、iU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的经度编号，lo0、…、lon、…、loU分别为栅格组grid0,……,grid n ,……,gridU中的经度编号的总个数，J i0,0、…、J in,n 、…、J iU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的经度编号i0、…、in、…、iU对应的经度坐标。
矩阵V0、……、V n 、……、V U可以表示为：
Figure DEST_PATH_IMAGE018
其中，矩阵V0、……、V n 、……、V U分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的纬度坐标，j0、…、jn、…、jU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的纬度编号，la0、…、lan、…、laU分别为栅格组grid0,……,grid n ,……,gridU中的纬度编号的总个数，V j0,0、…、V jn,n 、…、V jU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的纬度编号j0、…、jn、…、jU对应的纬度坐标。
矩阵G0、……、G n 、……、G U可以表示为：
Figure DEST_PATH_IMAGE019
其中，矩阵G0、……、G n 、……、G U分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的高度坐标，e0、…、en、…、eU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的高度编号，h0、…、hn、…、hU分别为栅格组grid0,……,grid n ,……,gridU中的高度编号的总个数，G e0,0、…、G en,n 、…、G eU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的高度编号e0、…、en、…、eU对应的高度坐标。
步骤102，分别计算各栅格组中各栅格与各雷达节点对应的方位、俯仰和距离。
步骤103，离线构建各栅格组与各通道距离单元对应的索引信息集。索引信息集包括栅格组中各栅格在对应通道下的索引信息，索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息。
示例性的，索引信息集表示为：
Figure DEST_PATH_IMAGE020
其中，I c0,……,I cn ,……,I cU分别表示栅格组grid0,……,grid n ,……,gridU对应的索引信息集，
Figure 916028DEST_PATH_IMAGE002
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射方位角度检索信息
Figure 988501DEST_PATH_IMAGE003
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收方位角度检索信息，
Figure 122810DEST_PATH_IMAGE004
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射俯仰角检索信息，
Figure 955768DEST_PATH_IMAGE005
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收俯仰角检索信息，I R0,……,I Rn ,……,I RU分别表示栅格组grid0,……,grid n ,……,gridU中栅格的距离检索信息。
步骤104，预处理各通道的回波数据，得到各通道的预处理结果序列。预处理包括匹配滤波处理和运动目标探测处理。
示例性的，预处理结果序列可以表示为：
Figure DEST_PATH_IMAGE021
。其中，S c (d)表示第c个通道中第d个距离单元的预处理结果序列。w为脉冲编号，W为脉冲个数。s c [d,w]为第c个通道中第d个距离单元的第w个脉冲的预处理结果序列。
步骤105，确定当前栅格组中存在目标的栅格：基于预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出当前栅格组中存在目标的栅格，其中，当前栅格组为与当前迭代次数相对应的栅格组。
示例性的，步骤105具体包括：
基于预处理结果序列以及与当前迭代次数对应的距离检索信息，分别计算当前栅格组中各栅格的各接收通道的单通道检测统计量；将同一栅格的各接收通道的单通道检测统计量进行累加，得到当前栅格组中各栅格分别对应的累加结果；将累加结果与预设检测器门限进行比较，选取累加结果大于预设检测器门限的栅格，作为当前栅格组中存在目标的栅格。
具体的，将当前栅格组记为grid u ，u表示当前迭代次数，则当前栅格组grid u 中各栅格的各接收通道的单通道检测统计量可以表示为
Figure DEST_PATH_IMAGE022
，其中，c∈[1, N]且c∈Z，N为接收通道个数，I Ru 为当前栅格组grid u 中栅格的距离检索信息，
Figure DEST_PATH_IMAGE023
为检测单元噪声。用C表示所有覆盖栅格中心的接收通道集合，则有∀c∈C。将当前栅格组grid u 中同一栅格的各接收通道的单通道检测统计量进行累加，得到当前栅格组grid u 中各栅格分别对应的累加结果K 1，则累加结果K 1可以表示为
Figure DEST_PATH_IMAGE024
。分别将当前栅格组grid u 中各栅格对应的累加结果K 1与预设检测器门限T进行比较，选取累加结果K 1大于预设检测器门限T的栅格，作为当前栅格组grid u 中存在目标的栅格。
步骤106，基于M/N准则进行逻辑判决：
将当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于预设单通道门限的栅格组成第一门限栅格集合。
具体的，将当前栅格组grid u 中存在目标的栅格的各接收通道的单通道检测统计量
Figure 705216DEST_PATH_IMAGE022
分别与预设单通道门限Tc进行比较，将单通道检测统计量
Figure 835459DEST_PATH_IMAGE022
大于预设单通道门限Tc的栅格分别记为p1,p2,…,pN2，并将p1,p2,…,pN2组成第一门限栅格集合P1，则有
Figure DEST_PATH_IMAGE025
，其中，N2为第一门限栅格集合P1中的栅格的数量。
将第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计第一门限栅格集合中各栅格对应的单通道检测统计量大于预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断第一门限栅格集合中各栅格对应的M/N是否大于预设单通道门限，将M/N大于预设单通道门限的栅格组成第二门限栅格集合。
具体的，将第一门限栅格集合P1中各栅格对应的接收通道的个数记为N，统计第一门限栅格集合P1中各栅格对应的单通道检测统计量大于预设单通道门限Tc的接收通道的个数，将统计出的接收通道的个数记为M，依据M/N逻辑判决准则，依次判断第一门限栅格集合P1中各栅格对应的M/N是否大于预设单通道门限Tc，将M/N大于预设单通道门限Tc的栅格组成第二门限栅格集合P2。
步骤107，序贯凝聚并镜像消除：对第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将点迹集合中的点迹按累加结果降序排序，将排序后的点迹集合中的第一个点迹存入最终点迹集合，并将第一个点迹从点迹集合中删除，得到新的点迹集合，其中，累加结果由与点迹相对应的栅格的各接收通道的单通道检测统计量累加得到。
具体的，对第二门限栅格集合P2进行序贯凝聚，得到对应的点迹集合Ps，将点迹集合Ps中的点迹按累加结果K 1降序排序，将排序后的点迹集合Ps中的第一个点迹存入最终点迹集合Pcl，并将该第一个点迹从点迹集合Ps中删除，得到新的点迹集合Ps。
示例性的，在将排序后的点迹集合中的第一个点迹存入最终点迹集合之后，所述方法还包括：根据当前栅格组对应的索引信息集查找第一个点迹的位置对应的各接收通道的预处理结果序列，将查找出的预处理结果序列置为0。通过将查找出的预处理结果序列置为0，可以有效减少后续迭代栅格检测过程中的计算量，从而提高检测效率。
具体的，在将排序后的点迹集合Ps中的第一个点迹存入最终点迹集合Pcl之后，还可以根据当前栅格组grid u 对应的索引信息集I cu 查找该第一个点迹的位置对应的各接收通道的预处理结果序列，将查找出的预处理结果序列S c (I Ru )置为0。
步骤108，判断镜像消除是否完成：判断新的点迹集合是否为空集：若新的点迹集合不为空集，则将新的点迹集合作为当前栅格组，回到确定当前栅格组中存在目标的栅格的步骤；若新的点迹集合为空集，则判断当前迭代次数是否达到预设迭代次数。
具体的，判断新的点迹集合Ps是否为空集：
若新的点迹集合Ps不为空集，则说明镜像消除未完成，将新的点迹集合Ps作为当前栅格组grid u 即令grid u =Ps，回到确定当前栅格组中存在目标的栅格的步骤即步骤105。示例性的，在新的点迹集合Ps不为空集时，还可以清除当前栅格组grid u 对应的索引信息集I cu 中不在新的点迹集合Ps中的栅格的索引信息。通过清除当前栅格组对应的索引信息集中不在新的点迹集合中的栅格的索引信息，可以有效减少后续迭代栅格检测过程中的计算量，从而提高检测效率。
若新的点迹集合Ps为空集，则说明镜像消除完成，执行步骤109。
步骤109，判断当前迭代次数是否达到预设迭代次数：若当前迭代次数未达到预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，回到确定当前栅格组中存在目标的栅格的步骤，其中，指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到预设迭代次数，则将最终点迹集合作为目标检测结果。
具体的，判断当前迭代次数u是否达到预设迭代次数U即u是否等于U：
若当前迭代次数u未达到预设迭代次数U即Ps为空集且u≠U，则将当前迭代次数加1即令u=u+1，清除当前栅格组grid u 中不在指定区域的栅格，回到确定当前栅格组中存在目标的栅格的步骤即步骤105。其中，指定区域根据上一个当前栅格组即grid u-1中存在目标的栅格确定。示例性的，在当前迭代次数u未达到预设迭代次数U即Ps为空集且u≠U时，还可以清除当前栅格组grid u 对应的索引信息集I cu 中不在指定区域的栅格的索引信息。通过清除当前栅格组grid u 对应的索引信息集I cu 中不在指定区域的栅格的索引信息，可以有效减少后续迭代栅格检测过程中的计算量，从而提高检测效率。
若当前迭代次数u达到预设迭代次数U即Ps为空集且u=U，则将最终点迹集合Pcl作为目标检测结果。
示例性的，指定区域可以通过以下步骤确定：
将最终点迹集合中各栅格的位置信息存储为位置集合，位置信息包括经度坐标、纬度坐标、高度坐标。
具体的，可以将最终点迹集合Pcl中各栅格的位置信息存储为位置集合Location，位置集合Location可以表示为Location=[(Jq,Vq,Gq) | q=1,2,…,Nloc]，其中，Jq为栅格q的经度坐标，Vq为栅格q的纬度坐标，Gq为栅格q的高度坐标，q=1,2,…,Nloc为最终点迹集合Pcl中的栅格编号，Nloc为最终点迹集合Pcl中的栅格数量。
基于位置集合、初始栅格尺寸、预设迭代栅格尺寸比例系数以及上一个当前迭代次数，在上一个当前栅格组中确定指定区域。
具体的，可以基于位置集合Location、初始栅格尺寸q1至q3、预设迭代栅格尺寸比例系数k以及上一个当前迭代次数u-1，在上一个当前栅格组grid u-1中确定指定区域。例如，指定区域可以为上一个当前栅格组grid u-1中经度处于[Jq-(3q1/2k u-1), Jq+(3q1/2k u-1)]范围、纬度处于[Vq-(3q2/2k u-1), Vq+(3q2/2k u-1)]范围、高度处于[Gq-(3q3/2k u-1), Gq+(3q3/2k u-1)]范围的区域，其中，u为当前迭代次数。
本公开实施方式相对于现有技术而言，离线建立多个尺寸的栅格以及栅格与各通道距离单元的索引关系，有效保证了迭代栅格检测的时效性；利用多个尺寸的栅格对三维探测空间进行迭代检测，在迭代检测过程中不断基于上一次迭代检测获得的存在目标的位置及其邻域范围缩小迭代检测范围直至达到预设迭代次数，有效提升了检测效率；并且，在迭代检测过程中，还可以不断利用尺寸更小的栅格在对应的迭代检测范围内进行检测，从而有效保证了检测精度。
为使本领域技术人员能够更好地理解上述实施方式，下面以一具体示例进行说明。
如图2所示，一种基于迭代栅格的分布式雷达目标检测方法，包括以下步骤：
步骤S01：将预设迭代次数记为U，预设迭代栅格尺寸比例系数记为k，变量u赋初值为0作为当前迭代次数。
步骤S02：离线构建三维空间栅格：分别利用尺寸为(q1×q2×q3)/(k3)0,……,(q1×q2×q3)/(k3) n ,……,(q1×q2×q3)/(k3)U的长方体的栅格，对三维探测空间进行均匀划分，得到对应的栅格组grid0,……,grid n ,……,gridU (n=0,1,2,…,U)。其中，q1代表初始栅格尺寸中栅格的经度维尺寸，q2代表初始栅格尺寸中栅格的纬度维尺寸，q3代表初始栅格尺寸中栅格的高度维尺寸，q1、q2、q3均远大于雷达距离分辨率，n=0,1,2,…,U表示栅格组的编号。使用矩阵J0、……、J n 、……、JU存储各尺寸栅格的经度坐标，使用矩阵V0、……、V n 、……、V U存储各尺寸栅格的纬度坐标，使用矩阵G0、……、G n 、……、G U存储各尺寸栅格的高度坐标。
矩阵J0、……、J n 、……、JU表示为：
Figure DEST_PATH_IMAGE026
其中，矩阵J0、……、J n 、……、JU分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的经度坐标，i0、…、in、…、iU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的经度编号，lo0、…、lon、…、loU分别为栅格组grid0,……,grid n ,……,gridU中的经度编号的总个数，J i0,0、…、J in,n 、…、J iU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的经度编号i0、…、in、…、iU对应的经度坐标。
矩阵V0、……、V n 、……、V U表示为：
Figure DEST_PATH_IMAGE027
其中，矩阵V0、……、V n 、……、V U分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的纬度坐标，j0、…、jn、…、jU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的纬度编号，la0、…、lan、…、laU分别为栅格组grid0,……,grid n ,……,gridU中的纬度编号的总个数，V j0,0、…、V jn,n 、…、V jU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的纬度编号j0、…、jn、…、jU对应的纬度坐标。
矩阵G0、……、G n 、……、G U表示为：
Figure DEST_PATH_IMAGE028
其中，矩阵G0、……、G n 、……、G U分别存储栅格组grid0,……,grid n ,……,gridU中各栅格的高度坐标，e0、…、en、…、eU分别为栅格组grid0,……,grid n ,……,gridU中各栅格的高度编号，h0、…、hn、…、hU分别为栅格组grid0,……,grid n ,……,gridU中的高度编号的总个数，G e0,0、…、G en,n 、…、G eU,U分别为栅格组grid0,……,grid n ,……,gridU中栅格的高度编号e0、…、en、…、eU对应的高度坐标。
步骤S03：分别计算各栅格组中各三维栅格与各雷达节点对应的方位、俯仰和距离。
步骤S04：离线构建各栅格组与各通道距离单元对应的索引信息集。即，离线构建各尺寸栅格对应的栅格组grid0,……,grid n ,……,gridU与数据距离单元的索引，得到对应的索引信息集I c0,……,I cn ,……,I cU，表示为：
Figure DEST_PATH_IMAGE029
其中，I c0,……,I cn ,……,I cU分别表示栅格组grid0,……,grid n ,……,gridU对应的索引信息集，
Figure 723692DEST_PATH_IMAGE002
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射方位角度检索信息，
Figure 306596DEST_PATH_IMAGE003
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收方位角度检索信息，
Figure 784982DEST_PATH_IMAGE004
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射俯仰角检索信息，
Figure 710344DEST_PATH_IMAGE005
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收俯仰角检索信息，I R0,……,I Rn ,……,I RU分别表示栅格组grid0,……,grid n ,……,gridU中栅格的距离检索信息。
步骤S05：预处理即匹配滤波处理和运动目标探测处理各通道的回波数据，得到各通道的预处理结果序列
Figure DEST_PATH_IMAGE030
。其中，S c (d)表示第c个通道中第d个距离单元的预处理结果序列。w为脉冲编号，W为脉冲个数。s c [d,w]为第c个通道中第d个距离单元的第w个脉冲的预处理结果序列。
步骤S06：计算当前栅格组grid u 中各栅格的单通道检测统计量，累加单通道检测统计量并进行累积量阈值检测：
计算当前栅格组grid u 中各栅格的各接收通道的单通道检测统计量
Figure DEST_PATH_IMAGE031
，其中，c∈[1, N]且c∈Z，N为接收通道个数，I Ru 为当前栅格组grid u 中栅格的距离检索信息，
Figure DEST_PATH_IMAGE032
为检测单元噪声。用C表示所有覆盖栅格中心的接收通道集合，则有∀c∈C。将当前栅格组grid u 中同一栅格的各接收通道的单通道检测统计量进行累加，得到当前栅格组grid u 中各栅格分别对应的累加结果
Figure DEST_PATH_IMAGE033
。分别将当前栅格组grid u 中各栅格对应的累加结果K 1与预设检测器门限T进行比较，选取累加结果K 1大于预设检测器门限T的栅格，作为当前栅格组grid u 中存在目标的栅格。
步骤S07：基于M/N准则进行逻辑判决：
（1）将当前栅格组grid u 中存在目标的栅格的各接收通道的单通道检测统计量
Figure 425271DEST_PATH_IMAGE031
分别与预设单通道门限Tc进行比较，统计单通道检测统计量
Figure 292208DEST_PATH_IMAGE031
大于预设单通道门限Tc的栅格，组成第一门限栅格集合
Figure DEST_PATH_IMAGE034
，其中，N2为第一门限栅格集合P1中的栅格的数量。
（2）统计第一门限栅格集合P1中各栅格对应的单通道检测统计量大于预设单通道门限Tc的接收通道的个数，记为M，依据M/N逻辑判决准则，依次判断第一门限栅格集合P1中各栅格对应的M/N是否大于预设单通道门限Tc，将M/N大于预设单通道门限Tc的栅格组成第二门限栅格集合P2。
步骤S08：序贯凝聚并镜像消除：
对第二门限栅格集合P2进行序贯凝聚，得到对应的点迹集合Ps，将点迹集合Ps中的点迹按累加结果K 1降序排序，将排序后的点迹集合Ps中的第一个点迹存入最终点迹集合Pcl，根据当前栅格组grid u 对应的索引信息集I cu 查找该第一个点迹的位置对应的各接收通道的预处理结果序列S c (I Ru )，将查找出的预处理结果序列S c (I Ru )置为0，并将该第一个点迹从点迹集合Ps中删除，得到新的点迹集合Ps。
步骤S09：判断镜像消除是否完成：
判断新的点迹集合Ps是否为空集：若新的点迹集合Ps不为空集，则令grid u =Ps，清除当前栅格组grid u 对应的索引信息集I cu 中不在新的点迹集合Ps中的栅格的索引信息，并回到步骤S06。若新的点迹集合Ps为空集，则执行步骤S010。
步骤S010：判断当前迭代次数u是否达到预设迭代次数U即u是否等于U：
若Ps为空集且u≠U，则：
（1）令u=u+1，统计最终点迹集合Pcl中各栅格的位置信息包括经度坐标、纬度坐标、高度坐标，存储为集合Location=[(Jq,Vq,Gq) | q=1,2,…,Nloc]，其中，Jq为栅格q的经度坐标，Vq为栅格q的纬度坐标，Gq为栅格q的高度坐标，q=1,2,…,Nloc为最终点迹集合Pcl中的栅格编号，Nloc为最终点迹集合Pcl中的栅格数量。
（2）基于当前的最终点迹集合Pcl中已存储的可能存在目标的栅格，参考位置集合Location存储的位置信息，依次在上一个当前栅格组grid u-1中查找经度处于[Jq-(3q1/2k u -1), Jq+(3q1/2k u-1)]范围、纬度处于[Vq-(3q2/2k u-1), Vq+(3q2/2k u-1)]范围、高度处于[Gq-(3q3/2k u-1), Gq+(3q3/2k u-1)]范围的所有区域作为指定区域，清除当前栅格组grid u 中不在指定区域的栅格，清除当前栅格组grid u 对应的索引信息集I cu 中不在指定区域的栅格的索引信息，并回到步骤S06。
若Ps为空集且u=U，则执行步骤S011。
步骤S011：将最终点迹集合Pcl作为目标检测结果，即，将最终点迹集合Pcl作为最终检测出的点迹集合并进行点迹输出。
本公开的另一个实施方式涉及一种基于迭代栅格的分布式雷达目标检测装置，如图3所示，包括：
第一构建模块301，用于离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与三维探测空间相对应的多个栅格组，其中，初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸，经度维尺寸、纬度维尺寸和高度维尺寸均远大于雷达距离分辨率；
计算模块302，用于分别计算各栅格组中各栅格与各雷达节点对应的方位、俯仰和距离；
第二构建模块303，用于离线构建各栅格组与各通道距离单元对应的索引信息集，索引信息集包括栅格组中各栅格在对应通道下的索引信息，索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息；
预处理模块304，用于预处理各通道的回波数据，得到各通道的预处理结果序列，预处理包括匹配滤波处理和运动目标探测处理；
确定模块305，用于确定当前栅格组中存在目标的栅格：基于预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出当前栅格组中存在目标的栅格，其中，当前栅格组为与当前迭代次数相对应的栅格组；
判决模块306，用于基于M/N准则进行逻辑判决：将当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于预设单通道门限的栅格组成第一门限栅格集合；将第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计第一门限栅格集合中各栅格对应的单通道检测统计量大于预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断第一门限栅格集合中各栅格对应的M/N是否大于预设单通道门限，将M/N大于预设单通道门限的栅格组成第二门限栅格集合；
消除模块307，用于序贯凝聚并镜像消除：对第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将点迹集合中的点迹按累加结果降序排序，将排序后的点迹集合中的第一个点迹存入最终点迹集合，并将第一个点迹从点迹集合中删除，得到新的点迹集合，其中，累加结果由与点迹相对应的栅格的各接收通道的单通道检测统计量累加得到；
第一判断模块308，用于判断镜像消除是否完成：判断新的点迹集合是否为空集：若新的点迹集合不为空集，则将新的点迹集合作为当前栅格组，重新触发确定模块305；若新的点迹集合为空集，则触发第二判断模块309；
第二判断模块309，用于判断当前迭代次数是否达到预设迭代次数：若当前迭代次数未达到预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，重新触发确定模块305，其中，指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到预设迭代次数，则将最终点迹集合作为目标检测结果。
本公开实施方式提供的基于迭代栅格的分布式雷达目标检测装置的具体实现方法，可以参见本公开实施方式提供的基于迭代栅格的分布式雷达目标检测方法所述，此处不再赘述。
本公开实施方式相对于现有技术而言，离线建立多个尺寸的栅格以及栅格与各通道距离单元的索引关系，有效保证了迭代栅格检测的时效性；利用多个尺寸的栅格对三维探测空间进行迭代检测，在迭代检测过程中不断基于上一次迭代检测获得的存在目标的位置及其邻域范围缩小迭代检测范围直至达到预设迭代次数，有效提升了检测效率；并且，在迭代检测过程中，还可以不断利用尺寸更小的栅格在对应的迭代检测范围内进行检测，从而有效保证了检测精度。
本公开的另一个实施方式涉及一种电子设备，如图4所示，包括：
至少一个处理器401；以及，
与至少一个处理器401通信连接的存储器402；其中，
存储器402存储有可被至少一个处理器401执行的指令，指令被至少一个处理器401执行，以使至少一个处理器401能够执行上述实施方式所述的基于迭代栅格的分布式雷达目标检测方法。
其中，存储器和处理器采用总线方式连接，总线可以包括任意数量的互联的总线和桥，总线将一个或多个处理器和存储器的各种电路连接在一起。总线还可以将诸如外围设备、稳压器和功率管理电路等之类的各种其他电路连接在一起，这些都是本领域所公知的，因此，本文不再对其进行进一步描述。总线接口在总线和收发机之间提供接口。收发机可以是一个元件，也可以是多个元件，比如多个接收器和发送器，提供用于在传输介质上与各种其他装置通信的单元。经处理器处理的数据通过天线在无线介质上进行传输，进一步，天线还接收数据并将数据传送给处理器。
处理器负责管理总线和通常的处理，还可以提供各种功能，包括定时，外围接口，电压调节、电源管理以及其他控制功能。而存储器可以被用于存储处理器在执行操作时所使用的数据。
本公开的另一个实施方式涉及一种计算机可读存储介质，存储有计算机程序，计算机程序被处理器执行时实现上述实施方式所述的基于迭代栅格的分布式雷达目标检测方法。
即，本领域技术人员可以理解，实现上述实施方式所述方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得一个设备（可以是单片机，芯片等）或处理器（processor）执行本公开各个实施方式所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（Read-Only Memory，ROM）、随机存取存储器（Random Access Memory，RAM）、磁碟或者光盘等各种可以存储程序代码的介质。
本领域的普通技术人员可以理解，上述各实施方式是实现本公开的具体实施方式，而在实际应用中，可以在形式上和细节上对其作各种改变，而不偏离本公开的精神和范围。

Claims (9)

1.一种基于迭代栅格的分布式雷达目标检测方法，其特征在于，所述方法包括以下步骤：
离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与所述三维探测空间相对应的多个栅格组，其中，所述初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸；
分别计算各所述栅格组中各栅格与各雷达节点对应的方位、俯仰和距离；
离线构建各所述栅格组与各通道距离单元对应的索引信息集，所述索引信息集包括所述栅格组中各栅格在对应通道下的索引信息，所述索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息；
预处理各通道的回波数据，得到各通道的预处理结果序列，所述预处理包括匹配滤波处理和运动目标探测处理；
确定当前栅格组中存在目标的栅格：基于所述预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出所述当前栅格组中存在目标的栅格，其中，所述当前栅格组为与当前迭代次数相对应的栅格组；
基于M/N准则进行逻辑判决：将所述当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于所述预设单通道门限的栅格组成第一门限栅格集合；将所述第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计所述第一门限栅格集合中各栅格对应的单通道检测统计量大于所述预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断所述第一门限栅格集合中各栅格对应的M/N是否大于所述预设单通道门限，将M/N大于所述预设单通道门限的栅格组成第二门限栅格集合；
序贯凝聚并镜像消除：对所述第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将所述点迹集合中的点迹按累加结果降序排序，将排序后的所述点迹集合中的第一个点迹存入最终点迹集合，并将所述第一个点迹从所述点迹集合中删除，得到新的点迹集合，其中，所述累加结果由与所述点迹相对应的栅格的各接收通道的单通道检测统计量累加得到；
判断镜像消除是否完成：判断所述新的点迹集合是否为空集：若所述新的点迹集合不为空集，则将所述新的点迹集合作为当前栅格组，回到确定当前栅格组中存在目标的栅格的步骤；若所述新的点迹集合为空集，则判断当前迭代次数是否达到所述预设迭代次数；
判断当前迭代次数是否达到所述预设迭代次数：若当前迭代次数未达到所述预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，回到确定当前栅格组中存在目标的栅格的步骤，其中，所述指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到所述预设迭代次数，则将最终点迹集合作为目标检测结果；
所述确定当前栅格组中存在目标的栅格，具体包括：
基于所述预处理结果序列以及与当前迭代次数对应的距离检索信息，分别计算所述当前栅格组中各栅格的各接收通道的单通道检测统计量；
将同一栅格的各接收通道的单通道检测统计量进行累加，得到所述当前栅格组中各栅格分别对应的累加结果；
将所述累加结果与预设检测器门限进行比较，选取所述累加结果大于所述预设检测器门限的栅格，作为所述当前栅格组中存在目标的栅格。
2.根据权利要求1所述的方法，其特征在于，所述离线构建三维空间栅格，具体包括：
分别利用尺寸为(q1×q2×q3)/(k3)0,……,(q1×q2×q3)/(k3) n ,……,(q1×q2×q3)/(k3)U的长方体的栅格，对所述三维探测空间进行均匀划分，得到对应的栅格组grid0,……,grid n ,……,gridU；
其中，q1代表所述初始栅格尺寸中栅格的经度维尺寸，q2代表所述初始栅格尺寸中栅格的纬度维尺寸，q3代表所述初始栅格尺寸中栅格的高度维尺寸，n=0,1,2,…,U表示栅格组的编号，U表示所述预设迭代次数，k表示所述预设迭代栅格尺寸比例系数。
3.根据权利要求2所述的方法，其特征在于，所述索引信息集表示为：
Figure 503781DEST_PATH_IMAGE001
其中，I c0,……,I cn ,……,I cU分别表示栅格组grid0,……,grid n ,……,gridU对应的索引信息集，
Figure DEST_PATH_IMAGE002
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射方位角度检索信息，
Figure DEST_PATH_IMAGE003
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收方位角度检索信息，
Figure DEST_PATH_IMAGE004
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的发射俯仰角检索信息，
Figure DEST_PATH_IMAGE005
分别表示栅格组grid0,……,grid n ,……,gridU中栅格的接收俯仰角检索信息，I R0,……,I Rn ,……,I RU分别表示栅格组grid0,……,grid n ,……,gridU中栅格的距离检索信息。
4.根据权利要求2所述的方法，其特征在于，所述指定区域通过以下步骤确定：
将最终点迹集合中各栅格的位置信息存储为位置集合，所述位置信息包括经度坐标、纬度坐标、高度坐标；
基于所述位置集合、所述初始栅格尺寸、所述预设迭代栅格尺寸比例系数以及上一个当前迭代次数，在上一个当前栅格组中确定所述指定区域。
5.根据权利要求4所述的方法，其特征在于，
所述位置集合表示为Location=[(Jq,Vq,Gq) | q=1,2,…,Nloc]，其中，Jq为栅格q的经度坐标，Vq为栅格q的纬度坐标，Gq为栅格q的高度坐标，q=1,2,…,Nloc为最终点迹集合中的栅格编号，Nloc为最终点迹集合中的栅格数量；
所述指定区域为上一个当前栅格组中经度处于[Jq-(3q1/2k u-1), Jq+(3q1/2k u-1)]范围、纬度处于[Vq-(3q2/2k u-1), Vq+(3q2/2k u-1)]范围、高度处于[Gq-(3q3/2k u-1), Gq+(3q3/2k u -1)]范围的区域，其中，u为当前迭代次数。
6.根据权利要求1至5任一项所述的方法，其特征在于，
在将排序后的所述点迹集合中的第一个点迹存入最终点迹集合之后，所述方法还包括：根据当前栅格组对应的索引信息集查找所述第一个点迹的位置对应的各接收通道的预处理结果序列，将查找出的预处理结果序列置为0；
在所述新的点迹集合不为空集时，所述方法还包括：清除当前栅格组对应的索引信息集中不在所述新的点迹集合中的栅格的索引信息；
在当前迭代次数未达到所述预设迭代次数时，所述方法还包括：清除当前栅格组对应的索引信息集中不在所述指定区域的栅格的索引信息。
7.一种基于迭代栅格的分布式雷达目标检测装置，其特征在于，所述装置包括：
第一构建模块，用于离线构建三维空间栅格：分别利用基于初始栅格尺寸、预设迭代栅格尺寸比例系数和预设迭代次数确定的不同尺寸长方体的栅格，对三维探测空间进行均匀划分，得到与所述三维探测空间相对应的多个栅格组，其中，所述初始栅格尺寸包括栅格的经度维尺寸、纬度维尺寸和高度维尺寸；
计算模块，用于分别计算各所述栅格组中各栅格与各雷达节点对应的方位、俯仰和距离；
第二构建模块，用于离线构建各所述栅格组与各通道距离单元对应的索引信息集，所述索引信息集包括所述栅格组中各栅格在对应通道下的索引信息，所述索引信息包括发射方位角度检索信息、接收方位角度检索信息、发射俯仰角检索信息、接收俯仰角检索信息、距离检索信息；
预处理模块，用于预处理各通道的回波数据，得到各通道的预处理结果序列，所述预处理包括匹配滤波处理和运动目标探测处理；
确定模块，用于确定当前栅格组中存在目标的栅格：基于所述预处理结果序列以及与当前迭代次数对应的距离检索信息，在当前栅格组进行栅格检测，确定出所述当前栅格组中存在目标的栅格，其中，所述当前栅格组为与当前迭代次数相对应的栅格组；
判决模块，用于基于M/N准则进行逻辑判决：将所述当前栅格组中存在目标的栅格的各接收通道的单通道检测统计量分别与预设单通道门限进行比较，将单通道检测统计量大于所述预设单通道门限的栅格组成第一门限栅格集合；将所述第一门限栅格集合中各栅格对应的接收通道的个数记为N，统计所述第一门限栅格集合中各栅格对应的单通道检测统计量大于所述预设单通道门限的接收通道的个数M，依据M/N逻辑判决准则，依次判断所述第一门限栅格集合中各栅格对应的M/N是否大于所述预设单通道门限，将M/N大于所述预设单通道门限的栅格组成第二门限栅格集合；
消除模块，用于序贯凝聚并镜像消除：对所述第二门限栅格集合进行序贯凝聚，得到对应的点迹集合，将所述点迹集合中的点迹按累加结果降序排序，将排序后的所述点迹集合中的第一个点迹存入最终点迹集合，并将所述第一个点迹从所述点迹集合中删除，得到新的点迹集合，其中，所述累加结果由与所述点迹相对应的栅格的各接收通道的单通道检测统计量累加得到；
第一判断模块，用于判断镜像消除是否完成：判断所述新的点迹集合是否为空集：若所述新的点迹集合不为空集，则将所述新的点迹集合作为当前栅格组，重新触发所述确定模块；若所述新的点迹集合为空集，则触发第二判断模块；
第二判断模块，用于判断当前迭代次数是否达到所述预设迭代次数：若当前迭代次数未达到所述预设迭代次数，则将当前迭代次数加1，清除当前栅格组中不在指定区域的栅格，重新触发所述确定模块，其中，所述指定区域根据上一个当前栅格组中存在目标的栅格确定；若当前迭代次数达到所述预设迭代次数，则将最终点迹集合作为目标检测结果；
所述确定当前栅格组中存在目标的栅格，具体包括：
基于所述预处理结果序列以及与当前迭代次数对应的距离检索信息，分别计算所述当前栅格组中各栅格的各接收通道的单通道检测统计量；
将同一栅格的各接收通道的单通道检测统计量进行累加，得到所述当前栅格组中各栅格分别对应的累加结果；
将所述累加结果与预设检测器门限进行比较，选取所述累加结果大于所述预设检测器门限的栅格，作为所述当前栅格组中存在目标的栅格。
8.一种电子设备，其特征在于，包括：
至少一个处理器；以及，
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1至6中任一项所述的基于迭代栅格的分布式雷达目标检测方法。
9.一种计算机可读存储介质，存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至6中任一项所述的基于迭代栅格的分布式雷达目标检测方法。
CN202210712106.9A 2022-06-22 2022-06-22 基于迭代栅格的分布式雷达目标检测方法及装置 Active CN114779225B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210712106.9A CN114779225B (zh) 2022-06-22 2022-06-22 基于迭代栅格的分布式雷达目标检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210712106.9A CN114779225B (zh) 2022-06-22 2022-06-22 基于迭代栅格的分布式雷达目标检测方法及装置

Publications (2)

Publication Number Publication Date
CN114779225A CN114779225A (zh) 2022-07-22
CN114779225B true CN114779225B (zh) 2022-09-13

Family

ID=82422366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210712106.9A Active CN114779225B (zh) 2022-06-22 2022-06-22 基于迭代栅格的分布式雷达目标检测方法及装置

Country Status (1)

Country Link
CN (1) CN114779225B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383339A (zh) * 2016-08-30 2017-02-08 电子科技大学 一种多站点雷达信号级联合检测的镜像目标抑制方法
CN107015205A (zh) * 2017-03-15 2017-08-04 电子科技大学 一种分布式mimo雷达检测的虚假目标消除方法
CN113569446A (zh) * 2021-07-02 2021-10-29 上海机电工程研究所 基于三维空间网格离散化的防空雷达优化部署方法及系统
DE102021003786A1 (de) * 2021-07-23 2021-11-04 Daimler Ag Verfahren zum Erzeugen einer Radarrasterkarte für ein Kraftfahrzeug mittels einer Radarvorrichtung, sowie Radarvorrichtung
CN113740849A (zh) * 2021-10-29 2021-12-03 中国电子科技集团公司信息科学研究院 一种多雷达自组织协同探测系统及方法
CN113835078A (zh) * 2021-11-30 2021-12-24 中国电子科技集团公司信息科学研究院 基于局部三维栅格的信号级联合检测方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113267754B (zh) * 2021-05-13 2023-03-03 中国人民解放军军事科学院战争研究院 一种基于三维网格的地形遮挡雷达探测范围计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383339A (zh) * 2016-08-30 2017-02-08 电子科技大学 一种多站点雷达信号级联合检测的镜像目标抑制方法
CN107015205A (zh) * 2017-03-15 2017-08-04 电子科技大学 一种分布式mimo雷达检测的虚假目标消除方法
CN113569446A (zh) * 2021-07-02 2021-10-29 上海机电工程研究所 基于三维空间网格离散化的防空雷达优化部署方法及系统
DE102021003786A1 (de) * 2021-07-23 2021-11-04 Daimler Ag Verfahren zum Erzeugen einer Radarrasterkarte für ein Kraftfahrzeug mittels einer Radarvorrichtung, sowie Radarvorrichtung
CN113740849A (zh) * 2021-10-29 2021-12-03 中国电子科技集团公司信息科学研究院 一种多雷达自组织协同探测系统及方法
CN113835078A (zh) * 2021-11-30 2021-12-24 中国电子科技集团公司信息科学研究院 基于局部三维栅格的信号级联合检测方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Grid cell based detection strategy for MIMO radar with widely separated subarrays";Chao, SY 等;《AEU-INTERNATIONAL JOURNAL OF ELECTRONICS AND COMMUNICATIONS》;20121230;第66卷(第9期);正文第741-751页 *
"Radar Target Detection via GAMP: A Sparse Recovery Strategy Off the Grid";Benzhou Jin 等;《IEEE Transactions on Vehicular Technology》;20211230;第70卷(第5期);正文第4153-4165页 *
"基于空间栅格划分的协同探测系统检测算法研究";李帅;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180215(第2期);正文第20-59页 *

Also Published As

Publication number Publication date
CN114779225A (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
CN110991311B (zh) 一种基于密集连接深度网络的目标检测方法
CN111327377B (zh) 场强预测方法、装置、设备及存储介质
CN107038717A (zh) 一种基于立体栅格自动分析3d点云配准误差的方法
CN108764168B (zh) 用于成像卫星在多障碍物海面搜索移动目标的方法及系统
CN109557532B (zh) 基于三维霍夫变换的检测前跟踪方法、雷达目标检测系统
JP7417745B2 (ja) 小包を検出するための方法、装置、計算装置、物流システム及び記憶媒体
CN101794437B (zh) 一种高光谱遥感影像异常探测方法
CN113835078B (zh) 基于局部三维栅格的信号级联合检测方法及装置
CN112070870B (zh) 点云地图评估方法、装置、计算机设备和存储介质
CN110795978B (zh) 路面点云数据提取方法、装置、存储介质及电子设备
CN108562867B (zh) 一种基于聚类的指纹定位方法及装置
CN114511679A (zh) 点云数据处理方法、装置、设备及存储介质
CN108803659B (zh) 基于魔方模型的多窗口启发式三维空间路径规划方法
CN114814807B (zh) 基于动态多尺度栅格的分布式雷达目标检测方法及装置
CN115410104A (zh) 一种获取飞行器的图像采集点的数据处理系统
KR101438377B1 (ko) 이동체의 위치 판단 장치 및 위치 판단 방법
CN110769369B (zh) 位置判断方法、判断装置、服务器和计算机可读存储介质
CN114386466B (zh) 一种用于脉冲星搜寻中候选体信号挖掘的并行的混合聚类方法
CN114529615A (zh) 雷达标定方法、装置及存储介质
CN114779225B (zh) 基于迭代栅格的分布式雷达目标检测方法及装置
CN108680897A (zh) 一种室内定位方法、装置、电子设备及存储介质
CN112097772B (zh) 机器人及其地图构建方法和装置
CN117705119A (zh) 大场景3d地图动态加载的导航定位方法、设备及存储介质
KR102260556B1 (ko) 전역 정보와 지역 정보 통합을 통한 딥러닝 기반 주차구획 검출 방법 및 장치
CN103092897A (zh) 一种用于点云数据处理的快速k近邻搜索方法

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