CN114720665A - 测土配方施肥土壤全氮异常值检测方法及装置 - Google Patents
测土配方施肥土壤全氮异常值检测方法及装置 Download PDFInfo
- Publication number
- CN114720665A CN114720665A CN202210650110.7A CN202210650110A CN114720665A CN 114720665 A CN114720665 A CN 114720665A CN 202210650110 A CN202210650110 A CN 202210650110A CN 114720665 A CN114720665 A CN 114720665A
- Authority
- CN
- China
- Prior art keywords
- total nitrogen
- soil
- data
- abnormal
- data set
- 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.)
- Granted
Links
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 645
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 324
- 239000002689 soil Substances 0.000 title claims abstract description 316
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 222
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000012360 testing method Methods 0.000 title claims abstract description 71
- 230000004720 fertilization Effects 0.000 title claims abstract description 60
- 238000001514 detection method Methods 0.000 claims abstract description 138
- 238000007619 statistical method Methods 0.000 claims abstract description 17
- 230000005856 abnormality Effects 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 9
- 238000007726 management method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 230000007613 environmental effect Effects 0.000 claims description 6
- 229910052698 phosphorus Inorganic materials 0.000 claims description 6
- 239000011574 phosphorus Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 5
- 239000011591 potassium Substances 0.000 claims description 5
- 229910052700 potassium Inorganic materials 0.000 claims description 5
- 238000013450 outlier detection Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000002085 persistent effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000012271 agricultural production Methods 0.000 claims 1
- 238000009472 formulation Methods 0.000 claims 1
- 239000000203 mixture Substances 0.000 claims 1
- 238000010801 machine learning Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000002028 Biomass Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000007418 data mining Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012732 spatial analysis Methods 0.000 description 3
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 229910002651 NO3 Inorganic materials 0.000 description 1
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000000618 nitrogen fertilizer Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000013107 unsupervised machine learning method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/245—Earth materials for agricultural purposes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2465—Query processing support for facilitating data mining operations in structured databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Medicinal Chemistry (AREA)
- Evolutionary Computation (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Food Science & Technology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computational Linguistics (AREA)
- Medical Informatics (AREA)
- Probability & Statistics with Applications (AREA)
- Fuzzy Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种测土配方施肥土壤全氮异常值检测方法及装置,涉及作物种植领域。该方法包括:获取待检测土壤全氮数据集;根据土壤全氮数据集,基于孤立森林检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;根据土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;根据第一全氮数据集和第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。该方法能够充分地挖掘数据在不同维度、不同层次上的异常性,有效地避免单一方法下造成的漏判、误判问题,极大地提升测土配方施肥土壤全氮数据的异常值检测准确度。
Description
技术领域
本发明涉及作物种植领域,尤其涉及一种测土配方施肥土壤全氮异常值检测方法及装置。
背景技术
测土配方施肥是破解农田肥料利用率偏低、农业面源污染的关键技术措施。测土配方施肥已经累积了数千万条的土壤测试数据,形成了规模庞大的测土配方数据库。其中,土壤全氮是测土配方土壤测试基础五项之一,是表征土壤养分的重要指标。某些地区近50%耕地面积的土壤缺氮,农业种植过程中氮肥施用量较大,因此,准确的土壤氮素数据能够为实际种植提供重要的施肥依据。测土配方施肥数据具有数据体量大、时间跨度长、人工填报汇总录入的特点,然而,因其地区田块、采集方式、测试方法、衡量标准等方面的差异,造成了数据来源复杂标准不统一、人工填报容易出错、数据缺乏校验等问题。因此,需要利用面向大数据的高效算法尤其是AI算法识别异常数据,实现阈值范围内的大数据全局、局部异常值高效捕获,提升数据可用性。
目前现有技术主要是针对局限特定领域数据的异常值检测,如面向交通、计算机、医疗等领域,这些类型的数据采集方式规范、阈值范围明确、时间关系清晰,具有明确的评价标准,容易通过算法检测。然而,农业中土壤要素受到作物、气象、管理等多方面复杂影响,时间序列检测等技术方法在呈现“离散采样、空间广泛、品类繁多”特点的测土配方施肥数据中无法开展应用,相关技术难以有效应用在农业大数据尤其是测土配方施肥数据的异常值检测。因此,需要针对农业测土配方施肥数据的群体特征、空间分布特点,提出针对性的数据异常值检测评价专有技术方法,通过深入挖掘数据多维度特征实现高效、准确的测土配方施肥土壤全氮数据异常值检测。
发明内容
针对现有技术存在的问题,本发明提供一种测土配方施肥土壤全氮异常值检测方法及装置。
本发明提供一种测土配方施肥土壤全氮异常值检测方法,包括:获取待检测土壤全氮数据,所述待检测土壤全氮数据中包括不同时期的土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;根据所述土壤全氮数据集,基于独立森林iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述根据所述待检测土壤全氮数据,确定土壤全氮数据集,包括:将历史气象数据、环境数据、作物特征数据、管理特征数据的时序数据,输入APSIM模型,得到与所述时序数据步长对应的不同时期的辅助土壤全氮数据;根据所述辅助土壤全氮数据和所述待检测土壤全氮数据,构建所述土壤全氮数据集;其中,所述辅助土壤全氮数据量大于待检测土壤全氮数据。
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,包括:根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree;遍历iForest中的每棵iTree,计算每个全氮数据在每棵树中的路径长度,根据路径长度计算每个全氮数据的异常分数;根据所述异常分数,确定每个土壤全氮数据为合理和异常的检测结果。
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree,包括:根据所述土壤全氮数据集随机采样子集,根据子集随机选取数据作为分裂值;对于子集剩下的数据进行随机选取,若选取的数据小于分裂值,则将数据放在左子树,反之放在右子树,重复递归地构造左子树和右子树,直至满足子集中只剩下一条数据或多条相同的数据或者树达到最大高度。
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述根据路径长度计算每个全氮数据的异常分数,包括,根据如下公式计算全氮数据的异常分数;
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,包括:采用皮尔森相关性算法分析多个辅助指标与土壤全氮的相关关系,选择具有强相关性的辅助指标进行异常判别;所述辅助指标包括土壤有效磷、有机质、速效钾和pH;对于全氮数据集和辅助指标数据集,计算K邻域的局部莫兰指数,得到每个数据点的相关性数值;若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况相同,则全氮数据点为合理值;若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况不同,则全氮数据点为异常;若全氮数据和辅助指标数据的相关性数值都存在异常但是异常点位或类型不一致,则全氮数据点为异常。
根据本发明提供的一种测土配方施肥土壤全氮异常值检测方法,所述计算K邻域的局部莫兰指数之前,还包括:对于每个全氮数据点周围邻近点数量,根据K值从小到大生成不同大小的邻域值;对每个邻域进行全氮数据探测,过程中根据局部莫兰指数产生异常值,确定异常值数量变化稳定的空间邻域,作为所述K邻域。
本发明还提供一种测土配方施肥土壤全氮异常值检测装置,包括:数据采集模块,用于获取待检测土壤全氮数据,所述待检测土壤全氮数据中包括不同时期的土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;第一检测模块,用于根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;第二检测模块,用于根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;综合处理模块,用于根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
本发明还提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述测土配方施肥土壤全氮异常值检测方法。
本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上述任一种所述测土配方施肥土壤全氮异常值检测方法。
本发明提供的测土配方施肥土壤全氮异常值检测方法及装置,通过基于机器学习的iForest检测模型与GIS空间统计分析相融合的模式,能够充分地挖掘数据在不同维度、不同层次上的异常性,有效地避免单一方法下所造成的漏判、误判问题,极大地提升测土配方施肥土壤全氮数据的异常值检测准确度。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的测土配方施肥土壤全氮异常值检测方法的流程示意图;
图2是本发明提供的测土配方施肥土壤全氮异常值检测装置的结构示意图;
图3是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前针对专有场景的异常值检测是数据挖掘领域的研究热点,根据霍金斯对离群值(Hawkins-Outlier)的定义,离群值是一种与其他观测值相差甚远的观测值,以至于引起人们怀疑它是由另一种机制产生的。因此,我们可以提出定义“异常”的两个假设:1)异常数据跟样本中大多数的数据不太一样;2)异常数据在整体数据样本中占比较小。此外,地理学第一定律证明:任何事物都是与其他事物相关的,越相近的事物关联越紧密。几乎所有空间数据都具有空间依赖性或空间自相关特征,空间单元与邻近单元之间针对同一属性值存在潜在的相互依赖的特性。基于以上假设,针对测土配方土壤全氮大数据的群体和空间特点,可以从以下三个方面深入挖掘:1)通过扩展数据集的体量,使得全局异常值更加“孤立”,更易于被算法识别;2)结合高效的机器学习检测算法,挖掘测土配方全氮数据特征,快速识别全局异常值;3)利用空间异常值识别方法,挖掘空间层面的突出异常,识别空间局部异常值。因此,通过数据体量扩容、数据群体全局识别、空间局部检测多种方法项组合,能够提升测土配方土壤全氮数据异常值识别能力。
本发明以测土配方施肥土壤全氮数据异常值精准检测为目标,探索结合农学机理模型、机器学习异常检测算法与GIS空间数据挖掘分析算法,形成全局与局部特征相结合的异常值检测新方法。提出一种基于模型扩展的土壤全氮异常值机器学习检测方法,识别具有统计显著性的全局、局部异常值,针对数据特征,提出异常值综合判定与效果验证方法,实现测土配方施肥土壤全氮数据的高效异常值识别检测。
下面结合图1-图3描述本发明的测土配方施肥土壤全氮异常值检测方法及装置。图1是本发明提供的测土配方施肥土壤全氮异常值检测方法的流程示意图,如图1所示,本发明提供测土配方施肥土壤全氮异常值检测方法,包括:
101、获取待检测土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;所述数据集中包括不同时期的土壤全氮数据。
根据待检测土壤全氮数据,构建待检测土壤全氮数据集中进行检测,土壤全氮数据集为不同时期的土壤全氮数据序列,从中需确定异常的全氮数据。
102、根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集。
本发明采用无监督机器学习方法筛选测土配方施肥土壤全氮全局异常数据。以步骤101获取的数据集为基础,通过带入机器学习异常检测算法,判别异常值。
本发明提出基于无监督机器学习异常值检测方法采用iForest算法设计,该方法利用异常值数据的两个特征:1)异常数据占数据集总体规模的比重较小;2)异常数据相比正常数据的属性值存在明显的差异。
本发明所提出的基于无监督学习的测土配方全氮全局数据异常值检测方法在于,在仅包含数值类型的训练集中,对数据进行递归的划分,直至隔离树iTree将每个数据与其他数据区别开来。全局异常值对隔离具有较强的敏感性,因此异常数据更接近于树的根节点,而正常数据离根节点较远,通过这种思路实现全局异常值监测。
103、根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集。
另一方面,本发明还采用基于GIS的空间统计分析异常识别方法,其思路在于:通
过利用GIS技术中的局部莫兰指数()统计量方法,识别出区域
内具有统计显著性的热点、冷点,进而判定出局部空间异常值。可选地,可结合如测土配方
基础五项中的另外四项指标作为辅助指标,开展土壤全氮数据局部时空自相关分析。
104、根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
本发明通过基于机器学习的无监督数据全局异常检测和基于GIS空间分析的局部异常值检测相结合的方法,分别挖掘大数据下的内部特征实现准确异常值判定。
进一步,整合机器学习和空间分析所得出的异常值判定结果,筛选所有数据中标
签为“”的数据,结合机器学习异常识别方法所得出的结果/(合理/异常)与空
间分析得出的结果/(合理/异常)。判定规则把包括:若结果组合为+,则数据
判定为正常值;若结果组合为+,则判定数据为异常值;若结果组合为+或+,则判定数据为疑似异常值。其中,正常值表示具有较准确性高置信度的数据,被认
为是准确可靠的数据;异常值表示明确具有异常特征的数据,被认为是需要被剔除和修正
的数据;疑似异常值数据表示在某个层面出现了异常特征,但又符合某种现实逻辑,根据实
际数据应用场景及数据准确性要求程度进行判别归类。
针对本发明所提出的测土配方土壤全氮数据异常值检测方法的效果识别,提出
“盲样法”进行验证,即在数据集中手动添加N个异常值(),进而根据识别结果验证
效果。正确率(accuracy)的计算方法可以表达为:
式中,TP表示被准确识别出的盲样异常值个数;TN表示未被识别出的盲样异常值个数。正确率越高,异常值检测方法效果越好。
本发明提供的测土配方施肥土壤全氮异常值检测方法,通过基于机器学习的iForest检测模型与GIS空间统计分析相融合的模式,能够充分地挖掘数据在不同维度、不同层次上的异常性,有效地避免单一方法下所造成的漏判、误判问题,极大地提升测土配方施肥土壤全氮数据的异常值检测准确度。
在一个实施例中,所述根据所述待检测土壤全氮数据,确定土壤全氮数据集,包括:将历史气象数据、环境数据、作物特征数据、管理特征数据的时序数据,输入APSIM模型,得到与所述时序数据步长对应的不同时期的辅助土壤全氮数据,根据所述辅助土壤全氮数据和所述待检测土壤全氮数据,构建所述土壤全氮数据集。
本发明提出采用农学机理模型构建扩展数据集为后续智能模型训练提供充足“素材”的思路,采用APSIM模型作为底层核心农学机理模型基础,其中,气候数据是APSIM模型的基础驱动要素,模型在逐日气象要素的驱动下,对作物生理生态、土壤剖面水肥动态和土壤侵蚀量等过程进行数值模拟,实现土壤氮素演化变化模拟。
本发明采用的农学机理模型APSIM所需的气象驱动要素包括:逐日太阳辐射量(MJ/m2)、逐日最高气温(℃)、逐日最低气温(℃)、逐日降水量(mm)、纬度(十进制度)、月平均气温(℃)和月均温变化(℃)等参数项。此外,需要初始化种植环境信息,包含初始土壤氮含量、初始土壤磷含量、初始土壤含水量、初始土壤pH等,为了进一步精准模拟土壤氮素的演变,需要记录种植管理过程信息包括:施肥情况、灌溉情况等。
针对土壤全氮的演化模拟,APSIM模型中氮在作物中的吸收规律取决于土壤的供氮能力、作物需氮量和作物的实际含氮量的变化。其中,针对土壤全氮的组分,土壤供氮的变化是模拟无机氮通过质流和扩散被作物吸收的总和,吸收氮素的公式表达为:
其中,作物需氮模拟作物吸收氮素与氮素迁移规律的过程表现,其模拟过程机理是作物生长发育的任何一天作物的需氮量是当日作需氮量与前日作物需氮亏缺量的加和,而作物当日需氮量由当日作物潜在生物量的增量与作物各部分器官的含氮百分数决定,潜在生物量的增量由作物的光合作物固定的有机物来决定。因此,作物需氮亏缺量是作物含氮量与含氮临界值之间的差,该过程可以表达为:
本发明所提出的数据集构建方法在于,通过所收集的长时序逐日历史气象、环境、
作物、管理数据,带入APSIM模型进行模拟,记录逐日的土壤氮素变化结果。其中,初始指标
通过随机选择方式选取,土壤全氮选取范围为[0.3g/kg, 6g/kg],有效磷选取范围为[3mg/
kg, 100mg/kg],有机质选取范围为[3g/kg, 45g/kg],速效钾选取范围为[30mg/kg,
300mg/kg],pH选取范围为[4, 9],土壤含水量选取范围为[40%, 80%]。根据本发明所提出
的基于模型的扩展方法,将模型计算所生成的逐日土壤全氮值进行记录。原始数据
标记数据类型为“Ori”,所有模型生成的扩展数据标记类型为“Gen”。
最终,将基于APSIM模型的土壤全氮数据作为辅助,和待检测的土壤全氮数据一起构建为土壤全氮数据集。其中的历史气象数据、环境数据、作物特征数据、管理特征数据的时序数据,可以为待检测土壤全氮数据相同情况下的历史气象数据、环境数据、作物特征数据、管理特征数据的时序数据。同时,辅助土壤全氮数据量大于待检测土壤全氮数据,如大于可以为两倍以上,或者十倍百倍。
本发明实施例的测土配方施肥土壤全氮异常值检测方法,通过对待检测土壤全氮数据进行大量扩充,可以避免仅依据少量待检测土壤全氮数据进行检测时准确度不高的问题,同时通过APSIM模型生成海量的可靠全氮数据,对待检测土壤全氮数据进行海里数据扩充,进一步提高检测精度。
在一个实施例中,所述根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,包括:根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree;遍历iForest中的每棵iTree,计算每个全氮数据在每棵树中的路径长度,根据路径长度计算每个全氮数据的异常分数;根据所述异常分数,确定每个土壤全氮数据为合理和异常的检测结果。
本发明所采用的方法通过构造一定数目的iTree来组成iForest。随机采样提取土
壤全氮数据集的子集来构造每棵iTree,以保证iTree的多样性。通过遍历iForest中的每
棵iTree,计算土壤全氮数据集中数据d在每棵树中的路径长度,然后根据其路径长度计
算d的异常分数,从而判断d是否异常。
在一个实施例中,所述根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree,包括:根据所述土壤全氮数据集随机采样子集,根据子集随机选取数据作为分裂值;对于子集剩下的数据进行随机选取,若选取的数据小于分裂值,则将数据放在左子树,反之放在右子树,重复递归地构造左子树和右子树,直至满足子集中只剩下一条数据或多条相同的数据或者树达到最大高度。
iForest算法的核心在于构建由iTree组成的森林(iForest)。其中,iTree定义为
令T是一棵二叉树,N是T的节点,若N是叶子节点,则称其为外部节点,若N是一个具有两个子
节点,则称其为内部节点。iTree从数据集中随机地选择属性A和
分裂值P,然后按照A的值(记为)对每个数据进行划分。如果,则将数据放在左子树,反之放在右子树。以此递归地构造左子树和右子树,直至满足下列条件之
一:1)D中只剩下一条数据或多条相同的数据;2)树达到最大高度。此外,路径长度定义为在
一棵iTree中,从根节点到外部节点所经历边的数目称为路径长度,记为。由于iTree与
二叉查找树的结构等价,因此包含d的叶子节点的路径长度等于二叉查找树中失败查询的
路径长度。给定数据集D,二叉查找树中失败查询的路径长度:
本发明所采用的方法通过构造一定数目的iTree来组成iForest。随机采样提取
的子集来构造每棵iTree,以保证的多样性。通过遍历iForest中的每棵iTree,计算
数据在每棵树中的路径长度,然后根据其路径长度计算的异常分数,从而判断是否
异常。
在一个实施例中,所述根据路径长度计算每个全氮数据的异常分数,包括,根据如下公式计算全氮数据的异常分数;
当时,S=0.5,即当所有数据均返回的S≈0.5时,全部样本中没
有明显的异常值;当时,S=1,即当数据返回的S非常接近于1时,它们是异常
值;当时,S=0,即当数据返回的S远小于0.5时,它们有很大的可能为正
常值。
本步骤所提出的基于无监督学习的机器学习全局异常值检测方法,通过上述步骤
完成对测土配方土壤全氮扩展数据集的异常值筛选,将每条数据分类为“合理”和“异常”两
类,筛选结果被记录入数据库,用于后续步骤的综合判别,其中,基于无监督机器学习异常
值检测判定为合理的点标记为,判定为异常的点标记为。
本发明的测土配方施肥土壤全氮异常值检测方法,通过采用高性能的无监督学习算法的异常值检测方法,实现面向大体量数据的全局异常值识别,立足大数据环境,进一步改进算法提升异常值识别效率,结合随机投影策略加速运行框架改进iForest算法,有效提升大体量、具体化场景下全局异常值检测的性能。
在一个实施例中,所述根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,包括:采用皮尔森相关性算法分析多个辅助指标与土壤全氮的相关关系,选择具有强相关性的辅助指标进行异常判别;所述辅助指标包括土壤有效磷、有机质、速效钾和pH;对于全氮数据集和辅助指标数据集,计算K邻域的局部莫兰指数,得到每个数据点的相关性数值;若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况相同,则全氮数据点为合理值;若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况不同,则全氮数据点为异常;若全氮数据和辅助指标数据的相关性数值都存在异常但是异常点位或类型不一致,则全氮数据点为异常。
本发明以空间自相关理论中的局部莫兰指数方法进行异常值判定,分析观测值的高值或低值聚集、高低或低高异常值分布。最终,基于相关性分析和辅助指标识别测土配方施肥土壤全氮数据异常点。
为提高莫兰指数判断异常数据后,检测数据异常或正常的准确度。本发明选择测
土配方土壤有效磷、有机质、速效钾和pH作为土壤全氮指标的辅助指标,采用皮尔森相关性
算法分析多个辅助指标与土壤全氮的相关关系,选择具有强相关性的辅助指标进行异常判
别:如果检测指标数据的异常值检测结果与辅助指标数据的异常值检测结果相同,即莫兰
指数计算得到的相关性数值均表现为高低“H-L”型或者低高“L-H”型,则认为该点属于合理
范畴非异常值;如果检测指标数据的异常值检测结果与辅助指标数据的异常值检测结果不
同,即检测指标存在异常而辅助指标不存在异常,或者两者都存在异常但是异常点位或类
型不一致,则认为该点明确存在异常。其中,基于GIS空间统计分析异常识别判定为合理的
点标记为,判定为异常的点标记为。
在一个实施例中,所述计算K邻域的局部莫兰指数之前,还包括:对于每个全氮数据点周围邻近点数量,根据K值从小到大生成不同大小的邻域值;对每个邻域进行全氮数据探测,过程中根据局部莫兰指数产生异常值,确定异常值数量变化稳定的空间邻域,作为所述K邻域。
首先要设定每个点周围制定的合理的K个临近点,该参数为后续计算莫兰指数的关键参数,过大和过小都会导致无法准确计算异常数据。考虑到本发明涉及的空间异常值分析面向的海量数据且分布区域广,作为可选实施例,设定K的取值范围为[10,M/5]的整数,M为全氮数据集的样本总数。
为了得到合适的K邻域,进一步提高莫兰指数查找异常数据的准确度。本发明实施例通过迭代方法寻找准确的K值。
K值由小到大逐一迭代,每次迭代分布带入局部莫兰指数对全部数据运算产生异
常值结果,通过持续迭代,能够产生连续、稳定数量的异常值判断结果的K值。最终,在生成
的空间权重矩阵中,临近的个点在矩阵中的结果为。局部莫兰指数依次利用各个
值生成的权重矩阵对空间点分布重新定义,生成不同权重矩阵对应下的异常值分布情况。
其中,局部莫兰指数计算表达式如下:
本发明所提出的方法立足大数据多维度数据挖掘分析,其本质决定了数据体量巨大、运算量庞大,为了提升针对测土配方土壤全氮异常值检测的效率,本发明能够有效提升模型识别的鲁棒性,避免过模型拟合所导致的误分类,配合iForest等基于树的模型所天然具备的并行计算能力,极大地提升土壤全氮的异常值识别效率。
下面对本发明提供的测土配方施肥土壤全氮异常值检测装置进行描述,下文描述的测土配方施肥土壤全氮异常值检测装置与上文描述的测土配方施肥土壤全氮异常值检测方法可相互对应参照。
图2是本发明提供的测土配方施肥土壤全氮异常值检测装置的结构示意图,如图2所示,该测土配方施肥土壤全氮异常值检测装置包括:数据采集模块201、第一检测模块202、第二检测模块203和综合处理模块204。其中,数据采集模块201用于获取待检测土壤全氮数据,所述待检测土壤全氮数据中包括不同时期的土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;第一检测模块202用于根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;第二检测模块203用于根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;综合处理模块204用于根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
本发明实施例提供的装置实施例是为了实现上述各方法实施例的,具体流程和详细内容请参照上述方法实施例,此处不再赘述。
本发明实施例所提供的测土配方施肥土壤全氮异常值检测装置,其实现原理及产生的技术效果和前述测土配方施肥土壤全氮异常值检测方法实施例相同,为简要描述,测土配方施肥土壤全氮异常值检测装置实施例部分未提及之处,可参考前述测土配方施肥土壤全氮异常值检测方法实施例中相应内容。
图3是本发明提供的电子设备的结构示意图,如图3所示,该电子设备可以包括:处理器(processor)301、通信接口(Communications Interface)302、存储器(memory)303和通信总线304,其中,处理器301,通信接口302,存储器303通过通信总线304完成相互间的通信。处理器301可以调用存储器303中的逻辑指令,以执行测土配方施肥土壤全氮异常值检测方法,该方法包括:获取待检测土壤全氮数据,所述待检测土壤全氮数据中包括不同时期的土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
此外,上述的存储器303中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的测土配方施肥土壤全氮异常值检测方法,该方法包括:获取待检测土壤全氮数据,所述待检测土壤全氮数据中包括不同时期的土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种测土配方施肥土壤全氮异常值检测方法,其特征在于,包括:
获取待检测土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;所述土壤全氮数据集中包括不同时期的土壤全氮数据;
根据所述土壤全氮数据集,基于独立森林iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;
根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;
根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
2.根据权利要求1所述的测土配方施肥土壤全氮异常值检测方法,其特征在于,所述根据所述待检测土壤全氮数据,确定土壤全氮数据集,包括:
将历史气象数据、环境数据、作物特征数据、管理特征数据的时序数据,输入农业生产系统模拟器APSIM模型,得到与所述时序数据步长对应的不同时期的辅助土壤全氮数据;
根据所述辅助土壤全氮数据和所述待检测土壤全氮数据,构建所述土壤全氮数据集;
其中,所述辅助土壤全氮数据量大于待检测土壤全氮数据。
3.根据权利要求1所述的测土配方施肥土壤全氮异常值检测方法,其特征在于,所述根据所述土壤全氮数据集,基于iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,包括:
根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree;
遍历iForest中的每棵iTree,计算每个全氮数据在每棵树中的路径长度,根据路径长度计算每个全氮数据的异常分数;
根据所述异常分数,确定每个土壤全氮数据为合理和异常的检测结果。
4.根据权利要求3所述的测土配方施肥土壤全氮异常值检测方法,其特征在于,所述根据所述土壤全氮数据集中随机采样的每个子集,构造iForest检测模型的每棵iTree,包括:
根据所述土壤全氮数据集随机采样子集,根据子集随机选取数据作为分裂值;
对于子集剩下的数据进行随机选取,若选取的数据小于分裂值,则将数据放在左子树,反之放在右子树,重复递归地构造左子树和右子树,直至满足子集中只剩下一条数据或多条相同的数据或者树达到最大高度。
6.根据权利要求1所述的测土配方施肥土壤全氮异常值检测方法,其特征在于,所述根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,包括:
采用皮尔森相关性算法分析多个辅助指标与土壤全氮的相关关系,选择具有强相关性的辅助指标进行异常判别;所述辅助指标包括土壤有效磷、有机质、速效钾和pH;
对于全氮数据集和辅助指标数据集,计算K邻域的局部莫兰指数,得到每个数据点的相关性数值;
若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况相同,则全氮数据点为合理值;若全氮数据的异常值检测结果与辅助指标数据的异常值检测结果的相关性数值高低或底高变化情况不同,则全氮数据点为异常;若全氮数据和辅助指标数据的相关性数值都存在异常但是异常点位或类型不一致,则全氮数据点为异常。
7.根据权利要求6所述的测土配方施肥土壤全氮异常值检测方法,其特征在于,所述计算K邻域的局部莫兰指数之前,还包括:
对于每个全氮数据点周围邻近点数量,根据K值从小到大生成不同大小的邻域值;
对每个邻域进行全氮数据探测,过程中根据局部莫兰指数产生异常值,确定异常值数量变化稳定的空间邻域,作为所述K邻域。
8.一种测土配方施肥土壤全氮异常值检测装置,其特征在于,包括:
数据采集模块,用于获取待检测土壤全氮数据,根据所述待检测土壤全氮数据确定土壤全氮数据集;所述土壤全氮数据集中包括不同时期的土壤全氮数据;
第一检测模块,用于根据所述土壤全氮数据集,基于独立森林iForest检测模型,得到每个土壤全氮数据为合理和异常的检测结果,作为第一全氮数据集;
第二检测模块,用于根据所述土壤全氮数据集,基于GIS空间统计分析的异常检测方法,确定每个土壤全氮数据为合理和异常的检测结果,作为第二全氮数据集;
综合处理模块,用于根据所述第一全氮数据集和所述第二全氮数据集中的合理和异常结果,确定土壤全氮数据集中所有数据的合理或异常检测结果。
9.一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7任一项所述测土配方施肥土壤全氮异常值检测方法。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述测土配方施肥土壤全氮异常值检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210650110.7A CN114720665B (zh) | 2022-06-10 | 2022-06-10 | 测土配方施肥土壤全氮异常值检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210650110.7A CN114720665B (zh) | 2022-06-10 | 2022-06-10 | 测土配方施肥土壤全氮异常值检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114720665A true CN114720665A (zh) | 2022-07-08 |
CN114720665B CN114720665B (zh) | 2022-09-16 |
Family
ID=82233076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210650110.7A Active CN114720665B (zh) | 2022-06-10 | 2022-06-10 | 测土配方施肥土壤全氮异常值检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114720665B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117571623A (zh) * | 2024-01-17 | 2024-02-20 | 北京市农林科学院智能装备技术研究中心 | 一种便携式测量土壤全磷的检测装置及方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884661A (zh) * | 2014-02-21 | 2014-06-25 | 浙江大学 | 基于土壤可见-近红外光谱库的土壤全氮实时检测方法 |
CN105550938A (zh) * | 2015-12-03 | 2016-05-04 | 中国农业大学 | 县域耕地质量评价成果异常值检验方法 |
CN107607486A (zh) * | 2017-09-25 | 2018-01-19 | 中国农业大学 | 一种土壤全氮检测方法及装置 |
CN107764976A (zh) * | 2017-09-28 | 2018-03-06 | 广州市健坤网络科技发展有限公司 | 土壤氮素快速诊断方法及在线监测系统 |
CN109345137A (zh) * | 2018-10-22 | 2019-02-15 | 广东精点数据科技股份有限公司 | 一种基于农业大数据的异常值检测方法 |
CN112365391A (zh) * | 2020-11-12 | 2021-02-12 | 河北工程大学 | 一种基于“国土调查”数据的用地多样性计量方法 |
CN112861770A (zh) * | 2021-03-02 | 2021-05-28 | 清华大学 | 一种水质污染预警方法及装置 |
CN113486134A (zh) * | 2021-07-20 | 2021-10-08 | 杭州鲁尔物联科技有限公司 | 降雨量异常检测方法、装置、计算机设备及存储介质 |
CN113902249A (zh) * | 2021-09-02 | 2022-01-07 | 北京市农林科学院信息技术研究中心 | 一种土壤重金属影响因素解析方法及装置 |
US20220044133A1 (en) * | 2020-08-07 | 2022-02-10 | Sap Se | Detection of anomalous data using machine learning |
CN114385639A (zh) * | 2022-01-13 | 2022-04-22 | 湖北中南鹏力海洋探测系统工程有限公司 | 一种近海地波雷达数据存储方法及合成方法 |
CN114486786A (zh) * | 2022-03-03 | 2022-05-13 | 上海园林绿化建设有限公司 | 土壤有机质测定方法及测定系统 |
-
2022
- 2022-06-10 CN CN202210650110.7A patent/CN114720665B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884661A (zh) * | 2014-02-21 | 2014-06-25 | 浙江大学 | 基于土壤可见-近红外光谱库的土壤全氮实时检测方法 |
CN105550938A (zh) * | 2015-12-03 | 2016-05-04 | 中国农业大学 | 县域耕地质量评价成果异常值检验方法 |
CN107607486A (zh) * | 2017-09-25 | 2018-01-19 | 中国农业大学 | 一种土壤全氮检测方法及装置 |
CN107764976A (zh) * | 2017-09-28 | 2018-03-06 | 广州市健坤网络科技发展有限公司 | 土壤氮素快速诊断方法及在线监测系统 |
CN109345137A (zh) * | 2018-10-22 | 2019-02-15 | 广东精点数据科技股份有限公司 | 一种基于农业大数据的异常值检测方法 |
US20220044133A1 (en) * | 2020-08-07 | 2022-02-10 | Sap Se | Detection of anomalous data using machine learning |
CN112365391A (zh) * | 2020-11-12 | 2021-02-12 | 河北工程大学 | 一种基于“国土调查”数据的用地多样性计量方法 |
CN112861770A (zh) * | 2021-03-02 | 2021-05-28 | 清华大学 | 一种水质污染预警方法及装置 |
CN113486134A (zh) * | 2021-07-20 | 2021-10-08 | 杭州鲁尔物联科技有限公司 | 降雨量异常检测方法、装置、计算机设备及存储介质 |
CN113902249A (zh) * | 2021-09-02 | 2022-01-07 | 北京市农林科学院信息技术研究中心 | 一种土壤重金属影响因素解析方法及装置 |
CN114385639A (zh) * | 2022-01-13 | 2022-04-22 | 湖北中南鹏力海洋探测系统工程有限公司 | 一种近海地波雷达数据存储方法及合成方法 |
CN114486786A (zh) * | 2022-03-03 | 2022-05-13 | 上海园林绿化建设有限公司 | 土壤有机质测定方法及测定系统 |
Non-Patent Citations (2)
Title |
---|
何为等: "基于LSTM的催化裂化装置NO_x排放预测模型及应用", 《西安石油大学学报(自然科学版)》 * |
厉彦玲等: "黄河三角洲典型地区耕地土壤养分空间预测", 《自然资源学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117571623A (zh) * | 2024-01-17 | 2024-02-20 | 北京市农林科学院智能装备技术研究中心 | 一种便携式测量土壤全磷的检测装置及方法 |
CN117571623B (zh) * | 2024-01-17 | 2024-04-16 | 北京市农林科学院智能装备技术研究中心 | 一种便携式测量土壤全磷的检测装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114720665B (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mao et al. | Remote sensing evaluation of CLM4 GPP for the period 2000–09 | |
Salvati et al. | The environmental “risky” region: identifying land degradation processes through integration of socio-economic and ecological indicators in a multivariate regionalization model | |
Whittaker et al. | Detection of overparameterization and overfitting in an automatic calibration of SWAT | |
CN113222242A (zh) | 一种森林火险等级预报方法、装置、计算机设备及存储介质 | |
CN114720665B (zh) | 测土配方施肥土壤全氮异常值检测方法及装置 | |
Sahu et al. | Testing some grouping methods to achieve a low error quantile estimate for high resolution (0.25 x 0.25) precipitation data | |
CN116129262A (zh) | 一种面向宜机化改造的耕地适宜性评价方法及系统 | |
Zhang et al. | Analysis of flash droughts in China using machine learning | |
CN113139717B (zh) | 作物苗情分级遥感监测方法及装置 | |
He et al. | Physics guided neural networks for time-aware fairness: an application in crop yield prediction | |
Toma et al. | Historical and future dynamics of land use land cover and its drivers in Ajora‐Woybo watershed, Omo‐Gibe basin, Ethiopia | |
Peña Quiñones et al. | Radius of influence of air temperature from automated weather stations installed in complex terrain | |
He et al. | Does partition matter? A new approach to modeling land use change | |
Kuriakose et al. | Indian Crop Yield Prediction using LSTM Deep Learning Networks | |
Vanarase et al. | Crop Prediction Using Data Mining and Machine Learning Techniques | |
Golian et al. | Evaluation of sub-selection methods for assessing climate change impacts on low-flow and hydrological drought conditions | |
Zhang et al. | Analysis of Flash Drought in China using Artificial Intelligence models | |
Chen et al. | Study on the impact of low‐temperature stress on winter wheat based on multi‐model coupling | |
Taylor et al. | Improving spatial predictions of Eucalypt plantation growth by combining interpretable machine learning with the 3-PG model | |
Komolafe et al. | Predictive Modeling for Land Suitability Assessment for Cassava Cultivation | |
CN117828312B (zh) | 一种流域水文环境的管理方法及相关设备 | |
Gao et al. | Ecological risk assessment of grassland vegetation change based on Bayesian model in Xilin Gol League, China | |
Huang et al. | Air quality monitoring of landscape architecture based on multi-sensor fusion | |
CN113344247B (zh) | 一种基于深度学习的电力设施选址预测方法与系统 | |
Subhashi et al. | Selecting Suitable Plants for a Given Area using Data Analysis Approaches |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |