利用 LiDAR 点云数据或多源遥感影像,结合 QGIS 和 Python 估算森林树高的完整流程。
导入 Sentinel-2 影像、Landsat 影像、DEM 或 LiDAR/CHM 数据,检查坐标系、分辨率和有效范围。
确保所有数据统一到同一坐标系,建议统一到研究区对应的 UTM 投影。
对 Sentinel-1/2 或 Landsat 数据进行必要的大气校正、几何校正、裁剪和云掩膜处理。
Sentinel-2 Level-2A 产品已完成大气校正,可重点检查云掩膜和研究区裁剪。
提取 NDVI、EVI、纹理特征、DEM、坡度、坡向等变量,构建样本点或像元级特征矩阵。
使用 rasterio、geopandas、rasterstats 可以批量提取样本点对应的栅格值。
以 CHM 或实测树高为标签,使用随机森林、XGBoost 或线性模型进行训练和验证。
至少报告 R²、RMSE、MAE,并区分训练集和验证集结果。
将树高预测结果导出为 GeoTIFF,在 QGIS 中进行专题制图和结果检查。
注意去除异常值和无效区域,避免把无 CHM 覆盖区域误判为低树高。