弱监督时序特征驱动的遥感耕地产品融合制图方法技术

本发明专利技术公开了弱监督时序特征驱动的遥感耕地产品融合制图方法，包括以下步骤：弱监督时序数据集构建；构造由弱监督损失改造后的utae时序网络模型；弱监督模型训练，得到初步的模型后再针对难样本进行补充训练；最后使用训练好的模型进行适当范围的耕地覆盖制图。本发明专利技术将弱监督学习思想应用于耕地产品融合制图任务，通过引入时序模型提高了产品融合精度。度。度。

【技术实现步骤摘要】
弱监督时序特征驱动的遥感耕地产品融合制图方法


[0001]本专利技术属于遥感耕地产品的制图
，尤其涉及弱监督时序特征驱动的遥感耕地产品融合制图方法。

技术介绍

[0002]耕地作为人类生存和社会发展的基本资源，其面积和分布受到我国政府的密切关注，要严守18亿亩耕地红线，及时准确的耕地制图至关重要。随着遥感技术的快速发展，大批不同时空分辨率，光谱分辨率的卫星传感器的升空为更快速更高精度的地表资源监测提供了机会。目前，有两类产品能够为耕地监测服务，耕地专题制图产品与全球土地覆盖地图。
[0003]对于耕地专题图制作，各国开展了一系列相关项目，常见的耕地制图产品有美国的CDL(Cropland Data Layer)作物分类图，加拿大的AFCC产品以及欧洲的Sen
‑
Agri，由于这些地图的制作往往需要制作足量的标签数据用于模型的学且，且需要大量的实地调查验证产品可靠性，因此产品覆盖范围只局限在上述项目实施地区，通常难以大范围的进行耕地专题地图的制作。
[0004]对于全球土地覆盖产品，不同国家以及组织发布了许多区域或全球范围的土地覆盖产品。例如，谷歌的10m近实时的全球地表覆盖项目：Dynamic World，Esri公司发布的2017
‑
2021为期5年的地表覆盖产品，清华大学的FROM
‑
GLC等。这些高空间分辨率的地表覆盖产品的出现为全球尺度下的耕地监测提供了机会，但是由于这类产品首要保障产品全球覆盖且综合保障所有土地覆盖类型整体的精度，因此为保证遥感数据的全球覆盖，常使用多个影像拼接而成单一时相数据或时间间隔较长的数据进行分类。而耕地这一覆盖类型的光学特征随季节改变有明显的变化，相同地区不同时间的耕地影像在纹理，色彩等视觉特征也会发生很大改变，因此使用单一时相的遥感影像制图得到的土地覆盖结果中耕地的精度并不高。
[0005]由于耕地的分布和面积一般不会随时间有剧烈的变化，为了低成本的进行高精度的全球耕地覆盖产品制作和更新，许多研究者尝试将多个土地覆盖产品根据一定的标准融合来生成新的精度更高的耕地专题图。下面从产品融合的数据源和产品融合方法两方面来介绍当前的耕地产品融合技术现状与技术不足之处。
[0006]产品融合数据源：遥感影像的空间分辨率限制了土地覆盖产品的制图精度，早期土地覆盖产品以较低的空间分辨率为主，随着欧洲哨兵系列卫星升空，我国的高分、资源系列卫星相继开放共享数据，土地覆盖产品的制作迈进了中、高空间分辨率时代。随着机器学习方法和遥感云计算平台的发展，高空间分辨率，高频次的全球土地覆盖制图成为可能，当前常用的几个全球土地覆盖产品如下：
[0007]谷歌的Dynamic World产品：谷歌公司在2022年发布，其空间分辨率为10m，覆盖全球范围，时间跨度为2015
‑
2022年，使用哨兵2号的1C级产品作为影像，该产品使用深度学习模型作为分类算法，以24000个全球各地标注图像为训练基础，凭借谷歌自身的云平台高强
度算力，做到了近实时的土地覆盖制图。
[0008]欧空局的World Cover产品：ESA在2021年发布的产品，WC产品发布了2020年与2021年两期，空间分辨率为10m，覆盖全球范围，WC使用哨兵2号和哨兵1号作为基础数据，综合使用了OSM等开源地理数据作为辅助，用随机森林分类树算法在100
×
100米网格中手工标记的像素上进行训练预测得到。
[0009]Esri的Land Cover产品：Esri公司2021年发布的产品，空间分辨为10m，时间跨度为2017
‑
2021，覆盖全球范围，使用哨兵2号作为模型输入数据，每期产品使用单时相数据。Esri的产品使用深度学习方法制作，采用的训练数据集与谷歌的产品相同，但在分类结束后根据先验对分类结果做了微调。
[0010]然而，一方面，这些产品的制作过程没有考虑到耕地的时空异质性，另一方面，由于不同数据之间使用的土地覆盖分类体系不同，分类方法不同，验证方法不同，不同产品之间在协同融合时必然会出现较大不一致性，不同的产品直接叠加可能将不同产品的错误分类结果带到融合产品从而导致融合精度降低。
[0011]产品融合技术方法：产品融合方法需要克服现有土地覆盖产品自身缺陷以及不同产品之间的不协调、不兼容的限制，才能融合得到更高精度的融合结果。当前的产品融合方法可以从学习式和非学习式两类产品融合方法来分类介绍。
[0012]非学习式：非学习式方法试图结合统计的数据融合方法来解决现有产品之间兼容性、可比性和准确性不够的问题。按照融合方法这些研究可以分为基于地理加权回归(GWR)的和基于融合决策规则的。GWR方法其特征是回归参数随着空间距离的变化而变化，因此能够较好的反应空间位置因素带来的影响。Linda See等人在非专利文献“Building a hybrid land cover map with crowdsourcing and geographically weighted regression，ISPRS Journal of Photogrammetry and Remote Sensing,103:48
‑
56，2015”通过GWR方法将GLC2000,MODIS全球土地覆盖产品，GlobCover三个产品融合，生成了全球300m产品图。
[0013]融合决策规则包含有多种模型，贝叶斯理论，Dempster
‑
Shafer证据理论。贝叶斯理论在已知先验概率分布的条件下可提高类别之间的可分性，A.P
é
rez
‑
Hoyos等在非专利文献：Integrating multiple land cover maps through amulti
‑
criteria analysis to improve agricultural monitoring in Africa,International Journal of Applied Earth Observation and Geoinformation，2020”中使用贝叶斯定理制作了融合的全球土地覆盖地图，给出了每个像素分类的置信度值，减少了因数据协调性和空间采样等因素产生的误差。
[0014]学习式：学习式的产品融合方法试图收集到不同产品高可信度的样本，将其作为标签输入到分类器中，并通过额外的遥感影像重新生成一幅精度更高的局部或全球的地图。由于产品之间质量参差不齐，通常需要一个样本筛选机制来选取可靠样本，常用的方法有基于规则筛选，如只保留多年数据均没有发生变化的样本，只保留中心像素以减少边缘效应的影响，或使用边缘缩减算子以去除边界像素，Liu等人在非专利文献“High
‑
resolution multi
‑
temporal mapping of global urban land usi本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.弱监督时序特征驱动的遥感耕地产品融合制图方法，其特征在于，包括以下步骤：弱监督时序数据集构建；构造由弱监督损失改造后的utae时序网络模型；弱监督模型训练，得到初步的模型后再针对难样本进行补充训练；最后使用训练好的模型进行适当范围的耕地覆盖制图。2.根据权利要求1所述的弱监督时序特征驱动的遥感耕地产品融合制图方法，其特征在于，所述弱监督时序数据集中的弱监督标签制作过程如下：对给定的三个数据源产品P
x
,P
y
,P
z
进行标签映射，将产品中除耕地外的地类转化为相同的背景像元，保留耕地像元，得到变换后的耕地专题产品P
x
′
,P
y
′
,P
z
′
；将三个产品进行栅格相加，得到栅格值范围0
‑
3的单一栅格产品A，对栅格A做标签映射，将栅格值为1的像素映射至2，栅格值为3的像素映射至1，得到最终弱监督标签L
w
，以上弱监督标签中栅格值为2的部分不参与模型的训练，栅格值为1的部分为耕地像素，0为背景像素；该步骤的公式化描述如下：L
w
＝f
″
(Sum(f
′
(P
x
),f
′
(P
y
),f
′
(P
z
)))其中P
x
,P
y
,P
z
代表原始产品，f'为第一步标签映射规则，Sum代表栅格相加，f
″
为第二步标签映射规则。3.根据权利要求1所述的弱监督时序特征驱动的遥感耕地产品融合制图方法，其特征在于，所述构造由弱监督损失改造后的utae时序网络模型，包括：步骤201：获取utae解码层最后一层的输出，即网络softmax层前一层的经过解码后的影像特征向量V
f
，此时特征向量尺寸为C*W*H，其中W,H为空间维，C为通道维；在空间维度上，以1％的比例随机选取特征像素获得数量为0.01*W*H的特征集合S＝{p1,p2,p3...p
n
}，n＝0.01
×
W
×
H，其中每个特征的尺寸为C*1*1；步骤202：计算空间维度相关性：对于步骤201中得到的特征集合S，对其中的每个元素做如下操作：取该特征p
i
在空间维度上相邻的8个像素，若处于图像边界上则只有3
‑
5个像素；计算该元素与其空间维度上相邻像元特征的欧式距离，并选取距离最近和距离最远的两个结果作为与该像素距离最近和最远的两个特征向量；步骤203：计算特征维度相关性：对于步骤201中得到的特征集合S，对其中的每个元素做如下操作：取该特征p
i
在通道维度上多个特征向量V'构成一个新的通道特征集合S
v
，其中每个向量的尺寸为1*1*1，依次计算S
v
中每个元素之间的余弦相似度，取其中的最小值作为通道特征上的最小相似度C
min
；步骤204：构造弱监督损失：对特征集合S中每个元素在步骤202与步骤203中对所做计算得到的每个元素在空间上最小和最大空间相似度距离d
min
和d
max
，在通道维度上计算得到的最小特征相似度C
min
三个特征值做求和操作，得到弱监督损失：其中α,β,γ分别为各个特征相似度的权重参数，用于调整各个相似度对最终损失的贡
献；步骤205：构造最终模型的完整损失：只在标签值0和1...

【专利技术属性】
技术研发人员：陶超，胡傲然，王宇泽，李海峰，
申请(专利权)人：中南大学，
类型：发明
国别省市：