前言：
星载激光雷达数据由于覆盖范围较广,已经成为大区域森林垂直结构参数定量反演的最佳选择之一，但激光光斑分布离散。因此,需要融合其他卫星遥感影像(如Landsat-8、Sentinel-2等)实现大区域高分辨率森林冠层高度制图。我们可以通过GEE平台读取这些数据并进行处理最终实现森林冠层高度的制图。
1、数据准备与处理
数据集：（1）GEDI L2A Raster Canopy Top Height；（2）Sentinel-2 Surface Reflectance；（3）ESA WorldCover 10m v100；（4）NASA SRTM Digital Elevation 30m
1.1 sentinel-2数据处理
影像去云、计算指数、合成影像
1.2 计算坡度
坡度大于30度得到的冠层高度相对误差较大
1.3 GEDI数据处理
去除低质量、坡度大于30度的波形数据
2、建立回归模型
2.1 图层叠加，仅保留森林区域
2.2图层重采样
因为GEDI水平测量精度为 ± 9米，这会导致gedi的树高值与sentinel-2像元不匹配，重采样到一个更大的像元
2.3获取训练数据
gedi数据过于离散，我们要保证采样在gedi有值的地方，所以掩膜掉这些空值。采样1000个点，过多的话计算量太大可能会报错。
2.4 设置随机森林分类器
2.5 散点图
观测值与估测值的散点图
2.6 生成森林冠层高度图
以上代码借鉴了这篇文章：https://spatialthoughts.com/2024/02/07/agb-regression-gee/
