基于 GEE 降尺度监测叶面积指数 LAI

目录

1 叶面积指数（LAI）

2 完整代码

3 运行结果

---

1 叶面积指数（LAI）

叶面积指数（LAI）定义为单位面积上所有叶子表面积的一半，是研究植被垂直分布的一个重要参数，它调控着许多生物物理过程，同时也为植被冠层表面最初的能量交换提供结构化定量信息，可直接作为量化植被冠层结构和生物特征的指标。因此叶面积指数监测对于监测农作物的生长态势和产量估算、维持森林系统的生态平衡和固碳能力、以及评估城市生态效力等方面具有十分重要的意义。

叶面积指数遥感监测方法：

1. 统计模型法：以光谱数据或其植被指数作为自变量建立的估算模型。

2. 光学模型法：基于植被的BRDF、建立在辐射传输模型基础上的一种模型，具有相当强的物理基础，不依赖于植被类型、背景环境，具有普适性。其辐射传输模型为冠层反射率与太阳和观测的天顶角、方位角，以及关于植被冠层的物理参数。但辐射传输方法复杂，无 法直接得到LAI，需要间接使用反演方法得到结果，常用的反演方法包括查找表法、神经网络法。

3. 激光雷达能量法：基于激光雷达传输方程，使用激光雷达全波形数据记录的回波能量信息计算植被间隙率，最后通过beer定律计算出叶面积指数。

降尺度（根据LAI与NDVI的关系降尺度）：

假定LAI与NDVI之间是线性关系，并且不同空间尺度也呈线性关系，而且大尺度NDVI像元值等于组成它的小尺度像元值的平均值，即：

代码逻辑为：

1. 准备modis LAI的500m数据、modis植被指数250米数据。
2. modis植被指数重采样500米，计算α。
3. 用modis植被指数250米乘以α，得到250米LAI。

2 完整代码

//定义研究区范围
var geometry = ee.FeatureCollection("projects/ee-zhangkanghnust/assets/HuNan");

Map.centerObject(geometry,7)
var startdate=ee.Date.fromYMD(2023, 7, 1);
var enddate=ee.Date.fromYMD(2023, 7, 20);

var datasetL = ee.ImageCollection('MODIS/061/MCD15A3H')
                  .filterBounds(geometry)
                  .filterDate(startdate,enddate)
var datasetN = ee.ImageCollection('MODIS/061/MOD13Q1')
                  .filterBounds(geometry)
                  .filterDate(startdate,enddate)
// 得到500米lai和250米ndvi
var LAI_composite=datasetL.median().clip(geometry)
var LAI500 = LAI_composite.select('Lai').multiply(0.1);
var dataNdvi = datasetN.first().clip(geometry).select('NDVI').multiply(0.0001)

// 重采样得到500米ndvi,注意这里设置的投影，投影很重要，否则无法显示     
var ndvi_proj= dataNdvi.projection()
var ndvi500 = dataNdvi.setDefaultProjection(ndvi_proj)
                 .resample('bilinear').reproject({
                      crs:ndvi_proj,
                      scale: 500
                    });
// 防止ndvi为0，无法除 
var epsilon = 1e-10;  
var ndvi500 = ndvi500.where(ndvi500.neq(0), ndvi500).where(ndvi500.eq(0), epsilon);  
// 得到比率  
var percent=LAI500.divide(ndvi500)
print("percent",percent)
// 得到250米的lai
var image250=ee.Image.constant(1).reproject({  
  crs: dataNdvi.projection().crs(),   
  scale: 250 
}).clip(geometry);
var percent250 = image250.multiply(percent)
var LAI250 = percent250.multiply(ndvi500)
var palettes = require('users/gena/packages:palettes')  
var palette01=palettes.colorbrewer.RdYlGn[11]
Map.addLayer(LAI250, {min:0,max:10,palette:palette01}, 'lai')

3 运行结果

研究区（湖南省）叶面积指数计算结果