劈窗算法最初是为反演海面温度开发的,具体地说是针对NOAA/AVHRR的4和5通道设计的,后来也被用来反演地表温度,这种算法较成熟,精度也高。劈窗算法以地表热辐射传导方程为基础,利用10~13μm 大气窗口内,两个相邻热红外通道(一般为10.5~11.5μm、11.5~12.5μm)对大气吸收作用的不同,通过两个通道测量值的各种组合来剔除大气的影响,进行大气和地表比辐射率的修正。
表达式为:TS= T4+ A(T4- T5) + B
其中:TS为地表真实温度,T4和 T5分别为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 技术流程图
注:(1)按照本流程反演出来的结果是SST。陆地上的值可以视为无效值,若要得到正确的陆表温度,需要加入海陆分离的步骤,以及城镇和自然表面的比辐射率计算。
(2)MODIS数据下载:nasa官网:http://modis.gsfc.nasa.gov/.
http://ladsweb.nascom.nasa.gov/data/search.html
下面详细介绍处理流程。
第一步:数据打开
Open->Open as->EOS->MODIS,选择MOD021KM.A2013185.0245.005.2013185094144.hdf文件,点击OK打开。分为三个数据集:热红外数据(20-36波段),可见光到短波红外的辐射率数据(1-19、26波段),可见光到短波红外的反射率数据