CN105550938B - 县域耕地质量评价成果异常值检验方法 - Google Patents

县域耕地质量评价成果异常值检验方法 Download PDF

Info

Publication number
CN105550938B
CN105550938B CN201511031521.4A CN201511031521A CN105550938B CN 105550938 B CN105550938 B CN 105550938B CN 201511031521 A CN201511031521 A CN 201511031521A CN 105550938 B CN105550938 B CN 105550938B
Authority
CN
China
Prior art keywords
abnormal value
spot
cultivated land
pattern spot
pattern
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201511031521.4A
Other languages
English (en)
Other versions
CN105550938A (zh
Inventor
杨永侠
王旭
施彦如
孙婷
孟丹
朱德海
杨建宇
汪念
岳彦利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Agricultural University
Original Assignee
China Agricultural University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Agricultural University filed Critical China Agricultural University
Priority to CN201511031521.4A priority Critical patent/CN105550938B/zh
Publication of CN105550938A publication Critical patent/CN105550938A/zh
Application granted granted Critical
Publication of CN105550938B publication Critical patent/CN105550938B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Husbandry (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种县域耕地质量评价成果异常值检验方法，包括根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵；根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑；统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑。本发明的县域耕地质量评价成果异常值检验方法操作简单、易于实现，可快速、准确地对县域耕地质量评价成果异常值进行检验。

Description

县域耕地质量评价成果异常值检验方法
技术领域
本发明涉及耕地质量调查评价领域，尤其涉及一种县域耕地质量评价成果异常值检验方法。
背景技术
随着我国耕地保护工作的加强，耕地质量评价已经成为每年的常态性工作。耕地质量评价成果在耕地综合生产能力核算、耕地占补平衡考核、土地利用总体规划修编以及土地综合整治等土地管理方面发挥了重要的作用。因此，对耕地质量评价成果的合理性、准确性进行检验具有十分重要的现实意义。
在现有技术中，对耕地质量评价成果的异常值进行检验主要通过重现计算过程，较为常见的异常值检验方法是空间自相关性分析方法。但目前所采用的空间自相关技术对异常值的选取仅仅依据单一变量(距离或相邻关系)进行判定，而未能考虑耕地面积对耕地图斑的影响，导致检验结果的合理性、准确性较低。
因此，如何针对上述缺陷提供一种合理、准确的县域耕地质量评价成果异常值检验方法成为目前急需解决的技术问题之一。
发明内容
为了解决上述技术问题，本发明提出了一种县域耕地质量评价成果异常值检验方法，该方法包括：
根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵；
根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑；
统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
优选地，所述根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵，包括：
根据空间权重矩阵的元素计算公式建立所述不同预设阈值距离下的空间权重矩阵；
所述空间权重矩阵的元素计算公式为
Figure BDA0000898477240000021
其中，Wij为所述空间权重矩阵W的元素，表示第j块耕地图斑对第i块耕地图斑的影响权重，Aij为在当前预设阈值距离内所述第j块耕地图斑对第i块耕地图斑的空间权重；
Figure BDA0000898477240000022
为所述第j块耕地图斑对所述第i块耕地图斑面积归一化后的值，maxi和mini分别为在所述第i块耕地图斑的所述当前预设阈值距离内的耕地图斑面积的最大值和最小值；Mj为所述第i块耕地图斑的所述当前预设阈值距离内所述第j块耕地的面积。
优选地，所述根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑，包括：
根据所述不同预设阈值距离，确定使得所述每块耕地图斑均具有至少一块相邻耕地图斑的最小预设阈值距离L1，以根据该最小预设阈值距离L1建立该最小预设阈值距离L1下的空间权重矩阵；
根据该最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，以根据该局部莫兰指数提取所述县域耕地质量评价成果图的疑似异常值图斑。
优选地，所述确定所述每块耕地图斑的局部莫兰指数的公式为
Figure BDA0000898477240000031
其中，Ii为所述第i块耕地图斑的局部莫兰指数，xi和xj分别为所述第i块耕地图斑和所述第j块耕地图斑的自然等指数值，
Figure BDA0000898477240000032
为所述第i块耕地图斑的自然等指数值的均值，
Figure BDA0000898477240000033
表示所述第i块耕地图斑的自然等指数值的方差，n为所述县域耕地质量评价成果图中的图斑总数；
若所述Ii的值为负值时，确定所述第i块耕地图斑为疑似异常值图斑。
优选地，所述根据所述最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，还包括：
利用Z检验方法检验所述每块耕地图斑的局部莫兰指数的置信度，以根据Z检验结果确定满足预设置信度要求的局部莫兰指数。
优选地，所述统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑，包括：
根据所述Z检验结果中所述每块疑似异常值图斑的Z得分以及所述疑似异常值图斑的局部莫兰指数确定所述每块疑似异常值图斑的聚集类型；
根据所述不同预设阈值距离下的增量空间自相关确定所述不同预设阈值距离下的莫兰指数，并根据该莫兰指数确定缓冲区距离；
根据所述缓冲区距离从每个所述疑似异常值图斑的边界分别向外侧建立所述每个所述疑似异常值图斑的缓冲区；
根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，并将该差值和所述每块耕地图斑的自然等指数标准差进行比较；
根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
优选地，所述根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，包括：
若所述疑似异常值图斑的聚集类型为高-低聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最大值与该疑似异常值图斑的自然等指数的差值；
若所述疑似异常值图斑的聚集类型为低-高聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最小值与该疑似异常值图斑的自然等指数的差值。
优选地，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑，包括：
若所述差值大于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为确定异常值图斑；
若所述差值小于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为非确定异常值图斑。
优选地，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑之后，该方法还包括：
判断每一块所述疑似异常值图斑的缓冲区中所述疑似异常值图斑的数量；
当所述疑似异常值图斑的缓冲区中存在另外一块疑似异常值图斑，且当所述另外一块疑似异常值图斑已确定为非确定异常值图斑时，在确定该缓冲区所属的疑似异常值图斑是确定异常值图斑或非确定异常值图斑时，将该块已确定的非确定异常值图斑加入到所述全部非疑似异常值图斑中；
当所述疑似异常值图斑的缓冲区中存在至少4块所述疑似异常值图斑时，将该缓冲区所述的疑似异常值图斑确定为非确定疑似异常值图斑。
本发明的县域耕地质量评价成果异常值检验方法操作简单、易于实现，可合理、准确地对县域耕地质量评价成果异常值进行检验。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出了本发明一个实施例的县域耕地质量评价成果异常值检验方法流程图；
图2示出了本发明另一实施例的根据异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1示出了本发明一个实施例的县域耕地质量评价成果异常值检验方法流程图。如图1所示，该方法包括：
S1：根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵；
S2：根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑；
S3：统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
本实施例的县域耕地质量评价成果异常值检验方法操作简单、易于实现，可合理、准确地对县域耕地质量评价成果异常值进行检验。
在上述实施例的基础上，步骤S1优选地还包括：
S11：根据空间权重矩阵的元素计算公式建立所述不同预设阈值距离下的空间权重矩阵；
所述空间权重矩阵的元素计算公式为
Figure BDA0000898477240000061
其中，Wij为所述空间权重矩阵W的元素，表示第j块耕地图斑对第i块耕地图斑的影响权重，Aij为在当前预设阈值距离内所述运用反距离权重法建立的第j块耕地图斑对第i块耕地图斑的空间权重；
Figure BDA0000898477240000062
为所述第j块耕地图斑对所述第i块耕地图斑面积归一化后的值，maxi和mini分别为在所述第i块耕地图斑的所述当前预设阈值距离内的耕地图斑面积的最大值和最小值；Mj为所述第i块耕地图斑的所述当前预设阈值距离内所述第j块耕地的面积优选地，建立上述权重矩阵的时候，所有的预设阈值距离都是以面状耕地单元的质心为标准建立的。
作为本实施例的优选，步骤S2可包括：
S21：根据所述不同预设阈值距离，确定使得所述每块耕地图斑均具有至少一块相邻耕地图斑的最小预设阈值距离L1，以根据该最小预设阈值距离L1建立该最小预设阈值距离L1下的空间权重矩阵；
S22：根据该最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，以根据该局部莫兰指数提取所述县域耕地质量评价成果图的疑似异常值图斑。
进一步地，步骤S22中确定所述每块耕地图斑的局部莫兰指数的公式为
Figure BDA0000898477240000071
其中，Wij为当前预设阈值距离(最小预设阈值距离L1)下的空间权重矩阵的元素，Ii为所述第i块耕地图斑的局部莫兰指数，xi和xj分别为所述第i块耕地图斑和所述第j块耕地图斑的自然等指数值，
Figure BDA0000898477240000072
为所述第i块耕地图斑的自然等指数值的均值，
Figure BDA0000898477240000073
表示所述第i块耕地图斑的自然等指数值的方差，n为所述县域耕地质量评价成果图中的图斑总数；
在此基础上，步骤S22中所述根据所述最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，还可包括：
S221：利用Z检验方法检验所述每块耕地图斑的局部莫兰指数的置信度，以根据Ii的值和Z检验结果确定疑似异常值图斑。
具体地，考虑耕地的复杂性，为了更加准确，可以使用|Z|大于1.65的置信度为90％去选择疑似异常值。即，若所述Ii的值为负值、Z检验的绝对值大于1.65时，确定所述第i块耕地图斑为疑似异常值图斑。
图2示出了本发明另一实施例的根据异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑的流程图。如图2所示，步骤S3可包括：
S31：根据所述Z检验结果中所述每块疑似异常值图斑的Z得分以及所述疑似异常值图斑的局部莫兰指数确定所述每块疑似异常值图斑的聚集类型；
S32：根据所述不同预设阈值距离下的增量空间自相关确定所述不同预设阈值距离下的莫兰指数，并根据该莫兰指数确定缓冲区距离；根据该莫兰指数确定缓冲区距离即以莫兰指数的最高值所对应的阈值距离为缓冲区距离。
S33：根据所述缓冲区距离从每个所述疑似异常值图斑的边界分别向外侧建立所述每个所述疑似异常值图斑的缓冲区；
S34：根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，并将该差值和所述每块耕地图斑的自然等指数标准差进行比较；
S35：根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
进一步地，S34中根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，可包括：
S341：若所述疑似异常值图斑的聚集类型为高-低聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最大值与该疑似异常值图斑的自然等指数的差值；
S342：若所述疑似异常值图斑的聚集类型为低-高聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最小值与该疑似异常值图斑的自然等指数的差值。
作为本实施例的优选，步骤S35可包括：
S351：若所述差值大于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为确定异常值图斑；
S352：若所述差值小于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为非确定异常值图斑。
进一步地，S35之后，该方法还包括：
S36：判断每一块所述疑似异常值图斑的缓冲区中所述疑似异常值图斑的数量；
S361：当所述疑似异常值图斑的缓冲区中存在另外一块疑似异常值图斑，且当所述另外一块疑似异常值图斑已确定为非确定异常值图斑时，在确定该缓冲区所属的疑似异常值图斑是确定异常值图斑或非确定异常值图斑时，将该块已确定的非确定异常值图斑加入到所述全部非疑似异常值图斑中；
S362：当所述疑似异常值图斑的缓冲区中存在至少4块所述疑似异常值图斑时，将该缓冲区所述的疑似异常值图斑确定为非确定疑似异常值图斑。
本实施例的县域耕地质量评价成果异常值检验方法，可以快速、准确判定耕地质量等指数异常值，通过综合考虑面积与距离建立空间权重矩阵，并利用空间自相关算法和统计学方法得到的判定标准，可得到合理的判定标准，能够有效避免由于主观判断制定标准所造成的偏误，进而可以更好地反映耕地质量的实际情况。
以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (9)

1.一种县域耕地质量评价成果异常值检验方法，其特征在于，包括：
根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵；
根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑；
统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
2.如权利要求1所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据县域耕地质量评价成果图中每块耕地图斑面积以及不同预设阈值距离建立所述不同预设阈值距离下的空间权重矩阵，包括：
根据空间权重矩阵的元素计算公式建立所述不同预设阈值距离下的空间权重矩阵；
所述空间权重矩阵的元素计算公式为
Figure FDA0002433054660000011
其中，Wij为所述空间权重矩阵W的元素，表示第j块耕地图斑对第i块耕地图斑的影响权重，Aij为在当前预设阈值距离内所述第j块耕地图斑对第i块耕地图斑的空间权重；
Figure FDA0002433054660000012
为所述第j块耕地图斑对所述第i块耕地图斑面积归一化后的值，maxi和mini分别为在所述第i块耕地图斑的所述当前预设阈值距离内的耕地图斑面积的最大值和最小值；Mj为所述第i块耕地图斑的所述当前预设阈值距离内所述第j块耕地的面积。
3.如权利要求2所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据所述不同预设阈值距离下的空间权重矩阵确定所述县域耕地质量评价成果图的疑似异常值图斑，包括：
根据所述不同预设阈值距离，确定使得所述每块耕地图斑均具有至少一块相邻耕地图斑的最小预设阈值距离L1，以根据该最小预设阈值距离L1建立该最小预设阈值距离L1下的空间权重矩阵；
根据该最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，以根据该局部莫兰指数提取所述县域耕地质量评价成果图的疑似异常值图斑。
4.如权利要求3所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述确定所述每块耕地图斑的局部莫兰指数的公式为
Figure FDA0002433054660000021
其中，Ii为所述第i块耕地图斑的局部莫兰指数，xi和xj分别为所述第i块耕地图斑和所述第j块耕地图斑的自然等指数值，
Figure FDA0002433054660000022
为所述第i块耕地图斑的自然等指数值的均值，
Figure FDA0002433054660000023
表示所述第i块耕地图斑的自然等指数值的方差，n为所述县域耕地质量评价成果图中的图斑总数；
若所述Ii的值为负值时，确定所述第i块耕地图斑为疑似异常值图斑。
5.如权利要求4所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据所述最小预设阈值距离L1下的空间权重矩阵确定所述每块耕地图斑的局部莫兰指数，还包括：
利用Z检验方法检验所述每块耕地图斑的局部莫兰指数的置信度，以根据Z检验结果确定满足预设置信度要求的局部莫兰指数。
6.如权利要求5所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述统计所述县域所有耕地图斑的自然等指数标准差，并根据该自然等指数标准差确定异常值判定标准，以根据该异常值判定标准判断所述疑似异常值图斑为确定异常值图斑或非确定异常值图斑，包括：
根据所述Z检验结果中所述每块疑似异常值图斑的Z得分以及所述疑似异常值图斑的局部莫兰指数确定所述每块疑似异常值图斑的聚集类型；
根据所述不同预设阈值距离下的增量空间自相关确定所述不同预设阈值距离下的莫兰指数，并根据该莫兰指数确定缓冲区距离；
根据所述缓冲区距离从每个所述疑似异常值图斑的边界分别向外侧建立所述每个所述疑似异常值图斑的缓冲区；
根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，并将该差值和所述每块耕地图斑的自然等指数标准差进行比较；
根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑。
7.如权利要求6所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据每个所述疑似异常值图斑的聚集类型，获取该疑似异常值图斑对应的缓冲区内的非疑似异常值图斑的自然等指数的最大值或最小值与该疑似异常值图斑的自然等指数的差值，包括：
若所述疑似异常值图斑的聚集类型为高-低聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最大值与该疑似异常值图斑的自然等指数的差值；
若所述疑似异常值图斑的聚集类型为低-高聚集类型，则获取该疑似异常值图斑对应的缓冲区内的全部非疑似异常值图斑的自然等指数的最小值与该疑似异常值图斑的自然等指数的差值。
8.如权利要求6所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑，包括：
若所述差值大于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为确定异常值图斑；
若所述差值小于所述每块耕地图斑的自然等指数标准差，则确定该疑似异常值图斑为非确定异常值图斑。
9.如权利要求6所述的县域耕地质量评价成果异常值检验方法，其特征在于，所述根据比较结果，确定该疑似异常值图斑为确定异常值图斑或非确定异常值图斑之后，该方法还包括：
判断每一块所述疑似异常值图斑的缓冲区中所述疑似异常值图斑的数量；
当所述疑似异常值图斑的缓冲区中存在另外一块疑似异常值图斑，且当所述另外一块疑似异常值图斑已确定为非确定异常值图斑时，在确定该缓冲区所属的疑似异常值图斑是确定异常值图斑或非确定异常值图斑时，将该块已确定的非确定异常值图斑加入到所述全部非疑似异常值图斑中；
当所述疑似异常值图斑的缓冲区中存在至少4块所述疑似异常值图斑时，将该缓冲区所述的疑似异常值图斑确定为非确定疑似异常值图斑。
CN201511031521.4A 2015-12-03 2015-12-31 县域耕地质量评价成果异常值检验方法 Expired - Fee Related CN105550938B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511031521.4A CN105550938B (zh) 2015-12-03 2015-12-31 县域耕地质量评价成果异常值检验方法

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN2015108818658 2015-12-03
CN201510881865 2015-12-03
CN201511031521.4A CN105550938B (zh) 2015-12-03 2015-12-31 县域耕地质量评价成果异常值检验方法

Publications (2)

Publication Number Publication Date
CN105550938A CN105550938A (zh) 2016-05-04
CN105550938B true CN105550938B (zh) 2020-08-11

Family

ID=55830116

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511031521.4A Expired - Fee Related CN105550938B (zh) 2015-12-03 2015-12-31 县域耕地质量评价成果异常值检验方法

Country Status (1)

Country Link
CN (1) CN105550938B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106682695A (zh) * 2016-12-28 2017-05-17 中国农业大学 一种基于支持向量机的县域耕地自然质量评价方法
CN107403004B (zh) * 2017-07-24 2020-07-24 邱超 一种基于地形数据的遥测雨量站点可疑数值检验方法
CN107862978A (zh) * 2017-11-01 2018-03-30 北京珍宝岛中药产业有限公司 一种生成中药材产地质量地图的方法和装置
CN108830438A (zh) * 2018-04-09 2018-11-16 中国农业大学 一种耕地产能指数异常检验方法及系统
CN109409748B (zh) * 2018-10-29 2021-07-20 广东友元国土信息工程有限公司 一种耕地质量评价指标关联性的检核方法及系统
AR121386A1 (es) * 2020-02-21 2022-06-01 Climate Corp Señalización de diferencias operacionales en implementos agrícolas
CN114720665B (zh) * 2022-06-10 2022-09-16 北京市农林科学院智能装备技术研究中心 测土配方施肥土壤全氮异常值检测方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008047387A2 (en) * 2006-10-17 2008-04-24 Speck System Limited A system for land record information management
CN101430735A (zh) * 2008-11-13 2009-05-13 中国农业大学 一种保护性耕作模式选择方法
CN102722577A (zh) * 2012-06-05 2012-10-10 中兴通讯股份有限公司 指标动态权重的确定方法及装置
CN103559389A (zh) * 2013-10-21 2014-02-05 中国农业大学 一种县域耕地等级变化监测样点布设的方法
CN104318066A (zh) * 2014-09-26 2015-01-28 武汉大学 一种自然地表特征表征方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008047387A2 (en) * 2006-10-17 2008-04-24 Speck System Limited A system for land record information management
CN101430735A (zh) * 2008-11-13 2009-05-13 中国农业大学 一种保护性耕作模式选择方法
CN102722577A (zh) * 2012-06-05 2012-10-10 中兴通讯股份有限公司 指标动态权重的确定方法及装置
CN103559389A (zh) * 2013-10-21 2014-02-05 中国农业大学 一种县域耕地等级变化监测样点布设的方法
CN104318066A (zh) * 2014-09-26 2015-01-28 武汉大学 一种自然地表特征表征方法

Also Published As

Publication number Publication date
CN105550938A (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105550938B (zh) 县域耕地质量评价成果异常值检验方法
US7107187B1 (en) Method for modeling system performance
CN104900551B (zh) 一种晶圆质量管控方法和装置
WO2020119372A1 (zh) 场强预测方法、装置、设备及存储介质
CN107798210B (zh) 一种多模型滑坡位移预测方法及其系统
WO2018184304A1 (zh) 一种网元健康状态的检测方法及设备
CN111967717A (zh) 一种基于信息熵值的数据质量评价方法
WO2021103577A1 (zh) 一种基于mdt的室内信号外泄分析方法及系统
CN109378823A (zh) 一种电压暂降水平的综合评估方法
Kock Common structural variation reduction in PLS-SEM: Replacement analytic composites and the one fourth rule
CN108510072A (zh) 一种基于混沌神经网络的河道流量监测数据质量控制方法
WO2018196538A1 (zh) 生产环境安全性的分析方法、装置和系统
CN114813757A (zh) 一种塑胶产品全角度外观检测方法及装置
CN107729293B (zh) 一种基于多元自适应回归的地理空间异常探测方法
CN115542236A (zh) 电能表运行误差估计方法及装置
CN113947504B (zh) 一种基于随机森林法的窃电分析方法及其系统
CN118113574A (zh) 一种服务器健康度评估方法、装置、设备及存储介质
CN110913407A (zh) 重叠覆盖的分析方法及装置
JP6192432B2 (ja) リスク計量システム
CN116932679A (zh) 一种基于gis技术多模式行业气象服务检验应用方法
TWI592810B (zh) 衡量二元資料於時間等級的群聚級別的無母數分析方法
CN108419247B (zh) 一种无线通信基站小区相关性计算方法及电子设备
CN116339262A (zh) 一种基于人工智能的数控加工生产质量监测系统
CN103852562B (zh) 判断样品检测数据可疑值的方法
CN113689526B (zh) 地图中无效区域的划分方法及装置、电子设备

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200811

Termination date: 20211231