[发明专利]毫米波雷达反射率逐层修复方法、系统及计算机存储介质

权利要求书

1.一种毫米波雷达反射率逐层修复方法，其特征在于，包括：

步骤S1、毫米波雷达反射率网格化处理；具体包括：

确定毫米波雷达的探测地点，收集探测过程中毫米波雷达探测的一次反射率数据R，根据R所对应的空间探测范围，设定分辨率为m×m的矩形网格，并将R按照设定的分辨率划分为x×y个网格，x为东西向的网格数，y为南北向的网格数；并统计网格化后的反射率数据R的最大值Rmax和最小值Rmin，并划分r段，每段跨度为：ΔR＝(Rmax-Rmin)/r；

步骤S2、窗口网格和抽吸层数设置；具体包括：

将反射率数据R记为第0层，设定整个修补过程的数据抽吸的层数L，并设定向下各层抽吸的窗口网格，所设置的窗口网格东西向和南北向至少有一个包含上一层大于或等于2的整数倍单位网格、且所设置的窗口网格在横移和竖移情况下都能均分上一层的全部单位网格；

步骤S3、毫米波雷达反射率逐层抽吸；具体包括：

根据各层相应的窗口网格对反射率数据R进行逐层抽吸以完成L层的抽吸计算；每抽吸一次，计算当前窗口网格的数值，分以下两种情况：

A、当窗口网格中没有缺失数据，则抽吸后的数据为窗口中所有值的平均值；

B、当抽吸窗口有缺失数据，则该抽吸计算后的网格也为缺失值；

步骤S4、L层雷达反射率数据蔓延补全；具体包括：

统计L层抽吸结果中缺失值个数为n、每一个缺失值周围网格数G和有效值个数I，并计算有效值比例为p＝I/G；按照p的值从高到低，排列数据补全的顺序，当两个或多个缺失网格对应的p值相同时，则按照缺失网格位置，从上到下，从左到右排序，记录排列后的顺序为p1,p2…pn；

对于p1对应的缺失网格，按照步骤S1中设定的反射率分段，计算其周围网格中有效值数据值的频率；当有效值频率中存在唯一最大值，则计算最大频率段中，所有值的平均值，并将该平均值作为缺失值的数值；当有效值频率中存在多个最大值，则计算多个最大频率段中，所有值的平均值，并将该平均值作为缺失值的数值；

补全p1对应的缺失网格后，重新统计L层中缺失值周围有效值的比例，并重新排序找出有效值比例最高的点，按照上述策略补全；以此类推，直到L层上所有缺失数据补全；补全后的缺失值分别记为：Rc1、Rc2、…、Rcn；

步骤S5、窗口网格内数据反向计算；具体包括：

根据步骤S4统计的L层缺失值个数n和补全的数据，以及L-1层上的窗口网格，反向计算L-1层中的缺失值；对于第i个缺失值：

C、当对应L层窗口网格中的数值均为缺失值，则L-1层中，每一个缺失值的值为Rci；

D、当对应L层窗口网格中的数值存在部分缺失值，将L层的补全值和L-1层的有效值对比，分下述三种情况：

D1、当L层的补全值均小于L-1层的有效值，则补全值为窗口内的缺失值；

D2、当L层的补全值均大于L-1层的有效值，则补全值就为窗口内的缺失值；

D3、当L层的补全值处于L-1层有效值的最大最小值之间，窗口内的缺失值计算公式如下：

其中，Rc_i^L-1为第L-1层的缺失值，S_k为窗口网格中L-1层上的有效值，共有t个；

重复上述反向计算步骤直至完成第0层雷达反射率的补全。

2.根据权利要求1所述的方法，其特征在于，各层抽吸的窗口网格大小一样。

3.根据权利要求1所述的方法，其特征在于，各层抽吸的窗口网格大小不一样。

4.一种毫米波雷达反射率逐层修复系统：包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至3任一所述方法的步骤。

5.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现上述权利要求1至3任一所述方法中的步骤。