新疆大学地理与遥感学院马旭团队在遥感顶刊提出一种基于Sentinel-2数据估算干旱区生物覆盖度的多角度新方法

题目：A New Multiangle Method for Estimating Fractional Biocrust Coverage From Sentinel-2 Data in Arid Areas

作者：Hui Sun, Xu Ma, Ying Liu, Guiyun Zhou, Jianli Ding, Lei Lu, Tiejun Wang, Qiuli Yang, Qingtai Shu, and Fei Zhang

期刊：IEEE Transactions on Geoscience and Remote Sensing (SCI Q1)

一作单位：College of Geography and Remote Sensing Sciences, Xinjiang University, Urumqi, 830046, China

文献链接：http://dx.doi.org/10.1109/tgrs.2024.3361249

导读

新疆大学地理与遥感科学学院马旭老师团队硕士研究生孙卉在《IEEE Transactions on Geoscience and Remote Sensing》发表了题为《A New Multiangle Method for Estimating Fractional Biocrust Coverage From Sentinel-2 Data in Arid Areas》的研究论文，本研究旨在解决干旱区生物结皮覆盖度遥感产品的缺失的问题，并进一步为研究干旱区的生物结皮时空变化对全球气候变化的响应机理提供方法参考。

摘要

生物结皮（biocrust）的时空分布可用于监测沙漠生态系统中的区域水资源分布情况。然而，由于高空间分辨率的生物结皮覆盖度遥感数据的缺失限制了该领域的科学研究。为了解决这个问题，本研究建立了一个生物结皮估算模型（Estimation model for biocrusts，EMBC），从空间分辨率为10米的Sentinel-2影像中绘制了沙漠地区的生物结皮覆盖度图。首先，本研究基于辐射传输理论开发了一个生物结皮覆盖指数（FBCI）。接下来，建立了一个关于FBCI的多角度计算方程，并使用线性核驱动模型的反算方法获得了像元二分模型的参数。最后，使用关于FBCI的像元二分模型来计算生物结皮覆盖度。本研究使用场地测量数据来验证模型，并将估算结果与随机森林模型和反向传播神经网络模型的估算结果进行比较。结果表明，EMBC估算的生物结皮覆盖度与场地测量数据高度一致（均方根误差（RMSE）=0.0774，系统偏差=-4.05%）。此外，EMBC和其他两个模型估计的生物结皮覆盖度在空间分布上具有很高的一致性（RMSE<0.0998）。本研究提出的EMBC可以准确估算沙漠中的生物结皮覆盖度，可以成为监测干旱环境中干旱情况的重要技术。

结果

图1 生物结皮估算模型（Estimation model  for biocrusts，EMBC）的流程图

图2 滑动窗口算法。（a）提取2×2滑动窗口中4个像元的反射率的过程；（b）基于2×2滑动窗口提取的4个像元的反射率求解对应滑动窗口范围内的F的过程；（c）设置滑动步长为1个像元，计算2×2滑动窗口中像元值的过程，（d）对Sentinel-2遥感影像进行迭代计算得到平均F影像的过程。

图3 计算FBCImax和FBCImin的方法流程图：（a）棋盘分块和蒙特卡洛算法执行过程；（b）使用最小二乘算法求解多角度计算方程组的过程。

图4 FBCI阈值统计。（a）实验室测量的各组分的NDVI与FBCI的阈值；（b）不同土地覆盖类型中的FBCI的阈值。FBCI为生物结皮覆盖指数；NDVI为归一化植被指数；1表示草地覆盖的盐渍化区域；2为沙漠灌木覆盖区域；3为稀疏植被覆盖的区域；BC表示存在生物结皮覆盖；non-BC表示不存在生物结皮覆盖。

图5 使用生物结皮估算模型（Estimation model  for biocrusts，EMBC）估算的（a）3月26日，（b）4月18日，（c）5月28日，（d）6月22日，（e）7月22日和（f）8月21日的生物结皮覆盖度影像。影像的空间分辨率为10米，白色区域代表被掩膜的农田和水体区域。

图6 使用生物结皮估算模型（Estimation model  for biocrusts，EMBC）对三种生物结皮在不同生长阶段的覆盖度估算比较结果。C_BC代表蓝藻结皮，L_BC代表地衣结皮，M_BC代表苔藓结皮。

图7 使用生物结皮估算模型（Estimation model  for biocrusts，EMBC）估计的三种生物结皮的覆盖度与实地测量数据的比较结果。（a）蓝藻结皮；（b）地衣结皮；（c）苔藓结皮。FBC_M代表实地测量的生物结皮覆盖度；FBC_EMBC代表使用EMBC估计的生物结皮覆盖度；R2代表决定系数；RMSE代表均方根误差；SD代表相对误差。

图8 使用不同的算法从Sentinel-2影像中估算生物结皮覆盖度。（a）生物结皮估算模型（Estimation model  for biocrusts，EMBC）；（b）随机森林算法（Random forest，RF）；（c）反向传播神经网络（Backpropagation neural network，BPNN）；草地覆盖（短命植被）的盐碱化区域：（d）EMBC；（e）RF；（f）BPNN；沙漠灌木覆盖区域：（g）EMBC；（h）RF；（i）BPNN；稀疏植被覆盖区域：（j）EMBC；（k）RF；（l）BPNN。影像的空间分辨率为10米，白色区域代表被掩膜的农田和水体区域。

图9 使用不同的算法从Sentinel-2影像中估算生物结皮覆盖度的密度散点图。（a）EMBC vs RF；（b）EMBC vs BPNN；（c）BPNN vs RF。FBC_EMBC代表使用EMBC估算的生物结皮覆盖度；FBC_RF代表使用RF估算的生物结皮覆盖度； FBC_BPNN代表使用BPNN估算的生物结皮覆盖度；R2代表决定系数；RMSE代表均方根误差。

图10 不同算法估算的生物结皮覆盖度与实地测量数据之间的比较结果： （a）生物结皮估算模型（Estimation model  for biocrusts，EMBC）；（b）随机森林算法（Random forest，RF）；（c）反向传播神经网络（Backpropagation neural  network，BPNN）。FBC_M代表实地测量数据；FBC_EMBC代表使用EMBC估算的生物结皮覆盖度；FBC_RF代表使用RF估算的生物结皮覆盖度； FBC_BPNN代表使用BPNN估算的生物结皮覆盖度；R2代表决定系数；RMSE代表均方根误差。

结论

本研究提出了一种名为生物结皮估算模型（Estimation model for biocrusts，EMBC）的新方法来估算生物结皮覆盖度。在EMBC中，本研究基于辐射传输理论建立了新的生物结皮指数（即生物结皮覆盖指数（Fractional biocrust cover index，FBCI））。FBCI反映了生物结皮光谱与生物结皮覆盖度之间的数学关系。同时，我们采用线性核驱动模型从Sentinel-2影像中获取多角度反射率数据，扩展了基于多角度的方法，使其能够应用于Sentinel-2影像。通过使用EMBC，我们绘制了研究区域的生物结皮覆盖度图。结果表明，在沙漠场景中，EMBC能够利用Sentinel-2影像估算出生物结皮覆盖度，并且估算结果与实地测量数据以及两组对比验证数据高度一致。本研究提出的方法为监测生物结皮时空变化提供了新的视角，对深入分析生物结皮的生态效应以及评估沙漠化风险具有重要意义。

作者团队

第一作者：孙卉，新疆大学地理与遥感科学学院2021级硕士研究生，主要从事生物结皮的定量遥感研究。

第一作者、通讯作者：马旭，新疆大学地理与遥感科学学院副教授，主要从事植被遥感反射率建模研究；植被结构及光学仪器测量算法研究；植被参数遥感反演研究；植被群落和气候变化研究。通讯作者邮箱：maxu2020@xju.edu.cn；课题组网站：http://www.ma-lab.cn/。