一种基于MODIS与Landsat数据的改进型地表温度日较差降尺度估算方法与流程

本发明涉及干旱、土壤水分遥感监测领域，具体涉及一种基于modis与landsat数据的改进型地表温度日较差降尺度估算方法。

背景技术：

目前遥感监测干旱和土壤水分是热点问题。不管是监测干旱还是土壤水分，对于大众服务来讲，均需要被告知地表最高温度和地表最低温度分别是多少。现有获取地表温度日较差的方法主要是依据modis数据，获取空间分辨率为1000m地表温度日较差，空间分辨率低，掩盖了地表温度日较差空间差异信息，仅适用于大尺度，地形简单，空间异质性小的区域。

技术实现要素：

为解决上述问题，本发明提供了一种基于modis与landsat数据的改进型地表温度日较差降尺度估算方法。

为实现上述目的，本发明采取的技术方案为：

一种基于modis与landsat数据的改进型地表温度日较差降尺度估算方法，基于土壤热传导理论建立土壤温度日较差模型，以landsat卫星数据和modismod11a2地表温度产品为数据源，获取空间分辨率为30m的地表温度日较差数据；具体包括如下计算步骤：

具体包括如下步骤

令z=0，=，可以如下式推导：

由：；

；

可得：

由：；

；

；

可得：；

其中，为土壤深度为zcm，时间为t时的地表温度，为地表温度变化周期函数，为通过日最低温（tmin）估算的地温日较差，通过日最高温（tmax）估算的地温日较差，为最终估算的地表温度日较差，ttm为由tm或tirs传感器数据计算得到的地表温度数据单位℃；tmax为modis数据获得的白天地表温度单位℃；tmin为modis数据获得的夜晚地表温度单位℃，t为地表温度日较差的待估算时刻，即landsat卫星过境时间；tmin为modis卫星夜晚过境时刻；tmax为modis卫星白天过境时刻；p是一个周期,通常为24h，p1为2倍的modis夜晚过境时刻与白天过境时刻的差值，单位：h，即2✖（tmax-tmin），w1、w2为时间权重；α是土壤热扩散率，单位cm²/h。

本发明具有以下有益效果：

将地表温度日较差数据空间分辨率由1000m提高至30m，可以获取更精细的地表温度日较差纹理信息。

附图说明

图1是由modis数据获得的空间分辨率为1000m地表温度日较差数据图2为基于modis与landsat数据的改进型地表温度日较差降尺度获得的空间分辨率为30m地表温度日较差数据。

具体实施方式

下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。

土壤温度日变化曲线与太阳辐射日变化曲线相似，土壤的热量是由于地面接受了太阳辐射增温后，向下传导给土壤，土壤才得以增温，因此温度和辐射二者有直接相关关系。随离地表的距离向下的增加，温度变化振幅越小，向下传输过程中，每层土壤都要留下一部分热量，所以越往下获得的热量越小，增温幅度越小，同时热量的传导需要一定的时间，每层达到最高或最低温的时间也越滞后。基于上述的土壤热传导理论建立土壤日温差模型，以landsat卫星数据和modismod11a2地表温度产品为数据源，获取空间分辨率为30m的地表温度日较差数据。

具体的，通过以下公式获取日较差数据：

；

其中，为最终估算的地表温度日较差，ttm为由tm或tirs传感器数据计算得到的地表温度数据单位℃；tmax为modis数据获得的白天地表温度单位℃；tmin为modis数据获得的夜晚地表温度单位℃，t为地表温度日较差的待估算时刻，即landsat卫星过境时间；tmin为modis卫星夜晚过境时刻；tmax为modis卫星白天过境时刻；p是一个周期,通常为24h，p1为2倍的modis夜晚过境时刻与白天过境时刻的差值，单位：h，即2✖（tmax-tmin），w1、w2为时间权重；α是土壤热扩散率，单位cm²/h。

由图1、图2可以看出基于modis与landsat数据的改进型地表温度日较差降尺度获得的地表温度日较差数据空间分辨率较高，并且能够更精细的反映出地表温度日较差的纹理信息。

本发明使用modis与landsat数据可以将地表温度日较差空间分辨率由1000m提升至30m，30m地表温度日较差可以反映出更精细的地表温度日较差纹理信息，使其在不受区域尺度和地形限制的同时，扩大了其适用性。在使用landsat数据获取的地表温度大于modis数据白天地表温度时，使用该方法降尺度得到的地表温度日较差不准确，此情况不建议使用该方法。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。