ENVI 系统默认对 MODIS 数据进行定标,所以通过 File->Open As-EOS->MODIS 打开的
是经过定标的数据集,分别是反射率(reflectance)、辐射率(radiance)以及发射率
(Emissive),热红外数据被自动定标为发射率,而有时候我们需要用到热红外波段的辐射
亮度数据(如计算亮度温度),所以需要读取原始的热红外数据 DN 值,进行手动定标。下
面介绍读取原始 DN 值和手动定标为辐射亮度值的过程。本文引用网上资源
第一步,原始热红外数据集的读取
打开 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 波段。
第二步,辐射亮度定标
打开/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
波段的辐射亮度。
同样的方法可以得到其他任意波段的辐射亮度。