县域耕地质量评价成果异常值检验方法
【技术领域】
[0001] 本发明涉及耕地质量调查评价领域，尤其涉及一种县域耕地质量评价成果异常值 检验方法。
【背景技术】
[0002] 随着我国耕地保护工作的加强，耕地质量评价已经成为每年的常态性工作。耕地 质量评价成果在耕地综合生产能力核算、耕地占补平衡考核、土地利用总体规划修编以及 土地综合整治等土地管理方面发挥了重要的作用。因此，对耕地质量评价成果的合理性、准 确性进行检验具有十分重要的现实意义。
[0003] 在现有技术中，对耕地质量评价成果的异常值进行检验主要通过重现计算过程， 较为常见的异常值检验方法是空间自相关性分析方法。但目前所采用的空间自相关技术对 异常值的选取仅仅依据单一变量(距离或相邻关系)进行判定，而未能考虑耕地面积对耕地 图斑的影响，导致检验结果的合理性、准确性较低。
[0004] 因此，如何针对上述缺陷提供一种合理、准确的县域耕地质量评价成果异常值检 验方法成为目前急需解决的技术问题之一。
【发明内容】
[0005] 为了解决上述技术问题，本发明提出了一种县域耕地质量评价成果异常值检验方 法，该方法包括：
[0006] 根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立 所述不同预设阈值距离下的空间权重矩阵；
[0007] 根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果 图的疑似异常值图斑；
[0008] 统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确 定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑 或非确定异常值图斑。
[0009] 优选地，所述根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈 值距离建立所述不同预设阈值距离下的空间权重矩阵，包括：
[0010] 根据空间权重矩阵的元素计算公式建立所述不同预设阈值距离下的空间权重矩 阵；
[0011] 所述空间权重矩阵的元素计算公式为
[0013]其中，为所述空间权重矩阵W的元素，表示第j块耕地图斑对第i块耕地图斑的影 响权重，为在当前预设阈值距离内所述第j块耕地图斑对第i块耕地图斑的空间权重；
，为所述第j块耕地图斑对所述第i块耕地图斑面积归一化 后的值，maxdPmim分别为在所述第i块耕地图斑的所述当前预设阈值距离内的耕地图斑面 积的最大值和最小值;M」为所述第i块耕地图斑的所述当前预设阈值距离内所述第j块耕地 的面积。
[0014] 优选地，所述根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质 量评价成果图的疑似异常值图斑，包括：
[0015] 根据所述不同预设阈值距离，确定使得所述每块耕地图斑均具有至少一块相邻耕 地图斑的最小预设阈值距离U，以根据该最小预设阈值距离1^建立该最小预设阈值距离U 下的空间权重矩阵；
[0016] 根据该最小预设阈值距离。下的空间权重矩阵确定所述每块耕地图斑的局部莫 兰指数，以根据该局部莫兰指数提取所述县域耕地质量评价成果图的疑似异常值图斑。
[0017] 优选地，所述确定所述每块耕地图斑的局部莫兰指数的公式为
[0019] 其中，^为所述第i块耕地图斑的局部莫兰指数，XdPXj分别为所述第i块耕地图斑 和所述第j块耕地图斑的自然等指数值，:?为所述第i块耕地图斑的自然等指数值的均值，
^表示所述第i块耕地图斑的自然等指数值的方差，η为所述县域耕 地质量评价成果图中的图斑总数；
[0020] 若所述的值为负值时，确定所述第i块耕地图斑为疑似异常值图斑。
[0021] 优选地，所述根据所述最小预设阈值距离1^下的空间权重矩阵确定所述每块耕地 图斑的局部莫兰指数，还包括：
[0022] 利用Z检验方法检验所述每块耕地图斑的局部莫兰指数的置信度，以根据Z检验结 果确定满足预设置信度要求的局部莫兰指数。
[0023] 优选地，所述统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等 指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确 定异常值图斑或非确定异常值图斑，包括：
[0024]根据所述Z检验结果中所述每块疑似异常值图斑的Z得分以及所述疑似异常值图 斑的局部莫兰指数确定所述每块疑似异常值图斑的聚集类型；
[0025]根据所述不同预设阈值距离下的增量空间自相关确定所述不同预设阈值距离下 的莫兰指数，并根据该莫兰指数确定缓冲区距离；
[0026]根据所述缓冲区距离从每个所述疑似异常值图斑的边界分别向外侧建立所述每 个所述疑似异常值图斑的缓冲区；
[0027]根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区 内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指 数的差值，并将该差值和所述每块耕地图斑的自然等指数标准差进行比较；
[0028] 根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
[0029] 优选地，所述根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑 对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图 斑的自然等指数的差值，包括：
[0030] 若所述疑似异常值图斑的聚集类型为高-低聚集类型，则获取该疑似异常值图斑 对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最大值与该疑似异常值图斑的 自然等指数的差值；
[0031] 若所述疑似异常值图斑的聚集类型为低-高聚集类型，则获取该疑似异常值图斑 对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最小值与该疑似异常值图斑的 自然等指数的差值。
[0032] 优选地，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定 异常值图斑，包括：
[0033]若所述差值大于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图 斑为确定异常值图斑；
[0034]若所述差值小于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图 斑为非确定异常值图斑。
[0035]优选地，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定 异常值图斑之后，该方法还包括：
[0036]判断每一块所述疑似异常值图斑的缓冲区中所述疑似异常值图斑的数量；
[0037] 当所述疑似异常值图斑的缓冲区中存在另外一块疑似异常值图斑，且当所述另外 一块疑似异常值图斑已确定为非确定异常值图斑时，在确定该缓冲区所属的疑似异常值图 斑是确定异常值图斑或非确定异常值图斑时，将该块已确定的非确定异常值图斑加入到所 述全部非疑似异常值图斑中；
[0038] 当所述疑似异常值图斑的缓冲区中存在至少4块所述疑似异常值图斑时，将该缓 冲区所述的疑似异常值图斑确定为非确定疑似异常值图斑。
[0039] 本发明的县域耕地质量评价成果异常值检验方法操作简单、易于实现，可合理、准 确地对县域耕地质量评价成果异常值进行检验。
【附图说明】
[0040] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图