一种广覆盖的土壤重金属遥感混合制图方法技术

本发明专利技术公开了一种广覆盖的土壤重金属遥感混合制图方法，通过获取目标区域内的土壤采样点重金属含量值和其他辅助数据，并进行预处理；并进行插值处理，进一步计算残差；采用三次样条函数插值拟合空间裸土范围内的残差函数，校正步骤5得到的遥感预测重金属结果；最终得到区域内裸土土壤重金属含量的预测值。本发明专利技术的优点在于：通过借助遥感数据获取的地面信息提高局部反演的精度，充分发挥插值方法对在结果异常值的控制，扩大土壤制图的空间范围和提高重金属制图的准确度。高重金属制图的准确度。高重金属制图的准确度。

【技术实现步骤摘要】
一种广覆盖的土壤重金属遥感混合制图方法


[0001]本专利技术涉及遥感和环境领域，具体涉及基于土壤中重金属含量的空间制图方法。

技术介绍

[0002]土壤重金属污染严重影响到国家粮食安全和人体健康，开展土壤重金属调查，实现大范围空间制图对经济社会的稳定发展具有重要意义。
[0003]目前利用地面样本点采样建立反演土壤重金属的空间预测性制图，都是通过已知点推测未知点或其他易于获取的量，建立相关模型，预测难于获取的量，达到制图目的。常见的土壤空间预测性制图方法主要有基于地统计学的空间插值法和基于遥感反演预测制图。
[0004]地统计学方法单纯的考虑了土壤重金属元素的空间变化影响，可以仅通过已有的地面采样点，实现空间制图，这也是早期土壤制图广泛采用的技术。但仅考虑采样点之间的空间关系，在较为复杂环境和采样点较少条件下得到的效果往往具有较大的偏差。后期发展的一些插值方法结合土壤发生学理论，通过易于获得的成土环境因素，如植被、地形及其派生因素，建立多变量的综合模型推断难于测定的土壤属性值，改进插值模型。该方法的运用和土壤形成的特殊性有关，在大范围复杂环境中的精度会受影响。基于遥感数据的反演制图充分利用了遥感影像可快速获取大范围空间的地物光谱信息，为大区域的土壤重金属研究制图提供数据支撑。尤其是高光谱技术的发展，可以实现“图谱合一”，为土壤属性的定量反演提供数据基础。地面采样点数据的质量和数量是两种方法的共同限制因素，通常由于地面采样工作量大，地形复杂和偏僻地区采样难度大等因素，采样点的数量往往成为土壤重金属高精度的反演的主要限制因素。
[0005]综上所述，目前土壤属性的空间预测制图技术中，地统统计学方法发展较为成熟且数据材料较为宽松，但土壤局部环境变异大和在大范围应用中精度受限，需要足够且均匀的地面采样点才能保证结果的精度；基于遥感数据的土壤反演是近些年研究较多的制图方法，大范围应用上有很大的优势，但容易受遥感影像数据本身和环境因素影响，模型的迁移性较差，需要大数据训练建模才能保证结果精度，以上两种单一方法均存在明显的局限性。

技术实现思路

[0006]本次专利技术针对目前采样点数量不足，大范围土壤重金属反演制图精度低等问题，提出一种融合空间插值和残差校正等技术，实现高精度、广覆盖的的遥感土壤重金属制图的方法。
[0007]为实现上述技术目的，本专利技术采用了以下技术方案：
[0008]一种广覆盖的土壤重金属遥感混合制图方法，包括以下步骤：
[0009]步骤1，获取目标区域内的土壤采样点重金属含量值和其他辅助数据，并进行预处理。辅助数据包括高程、经纬度、坡向和坡度数据。可以增加对目标区域的特殊性信息描述。
土壤地面采样点数据信息至少包含土壤采样点的位置和重金属含量值，所有的数据通过空间整合到统一大小的网格空间L中。对以上所有的面状数据进行空间重采样，统一到规则的网格M中，其中包含重金属属性信息共l个网格，空间范围记为L。网格为后续步骤处理的基本单元。
[0010]步骤2，采用协同克里金插值法将点状土壤重金属采样数据集P插值到覆盖整个目标区域的面状重金属栅格数据P
’
。协同克里格插值为已有技术，其原理是利用影响土壤重金属浓度的经纬度、海拔、坡度和坡向因子之间的相关性提高对未知点土壤重金属浓度值的预测精度。采用留一法交叉验证计算包含重金属属性信息数据集P的插值结果的残差数据集ε
L
＝ε1，ε2，...ε
j
...，ε
l
。
[0011]ε
jL
＝P
jL
‑
P
′
jL
[0012]式中，P
jL
为网格j处真实的土壤重金属含量值，P
’
jL
为插值计算得到的网格j土壤重金属含量的预测值。
[0013]步骤3，利用步骤2得到的残差信息提取采样点及其所在位置周围网格插值结果作为土壤重金属的模拟采样数据集。本步骤的目的是尽可能利用插值结果扩大原有土壤样本信息量。对插值残差大的，提取其周围取样距离小，反之则越大，全部网格点得到的模拟采样空间范围记为A。
[0014]步骤4，采用多时相遥感影像提取目标区域的裸土范围。裸土面积的提取方法为已有技术，如决策树分类法等。分别提取n景景影像的空间裸土范围：S1，S2，...S
i
...，S
n
。最终需要制图的空间范围全部影像裸土空间的范围的并集，记为S。
[0015]S＝∪(S1，S2，
…
S
i
…
，S
n
)
[0016]步骤5、6是为对每一景影像i的操作。
[0017]步骤5，通过空间匹配技术建立每景影像中裸土区域S
i
和模拟采样点空间A
i
的空间交集为建模数据集，并建立土壤重金属的反演模型，记其空间为AS
i
。通过遥感反演模型预测裸土区域S
i
以内的土壤重金属含量结果并采用留一交叉验证计算空间A
i
部分的数据残差
[0018]AS
i
＝A
i
∩S
i
[0019][0020]其中为步骤5反演模型得到AS
i
区域网格j的重金属浓度预测值，为步骤3中模拟样本网格的重金属含量值。
[0021]步骤6，采用三次样条函数插值对步骤5AS
i
空间网格中的残差数据ζ
jAS
进行拟合，得到覆盖空间范围S
i
面状残差函数ζ
’
j
。对步骤5中S
i
以内，A
i
范围以外空间B
i
的预测结果P
″
j
进一步校正，记为最终覆盖目标区域S
i
空间的重金属插值结果
[0022][0023][0024]步骤7，对n景影像数据重复步骤5、6，分别得到不同时相的和空间的土壤预测结
果，对每个网格预测得到的全部重金属浓度结果取均值，得到区域内裸土土壤重金属含量的预测值。
[0025]与现有技术相比，本专利技术具有如下有益效果：
[0026]本专利技术通过借助遥感数据获取的地面信息提高局部反演的精度，充分发挥插值方法对在结果异常值的控制，扩大土壤制图的空间范围和提高重金属制图的准确度。
附图说明
[0027]图1为本专利技术的步骤处理流程图。
具体实施方式
[0028]下面对本专利技术实施步骤做进一步描述。应当指出，对于本
的普通技术人员来说，在不脱离本专利技术原理的前提下，还可以对本专利技术进行若干改进和修饰，这些改进和修饰也落入本专利技术权利要求的保护范围内。
[0029]本专利技术技术方案可采用计算机软件实现处理。以下结合图1详述高精度土壤重金属空间预测和制图流程。...

【技术保护点】

【技术特征摘要】
1.一种广覆盖的土壤重金属遥感混合制图方法，其特征在于，包括以下步骤：步骤1，获取目标区域内的土壤采样点重金属含量值和其他辅助数据，并进行预处理；步骤2，采用协同克里金插值法将点状土壤重金属采样数据集P插值到覆盖整个目标区域的面状重金属栅格数据P
’
；步骤3，利用步骤2得到的残差信息提取土壤重金属的模拟采样数据集；步骤4，通过多时相遥感影像提取目标区域内的裸土面积；步骤5，建立每景影像土壤重金属反演模型并预测该景影像裸土区域全部的重金属浓度值，并计算残差；步骤6，采用三次样条函数插值拟合空间裸土范围内的残差函数，校正步骤5得到的遥感预测重金属结果；步骤7，对n景影像数据重复步骤5、6，分别得到不同时相的和空间的土壤预测结果，对每个网格预测得到的全部重金属浓度结果取均值，得到区域内裸土土壤重金属含量的预测值。2.根据权利要求1所述的一种广覆盖的土壤重金属遥感混合制图方法，其特征在于：在步骤3中构建的残差加权的模拟样本数据集。3.根据权利要求2所述的一种广覆盖的土壤重金属遥感混合制图方法，其特征在于：步骤3中构建的残差加权的模拟样本构建过程如下：S01,对步骤2得到残差数据集ε
L
进行归...

【专利技术属性】
技术研发人员：蒋曼，梁建，刘银年，孙德新，尹忠海，
申请(专利权)人：南通智能感知研究院，
类型：发明
国别省市：