专利名称：采用核密度估计的闪电密度分布及ng值估计方法
技术领域：
本发明属于GIS与气象信息技术交叉领域，地理信息系统(GIS)强大之处是能将空间地理分布状况及所具有的属性进行数字存储，建立数据管理系统，同时开发各种分析和处理功能，以便快速获取信息，并将处理结果以地图、图元及数据的形式表示出来。因此， 处理海量且庞杂的实时雷电数据信息，需要通过准确快捷的方法求得区域闪电密度分布及闪电NG值，为雷电防护的分析和决策提供了科学手段，大大提高了以多种形式对雷电数据资料的统计效率。
背景技术：
自然界雷电的发生具有很强的随机性、瞬时性和危险性，它与大气强对流系统密切相关。雷电监测网络主要由地面大气电场监测网、闪电定位监测网、地面自动观测站组成，这个网络包含3-5个电场仪、覆盖全省的9个闪电定位观测点和1-2个自动气象观测站。在当天雷电事件发生后，每个仪器都能即时地监测到的每次雷电的原始数据信息(包括序号、时间、经纬度、强度等信息)，并且以文本信息即时地更新传到服务器中，有时每天的雷电次数少则几条，而梅雨天的次数可能多达数万条。因此，当雷电的原始数据资料以文本记录形式存放时，空间庞大且信息冗杂，而且其灾害分布形势不具有直观性和鲜明性。地理信息系统(GIQ是20世纪60年代中期兴起的一门交叉边缘学科。GIS利用计算机建立地理数据库，将空间地理分布状况及所具有的属性进行数字存储，建立数据管理系统，同时开发各种分析和处理功能，以便快速获取信息，并将处理结果以地图、图元及数据的形式表示出来。因此，结合雷电监测网络的实时雷电数据信息，并通过GIS的强大的空间查询分析和图形可视化功能为雷电防护的分析和决策提供了科学手段，一方面能够利用GIS的地图及图层加载功能，实现雷电精确定位及资料回放显示效果；另一方面能够利用GIS的空间查询分析及专题图制作功能，实现雷电资料查询显示以及闪电密度分析功能，大大提高了以多种形式对雷电数据资料的统计效率。参考文献[1]陈家宏，樊灵孟.雷电定位系统和雷电参数[C]. The 3rd Asian Light-ning Protection Forum. Taisei Corp, Tokyo, Japan,2005.[2]冯民学，焦雪等.江苏省雷电分布特征分析[J].气象科学，2009，^K2) 246-250[3]刘岩，王振会，张慧良等.2007年浙江地区梅雨期闪电特征的分析[J].气象科学，2009，29 U) :225-228.[4]刘艳，马劲松.核密度估计法在西藏人口空间分布研究中的应用，西藏科技， 2007(4)，1498-1504[5]杨宜东，孙志辉，张净.基于核密度估计的分布数据流离群点检测.计算机研究与发展，2005，42 (9)[6]国网武汉高压研究院.福建省雷电参数研究[R].武汉国网武汉高压研究院，2006.[7]Krider E,Pifer A E,Vance D L. Lightuing direction finding system for forest fire detection.Ball Amer Metor SoC，1980,61:980-986.[8]B. W Silverman, Density Estimation for Statistics and Data Analysis, Momographs on Statistics and Applied Probability,1986.
发明内容
本发明目的是针对现有技术存在的缺陷提供一种，提供采用核密度估计的闪电密度分布及NG值统计方法，该方法有利于在雷电业务平台中进行闪电分布密度图及灾害指数NG的计算。本发明为实现上述目的，采用如下技术方案本发明采用核密度估计的闪电密度分布及NG值估计方法，包括以下步骤(1)读取闪电资料数据库中的时间、位置信息、属性信息，形成点空间元素集合；(2)确定平滑带宽，计算闪电密度计算尺度及密度函数，求出闪电密度集体；(3)将闪电点矢量图层按照密度值进行渲染，按照条件区域进行适量图层切割；(4)将矢量图层栅格化，并且按照上述条件区域进行栅格图层切割；(5)计算用户选择的区域面积及区域闪电次数，求出闪电NG值，最终密度图与数字综合展示。优选地，步骤(1)的实现方法是输入时间跨度条件后，读取闪电资料数据库中相应闪电的信息，形成对应的二维空间链表；再利用ArcGIS Engine 9. 2组件读取二维空间链表，将数据信息转成GIS地图元素，最终形成闪电点图层。优选地，步骤O)的实现方法是闪电集合S是闪电Si空间R上的离散点集，S = Is1, S2, ... , sn}，给定带宽bandwidth或平滑量τ，计算点密度式中，带宽参数τ >0，是一个以s为中心的圆的半径；δ τ (s)代表了密度计算的尺度，k是一个概率密度函数即核心kernel密度函数；并将点密度添加到相应闪电元素的字段 density 中，i = 1,2,…η。优选地，步骤C3)的实现方法是根据闪电元素字段density，选择色带、区间值参数，利用ArcGIS Engine 9. 2组件进行渲染，形成矢量面图层。优选地，步骤的实现方法是利用ArcGIS Engine 9. 2组件，设定分析环境参数，包括地理分析环境参数，通过执行IGeoAnalysisEnvironment接口实现；栅格分析环境参数，通过执行IRasterAnalysisEnvironment接口实现对矢量图层的栅格化。优选地，步骤(5)的实现方法是将所有闪电与条件区域作交集运算，求得该区域的闪电总次数，并计算条件区域的总面积，即可求出闪电NG值NG = M/Y/P其中，M为闪电总次数，Y为计算年限，P为区域面积。本发明可以对复杂且庞大的雷电灾害数据进行时域和空间分布特征分析，利用核密度估计法求得区域内闪电密度分布以及年平均地面落雷密度(NG)值，并借助强大的GIS 空间数据表现形式，了解区域内闪电的分布特征对确定闪电密度有重要意义，以便气象预报人员能够出色地完成闪电观测及分析工作。通常雷电地基监测网络由地面大气电场监测网、闪电定位监测网、地面自动观测站组成，但是粗糙的原始数据依然无法及时、便捷地提供给雷电防护人员观测分析。
图1:本发明原理图。
具体实施例方式
权利要求
1.一种采用核密度估计的闪电密度分布及NG值估计方法，其特性在于，包括以下步骤(1)读取闪电资料数据库中的时间、位置信息、属性信息，形成点空间元素集合；(2)确定平滑带宽，计算闪电密度计算尺度及密度函数，求出闪电密度集体；(3)将闪电点矢量图层按照密度值进行渲染，按照条件区域进行适量图层切割；(4)将矢量图层栅格化，并且按照上述条件区域进行栅格图层切割；(5)计算用户选择的区域面积及区域闪电次数，求出闪电NG值，最终密度图与数字综合展示。
2.根据权利要求1所述的采用核密度估计的闪电密度分布及NG值估计方法，其特征在于步骤(1)的实现方法是输入时间跨度条件后，读取闪电资料数据库中相应闪电的信息，形成对应的二维空间链表；再利用ArcGIS Engine 9. 2组件读取二维空间链表，将数据信息转成GIS地图元素，最终形成闪电点图层。
3.根据权利要求1所述的采用核密度估计的闪电密度分布及NG值估计方法，其特征在于步骤⑵的实现方法是闪电集合S是闪电Si空间R上的离散点集，S = Is1, s2,..., sn}，给定带宽bandwidth或平滑量τ，计算点密度1^As) τ VtJ式中，带宽参数τ >0，是一个以s为中心的圆的半径；δ τ(8)代表了密度计算的尺度，k是一个概率密度函数即核心kernel密度函数；并将点密度添加到相应闪电元素的字段 density 中，i = 1，2，…η。
4.根据权利要求1所述的采用核密度估计的闪电密度分布及NG值估计方法，其特征在于步骤(3)的实现方法是根据闪电元素字段density，选择色带、区间值参数，利用 ArcGIS Engine 9. 2组件进行渲染，形成矢量面图层。
5.根据权利要求1所述的采用核密度估计的闪电密度分布及NG值估计方法，其特征在于步骤(4)的实现方法是利用ArcGIS Engine 9. 2组件，设定分析环境参数，包括地理分析环境参数，通过执行IGeoAnalysisEnvironment接口实现；栅格分析环境参数，通过执行 ^asterAnalysisEnvironment接口实现对矢量图层的栅格化。
6.根据权利要求1所述的采用核密度估计的闪电密度分布及NG值估计方法，其特征在于步骤(5)的实现方法是将所有闪电与条件区域作交集运算，求得该区域的闪电总次数，并计算条件区域的总面积，即可求出闪电NG值NG = M/Y/P其中，M为闪电总次数，Y为计算年限，P为区域面积。
全文摘要
本发明公布了一种采用核密度估计的闪电密度分布及NG值估计方法。本发明本方法实现①结合地理信息系统(GIS)技术，构建和维护雷电空间数据库，并建立雷电实时定位监测系统来获取雷电发生的信息源，以空间要素形式实时地呈现闪电的定位分布；②利用核密度估计方法，求得区域闪电密度分布，再将用户选择的图形区域作为条件层作切割，形成区域密度分布图，并且能够消除边缘效应。③计算核密度估计的区域总面积(公里)以及年总闪电次数，求得对应区域的NG值，形成闪电密度图和数值结合的气象产品，最后通过共享网络来发布和共享雷电资料分析产品，通过以上闪电密度统计方法实现雷电业务服务的根本目标。
文档编号G06F19/00GK102254090SQ20111015303
公开日2011年11月23日 申请日期2011年6月8日 优先权日2011年6月8日
发明者侯荣涛, 史鑫明, 朱斌 申请人:南京信息工程大学