留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

利用Landsat TM影像进行地表温度像元分解

王斐 覃志豪 宋彩英

王斐, 覃志豪, 宋彩英. 利用Landsat TM影像进行地表温度像元分解[J]. 武汉大学学报 ● 信息科学版, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
引用本文: 王斐, 覃志豪, 宋彩英. 利用Landsat TM影像进行地表温度像元分解[J]. 武汉大学学报 ● 信息科学版, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
WANG Fei, QIN Zhihao, SONG Caiying. An Efficient Approach for Pixel Decomposition of Land Surface Temperature from Landsat TM Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
Citation: WANG Fei, QIN Zhihao, SONG Caiying. An Efficient Approach for Pixel Decomposition of Land Surface Temperature from Landsat TM Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604

利用Landsat TM影像进行地表温度像元分解

doi: 10.13203/j.whugis20140604
基金项目: 

国家自然科学基金 41471300

详细信息
    作者简介:

    王斐, 博士生, 主要从事热红外遥感与农业遥感研究。essiwf@163.com

    通讯作者: 覃志豪, 博士, 教授。zhihaoqin@163.com
  • 中图分类号: P237.9

An Efficient Approach for Pixel Decomposition of Land Surface Temperature from Landsat TM Data

Funds: 

The National Natural Science Foundation of China 41471300

More Information
    Author Bio:

    WANG Fei, PhD candidate, specializes in thermal remote sensing and agriculture remote sensing. E-mail: essiwf@163.com

    Corresponding author: QIN Zhihao, PhD, professor. E-mail: zhihaoqin@163.com
图(7)
计量
  • 文章访问数:  666
  • HTML全文浏览量:  9
  • PDF下载量:  325
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-12-28
  • 刊出日期:  2017-01-05

利用Landsat TM影像进行地表温度像元分解

doi: 10.13203/j.whugis20140604
    基金项目:

    国家自然科学基金 41471300

    作者简介:

    王斐, 博士生, 主要从事热红外遥感与农业遥感研究。essiwf@163.com

    通讯作者: 覃志豪, 博士, 教授。zhihaoqin@163.com
  • 中图分类号: P237.9

摘要: 提出了一种基于Landsat TM的地表温度二次像元分解方法,将地表温度的空间分辨率从120 m提高到30 m。首先,利用地表类型的线性统计模型(E-DisTrad)获取初次分解子像元的地表温度,计算得到初次分解子像元的辐亮度;然后,利用面向对象的图像分割方法获取二次分解子像元的权重,实现对地表温度的二次分解;最后,采用升尺度再分解的验证方法进行精度分析,并选取了北京市TM影像进行实例分析。实验结果表明,二次像元分解模型不仅能有效地提高地表温度的空间分辨率,反映出不同地表类型地表温度的空间差异性,而且保证了像元分解前后能量值的一致性,非常适合于复杂地表覆盖地区的热红外波段遥感影像数据的降尺度处理。

English Abstract

王斐, 覃志豪, 宋彩英. 利用Landsat TM影像进行地表温度像元分解[J]. 武汉大学学报 ● 信息科学版, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
引用本文: 王斐, 覃志豪, 宋彩英. 利用Landsat TM影像进行地表温度像元分解[J]. 武汉大学学报 ● 信息科学版, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
WANG Fei, QIN Zhihao, SONG Caiying. An Efficient Approach for Pixel Decomposition of Land Surface Temperature from Landsat TM Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
Citation: WANG Fei, QIN Zhihao, SONG Caiying. An Efficient Approach for Pixel Decomposition of Land Surface Temperature from Landsat TM Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(1): 116-122. doi: 10.13203/j.whugis20140604
  • 高空间分辨率的地表温度(land surface temperature,LST)数据对于城市热岛监测、土壤水分估算、地气水热交换等有着重要的作用[1-3]。目前主要搭载热红外波段传感器的遥感平台中,可见光/近红外波段空间分辨率都比热红外波段空间分辨率高[4],使得从热红外波段反演得到的地表温度空间分辨率也不高。因此,如何综合可见光/近红外波段的信息提高地表温度的空间分辨率一直是热红外遥感应用研究的重要方向[5-6]

    目前常用的分解方法是文献[7]提出的Distrad算法,该算法通过拟合地表温度与归一化植被指数(normalized difference vegetation index, NDVI)的线性统计关系LST-NDVI实现地表温度的降尺度。文献[8]基于LST-NDVI统计模型中干、湿边组成的三角形关系,提出了热红外波段像元分解的三角算法。TsHARP方法[9]提出用植被覆盖度来代替NDVI。但上述模型仅考虑了单一的植被指数,因此,其他不同的遥感指数也逐渐被引入。文献[10]分析了地表温度与归一化裸土指数的统计关系。文献[11]根据不同的地表类型,分别统计了地表温度与15种遥感指数之间的线性关系,提出了E-DisTrad算法,并分析得出在城市地表中不透水层覆盖度与地表温度的拟合关系最佳。上述模型利用地表温度与遥感指数的统计关系实现了地表温度像元分解,但是整个分解过程中没有考虑像元分解前后辐射能量值的一致性。

    本文以北京市Landsat TM影像为数据源,首先基于地表类型的线性统计模型(E-DisTrad)获得高空间分辨率的初次分解地表温度,计算得到初次分解子像元的辐亮度。然后,采用面向对象的图像分割方法对初次分解子像元辐亮度进行4×4规则分割,得到初次分解子像元的辐亮度与其分割后每个对象的辐亮度之间的权重。最后,根据上述权重实现对地表温度父像元的二次分解。由于该方法涉及地表温度两次分解,因此,本文将此方法称之为地表温度二次像元分解方法。

    • 本文使用Landsat TM数据,成像地区位于北京,过境时间是2009-07-20T10:42:26。对TM数据进行辐射定标、大气校正等预处理后,选取了北京市一个矩形区域作为本文的实验区,实验区TM假彩色合成影像(RGB432)如图 1所示。

      图  1  实验区Landsat TM假彩色合成图

      Figure 1.  TM False-Color Composite Image of Experimental Area

    • Landsat TM可见光/近红外波段的空间分辨率为30 m,热红外波段的空间分辨率为120 m,反演得到的地表温度的空间分辨率也是120 m,从而极大地限制了地表温度的应用。本文的研究目标是以高空间分辨的可见光/近红外波段为辅助数据,提出一种新的地表温度像元分解算法,以期将地表温度的空间分辨率提高到30 m。本文提出的二次像元分解方法的原理如图 2所示。

      图  2  二次像元分解方法原理图

      Figure 2.  Double-Step Pixel Decomposition Schematics

      利用120 m空间分辨率的地表温度数据Ts,对于任一个父像元,可以计算得到其分解前的辐亮度值(I)。以高空间分辨率的可见光/近红外波段数据为基础,以每个分解后子像元占其父像元辐射能量值的权重(Wi)为辅助数据,从而可以确定每个分解后子像元的辐亮度值(Ii),获得16个30 m空间分辨率的子像元(Tsi)地表温度。因此,如何获取子像元的权重是本文方法的重点。

      本文首先利用E-DisTrad方法分别统计地表温度与不同地表类型的遥感指数之间的线性关系,从而获得高空间分辨率的初次分解地表温度,计算得到初次分解子像元的辐亮度值,为二次分解子像元的权重的计算提供基础。因此,如何获取初次分解地表温度是本文方法的第一步。

    • 初次分解地表温度的计算是通过E-DisTrad方法进行的。水体的地表温度趋于稳定,故可以在分解时赋予一个合理的统计值,对植被和建筑用地两类地表类型,分别统计NDVI和不透水层(impervious surface area, ISA)覆盖度与地表温度的线性统计关系。

      利用可见光/近红外波段计算得到30 m空间分辨率的归一化植被指数NDVI30和不透水层指数ISA30。由于从热红外波段反演得到的地表温度空间分辨率为120 m,因此,需要将NDVI30和ISA30重采样到120 m,从而分别统计地表温度LST120与归一化植被指数NDVI120和不透水层指数ISA120的线性关系,得:

      $$ \begin{array}{l} {\rm{植被: }}{\rm{LS}}{{\rm{T}}_{120}} = a{\rm{NDV}}{{\rm{I}}_{120}} + b\\ 建筑用地:{\rm{LS}}{{\rm{T}}_{120}} = c{\rm{IS}}{{\rm{A}}_{120}} + d \end{array} $$ (1)

      式中,abcd分别为地表温度与遥感指数之间的线性拟合参数。LST120采用单窗算法反演,ISA120采用线性光谱分解模型计算,NDVI120采用TM波段运算获得,上述参数的计算详见参考文献[12-14]。

      E-DisTrad方法假设地表温度与遥感指数的拟合关系在不同空间分辨率上是一致的,因此,可以将上述线性关系直接应用于NDVI30和ISA30遥感指数,从而获得30 m空间分辨率的初次分解地表温度LST30,得:

      $$ \begin{array}{l} 植被:{\rm{LS}}{{\rm{T}}_{30}} = a{\rm{NDV}}{{\rm{I}}_{30}} + b\\ 建筑用地:{\rm{LS}}{{\rm{T}}_{30}} = c{\rm{IS}}{{\rm{A}}_{30}} + d \end{array} $$ (2)

      通过上述分解方法获取的初次分解地表温度并不能保证能量守恒,因此,以初次分解地表温度为基础,获取每个二次分解子像元占其父像元辐亮度的权重,从而得到每个二次分解子像元的辐亮度值,实现对地表温度的二次分解。

    • 根据120 m空间分辨率的地表温度,对于任一个父像元,根据Stenfan-Boltzman定律,可以计算得到辐亮度值I

      $$ I = \varepsilon \sigma {T^4} $$ (3)

      式中,I为父像元的辐射亮度值(W·m-2);ε为父像元的比辐射率;T为父像元的地表温度(K);σ为Stenfan-Boltzman常数,σ=5.67×10-8W·m-2·K-4。二次分解子像元权重计算方法见图 3

      图  3  权重计算方法图

      Figure 3.  Diagram of Calculation of Weight

      图 3所示,对于任一个30 m空间分辨率的初次分解子像元地表温度计算得到辐亮度值:

      $$ {I_k} = {\varepsilon _k}\sigma T_k^4 $$ (4)

      式中,Ik为第k个初次分解子像元的辐亮度值(W·m-2);εk为第k个初次分解子像元的比辐射率;Tk为第k个初次分解子像元的地表温度(K)。

      然后,根据面向对象的图像分割方法,对初次分解子像元的辐亮度进行4×4规则分割,使分割得到的每个对象与其分解前父像元是相对应的。则每个分割后对象的辐亮度值为:

      $$ I\prime = \sum\limits_{k = 1}^{16} {{I_k}{S_2}/{S_1}} $$ (5)

      式中,I′为分割后对象的辐亮度值(W·m-2);S1S2分别为分割后对象和初次分解子像元的面积(m2)。从而可以确定初次分解子像元的辐射能量与其分割后每个对象的辐射能量之间的权重,得到每个二次分解子像元的权重Wk

      $$ {W_k} = \frac{{{I_k}{S_2}}}{{I\prime {S_1}}} = {I_k}/\sum\limits_{k = 1}^{16} {{I_k}} $$ (6)

      最后,利用上述权重值实现对其分解前120 m空间分辨率的父像元的二次分解,得到每个二次分解子像元的辐亮度值(Ik):

      $$ {{I'}_k} = {W_k}I{S_1}/{S_2} = 16{W_k}I $$ (7)

      因此,二次分解子像元的地表温温度Tsk为:

      $$ {T_{{s_k}}} = {({{I'}_k}/\sigma {\varepsilon _k})^{0.25}} $$ (8)
    • 由于地表温度每天变化较快,很难找到实时的温度场来验证TM影像的分解结果,故本文借鉴文献[15]提出的“先升尺度再分解”验证方法,评价方法的优劣性。

      首先,利用单窗算法反演得到120 m空间分辨率的地表温度,可以认为该地表温度能够反映出实验区真实的地表温度空间分布,即为“真值”。然后,采用均值法将上述地表温度重采样到480 m。最后,分别采用E-DisTrad方法和二次像元分解方法将480 m空间分辨率的地表温度数据降尺度到120 m,与“真值”进行比较分析,评价两种分解方法的差异性。

    • 利用单窗算法反演得到的地表温度图像的空间分辨率是120 m,可以认为该图像代表真实的地表温度空间分布 Ts (见图 4(a))。图 4(b)4(c)分别表示E-DisTrad和二次像元分解两种地表温度图像分解模型的结果TeTw

      图  4  地表温度分解结果空间分布图

      Figure 4.  LST Spatial Distribution of Different Methods

      图 4可知,实验区内建筑用地的地表温度最高,植被次之,水体的地表温度最低,3种算法计算得到的地表温度具有相似的整体空间格局。高温区(≥310 K)主要集中于北京市城区的南部,次高温区(305~310 K)主要集中于中部,低温区(≤305 K)主要集中于水体和植被覆盖度较高的郊区。这与北京中心城区快速的城市化进程相吻合。与Ts(图 4(a))相比,TeTw(图 4(b)4(c))地表温度的空间分辨率都为30 m,相较于热红外波段的空间分辨率(120 m)提高了3倍。对于E-DisTrad模型,高温区和低温区的像元明显增多,地表温度趋于305~307 K的像元几乎没有,与分解前实验区地表温度的空间分布一致性较差。对于二次像元分解模型,地表温度的整体变化趋势更为平缓,但是,地表温度较低(302~305 K)的像元增多,多发现于郊区植被中。

    • 为了验证二次像元分解方法的合理性,本文采用“先升尺度再分解”的方法获得高空间分辨率的地表温度数据(见图 5)。其中,图 5(a)是TM6波段采用单窗算法反演得到的120 m空间分辨率的地表温度(Ts),作为两种方法对比的真值;图 5(b)5(c)分别是将120 m的热红外波段数据升尺度到480 m后,经过E-DisTrad方法和二次像元分解方法得到的120 m空间分辨率的地表温度结果(TeTw)。

      图  5  地表温度分解结果验证图

      Figure 5.  Validation of LST spatial Distribution for Different Methods

      Ts相比较,Te地表温度的统计值明显偏高,中心区域的地表温度普遍高于320 K,多被发现在植被和城市用地的混合地区,地表温度趋于305~310 K的像元几乎没有。Tw地表温度的统计值与观测值较为接近,地表温度的空间分布与观测结果更为一致,地表温度的变化趋势更为平缓。这与地表温度的频数统计曲线(见图 6(b))相一致。

      图  6  遥感指数与地表温度像元频数统计对比图

      Figure 6.  Frequent Comparison for Remote Sensing Indices and LST

      图 6(b)中可以看出,二次像元分解方法结果Tw(紫色)的地表温度频数统计曲线与单窗算法反演结果Ts(红色)更为接近,不同地表温度区间的像元数目也较为接近。与单窗算法反演结果相比,在302~310 K地表温度区间内,E-DisTrad方法降尺度得到的地表温度像元数目明显少于TM6反演结果,而在314 K和300 K处(图 6(b)AB)地表温度像元频数出现两个峰值,这是由于表征建筑用地的ISA遥感指数在0.7~0.9区间内的像元频数最大(图 6(a)A),从而导致城市地区地表温度初次分解结果多集中在314 K左右,高估了城市地表的温度。与此相同的是,表征自然表面的NDVI遥感指数在0.6~0.8区间内的像元频数也比较多(图 6(a)B),使得初次分解温度在300 K处出现另一个峰值。这两个峰值的出现,使得302~310 K地表温度区间内的像元数目远远低于实际真值。但是,通过二次像元分解方法分解后,上述情况得到明显改善,二次像元分解方法得到的不同地表温度区间内像元频数与单窗算法反演结果更为接近。需要特别指出的是,在312~314 K地表温度区间内,二次像元分解方法结果的像元数目增多,多发现于建筑用地和林草地等的混合区域。图 7给出了TeTwTs之间地表温度的散点图。从图 7中可以看出,TwTs之间地表温度的散点图拟合线的斜率更接近1,R2值更大,更加印证了基于二次像元分解模型降尺度得到的地表温度与观测结果更为接近。

      图  7  不同分解方法所得结果与TM6反演结果的对比散点图

      Figure 7.  Comparison of Scatters Between Different Decomposed Results and LST Retrieved from TM6

    • 本文提出了一种能量平衡角度的地表温度二次像元分解方法,并将该方法应用到具有复杂地表覆盖类型的北京市Landsat TM遥感影像中,实现了地表温度的分解,将空间分辨率提高到30 m,相较于热红外波段的空间分辨率(120 m)提高了3倍。

      同TM6波段直接利用单窗算法反演得到的地表温度相比,本文采用升尺度再降尺度的验证方法,讨论了E-DisTrad模型和二次像元分解模型分解后获得的地表温度的精度。分解后地表温度的空间分布、频数统计曲线和散点图分析表明,利用二次像元分解模型降尺度得到的地表温度不仅能准确地反映出实验区不同地表类型地表温度空间变化的差异性,而且有效地保证了分解前后影像的总体及局部能量的守恒,非常适合用于复杂地表覆盖地区的热红外波段影像数据的降尺度处理。

      需要指出的是,本文中初次分解地表温度的计算以E-DisTrad模型为基础,其他不同的地表温度像元分解方法也可以进行引入实验;子像元权重计算过程中采用面向对象的4×4规则分割,但是对象的定义不单单是规则的矩形,其他分割方法仍需进行进一步的研究。

参考文献 (15)

目录

    /

    返回文章
    返回