主权项 |
一种面向数字土壤制图的地形协同变量选取方法,其特征在于,包括如下步骤:步骤01.针对各个土壤样本点对应的各个地形因子变量V<sub>i,j</sub>进行预处理,使其符合正态分布;其中,1≤i≤m,1≤j≤n<sub>i</sub>,m为地形因子的种类数,n<sub>i</sub>为对应第i种地形因子的计算方法的数量;步骤02.针对经步骤01处理后的各个地形因子变量V<sub>i,j</sub>与各个土壤样本点的土壤属性数据S<sub>long,lati</sub>分别进行单地形因子多算法的皮尔逊相关性分析,分别获得各种地形因子对应其各种计算方法与该土壤样本点的土壤属性数据ΣS<sub>long,lati</sub>之间的相关性结果pr<sub>i,j</sub>,并分别构成各种地形因子对应其各种计算方法与土壤属性数据ΣS<sub>long,lati</sub>之间的相关性集合<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>pr</mi><mi>i</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>pr</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><msub><mi>pr</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><msub><mi>pr</mi><msub><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000493690910000012.GIF" wi="598" he="88" /></maths>步骤03.分别针对各种地形因子对应其各种计算方法与土壤属性数据ΣS<sub>long,lati</sub>之间的相关性集合<img file="FDA0000493690910000013.GIF" wi="578" he="84" />中的各个相关性结果pr<sub>i,j</sub>,分别采用如下公式(1)进行评价:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>LDTASET</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mo>|</mo><mi>λ</mi><mo>×</mo><mfrac><msub><mi>pr</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mfrac><mo>|</mo><mo>×</mo><mi>ξ</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000493690910000011.GIF" wi="1249" he="175" /></maths>针对各个pr<sub>i</sub>集合,将各个集合中的各LDTASET<sub>i,j</sub>按大小降序排列,由大至小取出前u个LDTASET<sub>i,j</sub>,分别获得各种地形因子分别对应其u种最优计算方法的地形因子变量V′<sub>i,j</sub>,并分别构成各种地形因子对应其最优计算方法的地形因子变量集合V<sub>i</sub>;其中,LDTASET<sub>i,j</sub>为各种地形因子对应其各种计算方法的评价结果;T<sub>i,j</sub>表示对应第i种地形因子的第j种计算方法的计算时间;ξ表示地形因子栅格值的分布是否符合正态分布,符合ξ为1,否则ξ为0;λ为膨胀系数,为预设计算精度控制系数;u为预设选取地形因子对应其最优算法的个数,且u≤min(n<sub>i</sub>);步骤04.针对中各种地形因子对应其最优计算方法的地形因子变量集合V<sub>i</sub>中的各个地形因子变量V′<sub>i,j</sub>进行线性变化,并进行降维操作;步骤05.针对步骤04中进行降维操作后获得的各种地形因子对应其最优计算方法的地形因子变量集合V<sub>i</sub>中的各个地形因子变量V′<sub>i,j</sub>与各个土壤样本点的土壤属性数据S<sub>long,lati</sub>分别进行多地形因子单算法的皮尔逊相关性分析,分别获得各种地形因子对应其最优计算方法与土壤属性数据ΣS<sub>long,lati</sub>之间的相关性结果pr′<sub>i,j</sub>;步骤06.针对步骤05中获得的各种地形因子对应其最优计算方法与土壤属性数据ΣS<sub>long,lati</sub>之间的相关性结果pr′<sub>i,j</sub>,分别采用如下公式(2)进行评价:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>MDTASET</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mo>|</mo><mfrac><msubsup><mi>pr</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>′</mo></msubsup><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mfrac><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000493690910000021.GIF" wi="1097" he="175" /></maths>其中,MDTASET<sub>i,j</sub>为该土壤样本中各种地形因子对应其最优计算方法的评价结果;步骤07.根据降维操作后获得的地形因子的种类数c,以及其各自对应的最优的计算方法的数量u,获得(c×u)种地形因子变量组合,针对该(c×u)种地形因子变量组合,根据如下公式(3):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>MLDTA</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>c</mi><mo>×</mo><mi>m</mi></mrow></munderover><mo>|</mo><msubsup><mi>pr</mi><mi>k</mi><mo>′</mo></msubsup><mo>-</mo><msubsup><mi>pr</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>′</mo></msubsup><mo>|</mo></mrow><mrow><msub><mi>MDTASET</mi><mi>k</mi></msub><mo>-</mo><msub><mi>MDTASET</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000493690910000022.GIF" wi="1310" he="206" /></maths>获得各种地形因子变量组合的评价结果MLDTA<sub>k</sub>,1≤k≤(c×u);取出MLDTA<sub>k</sub>的最大值,根据对应地形因子变量组合获取地形协同变量。 |