一种基于空间插值的电磁频谱分布预测和动态可视化方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106682234A(43)申请公布日2017.05.17(21)申请号201710033617.7(22)申请日2017.01.17(71)申请人北京工业大学地址100124北京市朝阳区平乐园100号(72)发明人牛项朋刘红杰尹良郭健赵光焰(74)专利代理机构北京思海天达知识产权代理有限公司11203代理人张慧(51)Int.Cl.G06F17/30(2006.01)G06Q10/04(2012.01)权利要求书2页说明书4页附图1页(54)发明名称一种基于空间插值的电磁频谱分布预测和动态可视化方法(57)摘要本发明公开一种涉及基于空间插值的电磁频谱分布预测和动态可视化方法，对来自网格化监测的数据按照时间粒度聚合，对每一个时间段内的监测数据使用Kriging空间插值算法，拟合监测站点间距离与电磁信号强度的统计关系，对未知区域的电磁频谱分布数据进行最优无偏插值预测。对每个时间段下的插值预测数据进行热力图渲染并按照时间先后顺序输出到GIS地图中，进而实现对监测区域电磁频谱分布的预测和动态可视化展示。采用本发明的技术方案，可以提高频谱监测数据利用效率。CN106682234ACN106682234A权利要求书1/2页1.一种基于空间插值的电磁频谱分布预测和动态可视化方法，其特征在于，包含以下步骤：第一步，读取网格监测数据，其包括频率、经度、维度、时间、场强；第二步，根据所选频段、门限、时间粒度对监测数据进行筛选压缩，将监测区域进行栅格化处理，栅格交点作为未知数据点；第三步，构建当前时间段ti每个时间段内的已知监测点数据Zij(Lng，Lat，E)；第四步，计算筛选后监测点两两之间的距离{dij}和半方差{γij}；第五步，对距离{dij}和半方差{γij}数据进行曲线拟合，获取拟合曲线函数γij＝f(dij)；第六步，根据拟合曲线重新计算已知监测点两两之间的半方差{γij}；第七步，任选一个未知点，利用拟合函数计算其与所有已知监测点之间的距离{dio}和半方差{γio＝f(dio)}；第八步，计算未知点到各个已知监测点的权重系数λi；第九步，计算未知点的信号强度得到未知点数据Zo(Lng，Lat，E)；第十步，重复第七步到第九步，直到完成所有未知点的数据计算，生成当前时间段ti内的未知点预测数据集{Zo}，并存储；第十一步，判断是否完成所有时间段的插值计算和结果存储，若完成则将所有时间段内的数据集{Zo}进行热力图渲染，并按照时间先后顺序依次展示在GIS地图上；若未完成，则重复第三步到第十步运算过程。2.如权利要求1所述的基于空间插值的电磁频谱分布预测和动态可视化方法，其特征在于，第四步中，监测点之间的距离dij(Zi，Zj)(km)计算涉及球面两点间距计算，经纬度参数需要做一定的转换，转换公式如下：dij[Zi(Lngi，Lati，Ei)，Zj(Lngi，Latj，Ej)]＝-1int((cos(sin(Lati/180×π)×sin(Latj/180×π)+cos(Latj/180×π)×cos(Latj/180×π)×cos(Lngi/180×π-Lngj/180×π))×180×60÷π)×100000)÷100000×1.852各监测点间的半方差γij计算：3.如权利要求1所述的基于空间插值的电磁频谱分布预测和动态可视化方法，其特征在于，第八步中采用Kriging空间插值算法进行曲线拟合，具体为：在Kriging插值中，权重系数λ与半方差{γij}、{γio}存在线性函数关系，具体关系如下：2CN106682234A权利要求书2/2页对矩阵求逆，即可解得：3CN106682234A说明书1/4页一种基于空间插值的电磁频谱分布预测和动态可视化方法技术领域[0001]本发明属于无线电监测技术领域，具体涉及一种基于空间插值的电磁频谱分布预测和动态可视化方法。背景技术[0002]在无线电管理机构的日常无线电监测中，网格化监测作为无线电监测精细化管理的有效手段，正在成为一种越来越普遍的监测方式。然而目前的监测方向集中在频段(频率)占用度的统计分析、单站点频谱监测、不明信号侦听解调和测向定位等方面。然而相对于监测产生的大量频谱监测数据，目前的监测手段对监测数据的利用率仍然显得偏低，对监测区域内的电磁频谱分布和变化情况依然停留在独立站点监测处理的阶段，对于站点以外的大面积区域中的电磁频谱连续分布和变化情况监测没有有效的方法。发明内容[0003]本发明的目的是针对无线电网格化监测中产生大量监测数据，却不能及时反映监测区域内电磁频谱分布和变化情况，数据低效利用的现状，提供一种基于Kriging插值算法的空间电磁态势预测及动态可视化方法。为网格化监测提供一种新的监