【WRF数据准备】基于GEE下载静态地理数据-叶面积指数LAI及绿色植被率Fpar

GEE的介绍可参见另一博客 【GEE学习第一期】GEE介绍、注册及基本使用，本博客主要介绍基于GEE平台下载所需静态地理数据。

准备：WRF所需静态地理数据（Static geographical data）

主要的所需静态地理数据如下：

• Albedo
• 叶面积指数LAI: Leaf area index
• 光合有效辐射分数Fpar: Green vegetation fraction
• 不透水面积ISA: Impervious Surface Area
• 土地利用Land use

GEE的相关静态地理数据介绍可参见另一博客-【GEE数据库】WRF常用数据集总结。

数据范围说明

在确定研究区域范围（此博客以粤港澳大湾区为例）后，可下载处理所需的静态地理数据。

经纬度范围：

• 网格1：经度（106 120.58） 纬度（16.626 29.748）
• 网格2：经度（108.43 118.15） 纬度（18.813 27.561）
• 网格3：经度（111.346 115.477） 纬度（21.486 24.402）

更新，采用Lambert投影确定的研究区范围如下：

网格经纬度范围：上：33 下：14 左：102 右：123，下载数据年份选择2020年。

基于GEE下载叶面积指数及绿色植被率

叶面积指数LAI（leaf area index） 为单位土地面积绿色叶片的单面面积总和，即：叶面积指数=叶片总面积/土地面积，它与植被的密度、结构（单层或复层）、树木的生物学特性（分枝角、叶着生角、耐荫性等）和环境条件（光照、水分、土壤营养状况）有关，是表示植被利用光能状况和冠层结构的一个综合指标。

叶面积指数的测定有直接方法、间接方法。由于直接方法具有一定的破坏性，耗时耗力，并且无法反映大面积、大范围内的植被LAI分布情况，间接测定方法尤其是其中的遥感法，以其具有的快速、实时、大面积、不受植被类型影响等优势，成为应用最为广泛的LAI监测手段。

随着遥感技术的飞速发展，越来越多高时间分辨率、高空间分辨率、多光谱、高光谱遥感卫星数据实现了LAI监测。

GEE数据集介绍

Google Earth Engine ——MYD15A2H V6 MODIS综合叶面积指数（LAI）和光合有效辐射分数（FPAR）产品是一个500米分辨率的8天综合数据集。

• 叶面积指数（LAI）
• 光合有效辐射吸收比率（fraction of absorbed photosynthetically active radiation, FAPAR/FPAR）是描述太阳光在冠层辐射传输过程中植被吸收比例的参量， 一般定义为植被对波长在400 nm—700 nm 间太阳辐射能量的吸收比率。

数据介绍官网-MYD15A2H.061: Aqua Leaf Area Index/FPAR 8-Day Global 500m

GEE中代码：

ee.ImageCollection('MODIS/061/MYD15A2H')

数据下载：LAI（叶面积指数）和Fpar（绿色植被率）

1、选取空间范围

//var geometry = ee.FeatureCollection('文件地址/Boundary');
var geometry = ee.Geometry.Rectangle([-102,34,-84,48],'EPSG:4326',false);
//var geometry = ee.Geometry.Rectangle([-94.9910,40.6107,-92.3703,42.57729],'EPSG:4326',false)

说明：EPSG:4326是WGS 84坐标系统，也称为全球地球坐标系统(GCS)，是一种广泛使用的地理坐标系统。 这个坐标系统使用经度(longitude)和纬度(latitude)来定义地球上的位置，其中经度的范围是-180到+180，纬度的范围是-90到+90。

2、选取数据集来源：类型、时间范围等

• ee.ImageCollection(‘MODIS/006/MCD15A3H’)：加载 MODIS 版本 006 的 MCD15A3H 数据集，这个数据集提供了“叶面积指数”（LAI）的信息。
• .select(‘Lai’)：从数据集中选择“Lai”波段，提取叶面积指数。
• .filterDate()：过滤图像集合，只保留在 2012 年 1 月 1 日到 2012 年 12 月 31 日之间的图像。

var LAI = ee.ImageCollection('MODIS/006/MCD15A3H').select('Lai').filterDate(ee.Date('2012-01-01'), ee.Date('2012-12-31'))var Fpar = ee.ImageCollection('MODIS/006/MCD15A3H').select('Fpar').filterDate(ee.Date('2012-01-01'), ee.Date('2012-12-31'))

3、创建一个每月的总数集合
这段代码用于处理 MODIS 数据集中的 LAI 和 Fpar，通过按月计算每个月的平均值

• ee.List.sequence(1, 12)：生成一个包含从 1 到 12 的整数列表，表示一年中的每个月。

var months = ee.List.sequence(1, 12);
//print(months);
//Output: [1,2,3,4,5,6,7,8,9,10,11,12]

• months.map(function(m) {…})：对每个月（m）进行迭代。
• LAI.filter(ee.Filter.calendarRange(m, m, ‘month’))：从 LAI 图像集合中过滤出当前月份的图像。
• .mean()：计算该月份图像的平均值。
• .set(‘month’, m)：为计算出的平均图像设置一个属性，标记其对应的月份。
• ee.ImageCollection.fromImages(…)：将所有月份的平均图像组合成一个新的图像集合。

var byMonth_LAI = ee.ImageCollection.fromImages(months.map(function(m) {return LAI.filter(ee.Filter.calendarRange(m, m, 'month')).mean().set('month', m);
}).flatten());

这部分代码与计算 LAI 平均值的部分类似，只是它处理的是 Fpar 图像集合。
通过相同的逻辑，过滤出每个月的图像并计算平均值，最终生成一个包含每个月 Fpar 平均值的图像集合。

var byMonth_Fpar = ee.ImageCollection.fromImages(months.map(function(m) {return Fpar.filter