[发明专利]一种土壤指标空间变化的判断方法

说明书

本发明属于土壤探测技术领域。公开了一种土壤指标空间变化的判断方法，包括如下步骤：确定进行变化探测的土壤指标C；将前期观测点数据C1(x，y)使用克里金插值方法得到区域数据，得到插值后的空间数据C1′(x，y)及其标准误的空间分布数据C1′E(x，y)；将后期观测点数据C2(x，y)使用克里金插值方法得到区域数据C2′(x，y)及其标准误的空间分布数据C2′E(x，y)；生成空间数据dC(x，y)，用以测度变化状况；其标准误为：使用下式计算zC(x，y)值，得到其空间数据：zC(x，y)＝dC(x，y)/dCE(x，y)取显著性水平α，“不变”区间的取值依据标准正态分布函数：[‑Z(α/2)，+Z(α/2)]。本发明的方法利用构造的z值，可以在给定显著性水平下，确定土壤指标是否发生了“增加”、“降低”或不变。

技术领域

本发明涉及土壤环境探测技术领域，具体涉及一种土壤指标空间变化的判断方法。

背景技术

本发明对于背景技术的描述属于与本发明相关的相关技术，仅仅是用于说明和便于理解本发明的发明内容，不应理解为申请人明确认为或推定申请人认为是本发明在首次提出申请的申请日的现有技术。

环境指标包括pH值、有机质等土壤基本性质指标，氮、磷等营养元素含量，也可以是镉、汞、砷等有害元素含量，或者碘、硒等健康元素含量；也可能是以上元素的某种化学提取态。区域土壤环境指标变化探测通常使用土壤地球化学调查等手段，通过采集土壤样品送实验室测试，获取土壤环境指标的空间分布状况。然而，当前的土壤环境管理还需要知道指标的变化状况：“增加”，“减少(或降低)”，还是“没变”，以便回答土壤环境状况“继续恶化”、“得到遏制”，还是“稳中向好”等问题。

通常使用重访获取“两期”(或更多期)调查或监测数据判断上述状况。然而，该方法面临以下障碍：观测数据有误差，点位无法严格重合，由“点”到“面”的数据推算带来不确定性。本发明针对该情景给出数据分析方案，克服以上障碍，得到科学结论。

发明内容

本发明实施例的目的是提供一种土壤指标空间变化的判断方法，本发明的方法利用构造的z值，可以在给定显著性水平下，确定土壤指标是否发生了“增加”、“降低”或不变。

本发明实施例的目的是通过如下技术方案实现的：

一种土壤指标空间变化的判断方法，包括如下步骤：

确定进行变化探测的土壤指标C；

将前期观测点数据C1(x,y)使用克里金插值方法得到区域数据，得到插值后的空间数据C1'(x,y)及其标准误的空间分布数据C1'E(x,y)；将后期观测点数据C2(x,y)使用克里金插值方法得到区域数据C2'(x,y)及其标准误的空间分布数据C2'E(x,y)；

生成空间数据dC(x,y)，用以测度变化状况；其标准误为：

使用下式计算zC(x,y)值，得到其空间数据：

zC(x,y)＝dC(x,y)/dCE(x,y)

取显著性水平α，“不变”区间的取值依据标准正态分布函数：[-Z(α/2),+Z(α/2)]。

进一步的，所述的土壤指标空间变化的判断方法，其特征在于，生成空间数据dC(x,y)的方法为使用dC(x,y)＝C2'(x,y)-C1'(x,y)运算生成空间数据dC(x,y)。

进一步的，所述的土壤指标空间变化的判断方法，其特征在于，生成空间数据dC(x,y)的方法为使用d'(x,y)＝C2'(x,y)/C1'(x,y)运算生成空间数据dC(x,y)。

进一步的，进行变化探测的土壤指标C为土壤观测指标或土壤观测指标的变换值。

进一步的，所述的土壤观测指标的变换值为对数变换。