买专利卖专利找龙图腾,真高效! 查专利查商标用IPTOP,全免费!专利年费监控用IP管家,真方便!
申请/专利权人:中国科学院地理科学与资源研究所
摘要:本发明提供一种城市洪涝灾害经济损失预测方法及装置,涉及遥感监测技术领域,所述方法包括:使用初始修正因子,计算目标区域内各个像元的水体指数值;将各个像元的水体指数值与实际的淹水区域分布进行对比评估,以得到评估结果;根据评估结果,调整初始修正因子的方向和幅度,以获得最终的修正因子;根据最终修正因子重新计算目标区域的水体指数,以得到最终的水体指数;根据目标区域灾前和灾后的遥感图像、最终的水体指数,提取淹水区域;将提取的淹水区域与农田分布图进行叠加分析,计算受灾农田面积。本发明可以使提取的淹水区域更加准确,并与农田分布图进行叠加分析,从而更准确地计算受灾农田面积。
主权项:1.一种城市洪涝灾害经济损失预测方法,其特征在于,所述方法包括:获取目标区域灾前和灾后的遥感图像,以及目标区域的地表特征、大气条件以及遥感数据的特性;根据目标区域的地表特征、大气条件以及遥感数据的特性,分析目标区域的特点,确定影响水体指数计算的因素;根据影响水体指数计算的因素,设定一个初始修正因子作为计算的起点;使用初始修正因子,计算目标区域内各个像元的水体指数值;将各个像元的水体指数值与实际的淹水区域分布进行对比评估,以得到评估结果;根据评估结果,调整初始修正因子的方向和幅度,以获得最终的修正因子;根据最终修正因子重新计算目标区域的水体指数,以得到最终的水体指数;根据目标区域灾前和灾后的遥感图像、最终的水体指数,提取淹水区域;将提取的淹水区域与农田分布图进行叠加分析,计算受灾农田面积,并通过云计算平台以及预设的经济数据,计算经济损失;使用初始修正因子,计算目标区域内各个像元的水体指数值,包括:确定目标区域在遥感影像中的位置,获取目标区域行列号范围;从第一个像元开始,按照行优先的顺序,逐像元读取对应的绿光波段和近红外波段的反射率值;对于每个像元,提取对应的绿光波段反射率值Gr,提取每个像元对应的近红外波段反射率值NI;根据绿光波段反射率值Gr、近红外波段反射率值NI以及初始修正因子,计算水体指数ND;将水体指数ND赋值给当前像元,作为当前像元的水体指数值;继续处理下一个像元,重复操作,直到目标区域内的所有像元都计算完毕;水体指数ND的计算公式为: ;其中,Fo表示初始修正因子,表示与地表高程相关的参数,表示与太阳天顶角相关的参数,表示与大气水汽含量相关的参数,表示地表高程的值,W表示大气水汽含量的值,S表示太阳光线与地表法线之间的夹角,表示太阳天顶角的余弦值;将各个像元的水体指数值与实际的淹水区域分布进行对比评估,以得到评估结果,包括:获取实际的淹水区域分布数据,将实际淹水区域分布数据转换为与水体指数计算结果相同的空间参考和分辨率;确定水体指数的阈值,将水体指数图按照确定的阈值进行二值化,得到水体和非水体的分类结果;将水体指数分类结果与实际淹水区域分布进行空间叠加分析,以得到空间叠加分析结果;根据空间叠加分析结果,计算评估指标;根据评估指标,分析水体指数分类结果的准确性,以得到评估结果;根据评估结果,调整初始修正因子的方向和幅度,以获得最终的修正因子,包括:计算每个像素的水体指数ND与实际淹水区域A之间的差异D,其中,Di=NDi-Ai,i表示像素索引,NDi是第i个像素的水体指数,Ai是第i个像素的实际淹水情况,Di表示第i个像素的水体指数与实际淹水情况之间的差异;根据水体指数ND与实际淹水区域A之间的差异D,计算加权标准差SD;根据加权标准差以及初始修正因子,计算最终的修正因子;加权标准差SD的计算公式为: ;其中,Di表示第i个像素的水体指数与实际淹水情况之间的差异;Wi表示第i个像素的权重,表示第i个像素所代表的地物类型Ti的影响系数,i和j表示索引,Dj表示第j个像素的水体指数与实际淹水情况之间的差异,表示第j个像素所代表的地物类型Tj的影响系数,Wj表示第j个像素的权重;最终的修正因子的计算公式为: F n =Fo±k×SD;其中,Fn表示最终的修正因子,k表示调整系数,如果差异D>0,则Fn=Fo-k×SD;如果差异D<0,则Fn=Fo+k×SD。
全文数据:
权利要求:
百度查询: 中国科学院地理科学与资源研究所 一种城市洪涝灾害经济损失预测方法及装置
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。