基于无人机观测的高光谱 BRDF 可表征平坦沙漠地表的光学特性：与实验室和卫星数据的综合对比研究

Hyperspectral BRDF based on UAV measurements can characterize optical properties of flat desert surfaces: A comprehensive comparison with laboratory and satellite data研究单位东北师范大学摘要沙漠是陆地生态系统的重要组成部分其光学反射特性对于理解地表辐射平衡至关重要。将高光谱与角度观测相结合被认为是表征沙漠地表光学反射特性的有效方法。然而由于测量配置的限制地面和卫星观测仪器难以在合适的尺度上同时获取沙漠地表的高光谱和角度反射率。为克服这一局限性本文采用搭载光谱仪的无人机UAV对中国西北地区六处沙漠区域开展了高光谱双向反射分布函数BRDF测量。通过与在相同观测几何条件下获得的实验室和卫星观测结果进行对比发现无人机测得的高光谱反射率在光谱形态上与实验室测量结果高度一致而无人机获取的角度反射率在半球空间内的 BRDF 表征与卫星观测结果具有良好的一致性。此外本文还将基于 BRDF 模型反演得到的角度反射率与 MODIS 在选定波长下的沙漠实测角度反射率进行了比较。结果表明基于无人机测得的高光谱 BRDF 构建的 BRDF 模型在光谱 BRDF 模拟方面表现出更优的性能。基于贝叶斯反演的不确定性分析进一步证实了无人机数据在有效表征沙漠地表光学特性方面的可靠性。上述结果表明基于无人机的光谱与角度联合测量在定量刻画高光谱 BRDF 方面具有较强的稳健性可为获取沙漠地表全面的原位光学特性提供有力支撑并可作为未来沙漠地表辐射平衡研究的重要验证数据来源。研究方法图1 无人机观测位置及中国沙漠分布。 绿色三角形表示无人机观测位置阿拉善ALS、敦煌DH、格尔木GEM、民勤MQ、乌海WH和小柴旦西XCDH_W表1 飞行试验总结表2 六个研究地点的质地特性图2 基于无人机的高光谱 BRDF 评估流程示意图图3 六个荒漠地点的野外图像。(A)阿拉善(ALS)(B)敦煌(Dh)(C)格尔木(GEM)(D)民勤(MQ)(E)乌海(WH)和(F)小柴旦西部Est(XCDH_W)图4 本研究所用六种沙漠样品的显微图像。(a) 阿拉善ALS(b) 敦煌DH(c) 格尔木GEM(d) 民勤MQ(e) 乌海WH(f) 小柴旦西XCDH_W数据采集本研究的无人机BRDF测量是通过无人机系统和地面系统协同完成。无人机系统包括A660型无人机、云台控制系统、地面站以及RS-8800光谱仪。A660无人机最大载重为10Kg续航时间最长可达60分钟可配备多种云台搭载不同载荷如高光谱成像仪、多光谱成像仪、光谱辐射计、红外热成像仪、激光雷达系统、倾斜摄影、可见光成像等遥感设备可满足多种航测需求。RS-8800光谱仪光谱范围350-2500nm可快速实现野外紫外-可见-近红外全波段光谱数据测量采样间隔1nmFOV25°。地面系统包括RS-8800光谱仪、自动伸缩支架和控制箱。图5 (a) 无人机飞行方案竖直线表示飞行高度水平线表示飞行半径黑点表示观测天顶角黄色水平条表示天底方向测量的视场. (b) BRDF 测量布设红色虚线表示飞行路径黄色点表示观测位置观测天顶角间隔为 10°观测方位角间隔为 30°. (c) 无人机系统组成无人机、光谱仪、云台和相机图6 (a) 地面测量系统由光谱仪、可伸缩支架和控制箱组成。光谱仪的光纤探头和参考面板安装在可伸缩支架上。(b) 测量过程中展开的参考面板研究结果图7 (a) 无人机在天底方向VZA 0°SZA 30°观测得到的六个沙漠站点的半球方向反射函数HDRF。(b) 实验室在天底方向VZA 0°SZA 30°观测得到的六个沙漠站点的双向反射因子BRF。(c) 六个沙漠站点在主平面内不同观测天顶角VZA条件下无人机高光谱 HDRF 与实验室高光谱 BRF 之间的光谱角制图SAM。表中每个数值表示对应的 SAM 值SAM 值越小表明两种观测结果在光谱形态上的相似性越高。“–”表示前向散射方向。太阳天顶角SZA为 30°。图8 在 858 nm 波长下不同太阳天顶角条件下无人机获取的 HDRF 在半球空间中的角度分布第一行为 SZA 30°第二行为 SZA 40°第三行为 SZA 50°。黑色五角星表示起始太阳天顶角SZA_start。每一列表示一个沙漠站点从左到右依次为阿拉善ALS、敦煌DH、格尔木GEM、民勤MQ、乌海WH和小柴旦西XCDH_W。图9 在 858 nm 波长下不同太阳天顶角条件第一行为 SZA 30°第二行为 SZA 40°第三行为 SZA 50°下无人机观测得到的各向异性因子ANIF在半球空间中的角度分布。黑色五角星表示起始太阳天顶角SZA_start。每一列表示一个沙漠站点从左到右依次为阿拉善ALS、敦煌DH、格尔木GEM、民勤MQ、乌海WH和小柴旦西XCDH_W。图10 在 858 nm 波长下不同太阳天顶角条件第一行为 SZA 30°第二行为 SZA 40°第三行为 SZA 50°下无人机观测得到的 HDRF 在半球空间中的角度分布。黑色五角星表示起始太阳天顶角SZA_start。每一列表示一个沙漠站点从左到右依次为阿拉善ALS、敦煌DH、格尔木GEM、民勤MQ、乌海WH和小柴旦西XCDH_W。图11 六个站点在三个太阳天顶角SZA 30°、40°、50°和 61 个观测天顶角VZA条件下无人机观测的 HDRF波长范围 400–2450 nm与卫星观测的七个反射波段 HDRF 对比。黑点表示卫星观测 HDRF 的中位值误差线的上端和下端分别表示卫星观测 HDRF 的最大值和最小值。图12 在 858 nm 波长、太阳天顶角 SZA 40° 条件下无人机观测 HDRF 与卫星观测 HDRF 之间的相对差RD在半球空间中的角度分布第一行及卫星观测 HDRF 在半球空间中的角度分布第二行。平均 |RD| 为 61 个观测天顶角VZA绝对相对差的平均值。黑色五角星表示起始太阳天顶角SZA_start。每一列表示一个沙漠站点从左到右依次为阿拉善ALS、敦煌DH、格尔木GEM、民勤MQ、乌海WH和小柴旦西XCDH_W图13 六个站点在三个太阳天顶角SZA、61 个观测天顶角VZA及 1653 个波长条件下RPV 模型a和 SOILSPECT 模型b模拟 HDRF 与无人机实测 HDRF 的对比样本量 N 181,4994图14 本研究中基于无人机的所有数据下RPV 模型和 SOILSPECT 模型各参数的平均不确定性图15 (a) ABL 站点的实地照片 (b) 无人机高光谱 BRDF 测量的现场场景 (c) 无人机在 50 米高度拍摄的天底视角图像图16 在 858 nm 波长下ABL 站点太阳天顶角为 30°、40° 和 50° 条件下无人机观测 HDRF 的半球空间角度分布(a–c) 和 (d–f) 分别表示不同观测或处理方式的 HDRF。黑色五角星表示起始太阳天顶角SZA_start图17 在 858 nm 波长下ABL 站点太阳天顶角为 30°、40° 和 50° 条件下无人机观测 HDRF 与卫星观测 HDRF 之间的相对差RD在半球空间中的角度分布。黑色五角星表示起始太阳天顶角SZA_start研究结论在本研究中我们通过将无人机UAV观测的高光谱 BRDF 与实验室光谱和卫星观测 BRDF 进行对比评估了无人机测量在表征沙漠地表光学特性方面的有效性。本研究利用搭载光谱仪的无人机在中国西北地区六个沙漠区域的不同太阳天顶角SZA下采集高光谱 BRDF 数据。结果表明无人机测得的高光谱反射率在光谱形态上与实验室测量高度一致光谱角制图 SAM 0.074这验证了无人机在不同照明-观测几何条件下表征沙漠地表高光谱反射特性的准确性。无人机观测与卫星观测 HDRF 之间的相对差较小平均 |RD| ~ 5%进一步证明了无人机测量在野外条件下表征沙漠地表角度反射特性的可靠性。基于无人机数据构建的 BRDF 模型模拟精度高相对均方根误差 RRMSE 4%反演得到的角度参数不确定性低平均不确定性 2%验证了无人机高光谱 BRDF 与沙漠地表光学特性之间的良好匹配。无人机测量弥补了实验室测量与卫星遥感在地块尺度约 1-30 m高光谱 BRDF 测量上的空白。未来可以进一步探索基于无人机的高光谱 BRDF 在估算沙漠属性方面的潜在应用并尝试将其与卫星遥感结合用于沙漠属性制图。总之无人机测量在表征沙漠地表光学特性方面具有显著效果为沙漠辐射平衡研究提供了新的技术手段和宝贵的数据支持。原文链接https://doi.org/10.1016/j.rse.2025.115228