CN102565554A - 一种生成三维雷电定位地闪点分布图的方法 - Google Patents
一种生成三维雷电定位地闪点分布图的方法 Download PDFInfo
- Publication number
- CN102565554A CN102565554A CN201210026360XA CN201210026360A CN102565554A CN 102565554 A CN102565554 A CN 102565554A CN 201210026360X A CN201210026360X A CN 201210026360XA CN 201210026360 A CN201210026360 A CN 201210026360A CN 102565554 A CN102565554 A CN 102565554A
- Authority
- CN
- China
- Prior art keywords
- lightning
- geographical
- basic unit
- point
- thunder
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种生成三维雷电定位地闪点分布图的方法,该方法首先利用现有的二维雷电定位地闪点分布图和海拔高程地图,通过栅格转换,生成新的三维雷电定位地闪点分布图的地理基层和雷电地闪点基层,然后统一地理基层和雷电地闪点基层的坐标系,最后对地理基层和雷电地闪点基层实施投影叠置分析,把地理基层对应位置的海拔高程数据传递给雷击地闪点,生成三维雷电定位地闪点分布图。该方法改进了现有传统二维雷电定位地闪点分布图的不足,使雷击地闪点获取了三维海拔属性,为进一步分析雷电活动规律和改进雷电参数算法打下了基础。
Description
技术领域
本发明属于防雷工程技术领域,具体地说涉及一种获得三维雷电定位地闪点分布图的方法。
背景技术
雷电定位技术是目前电力系统研究雷云对地放电,确定雷击地闪点分布的主流方法,该定位技术是采用雷电电磁场遥测技术实现对雷击故障点的定位和雷电流计算,记录并实时显示探测区域发生的地闪雷电的时间、位置、强度等参数,生成雷电定位地闪点分布图,对电力系统输电线路防雷保护设计及雷击线路跳闸后的现场巡查工作都很有帮助,但它也存在以下不足:现有的雷电定位系统是基于二维地理信息系统,反映的是雷击地闪点在平面维度上的信息,得到的雷电定位地闪点分布图是二维的,随之带来以下问题:(1)无法获得雷击地闪点与三维海拔高程之间的关系,从而无法研究雷电流参数随海拔变化的规律;(2)在计算地面落雷密度参数信息时,由于所用的是二维地图,计算用公式中的面积是实际地表面积在水平方向上的投影,导致计算出来的结果偏大,无法真实反映地面落雷密度参数。
综上所述,传统的雷电定位地闪点分布图不能满足精确计算雷电流参数和研究海拔因素对雷电规律影响的需求。在雷电定位领域,亟需一种可以统计三维地理信息下的雷电活动规律的方法。
发明内容
本发明的目的是针对传统雷电定位地闪点分布图只能实现雷击地闪点在二维地图的显示,无法反映雷击地闪点分布与三维海拔之间的关系,在山区等复杂地理环境下对输电线路防雷措施的指导意义受限的现状,通过对现有雷电定位系统获取的雷电地闪定位数据和反映三维地图的海拔高程数据的叠置加工,提出了一种生成三维雷电定位地闪点分布图的方法,以弥补传统二维雷电定位地闪点分布图的不足。
本发明的目的通过如下技术方案实现。
一种生成三维雷电定位地闪点分布图的方法,包括下列步骤:
步骤1:将现有的二维雷电定位地闪点分布图和海拔高程地图,通过栅格转换,生成新的三维雷电定位地闪点分布图的地理基层和雷电地闪点基层;
步骤2:统一地理基层和雷电地闪点基层的坐标系;
步骤3:对地理基层和雷电地闪点基层实施投影叠置分析,把地理基层对应位置的海拔高程数据传递给雷击地闪点,生成三维雷电定位地闪点分布图;
其中:
步骤1中把地图中地理实体模拟图形转换成计算机能够处理的栅格数据的方法为,将地理图形空间划分成规则的小单元即像元,地理空间位置由所述像元的行列号表示,地图数据的分辨率由像元的大小表示,地理空间实体由若干像元隐含描述;点地理实体用一个栅格单元表示;线状地理实体则用沿线走向的一组相邻栅格单元表示,每个栅格单元最多只有两个相邻单元在线上;面或区域地理实体用记有区域属性的相邻栅格单元的集合表示,每个栅格单元可有多于两个的相邻单元同属一个区域;地理地图数据直接记录属性的指针或属性本身,而所在位置则根据行列号转换为相应的坐标给出;每个栅格存贮单元的行列位置可以方便地根据其在数据文件中的记录位置得到,且行列坐标可以转为其他坐标系下的坐标;
步骤2中把地理基层和雷电地闪点基层的坐标系统一的方法为:
式中,X、Y、Z是中国西安80或北京54坐标系的点坐标;Xw、Yw、Zw是WGS-84坐标系下的点坐标,ζx、ζy、ζz是坐标系之间的旋转参数;X0、Y0、Z0是坐标系之间的平移参数;K是坐标系之间的尺度参数;其中七个转换参数可以通过地图区域的三个公共点求得;
步骤3中把地理基层和雷电地闪点基层实施投影的方法为:采用正解变换的方法,即通过建立地理基层投影变换为雷电地闪点基层投影的解析关系式,由雷电地闪点基层投影的数字化坐标x、y变换到地理基层投影的直角坐标X、Y;且步骤3中把地理基层和雷电地闪点基层实施叠置分析的方法为:通过对投影地理区域多重属性的模拟,寻找和确定同时具有几种地理属性的位置,按照确定的地理坐标,对其进行逻辑交运算,即:
式中,x是雷击地闪点的地理坐标;A,B分别是地理基层和雷电地闪点基层,通过逻辑交运算,使雷击地闪点得到地理基层的海拔高程属性。
本发明方法获得的三维雷电定位地闪图能包含传统雷电定位系统覆盖区域的雷电定位地闪点的三维地理坐标信息,具有以下特点和优点:
1、本发明提出的方法能够获取雷击地闪点的三维地理坐标信息,比传统雷电定位图多增加了地闪点的海拔信息;
2、本发明提出的方法能够运用于研究雷击地闪点随三维地理变化的规律;
3、本发明提出的方法有益于改善现有计算雷击地闪密度的算法,进一步精确雷击地闪密度的计算值。
4、本发明提出的方法算法简单可靠,可以快速准确地实现三维雷电地闪点分布图。
附图说明
图1是本发明方法的总体流程图;
图2是本发明实施例地理图形转换成栅格数据示意图;其中图2(a)是图形模拟表示的地理对象,其中图2(b)是地理对象对应的栅格数据模型;
图3是本发明实施例投影叠置分析示意图;
图4是本发明实施例三维雷电定位地闪点分布示意图(图中黑点代表雷击地闪点)。
具体实施方式
下面通过云南某地的三维雷电定位地闪点分布图实施例,结合附图,对本发明的技术方案做进一步说明。
如图1所示,一种三维雷电定位分布图的生成方法,包括下列步骤:
步骤1:将现有的二维雷电定位地闪点分布图和海拔高程地图,通过栅格转换,生成新的三维雷电定位地闪点分布图的地理基层和雷电地闪点基层;
步骤2:统一地理基层和雷电地闪点基层的坐标系;
步骤3:对地理基层和雷电地闪点基层实施投影叠置分析,把地理基层对应位置的海拔高程数据传递给雷击地闪点,生成三维雷电定位地闪点分布图;
其中:
步骤1中把地图中地理实体模拟图形转换成计算机能够处理的栅格数据的方法为,将地理图形空间划分成规则的小单元即像元,地理空间位置由所述像元的行列号表示,地图数据的分辨率由像元的大小表示,地理空间实体由若干像元隐含描述;点地理实体用一个栅格单元表示;线状地理实体则用沿线走向的一组相邻栅格单元表示,每个栅格单元最多只有两个相邻单元在线上;面或区域地理实体用记有区域属性的相邻栅格单元的集合表示,每个栅格单元可有多于两个的相邻单元同属一个区域;地理地图数据直接记录属性的指针或属性本身,而所在位置则根据行列号转换为相应的坐标给出;每个栅格存贮单元的行列位置可以方便地根据其在数据文件中的记录位置得到,且行列坐标可以转为其他坐标系下的坐标;
步骤2中把地理基层和雷电地闪点基层的坐标系统一的方法为:
式中,X、Y、Z是中国西安80或北京54坐标系的点坐标;Xw、Yw、Zw是WGS-84坐标系下的点坐标,ζx、ζy、ζz是坐标系之间的旋转参数;X0、Y0、Z0是坐标系之间的平移参数;K是坐标系之间的尺度参数;其中七个转换参数可以通过地图区域的三个公共点求得;
步骤3中把地理基层和雷电地闪点基层实施投影的方法为:采用正解变换的方法,即通过建立地理基层投影变换为雷电地闪点基层投影的解析关系式,由雷电地闪点基层投影的数字化坐标x、y变换到地理基层投影的直角坐标X、Y;且步骤3中把地理基层和雷电地闪点基层实施叠置分析的方法为:通过对投影地理区域多重属性的模拟,寻找和确定同时具有几种地理属性的位置,按照确定的地理坐标,对其进行逻辑交运算,即:
式中,x是雷击地闪点的地理坐标;A,B分别是地理基层和雷电地闪点基层。通过逻辑交运算,使雷击地闪点得到地理基层的海拔高程属性。
第一步:在目前正在使用的云南雷电定位系统数据库中提取云南某地的二维雷电定位地闪点分布图(20110701.shp),在中国科学技术数据库中提取云南某地的三维海拔高程地图(ASTGTM_N25E103W_DEM_UTM.img),利用本发明实现方法中的步骤1,将上述两地图通过栅格变换的方法转化成计算机可以处理的栅格数据。如图2(a)中表示了一个代码为1的地理点实体,一条代码为2的地理线实体,一个代码为3的地理面实体,图2(b)中则表示代码分别为1、2、3的点、线、面地理实体对应的转换后的栅格数据模型。经过转换,可以生成新的云南某地的三维雷电定位地闪点分布图的地理基层和雷电地闪点基层。
第二步:使用本发明实现方法中的步骤2把云南某地的三维海拔高程地图的WGS84坐标系转换成云南某地的二维雷电定位地闪点分布图中国西安80坐标系,保证两幅图的坐标起点,旋转参数和坐标尺度保持一致,以保证下一步叠置分析的精确性。
第三步:把云南某地的二维雷电定位地闪点分布基层投影叠置到对应的三维海拔高程地图之上如图3所示,进行叠置分析,通过逻辑交计算,把雷击地闪点对应位置的海拔数据属性从云南某地的三维地理基层传递给点P,依次类推,用批处理的方法使云南某地的雷击地闪点Pi(i=1,2...n,n代表雷击地闪点总数)全部获取海拔数据属性,从而可以得到此云南某地的三维雷电定位地闪点分布图。如图4所示。
以上对本发明实施方式进行了详细的介绍,本文中应用了具体实施例对本发明进行了阐述,以上实施例的说明只是用于帮助理解本发明的系统及方法。同时,对本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (1)
1. 一种生成三维雷电定位地闪点分布图的方法,其特征在于,所述方法包括下列步骤:
步骤1:将现有的二维雷电定位地闪点分布图和海拔高程地图,通过栅格转换,生成新的三维雷电定位地闪点分布图的地理基层和雷电地闪点基层;
步骤2:统一地理基层和雷电地闪点基层的坐标系;
步骤3:对地理基层和雷电地闪点基层实施投影叠置分析,把地理基层对应位置的海拔高程数据传递给雷击地闪点,生成三维雷电定位地闪点分布图;
其中:
步骤1中把地图中地理实体模拟图形转换成计算机能够处理的栅格数据的方法为,将地理图形空间划分成规则的小单元即像元,地理空间位置由所述像元的行列号表示,地图数据的分辨率由像元的大小表示,地理空间实体由若干像元隐含描述;点地理实体用一个栅格单元表示;线状地理实体则用沿线走向的一组相邻栅格单元表示,每个栅格单元最多只有两个相邻单元在线上;面或区域地理实体用记有区域属性的相邻栅格单元的集合表示,每个栅格单元可有多于两个的相邻单元同属一个区域;地理地图数据直接记录属性的指针或属性本身,而所在位置则根据行列号转换为相应的坐标给出;每个栅格存贮单元的行列位置可以方便地根据其在数据文件中的记录位置得到,且行列坐标可以转为其他坐标系下的坐标;
步骤2中把地理基层和雷电地闪点基层的坐标系统一的方法为:
式中,X、Y、Z是中国西安80或北京54坐标系的点坐标;Xw、Yw、Zw是WGS-84坐标系下的点坐标,ζx、ζy、ζz是坐标系之间的旋转参数;X0、Y0、Z0是坐标系之间的平移参数;K是坐标系之间的尺度参数;其中七个转换参数可以通过地图区域的三个公共点求得;
步骤3中把地理基层和雷电地闪点基层实施投影的方法为:采用正解变换的方法,即通过建立地理基层投影变换为雷电地闪点基层投影的解析关系式,由雷电地闪点基层投影的数字化坐标x、y变换到地理基层投影的直角坐标X、Y;且步骤3中把地理基层和雷电地闪点基层实施叠置分析的方法为:通过对投影地理区域多重属性的模拟,寻找和确定同时具有几种地理属性的位置,按照确定的地理坐标,对其进行逻辑交运算,即:
式中,x是雷击地闪点的地理坐标;A,B分别是地理基层和雷电地闪点基层,通过逻辑交运算,使雷击地闪点得到地理基层的海拔高程属性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210026360.XA CN102565554B (zh) | 2012-02-07 | 2012-02-07 | 一种生成三维雷电定位地闪点分布图的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210026360.XA CN102565554B (zh) | 2012-02-07 | 2012-02-07 | 一种生成三维雷电定位地闪点分布图的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102565554A true CN102565554A (zh) | 2012-07-11 |
CN102565554B CN102565554B (zh) | 2014-12-31 |
Family
ID=46411482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210026360.XA Active CN102565554B (zh) | 2012-02-07 | 2012-02-07 | 一种生成三维雷电定位地闪点分布图的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102565554B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258140A (zh) * | 2013-05-29 | 2013-08-21 | 广东电网公司江门供电局 | 用于雷暴过程还原、雷电通道分析的方法及其系统 |
CN103605790A (zh) * | 2013-12-04 | 2014-02-26 | 武汉聪源科技有限责任公司 | 一种用于引雷设施选址分析的方法及其系统 |
CN107765102A (zh) * | 2017-09-13 | 2018-03-06 | 广东电网有限责任公司江门供电局 | 一种雷电参数反演计算优化方法 |
CN110221180A (zh) * | 2019-07-01 | 2019-09-10 | 国网四川省电力公司电力科学研究院 | 一种10kV配电线路雷击故障识别与定位方法 |
CN110426562A (zh) * | 2019-06-24 | 2019-11-08 | 乐山师范学院 | 基于分层搜索和距离空间投影的高精度闪电三维定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101621191A (zh) * | 2009-08-06 | 2010-01-06 | 重庆市电力公司超高压局 | 超高压输电线路三维可视化处理方法及系统 |
CN101650444A (zh) * | 2008-08-11 | 2010-02-17 | 华东电力试验研究院有限公司 | 落雷密度统计方法 |
CN101872492A (zh) * | 2010-06-09 | 2010-10-27 | 中国科学院深圳先进技术研究院 | 三维仿真城市的多角度地图实现方法 |
CN101908062A (zh) * | 2010-07-02 | 2010-12-08 | 中国科学院计算技术研究所 | Gis空间谓词判断方法及其系统 |
CN102096713A (zh) * | 2011-01-29 | 2011-06-15 | 广州都市圈网络科技有限公司 | 一种基于网格化二三维地图匹配方法及系统 |
-
2012
- 2012-02-07 CN CN201210026360.XA patent/CN102565554B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650444A (zh) * | 2008-08-11 | 2010-02-17 | 华东电力试验研究院有限公司 | 落雷密度统计方法 |
CN101621191A (zh) * | 2009-08-06 | 2010-01-06 | 重庆市电力公司超高压局 | 超高压输电线路三维可视化处理方法及系统 |
CN101872492A (zh) * | 2010-06-09 | 2010-10-27 | 中国科学院深圳先进技术研究院 | 三维仿真城市的多角度地图实现方法 |
CN101908062A (zh) * | 2010-07-02 | 2010-12-08 | 中国科学院计算技术研究所 | Gis空间谓词判断方法及其系统 |
CN102096713A (zh) * | 2011-01-29 | 2011-06-15 | 广州都市圈网络科技有限公司 | 一种基于网格化二三维地图匹配方法及系统 |
Non-Patent Citations (3)
Title |
---|
李永福等: "基于雷电定位数据的雷电流参数岁海拔变化规律", 《高电压技术》 * |
童雪芳等: "雷电定位系统地闪密度分布图与雷击故障相关性分析", 《高电压技术》 * |
迟文学等: "雷电监测数据三维可视化技术研究", 《计算机工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258140A (zh) * | 2013-05-29 | 2013-08-21 | 广东电网公司江门供电局 | 用于雷暴过程还原、雷电通道分析的方法及其系统 |
CN103258140B (zh) * | 2013-05-29 | 2015-04-01 | 广东电网公司江门供电局 | 用于雷暴过程还原、雷电通道分析的方法及其系统 |
CN103605790A (zh) * | 2013-12-04 | 2014-02-26 | 武汉聪源科技有限责任公司 | 一种用于引雷设施选址分析的方法及其系统 |
CN103605790B (zh) * | 2013-12-04 | 2017-02-08 | 武汉聪源科技有限责任公司 | 一种用于引雷设施选址分析的方法及其系统 |
CN107765102A (zh) * | 2017-09-13 | 2018-03-06 | 广东电网有限责任公司江门供电局 | 一种雷电参数反演计算优化方法 |
CN107765102B (zh) * | 2017-09-13 | 2019-11-22 | 广东电网有限责任公司江门供电局 | 一种雷电参数反演计算优化方法 |
CN110426562A (zh) * | 2019-06-24 | 2019-11-08 | 乐山师范学院 | 基于分层搜索和距离空间投影的高精度闪电三维定位方法 |
CN110426562B (zh) * | 2019-06-24 | 2021-06-04 | 乐山师范学院 | 基于分层搜索和距离空间投影的高精度闪电三维定位方法 |
CN110221180A (zh) * | 2019-07-01 | 2019-09-10 | 国网四川省电力公司电力科学研究院 | 一种10kV配电线路雷击故障识别与定位方法 |
CN110221180B (zh) * | 2019-07-01 | 2021-07-06 | 国网四川省电力公司电力科学研究院 | 一种10kV配电线路雷击故障识别与定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102565554B (zh) | 2014-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110675494B (zh) | 一种风电场风机位置三维可视化优选方法 | |
Gao et al. | Three-dimensional propagation characteristics of the upward connecting leaders in six negative tall-object flashes in Guangzhou | |
CN103559360B (zh) | 一种三维架空电力线模型的构建方法及系统 | |
CN103196426A (zh) | 一种全站仪联合三维激光扫描仪用于建筑测量的方法 | |
CN102662179A (zh) | 基于机载激光雷达的三维优化选线方法 | |
CN102565554B (zh) | 一种生成三维雷电定位地闪点分布图的方法 | |
CN104835202A (zh) | 一种三维虚拟场景快速构建方法 | |
JP5020106B2 (ja) | 風力発電量予測システム、方法及びプログラム | |
CN102706323B (zh) | 基于机载激光雷达数据的塔基断面提取方法 | |
CN102750359A (zh) | 基于输电线路三维模型的空间信息提取与格式转换方法 | |
CN105354882A (zh) | 一种大数据架构的大空间范围输电三维全景展示平台构建方法 | |
CN105913140A (zh) | 一种基于建筑信息模型的发电预测方法和装置 | |
CN104217010B (zh) | 一种基于分段地形的导线跳闸率计算及评估方法 | |
CN114170524A (zh) | 一种高时-空卫星数据融合的单日潮滩数字地形构建方法 | |
CN111540046B (zh) | 一种基于虚拟现实的输电通道仿真方法 | |
Zhou et al. | Smart grid data mining and visualization | |
CN109472862A (zh) | 一种变电站三维建模系统 | |
CN106844968A (zh) | 基于环境敏感目标自动识别的输变电工程电磁环境预测方法 | |
CN112541620B (zh) | 一种预测精度及效率高的台风风暴增水预测方法及系统 | |
CN112862921B (zh) | 一种电网分布图像绘制方法 | |
Su et al. | Automatic multi-source data fusion technique of powerline corridor using UAV Lidar | |
CN103235290B (zh) | 一种基于地理空间点阵的雷达探测数据处理方法 | |
Schuffert | An automatic data driven approach to derive photovoltaic-suitable roof surfaces from ALS data | |
CN102840850B (zh) | 一种输电线路实际防雷保护角的测定方法 | |
Hu et al. | A 3D Modeling Method for Substations Based on Oblique Photography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220624 Address after: 650217 yundian science and Technology Park, No.105, Yunda West Road, economic development zone, Kunming City, Yunnan Province Patentee after: YUNNAN ELECTRIC POWER TEST & RESEARCH INSTITUTE (GROUP) Co.,Ltd. Address before: 650217 Yunnan Kunming Economic and Technological Development Zone Yunda West Road Middle Yundian Science Park Patentee before: YUN NAN ELECTRIC TEST & RESEARCH INSTITUTE GROUP Co.,Ltd. ELECTRIC INSTITUTE |