一种基于卫星遥感数据的植被温度和裸地温度估算方法技术

本发明专利技术公开了一种基于卫星遥感数据的植被温度和裸地温度估算方法，包括以下步骤：步骤一、准备卫星遥感数据和地面气象观测数据；步骤二、根据步骤一的数据计算植被覆盖度；步骤三、根据步骤一的数据分别计算裸地反照率、裸地发射率和植被发射率；步骤四、计算裸地温度和植被温度。本发明专利技术以地表能量平衡原理为基础，不依赖于其他像元信息，应用目标像元的卫星遥感数据，包括但不限于地表温度、地表发射率、植被指数、植被类型和地面气象观测数据为输入，计算出目标像元的植被温度和裸地温度，具有空间独立性，提高了估算精度。

【技术实现步骤摘要】
一种基于卫星遥感数据的植被温度和裸地温度估算方法
本专利技术涉及一种植被温度和裸地温度估算方法，属于卫星遥感
，尤其涉及一种基于卫星遥感数据的植被温度和裸地温度估算方法。
技术介绍
卫星遥感反演的地表温度(LandSurfaceTemperature(LST)或TL)是以像元尺度表征的，仅代表像元所覆盖区域的平均陆地表面温度，无法区分不同植被覆盖度的像元内的植被温度和裸地温度。地表温度(LandSurfaceTemperature(LST)或TL)表征了地球陆地表面的温度状态，是地表能量平衡和水循环中的关键变量，在地-气相互作用和气候变化研究中扮演着十分重要的角色，对气象、水文、农业、生态和环境等领域均有着重要的应用价值。卫星遥感反演的地表温度由于具有固定的分辨率，如NASAAVHRR-LST约为1km×1km，MODISLST约为500m×500m，Landsat-LST约为30m×30m等，其地表温度表征的是像元所覆盖区域的平均陆地表面温度。对于像元所覆盖区域的均为植被时，该地表温度则代表了植被温度；对于像元所覆盖区域的均为裸地时，该地表温度则代表了裸地温度；但是，对于像元所覆盖区域有一部分区域为植被，而其他区域为裸地时，则无法直接区分该像元内的植被温度和裸地温度。由于长时间序列的卫星遥感观测数据(如NASAAVHRR-LST自1982年开始)的空间分辨率较为粗糙，像元内的植被覆盖度发生了显著的变化。因此，为得到更为精确的植被动态变化对环境和气候的影响信息，需要将遥感反演的地表温度进一步区分出植被温度和裸地温度。现有技术在估算植被温度或裸地温度时，主要采用邻近像元替代或空间插值法，如最近邻点插值法(NearestNeighbor)，自然邻点插值法(NaturalNeighbor)，克里金插值法(Kriging)等。例如，估算某一像元的植被温度时，用该像元附近的完全植被覆盖区的地表温度替代，或者进行空间插值得到该目标像元的植被温度。这类方法主要特点是计算简单，且目标像元的计算值高度依赖于邻近像元的有效个数及其和目标像元的相关性。因而，对于空间一致性越好的区域或者植被覆盖度高的区域，现有的空间插值方法得到的植被温度估算精度越高。但对于空间异质性越高的区域(如山地、丘陵以及植被类型过渡地带等)、或者植被覆盖度中等或越低的区域(如城市下垫面、荒漠、草地、稀疏森林等)，现有的空间插值方法得到的计算值误差就越大，导致无法得到可接受的植被温度和裸地温度估算值。现有的空间插值方法是一种统计算法，没有物理机制支撑，对于空间异质性高的区域(如山地、丘陵以及植被类型过渡地带等)、或者植被覆盖度中等或低的区域(如城市下垫面、荒漠、草地、稀疏森林等)，现有的空间插值方法得到的计算值误差较大，无法得到可接受的植被温度和裸地温度。本专利技术针对现有的计算方法没有物理机制支撑，无法适用于空间异质性高的区域等缺点，提出了一种以地表能量平衡原理为基础的，不依赖于其他像元信息，应用目标像元的卫星遥感数据(包括地表温度、地表发射率、植被指数、植被类型等)和地面气象观测数据为输入，计算出目标像元的植被温度和裸地温度的估算方法。
技术实现思路
为了解决上述技术所存在的不足之处，本专利技术提供了一种基于卫星遥感数据的植被温度和裸地温度估算方法。为了解决以上技术问题，本专利技术采用的技术方案是：一种基于卫星遥感数据的植被温度和裸地温度估算方法，包括以下步骤：步骤一、准备卫星遥感数据和地面气象观测数据；步骤二、根据步骤一的数据计算植被覆盖度；步骤三、根据步骤一的数据分别计算裸地反照率、裸地发射率和植被发射率；步骤四、计算裸地温度和植被温度。进一步地，卫星遥感数据包括但不限于卫星遥感的陆地地表温度TL、地表反照率αL、地表发射率εL、归一化植被指数NDVI和叶面积指数LAI、土地利用、植被类型数据；地面气象观测数据包括但不限于下行短波辐射SW↓和下行长波辐射LW↓数据。进一步地，植被覆盖度fV由遥感的归一化植被指数NDVI或叶面积指数LAI计算得到。进一步地，植被覆盖度fV的具体计算过程如下：其中，n通常取值为2，NDVImax和NDVImin为NDVI的最大和最小值；或者也可以利用LAI估算fV，其中，LAIref为参考叶面积指数，对于森林和灌木其取值为2.5，而对于草地、湿地和农田其取值为1.4；e为自然常数。进一步地，裸地反照率αS是从遥感地表反照率αL数据中找出同纬度或附件区域中无植被覆盖区的地表反照率αL的值；裸地发射率εS是从遥感地表发射率εL数据中找出同纬度或附件区域中无植被覆盖区的地表发射率εL的值；植被发射率εV是根据植被类型分布图，由遥感地表发射率εL数据中找出相同植被类型的完全植被覆盖区的地表发射率εL的值。进一步地，裸地温度TS是由气象观测的下行短波辐射SW↓、下行长波辐射LW↓、卫星遥感的陆地地表温度TL、地表反照率αL、地表发射率εL、叶面积指数LAI、裸地反照率αS、裸地发射率εS、植被发射率εV和参数k计算得到的，其具体计算公式如下：其中，σ＝5.67×10-8Wm-2K-4为斯蒂芬-波尔茨曼常数，k为消减系数,0＜k＜1，可根据不同植被类型估算最优参数值。进一步地，植被温度TV的具体计算公式如下：其中，TS为裸地温度，TL为卫星遥感的陆地地表温度，εL为地表发射率、fV为植被覆盖度，εV为植被发射率，εS为裸土发射率。本专利技术以地表能量平衡原理为基础，不依赖于其他像元信息，应用目标像元的卫星遥感数据，包括但不限于地表温度、地表发射率、植被指数、植被类型和地面气象观测数据为输入，计算出目标像元的植被温度和裸地温度，具有空间独立性，提高了估算精度。附图说明图1为本专利技术的整体流程示意图。图2为地表温度、植被温度和裸地温度随叶面积指数增加的变化关系图。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步详细的说明。图1所示的一种基于卫星遥感数据的植被温度和裸地温度估算方法，包括以下步骤：步骤一、准备卫星遥感数据和地面气象观测数据；步骤二、根据步骤一的数据计算植被覆盖度；步骤三、根据步骤一的数据分别计算裸地反照率、裸地发射率和植被发射率；步骤四、计算裸地温度和植被温度。步骤一中，卫星遥感数据包括但不限于卫星遥感的陆地地表温度TL、地表反照率αL、地表发射率εL、归一化植被指数NDVI和叶面积指数LAI、土地利用、植被类型数据；地面气象观测数据包括但不限于下行短波辐射SW↓和下行长波辐射LW↓数据。步骤二中的植被覆盖度fV由遥感的归一化植被指数NDVI或叶面积指数LAI计算得到；步骤三中的裸地反照率αS是从遥感地表反照率αL数据中找出同纬度或附件区域中无植被覆盖区的地表反照率αL的值；裸地发射率εS是从遥感地表发射率εL数据中找出同纬度或附件区域中无植被覆盖区的地表本文档来自技高网...

【技术保护点】
1.一种基于卫星遥感数据的植被温度和裸地温度估算方法，其特征在于：包括以下步骤：/n步骤一、准备卫星遥感数据和地面气象观测数据；/n步骤二、根据步骤一的数据计算植被覆盖度；/n步骤三、根据步骤一的数据分别计算裸地反照率、裸地发射率和植被发射率；/n步骤四、计算裸地温度和植被温度。/n

【技术特征摘要】
1.一种基于卫星遥感数据的植被温度和裸地温度估算方法，其特征在于：包括以下步骤：
步骤一、准备卫星遥感数据和地面气象观测数据；
步骤二、根据步骤一的数据计算植被覆盖度；
步骤三、根据步骤一的数据分别计算裸地反照率、裸地发射率和植被发射率；
步骤四、计算裸地温度和植被温度。


2.根据权利要求1所述的基于卫星遥感数据的植被温度和裸地温度估算方法，其特征在于：所述卫星遥感数据包括但不限于卫星遥感的陆地地表温度TL、地表反照率αL、地表发射率εL、归一化植被指数NDVI和叶面积指数LAI、土地利用、植被类型数据；地面气象观测数据包括但不限于下行短波辐射SW↓和下行长波辐射LW↓数据。


3.根据权利要求1所述的基于卫星遥感数据的植被温度和裸地温度估算方法，其特征在于：所述植被覆盖度fV由遥感的归一化植被指数NDVI或叶面积指数LAI计算得到。


4.根据权利要求3所述的基于卫星遥感数据的植被温度和裸地温度估算方法，其特征在于：所述植被覆盖度fV的具体计算过程如下：



其中，n通常取值为2，NDVImax和NDVImin为NDVI的最大和最小值；
或者也可以利用LAI估算fV，



其中，LAIref为参考叶面积指数，对于森林和灌木其取值为2.5，而对于草地、湿地和农田其取值为1....

【专利技术属性】
技术研发人员：张选泽，张永强，刘昌明，
申请(专利权)人：中国科学院地理科学与资源研究所，
类型：发明
国别省市：北京;11