使用遥感数据计算蒸散发所需文件

在计算之前一定要注意:

1.坐标系转换;
2.map尺寸;
3.输入文件中地表温度文件为开氏温度。

使用软件或代码均需要以下map:

1.反照率albedo

采用梁顺林2001年发表的转换公式:

A=0.356a2+0.130a4+0.373a5+0.085a6+0.072a7-0.0018

A=0.356b(Blue)+0.130b(Red)+0.373b(NIR)+0.085b(SWIR1)+0.072b(SWIR2)-0.0018

首先需要将landsat波段反射率数值除以10000,然后带入公式。即:

A=(0.356b2+0.130b4+0.373b5+0.085b6+0.072*b7)/10000-0.0018

2.归一化植被指数NDVI

NDVI = (NIR-R)/(NIR+R)

(在聚类分析时已获得NDVI map)

3.地表温度LST

参考基于单波段算法的地表温度反演教程

4.发射率emissivity(比辐射率)

在反演地表温度过程中已获得发射率map。

5.DEM

地理空间数据云下载(注意坐标系的投影转换以及图像的尺寸裁剪)。