MODIS数据火点提取方法

MODIS数据火点提取
【数据准备】
所需数据：MODIS02数据
注意下载1KM的数据（因为需要用到红外波段，而红外波段的空间分辨率较低）
网站：https://ladsweb.modaps.eosdis.nasa.gov/search/imageViewer

这里我们需要用的是Terra数据，点击Select a Sensor - MODIS:Terra,选择MOD021KM。

选择所需的TIME时间范围，点击ADD DATE。

点击LOCATION位置范围，右侧有不同的位置选择方式，标明需要哪个位置范围的影像。

点击IMAGES，查看筛选结果（右下角标明影像的地理范围），点击右侧对话框的Files进行下载。

【计算亮度温度值】
ENVI系统默认对MODIS数据进行定标，所以通过File->Open As-EOS->MODIS打开的是经过定标的数据集，分别是反射率（reflectance）、辐射率（radiance）以及发射率（Emissive），热红外数据被自动定标为发射率，而有时候我们需要用到热红外波段的辐射亮度数据（如计算亮度温度），所以需要读取原始的热红外数据DN值，进行手动定标。下面介绍读取原始DN值和手动定标为辐射亮度值的过程。
（1）原始热红外数据集的读取
ENVI EX打开File->Open As->Generic formats->HDF4，选择hdf文件，选择发射率数据集，点击OK，选择数据存储方式BSQ，点击OK。

打开的是热红外原始数据集，第20-36波段，共16个波段，分别是：20、21、22、23、24、25、27、28、29、30、31、32、33、34、35、36波段。

（2）辐射亮度定标
打开/Raster Management/Data-Specific Utilities/View HDF Dataset Attributes，选择原始的hdf数据，点击OK，选择相应的热红外数据集：Earth View 1KM Emissive Bands Scaled Integers，点击OK，得到一个文件说明面板：

radiance_scales，和radiance_offset这两项参数代表波段的增益和偏移量，是辐射定标的系数。比如要计算31波段的辐射亮度，读取到scales为0.00084002，offsets为1577.33972168，带入MODIS辐射定标的通用公式：Radiance=scales*(DN-offsets)，即可以得到该波段的辐射亮度，单位见说明面板下方。单位容易换算错，一定要注意。
打开/Band Algebra/Band Math工具，输入公式：0.00084002*(B31-1577.33972168)，（B31是第31波段的DN值），点击OK，选择第31波段数据为B31，设置路径和文件名，点击OK。得到的结果就是31波段的辐射亮度。

同样的方法可以得到其他任意波段的辐射亮度。算出21、31波段的辐射亮度。
（3）计算波段21、31的亮度温度值
亮温是指辐射出与观测物体相等辐射能量的黑体温度，可以根据普朗克公式计算得到，MODIS第i波段的亮度温度Ti由下式计算：

其中 i是波段i的中心波长。
针对MODIS的第31,21波段,其值可分别取31=11.28,21=3.959m；
C1和C2是常量,C1=1.19104356×10-16W·m2,C2=1.4387685×104m·K.
由于需要注意单位转化问题，为了便于计算，可将该式进行简化，设K2=C2/ i,K1=C1/ i^5。则上式可转化为：

对于31波段，K1=729.54163436, K2=1304.41387;
对于21波段,K1=122461.38, K2=3634.17;
【火点提取】
采用如下筛选规则：
T4>Mean(T4)+3*StdDev(T4)
AND
T4-T11>Mean(T4-T11)+3*StdDev(T4-T11)
由于结果中火点像素数很少，不便观察，因此对结果做了一个膨胀（在ENVI CLASSIC中点击Filter-Convolutions and Morphology-Morphology-Dilate），得到结果如图，红色表示火点：

并计算CDI归一化云检测指数，进行去云处理：

CDI大于0即为云。
提取结果如图：白色为云和水体。

去云后的火点提取数据为：
白色为从上一步火点中剔除的云，红色为火点

下载MOD14温度异常数据并进行对比（左图为我们的结果，右图为modis温度异常数据图）

可以看到火点与温度异常点分布较为一致，但是仍有一些温度正常但被划分为火点的像素。

火点提取结果可与https://worldview.earthdata.nasa.gov网站上的火点数据进行对比
参考资料：
MODIS产品功能解释：https://wenku.baidu.com/view/6fd329dcf524ccbff0218440.html
MODIS辐射亮度计算：http://blog.sina.com.cn/s/blog_764b1e9d0102w1xy.html
MODIS云检测：https://wenku.baidu.com/view/f419031b964bcf84b9d57b5d.html
MODIS数据反演地表温度：https://wenku.baidu.com/view/8d6b9c66caaedd3383c4d3f4.html
MODIS林火识别算法的验证分析 ,周小成汪小钦

MODIS数据火点提取方法相关推荐

使用MODIS数据筛选火点
流程图 1.数据源1)MODIS021KM (2)GF-4号数据(GF-4PMS) (3)MOD13(植被指数) (4)MOD14火灾产品在ENVI软件中打开modis数据在这里插入图片描述 2. ...
基于Modis数据监测森林火灾
2月10日凌晨,四川省凉山州木里县三桷桠乡里铺村和高房子村交界处发生森林火灾.坐标为101°29′21.47″E,28°04′20.40″,明火于12日9时26分全部扑灭. 为对此次灾害进行监测分析, ...
Python应用实战案例-Python使用MODIS数据实现温度植被干旱指数TVDI的计算
1.数据下载数据及代码参见温度植被干旱指数TVDI 采用的数据为MODIS植被指数产品MOD13A3.地表温度产品MOD11A2以及SRTM DEM产品. MODIS数据来源于美国航空航天局(Nat ...
modis数据介绍_【更新90篇】地理数据科学技术文章合集，欢迎大家点赞、在看、转发三连！...
最近因微信公众号修改了规则, 不将公众号标星可能无法及时收到推送, 还请大家把"走天涯徐小洋地理数据科学"加个星标. 公众号技术文章合集目录,目前已整理以下8个部分,共计90篇文章 ...
用python批量下载modis数据的速度怎么样_MODIS数据的简介和下载（五）——应用密钥的Python脚本下载...
接上文,由于NASA官方弃用了FTP改用HTTPS服务,所以关于MODIS数据的下载方式有所改变.完整系列博客可以参照以下链接.本篇主要接着上一篇没有讲完的应用密钥的脚本下载介绍. 1 官方教程与说明 ...
MODIS数据的简介和下载（三）——MODIS数据下载方式（基于MODIS Web Service）
版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/ESA_DSQ/article/details/70176214 这是MODIS数据的简介和下载的最后 ...
用MODIS数据借助MATLAB绘制世界植被分布图
目录下载数据下载文件信息查看软件HDFview(可选) 用MATLAB导入并读取数据下载数据可以下载数据的网址有很多,本教程以美国NASA的MODIS数据为例,下载链接在这里选择想要下载数据 ...
关于批量下载MODIS数据的坑
关于批量下载MODIS数据的坑与正确姿势 1. 第一个坑 1.1 Ladsweb的坑由于之前做项目的时候去Ladsweb官网下载数据的时候还好好的,于是经过被批评改正之后回来就信心满满地去Ladsw ...
基于MODIS数据的NDVI与LST相关性分析（IDL代码实现）
1 数据预处理 (1)数据提取我们可以选取2018年5月初华东地区MODIS中的MOD11A2和MOD13A2的16天合成LST和NDVI产品数据,下载地址:MODIS数据下载网站下载数据需要注册 ...
MODIS数据下载，HEG批量拼接处理，arcmap批量掩膜提取，python栅格计算
前言: 2022.12.30修改:补Python编程基础的过程中发现使用Print也可以向文件里写东西,查看写过的代码,然后发现HEG批量拼接处理代码中有的地方在复制粘贴至博客的时候出现了缩进错误,评 ...

MODIS数据火点提取方法

MODIS数据火点提取方法相关推荐

最新文章

热门文章