请选择 进入手机版 | 继续访问电脑版

利用卫星影像评估火灾后迹地

[复制链接]
蝶蝶已蝶已蝶蝶 发表于 2021-1-1 10:33:17 | 显示全部楼层 |阅读模式 打印 上一主题 下一主题
2020年3月28日19时30分,四川省凉山州木里县乔瓦镇锄头湾村与项脚蒙古族乡项脚村交界处发生森林火灾。火灾发生后,我们需要努力丈量火灾后迹地,为森林重建和植被演替创造基线。借助火灾前后的Landsat8影像,使用归一化燃烧指数盘算损失,然后将数据作为要素类发布以与其他人共享。
比力波段组合

首先,通过usgs获取 Landsat 8 影像初步检察火灾后迹地。相识研究区域后,可以更改影像的波段组合以便更好地检察燃烧区域。然后,可创建自界说组合来突出火灾后迹地。
打开工程

在开始分析之前,可下载含有任务原始数据的工程包并将其打开。

  • 转至Muli_ForestFire.zipMuli_ForestFire.z01,可下载本人分享的舆图包,由于上传限制,故分包压缩,因此只需下载Muli_ForestFire.zip和Muli_ForestFire.z01,放在同一目次下解压即可。
    Muli_Forestfire是一个舆图包。舆图包中包含用于分析的原始数据。
  • 浏览至下载的Muli_ForestFire.mpkx,双击后,ArcGIS Pro将自动导入。
    舆图包随即打开,并在内容窗格中显示两个图层。这两个图层为影像图层。

增强影像

原始下载的两个Landsat 8影像图层,分别拍摄于2019年5月21日和2020年5月7日。箭头指示区域为火灾发生区域。


  • 从USGS下载的原始图像载入时,图像大概很暗,难以看清,可以通过调解亮度、对比度和Gamma以更好地检察图像。亮度可决定图像是较亮照旧较暗。对比度可决定各个要素之间的可分辨水平。Gamma 可决定影像检测光的方式与实在际发光方式之间的关系。增加上述三项的值会提高图像的可见性。

  • 在应用程序顶部的功能区上,单击外观选项卡。在增强群组中,将图层亮度调至 20,图层对比度调至25,图层Gamma增加到1.8。按Enter键后,图像会立即更改。

  • 这里重点说一下gamma值
    Gamma 指的是栅格数据集中等灰度值之间的对比度。Gamma 对栅格数据集的玄色或白色值不会产生影响,而仅影响中间值。通过应用 gamma 校正,可以控制栅格数据集的总体亮度。别的,gamma 系数不仅会改变影像的亮度,还会改变红绿蓝之间的比率。
    小于 1 的 Gamma 值将低落较暗区域的对比度并增加较亮区域的对比度。这将使影像变暗,但不会低落影像的暗中或豁亮区域的饱和度。这有助于出现较亮要素的细节,如修建物顶部。相反,大于 1 的 gamma 值将增加较暗区域的对比度,如修建物阴影。使用高程数据时,大于 1 的 Gamma 值也将有助于出现较低高程区域的细节。
    在下图中,可以看到调解用于显示栅格数据集的 Gamma 值的效果:

  • 颠末增强后,可以更细致地检察情况。在2019年5月21日的多光谱影像中,可以清楚的看到山地、植被、河道。同比2020年5月7日火灾发生后的多光谱影像,可以或许发现植被类型和山地坡度对火灾造成的影响,尤其是火势伸张的速度。

  • 2019年和2020年该区域影像还有一处显着的区别是2019年影像有一大片云层遮盖在火灾发生区域的西北角。虽然2020年5月7日的影像可以清楚的看到迹地,但其准确边界不易确定。接下来,我们通过更改影像波段组合的方式来突出显示灾区。
TIPS:通过创建并添加栅格至镶嵌数据集,可以直接检察融合后的多光谱影像,提高地物的判读识别精度
检察差别波段组合

Landsat 影像可丈量电磁光谱的波长范围,包罗一些人眼不可见的波长。这些范围叫做光谱波段。
数值名称该波段的最佳显示1沿海气溶胶浅海、细粉尘颗粒2蓝深海、大气3绿植被4红人造物体、土壤、植被5近红外海岸线、植被6短波红外 1云层穿透度、土壤和植被湿度7短波红外 2改善的云层穿透度、土壤和植被湿度8全色好坏影像、更清晰的细节9卷云卷云10热红外 1热谱舆图、预计土壤湿度11热红外 2改善的热谱舆图、预计土壤湿度在应用程序顶部的功能区上,单击外观选项卡。可以通过波段组合模板,选择符合组合来突出显示感兴趣区域。


  • 现在影像使用的“自然色”波段组合将三个波段举行联合以大抵生成人眼可见的影像。

  • 外观选项卡的渲染组中,单击波段组合并选择彩色红外

    红色和绿色波段突出显示了植被,这对于检察火灾很有资助,因为未遭受火灾的茂密植被区域和植被已遭火灾摧毁的区域有着显着对比。使用另一波段来突出显示植被,如近红外线(波段 5)可提高对比度。在该图像中,植被显示为红色。两个火灾区域显示为暗棕色。与原始图像相比,火灾显示更加清晰
  • 外观选项卡的渲染组中,单击波段组合并选择陆地/水域交界面

    该组合的主要目的是为了刻画陆地和水域,同时穿透霾雾(或是火灾发生时产生的烟雾)。本影像中,火灾周围看不到烟雾,边界变得更加清晰。但是,燃烧区域显示为橙色而其周边山坡显示为黄色。这使得已伸张至山脉中的迹地变得更加难以分辨。
  • 再次单击波段组合,并选择植被分析

    该组合将彩色红外组合的植被突出显示与陆地/水域交界面组合的霾雾穿透度联合起来。火灾迹地与山地的肴杂部门已不像之前的组合那么严重。
创建自界说波段组合

以上是使用预配置的波段组合,现在我们选择自己来创建自界说的波段组合,改进“植被分析”波段组合。

  • 选中要变更的图层,“植被分析”波段组合使用短波红外 1 波段来淘汰霾雾,并使用近红外和红色波段来突出显示植被。从短波红外 1 波段切换至短波红外 2 波段可改善霾雾(或云层)穿透度。

  • 波段自动变革的同时,舆图上的图像也会随之变革。由于变革发生在两个相似波长之间,因此舆图上的区别并不十分显着,但确实有所变革。下一步,用蓝色波段替换红色波段。

    舆图上的图像发生轻微变革。但已经可以清晰的展示火灾发生区域。而且此种波段组合能较好的淘汰透霾雾。

要想以后能轻松访问此波段组合,可以将其生存为“外观”选项卡中的默认组合。
只管此波段组合更适合生成准确的火灾后迹地边界,但还需要对其举行一些目视解译以确定哪些区域已燃烧,哪些区域未燃烧。接下来,使用燃烧指数的数学公式来量化盘算燃烧面积,此公式可提供受灾区域的更准确丈量。
盘算燃烧指数

归一化燃烧比率以量化形式来确定燃烧面积,是使用数学方式对近红外波段和短波红外 2 波段(分别为波段 5 与波段 7)举行比力,以确定火伤严重性。然后,比力 2019 年与 2020 年影像的 NBR 以盘算 NBR 的更改,仅显示在两张图像拍摄日期之间发生燃烧的区域。
需要注意的是,NBR盘算前应将影像更正为反射率值。因此需要大气校正。
盘算归一化燃烧比率

需要对 NBR 举行两次盘算:一次盘算 2019年图像,另一次盘算 2020 年图像。要举行盘算,可以使用栅格盘算器地理处理惩罚工具,并使用以下等式:
NBR =(NIR - SWIR2)/(NIR + SWIR2)
使用ArcGIS Pro影像指数模板盘算方法

ArcGIS Pro2.5版本开始逐渐丰富了影像指数模板,可以通过该方式快速生产需要的成果。

  • 在功能区上,单击影像选项卡,然后单击指数窗格,单击NBR

  • 选择对应的波段,近红外为Band_5,短波红外为Band_7,单击确定

  • 得到NBR盘算效果。

使用栅格盘算器盘算方法

使用栅格盘算器需要将波段5(近红外)和波段7(短波红外2)加载到内容窗格中。

4. 在目次窗格中,找到影像.tif并展开,右键单击波段5(NearInfrared )并选择添加至当前舆图
波段随即添加到舆图。看上去像是好坏图像,因为该图像未作为 RGB 合成的一部门举行显示。单波段图像图层通常看似好坏图像。
同理添加波段7(ShortWaveInfrared_2)

5. 单击功能区上的分析选项卡,然后单击工具

地理处理惩罚窗格中,单击搜索框并输入栅格盘算器。从效果列表中,单击栅格盘算器(空间分析工具)栅格盘算器(图像分析工具)(两者都可以)

6. 在栅格下,双击 Band_5,将其添加到表达式框中,然后在运算符下,双击减号将其添加至表达式。然后双击Band_7。


  • 完成的表达式 (“Band5” - “Band7”) / ((“Band5” + “Band7”)。

  • 单击运行,得到20200507的NBR成果。同理完成20190521的NBR成果。

    TIPS:有时NBR大概NDVI等归一化指数盘算的效果大概出现[-1,1]以外的异常值出现,这与大气校正时所选的模子有一定关系,这些像元通常少少,不敷以影响。
盘算NBR的变革

再一次使用栅格盘算器工具来盘算两图像之间的 NBR 变革。使用此工具,将移除 2019 至 2020 年间未燃烧区域的值,仅显示燃烧区域并移除图像上的所有其他内容。

  • 将盘算的NBR效果加载到内容窗格中,点击栅格盘算器工具。在“舆图代数”表达式框中,创建表达式 “NBR_20190521_ROI_T.tif” - “NBR_20200507_ROI_T.tif”,单击运行。

  • 得到的变革效果可以看到,火灾位置看起来险些为纯白色,与周围的灰色和玄色区域形成显着对比。

  • 打开符号系统窗格。将色带更改为绿色到红色方案。可以沿着黄色边界举行区划,得到受灾的全部区域。

    根据差值的效果,我们可以对其举行阈值分割,能给得到差别受伤水平的植被分类。这里不再详述。

来源:https://blog.csdn.net/weixin_43809925/article/details/111909097
免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

发布主题

专注素材教程免费分享
全国免费热线电话

18768367769

周一至周日9:00-23:00

反馈建议

27428564@qq.com 在线QQ咨询

扫描二维码关注我们

Powered by Discuz! X3.4© 2001-2013 Comsenz Inc.( 蜀ICP备2021001884号-1 )