[发明专利]一种面向数字土壤制图的地形协同变量选取方法

权利要求书

1.一种面向数字土壤制图的地形协同变量选取方法，其特征在于，包括如下步骤：

步骤01.针对各个土壤样本点对应的各个地形因子变量V_i,j进行预处理，使其符合正态分布；其中，1≤i≤m，1≤j≤n_i，m为地形因子的种类数，n_i为对应第i种地形因子的计算方法的数量；

步骤02.针对经步骤01处理后的各个地形因子变量V_i,j与各个土壤样本点的土壤属性数据S_long,lati分别进行单地形因子多算法的皮尔逊相关性分析，分别获得各种地形因子对应其各种计算方法与该土壤样本点的土壤属性数据ΣS_long,lati之间的相关性结果pr_i,j，并分别构成各种地形因子对应其各种计算方法与土壤属性数据ΣS_long,lati之间的相关性集合pri=(pri,1,···,pri,j,···,pri,ni);]]>

步骤03.分别针对各种地形因子对应其各种计算方法与土壤属性数据ΣS_long,lati之间的相关性集合中的各个相关性结果pr_i,j，分别采用如下公式（1）进行评价：

LDTASETi,j=|λ×pri,jTi,j|×ξ---(1)]]>

针对各个pr_i集合，将各个集合中的各LDTASET_i,j按大小降序排列，由大至小取出前u个LDTASET_i,j，分别获得各种地形因子分别对应其u种最优计算方法的地形因子变量V′_i,j，并分别构成各种地形因子对应其最优计算方法的地形因子变量集合V_i；其中，LDTASET_i,j为各种地形因子对应其各种计算方法的评价结果；T_i,j表示对应第i种地形因子的第j种计算方法的计算时间；ξ表示地形因子栅格值的分布是否符合正态分布，符合ξ为1，否则ξ为0；λ为膨胀系数，为预设计算精度控制系数；u为预设选取地形因子对应其最优算法的个数，且u≤min(n_i)；

步骤04.针对中各种地形因子对应其最优计算方法的地形因子变量集合V_i中的各个地形因子变量V′_i,j进行线性变化，并进行降维操作；

步骤05.针对步骤04中进行降维操作后获得的各种地形因子对应其最优计算方法的地形因子变量集合V_i中的各个地形因子变量V′_i,j与各个土壤样本点的土壤属性数据S_long,lati分别进行多地形因子单算法的皮尔逊相关性分析，分别获得各种地形因子对应其最优计算方法与土壤属性数据ΣS_long,lati之间的相关性结果pr′_i,j；

步骤06.针对步骤05中获得的各种地形因子对应其最优计算方法与土壤属性数据ΣS_long,lati之间的相关性结果pr′_i,j，分别采用如下公式（2）进行评价：

MDTASETi,j=|pri,j′Ti,j|---(2)]]>

其中，MDTASET_i,j为该土壤样本中各种地形因子对应其最优计算方法的评价结果；

步骤07.根据降维操作后获得的地形因子的种类数c，以及其各自对应的最优的计算方法的数量u，获得(c×u)种地形因子变量组合，针对该(c×u)种地形因子变量组合，根据如下公式（3）：

MLDTAk=Σk=1c×m|prk′-prk-1′|MDTASETk-MDTASETk-1---(3)]]>

获得各种地形因子变量组合的评价结果MLDTA_k，1≤k≤(c×u)；取出MLDTA_k的最大值，根据对应地形因子变量组合获取地形协同变量。