基于MODIS影像下的温度反演研究

摘要： 利用遥感影像获取地表温度信息，对于检测地表状况和生态环境变化有重要意义。在对地观测卫星MODIS的16个热红外波段中，其中有6个波段是针对地表温度来设计的，其中第31（10.78～11.28um）和32（11.77～12.27um）与TM的第6波段区间基本相同。运用单窗算法分别对这2个波段建立方程，然后通过对Planck方程的辐射率和温度计算模拟，发现温度和辐射率近似线性相关，本文通过对MODIS遥感影像进行数据读取、辐射定标、反演温度、区域裁剪预处理，得到了研究区地表温度的空间分布变化，最后运用单窗算法对31通道的影像反演了研究区的温度分布。关键词： MODIS；地表温度；单窗算法

中图分类号：TP751 文献标识码：A 文章编号：1006-4311（2014）24-0224-02

0 引言

地表温度[1-2]是研究区域地表能量平衡和资源环境变化的重要参数之一，它在众多领域有着广泛的需求，因此准确获取其分布的时空状况具有非常重要的意义。要想获得大区域的地表温度，就需要通过遥感手段，因为温度的遥感反演是遥感技术的关键环节之一，具有非常重要的现实意义。

MODIS（中分辨率成像光谱仪）[3]是新一代地球观测传感器，搭载于NASA的地球观测系统（EOS）TERRA和AQUA两颗卫星之上，其数据在全世界范围都可免费接收。MODIS数据采用HDF格式进行存储、。HDF的数据结构是一种分层式数据管理结构，存储了MODIS数据以及相应元数据，提供了MODIS数据辐射校正所需的相应卫星参数。

1 研究区概况

会泽县位于云南省东北部、金沙江东岸、曲靖市西北部，全县占地面积6077平方公里，地处东经103°03′～ 103°55′、北纬25°48′～27°04′之间。而县城所在地金钟镇，海拔达到2120米，距离曲靖市区有187公里，距省会昆明205公里。

2 数据来源

本文采用2010年1月8日12时25分（本地时间）Terra（AM）卫星MODIS数据，影像覆盖范围为北纬26°～28°N，东经102°～103°E。选取云南省会泽县为研究区域，没有云层覆盖，无雪覆盖，质量良好。本次试验选取MODIS的第31波段（空间分辨率为1km）为试验对象。

3 数据处理

3.1 单窗算法反演温度 单窗算法是覃志豪[4-6]针对Landsat TM 影像在只有一个热红外波段TM6的时候提出的，该算法最大的优点就是大气参数少，但是精度高。因此，在对MODIS数据利用ENVI自带的GEOREFERENCE MODIS模块作几何校正，去除蝴蝶结效应的基础上，运用单窗算法对MODIS的31通道影像反演了会泽地区的温度分布。

如果大气透射率τ6，大气平均作用温度Ta和地表辐射率ε6已知，则可用如下单窗算法从像元的亮点温度值T6中推算该像元的实际地表温度TS：

3.2 亮度温度的计算

我们可以使用Planck方程求解出星上亮度温度。对31波段应用Planck方程，并化简可得31波段的亮温计算公式：

本文用到的MODIS影像是摄取于1月上旬，对会泽地区而言，通过利用上述方法对31波段进行地面温度反演发现，31波段的星上亮度温度非常接近。通过用单通道方法反演，第31波段的最高温度为22.85℃，平均温度为8.15℃。符合影像摄取时的当地实际温度范围，这种大面积温度分布规律在TM影像上呈现不出来，但是反映大区域的温度分布差异是MODIS影像的一个优势。如图2。

如图3，散点图中可以清晰的看出研究区的温度区间为-8℃-22℃之间，为我们研究该区域给出了直观的分析。

由图4会泽区空间温度分布图看出，温度的地域分布很有特点，总体的特征是随着维度的地表温度增大的，而温度大致还可分为3个区间：11℃。这种地温产生的的原因主要是由一些云块引起的异常，是因为较浓厚的云层遮蔽，使得传感器接受到地面极低的热辐射量，从而反演出的温度值极低。

4 结论

本文利用基于单窗算法的地表温度反演算法，通过采用MODIS的第31波段提取地表比辐射率，结合热红外波段的亮温信息，反演了会泽地区的地表温度。结果表明：反演精度比较高。同时对MODIS的第31通道的参数进行了计算，对会泽地区的整个MODIS用该方法进行了反演，并对会泽地区反演的地表温度进行了分析。结果表明反演的地表温度比较合理。该方法适宜于该地区。

参考文献：

[1]Becker，F，Li，Z，L..Towards a local split window method over land surface[J]. Intermati-onalJourmal of Remote Sensing，1990（11）：369-393.

[2]LijuanZhang，XingpingWen，JunWang，etal. Applied Materiala and Manufacturing[C]. Trans Tech Publications， Vols.2014，889-890：1634-1637.

[3]覃志豪.Li Wenjuan.张名华等.单窗算法的大气参数估计方法[J].国土资源遥感，2003（56）.