火点识别简单实现

大家好，欢迎来到IT知识分享网。

一、实验数据下载

以澳大利亚数据下载为例，在NASA网站上具体下载方法如下：

网址：Find Data – LAADS DAAC

①选择Terra传感器，MOD021KM卫星数据，选择卫星产品

②根据需要设置卫星产品的时间

③选择下载的区域，可以根据国家（Countries）格网（Tiles）等选择

④如果只下载一副影像，直接点击下载即可，下载多幅影像可以用wget批量下载

二、火点识别

1、火点识别思路

       基本思路：假设一个或多个光谱通道的亮度温度超过预先设定的阈值，那么该像素就是火。因此使用阈值法。温度异常：亮温法——>数值提取

T(4微米)21/22波段>360K(夜间>330K)

其他方法：目视解译：选用波段CH7、2、1（或者  21、2、1）进行RGB彩色合成

鲜红色：正在燃烧的火点；暗红色：过火点；绿色：未燃烧的植被；深蓝色：水体；灰色：烟雾或云；黄褐色：裸露土地

2、火点识别ENVI实现

①辐射亮度计算

ENVI系统默认对MODIS数据进行定标，所以通过File->Open As-EOS->MODIS打开的是经过定标的数据集，分别是反射率（reflectance）、辐射率（radiance）以及发射率（Emissive），热红外数据被自动定标为发射率，需要读取原始的热红外数据DN值，进行手动定标。

【读取原始的热红外数据DN值工具：

Toolbox——Raster Management——Data-Specific Utilities——View HDF Dataset Attributes】

 读取原始的热红外数据DN值工具与结果

【利用BandMath工具进行辐射亮度定标计算：0.00006922*(b22-2730.)】

辐射亮度定标计算

辐射亮度结果图

注：如何打开热红外发射数据集？

①MODIS数据直接打开，只包含反射数据集：Open As——Optical Sensors——EOS——MODIS——选择下载的hdf数据，即可打开MODIS数据

 反射率波段7、2、1彩色合成

 ②打开热红外发射波段集：Open As——Generic Formats——HDF4——选择下载的hdf数据——Earth View 1KM Emissive Bands Scaled Integers（1KM地球视图发射波段缩放整数）——选择BSQ格式，即可打开热红外数据集

②亮温计算

亮温是指辐射出与观测物体相等辐射能量的黑体温度，可以根据普朗克公式计算得到，根据前人研究的结果，可以用如下公式计算22波段的亮温：3634.17/alog(1+.38/b22)

【BandMath工具：Toolbox——Band Algebra——Band Math】

 亮温计算方法 亮温计算结果图

③火点识别与提取

利用掩膜方法提取火点，将亮温大于310的地方定义为火点（选择的310作为临界点，是希望效果明显点）。

【掩膜工具：Toolbox——Raster Management——Masking——Build Mask】

掩膜提取

火点提取结果图

1. （学习中的小菜鸟，感谢指出不足，但是请轻喷，谢谢♥）