CN113111799A - 一种基于集合经验模态分解的耕地土壤肥力水平监测方法 - Google Patents
一种基于集合经验模态分解的耕地土壤肥力水平监测方法 Download PDFInfo
- Publication number
- CN113111799A CN113111799A CN202110418313.9A CN202110418313A CN113111799A CN 113111799 A CN113111799 A CN 113111799A CN 202110418313 A CN202110418313 A CN 202110418313A CN 113111799 A CN113111799 A CN 113111799A
- Authority
- CN
- China
- Prior art keywords
- soil fertility
- soil
- stress
- nutrient
- empirical mode
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 134
- 230000035558 fertility Effects 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 65
- 238000012544 monitoring process Methods 0.000 title claims abstract description 27
- 230000000694 effects Effects 0.000 claims abstract description 47
- 230000012010 growth Effects 0.000 claims abstract description 42
- 238000012216 screening Methods 0.000 claims abstract description 31
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 230000003595 spectral effect Effects 0.000 claims abstract description 16
- 239000002131 composite material Substances 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 230000035882 stress Effects 0.000 claims description 94
- 235000015097 nutrients Nutrition 0.000 claims description 61
- 241000209140 Triticum Species 0.000 claims description 39
- 235000021307 Triticum Nutrition 0.000 claims description 39
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 22
- 230000014509 gene expression Effects 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 20
- 241000607479 Yersinia pestis Species 0.000 claims description 14
- 229910001385 heavy metal Inorganic materials 0.000 claims description 13
- 241000196324 Embryophyta Species 0.000 claims description 11
- 241000238631 Hexapoda Species 0.000 claims description 11
- 201000010099 disease Diseases 0.000 claims description 11
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 11
- 229910052757 nitrogen Inorganic materials 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 8
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 claims description 8
- 229910052698 phosphorus Inorganic materials 0.000 claims description 8
- 239000011574 phosphorus Substances 0.000 claims description 8
- 229910052700 potassium Inorganic materials 0.000 claims description 8
- 239000011591 potassium Substances 0.000 claims description 8
- 238000000513 principal component analysis Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 7
- 230000006353 environmental stress Effects 0.000 claims description 4
- 239000005416 organic matter Substances 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 abstract description 2
- 230000009471 action Effects 0.000 description 11
- 238000001556 precipitation Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 8
- 230000001932 seasonal effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 208000005156 Dehydration Diseases 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000007774 longterm Effects 0.000 description 4
- 238000010219 correlation analysis Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 239000004016 soil organic matter Substances 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000003973 irrigation Methods 0.000 description 2
- 230000002262 irrigation Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000006355 external stress Effects 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 235000021049 nutrient content Nutrition 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008635 plant growth Effects 0.000 description 1
- 230000010152 pollination Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005527 soil sampling Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/58—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
- G06F16/5866—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using information manually generated, e.g. tags, keywords, comments, manually generated location and time information
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Library & Information Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Multimedia (AREA)
- Game Theory and Decision Science (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Cultivation Of Plants (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种基于集合经验模态分解的耕地土壤肥力水平监测方法,包括:通过高分一号遥感探测器采集耕地作物的光谱影像作为数据源;对数据源的光谱影像进行图像预处理,提取NDVI长时间序列;基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中;结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取。本发明以EEMD方法为主体进行土壤肥力胁迫作用甄别提取,能够有效获取单一胁迫因素对于耕地作物长势的影响,提高了土壤肥力水平监测效率。
Description
技术领域
本发明涉及土壤肥力水平监测技术领域,特别是涉及一种基于集合经验模态分解的耕地土壤肥力水平监测方法。
背景技术
传统非遥感领域的土壤肥力水平监测方法普遍采用试验田的手段,通过定期进行土壤采样化验及农作物产量测算等途径对土壤肥力水平变化进行监测。这种方法耗费的人力物力资源多,且监测结果受空间和时间的约束较大,难以实现大范围动态监测。
相比之下,遥感技术在便捷性及时效性等方面都有着极大的提高。王之杰等利用不同氮素水平下冬小麦冠层叶片的光谱反射特点差异,对冬小麦冠层氮素分布对其籽粒蛋白质含量的影响进行了研究;田玉珂等利用TM遥感数据对土壤肥力各项指标的预测模型进行蒙特卡罗模拟修正,获得了有机质、全氮、碱解氮、有效磷以及速效钾的最佳修正光谱参数;王磊等利用无人机成像高光谱数据对土壤有机质(SOM)的敏感光谱特征进行了分析,提出了一种土壤SOM快速、准确、无损估测的监测方法,为遥感领域精准农业的发展以及成像高光谱技术的应用提供了一定的参考作用。
目前,这些基于遥感技术的土壤肥力水平监测方法虽然能够满足很多情况下的需求,但其仍存在着几个层面上的缺陷与不足:一是土壤中部分速效养分的含量差异在植被光谱特征上的显著性不足,导致其最终的反演精度受限;二是土壤表层与深层养分差异较大,且深层土壤的光谱数据在遥感影像中获取困难,导致数据准确性不足等问题;另外,一般农作物长势的遥感反演结果均包含着多种胁迫因素的共同影响,如何对土壤肥力胁迫进行有效而准确的甄别也是遥感技术在土壤肥力领域需要解决的重要难题。而农作物长势受到多种胁迫因素的影响较为复杂,单一胁迫作用甄别结果容易受到其他胁迫因素的影响;此外,影响土壤肥力水平的养分指标较多,遥感光谱数据对单一养分指标含量差异的响应不够敏感。
发明内容
本发明的目的是提供一种基于集合经验模态分解的耕地土壤肥力水平监测方法,在EEMD分解过程中对不同于土壤肥力时间尺度的胁迫作用进行的提取,并将复杂的土壤肥力系统用一个综合评价数值进行表达,使遥感光谱数据对于其数值差异的响应程度更为明显。
为实现上述目的,本发明提供了如下方案:
一种基于集合经验模态分解的耕地土壤肥力水平监测方法,包括如下步骤:
步骤1,通过高分一号遥感探测器采集耕地作物的光谱影像作为数据源;
步骤2,对数据源的光谱影像进行图像预处理,并将归一化植被指数NDVI作为耕地作物长势的反演参数,提取NDVI长时间序列;
步骤3,基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中;
步骤4,结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取。
进一步的,所述方法还包括:结合测得的耕地土壤养分指标数据,从统计学角度对甄别提取结果进行验证,验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性。
进一步的,所述步骤3,基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中,具体包括:
采用集合经验模态分解EEMD将NDVI长时间序列按照年内、年间、年际三个时间尺度进行分解,得到多个IMF分量;
将水分、病虫害、土壤肥力以及重金属污染多种胁迫因素分解到不同的IMF分量中。
进一步的,采用集合经验模态分解EEMD将NDVI长时间序列按照年内、年间、年际三个时间尺度进行分解,得到多个IMF分量,其中,所述集合经验模态分解EEMD进行分解的过程包括:
对一组NDVI长时间序列的原始序列,加入高斯白噪声;
采用三次样条插值法得到加入高斯白噪声后的原始序列的上下包络线及均值;
减去原始序列的均值,得到新序列;
基于新序列得到一个IMF分量;
重复上述过程,进行多次分解,得到多个IMF分量。
进一步的,所述步骤4中,结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取,具体包括:
提取出代表水分、病虫害这类瞬时胁迫作用的年内波动组分序列;
提取代表耕地作物固有生长趋势的年间波动组分序列;
提取代表土壤肥力以及重金属污染这类环境胁迫作用的年际波动组分序列。
进一步的,所述结合测得的耕地土壤养分指标数据,从统计学角度对甄别提取结果进行验证,验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性,具体包括:
在SPSS软件环境中对有机质、全氮、有效磷、速效钾四项养分指标数据测定结果进行相关性检验,得到相关性矩;
根据相关型矩阵对四项养分指标数据进行主成分分析,得到土壤养分总体水平关于四项养分指标的线性表达式;
将甄别提取结果与所述土壤养分总体水平关于四项养分指标的线性表达式采用最小二乘法进行线性拟合,基于拟合结果验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性。
进一步的,根据相关型矩阵对四项养分指标数据进行主成分分析,得到土壤养分总体水平关于四项养分指标的线性表达式,具体包括:
根据相关型矩阵对四项养分指标数据进行主成分分析,选择前三项主成分作为评价土壤肥力水平的一组新变量,如下式所示:
F=0.5759F1+0.3172F2+0.0985F3+Δ(F)
式中,F表示土壤养分总体水平;Fi表示第i个主成分,Δ(F)表示余项;
基于新变量新变量关于四项养分指标的线性表达式,得到土壤养分总体水平关于四种土壤养分指标的线性表达式;
新变量关于四项养分指标的线性表达式如下所示:
F1=0.5623X1+0.3219X2+0.0264X3+0.5444X4;
F2=-0.1352X1-0.2097X2-0.8719X3+0.4214X4;
F3=0.7921X1-0.3073X2-0.2684X3-0.4540X4;
式中,X1表示有机质含量;X2表示全氮含量;X3表示有效磷含量;X4表示速效钾含量;
土壤养分总体水平关于四种土壤养分指标的线性表达式如下:
F=0.359X1+0.288X2+0.149X3+0.202X4
进一步的,所述耕地作物选择冬小麦。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的基于集合经验模态分解的耕地土壤肥力水平监测方法,通过集合经验模态分解方法—EEMD,对冬小麦NDVI长时间序列进行分解,获取年内、年间、年际三个时间尺度的胁迫作用曲线,并结合野外实测数据及统计学方法确定了代表土壤肥力胁迫作用的组分,实现土壤肥力胁迫的甄别与提取;其中,EEMD方法依照序列本身的时间尺度进行分解,无需借助任何基函数,相较于传统时域分析方法,EEMD在处理非线性非平稳序列过程中具有明显的优势,理论上能够实现任何信号的分解,同时,不同胁迫作用对于冬小麦长势影响也存在时间尺度上的差异,因此,以EEMD方法为主体的土壤肥力胁迫作用甄别提取方法能够有效获取单一胁迫因素对于冬小麦长势的影响,对土壤肥力水平的监测有很大的帮助。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例基于集合经验模态分解的耕地土壤肥力水平监测方法的流程示意图;
图2(a)为本发明实施例加入高斯白噪声的原始序列曲线图;
图2(b)为本发明实施例三次样条插值法得到的序列上下包络线及均值曲线图;
图2(c)为本发明实施例原始序列减去均值后的新序列的曲线图;
图2(d)为本发明实施例得到的第一个IMF分量的曲线图;
图2(e)为本发明实施例EEMD单次分解结果的曲线图;
图3为本发明实施例年内波动组分示意图;
图4为本发明实施例年间波动组分示意图;
图5(a)为本发明实施例季节平均模型示意图;
图5(b)为本发明实施例季节项分量示意图;
图6为本发明实施例年际波动组分示意图;
图7(a)为本发明实施例大厂县2015年降水量示意图;
图7(b)为本发明实施例大厂县2016年降水量示意图;
图7(c)为本发明实施例大厂县2017年降水量示意图;
图7(d)为本发明实施例大厂县2018年降水量示意图;
图7(e)为本发明实施例大厂县2019年降水量示意图;
图8(a)为本发明实施例陈府镇土壤肥力胁迫作用水平曲线图;
图8(b)为本发明实施例东马各庄土壤肥力胁迫作用水平曲线图;
图8(c)为本发明实施例兰庄户村土壤肥力胁迫作用水平曲线图;
图8(d)为本发明实施例梁庄村土壤肥力胁迫作用水平曲线图;
图8(e)为本发明实施例荣马坊村土壤肥力胁迫作用水平曲线图;
图8(f)为本发明实施例王必屯村土壤肥力胁迫作用水平曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于集合经验模态分解的耕地土壤肥力水平监测方法,在EEMD分解过程中对不同于土壤肥力时间尺度的胁迫作用进行的提取,并将复杂的土壤肥力系统用一个综合评价数值进行表达,使遥感光谱数据对于其数值差异的响应程度更为明显。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明实施例提供的基于集合经验模态分解的耕地土壤肥力水平监测方法,包括如下步骤:
步骤1,通过高分一号遥感探测器采集耕地作物的光谱影像作为数据源;
步骤2,对数据源的光谱影像进行图像预处理,并将归一化植被指数NDVI作为耕地作物长势的反演参数,提取NDVI长时间序列;
步骤3,基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中;具体包括:
采用集合经验模态分解EEMD将NDVI长时间序列按照年内、年间、年际三个时间尺度进行分解,得到多个IMF分量;
将水分、病虫害、土壤肥力以及重金属污染多种胁迫因素分解到不同的IMF分量中。
其中,所述集合经验模态分解EEMD进行分解的过程包括:
对一组NDVI长时间序列的原始序列,加入高斯白噪声;
采用三次样条插值法得到加入高斯白噪声后的原始序列的上下包络线及均值;
减去原始序列的均值,得到新序列;
基于新序列得到一个IMF分量;
重复上述过程,进行多次分解,得到多个IMF分量。
步骤4,结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取,具体包括:
提取出代表水分、病虫害这类瞬时胁迫作用的年内波动组分序列;
提取代表耕地作物固有生长趋势的年间波动组分序列;
提取代表土壤肥力以及重金属污染这类环境胁迫作用的年际波动组分序列。
所述方法还包括:结合测得的耕地土壤养分指标数据,从统计学角度对甄别提取结果进行验证,验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性,具体包括:
在SPSS软件环境中对有机质、全氮、有效磷、速效钾四项养分指标数据测定结果进行相关性检验,得到相关性矩;
根据相关型矩阵对四项养分指标数据进行主成分分析,得到土壤养分总体水平关于四项养分指标的线性表达式;
将甄别提取结果与所述土壤养分总体水平关于四项养分指标的线性表达式采用最小二乘法进行线性拟合,基于拟合结果验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性。
本发明实施例中以冬小麦为目标作物,从冬小麦长势长时间序列入手,选取高分一号16m分辨率高光谱影像作为数据源,归一化植被指数NDVI作为冬小麦长势的反演参数,基于一种时域信号分析方法——集合经验模态分解EEMD,对长时间序列进行不同时间尺度的分解,将水分、病虫害、土壤肥力以及重金属污染等多种胁迫因素分解到不同频率的IMF分量中。结合相关理论以及定期开展的野外实测工作数据结果,对不同的IMF分量进行甄别,提取出代表了水分、病虫害等瞬时胁迫的年内波动组分序列、代表冬小麦固有生长趋势的年间波动组分序列以及代表土壤肥力胁迫作用的年际波动组分序列,有效实现了土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取。结合测得的土壤养分指标数据,从统计学角度对分析结果进行验证,证明了此方法对土壤肥力胁迫进行甄别提取的可行性与准确性,并进一步确定了基于遥感影像数据参数的土壤肥力水平定量评价模型。
冬小麦长势时间序列构建:获取若干典型样区自2015年以来的高分一号卫星影像数据,选择NDVI作为冬小麦长势的反演参数,获取复合胁迫因素作用下的冬小麦长势长时间序列曲线。
EEMD分解过程实现:对于一组NDVI长时间序列,其在EEMD算法中的单次分解流程如图2(a)至图2(e)所示。其中三次样条插值法求信号上下包络线的过程可用spline函数实现,为与实际数据时间尺度相匹配,差值过程的采样频率设置为12,采样周期为1/12,即每一年内进行12次采样。
将EEMD分解过程重复N次即得到最终EEMD的分解结果。通常情况下,综合考虑对计算效率以及结果精度的影响,N的次数一般设置在50至100次之间。在这一区间内N的数值变化对结果精度的影响较小,为保证研究的科学性,本实验中设置N为100次。
在EEMD分解的整体流程中,添加高斯白噪声信号的强弱对最终结果的精度起到了决定性的影响。为选择合适的数值,研究中规定μ为加入的高斯白噪声信号标准差与原始序列间的比值,μ的数值过小则加入的噪声信号对原始序列的影响过小从而难以消除模态混叠的影响;μ过大则会影响分解结果的真实性。在一般的EEMD分解研究中,μ的数值通常设置在0.1-0.4之间,该数值的选取应结合实际研究内容考虑。
表1不同μ值下的年间波动组分相关系数
对于冬小麦而言,其不同年份间的固有生长趋势应大致接近,后文中将提到分解结果中的年间波动组分代表了冬小麦的固有生长趋势,因此通过计算比较不同μ值下各年份间年间波动组分的相关系数即可确定本实验中最佳的μ值。表1给出了不同μ值下EEMD分解结果年间波动组分各年份间的相关系数,可以看出,当μ=0.3时,平均相关系数达到最大值0.849,因此本研究中选取μ=0.3。
EEMD分解结果的组合:经过EEMD分解的NDVI时间序列根据频率差异被分解为两项周期大于年的组分、一项周期约等于年的组分、一项周期小于年的组分、一项无规则的组分以及一个残差项,根据研究需求对这些不同的分量进行组合得到三个不同时间尺度上的组分分量。通过后续甄别工作能够将这些组分归纳为以时间尺度分类的年际胁迫、年间胁迫以及年内胁迫(瞬时胁迫),同时根据不用胁迫因素的实际作用类型,可对胁迫因素进行有效的甄别。
①年内波动组合的合成:EEMD分解结果中的分量IMF1以及IMF2的频率较高,二者的周期均小于1年,因此可以判定IMF1以及IMF2为信号的年内波动组分,这一部分分量代表着植物生长周期中受到的瞬时胁迫影响作用以及信号处理过程中产生的误差结果。其中,由于信号处理过程产生的误差组分为噪声信号,这一组分不应参与模型最终分解结果的合成。因此,IMF2代表了冬小麦长势长时间序列分解结果中的年内波动组分,如图3所示。
②年间波动组分的合成:年间波动组分又可称为长时间序列分解结果中的季节项,其典型特征是信号组分具有一个明显的周期。根据农作物的实际物候参数不同,这个周期也会出现一定的变化。通常情况下,这个时间周期的长度是1年。但由于实际生产过程中耕地的利用方式以及农作物的耕作方式的影响,这个年间波动组分会呈现出周期在1年左右变动的趋势。为确定研究区域内冬小麦长势的年间波动组分,在此引入季节平均模型的方法以确定冬小麦长势的年间波动组分,其流程如下:
将长时间序列所包含年份中的每年同月的数据进行平均值求取,并重复6年以获取平均季节数据;
此平均季节数据代表了冬小麦长时间序列中所有周期特征的平均波动结果,能够反映出不同年份间冬小麦生长所共有的周期特征。利用EEMD算法对平均季节数据进行分解,得到代表该基本波动的IMF分量,如图5(a)和图5(b)所示;
将原始序列EEMD分解结果中所有与上一步中得到的IMF分量周期一致的分量进行组合,即可得到冬小麦长势的年间波动组分,如图4所示。
③年际波动组分的合成:年际波动组分又称为长时间序列分解结果的趋势项,它代表了信号在一个长时间尺度上的振动趋势。对于农作物而言,趋势项代表的意义是其受到时间尺度大于季节项的胁迫因素影响时所产生的响应。因此,将振动频率小于年频率的若干IMF分量进行合成,得到冬小麦长势长时间序列分解结果中的年际波动组分,如图6所示。
土壤肥力胁迫甄别与提取:
①年内波动组分:这部分是EEMD将原始序列分解得到的高频分量。高频分量的时间分辨率较高,频率分辨率较低,能够对原始序列的高频振荡部分进行很好地表现,可以检测信号夹带的瞬时异常现象。这些信息的特征表现为年内波动杂乱,振幅规律性不强。对相邻年份同一月份的信号数值进行相关性分析能够较好地对这一特征进行验证,在此选择3月、4月、5月、6月四个冬小麦生长较为旺盛的月份进行相关性分析,结果如表2所示。可以看出,相邻年份间的同一月份数值相关性较小,表明曲线的规律性较弱,符合上述猜想。因此,可以认为冬小麦长势长时间序列的年内波动组分代表了冬小麦受各种瞬时胁迫因素胁迫的作用结果。
表2年内波动组分相关性分析结果
相对而言,水分胁迫、病虫害胁迫等胁迫对于农作物长势的影响相较于土壤肥力、重金属污染等环境因素的响应更快,影响程度更大。而年内波动组分与其他组分相比,在曲线上表现为振动频率更大且振幅更大。这一特点能够反映出水分胁迫、病虫害胁迫的胁迫信息。因此,假定信号的高频分量代表了冬小麦受水分、病虫害等瞬时胁迫作用的响应结果。
图7(a)至图7(e)及表3给出了廊坊市大厂县2015年-2019年的年降水量。其中,2018年的年降水量远高于其他年份。根据冬小麦的耕作技术,较多的降水和连续阴雨天气会对冬小麦的长势和最终产量造成一定的影响,如影响扬花授粉、容易诱发病虫害、降低小麦粒重等。此外,从图中可以看出,大厂县降水主要集中于每年的7-9月份,结合该地区冬小麦耕作“春灌1水”的灌溉方式,水分胁迫作用在一年应在上半年和下半年各出现一个峰值。代表年内波动组分的曲线中,2018年的曲线振幅明显小于其他年份,并且曲线在每年的4月份以及9月份左右出现了两个波峰。这些特性均符合水分胁迫对冬小麦长势的影响,因此,可以认为信号的高频分量也就是年内波动组分能够代表冬小麦受水分、病虫害等瞬时胁迫作用影响的结果,将这些胁迫称为年内胁迫。
表3大厂县2015-2019年降水量
②年间波动组分:年间波动组分反映了农作物的固有生长趋势,可以近似为一个稳定的周期波动信号,如图5所示。对于试验区内的冬小麦,其在一年中有一个高峰振荡,出现在每年的5月至6月左右。理论上讲,农作物在不受外界胁迫因素影响时,不同年份的年间胁迫组分应呈现出相似的趋势,但在实际研究中,耕地作物类型的变更以及实验过程中数据获取的误差会导致年间波动组分在不用年份存在一定的差异,这种差异主要体现在波峰与波谷的数值差异上,但在时间尺度上峰值对应的时间点在不同年份内大体是一致的。
③年际波动组分:根据上文的研究结果,冬小麦长势的长时间序列由年内波动组分、年间波动组分、年内波动组分构成。在剔除代表水分胁迫、病虫害胁迫等瞬时胁迫的年内波动组分以及代表冬小麦固有生长趋势的年间波动组分后,剩余的年际波动组分可代表持续的环境胁迫对农作物长势的影响,包括土壤肥力胁迫、重金属胁迫等因素。从图6中的曲线可以看出,这种胁迫是长期而平缓的。土壤肥力胁迫的作用形式主体为农作物对土壤中养分的吸收,土壤中各种养分的含量通过不同途径的生物界循环在整体上是保持一定幅度内的平衡的,这一特征符合年际波动组分曲线平稳且振幅较小的特点。因此,可以认为年际波动组分中包含着土壤肥力胁迫作用的相应结果。同时,根据野外实测工作中利用手持XRF分析仪对几个典型样区耕地土壤中重金属元素含量的测定结果,如表4所示,所有样区耕地土壤中的重金属元素含量均在正常水平范围内,可以认为该地区受到重金属胁迫作用的影响很弱。因此,年际波动组分能够较好地反应研究区域内土壤肥力胁迫作用对于农作物长势的影响。如图8(a)至图8(f)给出了表4中几个样区的土壤肥力胁迫作用水平曲线。
表4土壤重金属污染元素含量测定结果(2020.12)
分析结果验证与土壤肥力水平综合评价模型构建:选择有机质、全氮、有效磷、速效钾四项典型养分指标进行含量测定,基于测定结果构建土壤肥力综合水平的定量评价模型。
①在SPSS软件环境中对表5所示四项养分指标测定结果进行相关性检验,得到相关性矩阵如表6所示;
表5大厂县部分地区土壤养分含量测定结果
表6土壤养分指标相关性检验结果
表7主成分分析结果
②根据相关型矩阵对四项养分指标数据进行主成分分析,结果如表7所示。前三项主成分的累计贡献率打到99.16%,基本能够代表总体水平。因此选择前三项主成分作为评价土壤肥力水平的一组新变量,如下式所示:
F=0.5759F1+0.3172F2+0.0985F3+Δ(F)
式中,F表示土壤养分总体水平;Δ(F)表示余项。
此外,新变量关于四项养分指标的线性表达式如下:
F1=0.5623X1+0.3219X2+0.0264X3+0.5444X4;
F2=-0.1352X1-0.2097X2-0.8719X3+0.4214X4;
F3=0.7921X1-0.3073X2-0.2684X3-0.4540X4.
式中,Fi表示第i个主成分;X1表示有机质含量;X2表示全氮含量;X3表示有效磷含量;X4表示速效钾含量。
将两组表达式相结合,即可得到土壤养分总体水平关于四种土壤养分指标的线性表达式,如下式所示:
F=0.359X1+0.288X2+0.149X3+0.202X4
基于上式得到几个典型样区的土壤肥力综合水平数值如表8所示。
表8大厂县样区土壤肥力水平
③与EEMD分解甄别结果进行拟合。表9列出了几个典型样区土壤肥力水平与EEMD分解的土壤肥力胁迫组分数据,利用最小二乘法对二者进行拟合,结果如表10所示。可以看出,虽然IMF11-IMF10比IMF10更能反映冬小麦长势的变化情况,但其数值变化幅度过小,导致拟合过程中的误差被大幅度提高。因此,选择单时相EEMD分解年际波动组分数据更适合对土壤肥力水平进行拟合,结果见下式:
F=0.1559+7669.2IMF-6.5*105*IMF2+2*105*IMF3
利用最小二乘法对数据进行线性拟合的方法,在实现过程上较为简单,且拟合结果的相关系数达到了0.9741,能够较好地反映出原始数据的变化水平。因此,此过程实现了土壤肥力水平基于遥感影像数据的定量评价。
表9拟合数据
表10几种不同拟合方式的拟合结果
本发明提供的基于集合经验模态分解的耕地土壤肥力水平监测方法,通过集合经验模态分解方法—EEMD,对冬小麦NDVI长时间序列进行分解,获取年内、年间、年际三个时间尺度的胁迫作用曲线,并结合野外实测数据及统计学方法确定了代表土壤肥力胁迫作用的组分,实现土壤肥力胁迫的甄别与提取;本方法在EEMD分解过程中对不同于土壤肥力时间尺度的胁迫作用进行提取和甄别,并将复杂的土壤肥力系统用一个综合评价数值进行表达,使遥感光谱数据对于其数值差异的响应程度更为明显。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,包括以如下步骤:
步骤1,通过高分一号遥感探测器采集耕地作物的光谱影像作为数据源;
步骤2,对数据源的光谱影像进行图像预处理,并将归一化植被指数NDVI作为耕地作物长势的反演参数,提取NDVI长时间序列;
步骤3,基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中;
步骤4,结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取。
2.根据权利要求1所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,所述方法还包括:结合测得的耕地土壤养分指标数据,从统计学角度对甄别提取结果进行验证,验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性。
3.根据权利要求1所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,所述步骤3,基于集合经验模态分解EEMD方法,对NDVI长时间序列进行不同时间尺度的分解,得到多个IMF分量,并将影响耕地作物长势的多种胁迫因素分解到不同的IMF分量中,具体包括:
采用集合经验模态分解EEMD将NDVI长时间序列按照年内、年间、年际三个时间尺度进行分解,得到多个IMF分量;
将水分、病虫害、土壤肥力以及重金属污染多种胁迫因素分解到不同的IMF分量中。
4.根据权利要求3所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,采用集合经验模态分解EEMD将NDVI长时间序列按照年内、年间、年际三个时间尺度进行分解,得到多个IMF分量,其中,所述集合经验模态分解EEMD进行分解的过程包括:
对一组NDVI长时间序列的原始序列,加入高斯白噪声;
采用三次样条插值法得到加入高斯白噪声后的原始序列的上下包络线及均值;
减去原始序列的均值,得到新序列;
基于新序列得到一个IMF分量;
重复上述过程,进行多次分解,得到多个IMF分量。
5.根据权利要求3所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,所述步骤4中,结合野外实测数据和历史农情数据,对不同的IMF分量进行甄别,确定代表土壤肥力胁迫作用的组分序列,实现土壤肥力胁迫作用从多种复合胁迫作用结果中的甄别提取,具体包括:
提取出代表水分、病虫害这类瞬时胁迫作用的年内波动组分序列;
提取代表耕地作物固有生长趋势的年间波动组分序列;
提取代表土壤肥力以及重金属污染这类环境胁迫作用的年际波动组分序列。
6.根据权利要求2所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,所述结合测得的耕地土壤养分指标数据,从统计学角度对甄别提取结果进行验证,验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性,具体包括:
在SPSS软件环境中对有机质、全氮、有效磷、速效钾四项养分指标数据测定结果进行相关性检验,得到相关性矩;
根据相关型矩阵对四项养分指标数据进行主成分分析,得到土壤养分总体水平关于四项养分指标的线性表达式;
将甄别提取结果与所述土壤养分总体水平关于四项养分指标的线性表达式采用最小二乘法进行线性拟合,基于拟合结果验证所述方法对土壤肥力胁迫作用进行甄别提取的可行性与准确性。
7.根据权利要求6所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,根据相关型矩阵对四项养分指标数据进行主成分分析,得到土壤养分总体水平关于四项养分指标的线性表达式,具体包括:
根据相关型矩阵对四项养分指标数据进行主成分分析,选择前三项主成分作为评价土壤肥力水平的一组新变量,如下式所示:
F=0.5759F1+0.3172F2+0.0985F3+Δ(F)
式中,F表示土壤养分总体水平;Fi表示第i个主成分,Δ(F)表示余项;
基于新变量新变量关于四项养分指标的线性表达式,得到土壤养分总体水平关于四种土壤养分指标的线性表达式;
新变量关于四项养分指标的线性表达式如下所示:
F1=0.5623X1+0.3219X2+0.0264X3+0.5444X4;
F2=-0.1352X1-0.2097X2-0.8719X3+0.4214X4;
F3=0.7921X1-0.3073X2-0.2684X3-0.4540X4;
式中,X1表示有机质含量;X2表示全氮含量;X3表示有效磷含量;X4表示速效钾含量;
土壤养分总体水平关于四种土壤养分指标的线性表达式如下:
F=0.359X1+0.288X2+0.149X3+0.202X4。
8.根据权利要求1所述的基于集合经验模态分解的耕地土壤肥力水平监测方法,其特征在于,所述耕地作物选择冬小麦。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110418313.9A CN113111799B (zh) | 2021-04-19 | 2021-04-19 | 一种基于集合经验模态分解的耕地土壤肥力水平监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110418313.9A CN113111799B (zh) | 2021-04-19 | 2021-04-19 | 一种基于集合经验模态分解的耕地土壤肥力水平监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113111799A true CN113111799A (zh) | 2021-07-13 |
CN113111799B CN113111799B (zh) | 2024-01-30 |
Family
ID=76718311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110418313.9A Active CN113111799B (zh) | 2021-04-19 | 2021-04-19 | 一种基于集合经验模态分解的耕地土壤肥力水平监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113111799B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116663779A (zh) * | 2023-07-27 | 2023-08-29 | 北京市农林科学院信息技术研究中心 | 耕地土壤多深度肥力点面转换方法及装置 |
WO2024099475A1 (zh) * | 2022-11-07 | 2024-05-16 | 北华航天工业学院 | 基于长时序ndvi的土壤重金属胁迫甄别方法 |
WO2024098587A1 (zh) * | 2022-11-10 | 2024-05-16 | 北华航天工业学院 | 基于ceemdan的农田土壤水分胁迫甄别与监测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073869A (zh) * | 2010-12-27 | 2011-05-25 | 中国农业大学 | 基于点位数据和遥感影像数据的区域耕地质量监测方法 |
CN106918816A (zh) * | 2017-03-21 | 2017-07-04 | 中国科学院遥感与数字地球研究所 | 农作物长势监测系统及其数据处理方法和模块 |
US20170228475A1 (en) * | 2016-02-05 | 2017-08-10 | The Climate Corporation | Modeling trends in crop yields |
CN108388866A (zh) * | 2018-02-27 | 2018-08-10 | 海南热带海洋学院 | 一种植被单调变化趋势检测方法及相关装置 |
CN109635731A (zh) * | 2018-12-12 | 2019-04-16 | 中国科学院深圳先进技术研究院 | 一种识别有效耕地的方法及装置、存储介质及处理器 |
CN110046415A (zh) * | 2019-04-08 | 2019-07-23 | 中国科学院南京地理与湖泊研究所 | 一种时空精细化的土壤有机质含量遥感动态反演方法 |
CN110398466A (zh) * | 2019-08-05 | 2019-11-01 | 北京绿土科技有限公司 | 基于遥感反演的农作物生长状态监测方法 |
US20210027429A1 (en) * | 2019-07-26 | 2021-01-28 | Zhejiang University Of Technology | Noise detection method for time-series vegetation index derived from remote sensing images |
-
2021
- 2021-04-19 CN CN202110418313.9A patent/CN113111799B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073869A (zh) * | 2010-12-27 | 2011-05-25 | 中国农业大学 | 基于点位数据和遥感影像数据的区域耕地质量监测方法 |
US20170228475A1 (en) * | 2016-02-05 | 2017-08-10 | The Climate Corporation | Modeling trends in crop yields |
CN106918816A (zh) * | 2017-03-21 | 2017-07-04 | 中国科学院遥感与数字地球研究所 | 农作物长势监测系统及其数据处理方法和模块 |
CN108388866A (zh) * | 2018-02-27 | 2018-08-10 | 海南热带海洋学院 | 一种植被单调变化趋势检测方法及相关装置 |
CN109635731A (zh) * | 2018-12-12 | 2019-04-16 | 中国科学院深圳先进技术研究院 | 一种识别有效耕地的方法及装置、存储介质及处理器 |
CN110046415A (zh) * | 2019-04-08 | 2019-07-23 | 中国科学院南京地理与湖泊研究所 | 一种时空精细化的土壤有机质含量遥感动态反演方法 |
US20210027429A1 (en) * | 2019-07-26 | 2021-01-28 | Zhejiang University Of Technology | Noise detection method for time-series vegetation index derived from remote sensing images |
CN110398466A (zh) * | 2019-08-05 | 2019-11-01 | 北京绿土科技有限公司 | 基于遥感反演的农作物生长状态监测方法 |
Non-Patent Citations (5)
Title |
---|
PIETER HAWINKEL: "A time series processing tool to extract climate-driven interannual vegetation dynamics using Ensemble Empirical Mode Decomposition (EEMD)", 《REMOTE SENSING OF ENVIRONMENT》, vol. 169, pages 375 - 389 * |
刘帝: "基于时空谱分析的黄淮海平原北部土壤肥力胁迫甄别与监测研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, pages 028 - 112 * |
李旭青等: "基于集合经验模态分解算法的土壤肥力胁迫甄别与监测", 《农业工程学报》, pages 137 - 146 * |
李燕丽;潘贤章;周睿;王昌昆;刘娅;石荣杰;陈冬峰;赵其国;: "长期土壤肥力因子变化及其与植被指数耦合关系", 生态学杂志, no. 03, pages 536 - 541 * |
范倩倩;赵安周;王金杰;张安兵;裴韬;: "1982―2015年黄土高原NDVI时空演变及其对气候变化的季节响应", 生态学杂志, no. 05, pages 1664 - 1675 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024099475A1 (zh) * | 2022-11-07 | 2024-05-16 | 北华航天工业学院 | 基于长时序ndvi的土壤重金属胁迫甄别方法 |
WO2024098587A1 (zh) * | 2022-11-10 | 2024-05-16 | 北华航天工业学院 | 基于ceemdan的农田土壤水分胁迫甄别与监测方法 |
CN116663779A (zh) * | 2023-07-27 | 2023-08-29 | 北京市农林科学院信息技术研究中心 | 耕地土壤多深度肥力点面转换方法及装置 |
CN116663779B (zh) * | 2023-07-27 | 2023-10-20 | 北京市农林科学院信息技术研究中心 | 耕地土壤多深度肥力点面转换方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113111799B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110082300B (zh) | 基于光谱参量的冬小麦冠层含水率监测模型建立方法 | |
CN113111799B (zh) | 一种基于集合经验模态分解的耕地土壤肥力水平监测方法 | |
Li et al. | Evaluating chlorophyll density in winter oilseed rape (Brassica napus L.) using canopy hyperspectral red-edge parameters | |
Kowalik et al. | Yield estimation using SPOT-VEGETATION products: A case study of wheat in European countries | |
Lee et al. | Estimating chemical properties of Florida soils using spectral reflectance | |
CN108985260A (zh) | 一种遥感及气象一体化水稻估产方法 | |
CN109711102B (zh) | 一种作物灾害损失快速评估方法 | |
CN108169161A (zh) | 一种基于改进型modis指数的玉米种植区域土壤湿度评估方法 | |
CN111044516B (zh) | 一种水稻叶绿素含量遥感估测方法 | |
Hu et al. | Evapotranspiration partitioning and its interannual variability over a winter wheat-summer maize rotation system in the North China Plain | |
Peng et al. | Winter wheat canopy water content monitoring based on spectral transforms and “three-edge” parameters | |
Chen et al. | Including soil water stress in process-based ecosystem models by scaling down maximum carboxylation rate using accumulated soil water deficit | |
Mittermayer et al. | Analysis of site-specific N balances in heterogeneous croplands using digital methods | |
Tian et al. | Extracting red edge position parameters from ground-and space-based hyperspectral data for estimation of canopy leaf nitrogen concentration in rice | |
Liu et al. | Winter wheat yield estimation based on assimilated Sentinel-2 images with the CERES-Wheat model | |
Wang et al. | Improving remote sensing of salinity on topsoil with crop residues using novel indices of optical and microwave bands | |
Sun et al. | Monitoring leaf nitrogen accumulation with optimized spectral index in winter wheat under different irrigation regimes | |
CN111537510B (zh) | 一种基于空间信息技术的农田防护林防风效应计量方法 | |
An et al. | The optimal soil water content models based on crop-LAI and hyperspectral data of winter wheat | |
US20240167947A1 (en) | CEEMDAN-based method for screening and monitoring soil moisture stress in agricultural fields | |
CN109765247B (zh) | 一种不同生育阶段麦类作物覆盖区域土壤水分反演方法 | |
Feng et al. | Monitoring leaf pigment status with hyperspectral remote sensing in wheat | |
Tola et al. | Determining soil organic carbon concentration in agricultural fields using a handheld spectroradiometer: Implication for soil fertility measurement | |
CN116128353A (zh) | 一种基于遥感生物量反演的优势树种BVOCs排放估算方法 | |
Li et al. | Developing a USLE cover and management factor (C) for forested regions of southern China |
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 |