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

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阈值（如NDVI<0.1且热红外波段31亮温>300K判定为云/裸土干扰），生成掩膜文件；再通过BasicTools→Masking→ApplyMask，将掩膜应用到辐射亮度影像。2.大气校正：劈窗算法（简化版，无需大气廓线）ENVI无内置MODTRAN/6S直接调用接口，推荐用劈窗算法工具（需加载ENVI扩展工具包，如ENVI_Split_Window_Tool）：下载并安装MODIS劈窗算法扩展工具（可从ESRI官网或开源社区获取）；打开工具，输入以下参数：参数项取值/说明波段31辐射亮度影像预处理后的MODISBand31数据波段32辐射亮度影像预处理后的MODISBand32数据大气水汽含量可从MOD07产品提取（或经验值：热带≈5g/cm²，温带≈2g/cm²）观测天顶角从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))\)（\(\lambda\)=11.03μm，Band31中心波长；\(c_1=3.7418e-16\)，\(c_2=1.4388e-2\)）；运行后生成地表温度影像，可通过ColorMapping调整色标，直观查看温度分布。5.结果验证与导出验证：导入实测温度点数据（ENVI→

File→OpenAs→Vector→PointFile），用BasicTools→Statistics→ComputeStatistics计算反演值与实测值的RMSE；导出：File→Export→ExporttoGeoTIFF，将温度影像转换为通用格式，用于ArcGIS等软件后续分析。二、IDL编程实现（自动化，适合批量处理）IDL（InteractiveDataLanguage）是MODIS数据处理的经典编程工具，可编写脚本实现批量反演，核心代码框架如下（聚焦劈窗算法，可直接修改参数运行）：1.核心代码框架（MODIS劈窗算法反演LST）idl;=====================1.读取MODIS数据=====================;读取MOD021KMHDF数据（辐射亮度）file_mod02='MOD021KM.A2025001.0000.061.2025002020558.hdf';替换为你的数据路径hdf_id=HDF_OPEN(file_mod02,/READ);提取Band31和Band32的辐射亮度（DN值）data31=HDFSD_GETDATA(hdf_id,'EV_1KM_Emissive',0,10);Band31对应索引10data32=HDFSD_GETDATA(hdf_id,'EV_1KM_Emissive',0,11);Band32对应索引11;提取定标系数（gain/offset）gain31=HDFATTR_GET(hdf_id,'EV_1KM_Emissive','gain',10)offset31=HDFATTR_GET(hdf_id,'EV_1KM_Emissive','offset',10)gain32=HDFATTR_GET(hdf_id,'EV_1KM_Emissive','gain',11)offset32=HDFATTR_GET(hdf_id,'EV_1KM_Emissive','offset',11)HDF_CLOSE,hdf_id;=====================2.辐射定标：DN→辐射亮度=====================L31=data31*gain31+offset31;Band31表观辐射亮度(W·m^-2·sr^-1·μm^-1)L32=data32*gain32+offset32;Band32表观辐射亮度;=====================3.计算星上亮温T31/T32（普朗克定律反推）=====================c1=3.7418e-16;普朗克常数1c2=1.4388e-2;普朗克常数2lambda31=11.03e-6;Band31中心波长(m)lambda32=11.95e-6;Band32中心波长(m);亮温计算公式（黑体辐射反推）T31=c2/(lambda31*ALOG(c1/(lambda31^5*L31)+1))T32=c2/(lambda32*ALOG(c1/(lambda32^5*L32)+1));=====================4.劈窗算法计算地表温度=====================;经验系数（适配中纬度、大气水汽含量2g/cm²，可根据区域调整）a0=0.0a1=1.019a2=-0.019a3=1.35Ts=a0+a1*T31+a2*T32+a3*(T31-T32);地表温度（K）Ts_C=Ts-273.15;转换为摄氏度;=====================5.掩膜异常值+输出结果=====================;掩膜云/雪（Ts_C>50或Ts_C<-20判定为异常）mask=(Ts_CGT50)OR(Ts_CLT-20)Ts_C[mask]=!VALUES.F_NAN;异常值设为NaN;保存为ENVI格式ENVI_WRITE_ENVI_FILE,Ts_C,$FILENAME='MODIS_LST_2025001.dat',$INTERLEAVE=0,$;BSQ格式MAP_INFO=ENVI_GET_MAP_INFO(FILE=file_mod02),$;继承地理信息DESCRIPTION='MODIS地表温度（℃）_劈窗算法'PRINT,'反演完成！地表温度范围：',MIN(Ts_C,/NAN),'~',MAX(Ts_C,/NAN),'℃'2.关键代码说明数据读取：HDFSD_GETDATA函数用于提取HDF文件中的波段数据，MODIS热红外波段索引可参考MODIS官方文档（Band31=10，Band32=11）；亮温计算：通过普朗克定律反推星上亮温（传感器接收到的黑体温度），是劈窗算法的基础；经验系数调整：a0-a3需根据研究区大气条件调整（如热带地区a3取1.5，高纬度取1.2），可参考MODIS劈窗算法相关文献；批量处理：可嵌套FILE_SEARCH函数遍历文件夹下所有MODIS数据，实现自动化批量反演。三、工具选择与注意事项工具优势适用场景注意事项ENVI可视化操作，无需编程单景数据处理、新手入门、结果可视化