CN111855592A - 水稻单位面积内地上部干物重遥感反演模型和方法 - Google Patents
水稻单位面积内地上部干物重遥感反演模型和方法 Download PDFInfo
- Publication number
- CN111855592A CN111855592A CN202010772402.9A CN202010772402A CN111855592A CN 111855592 A CN111855592 A CN 111855592A CN 202010772402 A CN202010772402 A CN 202010772402A CN 111855592 A CN111855592 A CN 111855592A
- Authority
- CN
- China
- Prior art keywords
- rice
- unit area
- dry matter
- model
- matter weight
- 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.)
- Pending
Links
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 191
- 235000009566 rice Nutrition 0.000 title claims abstract description 191
- 238000000034 method Methods 0.000 title claims abstract description 54
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 190
- 238000002310 reflectometry Methods 0.000 claims description 36
- 238000012549 training Methods 0.000 claims description 31
- 238000005070 sampling Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000005259 measurement Methods 0.000 claims description 19
- 239000000523 sample Substances 0.000 claims description 15
- 238000001035 drying Methods 0.000 claims description 8
- 241000196324 Embryophyta Species 0.000 claims description 7
- 102000004190 Enzymes Human genes 0.000 claims description 7
- 108090000790 Enzymes Proteins 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 7
- 239000012634 fragment Substances 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 abstract description 26
- 230000000694 effects Effects 0.000 abstract description 12
- 238000013461 design Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 8
- 238000003066 decision tree Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 230000008901 benefit Effects 0.000 description 7
- 238000007689 inspection Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 239000003337 fertilizer Substances 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000009849 deactivation Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000000985 reflectance spectrum Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000029553 photosynthesis Effects 0.000 description 2
- 238000010672 photosynthesis Methods 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241001483116 Neopicrorhiza scrophulariiflora Species 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003851 biochemical process Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 230000004962 physiological condition Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G17/00—Apparatus for or methods of weighing material of special form or property
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/55—Specular reflectivity
-
- 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/0098—Plants or trees
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Wood Science & Technology (AREA)
- Algebra (AREA)
- Botany (AREA)
- Databases & Information Systems (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种水稻单位面积内地上部干物重遥感反演模型，水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，并进一步提供了该极端随机树模型的模型参数。还提供了水稻单位面积内地上部干物重遥感反演方法。本发明的水稻单位面积内地上部干物重遥感反演模型能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演模型精度，设计巧妙，计算简便，易于实现，成本低，适于大规模推广应用。
Description
技术领域
本发明涉及农业遥感技术领域，特别涉及水稻单位面积内地上部干物重测量技术领域，具体是指一种水稻单位面积内地上部干物重遥感反演模型和方法。
背景技术
水稻地上部干物质的积累是水稻净光合作用的重要标志，受水稻光合能力、环境温度以及肥水等因素的影响，反映了水稻生理、长势以及肥水状况(胡继超,姜东,曹卫星,等.短期干旱对水稻叶水势、光合作用及干物质分配的影响[J].应用生态学报,2004,15(1):63-67)。
监测水稻单位面积内地上部干物重，不仅可以保证水稻生产的产量和品质，同时能够动态管理水稻的水肥施用，减少水稻生产中水肥的使用量，从而产生显著的经济和社会效益(薛利红,曹卫星,罗卫红.基于冠层反射光谱的水稻产量预测模型[J].遥感学报,2005,9(001):100-105)。传统的水稻单位面积内地上部干物重的监测主要采用破坏性采样的方法，需要在室内进行烘干称重，费时费力，无法及时的获取水稻的单位面积内地上部干物重。
水稻的生理生化过程中，体内某些特定物质和细胞结构的变化，导致了水稻反射光谱的改变。因此，可以用光谱的变化来获取水稻生长信息(周冬琴.基于冠层反射光谱的水稻氮素营养与籽粒品质监测[D].南京农业大学,2007)。目前，使用高光谱来监测水稻的生长状态已经在作物生长监测中得到了应用。随着光谱技术的发展和普及，使用光谱可以快捷迅速的获取水稻单位面积内地上部干物重信息已经成为越来越多水稻生产从业者和研究人员的共识。最常用的方式是使用便携式全波段光谱仪获取水稻生长信息，选择能够反映地上部干物质积累量(单位面积内地上部干物重)的特征波段构建反演模型。在构建水稻单位面积内地上部干物重反演模型的过程中，全波段光谱仪所测定的光谱范围涵盖了350nm～2500nm，水稻组分复杂，组分光谱特征波段部分重叠，水稻单位面积内地上部干物重特征光谱的确定困难，同时高光谱数据的快速处理成为基于高光谱数据估测水稻单位面积内地上部干物重的亟待解决的技术问题。
因此，希望提供一种水稻单位面积内地上部干物重遥感反演模型，其能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演模型精度。
发明内容
为了克服上述现有技术中的缺点，本发明的一个目的在于提供一种水稻单位面积内地上部干物重遥感反演模型，其能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演模型精度，适于大规模推广应用。
本发明的另一目的在于提供一种水稻单位面积内地上部干物重遥感反演模型，其设计巧妙，计算简便，易于实现，成本低，适于大规模推广应用。
本发明的另一目的在于提供一种水稻单位面积内地上部干物重遥感反演方法，其能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演精度，适于大规模推广应用。
本发明的另一目的在于提供一种水稻单位面积内地上部干物重遥感反演方法，其设计巧妙，操作简单方便，成本低，适于大规模推广应用。
为达到以上目的，在本发明的第一方面，提供一种水稻单位面积内地上部干物重遥感反演模型，其特点是，所述水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，所述极端随机树模型的模型参数为：'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
较佳地，所述极端随机树模型采用水稻的数据集训练而成，所述数据集包括所述水稻的m个样点的冠层反射率和单位面积内地上部干物重，m个所述样点均匀分布在水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率。
更佳地，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
在本发明的第二方面，提供了一种水稻单位面积内地上部干物重遥感反演方法，其特点是，包括以下步骤：
(1)测量水稻的冠层反射率；
(2)测量所述水稻的单位面积内地上部干物重：
(3)以所述冠层反射率为输入数据，采用Python语言的极端随机树模型进行计算，获得反演值，根据所述反演值与所述单位面积内地上部干物重计算决定系数R2，改变所述极端随机树模型的模型参数的取值，R2的变化越大，说明所述模型参数的重要性越大，将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵；
(4)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为输出结果，训练所述极端随机树模型，根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优，获得所述模型参数的调优值；
(5)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为所述输出结果，采用所述的模型参数的调优值，训练所述极端随机树模型，待所述极端随机树模型训练结束后，获得水稻单位面积内地上部干物重遥感反演模型，使用save方法保存所述水稻单位面积内地上部干物重遥感反演模型，如果需要使用所述水稻单位面积内地上部干物重遥感反演模型，使用load方法加载所述水稻单位面积内地上部干物重遥感反演模型使用。
较佳地，在所述步骤(1)中，所述测量采用高光谱辐射仪进行，所述测量的时间为10:00～14:00，所述高光谱辐射仪采用视场角为25°的镜头，所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米，所述传感器探头的地面视场范围直径为0.44米，所述传感器探头迎向阳光，所述测量采用标准板进行校正，所述标准板是反射率为95％～99％的标准白板。
较佳地，在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体包括：
采集所述水稻的地上部分植株，杀青，烘干，至恒重后测定重量，获得地上部干物重，根据采样覆盖面积，获得所述单位面积内地上部干物重。
更佳地，在所述步骤(2)中，所述杀青的温度为105℃，所述杀青的时间为20分钟～30分钟，所述烘干的温度为80℃～90℃。
较佳地，在所述步骤(3)中，所述模型参数调优秩次矩阵为：
Params＝{'min_impurity_split','ccp_alpha','max_depth','max_features','min_impurity_decrease','min_samples_leaf','min_samples_split','min_weight_fraction_leaf','splitter','criterion','max_leaf_nodes'}。
更佳地，在所述步骤(4)中，所述的模型参数的调优值为：
'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
较佳地，在所述步骤(1)中，所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率，m个所述样点均匀分布在所述水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率；在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体为测量该m个所述样点的所述单位面积内地上部干物重。
更佳地，在所述步骤(1)中，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
本发明的有益效果主要在于：
1、本发明的水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，极端随机树模型的模型参数为：'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None，对该模型进行检验，R2在0.85以上，因此，其能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演模型精度，适于大规模推广应用。
2、本发明的水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，极端随机树模型的模型参数为：'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None，对该模型进行检验，R2在0.85以上，因此，其设计巧妙，计算简便，易于实现，成本低，适于大规模推广应用。
3、本发明的水稻单位面积内地上部干物重遥感反演方法，包括：测量水稻的冠层反射率；测量水稻的单位面积内地上部干物重：以冠层反射率为输入数据，采用Python语言的极端随机树模型进行计算，以决定系数R2构建模型参数调优秩次矩阵；以冠层反射率为输入数据，以单位面积内地上部干物重为输出结果，训练极端随机树模型，根据模型参数调优秩次矩阵依次对模型参数进行调优，获得模型参数的调优值；以冠层反射率为输入数据，以单位面积内地上部干物重为输出结果，采用模型参数的调优值，训练极端随机树模型，获得水稻单位面积内地上部干物重遥感反演模型，对该模型进行检验，R2在0.85以上，因此，其能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演精度，适于大规模推广应用。
4、本发明的水稻单位面积内地上部干物重遥感反演方法，包括：测量水稻的冠层反射率；测量水稻的单位面积内地上部干物重：以冠层反射率为输入数据，采用Python语言的极端随机树模型进行计算，以决定系数R2构建模型参数调优秩次矩阵；以冠层反射率为输入数据，以单位面积内地上部干物重为输出结果，训练极端随机树模型，根据模型参数调优秩次矩阵依次对模型参数进行调优，获得模型参数的调优值；以冠层反射率为输入数据，以单位面积内地上部干物重为输出结果，采用模型参数的调优值，训练极端随机树模型，获得水稻单位面积内地上部干物重遥感反演模型，对该模型进行检验，R2在0.85以上，因此，其设计巧妙，操作简单方便，成本低，适于大规模推广应用。
本发明的这些和其它目的、特点和优势，通过下述的详细说明、附图和权利要求得以充分体现，并可通过所附权利要求中特地指出的手段、装置和它们的组合得以实现。
附图说明
图1是本发明的水稻单位面积内地上部干物重遥感反演方法的一具体实施例的流程示意图。
图2是图1所示的具体实施例的模型构建流程示意图。
图3是图1所示的具体实施例的模型检验结果示意图，其中实测值和预测值的单位均为g/m2。
具体实施方式
本发明人针对基于高光谱估测水稻单位面积内地上部干物重的需求，克服了因水稻组分复杂而导致的水稻单位面积内地上部干物重特征波段难以确定以及高光谱数据特征波段筛选费时费力的困难，提出了一种水稻单位面积内地上部干物重遥感反演模型，所述水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，所述极端随机树模型的模型参数为：'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
所述极端随机树模型可以采用任何合适的数据集训练而成，较佳地，所述极端随机树模型采用水稻的数据集训练而成，所述数据集包括所述水稻的m个样点的冠层反射率和单位面积内地上部干物重，m个所述样点均匀分布在水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率。所述水稻种植区域可以是多个生态点、多个品种的水稻种植区域。
所述m和所述n为正整数，可以根据需要确定，更佳地，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
本发明还提供了一种水稻单位面积内地上部干物重遥感反演方法，包括以下步骤：
(1)测量水稻的冠层反射率；
(2)测量所述水稻的单位面积内地上部干物重：
(3)以所述冠层反射率为输入数据，采用Python语言的极端随机树模型进行计算，获得反演值，根据所述反演值与所述单位面积内地上部干物重计算决定系数R2，改变所述极端随机树模型的模型参数的取值，R2的变化越大，说明所述模型参数的重要性越大，将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵；
(4)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为输出结果，训练所述极端随机树模型，根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优，获得所述模型参数的调优值；
(5)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为所述输出结果，采用所述的模型参数的调优值，训练所述极端随机树模型，待所述极端随机树模型训练结束后，获得水稻单位面积内地上部干物重遥感反演模型，使用save方法保存所述水稻单位面积内地上部干物重遥感反演模型，如果需要使用所述水稻单位面积内地上部干物重遥感反演模型，使用load方法加载所述水稻单位面积内地上部干物重遥感反演模型使用。
在所述步骤(1)中，所述测量可以采用任何合适的光谱仪和方法进行，较佳地，在所述步骤(1)中，所述测量采用高光谱辐射仪进行，所述测量的时间为10:00～14:00，所述高光谱辐射仪采用视场角为25°的镜头，所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米，所述传感器探头的地面视场范围直径为0.44米，所述传感器探头迎向阳光，所述测量采用标准板进行校正，所述标准板是反射率为95％～99％的标准白板。
在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤可以具体包括任何合适的方法，较佳地，在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体包括：
采集所述水稻的地上部分植株，杀青，烘干，至恒重后测定重量，获得地上部干物重，根据采样覆盖面积，获得所述单位面积内地上部干物重。
在所述步骤(2)中，所述杀青和所述烘干可以采用任何合适的条件，更佳地，在所述步骤(2)中，所述杀青的温度为105℃，所述杀青的时间为20分钟～30分钟，所述烘干的温度为80℃～90℃。
在所述步骤(3)中，所述模型参数调优秩次矩阵根据决定系数R2确定，较佳地，在所述步骤(3)中，所述模型参数调优秩次矩阵为：
Params＝{'min_impurity_split','ccp_alpha','max_depth','max_features','min_impurity_decrease','min_samples_leaf','min_samples_split','min_weight_fraction_leaf','splitter','criterion','max_leaf_nodes'}。
在所述步骤(4)中，所述的模型参数的调优值根据所述模型参数调优秩次矩阵依次确定，更佳地，在所述步骤(4)中，所述的模型参数的调优值为：
'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
为了使得所述水稻单位面积内地上部干物重遥感反演模型的精度更高，可以选取水稻种植区域的多个样点，测定多个样点的多个特征波段的冠层反射率以及多个样点的单位面积内地上部干物重，较佳地，在所述步骤(1)中，所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率，m个所述样点均匀分布在所述水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率；在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体为测量该m个所述样点的所述单位面积内地上部干物重。
在所述步骤(1)中，所述m和所述n为正整数，可以根据需要确定，更佳地，在所述步骤(1)中，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
下面结合具体实施例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。
实施例
本实施例的水稻单位面积内地上部干物重遥感反演方法基于实测高光谱数据，采用了水稻种植区域(江苏省淮安市农业科学研究院淮安地区稻麦种植基地，水稻品种为淮稻5号，采样时期为水稻拔节期)采集的水稻冠层反射率光谱数据以及水稻单位面积内地上部干物重数据，共48个采样点，这些采样点均匀分布并完全覆盖水稻种植区域整个区域。48个采样点数据采用随机的方法分为两部分，其中37个采样点的数据用于模型构建，11个采样点数据用于模型检验。水稻单位面积内地上部干物重遥感反演方法的流程如图1所示，包括如下步骤：
1.光谱测定。
水稻冠层光谱测定使用美国ASD生产的FieldSpec Pro便携式野外高光谱辐射仪，选择在天气晴朗、无风或者风速很小时进行，时间范围10:00～14:00，采样测试人员身着深色服饰，减少对光谱仪的影响或者干扰。采样时，选用视场角为25°的镜头，传感器探头垂直指向测量目标---冠层，距离冠层顶层垂直高度约1米，地面视场范围直径为0.44米，迎向阳光，取测量10次的反射光谱平均值，作为该采样点的光谱数据。测量过程中，对每个采样点的测量前后进行标准白板校正。如果测试过程中，环境光场分布发生变化，亦进行标准白板校正，本例使用的标准白板的反射率为99％。测定的光谱数据使用FieldSpec Pro便携式野外高光谱辐射仪随机软件RS3或者ViewSpec Pro软件检查，剔除异常光谱文件，对光谱数据进行插值计算，得到范围为350nm～2500nm，分辨率为1nm的光谱数据，计算光谱平行采样光谱的平均值，最后导出光谱数据并保存为ASCII文件。
2.水稻单位面积内地上部干物重测定
收集均匀分布于每个采样点的光谱测量视场的水稻地上部分植株，数量为6株，用吸水纸包好，带回实验室，先于105℃下杀青20分钟，然后在85℃下烘干，至恒重后测定重量，所得数据即为水稻地上部干物重，根据采样覆盖面积，获得水稻单位面积内地上部干物重，单位为g/m2。
3.模型构建
模型构建采用Python语言的极端随机树模型进行构建，请参见图2所示，模型构建主要包括以下步骤：
3.1数据校验
对获取的水稻冠层反射率数据进行校验，剔除异常的整条光谱曲线数据。本发明中异常光谱是指相邻光谱变化超过100％，且包含空值、负值的光谱值。
3.2数据的预处理
对校验后的水稻冠层反射率数据以及水稻单位面积内地上部干物重数据进行预处理，包括去除包含缺失值和空值的成对水稻冠层反射率数据和水稻单位面积内地上部干物重数据。
3.3数据集的划分
为保证模型训练和反演结果的合理评估，使用随机的方法，将整个数据集(48组)划分为两部分，80％(37组)的数据用于模型训练，20％(11组)的数据用于训练后的效果评估。
3.4训练数据集的划分
为保证模型训练的效果，使用随机的方法，在每次模型训练迭代时，将训练数据集划分为5个部分，对模型进行训练。
3.5模型参数调优秩次矩阵的构建
本发明中，模型训练过程模型参数的调优非常重要，为保证尽量获取最佳的模型调优，使用试错法进行模型参数调优。本发明使用决定系数R2(R2越接近于1越好)作为检验参数，构建了评估模型参数权重的参数秩次矩阵。根据训练数据集首先使用模型参数的默认值进行计算，获取反演值，根据反演值和单位面积内地上部干物重，计算决定系数R2，然后改变模型参数的取值，R2的变化越大，说明模型参数的重要性越大，将模型参数按照重要性的大小由大到小进行排列构建模型参数调优秩次矩阵，用于后续计算。
根据训练数据集中冠层反射率数据和对应的单位面积内地上部干物重数据，计算得到的模型参数调优秩次矩阵为：
Params＝{'min_impurity_split','ccp_alpha','max_depth','max_features','min_impurity_decrease','min_samples_leaf','min_samples_split','min_weight_fraction_leaf','splitter','criterion','max_leaf_nodes'}。
其中max_leaf_nodes变化不会引起模型的精度变化。
3.6模型构建
依据获取的模型参数调优秩次矩阵，将建模使用的数据，包括实测的冠层反射率数据和对应的实测的单位面积内地上部干物重数据，以实测的冠层反射率数据为输入数据，实测的单位面积内地上部干物重数据为输出结果，训练极端随机树模型，根据模型参数调优秩次矩阵依次对模型参数进行调优，获得该模型完整参数及取值，如下：
'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
待模型训练结束后，使用save方法保存模型，如需使用，运行load方法加载使用。
对于包含m个样本以及n个特征波段的数据集，Python语言的极端随机树模型的模型构建计算过程如下：
(1)使用所有的训练样本构建多个决策树；
(2)构建决策树时，根据评估分数使用分数最好的特征来构建决策树；
(3)在构建决策树的时候，在特征经验范围内均匀随机生成决策树的分叉值，在所有随机的划分点中，选择其中分数最高的作为结点的划分点，不限制决策树的深度；
(4)在训练结束之后，对未知样本x的预测可以通过对x上所有单个回归树的预测求平均来实现：
其中，
为最终的预测值，B为构建的决策树个数，fb为构建的单个决策树，x为样本数据。
3.7模型检验
使用构建模型之外的11个采样点高光谱数据输入模型，使用调优后的模型参数进行计算，得到预测值，分析预测值与实测值的关系，结果如图3所示，模型的R2为0.8786。而使用默认参数的模型R2为0。
本案例训练数据和检验数据的随机划分以及模型的构建、训练和检验使用了美国MathWorks公司开发的Matlab软件(version:R2020a 9.8.0.1380330)以及Python(version:3.7.0)，通过Matlab软件调用Python的极端随机树模型。
因此，本发明提出了一种新的、基于实测高光谱遥感数据的水稻单位面积内地上部干物重遥感反演模型，基于实测水稻冠层反射率数据和实地采集的水稻单位面积内地上部干物重数据，可快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，并通过构建模型参数调优秩次矩阵，使用试错法对模型参数调优，有效降低了线性模型过拟合的现象，大大提高了水稻单位面积内地上部干物重反演的精度，适用于不同生态区、不同品种以及主要生育期的水稻单位面积内地上部干物重的定量反演，从而获得水稻生理状态以及水肥供给的状态，提高了水稻栽培种植过程中的生长信息获取效率，并为水稻生产中水分肥料运筹提供了基础科学数据。
与现有技术相比，本发明具有以下优点：
(1)本发明使用的极端随机树模型(Extremely randomized trees model,ET)适用于基于高光谱的水稻单位面积内地上部干物重的反演，在综合考虑高光谱350～2500nm波段范围信息的基础上，考虑了水稻体内各种物质组成和细胞结构的光学特性，特别是复杂的组分对水稻单位面积内地上部干物重特征波段的影响和叠加效应，充分利用遥感数据中不同波段包含的水稻单位面积内地上部干物重信息，进行水稻单位面积内地上部干物重的反演；
(2)使用了极端随机树的机器学习算法，构建350～2500nm反射率与水稻单位面积内地上部干物重的模型，可有效减少使用线性回归等模型出现的过拟合现象，提高基于高光谱信息的水稻单位面积内地上部干物重反演的速度和效率；
(3)本发明充分考虑了模型训练和模型检验的独立性，使用随机分割的方法划分训练数据集和检验数据集，训练数据集只用于模型训练，检验数据集仅用于模型检验，以保证模型效果检验的合理性。
(4)由于模型的参数调优对于模型的计算精度非常重要，因此本发明中构建了模型参数秩次矩阵，以决定系数R2为评价参数，使用试错法进行模型参数调优，在保证参数调优效果的基础上，大大提高了模型训练和参数调优的速度。
(5)本发明提出的水稻地上部反演方法计算简便，适用于不同生态区、不同品种以及不同生育期的水稻单位面积内地上部干物重的遥感定量反演，可准确的反演水稻单位面积内地上部干物重，快速的获取水稻生理状况和长势等信息，同时为水稻种植栽培的水分肥料运筹管理提供科学数据。
综上，本发明的水稻单位面积内地上部干物重遥感反演模型能够快速、准确的获取水稻单位面积内地上部干物重信息，克服水稻组分复杂带来的光谱叠加效应造成的水稻单位面积内地上部干物重特征波段难以确定的困难，大大提高水稻单位面积内地上部干物重反演模型精度，设计巧妙，计算简便，易于实现，成本低，适于大规模推广应用。
由此可见，本发明的目的已经完整并有效的予以实现。本发明的功能及结构原理已在实施例中予以展示和说明，在不背离所述原理下，实施方式可作任意修改。所以，本发明包括了基于权利要求精神及权利要求范围的所有变形实施方式。
Claims (11)
1.一种水稻单位面积内地上部干物重遥感反演模型，其特征在于，所述水稻单位面积内地上部干物重遥感反演模型为Python语言的极端随机树模型，所述极端随机树模型的模型参数为：'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
2.如权利要求1所述的水稻单位面积内地上部干物重遥感反演模型，其特征在于，所述极端随机树模型采用水稻的数据集训练而成，所述数据集包括所述水稻的m个样点的冠层反射率和单位面积内地上部干物重，m个所述样点均匀分布在水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率。
3.如权利要求2所述的水稻单位面积内地上部干物重遥感反演模型，其特征在于，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
4.一种水稻单位面积内地上部干物重遥感反演方法，其特征在于，包括以下步骤：
(1)测量水稻的冠层反射率；
(2)测量所述水稻的单位面积内地上部干物重：
(3)以所述冠层反射率为输入数据，采用Python语言的极端随机树模型进行计算，获得反演值，根据所述反演值与所述单位面积内地上部干物重计算决定系数R2，改变所述极端随机树模型的模型参数的取值，R2的变化越大，说明所述模型参数的重要性越大，将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵；
(4)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为输出结果，训练所述极端随机树模型，根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优，获得所述模型参数的调优值；
(5)以所述冠层反射率为所述输入数据，以所述单位面积内地上部干物重为所述输出结果，采用所述的模型参数的调优值，训练所述极端随机树模型，待所述极端随机树模型训练结束后，获得水稻单位面积内地上部干物重遥感反演模型，使用save方法保存所述水稻单位面积内地上部干物重遥感反演模型，如果需要使用所述水稻单位面积内地上部干物重遥感反演模型，使用load方法加载所述水稻单位面积内地上部干物重遥感反演模型使用。
5.如权利要求4所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(1)中，所述测量采用高光谱辐射仪进行，所述测量的时间为10:00～14:00，所述高光谱辐射仪采用视场角为25°的镜头，所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米，所述传感器探头的地面视场范围直径为0.44米，所述传感器探头迎向阳光，所述测量采用标准板进行校正，所述标准板是反射率为95％～99％的标准白板。
6.如权利要求4所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体包括：
采集所述水稻的地上部分植株，杀青，烘干，至恒重后测定重量，获得地上部干物重，根据采样覆盖面积，获得所述单位面积内地上部干物重。
7.如权利要求6所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(2)中，所述杀青的温度为105℃，所述杀青的时间为20分钟～30分钟，所述烘干的温度为80℃～90℃。
8.如权利要求4所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(3)中，所述模型参数调优秩次矩阵为：
Params＝{'min_impurity_split','ccp_alpha','max_depth','max_features','min_impurity_decrease','min_samples_leaf','min_samples_split','min_weight_fraction_leaf','splitter','criterion','max_leaf_nodes'}。
9.如权利要求8所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(4)中，所述的模型参数的调优值为：
'min_impurity_split':0.0046954861925470655,'ccp_alpha':0.00047054761925470663,'max_depth':45,'max_features':'sqrt','min_impurity_decrease':0.0009710127579306127,'min_samples_leaf':1,'min_samples_split':2,'min_weight_fraction_leaf':0.010209198720162859,'splitter':'random','criterion':'mse','max_leaf_nodes':None。
10.如权利要求4所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(1)中，所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率，m个所述样点均匀分布在所述水稻种植区域，所述冠层反射率为n个特征波段的冠层反射率；在所述步骤(2)中，所述的测量所述水稻的单位面积内地上部干物重的步骤具体为测量该m个所述样点的所述单位面积内地上部干物重。
11.如权利要求10所述的水稻单位面积内地上部干物重遥感反演方法，其特征在于，在所述步骤(1)中，所述m为37，所述n个特征波段为2151个特征波段，所述2151个特征波段为从350nm波段至2500nm波段。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010772402.9A CN111855592A (zh) | 2020-08-04 | 2020-08-04 | 水稻单位面积内地上部干物重遥感反演模型和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010772402.9A CN111855592A (zh) | 2020-08-04 | 2020-08-04 | 水稻单位面积内地上部干物重遥感反演模型和方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111855592A true CN111855592A (zh) | 2020-10-30 |
Family
ID=72953156
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010772402.9A Pending CN111855592A (zh) | 2020-08-04 | 2020-08-04 | 水稻单位面积内地上部干物重遥感反演模型和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111855592A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112270132A (zh) * | 2020-11-11 | 2021-01-26 | 仲恺农业工程学院 | 基于高斯过程回归算法的水稻单位面积内地上部干物重遥感反演模型和方法 |
| CN112362659A (zh) * | 2020-11-26 | 2021-02-12 | 淮阴师范学院 | 基于贝叶斯岭回归算法的水稻地上部糖氮比遥感反演模型和方法 |
| CN112630161A (zh) * | 2020-12-01 | 2021-04-09 | 淮阴师范学院 | 基于k近邻回归算法的水稻单位面积地上部总有机碳含量遥感反演模型和方法 |
| CN112666091A (zh) * | 2020-12-01 | 2021-04-16 | 淮阴师范学院 | 基于随机梯度下降回归算法的水稻单位面积地上部可溶性糖积累量遥感反演模型和方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103048276A (zh) * | 2012-12-14 | 2013-04-17 | 北京农业信息技术研究中心 | 一种用于检测作物冠层叶片碳氮比的光谱指数构造方法 |
| CN104502283A (zh) * | 2014-12-15 | 2015-04-08 | 南京农业大学 | 估测大豆产量和地上部干重的一套两波段高光谱指数和预测模型 |
| CN106290197A (zh) * | 2016-09-06 | 2017-01-04 | 西北农林科技大学 | 水稻叶片全氮含量高光谱估测及估测模型构建方法 |
| CN107024439A (zh) * | 2017-03-23 | 2017-08-08 | 西北农林科技大学 | 一种水稻不同生育期叶绿素含量高光谱估测方法 |
| CN109754127A (zh) * | 2019-01-31 | 2019-05-14 | 浙江大学 | 基于无人机成像高光谱的水稻籽粒直链淀粉含量估测方法 |
| CN110376167A (zh) * | 2019-07-29 | 2019-10-25 | 北京麦飞科技有限公司 | 基于无人机高光谱的水稻叶片氮含量监测方法 |
-
2020
- 2020-08-04 CN CN202010772402.9A patent/CN111855592A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103048276A (zh) * | 2012-12-14 | 2013-04-17 | 北京农业信息技术研究中心 | 一种用于检测作物冠层叶片碳氮比的光谱指数构造方法 |
| CN104502283A (zh) * | 2014-12-15 | 2015-04-08 | 南京农业大学 | 估测大豆产量和地上部干重的一套两波段高光谱指数和预测模型 |
| CN106290197A (zh) * | 2016-09-06 | 2017-01-04 | 西北农林科技大学 | 水稻叶片全氮含量高光谱估测及估测模型构建方法 |
| CN107024439A (zh) * | 2017-03-23 | 2017-08-08 | 西北农林科技大学 | 一种水稻不同生育期叶绿素含量高光谱估测方法 |
| CN109754127A (zh) * | 2019-01-31 | 2019-05-14 | 浙江大学 | 基于无人机成像高光谱的水稻籽粒直链淀粉含量估测方法 |
| CN110376167A (zh) * | 2019-07-29 | 2019-10-25 | 北京麦飞科技有限公司 | 基于无人机高光谱的水稻叶片氮含量监测方法 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112270132A (zh) * | 2020-11-11 | 2021-01-26 | 仲恺农业工程学院 | 基于高斯过程回归算法的水稻单位面积内地上部干物重遥感反演模型和方法 |
| CN112362659A (zh) * | 2020-11-26 | 2021-02-12 | 淮阴师范学院 | 基于贝叶斯岭回归算法的水稻地上部糖氮比遥感反演模型和方法 |
| CN112630161A (zh) * | 2020-12-01 | 2021-04-09 | 淮阴师范学院 | 基于k近邻回归算法的水稻单位面积地上部总有机碳含量遥感反演模型和方法 |
| CN112666091A (zh) * | 2020-12-01 | 2021-04-16 | 淮阴师范学院 | 基于随机梯度下降回归算法的水稻单位面积地上部可溶性糖积累量遥感反演模型和方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111855589A (zh) | 水稻叶片氮积累量遥感反演模型和方法 | |
| CN111855590A (zh) | 水稻叶片淀粉积累量遥感反演模型和方法 | |
| CN111829965A (zh) | 水稻地上部淀粉积累量遥感反演模型和方法 | |
| CN111855591A (zh) | 水稻地上部碳氮比遥感反演模型和方法 | |
| CN111855592A (zh) | 水稻单位面积内地上部干物重遥感反演模型和方法 | |
| CN111855593A (zh) | 水稻叶片淀粉含量遥感反演模型和方法 | |
| CN113268923B (zh) | 一种基于模拟多光谱的夏玉米产量估算方法 | |
| CN109187441B (zh) | 基于冠层光谱信息的夏玉米含氮量监测模型的构建方法 | |
| CN107505271B (zh) | 基于氮素组分辐射传输模型的植株氮素估算方法和系统 | |
| CN102788752B (zh) | 基于光谱技术的农作物内部信息无损检测装置及方法 | |
| CN107796764B (zh) | 一种基于三波段植被指数的小麦叶面积指数估算模型的构建方法 | |
| CN105445229B (zh) | 基于氮平衡光谱法的小麦春季氮肥施用方法及其追氮量模型的构建方法 | |
| CN110567892B (zh) | 一种基于临界氮浓度的夏玉米氮素高光谱预测方法 | |
| CN110189793B (zh) | 基于高光谱的小麦氮肥生理利用率估测模型构建及不同氮效率小麦品种分类 | |
| CN112816618A (zh) | 一种氮高效小麦品种的筛选方法 | |
| CN112270131A (zh) | 基于ard回归算法的水稻叶面积指数遥感反演模型和方法 | |
| CN112362812A (zh) | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 | |
| CN109214591B (zh) | 一种木本植物地上生物量预测方法及系统 | |
| CN105608296B (zh) | 一种基于荔枝冠层光谱的叶片钾浓度反演方法 | |
| CN113049499A (zh) | 水体总氮浓度间接遥感反演方法、存储介质和终端设备 | |
| CN116385784A (zh) | 一种镉胁迫下水稻叶绿素含量测算方法及系统 | |
| CN112632847A (zh) | 基于XGBoost回归算法的水稻叶片淀粉含量遥感反演模型和方法 | |
| CN112525831A (zh) | 基于CatBoost回归算法的水稻地上部蛋白氮积累量遥感反演模型和方法 | |
| CN112630167A (zh) | 基于支持向量机回归算法的水稻叶片碳氮比遥感反演模型和方法 | |
| CN112270130A (zh) | 基于LightGBM回归算法的水稻叶片可溶性糖积累量遥感反演模型和方法 |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201030 |
|
| RJ01 | Rejection of invention patent application after publication |