CN109241217B - 一种基于地理网格的灾情数据空间统计、查询及展示方法 - Google Patents
一种基于地理网格的灾情数据空间统计、查询及展示方法 Download PDFInfo
- Publication number
- CN109241217B CN109241217B CN201810924799.1A CN201810924799A CN109241217B CN 109241217 B CN109241217 B CN 109241217B CN 201810924799 A CN201810924799 A CN 201810924799A CN 109241217 B CN109241217 B CN 109241217B
- Authority
- CN
- China
- Prior art keywords
- disaster
- geosot
- area
- data
- grids
- 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.)
- Active
Links
Images
Landscapes
- Alarm Systems (AREA)
Abstract
本发明公开了一种基于地理网格的灾情数据空间统计、查询及展示方法,方案具体为:获取受灾区域的经纬度范围以及灾情数据,灾情数据包含灾情发生的地理位置信息;使用预设统计层级的GeoSOT网格对受灾区域进行覆盖,获得受灾区域的GeoSOT网格集合;计算受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格的GeoSOT编码作为区位编码;建立灾情统计表,灾情统计表中表项对应受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格,将灾情数据分配至灾情发生的地理位置信息所在受灾区域GeoSOT网格对应表项中进行存储,每个表项以对应受灾区域GeoSOT网格的区位编码作为索引号。根据设定的查询条件展示层级,逐级汇总灾情数据,返回并分类浏览相应的灾情数据或其统计结果。
Description
技术领域
本发明涉及灾情数据统计技术领域,具体涉及一种基于地理网格的灾情数据空间统计、查询及展示方法。
背景技术
灾情统计是按照民政部《自然灾害情况统计制度》,对灾区所发生的自然灾害及其破坏损失情况进行调查分析和总结整理,涵盖灾害发生时间、灾害种类、受灾范围、灾害造成的损失以及救灾工作开展情况等众多方面,是制定减灾救灾和重建方案以及社会经济发展规划的重要依据,也是灾害管理的基础资料。因此,为及时、准确、全面、客观地掌握灾情和救灾工作情况,灾情统计是必不可少的工作环节。
目前,我国灾情统计主要是由各级民政部门负责组织、协调和管理。统计人员首先以县(市、区)为基本统计单元,部分省份已延伸至乡镇记录数字、文本、图片等格式的灾情数据,通过电话、传真、电子邮件、国家自然灾害灾情报送系统等形式上报,逐级汇总直至中央。
这种以行政区划为单位进行各类灾情数据统计的方式存在以下不足:
(1)受灾区域的空间区位多以文字形式进行描述,不涉及具体坐标,导致了受灾区域范围界定过程中存在误差,直接影响统计人员上报的灾情数据精度;
(2)主要形成县、市、省、中央等四级行政区划层面上的数据汇总与统计,尺度较少且不连续;
(3)表格、文档、影像等不同格式的灾情数据分开存储,忽略了其内在逻辑联系,无法进行关联检索及展示。
如何实现受灾区域精确定位,多源、多尺度灾情数据统计,是目前灾情数据统计领域亟待解决的问题。
发明内容
有鉴于此,本发明提供了一种基于地理网格的灾情数据空间统计、查询及展示方法,能够实现以地理范围连续的地理网格对灾情数据进行存储、查询和展示,能够将灾情数据对标地理位置,且可以实现多源、多尺度的灾情数据的统计、查询以及展示。
本发明一个实施例提供了一种基于地理网格的灾情数据空间统计方法,该方法包括:
步骤1,获取受灾区域的经纬度范围以及灾情数据,灾情数据包含灾情发生的地理位置信息。
步骤2,使用预设统计层级的GeoSOT网格对受灾区域进行覆盖,获得受灾区域的GeoSOT网格集合。
步骤3,计算受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格的GeoSOT编码作为区位编码。
步骤4,建立灾情统计表,灾情统计表中表项对应受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格,将灾情数据分配至灾情发生的地理位置信息所在受灾区域GeoSOT网格对应表项中进行存储,每个表项以对应受灾区域GeoSOT网格的区位编码作为索引号。
进一步地,获取受灾区域的经纬度范围,具体为:获取受灾区域边界,每隔设定距离设定采样点,采集获取采样点的经纬度信息作为受灾区域的经纬度范围。
进一步地,灾情数据具体包含灾情附件以及灾情发生的地理位置信息;灾情附件用于描述灾情,包含灾情图片、灾情视频以及灾情描述文字;灾情发生的地理位置信息为灾情附件所描述灾情发生的地理位置。
进一步地,预设统计层级的GeoSOT网格为第20层级GeoSOT网格。
本发明另一个实施例还提供了一种基于地理网格的灾情数据查询及展示方法,该方法包括如下步骤:
step1、依据本发明第一个实施例速提供的灾情数据空间统计方法建立灾情统计表。
step2,预先设置查询条件以及展示层级。
查询条件包括待查询区域的经纬度范围,使用设定层级的GeoSOT网格对待查询区域进行覆盖,获得待查询区域的GeoSOT网格集合。
展示层级大于预设统计层级。
step3,计算待查询区域的GeoSOT网格集合中各待查询区域GeoSOT网格的GeoSOT编码作为查询号;
step4,将查询号与灾情统计表中查询号进行匹配查询,获取索引号与查询号一致的表项中所存储的灾情数据作为待查询灾情数据。
step5,将待查询灾情数据按照所属受灾区域GeoSOT网格进行汇总,汇总至展示层级的GeoSOT网格,对展示层级的GeoSOT网格对应待查询灾情数据进行展示。
有益效果:
本发明所提供基于地理网格的灾情数据空间统计、查询及展示方法,是在目前已有的全球剖分网格GeoSOT网格的基础上,以预设统计层级的GeoSOT网格对灾情数据进行区域划分以及区位编码计算,从而能够以相应的区位编码作为索引号实现灾情数据统计表的建立,该灾情数据统计表实现了地理范围连续的灾情数据的统计存储,并且通过索引号实现了更为方便的检索,能够提供多源、多尺度的灾情数据的统计、查询以及展示。
附图说明
图1为本发明实施例所提供的基于地理网格的灾情数据空间统计方法流程图;
图2为本发明实施例所提供的基于地理网格的灾情数据查询及展示方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明所使用的GeoSOT剖分网格,其网格划分方案及编码方案参见公开号为CN102609525的专利申请。基于GeoSOT剖分和编码方案,本发明的基本思想是:使预设统计层级的GeoSOT网格对受灾区域进行无缝无重叠覆盖,得到受灾区域的GeoSOT网格集合;计算受灾区域的GeoSOT网格集合中各网格的编码;将各类灾情数据分配到对应受灾区域的GeoSOT网格,存储时以网格编码标识;根据设定的查询条件展示层级,逐级汇总灾情数据,返回并分类浏览相应的灾情数据或其统计结果。
具体地,本发明其中一个实施例提供了一种基于地理网格的灾情数据空间统计方法,其流程如图1所示,该方法包括:
步骤1,获取受灾区域的经纬度范围以及灾情数据,灾情数据包含灾情发生的地理位置信息。
要获取受灾区域的经纬度范围,可以采用如下具体方式,即:获取受灾区域边界,每隔设定距离设定采样点,该设定距离根据需求,实际实施可以根据统计对象属性和覆盖范围进行设定,例如本发明实施例中可以设置该设定距离为50米,采集获取采样点的经纬度信息作为受灾区域的经纬度范围,本发明实施例中,经纬度信息的采集可以使用GPS等定位设备。
例如,在本发明的一个具体实例中,围绕受灾区域一的边界共进行了12次采样,经纬度坐标依次为:
①(31.459961°N,103.625000°E)
②(31.460449°N,103.624512°E)
③(31.460937°N,103.625000°E)
④(31.461425°N,103.625000°E)
⑤(31.461425°N,103.625488°E)
⑥(31.461425°N,103.625977°E)
⑦(31.461425°N,103.626465°E)
⑧(31.460937°N,103.626465°E)
⑨(31.460449°N,103.626953°E)
⑩(31.459961°N,103.626465°E)
本发明实施例中,为了实现多源、多尺度的灾情数据统计,灾情数据具体包含灾情附件以及灾情发生的地理位置信息;灾情附件用于描述灾情,包含灾情图片、灾情视频以及灾情描述文字;灾情发生的地理位置信息为灾情附件所描述灾情发生的地理位置。
步骤2,使用预设统计层级的GeoSOT网格对受灾区域进行覆盖,获得受灾区域的GeoSOT网格集合。本发明实施例中使用预设统计层级的GeoSOT网格对受灾区域进行覆盖时,应当实现无缝无叠覆盖,即所获得的受灾区域的GeoSOT网格集合应当是无缝隙无重叠的网格集合。
本发明考虑到第20层级GeoSOT网格大小为2″×2″,在中国区域可精确到50米,根据受灾区域面积特点,将第20层级GeoSOT网格设置为最精细地理网格即可满足统计要求。因此本发明实施例中,预设统计层级的GeoSOT网格为第20层级GeoSOT网格,利用第20层级GeoSOT网格对受灾区域进行无缝无叠覆盖。
步骤3,计算受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格的GeoSOT编码作为区位编码。
GeoSOT网格编码方案是基于经纬度坐标空间进行定义的,共有四进制一维、二进制一维、二进制二维和十进制二维四种形式,编码的长度由网格层级决定。
以上述采样点的第一点(31.459961°N,103.625000°E)为例,根据其经纬度,计算其所在第20层级GeoSOT网格的四进制一维编码过程如下:
纬度第1级:31÷256=0余31
经度第1级:103÷256=0余103
则该GeoSOT网格在第1层对应的二进制一维编码为00,将00转化为四进制,进一步可求得该GeoSOT网格在第1层对应的四进制一维编码为0;
纬度第2级:31÷128=0余31
经度第2级:103÷128=0余103
则该GeoSOT网格在第2层对应的二进制一维编码为00,将00转化为四进制,进一步可求得该GeoSOT网格在第2层对应的四进制一维编码为0;
纬度第3级:31÷64=0余31
经度第3级:103÷64=1余39
则该GeoSOT网格在第3层对应的二进制一维编码为01,将01转化为四进制,进一步可求得该GeoSOT网格在第3层对应的四进制一维编码为1;
纬度第4级:31÷32=0余31
经度第4级:39÷32=1余7
则该GeoSOT网格在第4层对应的二进制一维编码为01,将01转化为四进制,进一步可求得该GeoSOT网格在第4层对应的四进制一维编码为1;
纬度第5级:31÷16=1余15
经度第5级:7÷16=0余7
则该GeoSOT网格在第4层对应的二进制一维编码为10,将10转化为四进制,进一步可求得该GeoSOT网格在第4层对应的四进制一维编码为2;
同理可求得该GeoSOT网格在第5-20层对应的四进制一维编码。
将该GeoSOT网格在1-20层的四进制一维编码顺序组合,得到该GeoSOT网格的四进制一维编码为00112233312320202220。
步骤4,建立灾情统计表,灾情统计表中表项对应受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格,将灾情数据分配至灾情发生的地理位置信息所在受灾区域GeoSOT网格对应表项中进行存储,每个表项以对应受灾区域GeoSOT网格的区位编码作为索引号。
例如,若灾情数据可以按照“受灾区域GeoSOT网格区位编码+灾害发生年份+月份+日+灾害种类代码”的形式将每一网格的灾情数据录入数据库,即“001122333123202022202008 05 12 09”作为第一个GeoSOT网格的信息,同理录入剩余17个GeoSOT网格信息。
若受灾区域GeoSOT网格受灾严重,为更清晰的说明受灾情况,对该受灾区域GeoSOT网格所代表的实际区域拍摄图片photo.jpg,则该图片在数据库中按照“受灾区域GeoSOT网格区位编码+图片”的格式存储,即“00112233312320202220photo.jpg”。
其余格式的灾情数据在数据库中的存储方式依此类推,均需以受灾区域GeoSOT网格的区位编码进行标识。
本发明第二个实施例提供了一种基于地理网格的灾情数据查询及展示方法,如图2所示,该方法包括如下步骤:
step1、依据上述灾情数据空间统计方法建立灾情统计表;
step2,预先设置查询条件以及展示层级;
查询条件包括待查询区域的经纬度范围,使用设定层级的GeoSOT网格对待查询区域进行无缝无叠覆盖,获得待查询区域的GeoSOT网格集合;
展示层级大于预设统计层级;
step3,计算待查询区域的GeoSOT网格集合中各待查询区域GeoSOT网格的GeoSOT编码作为查询号。
step4,将查询号与灾情统计表中查询号进行匹配查询,获取索引号与查询号一致的表项中所存储的灾情数据作为待查询灾情数据。
step,将待查询灾情数据按照所属受灾区域GeoSOT网格进行汇总,汇总至展示层级的GeoSOT网格,对展示层级的GeoSOT网格对应待查询灾情数据进行展示。
查询某一区域灾害情况,只需输入该待查询区域的网格区位编码即可搜索出与该编码一致的各类灾情数据,逐一浏览。例如,输入0011223331232020222,可搜索到2条数据:“0011223331232020222 2008 05 12 09”以及“0011223331232020222photo.jpg”。
根据年份、灾害种类等条件查询受灾区域,并根据搜索到的灾害次数对网格进行渲染,颜色深代表灾害多发区,应予以重点关注。例如,第18层级的GeoSOT网格001122333123202023共发生了2次灾害,而18层级上的其余GeoSOT网格只发生了1次灾害。
本发明一个实例中,以第18层级的GeoSOT网格作为展示层级,需进行两次灾情数据的逐级汇总。根据GeoSOT剖分方案,第20层级的网格00112233312320202220与其正上方、右上角、正右方的三个邻近网格共同构成第19级的GeoSOT网格0011223331232020222,将“0011223331232020222 2008 05 12 09”以及“0011223331232020222photo.jpg”2条灾情数据存入数据库,同理合并受灾区域一中的剩余网格,最终受灾区域一由6个第19层级的GeoSOT网格表示。第19层级向第18层级的灾情数据汇总方式与此相同,最终由4个第18层级的GeoSOT网格表示,将相应数据存入数据库。
综上,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于地理网格的灾情数据空间统计方法,其特征在于,该方法包括:
步骤1,获取受灾区域的经纬度范围以及灾情数据,所述灾情数据包含灾情发生的地理位置信息;
步骤2,使用预设统计层级的GeoSOT网格对所述受灾区域进行覆盖,获得所述受灾区域的GeoSOT网格集合;
步骤3,计算所述受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格的GeoSOT编码作为区位编码;
步骤4,建立灾情统计表,所述灾情统计表中表项对应所述受灾区域的GeoSOT网格集合中各受灾区域GeoSOT网格,将所述灾情数据分配至灾情发生的地理位置信息所在受灾区域GeoSOT网格对应表项中进行存储,每个表项以对应受灾区域GeoSOT网格的区位编码作为索引号。
2.如权利要求1所述的方法,其特征在于,所述获取受灾区域的经纬度范围,具体为:
获取受灾区域边界,每隔设定距离设定采样点,采集获取所述采样点的经纬度信息作为所述受灾区域的经纬度范围。
3.如权利要求1所述的方法,其特征在于,所述灾情数据具体包含灾情附件以及灾情发生的地理位置信息;
所述灾情附件用于描述灾情,包含灾情图片、灾情视频以及灾情描述文字;
所述灾情发生的地理位置信息为所述灾情附件所描述灾情发生的地理位置。
4.如权利要求1所述的方法,其特征在于,所述预设统计层级的GeoSOT网格为第20层级GeoSOT网格。
5.一种基于地理网格的灾情数据查询及展示方法,其特征在于,该方法包括如下步骤:
step1、依据如权利要求1所述的灾情数据空间统计方法建立灾情统计表;
step2,预先设置查询条件以及展示层级;
所述查询条件包括待查询区域的经纬度范围,使用所述预设统计层级的GeoSOT网格对所述待查询区域进行覆盖,获得所述待查询区域的GeoSOT网格集合;
所述展示层级大于所述预设统计层级;
step3,计算所述待查询区域的GeoSOT网格集合中各待查询区域GeoSOT网格的GeoSOT编码作为查询号;
step4,将所述查询号与所述灾情统计表中查询号进行匹配查询,获取索引号与查询号一致的表项中所存储的灾情数据作为待查询灾情数据;
step5,将所述待查询灾情数据按照所属受灾区域GeoSOT网格进行汇总,汇总至所述展示层级的GeoSOT网格,对展示层级的GeoSOT网格对应待查询灾情数据进行展示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810924799.1A CN109241217B (zh) | 2018-08-14 | 2018-08-14 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810924799.1A CN109241217B (zh) | 2018-08-14 | 2018-08-14 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109241217A CN109241217A (zh) | 2019-01-18 |
CN109241217B true CN109241217B (zh) | 2021-08-20 |
Family
ID=65071154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810924799.1A Active CN109241217B (zh) | 2018-08-14 | 2018-08-14 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109241217B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110296703A (zh) * | 2019-06-21 | 2019-10-01 | 中国人民解放军陆军工程大学 | 一种应用于大规模无人机群系统中的地理位置编码方法 |
CN110473251B (zh) * | 2019-08-14 | 2022-03-25 | 重庆市地理信息和遥感应用中心(重庆市测绘产品质量检验测试中心) | 基于网格空间索引的自定义范围空间数据面积统计方法 |
CN113656398A (zh) * | 2021-07-06 | 2021-11-16 | 中铁建港航局集团勘察设计院有限公司 | 面向航道信息数据内容的网格动态显示方法、系统及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609525A (zh) * | 2012-02-10 | 2012-07-25 | 北京大学 | 一种统一现有经纬度剖分网格的方法 |
CN103108189A (zh) * | 2013-03-01 | 2013-05-15 | 北京大学 | 一种用于北斗遥感影像数据传输方法 |
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
CN104166695A (zh) * | 2014-07-31 | 2014-11-26 | 北京大学 | 一种面向遥感数据内容的剖分、查询及动态显示方法 |
CN107153711A (zh) * | 2017-05-19 | 2017-09-12 | 北京旋极伏羲大数据技术有限公司 | 地理信息数据处理方法及装置 |
CN107404659A (zh) * | 2016-07-01 | 2017-11-28 | 北京都在哪网讯科技有限公司 | 网络视频显示方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150006569A1 (en) * | 2013-06-28 | 2015-01-01 | James Warner | Method and apparatus for preparing disaster alert notice report |
-
2018
- 2018-08-14 CN CN201810924799.1A patent/CN109241217B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609525A (zh) * | 2012-02-10 | 2012-07-25 | 北京大学 | 一种统一现有经纬度剖分网格的方法 |
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
CN103108189A (zh) * | 2013-03-01 | 2013-05-15 | 北京大学 | 一种用于北斗遥感影像数据传输方法 |
CN104166695A (zh) * | 2014-07-31 | 2014-11-26 | 北京大学 | 一种面向遥感数据内容的剖分、查询及动态显示方法 |
CN107404659A (zh) * | 2016-07-01 | 2017-11-28 | 北京都在哪网讯科技有限公司 | 网络视频显示方法及装置 |
CN107153711A (zh) * | 2017-05-19 | 2017-09-12 | 北京旋极伏羲大数据技术有限公司 | 地理信息数据处理方法及装置 |
Non-Patent Citations (1)
Title |
---|
基于GeoSOT编码的多元灾害数据一体化组织管理方法研究;廖永丰等;《地理与地理信息科学》;20130930;第29卷(第5期);第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109241217A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109241217B (zh) | 一种基于地理网格的灾情数据空间统计、查询及展示方法 | |
Martin | Optimizing census geography: the separation of collection and output geographies | |
CN101452072B (zh) | 一种用于土地监测的电子信息化系统及其方法 | |
US8799074B2 (en) | Methods and apparatus for performing enhanced queries for items such as maps and geospatial data | |
CN112860689B (zh) | 一种基于时空网格数据的平台系统 | |
CN106779417A (zh) | 工程勘察信息数字化采集、管理及综合应用方法 | |
CN107391514A (zh) | 一种数据普查成果输出方法 | |
CN104766168A (zh) | 一种地震三维可视化平台 | |
CN111221867B (zh) | 一种保护性建筑信息管理系统 | |
CN111488634A (zh) | 基于3dgis+bim的铁路勘察设计交付成果管理方法 | |
CN113742505B (zh) | 一种海量合成孔径雷达干涉测量(InSAR)数据在线可视化方法 | |
CN117557681B (zh) | 一种基于多源测绘数据的高精度地形图生成方法及装置 | |
CN113222540A (zh) | 一种生态保护红线的管理系统 | |
CN202533938U (zh) | 基于触摸查询机/移动设备的土地定级估价处理系统 | |
CN111382165B (zh) | 一种移动国土管理系统 | |
Kyung et al. | Spatial data warehouse design and spatial OLAP implementation for decision making of geospatial data update | |
CN113239076A (zh) | 一种基于三维影像的地理信息查询管理平台 | |
CN104964670A (zh) | 一种利用数字设备进行无纸化摄影测量的方法 | |
CN111309842A (zh) | 一种地理测绘专用仪、测绘方法及应用 | |
Saleh | Geographical Information System (GIS)-based analysis of road traffic accident blackspots in Federal Capital Territory (FCT) | |
CN113628060B (zh) | 一种基于WebGIS的红树林信息管理系统 | |
Scholar | Development of Cadastral Information System for Effective and Sustainable Land Management Zakari, Danladi; Surv. Ahmed, Babayo; Adamu Makama Pindiga; & Abubakar Abdulkadir | |
Stroh et al. | Establishing classification and hierarchy in populated place labeling for multiscale mapping for The National Map | |
Griffin | Data quality issues affecting GIS use for environmental problem-solving | |
Hassan et al. | Domination of Technological Tools in Cartographic Advancement: Perspective of Spatial Analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |