MODIS 温度反演常用工具（如 ENVI、IDL 代码）的操作要点.docxVIP

MODIS温度反演常用工具（如ENVI、IDL代码）的操作要点

以下是MODIS温度反演中主流工具（ENVI、IDL）?的实操要点，结合前文的反演步骤，聚焦关键操作流程和核心代码/参数设置，兼顾新手易上手性和实用性：

一、ENVI操作流程（可视化界面，适合非编程用户）

ENVI是遥感数据处理的主流软件，针对MODIS温度反演的核心操作围绕“数据预处理→大气校正→比辐射率计算→温度反演”展开，推荐使用ENVI5.3及以上版本（内置MODIS工具包）。

1.前期准备：数据加载与预处理

（1）MODIS数据拼接与格式转换

MODISL1B数据为HDF格式，需先转换为ENVI可识别的格式：

打开ENVI→工具箱→?RasterManagement→ModisConversionTools→ModisHDFtoENVI；

选择下载的MOD021KM（辐射亮度数据）和MOD03（地理定位数据），勾选需要的波段（重点选热红外波段31、32，以及可见光波段用于NDVI计算）；

输出路径设置后，点击OK，生成ENVI格式的辐射亮度影像（含地理信息）。

（2）辐射定标（自动完成）

上述转换过程中，ENVI会自动读取MODIS头文件中的gain/offset系数，将DN值转换为表观辐射亮度?\(L_{\lambda}\)，无需手动计算；若需验证，可在BasicTools→BandMath中输入公式：b1*gain+offset（gain/offset从MODIS元数据中提取）。

（3）云检测与掩膜

方法1：利用MOD35云产品（HDF格式），同样通过ModisHDFtoENVI转换为ENVI格式，将云像元（值为1/2）设为掩膜；

方法2：阈值法手动掩膜：工具箱→?BandMath，输入NDVI阈值（如NDVI0.1且热红外波段31亮温300K判定为云/裸土干扰），生成掩膜文件；再通过BasicTools→Masking→ApplyMask，将掩膜应用到辐射亮度影像。

2.大气校正：劈窗算法（简化版，无需大气廓线）

ENVI无内置MODTRAN/6S直接调用接口，推荐用劈窗算法工具（需加载ENVI扩展工具包，如ENVI_Split_Window_Tool）：

下载并安装MODIS劈窗算法扩展工具（可从ESRI官网或开源社区获取）；

打开工具，输入以下参数：

参数项

取值/说明

波段31辐射亮度影像

预处理后的MODISBand31数据

波段32辐射亮度影像

预处理后的MODISBand32数据

大气水汽含量

可从MOD07产品提取（或经验值：热带≈5g/cm2，温带≈2g/cm2）

观测天顶角

从MOD03数据中提取（影像头文件可查看）

点击运行，工具自动计算波段31/32的亮温，并代入劈窗算法公式（如前文\(T_s=a_0+a_1T_{31}+a_2T_{32}+a_3(T_{31}-T_{32})\)），输出地表温度影像（单位：K，可手动转换为℃：T(℃)=T(K)-273.15）。

3.比辐射率计算（植被-裸土混合像元）

第一步：计算NDVI工具箱→?Vegetation→NDVI，选择MODIS红波段（Band1）和近红外波段（Band2），输出NDVI影像；

第二步：计算植被覆盖度\(F_v\)BandMath输入公式：(b1-NDVI_min)/(NDVI_max-NDVI_min)（NDVI_min=0.05，NDVI_max=0.7，经验值），生成\(F_v\)影像；

第三步：线性混合模型计算比辐射率BandMath输入公式：b1*0.986+(1-b1)*0.975+0.001（b1为\(F_v\)影像，0.986=植被比辐射率，0.975=裸土比辐射率，0.001为修正项\(d\varepsilon\)），输出波段31的比辐射率影像。

4.温度反演（普朗克定律迭代）

若不用劈窗算法，需手动基于普朗克定律计算：

BandMath中输入普朗克定律反演公式（迭代求解需自定义函数，ENVI中可通过ENVIModeler实现）：核心公式（简化版，已知\(L_{s,\lambda}\)和\(\varepsilon_{\lambda}\)）：\(T_s=c_2/(\lambda\cdotln(c_1/(\lambda^5\cdotL_{s,\lambda}/\varepsilon_{\lambda})+1))\