RPI-GMM:一种基于结构且与物候无关的新算法,用于绘制最新10米分辨率的国家级橡胶种植园分布图

《Remote Sensing of Environment》:RPI-GMM: A novel structure-based and phenology-independent algorithm for mapping latest 10-m resolution national-level rubber plantations

【字体: 时间:2026年01月16日 来源:Remote Sensing of Environment 11.4

编辑推荐:

  橡胶种植园遥感监测中提出SAR光学融合的RPI-GMM框架,通过双极化SAR回波与SWIR波段构建结构敏感指数,结合无监督高斯混合模型实现物候无关分类,在老挝验证精度达87%,揭示70%种植园位于边境山地等新特征。

  
肖驰伟|岳子龙|冯志明|董金伟|卢 Juliet|Khin Htet Htet Pyone|Khampheng Boudmyxay
中国科学院地理科学与自然资源研究所,北京,中国

摘要

准确且最新的橡胶种植园地图对于生态环境和社会经济影响评估以及可持续的农林业管理非常有益。然而,现有的基于遥感的技术主要依赖于时间序列光学数据中的物候信号来识别橡胶种植园,这些方法受到持续云层覆盖、区域物候变化或不一致性以及高数据需求的限制。为了解决这些问题,我们提出了一种创新的、不依赖于物候的框架,该框架将橡胶种植园指数(RPI)与无监督高斯混合模型(GMM)分类器相结合。RPI是一个基于双极化Sentinel-1 SAR后向散射(VV/VH)和Sentinel-2短波红外反射率(第11波段)的结构敏感指数,能够捕捉种植园的规律性和冠层湿度特征。我们在六个不同类型的橡胶种植园样本区域评估了RPI-GMM框架的效果,这些区域在物候、地形和种植园结构上存在差异。结果显示,在物候强烈和较弱的情况下,以及在不同类型的地形中,分类准确率均超过0.87。我们的RPI-GMM方法在老挝的应用中取得了87.0%的总体准确率,并估计出2024年共有234,206公顷的橡胶种植园。空间分析显示,大约70%的橡胶种植园位于靠近中国和越南的边境地区,90%的种植园位于海拔1000米以下,80%的种植园位于坡度在3°到16°之间的区域。值得注意的是,我们这种简单且集成的RPI-GMM方法不需要时间序列数据或标记数据,确保了方法的稳健性、成本效益和可转移性。研究结果强调了基于结构的SAR-光学融合技术在未来全球或热带地区监测树木种植园动态方面的价值,并支持其在农林业管理中的更广泛应用。

引言

天然橡胶(Hevea brasiliensis)是一种有价值的经济作物和出口商品,其种植面积在过去三十年中迅速扩大,尤其是在东南亚(SEA)(Wang和Zhang,2025)。在SEA,从1993年到2016年,大约有410万公顷的天然森林被转换为橡胶种植园,其中许多位于生物多样性热点地区(Wang等人,2023)。虽然这种扩张促进了区域经济增长并改善了小农户的生活,但也引发了严重的生态环境问题,包括森林砍伐、生物多样性丧失、碳封存减少(Grass等人,2020)以及水文平衡问题(Fox等人,2014;Ziegler等人,2009)。因此,迫切需要可靠的监测方法来生成高分辨率的橡胶种植园地图。同时,一些国际框架(如欧盟森林砍伐法规(EUDR)直接针对由单一作物种植园扩张引起的森林砍伐问题,对新型和可扩展的遥感技术提出了越来越高的需求,以监测空间上准确且最新的信息(Azizan等人,2021b;Ye等人,2018;Yue和Xiao,2025)。
到目前为止,大多数现有的橡胶种植园遥感映射方法都依赖于物候特征(Kou等人,2015;Chen等人,2018;Xiao等人,2020a)。与常绿热带森林相比,橡胶树表现出明显的季节性落叶行为,在植被指数(如NDVI和RESI)中产生了特征性的时间变化(Abd Razak等人,2018;Fan等人,2015;Xiao等人,2020a)。这些物候差异已被广泛用于分类(Chen等人,2016;Xiao等人,2020b)。早期基于MODIS和Landsat图像的研究使用了最大值合成、归一化差异指数和时间差分等技术来描述区域尺度的橡胶季节性(Chen等人,2018;Li和Fox,2012)。随着高分辨率传感器(如Sentinel-2)和Google Earth Engine平台的出现,基于物候的映射的准确性和自动化水平显著提高,使用了样条拟合和随机森林等模型(Liu等人,2020;Shen等人,2023;Zhai等人,2024)。这些方法在物候稳定的地区(如中国西南部的西双版纳、越南中南部和泰国北部)取得了相当大的成功(B. Chen等人,2025a;Dong等人,2012;Xiao等人,2021)。然而,在许多热带环境中,基于物候的方法仍面临重大限制(Yue和Xiao,2025)。
重要的是,橡胶种植园在纬度、海拔和气候梯度上表现出多样的物候行为或不一致性(Azizan等人,2022)。例如,在印度尼西亚和马来西亚南部等赤道地区,橡胶树的落叶行为可能较为微妙或不一致,这使得与周围植被的物候区分变得困难(Azizan等人,2021a)。此外,热带地区的持续云层覆盖经常阻碍高质量光学时间序列数据的获取,影响物候曲线的可靠性(Azizan等人,2021b)。尽管基于物候的方法在某些地区取得了成功,但由于环境的复杂性和物候行为的多样性,其在热带地区的准确性仍然不确定(Hoang等人,2025;Sheil等人,2025)。因此,开发一种可扩展的、不依赖于物候的大面积橡胶种植园检测方法仍然是一个重要的研究挑战。
为了克服这些挑战,最近的研究开始探索利用结构特征进行种植园分类的潜力(Cheng等人,2023)。例如,在油棕种植园的映射中,规则的行种植在光学图像中产生了强烈的线性纹理,可以使用灰度共生矩阵(GLCM)或基于对象的图像分析等技术来捕捉(Fujiki等人,2016;Zhang等人,2020)。然而,特别是由小农户经营的橡胶种植园,在种植密度、年龄结构和边界清晰度方面往往更加异质,导致空间纹理较弱且不均匀(Cui等人,2022;Jia等人,2020)。因此,传统的基于纹理的方法在橡胶种植园映射中表现不佳。合成孔径雷达(SAR)具有全天候和全天成像的能力,在多云的热带地区提供了一个有前景的替代方案(Hashemi等人,2024)。迄今为止,基于SAR的橡胶映射主要利用极化后向散射信号(特别是双极化VV/VH数据)来区分森林和非森林(Chen等人,2020;Chen等人,2023;le Maire等人,2025;Trisasongko,2021)。这些方法在物候稳定的地区(如越南中南部和泰国北部)取得了显著的成功(B. Chen等人,2025a;Dong等人,2012;Xiao等人,2021)。然而,在许多热带环境中,基于物候的方法仍面临重大限制(Yue和Xiao,2025)。
重要的是,橡胶种植园在纬度、海拔和气候梯度上表现出多样的物候行为或不一致性(Azizan等人,2022)。例如,在印度尼西亚和马来西亚南部等赤道地区,橡胶树的落叶行为可能由于落叶时间的变化而显得微妙或不一致,使得与周围植被的物候区分变得困难(Azizan等人,2021a)。此外,热带地区的持续云层覆盖经常阻碍高质量光学时间序列数据的获取,影响物候曲线的可靠性(Azizan等人,2021b)。尽管基于物候的方法在某些地区取得了成功,但由于环境的复杂性和物候行为的多样性,其在热带地区的准确性仍不确定(Hoang等人,2025;Sheil等人,2025)。因此,开发一种可扩展的、不依赖于物候的大面积橡胶种植园检测方法仍然是一个重要的研究挑战。
为了克服这些挑战,最近的研究开始探索利用结构特征进行种植园分类的潜力(Cheng等人,2023)。例如,在油棕种植园的映射中,规则的行种植在光学图像中产生了强烈的线性纹理,可以使用灰度共生矩阵(GLCM)或基于对象的图像分析等技术来捕捉(Fujiki等人,2016;Zhang等人,2020)。然而,特别是由小农户经营的橡胶种植园,在种植密度、年龄结构和边界清晰度方面往往更加异质,导致空间纹理较弱且不均匀(Cui等人,2022;Jia等人,2020)。因此,传统的基于纹理的方法在橡胶种植园映射中往往表现不佳。合成孔径雷达(SAR)具有全天候和全天成像的能力,在多云的热带地区提供了一个有前景的替代方案(Hashemi等人,2024)。迄今为止,基于SAR的橡胶映射主要利用极化后向散射信号(特别是双极化VV/VH数据)来区分森林和非森林(Chen等人,2020;Chen等人,2023;le Maire等人,2025;Trisasongko,2021)。这些方法可以有效识别密集的自然森林,但由于橡胶种植园在种植密度、树龄和边界定义上的变化,通常无法直接将其与其他基于树木的土地覆盖类型(如油棕种植园)区分开来。与油棕种植园的强健且均匀的后向散射信号相比(Descals等人,2024),橡胶种植园在SAR图像中显示出更复杂和不规则的结构特征,从而使得提取独特模式变得复杂。到目前为止,还没有一种通用框架专门针对橡胶种植园的结构特征进行提取。
尽管有一些研究尝试使用景观指标(如斑块大小、冠层密度和边界复杂性)来区分小农户和工业橡胶系统(Hazir和Muda,2020;Kou等人,2015;le Maire等人,2025),但这些方法高度依赖于空间分辨率和当地的种植实践。使用较低分辨率(如500米或250米)的研究往往无法充分了解小农户系统的细节,因为它们无法捕捉到这些种植园的细微异质性(Dong等人,2012;Li和Fox,2012)。目前,大多数研究依赖于30米的分辨率,但即使在这种分辨率下,仍然缺乏一种稳健且可转移的结构基础方法,限制了它们在异质和破碎的热带地区的应用(B. Chen等人,2025a;Xiao等人,2021)。特别是,在平坦和山区之间进行区分时,橡胶种植园通常规模较小、形状不规则且分布稀疏,使用现有模型难以检测(Chen等人,2023)。这些限制在地形复杂和生物多样性丰富的地区尤为明显。此外,最近在物候研究方面的进展(Yue和Xiao,2025)表明,结合时间动态和更精细的空间特征可能显著提高在复杂景观中检测和分类橡胶种植园的能力。
在这里,我们提出了一种新的框架,该框架结合了来自Sentinel-1 SAR的结构敏感特征和来自Sentinel-2 SWIR数据的冠层湿度信息,以准确识别橡胶种植园。具体来说,我们通过结合Sentinel-1双极化后向散射(VV/VH)和Sentinel-2短波红外反射率(第11波段)构建了一个橡胶种植园指数(RPI),并使用高斯混合模型(GMM)进行无监督分类。这种方法利用SAR揭示的结构信息和SWIR的冠层湿度敏感性,以成本高效、不依赖于物候且空间可转移的方式识别橡胶种植园。我们在多个具有不同地形、物候和栽培系统的代表性热带地区评估了RPI-GMM方法的效果。最后,我们制作了老挝首个全国范围的10米分辨率橡胶种植园地图。实验结果证明了我们的RPI-GMM算法具有高分类准确性,为树木种植园动态的监测提供了可扩展且稳健的解决方案,支持可持续的农林业管理以及经济和商业运营实践。

方法片段

Sentinel-1合成

Sentinel-1是由欧洲空间局(ESA)发射的C波段SAR卫星,提供10米空间分辨率的全天候后向散射数据,重访周期约为六天(Torres等人,2012)。我们使用了地面距离检测(GRD)产品并提取了VV和VH极化通道。按照合成工作流程(Descals等人,2021),预处理包括轨道分离、辐射校准、热噪声去除、地形校正和入射角校正

方法论

RPI-GMM算法的工作流程如图2所示。它包括三个主要阶段:(1)RPI的开发,(2)应用高斯混合模型(GMM)进行无监督分类,以及(3)结果地图的准确性评估。这些步骤共同构成了本研究的方法论框架。

在不同物候条件下的RPI可分性

选择了两个代表性样本站点来检验RPI在不同物候条件下的可分性(图3):S1位于老挝北部,那里的橡胶种植园表现出强烈的落叶物候;S2位于印度尼西亚南部,那里的物候信号较弱,橡胶与其他植被类型更难以区分。
在S1或老挝北部,橡胶种植园显示出明显的单峰分布,RPI值集中在0.8以上

从特定于橡胶的RPI到通用的雷达种植园指数

在种植园映射领域,传统方法严重依赖于橡胶树的落叶物候,利用时间序列光学数据来跟踪落叶和重新绿化阶段以进行物种级别的区分(Abd Razak等人,2018;Dong等人,2012;Ye等人,2018)。尽管在某些情况下有效,但这些方法在热带地区受到持续云层覆盖、年际气候变化和物候表达不一致性的限制

结论

本研究提出了一种基于结构且不依赖于物候的橡胶种植园映射框架,结合了Sentinel-1 SAR和Sentinel-2 SWIR数据以及GMM分类器。所提出的RPI有效地捕捉了橡胶种植园独特的结构规律性和冠层湿度特征,使得在各种环境条件下都能进行准确分类。我们在六个具有不同物候条件的代表性热带地区进行了评估

CRediT作者贡献声明

肖驰伟:撰写 – 审稿与编辑,撰写 – 原稿,可视化,验证,资源管理,方法论,研究调查,资金获取,正式分析,概念化。岳子龙:撰写 – 审稿与编辑,撰写 – 原稿,可视化,验证,软件开发,方法论,正式分析,数据管理。冯志明:撰写 – 审稿与编辑,可视化,资金获取。董金伟:撰写 – 审稿与编辑,验证。

利益冲突声明

作者声明他们没有已知的财务利益或个人关系可能影响本文的报告。

致谢

我们衷心感谢三位匿名审稿人的建设性评论和建议,这些评论和建议极大地提高了我们工作的质量。本工作得到了国家自然科学基金(42371280、42130508和42001226)、中国-东盟卫星遥感应用重点实验室(KLCARS-2025-G03)以及研究所卓越青年科学家计划的支持
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

知名企业招聘

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号