[发明专利]一种基于协同克里金插值法的土壤锰元素含量预测方法有效
申请号: | 201711328434.4 | 申请日: | 2017-12-13 |
公开(公告)号: | CN108088981B | 公开(公告)日: | 2021-07-02 |
发明(设计)人: | 赵晋陵;梁栋;姜燕;高文武;黄林生;阮莉敏;徐超;翁士状;张东彦;洪琪 | 申请(专利权)人: | 安徽大学 |
主分类号: | G01N33/24 | 分类号: | G01N33/24 |
代理公司: | 北京劲创知识产权代理事务所(普通合伙) 11589 | 代理人: | 张铁兰 |
地址: | 230601 安徽省合肥市经*** | 国省代码: | 安徽;34 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种基于协同克里金插值法的土壤锰元素含量预测方法,采用协同克里金插值法来预测其含量。由于土壤中锰元素受其它较多因素的影响,所以本发明采用主成分分析法将这些影响因素组合成一个综合因子,进而将这个综合因子作为辅助变量对主变量锰元素进行协同克里金插值。首先收集了北京市房山区20个乡镇果园地、水浇地、菜地、灌溉水田和旱地这五种耕地类型的土壤重金属元素数据;然后使用方差分析法分析了这五种耕地类型对土壤中锰元素含量的影响;最后在不同的耕地类型下,运用协同克里金插值方法预测土壤中锰元素的含量,并根据标准平均值误差、均方根误差、标准化均方根对插值结果的精度进行分析。 | ||
搜索关键词: | 一种 基于 协同 克里金插值法 土壤 元素 含量 预测 方法 | ||
【主权项】:
1.一种基于协同克里金插值法的土壤锰元素含量预测方法,其特征在于,包括以下步骤:步骤1、方差分析为了探讨耕地类型对土壤重金属含量的影响,使用方差分析法研究了北京市房山区果园地、菜地、水浇地、灌溉水田和旱地五种耕地类型对土壤重金属的影响;五种耕地类型为自变量,是分类型变量;重金属锌、铁、铜、锰、硼、硫的含量为因变量,是数值型变量;步骤2、协同克里金插值法使用协同克里金插值法Cokriging预测了土壤重金属锰的含量;为了解决协调克里金插值时土壤重金属锌、铁、锰这三个辅助变量的权重问题,使用主成分分析法PCA将这三个存在相关性的变量通过降维技术转换为一个综合因子;然后利用综合因子作为辅助变量对主变量锰进行协同克里金插值;根据主成分分析原理,将土壤中重金属Zn、Fe、Cu作为一个综合土壤重金属变量,主成分构成公式如下:Fi =a1i Zn+a2i Fe+a3i Cu式中,Fi 为第i主成分变量,aji (j=1,2,3)分别表示为第i主成分对应原始变量Zn、Fe、Cu的系数,度量了相应原始变量对Fi 的重要性;Zn表示Zn的含量,Fe表示Fe的含量,Cu表示Cu的含量;综合因子的构成公式如下:Y=∑λi Fi 式中,Y为综合因子,Fi 为第i主成分,λi 为主成分的贡献率;协同克里金法插值的计算公式如下: Z * ( l , k ) 0 = Σ i = 1 n λ i Z 1 ( l , k ) 1 i + Σ j = 1 m λ j Z 2 ( l , k ) 2 j ]]> 式中,Z* (l,k)0 为(l,k)0 处土壤重金属Mn含量的估值,n和m为主变量和辅助变量的采样点数,Z1 (l,k)1i 为各点土壤重金属Mn的含量,λi 为赋予各点土壤重金属Mn含量的一组权重系数,且∑λi =1;Z2 (l,k)2j 为各点土壤综合重金属含量,λi 为赋予各点综合土壤重金属含量的一组权重系数;半方差函数公式如下: γ 12 ( h l , h k ) = 1 2 E { [ Z 1 ( l + h l , k + h k ) - Z 1 ( l , k ) ] [ Z 2 ( l + h l , k + h k ) - Z 2 ( l , k ) ] } ]]>
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于安徽大学,未经安徽大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201711328434.4/,转载请声明来源钻瓜专利网。