一种基于矢量引导的洪灾灾情评估方法与流程

1.本发明涉及灾情评估技术领域，尤其涉及一种基于矢量引导的洪灾灾情评估方法。


背景技术：

2.洪水灾害是我国面临的重大突发性自然灾害之一，将现代信息技术尤其是以“3s”为核心的空间信息技术运用于洪水灾害的动态监测、灾情评估、救灾决策支持等领域具有重大的意义，如何利用现代空间信息技术为洪水灾害监测、救援、评估等提供更科学、快速的辅助决策支持一直是急需解决的重大问题，加强对洪水灾害损失评估的研究，不仅能够及时的了解洪涝灾害发生的几率，而且也对洪涝灾害引起的损失加以预测，为防洪减灾部门救灾、抗灾计划的制定提供数据支持；
3.至今为止，还没有一套完善的传感器或系统能够从遥感影像稳定地、全自动地提取目标地物并做出评估，现有的一些洪灾灾情评估方法大都存在实施较为耗时，时效性不高的问题，且评估过程中所需数据众多，分析较为复杂，分析难度较高，另外现有的方法在评估精度上也有所不足，为减灾救灾服务提供更加快速和准确的决策支持，因此，本发明提出一种基于矢量引导的洪灾灾情评估方法以解决现有技术中存在的问题。


技术实现要素：

4.针对上述问题，本发明的目的在于提出一种基于矢量引导的洪灾灾情评估方法，该方法通过洪水期遥感影像得到水体范围，然后与地理矢量数据进行匹配可以快速定位至洪灾发生地，即可进行灾损评估得到灾损信息，在克服处理资源受限的条件下，实现灾害信息的快速提取，实施简单快捷，所需成本低，利用地理矢量信息的辅助功能，充分深挖不同数据内部和数据之间的协调能力，实现洪涝灾害评估。
5.为了实现本发明的目的，本发明通过以下技术方案实现：一种基于矢量引导的洪灾灾情评估方法，包括以下步骤：
6.步骤一：水体提取
7.基于待评估区域的遥感卫星遥感影像，利用归一化差异水体指数ndwi对待评估区域进行水体提取，其中ndwi定义如下：
[0008][0009]
步骤二：矢量引导
[0010]
根据步骤一，先将基于ndvi提取的水体转为矢量然后导出成图片，然后将待评估区域的水体矢量数据导出成图片，其中待评估区域的水体矢量数据图片作为被匹配模板；
[0011]
步骤三：确定洪水淹没范围
[0012]
根据步骤二，先对基于ndvi提取的水体矢量图片进行模板匹配，确定具体的灾害发生地，再将提取的水体与灾害发生地灾害发生前的卫星遥感影像进行叠加，确定洪水淹
没区域信息；
[0013]
步骤四：估算洪灾损失
[0014]
根据步骤三，根据洪水淹没区域信息对洪灾发生地的洪灾损失进行计算评估，得到灾损信息，公式如下：
[0015][0016]
其中，s为灾害发生地的直接经济损失，p
ab
是a种淹没范围内第b类财产值，t
ab
是第a种淹没级别下的第b类财产损失率，m为财产类型的总数，n是洪水淹没级别。
[0017]
进一步改进在于：所述步骤一中，所述水体的提取方法还包括基于阈值的多波段谱间关系法和阴影水体指数法，其中基于阈值的多波段谱间关系法主要利用影像多波段的综合优势进行水体信息的提取，包括谱间关系法与比值法，阴影水体指数法是在现有的多波段谱间关系法基础上，根据水体的光谱特性进行构建。
[0018]
进一步改进在于：所述步骤三中，所述洪水淹没区域信息包括洪水淹没范围和淹没区域土地利用类型，所述灾害发生地灾害发生前的影像为从全国地理信息资源目录服务上下载的待评估区域的土地利用数据影像。
[0019]
进一步改进在于：所述步骤三中，所述模板匹配具体为在一副大图像中搜寻目标，已知该图中存在要找的目标，且该目标同模板有相同的尺寸、方向和图像元素，通过ssda算法在图中找到目标，确定其坐标位置，模板图像记为t，待搜索图像记为f，计算图像f(x,y)在像素(u,v)的非相似度m(u,v)为：
[0020][0021]
以m(u,v)作为匹配尺度，式中(u,v)表示模板与图像重合部分的左上角像素坐标，模板的大小为：n
×
m。
[0022]
进一步改进在于：所述步骤四中，所述洪灾损失包括农作物损失、畜牧业损失和居民财产损失，所述居民财产损失包括房屋建筑损失和其他固定资产损失。
[0023]
进一步改进在于：所述农作物损失的估算公式如下：
[0024][0025]
其中，r
农
为洪灾农作物经济损失，η
ij
为第j级淹没水深下的第i类农作物洪灾损失率，w
ij
为第j级淹没水深下的第i类农作物正常年产值，n为农作物种类数，m为淹没水深等级数，r
存
为被洪水冲走的农作物或因受淹霉变造成的损失。
[0026]
进一步改进在于：所述畜牧业损失的估算公式如下：
[0027][0028]
其中，na为牲畜死亡的种类数，gt为死亡牲畜的活体重量，vt为死亡牲畜的批发价格。
[0029]
进一步改进在于：所述房屋建筑损失的估算公式如下：
[0030]
l建筑＝a淹
×c×v×r[0031]
其中，l建筑为房屋损失价值，a淹为对应的农村居民点或城镇淹没面积，c为农村居民点或城镇土地利用面积系数，v为单位面积上的房屋价值，r为损失率。
[0032]
进一步改进在于：所述其他固定资产损失的计算公式如下：
[0033]
l固资＝(a淹/h)
×c×
v固
×r[0034]
其中，l固资为其他固定资产损失，a淹为对应的农村居民点或城镇淹没面积，h为每户家庭房屋面积，c为农村居民点或城镇土地利用面积系数，v固为每户家庭的其他固定资产价值，r为损失率。
[0035]
本发明的有益效果为：本发明通过洪水期遥感影像得到水体范围，然后与地理矢量数据进行匹配可以快速定位至洪灾发生地，即可进行灾损评估得到灾损信息，在克服处理资源受限的条件下，实现灾害信息的快速提取，实施简单快捷，所需成本低，利用地理矢量信息的辅助功能，充分深挖不同数据内部和数据之间的协调能力，实现洪涝灾害评估，评估过程所需数据较少，分析较为简便，相比现有的评估方法效率更高，更为精确，基于0级或者1级数据进行处理，可以为减灾救灾服务提供快速的决策支持。
附图说明
[0036]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0037]
图1是本发明的方法流程图；
[0038]
图2是本发明实施例二中的研究区遥感影像图；
[0039]
图3是本发明实施例二中基于ndwi的水体提取图；
[0040]
图4是本发明实施例二中的待匹配提取水体图；
[0041]
图5是本发明实施例二中的水体匹配模板图；
[0042]
图6是本发明实施例二中的提取水体与遥感影像叠加图。
具体实施方式
[0043]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0044]
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”、“第四”等仅用于描述目的，而不能理解为指示或暗示相对重要性。
[0045]
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相
连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0046]
实施例一
[0047]
参见图1，本实施例提供了一种基于矢量引导的洪灾灾情评估方法，包括以下步骤：
[0048]
步骤一：水体提取
[0049]
基于待评估区域的遥感卫星遥感影像，利用归一化差异水体指数ndwi对待评估区域进行水体提取，其中ndwi定义如下：
[0050][0051]
水体的提取方法还包括基于阈值的多波段谱间关系法和阴影水体指数法，其中基于阈值的多波段谱间关系法主要利用影像多波段的综合优势进行水体信息的提取，包括谱间关系法与比值法，阴影水体指数法是在现有的多波段谱间关系法基础上，根据水体的光谱特性进行构建，能较好区分水体和阴影，能削弱积雪和山体裸地的影响，定义为：
[0052]
swi＝蓝波段+绿波段-近红外波段≥c(水体)
[0053]
swi＝蓝波段+绿波段-近红外波段＜c(阴影)
[0054]
式中c为经验性阈值；
[0055]
步骤二：矢量引导
[0056]
根据步骤一，先将基于ndvi提取的水体转为矢量然后导出成图片，然后将待评估区域的水体矢量数据导出成图片，其中待评估区域的水体矢量数据图片作为被匹配模板；
[0057]
步骤三：确定洪水淹没范围
[0058]
根据步骤二，先对基于ndvi提取的水体矢量图片进行模板匹配，确定具体的灾害发生地，模板匹配具体为在一副大图像中搜寻目标，已知该图中存在要找的目标，且该目标同模板有相同的尺寸、方向和图像元素，通过ssda算法在图中找到目标，确定其坐标位置，模板图像记为t，待搜索图像记为f，计算图像f(x,y)在像素(u,v)的非相似度m(u,v)为：
[0059][0060]
以m(u,v)作为匹配尺度，式中(u,v)表示的不是模板与图像重合部分的中心坐标，而是重合部分左上角像素坐标，模板的大小为：n
×
m；再将提取的水体与灾害发生地灾害发生前的卫星遥感影像进行叠加，确定洪水淹没区域信息，其中洪水淹没区域信息包括洪水淹没范围和淹没区域土地利用类型，所述灾害发生地灾害发生前的影像为从全国地理信息资源目录服务上下载的待评估区域的土地利用数据影像；
[0061]
步骤四：估算洪灾损失
[0062]
根据步骤三，根据洪水淹没区域信息对洪灾发生地的洪灾损失进行计算评估，得到灾损信息，公式如下：
[0063][0064]
其中，s为灾害发生地的直接经济损失，p
ab
是a种淹没范围内第b类财产值，t
ab
是第a种淹没级别下的第b类财产损失率，m为财产类型的总数，n是洪水淹没级别，其中洪灾损失包括农作物损失、畜牧业损失和居民财产损失，居民财产损失包括房屋建筑损失和其他固定资产损失，所述农作物损失的估算公式如下：
[0065][0066]
其中，r
农
为洪灾农作物经济损失，η
ij
为第j级淹没水深下的第i类农作物洪灾损失率，w
ij
为第j级淹没水深下的第i类农作物正常年产值，n为农作物种类数，m为淹没水深等级数，r
存
为被洪水冲走的农作物或因受淹霉变造成的损失；
[0067]
畜牧业损失的估算公式如下：
[0068][0069]
其中，na为牲畜死亡的种类数，gt为死亡牲畜的活体重量，vt为死亡牲畜的批发价格；
[0070]
房屋建筑损失的估算公式如下：
[0071]
l建筑＝a淹
×c×v×r[0072]
其中，l建筑为房屋损失价值，a淹为对应的农村居民点或城镇淹没面积，c为农村居民点或城镇土地利用面积系数，v为单位面积上的房屋价值，r为损失率；
[0073]
其他固定资产损失的计算公式如下：
[0074]
l固资＝(a淹/h)
×c×
v固
×r[0075]
其中，l固资为其他固定资产损失，a淹为对应的农村居民点或城镇淹没面积，h为每户家庭房屋面积，c为农村居民点或城镇土地利用面积系数，v固为每户家庭的其他固定资产价值，r为损失率。
[0076]
实施例二
[0077]
以洪涝灾害频发的四川省乐山市市中区作为研究对象，基于乐山市的gf1号的遥感影像，如图2所示，利用归一化差异水体指数ndwi进行水体提取，如图3所示；
[0078][0079]
将基于ndvi提取的水体转为矢量然后导出成图片，如图4所示，然后将乐山市的水体矢量数据导出成图片，如图5所示，然后对其进行模板匹配；
[0080]
经过基于矢量数据引导的模板匹配后，可以确定灾害发生地为乐山市市中区，将提取的水体与发生灾害前的影像叠加，如图6所示；
[0081]
将提取的水体与从全国地理信息资源目录服务上下载的乐山市土地利用数据叠加，然后得到洪水淹没范围的信息；
[0082]
通过对水体淹没区域的土地利用类型以及灾损区域定位进行分析，可以得到灾害
地点、淹没区域土地类型、淹没范围以及损毁价值等灾损信息，如下表1所示：
[0083]
表1灾损信息表
[0084][0085]
以示范应用为例统计各个阶段的运行时间，如表2所示，对于研究区乐山市中区进行重点区域洪水快速提取大约需要22分钟，可以为减灾救灾服务提供快速的决策支持。
[0086]
表2系统运行效率
[0087][0088]
基于矢量数据引导的洪水淹没检测方法还有其它实现方法，如利用基于正切空间的面状矢量要素形状相似性度量模型进行矢量匹配也可得到灾害发生地。
[0089]
(1)基于douglas-peucker的复杂面要素形状简化
[0090]
定义:设简化距离阈值为t；c为实平面上的闭合多边形，p0，p1，
…
，pn为该闭合多边形上顶点，并沿顺时针方向分布，计算出与p0距离最长的顶点pk(如有几个最长值取k值最小的pk)，连接p0与pk，利用直线p0pk分别对两段复合线p0－pk和pk－pn上的节点计算到直线p0pk的距离di，选取其中距离最大的点pj，如果di大于限差阈值，则保留点，反之剔除该点，利用保留的最大距离点pj将原复合线分为两段，并用同样的方法对位于它们之间的节点进行检测，重复此操作，直至节点到两端点连线的距离最大值小于限差阈值为止。
[0091]
(2)面要素形状描述方法
[0092]
假设多边形的某一顶点作为参考点p0，θ1表示起始边p0p1的方位角，φ1表示从起始边p0p1到p1p2的转角，φk表示沿着p
k－1
pk到pkp
k-1
的转角，多边形的正切空间形状描述函数为θ(l)，x轴代表从起点p0沿着多边形周边到多边形上各点pk的归一化距离y轴代表各点沿着周边的转角(以顺时针为正方向)的累加θk＝θ
k－1
+φ
k+1
(k＝3，4，
…
，n)，由于不同起始点、不同的方向所得到的面实体正切空间函数不同，在源匹配多边形中，定义p0(x0，y0)为起始点，其中x0＝max{x|(xi，yi)∈a}，y0＝min{y|x＝x0，(xi，yi)∈a}。从起始点p0沿着多边形顺时针旋转方向为正方向，p0p1为起始方向线。以p0为原点、距离阈值t为半径搜索包含在圆内的目标多边形b的节点集合q，选取与p0p1的方位角θ1差异最小的结点为目标匹配多边形的起始点p0'。
[0093]
(3)面要素形状匹配方法
[0094]
通过以上的形状描述方法对多边形a、b要素进行形状描述，将其形状化为正切空间表达式，分别为θa(l)和θb(l)，其s为x轴坐标，θ为y轴坐标。通过计算两个矢量要素之间的形状匹配距离来确定它们的相似性，进而判断两要素是否匹配。定义其匹配距离为
[0095][0096]
其中d
ab
的值越趋近于1，表示多边形a和b的形状越相似，匹配的程度越好。
[0097]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。