基于剖分网格的雷达探测能力计算方法、装置及设备

1.本技术涉及雷达探测技术领域，特别是涉及一种基于剖分网格的雷达探测能力计算方法、装置及设备。


背景技术：

2.雷达通过发射电磁波对目标进行照射并接受其回波来获得目标的位置信息。传统雷达探测能力计算的研究重点是水平方向和高度方向的二维模拟，这样的模型可视化雷达探测区域很难真实、高效的展示雷达的探测范围。
3.随着计算机性能的提升，并行化技术的快速发展，雷达三维探测能力可视化已成为研究热点。传统的雷达探测能力三维建模通过分析雷达探测范围边界与数字地形高程的相对位置关系构建雷达探测范围的三维模型。这样的处理流程算法复杂，计算量大，当雷达的参数发生变化时，模型很难实现快速响应。因此，急需提出一种三维空间下雷达探测能力的快速计算方法来进行有效的建模。
4.随着地理信息采集技术的发展，地理数据的规模不断扩大，要求gis算法能够快速访问较大规模的数据，地理要素数据作为一种重要的地理信息类型，在国土、地质、水文等领域起到关键作用，这些应用迫切需要提高地理要素数据访问的效率。


技术实现要素：

5.基于此，有必要针对上述技术问题，提供一种能够真实、高效的展示雷达探测范围的基于剖分网格的雷达探测能力计算方法、装置及设备。
6.一种基于剖分网格的雷达探测能力计算方法，所述方法包括：获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
7.在其中一实施例中，在计算各三维网格节点上的雷达功率密度之前，还包括：根据雷达的类型以及相关参数构建0-360
°
各角度对应的方向性系数对照表；在计算各三维节点上的雷达功率密度时，根据各三维网格节点相对于雷达的不同角度直接调用所述方向性系数对照表进行查表得到对应的方向性系数。
8.在其中一实施例中，在计算各三维网格节点上的雷达功率密度时：根据各三维网格节点的位置坐标以及雷达的位置计算得到该雷达针对各所述三
维网格节点的方位角、俯仰角以及距离；根据所述方位角、俯仰角以及距离进行计算得到对应三维网格节点的雷达功率密度。
9.在其中一实施例中，在计算各三维网格节点上的雷达功率密度时，采用以下公式：其中，式中，表示天线增益，表示在方位角为以及俯仰角为下雷达功率密度，表示雷达的发射功率，表示雷达的天线单元与目标三维网格节点之间的距离，表示距离r处的大气损耗，表示在方位角为以及俯仰角为下天线增益，为方向性系数，可由所述方向性系数对照表查表得到，表示方位角方向的方向性系数，表示俯仰角方向的方向性系数。
10.在其中一实施例中，当雷达为可移动雷达，则在计算处于该雷达所能作用的三维空间周边的边缘三维网格节点的雷达功率密度时：获取雷达的初始姿态以及移动范围；在所述移动范围内针对各所述边缘三维网格节点分别计算最优的雷达姿态；根据最优的雷达姿态时雷达所在的位置、方位角度以及俯仰角度与对应的边缘三维网格节点的位置坐标进行计算得到边缘三维网格节点的雷达功率密度。
11.在其中一实施例中，所述方法还包括：获取雷达的散射截面积；根据所述数据库中各三维网格节点的雷达功率密度以及散射截面积进行计算，并根据计算结果判断各三维网格节点是否能被检测到，以此得到雷达的感知范围；根据所述感知范围确定该雷达的探测能力。
12.在其中一实施例中，所述雷达探测能力计算方法适用雷达类型包括：缝阵雷达、相控阵雷达以及八木天线雷达。
13.在其中一实施例中，所述数据库采用clickhouse。
14.本技术还提供了一种基于剖分网格的雷达探测能力计算装置，所述装置包括：雷达作用三维空间得到模块，用于获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；全球三维网格划分模块，用于根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；雷达功率密度并行计算模块，用于将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网
格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；雷达探测能力确定模块，用于根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
15.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
16.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
17.上述基于剖分网格的雷达探测能力计算方法、装置及设备，通过采用geosot-3d全球剖分网格编码对雷达的所能作用范围进行剖分得到多个三维网格，这样能从低精度到高精度的对雷达所能作用的空间进行划分，再计算雷达在各网格区域的功率密度以真实展现雷达的探测范围，而在进行雷达功率密度的计算时，通过将多个三维网格划分为多个计算区域，再分配给对应进程进行并行计算，这样可高效的计算出所有三维网格区域内的功率密度，以达到快速计算雷达探测能力的目的。
附图说明
18.图1为一个实施例中雷达探测能力计算方法的流程示意图；图2为一个实施例中geosot-3d网络多级剖分示意图：（a）为一级剖分示意图，（b）为二级剖分示意图，（c）为三级剖分示意图；图3为一个实施例中按照空间分配计算任务情况示意图；
图4为一个实施例中并行算法计算时间与加速比的示意图：（a）为并行算法计算时间的示意图，（b）为并行算法计算加速比的示意图；图5为一个实施例中基于雷达探测能力计算方法的系统架构示意图；图6为一个实施例中雷达探测能力计算方法的模型示意图图7为一个实施例中雷达探测能力计算装置的结构框图；图8为一个实施例中计算机设备的内部结构图。
具体实施方式
19.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
20.如图1所示，提供了一种基于剖分网格的雷达探测能力计算方法，包括以下步骤：步骤s100，获取雷达的位置，根据雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；步骤s110，根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；步骤s120，将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；步骤s130，根据数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
21.在本实施例中，为了计算雷达在作用范围内的探测能力，采用了geosot-3d全球剖分网格编码对雷达作用范围空间进行剖分，也就是将雷达作用范围空间进行分割以多个三维立方体网格的形式进行表达，并且通过划分粒度的设置可确定各三维网格的大小，也就是说当划分粒度越小时，雷达探测范围被划分为更细小的网格，计算精度越高，可根据实际情况进行调节。而雷达的探测能力是通过计算各三维网格区域内的功率密度确定的，当完成雷达作用范围空间内所有三维网格的功率密度计算后，可得到该空间内雷达的功率密度分布以此来判断雷达的探测能力。而在进行功率密度计算时，本方法中采用了并行计算方法以提高计算效率。
22.在这里简要对geosot-3d全球剖分网格编码进行介绍：geosot-3d编码源自北京大学提出的基于2n一维整型数组地理坐标的地球剖分网格geosot，geosot采用经纬度剖分，具有边界不重叠、网格正交、经纬一致、与传统数据规格兼容性好等特点，是一种地球表面的剖分与编码方法。
23.geosot-3d在继承了这些特性的基础上，将geosot二维网格在高度维上进行扩展，将对自地心至地球外围50000km的地球立体空间，沿经线、纬线和高程方向进行八叉树剖分。为了使得剖分体元大小为二进制整型， geosot-3d将地球空间扩展为512
°×
512
°×
512
°
，将1
°
扩展为64
′
，将1
′
扩展为64
″
，实现整度、整分、整秒的递归八叉树剖分。geosot-3d
网络多级剖分示意图2所示，图2中（a）为一级剖分示意图，（b）为二级剖分示意图，（c）为三级剖分示意图，其中，图2（b）中编号0-7分别表示二级剖分形成的第0-7级等分体块，图2（c）中编号01-07分别表示第0级等分体块经过三级剖分形成的第01-07级等分体块。
24.geosot-3d立体网格分为32级，多级剖分后，将三维空间剖分成不同粒度的立体网格，第1-9级是度级网格，第10-15级网格为分级网格、第16-21级为秒级网格，第22-32级为秒以下网格。
25.平面投影类地图不能表达三维数据，可视化效果也不好，因此在本方法中使用geosot-3d对雷达作用范围进行划分。geosot-3d具有全球尺度、不同粒度、支持真三维空间剖分的特点，很适合对雷达作用范围进行表达。
26.在步骤s100中，首先获取雷达的位置，在这里已经预设了雷达类型，以及该雷达的相关参数。本方法可以适用于多种雷达包括但不限于：缝阵雷达、相控阵雷达以及八木天线雷达。
27.而与雷达相关的参数包括雷达阵元个数、阵元间距或y方向尺寸、工作频率。
28.在得到雷达位置后，根据雷达中心点以及雷达的作用范围寻找到雷达所能作用的三维空间。
29.在步骤s110中，根据预设的geosot-3d的划分粒度将雷达作用空间剖分为三维网格，剖分之后的格点是具有经度、纬度、高程位置信息的真三维节点。节点与geosot-3d编码一一对应，将雷达作用空间数据离散化，统一了空间编码，方便了雷达功率密度在空间上的表达。
30.在对各三维网格节点进行计算时，若给定经纬度与层级计算剖分网格编码，则首先根据层级确定经纬度所在网格的经度、纬度、高程序号，然后按照geosot编码规则，考虑10、16、22级的跳进，得出在一个方向上的编码。最后按照morton码的编码规则，将各维度编码组合成三维geosot剖分编码。
31.若给定剖分网格编码确认经纬度与高程，则首先按照morton码的规则对编码进行反解，得出各维度的编码，然后根据编码确认网格的序号，最后得出经纬度与高程坐标，也就是各三维网格节点的位置坐标以及与各位置坐标一一对应的网格编码。
32.在步骤s120中，计算各三维网格区域的雷达功率密度时，根据各三维网格节点的位置坐标以及雷达的位置计算得到该雷达针对各三维网格节点的方位角、俯仰角以及距离，并根据方位角、俯仰角以及距离进行计算得到对应三维网格节点的雷达功率密度。
33.在本实施例中，在辐射强度最大的方向上，考虑对特定方向的辐射功率的影响，为方位角，为俯仰角，距离r处的功率密度（峰值发射功率密度）：（1）其中，（2）在公式（1）中，表示天线增益，表示在方位角为以及俯仰角为
下雷达功率密度，表示雷达的发射功率，表示雷达的天线单元与目标三维网格节点之间的距离，表示距离r处的大气损耗，表示在方位角为以及俯仰角为下天线增益；在公式（2）中，为方向性系数，表示方位角方向的方向性系数，表示俯仰角方向的方向性系数。
34.通常情况下，雷达天线的视线方向为天线增益的最大方向，天线最大增益。也就是方位角和俯仰角均为0的时候。
35.而对于方向性系数的计算，根据雷达天线类别不同，选择的计算公式也是不同，对于孔径天线采用：（3）在公式（3）中，为y方向矩形孔径大小。
36.对于陈列天线采用：（4）在对雷达功率密度的计算过程中最耗时的部分是关于sin的计算，因此，为了加快计算效率，在进行预处理阶段先构建方向性系数对照表，并将该表保存在内存中，在进行计算时直接通过查表的方式得到方向性系数，以减少计算时间。
37.在本实施例中，构建方向性系数对照表时，根据雷达的类型以及相关参数构建0-360
°
各角度对应的方向性系数对照表。
38.具体的，根据雷达的类型以及相关参数依据不同的方向性系数计算公式，例如公式（3）-（4）计算出0-360
°
各角度对应的方向性系数，且精度值可调，在一些实施例中，精度值为。这样在计算各三维节点上的雷达功率密度时，只需要根据各三维网格节点相对于雷达的不同角度直接调用方向性系数对照表依照角度进行查表得到对应的方向性系数，不需要再次进行计算，节省了大量的计算时间。
39.因为雷达作用空间范围中可能划分有上百万个三维网格，也就意味着需要针对这百万个三维网格均要进行一次计算，其中计算较为耗时的sin计算，也就是方向性系数的计算在本实施例中采用了查表的方式以减少计算时间，在本方法中为了进一步的提高计算效率，在对各三维网格区域的雷达功率密度进行计算时，采用了mpi（high performance computing,信息传递接口）并行计算框架。
40.mpi并行计算框架的优点是多进程并行，与openmp等共享内存的计算框架相比，进程空间独立，可以有效避免内存与缓存的竞争，在无明显io与数据交换的前提下实现近线性加速比。同时，mpi算法可以直接跨节点运行。
41.在本实施例中，算法采用几个措施避免mpi进程的启动开销与进程间的数据交互。
42.首先根据算法主要使用场景，在预处理阶段进行方向性系数对照表的计算，且各进程进行全量计算，使得在载具位置（也就是雷达位置）与姿态更新时，可以以最小的代价实现计算结果的更新。
43.其次，按照空间网格对计算任务进行分割，各节点分别计算，不进行交互。
44.最后，采用clickhouse作为数据库，利用clickhouse的并发优势，每个进程计算完成后直接写入，避免数据汇总带来的通信与同步开销。
45.具体的，在对各三维网格节点的雷达功率密度进行计算时，按照空间将各三维网格节点划分为多个计算区域，也就是多个计算任务，并分配到不同的进程，由多个进程并行计算各自的分配到的计算任务。
46.以四个进程为例，空间上划分计算区域的情况如图3所示，不同的深浅表示不同进程被分配到的计算区域，而各计算区域内部由预设给定的划分粒度划分了对应的更加细小的网络，不同进程的计算互不干扰，在各进程计算结束后向0号进程发送一个确认信息，表示已经计算完成。
47.在进行计算时，部分网格计算快捷，因此在各进程内部，首先全量执行划分网络的算法，以雷达为中心构建geosot-3d网格。然后对各自所分配到的网格进行并行计算，结合clickhouse并发写入的特点，在计算出三维网格对应的结果后直接写入数据库，节省了进程通信的时间。当所有进程执行完毕，计算结果全部写入clickhouse后，0号进程收到确认信息，结束并行算法。
48.如图4所示，（a）为并行算法计算时间的示意图，（b）为并行算法计算加速比的示意图，可见计算时间随进程数不断减小，而加速度比几乎与进程数一致，即获得了近线性的加速比。
49.如图5所示，结合上述方法，本技术还提供了一种相适应的雷达探测能力计算系统，系统根据输入的雷达位置和三维网格参数相应的调用网格划分接口和雷达功率密度计算接口，采用多进程并行的方式计算个网格的结果数据，并直接写入数据库clickhouse供查询和验证。
50.从图5可以看出，通用雷达服务是在服务器启动时提供的雷达探测能力的计算接口。其中预设了特定雷达的参数，为特定雷达提供了各自的计算接口，即为某型雷达服务。载具的参数代表了雷达的初始位置和姿态，用户开启一个雷达的客户端后，可以根据雷达自身参数进行位置及姿态的调整，可计算调整后的新的雷达探测范围及在空间网格内的功率密度。
51.在本实施例中，当雷达为固定雷达时，针对各三维网格节点的位置，方位角及俯仰角度是不变的，也就是说各节点上雷达频率密度是不会改变。而当雷达为可移动时，对各三维网格节点发射探测信号时，由于雷达姿态的不同针对同一个三维网格节点计算得到的雷达功率密度也是不相同的。而实际上，在计算雷达功率密度时，希望可以在雷达辐射强度最大的方向上进行计算，此时雷达天线增益为最大g（0,0），也就是雷达与三维网格节点之间的方位角与俯仰角均为零的时候。这样处在雷达移动范围内的三维网格节点均可在雷达天线最大时计算雷达功率密度。
52.而针对雷达所能作用的三维空间周边的边缘网格节点的雷达功率密度，需要获取
雷达的初始姿态以及移动范围，在移动范围内针对各边缘三维网格节点分别计算最优的雷达姿态。根据最优的雷达姿态时雷达所在的位置、方位角度以及俯仰角度与对应的边缘三维网格节点的位置坐标进行计算得到边缘三维网格节点的雷达功率密度。
53.具体的，针对各边缘三维网格节点，最优的雷达姿态即为雷达与边缘三维网格节点之间呈现最小角度时。
54.在本实施例中，若可以获取到雷达散射截面积（rcs），还可以根据数据库中各三维网格节点的雷达功率密度以及散射截面积进行计算，并根据计算结果判断各三维网格节点是否能被检测到，以此得到雷达的感知范围，根据感知范围确定该雷达的探测能力。
55.如图6所示，本技术还给出了一种本方法中雷达频率密度计算模型。
56.上述基于剖分网格的雷达探测能力计算方法中，结合了geosot-3d的剖分方法，将雷达探测范围离散化，然后对各三维网格节点采用并行化的方法快速高效地计算出雷达在该点的功率密度，建立了一个基于三维剖分网格的，当雷达位置、方向变化时，能够及时响应的快速雷达探测能力计算系统。实验结果显示，采取12进程并行运行时，计算雷达作用范围内100万个位置点的雷达功率密度的平均总耗时约为40ms。
57.本技术提供了一套耦合性低、易用性强，支持多种类型雷达，能根据载具和雷达变化及时响应的雷达探测能力快速计算接口。
58.应该理解的是，虽然图1流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
59.在一个实施例中，如图7所示，提供了一种基于剖分网格的雷达探测能力计算装置，包括：雷达位置获取模块200、全球三维网格划分模块210、雷达功率密度并行计算模块220以及雷达探测能力确定模块230，其中：雷达位置获取模块200，用于获取雷达的位置；全球三维网格划分模块210，用于根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间，根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；雷达功率密度并行计算模块220，用于将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；雷达探测能力确定模块230，用于根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
60.关于基于剖分网格的雷达探测能力计算装置的具体限定可以参见上文中对于基于剖分网格的雷达探测能力计算方法的限定，在此不再赘述。上述基于剖分网格的雷达探测能力计算装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块
可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
61.在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储各三维网格节点的雷达功率密度数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于剖分网格的雷达探测能力计算方法。
62.本领域技术人员可以理解，图8中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
63.在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
64.在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：获取雷达的位置，根据所述雷达位置的中心点及作用范围得到该雷达所能作用的三维空间；根据geosot-3d在预设的划分粒度下将所述三维空间剖分为多个三维网格，并计算各三维网格节点的位置坐标以及网格编码；将各三维网格节点按照空间位置划分为预设数量的区域，将各区域分别匹配至不同的进程，由多个进程并行计算匹配区域内所有三维网格节点上的雷达功率密度，并在各进程完成计算后将各三维网格节点上的计算结果以及对应网格编码直接写入数据库中；根据所述数据库中的所有三维网格节点的雷达功率密度以及对应的网格编码得到所述雷达作用范围内的功率密度分布，根据所述功率密度分布确定该雷达的探测能力。
65.本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可
包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器（rom）、可编程rom（prom）、电可编程rom（eprom）、电可擦除可编程rom（eeprom）或闪存。易失性存储器可包括随机存取存储器（ram）或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram（sram）、动态ram（dram）、同步dram（sdram）、双数据率sdram（ddrsdram）、增强型sdram（esdram）、同步链路（synchlink） dram（sldram）、存储器总线（rambus）直接ram（rdram）、直接存储器总线动态ram（drdram）、以及存储器总线动态ram（rdram）等。
66.以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
67.以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。