一种广覆盖的土壤重金属遥感混合制图方法与流程

技术特征：
1.一种广覆盖的土壤重金属遥感混合制图方法，其特征在于，包括以下步骤：步骤1，获取目标区域内的土壤采样点重金属含量值和其他辅助数据，并进行预处理；步骤2，采用协同克里金插值法将点状土壤重金属采样数据集p插值到覆盖整个目标区域的面状重金属栅格数据p’；步骤3，利用步骤2得到的残差信息提取土壤重金属的模拟采样数据集；步骤4，通过多时相遥感影像提取目标区域内的裸土面积；步骤5，建立每景影像土壤重金属反演模型并预测该景影像裸土区域全部的重金属浓度值，并计算残差；步骤6，采用三次样条函数插值拟合空间裸土范围内的残差函数，校正步骤5得到的遥感预测重金属结果；步骤7，对n景影像数据重复步骤5、6，分别得到不同时相的和空间的土壤预测结果，对每个网格预测得到的全部重金属浓度结果取均值，得到区域内裸土土壤重金属含量的预测值。2.根据权利要求1所述的一种广覆盖的土壤重金属遥感混合制图方法，其特征在于：在步骤3中构建的残差加权的模拟样本数据集。3.根据权利要求2所述的一种广覆盖的土壤重金属遥感混合制图方法，其特征在于：步骤3中构建的残差加权的模拟样本构建过程如下：s01,对步骤2得到残差数据集ε
l
进行归一化，利用归一化的结果对设定基准取样距离d进行反残差加权，将各插值点所在的中心网格d
jl
距离半径内的插值数据作为模拟采样点，全部网格点得到的模拟采样空间范围记为a；所述模拟采样点包含插值反演的重金属浓度信息和位置信息；残差归一化公式如下为：其中ε
mah
和ε
min
分别为残差集ε的最大和最小值，残差加权公式为：其中d为基准取样距离。4.根据权利要求1所述的一种广覆盖的土壤重金属遥感混合制图方法，其特征在于：在步骤6采用三次样条函数插值对步骤5as
i
空间网格中的残差数据ζ
jas
进行拟合，得到覆盖空间范围s
i
面状残差函数ζ
′
j
，对步骤5中s
i
以内，a
i
范围以外空间b
i
的预测结果p
″
j
进一步校正，记为最终覆盖目标区域s
i
空间的重金属插值结果其中拟合公式如下：其中拟合公式如下：

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


技术研发人员：蒋曼 梁建 刘银年 孙德新 尹忠海
受保护的技术使用者：南通智能感知研究院
技术研发日：2021.07.26
技术公布日：2021/11/4