envi反演水质参数_ENVI5.2中基于MODIS数据的海表温度反演

劈窗算法最初是为反演海面温度开发的，具体地说是针对NOAA/AVHRR的4和5通道设计的，后来也被用来反演地表温度，这种算法较成熟，精度也高。劈窗算法以地表热辐射传导方程为基础，利用10~13μm

大气窗口内，两个相邻热红外通道(一般为10.5～11.5μm、11.5～12.5μm)对大气吸收作用的不同，通过两个通道测量值的各种组合来剔除大气的影响，进行大气和地表比辐射率的修正。表达式为：

T S= T 4+ A (T

4- T 5) + B

其中：T S为地表真实温度，T 4和 T

5分别为AVHRR的4和5通道，A和B为常量。

AVHRR 的通道4(10.15~11.13μm) 和通道5 (11.15~12.15μm)

恰与MODIS 的第31 波段 (10.178~ 11.128μm) 和32 波段 (11.177~ 12.127μm )

的中心波长相对应, 可将MODIS 的31、32 波段资料,

用于劈窗算法进行地表温度计算。很多学者对这个算法进行了推演，得到很多新的算法，如覃志豪、毛克彪等人。本文就是使用其他学者推演的算法。

利用MODIS数据劈窗算法反演海表温度技术流程如下图：

图：技术流程图

注：(1)按照本流程反演出来的结果是SST。陆地上的值可以视为无效值，若要得到正确的陆表温度，需要加入海陆分离的步骤，以及城镇和自然表面的比辐射率计算。

下面详细介绍处理流程。

第一步：读取原始数据

打开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这两项参数代