快速提取过火区域,卫星影像评估西昌山火方法初探
转自易智瑞中国
2020年3月30日15时,四川凉山州西昌市突发森林火灾,火势迅速蔓延,危及西昌市区。截至3月31日16时,火灾过火面积1000余公顷,毁坏面积为80余公顷。更加令人心痛的是,该起山火造成19名地方扑火人员牺牲。
易智瑞公司十分关注火情的发展,并在第一时间派出了无人机专家队伍前往火灾现场驰援:四川凉山森林火灾 易智瑞无人机专家赶赴现场驰援。
关于这场山火的损失评估,目前媒体的正式报道还停留在3月31日的数据。如何快速评估这场山火的严重性?小编尝试对公开可获取的卫星影像数据进行分析,来快速提取过火区域。以下是具体的分析过程,考虑不周之处敬请指正:
1、数据获取
在数据选择上,考虑到Landsat8卫星的重访周期为16天,哨兵2的为10天,A星和B星结合起来也就5天。相比较而言,哨兵2的数据比较充足,且波段融合后的影像精度可达10m,比Landsat8全色波段融合后15m更好,于是选用哨兵2卫星的数据进行西昌山火的过火区研究。
哨兵2数据的获取有很多途径,推荐从哥白尼数据中心(Copernicus Open Access Hub)去获取,具体操作这里就不介绍了。
地址为:https://scihub.copernicus.eu/
这里下载的是3月25日A星L1C级和4月9日B星L1C级卫星数据,能够充分的观测山火前后的状态。哨兵2卫星数据提供L1C和L2A级产品,L1C级进行了几何精校正而没有进行大气校正和辐射定标,L2A是L1C进行上述操作后的产品。但这里没直接用L2A的,用的是对L1C级产品进行预处理后的L2A级产品。
哨兵2产品类型
2、数据预处理
整个预处理是利用欧空局的SNAP进行处理的,同时利用Sen2Cor插件进行大气校正和辐射定标。Sen2Cor可以在SNAP中直接安装,也可以下载个单机版,再命令行执行,具体的安装及命令行使用方法就不展开说了。
附带下载地址:http://step.esa.int/main/third-party-plugins-2/sen2cor/
(1) 大气校正和辐射定标
这里用的是SNAP中的Sen2Cor插件,该插件需要在软件中安装。大气校正和辐射定标操作很简单,这里注意一点Resolution的选择,可以按照需要对不同精度的波段进行大气校正和辐射定标,如下图:
大气校正和辐射定标
(2) 重采样
这里的重采样就是将对应的波段数据重采样为想要的精度,当然我们希望将60m和20m分辨率的波段重采样为10m的,哨兵数据的可见光波段都是10m的,一步操作就能完成重采样,如下图所示有两种方法:
重采样方法一
重采样方法二
最后看看重采样后的结果,精度很高:
(3)选择波段
哨兵2数据经过大气校正后多了好多波段,比如气溶胶厚度云检测等等,但这里不需要那么多,所以可以选取需要的波段进行保存,然后进行后续处理,试了一下全部输出要20GB左右,还是选择常用的波段输出。
波段选择
3、过火区提取
将波段提取后的结果保存为TIFF后,这里选用ArcGIS Pro进行过火区提取。ArcGIS Pro的影像处理功能也很强大,对于大范围的影像数据,进行波段运算等相关操作只需一眨眼功夫就可以看到结果。
对于过火区的提取,经过查找相关论文后选择了三种方法进行提取和对比,下面简单说一下具体的提取方法。
在进行过火区提取前,手动绘制了泸山的区域,在Pro中进行泸山区域的影像裁剪工作,方便后面的提取。
裁剪
泸山片区
假彩色(843)合成山火前后对比,假彩色合成植被凸显为红色
(1)归一化植被指数(NDVI)法
NDVI = (NIR - R)/(NIR + R)
dNDVI = NDVI1 - NDVI2
这里对山火前和山火后的数据进行NDVI计算,用山火前的减去山火后的NDVI,从中提取过火区,然后利用阈值法对dNDVI进行处理,从而提取过火区。公式中NDVI1是山火前的NDVI,NDVI2是山火后的NDVI,后面也是如此。Pro中的相关操作如下图:
NDVI计算
dNDVI计算
阈值分割
(2)归一化火烧指数法(NBR)
NBR = (NIR - SWIR)/(NIR + SWIR)
dNBR = NBR1 - NBR2
NBR是利用近红外和短波红外波段构建的植被指数,利用过火区光谱反射率在短波红外波段上升和在近红外波段下降的特征,能够将过火区与其他地物分离。计算过程和上面的一样,只不过波段计算时需要自定义,如下图:
NBR计算
(3)过火区识别指数(BAI)法
BAI=1/((0.1-R)(0.1-R)+(0.06-NIR)(0.06-NIR))
dBAI = BAI1 - BAI2
BAI是基于地物像元值与参考光谱值的距离构建的,通过参考相关文献,将近红外和红光波段反射率的参考值分别为设为0.06和0.10,能够强调过火区的木炭信号,计算步骤和NBR的相同。
4、结果对比
首先,看一下dNDVI、dNBR、dBAI的计算结果,从结果看出他们都能较好地提取过火区,而dNBR计算后的结果更容易区分过火区。
接下来,是阈值分割后的结果,对于阈值分割,其阈值的选取尤为重要,由于对具体的山火区域不太了解,所以阈值选取也会造成误差,不过此处仅为研究。下图为阈值分割结果。
对阈值分割过的结果,会出现很多细小零散的区域,这里采用的是形态学算子进行剔除细小图斑,于是在ArcMap中进行形态学处理,这里采用的是闭运算,在矢量化模块中进行处理,如下图:
闭运算前后对比
最后,利用ArcGIS Pro中的栅格转矢量将闭运算后的数据转换为矢量数据,用字段计算器计算其面积,获取过火区的面积,此处没考虑地形。过火区及其面积如下图所示:
5、最后的验证
提取后发现,结果与官方3月31日发布的数据有点出入,下面为3月31日官方发布的过火区:
这是4月9日自己提取后的:
左上角差距倒不大,差距主要在右下角,需要验证是否提取有误。于是,找了右下角区域处的一座山,卧云山,如下图:
然后查找山火期间这座山有没有相关报道,结果如下图显示,4月1日消防官兵正在扑灭此山起火点上的余火。上面官方发布过火区是在3月31日,当时山火还没结束,可见提取结果具有可信度。
这里写自定义目录标题
欢迎使用Markdown编辑器
你好! 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章,了解一下Markdown的基本语法知识。
新的改变
我们对Markdown编辑器进行了一些功能拓展与语法支持,除了标准的Markdown编辑器功能,我们增加了如下几点新功能,帮助你用它写博客:
- 全新的界面设计 ,将会带来全新的写作体验;
- 在创作中心设置你喜爱的代码高亮样式,Markdown 将代码片显示选择的高亮样式 进行展示;
- 增加了 图片拖拽 功能,你可以将本地的图片直接拖拽到编辑区域直接展示;
- 全新的 KaTeX数学公式 语法;
- 增加了支持甘特图的mermaid语法1 功能;
- 增加了 多屏幕编辑 Markdown文章功能;
- 增加了 焦点写作模式、预览模式、简洁写作模式、左右区域同步滚轮设置 等功能,功能按钮位于编辑区域与预览区域中间;
- 增加了 检查列表 功能。
功能快捷键
撤销:Ctrl/Command + Z
重做:Ctrl/Command + Y
加粗:Ctrl/Command + B
斜体:Ctrl/Command + I
标题:Ctrl/Command + Shift + H
无序列表:Ctrl/Command + Shift + U
有序列表:Ctrl/Command + Shift + O
检查列表:Ctrl/Command + Shift + C
插入代码:Ctrl/Command + Shift + K
插入链接:Ctrl/Command + Shift + L
插入图片:Ctrl/Command + Shift + G
查找:Ctrl/Command + F
替换:Ctrl/Command + G
合理的创建标题,有助于目录的生成
直接输入1次#,并按下space后,将生成1级标题。
输入2次#,并按下space后,将生成2级标题。
以此类推,我们支持6级标题。有助于使用TOC
语法后生成一个完美的目录。
如何改变文本的样式
强调文本 强调文本
加粗文本 加粗文本
标记文本
删除文本
引用文本
H2O is是液体。
210 运算结果是 1024.
插入链接与图片
链接: link.
图片:
带尺寸的图片:
居中的图片:
居中并且带尺寸的图片:
当然,我们为了让用户更加便捷,我们增加了图片拖拽功能。
如何插入一段漂亮的代码片
去博客设置页面,选择一款你喜欢的代码片高亮样式,下面展示同样高亮的 代码片
.
// An highlighted block
var foo = 'bar';
生成一个适合你的列表
- 项目
- 项目
- 项目
- 项目
- 项目1
- 项目2
- 项目3
- 计划任务
- 完成任务
创建一个表格
一个简单的表格是这么创建的:
项目 | Value |
---|---|
电脑 | $1600 |
手机 | $12 |
导管 | $1 |
设定内容居中、居左、居右
使用:---------:
居中
使用:----------
居左
使用----------:
居右
第一列 | 第二列 | 第三列 |
---|---|---|
第一列文本居中 | 第二列文本居右 | 第三列文本居左 |
SmartyPants
SmartyPants将ASCII标点字符转换为“智能”印刷标点HTML实体。例如:
TYPE | ASCII | HTML |
---|---|---|
Single backticks | 'Isn't this fun?' | ‘Isn’t this fun?’ |
Quotes | "Isn't this fun?" | “Isn’t this fun?” |
Dashes | -- is en-dash, --- is em-dash | – is en-dash, — is em-dash |
创建一个自定义列表
-
Markdown
- Text-to- HTML conversion tool Authors
- John
- Luke
如何创建一个注脚
一个具有注脚的文本。2
注释也是必不可少的
Markdown将文本转换为 HTML。
KaTeX数学公式
您可以使用渲染LaTeX数学表达式 KaTeX:
Gamma公式展示 Γ ( n ) = ( n − 1 ) ! ∀ n ∈ N \Gamma(n) = (n-1)!\quad\forall n\in\mathbb N Γ(n)=(n−1)!∀n∈N 是通过欧拉积分
Γ ( z ) = ∫ 0 ∞ t z − 1 e − t d t . \Gamma(z) = \int_0^\infty t^{z-1}e^{-t}dt\,. Γ(z)=∫0∞tz−1e−tdt.
你可以找到更多关于的信息 LaTeX 数学表达式here.
新的甘特图功能,丰富你的文章
- 关于 甘特图 语法,参考 这儿,
UML 图表
可以使用UML图表进行渲染。 Mermaid. 例如下面产生的一个序列图:
这将产生一个流程图。:
- 关于 Mermaid 语法,参考 这儿,
FLowchart流程图
我们依旧会支持flowchart的流程图:
- 关于 Flowchart流程图 语法,参考 这儿.
导出与导入
导出
如果你想尝试使用此编辑器, 你可以在此篇文章任意编辑。当你完成了一篇文章的写作, 在上方工具栏找到 文章导出 ,生成一个.md文件或者.html文件进行本地保存。
导入
如果你想加载一篇你写过的.md文件,在上方工具栏可以选择导入功能进行对应扩展名的文件导入,
继续你的创作。
注脚的解释 ↩︎