[发明专利]获取田块尺度的农田蒸散量的方法

权利要求书

1.一种获取田块尺度的农田蒸散量的方法，其特征在于，包括以下步骤：

S1、对原始空间分辨率的热红外遥感数据与原始空间分辨率的可见光遥感数据进行空间配准；

S2、利用经过空间配准后的热红外遥感数据和可见光遥感数据计算子像元温度；

S3、利用原始空间分辨率的可见光遥感数据计算农田参数，所述农田参数包括反照率、比辐射率和植被覆盖度；

S4、根据所述反照率、比辐射率及子像元温度计算农田净辐射通量；

S5、利用所述农田净辐射通量和植被覆盖度计算土壤热通量；

S6、计算农田显热通量H，并根据所述农田显热通量、农田净辐射通量和土壤热通量以及农田能量平衡方程计算得到田块尺度的农田蒸散量。

2.如权利要求1所述的方法，其特征在于，步骤S1具体为：对原始空间分辨率的热红外遥感数据与原始空间分辨率的可见光遥感数据的空间坐标进行转换，以使得这两种数据的像素在空间位置上一致。

3.如权利要求2所述的方法，其特征在于，进一步地，步骤S1具体为：在原始空间分辨率的可见光遥感数据中选择图像控制点，并根据图像控制点求解热红外遥感数据与可见光遥感数据像素之间坐标转换的二次多项式系数，然后利用所述二次多项式系数将原始空间分辨率的热红外遥感数据像素进行空间坐标转换，转换到原始空间分辨率的可见光遥感数据的地理投影空间，得到经空间坐标转换后的热红外遥感数据Tr，转换公式如下：

X＝a₀+a₁x+a₂y+a₃x²+a₄y²+a₅xy                         (1)

Y＝b₀+b₁x+b₂y+b₃x²+b₄y²+b₅xy                       (2)

其中：x，y分别为热红外遥感数据像素的横坐标、纵坐标；X，Y分别为可见光遥感数据像素的地理空间横、纵坐标；a₀，a₁，a₂...a₅、b₀，b₁，b₂...b₅分别为横、纵坐标转换二次多项式系数。

4.如权利要求1所述的方法，其特征在于，步骤S2具体包括：

S21、对原始空间分辨率的可见光遥感数据进行重采样，以使其空间分辨率与原始空间分辨率的热红外遥感数据相同；

S22、利用重采样后的可见光遥感数据计算植被指数NDVI_L，并建立该植被指数与原始空间分辨率的热红外遥感数据间的模型，即式(3)、(4)，并求解出模型系数a、b、c：

T′＝a+bNDVI_L+cNDVI_L²                   (3)

NDVI_L＝(LB_nir-LB_red)/(LB_nir+LB_red)      (4)

T′指原始空间分辨率的热红外遥感数据对应的像元温度，LB_nir和LB_red分别代表原始空间分辨率的可见光遥感数据中的近红外波段和红光波段；

S23、利用原始空间分辨率的可见光遥感数据计算植被指数NDVI，并将植被指数NDVI代入公式(3)计算得到高分辨率温度Ts；

S24、计算热红外遥感温度Tr与T′之差，并将差值作为误差项补充到Ts中得到子像元温度Tsub，即根据式(5)计算子像元温度：

Tsub＝Ts+(Tr-T′)                       (5)

得到子像元温度的空间分辨率与原始空间分辨率的可见光遥感数据的空间分辨率相同。