关键词:
地表反照率
平坦地表
山区地表
哨兵二号
机器学习
遥感反演
摘要:
地表反照率定义为地表反射太阳短波辐射与入射短波辐射之比,是地球表层能量循环的关键参数之一,驱动地表辐射强迫变化,进而对区域地表气候变化产生影响。遥感卫星为全球尺度地表反照率制图提供了关键数据源。在过去几十年中,利用遥感卫星数据已经生产了区域/全球覆盖范围的百米至千米空间尺度的地表反照率遥感产品,如500m中分辨率成像光谱仪(MODIS)MCD43A3产品、0.5°多角度成像光谱仪(MISR)反照率产品和1000m全球陆地表面卫星(GLASS)反照率产品等。这些反照率产品已被应用于提供区域/全球范围的反照率观测,供气候变化监测和生物圈模型应用、辐射强迫估算、物候学研究等,总体精度在0.02-0.05之间。对于反照率产品的应用,如森林管理中特征斑块反照率分布、全球农业的典型农田大小或城市地区的主要不透水表面范围反照率的变化和分布等,这种粗尺度的反照率产品并不十分适合。此外,地表反照率遥感产品受地形效应影响严重,已有的反照率产品都没有考虑到地形影响,使得这些产品在山区的误差较大,难以获得满意的结果。高分辨率反照率具有更高的空间分辨率(从10m到100m)的优势,可以填补反照率产品从全球范围到区域尺度应用的空白,特别有助于量化地块尺度上土地利用变化引起地表反照率变化对农田生态系统干扰的影响。Sentinel-2卫星可以提供可见光至近红外波段精细尺度的观测值,空间分辨率为10m,时间分辨率约为5天,为构建高分辨率的地表反照率产品提供了重要数据源。
针对反照率产品空间分辨率低,没有考虑地形影响这两个问题,本论文以构建高分辨率地表反照率遥感产品为目标。首先,通过比较四种数据驱动的机器学习算法(人工神经网络(ANN)、K近邻算法(KNN)、随机森林(RF)和XGBoost(XGBT))的性能来评估最优的机器学习反演地表反照率算法;进而,利用Sentinel-2卫星数据,分平地和山地两类地形特征情况,反演高分辨率地表反照率遥感产品;最后,利用最优的反演模型,借助于Google Earth Engine(GEE)平台,生成了青藏高原2019年高分辨率地表反照率产品,并将其重采样到500m后与MODIS反照率产品进行比较,评价其时空一致性,同时也比较了Sentinel-2和MODIS地表反照率在冰雪地表的精准体现。本文的主要结论如下:
(1)在平坦地表,机器学习模型(ANN、KNN、RF、XGBT)在反演地表反照率方面具有很好的性能。通过测试数据与验证数据的检验,四个模型的精度都非常高,Bias范围为-0.0024~0.0018,R2高于0.99,RMSE范围在0.0055~0.0126。通过与地面站点反照率测量值相比,四类算法反演的Sentinel-2高分辨率反照率显示出很高的精度,Bias范围为0.0007~0.0113,R2高于0.9465,RMSE范围为0.0308~0.0335。其中,RF模型反演精度优于ANN、KNN和XGBT模型,Bias是0.0027,R2为0.9472,RMSE为0.0308。
(2)在山区表面,机器学习模型(ANN、KNN、RF、XGBT)在反演地表反照率方面也具有很好的性能。通过测试数据在不同坡度范围内进行检验,四个模型具有很高的精度。整体的Bias范围在-0.002~0.0023,R2高于0.96,RMSE范围在0.0048~0.0092。通过与地面站点反照率测量值相比,四类算法反演的Sentinel-2高分辨率反照率显示出很高的精度,Bias范围为-0.002~-0.0009,R2高于0.5871,RMSE范围为0.0254~0.0272。其中,RF模型反演精度优于ANN、KNN和XGBT模型,Bias是-0.002,R2为0.6162,RMSE为0.0254。
(3)青藏高原区域,Sentinel-2地表反照率与MODIS产品有很高的一致性,黑空反照率的Bias为0.003,R2为0.9387,RMSE为0.0282。地表反照率呈现季节性变化,冬季反照率值最高,春季开始降低,夏季将至最低值,秋季反照率值增加,冬季反照率增加到最高。Sentinel-2地表反照率在冰雪地表上有很高的准确性,可以很好的表现出冰雪地表的地表反照率分布并反映其变化。