一、数据简介
全球土地蒸发阿姆斯特丹模型 (GLEAM)是一套算法,可分别估计土地蒸发的不同组成部分(通常称为“蒸发蒸腾”):蒸腾、裸土蒸发、截留损失、开放水域蒸发和升华。
此外,GLEAM 还提供表面和根部土壤湿度、潜在蒸发和蒸发应力条件。该方法的基本原理是最大限度地恢复当前气候和环境变量卫星观测中包含的蒸发信息。
GLEAM 中的普里斯特利和泰勒方程根据对地表净辐射和近地表空气温度的观测来计算潜在蒸发量。基于微波植被光学深度 (VOD) 的观测和根区土壤湿度的估计,使用乘法蒸发应力因子将裸土、高冠层和短冠层的土地部分的潜在蒸发估计值转换为实际蒸发量。后者是使用多层流水平衡来计算的。
为了尝试纠正随机强迫误差,对表面土壤湿度的观测也被吸收到土壤剖面中。拦截损失是在 GLEAM 中使用 Gash 分析模型单独计算的。
最后,对水体和被冰和/或雪覆盖的区域的实际蒸发量的估计是基于修正的普里斯特利和泰勒方程。
主要特征:考虑土壤水分对蒸发的约束;森林拦截的详细参数化;广泛使用微波观测,在多云条件下具有优势。
二、数据下载
官网下载:www.gleam.eu
ftp下载:SFTP
Host: sftp://hydras.ugent.be
Port: 2225
Username: gleamuser
Password: v37_GLEAM2023#aw
三、数据分享
我已经下载全球1982~2022年的月和年尺度的数据,和各位免费分享。GLEAM全球1982~2022年数据免费分享 (qq.com)
四、数据处理
nc文件需要转换为tif,这是python代码。
# -*- encoding:utf-8 -*-
import arcpy
import os
# Set local variables
arcpy.CheckOutExtension("ImageAnalyst") # 检查许可
arcpy.CheckOutExtension("spatial") # 检查许可
# path = r"F:\modis_process\test\aod"
path = r'E:\EVI'
out_file = r"E:\EVI"
variable = "etp"
XDimension = "lon"
YDimension = "lat"
# Execute MakeNetCDFRasterLayer
dirs = os.listdir(path)
print(dirs)
dir_path = []
for dir in dirs:
path_1 = path + "//" + dir
dir_path.append(path_1)
print(dir_path)
i = 0
for nc in dir_path:
print(nc)
nc_file = dirs[i].strip(".nc")
print(nc_file)
arcpy.MakeNetCDFRasterLayer_md(nc, variable, XDimension, YDimension, nc_file)
outNC = nc.strip(".nc")
out_path = outNC + ".tif"
copy_raster = arcpy.CopyRaster_management(nc_file, out_path, background_value="", nodata_value="")
print(out_path)
i += 1
print("over")