归一化水指数模型应用洪灾监测技术方案设计Word下载.docx
《归一化水指数模型应用洪灾监测技术方案设计Word下载.docx》由会员分享,可在线阅读,更多相关《归一化水指数模型应用洪灾监测技术方案设计Word下载.docx(7页珍藏版)》请在冰豆网上搜索。
![归一化水指数模型应用洪灾监测技术方案设计Word下载.docx](https://file1.bdocx.com/fileroot1/2022-11/24/5e0eaad9-4694-4590-8da5-938977097816/5e0eaad9-4694-4590-8da5-9389770978161.gif)
应用遥感技术监测洪涝灾害的技术方法和流程有很多,下面介绍一种较为实用和简单的方法。
归一化水指数模型
Mcfeeters提出的归一化差异水体指数(简称归一化水指数)NDWI(NormalizedDifferenceWaterIndex)便是其中之一,其公式如下:
NDWI=(BGreen-BNIR)/(BGreen+BNIR)
式中,BGreen代表绿光波段;
BNIR代表近红外波段。
模型的基本原理:
由于水体的反射从可见光到中红外波段逐渐减弱,在近红外和中红外波长范围内吸收性最强,几乎无反射。
因此用可见光波段和近红外波段的反差构成的NDWI可以突出影像中的水体信息(水体的NDWI值大)。
另外由于植被在近红外波段的反射率一般最强,因此采用绿光波段与近红外波段的比值可以最大程度地抑制植被的信息,从而达到突出水体信息的目的。
这也是Mcfeeters选用近红外波段构建该指数的原因。
Mcfeeters模型考虑到植被的影响,后来很多学者又将土壤、建筑物等因素考虑进去,提出一些改进型的归一化水指数(徐涵秋):
MNDWI=(BGreen–BSWIR)/(BGreen+BSWIR)
BSWIR表示中红外波段。
Wilson等又提出(BNIR–BSWIR)/(BNIR+BSWIR)用于研究植被的水分及土壤湿度等。
技术流程
基本原理与流程:
分别计算洪灾前影像和洪灾后影像的NDWI,对两个时相的NDWI做差值运算,选择阈值对NDWI差值图像进行分类,如NDWI差值为正表示水体增加,可认为洪水淹没区。
详细流程如下图:
图1监测流程图
图像配准:
两个时相影像精确的空间配准是整个流程的前提。
NDWI计算:
可根据数据情况(如:
是否包括中红外波段)和主要地表覆盖类型(植被、建筑物、裸土等)选择不同的计算模型。
分类阈值确定:
理想情况是NDWI差值大于0为水体淹没区,由于其他地表覆盖类型的干扰,这个阈值往往不为0。
可选择一些算法自动计算,再结合目视预览效果进行调节。
图像分类及后处理:
选择大于“分类阈值”的区域作为输出结果。
最初结果中往往包括一些“小图斑”,选择性的对这些小图斑进行剔除。
为了满足应急救灾的需求,整个流程应该在统一界面下流程化完成。
详细实现步骤
直接使用ENVIEX下的遥感变化监测工具ImageChange完成整个流程。
ImageChange工具采用向导式操作,详细过程如下:
(1)启动ENVIEX,并打开两个时相的影像文件。
(2)双击Workflow->
ImageChange工具,分别为Time1File选择灾前影像,Time2File选择灾后影像,单击Next进入ImageCo-Registration面板。
图2选择前后时相的影像
(3)这一步选择是否执行图像配准。
如果图像没有精确配准,可选择自动选择同名点执行图像配准。
单击Next进入ChangeMethodChoice面板
(4)在ChangeMethodChoice面板中,选择ImageDifference,单击Next进入ImageDifference面板。
图3图像配准选项
(5)在ImageDifference面板中,选择DifferenceofFeatureIndex,之后选择WaterIndex(NDWI)。
(6)在Advanced选项中,默认为Band1和Band2选择了相应的波段。
这里计算公式为:
(Band2-Band1)/(Band2+Band1)。
根据数据的实际情况选择波段。
(7)单击Next进入ThresholdingorExport面板,直接选择ApplyThresholding后,单击Next进入ChangeThresholding面板。
图4计算NDWI及差值
(8)在ChangeThresholding面板中,选择感兴趣的变化:
IncreaseOnly(增加)。
提供了四种自动计算阈值方法,选择Otsu’s,勾选Preview预览结果分类结果。
在Manual选项中可以手动调整阈值。
(9)单击Next进入ClearUp面板。
图5确定分类阈值
(10)在Cleanup面板中,提供平滑滤波和像素聚合处理,将结果中的“小图斑”去除或者合并。
(11)单击Next进入Export面板,结果可输出为矢量或者单波段的分类结果。
图6分类后处理
图7输出结果