CN112362812A - 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 - Google Patents
基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 Download PDFInfo
- Publication number
- CN112362812A CN112362812A CN202011249318.5A CN202011249318A CN112362812A CN 112362812 A CN112362812 A CN 112362812A CN 202011249318 A CN202011249318 A CN 202011249318A CN 112362812 A CN112362812 A CN 112362812A
- Authority
- CN
- China
- Prior art keywords
- chlorophyll
- rice
- carotenoid content
- content ratio
- model
- 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 200
- 235000009566 rice Nutrition 0.000 title claims abstract description 200
- 235000021466 carotenoid Nutrition 0.000 title claims abstract description 182
- 229930002875 chlorophyll Natural products 0.000 title claims abstract description 182
- 235000019804 chlorophyll Nutrition 0.000 title claims abstract description 182
- -1 chlorophyll carotenoid Chemical class 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 53
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 200
- 238000002310 reflectometry Methods 0.000 claims description 42
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 39
- 150000001747 carotenoids Chemical class 0.000 claims description 34
- 238000012549 training Methods 0.000 claims description 29
- 239000000523 sample Substances 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000005070 sampling Methods 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 19
- 238000002835 absorbance Methods 0.000 claims description 10
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 229930002868 chlorophyll a Natural products 0.000 claims description 5
- 229930002869 chlorophyll b Natural products 0.000 claims description 5
- NSMUHPMZFPKNMZ-VBYMZDBQSA-M chlorophyll b Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C=O)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 NSMUHPMZFPKNMZ-VBYMZDBQSA-M 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 4
- 238000000870 ultraviolet spectroscopy Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 abstract description 26
- 230000000694 effects Effects 0.000 abstract description 13
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 16
- 230000012010 growth Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 8
- 229910052757 nitrogen Inorganic materials 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 238000007689 inspection Methods 0.000 description 6
- 230000035764 nutrition Effects 0.000 description 6
- 235000016709 nutrition Nutrition 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 230000035790 physiological processes and functions Effects 0.000 description 4
- 239000000049 pigment Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000032683 aging Effects 0.000 description 3
- 239000003337 fertilizer Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000618 nitrogen fertilizer Substances 0.000 description 3
- 230000029553 photosynthesis Effects 0.000 description 3
- 238000010672 photosynthesis Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000000985 reflectance spectrum Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 240000004859 Gamochaeta purpurea 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
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000003851 biochemical process Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 210000003763 chloroplast Anatomy 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion 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
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003711 photoprotective effect Effects 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
- 238000011155 quantitative monitoring Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035882 stress Effects 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/0098—Plants or trees
-
- 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
- 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
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/33—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/01—Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Crystallography & Structural Chemistry (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Databases & Information Systems (AREA)
- Botany (AREA)
- Wood Science & Technology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的Lars回归模型,并进一步提供了该Lars回归模型的模型参数。还提供了基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法。本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演模型精度。
Description
技术领域
本发明涉及农业遥感技术领域,特别涉及水稻叶片叶绿素类胡萝卜素含量比值测量技术领域,具体是指一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法。
背景技术
水稻叶片中叶绿体是水稻进行光合作用的主要场所。光合作用中吸收光能的色素主要包括叶绿素以及类胡萝卜素,是水稻与外界进行物质交换和能量转换的物质基础(杨杰.基于叶片高光谱指数的水稻氮素及色素含量监测研究[D]. 南京农业大学,2009)。其中,叶绿素(包含叶绿素a和叶绿素b)是主要的吸收光能物质,直接影响了水稻光和作用的光能利用;类胡萝卜素具有吸收和传递光能以及光保护功能。叶绿素和类胡萝卜素的含量,尤其是叶绿素和类胡萝卜素含量比值反映了水稻营养胁迫,尤其是氮素营养状况,同时也反映了光合能力以及衰老进程等一系列重要的生理状态和生长情况(周贤锋.色素含量比值进行作物氮素营养状况诊断方法研究[D].2017)。
定量监测水稻叶片叶绿素类胡萝卜素含量比值,不仅可以了解水稻光合作用以及叶片衰老进程等重要的生理状态,同时能够反映水稻氮素状况,为动态管理水稻的氮肥施用,减少水稻生产中氮肥的使用量,缓解大量氮肥施用带来的环境问题,从而产生显著的经济和社会效益(陈青春,田永超,姚霞,等.基于冠层反射光谱的水稻追氮调控效应研究[J].中国农业科学,2010, 43(20):4149-4157)。传统的水稻叶片叶绿素类胡萝卜素含量比值的监测主要采用破坏性采样的方法,需要在室内进行测定,费时费力,且时效性差,无法及时的获取水稻叶片色素含量,不利于推广应用。
在水稻的生理生化过程中,水稻植株体内某些特定物质和细胞结构的变化,导致了水稻反射光谱的改变。因此,可以用光谱的变化来获取水稻叶片叶绿素类胡萝卜素含量比值等水稻生长信息(周冬琴.基于冠层反射光谱的水稻氮素营养与籽粒品质监测[D].南京农业大学,2007)。目前,使用高光谱来监测水稻的生长状态以及生理参数已经在作物生长监测中得到了应用。随着光谱技术的发展和普及,高光谱数据可以快捷迅速的获取水稻叶片叶绿素类胡萝卜素含量比值信息已经成为越来越多水稻生产从业者和研究人员的共识。最常用的方式是使用便携式全波段光谱仪获取水稻生长信息,选择能够反映叶片叶绿素类胡萝卜素含量比值的特征波段构建反演模型。在构建水稻叶片叶绿素类胡萝卜素含量比值反演模型的过程中,全波段光谱仪所测定的光谱范围涵盖了 350nm~2500nm,但是由于水稻组分复杂,组分光谱特征波段部分重叠,水稻叶片叶绿素类胡萝卜素含量比值特征光谱的确定困难,同时高光谱数据的快速处理成为基于高光谱数据估测水稻叶片叶绿素类胡萝卜素含量比值的亟待解决的技术问题。
因此,希望提供一种水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演模型精度。
发明内容
为了克服上述现有技术中的缺点,本发明的一个目的在于提供一种基于 Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演模型精度,适于大规模推广应用。
本发明的另一目的在于提供一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其设计巧妙,计算简便,易于实现,成本低,适于大规模推广应用。
本发明的另一目的在于提供一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演精度,适于大规模推广应用。
本发明的另一目的在于提供一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其设计巧妙,操作简单方便,成本低,适于大规模推广应用。
为达到以上目的,在本发明的第一方面,提供一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其特点是,所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的Lars回归模型,所述Lars回归模型的模型参数为:'eps'=0.13382952895927658, 'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True, 'jitter'=None。
较佳地,所述Lars回归模型采用水稻的数据集训练而成,所述数据集包括所述水稻的m个样点的冠层反射率和叶片叶绿素类胡萝卜素含量比值,m个所述样点均匀分布在水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率。
更佳地,所述m为36,所述n个特征波段为2151个特征波段,所述2151个特征波段为从350nm波段至2500nm波段。
在本发明的第二方面,提供了一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特点是,包括以下步骤:
(1)测量水稻的冠层反射率;
(2)测量所述水稻的叶片叶绿素类胡萝卜素含量比值:测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量,将所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量进行比值计算,获得所述叶片叶绿素类胡萝卜素含量比值;
(3)以所述冠层反射率为输入数据,采用Python语言的Lars回归模型进行计算,获得反演值,根据所述反演值与所述叶片叶绿素类胡萝卜素含量比值计算决定系数R2,改变所述Lars回归模型的模型参数的取值,R2的变化越大,说明所述模型参数的重要性越大,将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵;
(4)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为输出结果,训练所述Lars回归模型,根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优,获得所述模型参数的调优值;
(5)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为所述输出结果,采用所述的模型参数的调优值,训练所述Lars回归模型,待所述Lars回归模型训练结束后,获得基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用save方法保存所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,如果需要使用所述的基于 Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用load方法加载所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型使用。
较佳地,在所述步骤(1)中,所述测量采用高光谱辐射仪进行,所述测量的时间为10:00~14:00,所述高光谱辐射仪采用视场角为25°的镜头,所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米,所述传感器探头的地面视场范围直径为0.44米,所述传感器探头迎向阳光,所述测量采用标准板进行校正,所述标准板是反射率为95%~ 99%的标准白板。
较佳地,在所述步骤(2)中,所述的测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量的步骤具体包括:
采集所述的水稻的叶片,加入叶绿素类胡萝卜素提取液于暗处提取,测定所述叶绿素类胡萝卜素提取液在470nm、645nm和663nm处的吸光度A470、A645和A663,使用下列公式计算叶绿素a浓度Chla、叶绿素b浓度Chlb、叶绿素浓度Chl 和类胡萝卜素浓度Car:
Chla(mg/L)=12.21×A663-2.81×A645,
Chlb(mg/L)=20.13×A645-5.03×A663,
Chl=Chla+Chlb,
Car(mg/L)=(1000×A470–3.27×Chla–104×Chlb)/299,
然后根据取样分布面积把所述叶绿素浓度Chl和所述类胡萝卜素浓度Car转换为所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量。
更佳地,在所述步骤(2)中,所述叶绿素类胡萝卜素提取液是无水乙醇和丙酮1:1混合液,所述吸光度采用紫外可见分光光度计测定。
较佳地,在所述步骤(3)中,所述模型参数调优秩次矩阵为:
Params={'eps','n_nonzero_coefs','normalize','precompute','fit_intercept', 'jitter'}。
更佳地,在所述步骤(4)中,所述的模型参数的调优值为:
'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True,'jitter'=None。
较佳地,在所述步骤(1)中,所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率,m个所述样点均匀分布在所述水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率;在所述步骤(2) 中,所述的测量所述水稻的叶片叶绿素类胡萝卜素含量比值的步骤具体为测量该 m个所述样点的所述叶片叶绿素类胡萝卜素含量比值。
更佳地,在所述步骤(1)中,所述m为36,所述n个特征波段为2151 个特征波段,所述2151个特征波段为从350nm波段至2500nm波段。
本发明的有益效果主要在于:
1、本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的Lars回归模型,Lars回归模型的模型参数为: 'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False, 'precompute'='auto','fit_intercept'=True,'jitter'=None,对该模型进行检验,R2在 0.85以上,因此,其能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演模型精度,适于大规模推广应用。
2、本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的Lars回归模型,Lars回归模型的模型参数为: 'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False, 'precompute'='auto','fit_intercept'=True,'jitter'=None,对该模型进行检验,R2在 0.85以上,因此,其设计巧妙,计算简便,易于实现,成本低,适于大规模推广应用。
3、本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,包括:测量水稻的冠层反射率;测量水稻的叶片叶绿素类胡萝卜素含量比值:以冠层反射率为输入数据,采用Python语言的Lars回归模型进行计算,以决定系数R2构建模型参数调优秩次矩阵;以冠层反射率为输入数据,以叶片叶绿素类胡萝卜素含量比值为输出结果,训练Lars回归模型,根据模型参数调优秩次矩阵依次对模型参数进行调优,获得模型参数的调优值;以冠层反射率为输入数据,以叶片叶绿素类胡萝卜素含量比值为输出结果,采用模型参数的调优值,训练Lars回归模型,获得基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,对该模型进行检验,R2在0.85以上,因此,其能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演精度,适于大规模推广应用。
4、本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,包括:测量水稻的冠层反射率;测量水稻的叶片叶绿素类胡萝卜素含量比值:以冠层反射率为输入数据,采用Python语言的Lars回归模型进行计算,以决定系数R2构建模型参数调优秩次矩阵;以冠层反射率为输入数据,以叶片叶绿素类胡萝卜素含量比值为输出结果,训练Lars回归模型,根据模型参数调优秩次矩阵依次对模型参数进行调优,获得模型参数的调优值;以冠层反射率为输入数据,以叶片叶绿素类胡萝卜素含量比值为输出结果,采用模型参数的调优值,训练Lars回归模型,获得基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,对该模型进行检验,R2在0.85以上,因此,其设计巧妙,操作简单方便,成本低,适于大规模推广应用。
本发明的这些和其它目的、特点和优势,通过下述的详细说明、附图和权利要求得以充分体现,并可通过所附权利要求中特地指出的手段、装置和它们的组合得以实现。
附图说明
图1是本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法的一具体实施例的流程示意图。
图2是图1所示的具体实施例的模型构建流程示意图。
图3是图1所示的具体实施例的模型检验结果示意图,其中实测值和预测值均为叶片叶绿素类胡萝卜素含量比值,无量纲。
具体实施方式
本发明人针对基于高光谱估测水稻叶片叶绿素类胡萝卜素含量比值的需求,克服了因水稻组分复杂而导致的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定以及高光谱数据特征波段筛选费时费力的困难,提出了一种基于 Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,所述的基于 Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的 Lars回归模型,所述Lars回归模型的模型参数为:'eps'=0.13382952895927658, 'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True, 'jitter'=None。
所述Lars回归模型可以采用任何合适的数据集训练而成,较佳地,所述Lars 回归模型采用水稻的数据集训练而成,所述数据集包括所述水稻的m个样点的冠层反射率和叶片叶绿素类胡萝卜素含量比值,m个所述样点均匀分布在水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率。所述水稻种植区域可以是多个生态点、多个品种的水稻种植区域。
所述m和所述n为正整数,可以根据需要确定,更佳地,所述m为36,所述n 个特征波段为2151个特征波段,所述2151个特征波段为从350nm波段至2500nm 波段。
本发明还提供了一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,包括以下步骤:
(1)测量水稻的冠层反射率;
(2)测量所述水稻的叶片叶绿素类胡萝卜素含量比值:测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量,将所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量进行比值计算,获得所述叶片叶绿素类胡萝卜素含量比值;
(3)以所述冠层反射率为输入数据,采用Python语言的Lars回归模型进行计算,获得反演值,根据所述反演值与所述叶片叶绿素类胡萝卜素含量比值计算决定系数R2,改变所述Lars回归模型的模型参数的取值,R2的变化越大,说明所述模型参数的重要性越大,将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵;
(4)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为输出结果,训练所述Lars回归模型,根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优,获得所述模型参数的调优值;
(5)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为所述输出结果,采用所述的模型参数的调优值,训练所述Lars回归模型,待所述Lars回归模型训练结束后,获得基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用save方法保存所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,如果需要使用所述的基于 Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用load方法加载所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型使用。
在所述步骤(1)中,所述测量可以采用任何合适的光谱仪和方法进行,较佳地,在所述步骤(1)中,所述测量采用高光谱辐射仪进行,所述测量的时间为10:00~14:00,所述高光谱辐射仪采用视场角为25°的镜头,所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米,所述传感器探头的地面视场范围直径为0.44米,所述传感器探头迎向阳光,所述测量采用标准板进行校正,所述标准板是反射率为95%~99%的标准白板。
在所述步骤(2)中,所述的测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量的步骤可以具体包括任何合适的方法,较佳地,
在所述步骤(2)中,所述的测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量的步骤具体包括:
采集所述的水稻的叶片,加入叶绿素类胡萝卜素提取液于暗处提取,测定所述叶绿素类胡萝卜素提取液在470nm、645nm和663nm处的吸光度A470、A645和A663,使用下列公式计算叶绿素a浓度Chla、叶绿素b浓度Chlb、叶绿素浓度Chl 和类胡萝卜素浓度Car:
Chla(mg/L)=12.21×A663-2.81×A645,
Chlb(mg/L)=20.13×A645-5.03×A663,
Chl=Chla+Chlb,
Car(mg/L)=(1000×A470–3.27×Chla–104×Chlb)/299,
然后根据取样分布面积把所述叶绿素浓度Chl和所述类胡萝卜素浓度Car转换为所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量。
在所述步骤(2)中,所述叶绿素类胡萝卜素提取液可以是任何合适的叶绿素类胡萝卜素提取液,所述吸光度可以采用任何合适的仪器进行测定,更佳地,在所述步骤(2)中,所述叶绿素类胡萝卜素提取液是无水乙醇和丙酮1:1混合液,所述吸光度采用紫外可见分光光度计测定。
在所述步骤(3)中,所述模型参数调优秩次矩阵根据决定系数R2确定,较佳地,在所述步骤(3)中,所述模型参数调优秩次矩阵为:
Params={'eps','n_nonzero_coefs','normalize','precompute','fit_intercept', 'jitter'}。
在所述步骤(4)中,所述的模型参数的调优值根据所述模型参数调优秩次矩阵依次确定,更佳地,在所述步骤(4)中,所述的模型参数的调优值为:
'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True,'jitter'=None。
为了使得所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型的精度更高,可以选取水稻种植区域的多个样点,测定多个样点的多个特征波段的冠层反射率以及多个样点的叶片叶绿素类胡萝卜素含量比值,较佳地,在所述步骤(1)中,所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率,m个所述样点均匀分布在所述水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率;在所述步骤(2)中,所述的测量所述水稻的叶片叶绿素类胡萝卜素含量比值的步骤具体为测量该m 个所述样点的所述叶片叶绿素类胡萝卜素含量比值。
在所述步骤(1)中,所述m和所述n为正整数,可以根据需要确定,更佳地,在所述步骤(1)中,所述m为36,所述n个特征波段为2151个特征波段,所述2151个特征波段为从350nm波段至2500nm波段。
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。
实施例
本实施例的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法基于实测高光谱数据,采用了水稻种植区域(江苏省淮安市农业科学研究院淮安地区稻麦种植基地,水稻品种为淮稻5号,采样时期为水稻拔节期)采集的水稻冠层反射率光谱数据以及水稻叶片叶绿素类胡萝卜素含量比值数据,共 48个采样点,这些采样点均匀分布并完全覆盖水稻种植区域整个区域。48个采样点数据采用随机的方法分为两部分,其中36个采样点的数据用于模型构建,12个采样点数据用于模型检验。基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法的流程如图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 株,用吸水纸包好,带回实验室,分离叶片,取约0.1克叶片鲜样,用无水乙醇和丙酮1:1混合液于暗处提取,使用日立(HITACHI)U-2800型紫外可见分光光度计测定470nm处的吸光度A470、645nm处的吸光度A645和663nm处的吸光度A663,然后使用下列公式计算叶片的叶绿素a浓度Chla、叶绿素b浓度Chlb、叶绿素浓度Chl和类胡萝卜素浓度Car:
Chla(mg/L)=12.21×A663-2.81×A645
Chlb(mg/L)=20.13×A645-5.03×A663
Chl=Chla+Chlb
Car(mg/L)=(1000×A470–3.27×Chla–104×Chlb)/299
最后,根据取样分布面积把叶绿素浓度和类胡萝卜素浓度转换为单位土地面积叶绿素含量Chlsoil和单位土地面积类胡萝卜素含量Carsoil,计算单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量的比值(CCR),即为叶片叶绿素类胡萝卜素含量比值:
CCR=Chlsoil/Carsoil
3.模型构建
模型构建采用Python语言的Lars回归模型进行构建,请参见图2所示,模型构建主要包括以下步骤:
3.1数据校验
对获取的水稻冠层反射率数据进行校验,剔除异常的整条光谱曲线数据。本发明中异常光谱是指相邻光谱变化超过100%,且包含空值、负值的光谱值。
3.2数据的预处理
对校验后的水稻冠层反射率数据以及水稻叶片叶绿素类胡萝卜素含量比值数据进行预处理,包括去除包含缺失值和空值的成对水稻冠层反射率数据和水稻叶片叶绿素类胡萝卜素含量比值数据。
3.3数据集的划分
为保证模型训练和反演结果的合理评估,使用随机的方法,将整个数据集(48 组)划分为两部分,75%(36组)的数据用于模型训练,25%(12组)的数据用于训练后的效果评估。
3.4训练数据集的划分
为保证模型训练的效果,使用随机的方法,在每次模型训练迭代时,将训练数据集划分为5个部分,对模型进行训练。
3.5模型参数调优秩次矩阵的构建
本发明中,模型训练过程模型参数的调优非常重要,为保证尽量获取最佳的模型调优,使用试错法进行模型参数调优。本发明使用决定系数R2(R2越接近于1 越好)作为检验参数,构建了评估模型参数权重的参数秩次矩阵。根据训练数据集首先使用模型参数的默认值进行计算,获取反演值,根据反演值和叶片叶绿素类胡萝卜素含量比值,计算决定系数R2,然后改变模型参数的取值,R2的变化越大,说明模型参数的重要性越大,将模型参数按照重要性的大小由大到小进行排列构建模型参数调优秩次矩阵,用于后续计算。
根据训练数据集中冠层反射率数据和对应的叶片叶绿素类胡萝卜素含量比值数据,计算得到的模型参数调优秩次矩阵为:
Params={'eps','n_nonzero_coefs','normalize','precompute','fit_intercept', 'jitter'}。
3.6模型构建
依据获取的模型参数调优秩次矩阵,将建模使用的数据,包括实测的冠层反射率数据和对应的实测的叶片叶绿素类胡萝卜素含量比值数据,以实测的冠层反射率数据为输入数据,实测的叶片叶绿素类胡萝卜素含量比值数据为输出结果,训练Lars回归模型,根据模型参数调优秩次矩阵依次对模型参数进行调优,获得该模型完整参数及取值,如下:
'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True,'jitter'=None。
待模型训练结束后,使用save方法保存模型,如需使用,运行load方法加载使用。
3.7模型检验
使用构建模型之外的12个采样点高光谱数据输入模型,使用调优后的模型参数进行计算,得到预测值,分析预测值与实测值(叶片叶绿素类胡萝卜素含量比值) 的关系,结果如图3所示,模型的R2为0.8984。而使用默认参数的模型R2为0。
本案例训练数据和检验数据的随机划分以及模型的构建、训练和检验使用了美国MathWorks公司开发的Matlab软件(version:R2020a 9.8.0.1380330)以及Python(version:3.7.0),通过Matlab软件调用Python的Lars回归模型。
因此,本发明提出了一种新的、基于实测高光谱遥感数据的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,基于实测水稻冠层反射率数据和实地采集的水稻叶片叶绿素类胡萝卜素含量比值数据,可快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,并通过构建模型参数调优秩次矩阵,使用试错法对模型参数调优,有效降低了线性模型过拟合的现象,大大提高了水稻叶片叶绿素类胡萝卜素含量比值反演的精度,适用于不同生态区、不同品种以及主要生育期的水稻叶片叶绿素类胡萝卜素含量比值的定量反演,从而获得水稻氮素营养、生理状态以及水肥供给的状态,提高了水稻栽培种植过程中的生长信息获取效率,并为水稻生产中水分肥料运筹提供了基础科学数据。
与现有技术相比,本发明具有以下优点:
(1)本发明使用的Lars回归模型(Least Angle Regression Model,LARM,最小角回归模型)适用于基于高光谱的水稻叶片叶绿素类胡萝卜素含量比值的反演,在综合考虑高光谱350~2500nm波段范围信息的基础上,考虑了水稻体内各种物质组成和细胞结构的光学特性,特别是复杂的组分对水稻叶片叶绿素类胡萝卜素含量比值特征波段的影响和叠加效应,充分利用遥感数据中不同波段包含的水稻叶片叶绿素类胡萝卜素含量比值信息,进行水稻叶片叶绿素类胡萝卜素含量比值的反演;
(2)使用了最小角回归机器学习算法,构建350~2500nm反射率与水稻叶片叶绿素类胡萝卜素含量比值模型,通过构建不同的分段线性的函数求解路径,提高了模型构建以及训练的效率,有效的利用了高光谱分辨率的地面反射率所包含叶绿素类胡萝卜素含量比值信息,也提高基于高光谱信息的水稻叶片叶绿素类胡萝卜素含量比值的反演速度和效率;
(3)本发明充分考虑了模型训练和模型检验的独立性,使用随机分割的方法划分训练数据集和检验数据集,训练数据集只用于模型训练,检验数据集仅用于模型检验,以保证模型效果检验的合理性。
(4)由于模型的参数调优对于模型的计算精度非常重要,因此本发明中构建了模型参数秩次矩阵,以决定系数R2为评价参数,使用试错法进行模型参数调优,在保证参数调优效果的基础上,大大提高了模型训练和参数调优的速度。
(5)本发明提出的水稻叶片叶绿素类胡萝卜素含量比值反演方法计算简便,适用于不同生态区、不同品种以及不同生育期的水稻叶片叶绿素类胡萝卜素含量比值的遥感定量反演,可准确的反演水稻叶片叶绿素类胡萝卜素含量比值,快速的获取水稻氮素营养、生理状况和长势等信息,同时为水稻种植栽培的水分肥料运筹管理提供科学数据。
综上,本发明的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型能够快速、准确的获取水稻叶片叶绿素类胡萝卜素含量比值信息,克服水稻组分复杂带来的光谱叠加效应造成的水稻叶片叶绿素类胡萝卜素含量比值特征波段难以确定的困难,大大提高水稻叶片叶绿素类胡萝卜素含量比值反演模型精度,设计巧妙,计算简便,易于实现,成本低,适于大规模推广应用。
由此可见,本发明的目的已经完整并有效的予以实现。本发明的功能及结构原理已在实施例中予以展示和说明,在不背离所述原理下,实施方式可作任意修改。所以,本发明包括了基于权利要求精神及权利要求范围的所有变形实施方式。
Claims (11)
1.一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其特征在于,所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型为Python语言的Lars回归模型,所述Lars回归模型的模型参数为:'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True,'jitter'=None。
2.如权利要求1所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其特征在于,所述Lars回归模型采用水稻的数据集训练而成,所述数据集包括所述水稻的m个样点的冠层反射率和叶片叶绿素类胡萝卜素含量比值,m个所述样点均匀分布在水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率。
3.如权利要求2所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,其特征在于,所述m为36,所述n个特征波段为2151个特征波段,所述2151个特征波段为从350nm波段至2500nm波段。
4.一种基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,包括以下步骤:
(1)测量水稻的冠层反射率;
(2)测量所述水稻的叶片叶绿素类胡萝卜素含量比值:测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量,将所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量进行比值计算,获得所述叶片叶绿素类胡萝卜素含量比值;
(3)以所述冠层反射率为输入数据,采用Python语言的Lars回归模型进行计算,获得反演值,根据所述反演值与所述叶片叶绿素类胡萝卜素含量比值计算决定系数R2,改变所述Lars回归模型的模型参数的取值,R2的变化越大,说明所述模型参数的重要性越大,将所述模型参数按照所述重要性的大小由大到小进行排列构建模型参数调优秩次矩阵;
(4)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为输出结果,训练所述Lars回归模型,根据所述模型参数调优秩次矩阵依次对所述模型参数进行调优,获得所述模型参数的调优值;
(5)以所述冠层反射率为所述输入数据,以所述叶片叶绿素类胡萝卜素含量比值为所述输出结果,采用所述的模型参数的调优值,训练所述Lars回归模型,待所述Lars回归模型训练结束后,获得基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用save方法保存所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,如果需要使用所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型,使用load方法加载所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型使用。
5.如权利要求4所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(1)中,所述测量采用高光谱辐射仪进行,所述测量的时间为10:00~14:00,所述高光谱辐射仪采用视场角为25°的镜头,所述便携式野外高光谱辐射仪的传感器探头垂直指向所述水稻的冠层并距离所述冠层的顶层的垂直高度为1米,所述传感器探头的地面视场范围直径为0.44米,所述传感器探头迎向阳光,所述测量采用标准板进行校正,所述标准板是反射率为95%~99%的标准白板。
6.如权利要求4所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(2)中,所述的测量所述水稻的叶片的单位土地面积叶绿素含量和单位土地面积类胡萝卜素含量的步骤具体包括:
采集所述的水稻的叶片,加入叶绿素类胡萝卜素提取液于暗处提取,测定所述叶绿素类胡萝卜素提取液在470nm、645nm和663nm处的吸光度A470、A645和A663,使用下列公式计算叶绿素a浓度Chla、叶绿素b浓度Chlb、叶绿素浓度Chl和类胡萝卜素浓度Car:
Chla(mg/L)=12.21×A663-2.81×A645,
Chlb(mg/L)=20.13×A645-5.03×A663,
Chl=Chla+Chlb,
Car(mg/L)=(1000×A470–3.27×Chla–104×Chlb)/299,
然后根据取样分布面积把所述叶绿素浓度Chl和所述类胡萝卜素浓度Car转换为所述单位土地面积叶绿素含量和所述单位土地面积类胡萝卜素含量。
7.如权利要求6所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(2)中,所述叶绿素类胡萝卜素提取液是无水乙醇和丙酮1:1混合液,所述吸光度采用紫外可见分光光度计测定。
8.如权利要求4所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(3)中,所述模型参数调优秩次矩阵为:
Params={'eps','n_nonzero_coefs','normalize','precompute','fit_intercept','jitter'}。
9.如权利要求8所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(4)中,所述的模型参数的调优值为:
'eps'=0.13382952895927658,'n_nonzero_coefs'=126,'normalize'=False,'precompute'='auto','fit_intercept'=True,'jitter'=None。
10.如权利要求4所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(1)中,所述的测量水稻的冠层反射率的步骤具体为测量水稻种植区域的m个样点的所述冠层反射率,m个所述样点均匀分布在所述水稻种植区域,所述冠层反射率为n个特征波段的冠层反射率;在所述步骤(2)中,所述的测量所述水稻的叶片叶绿素类胡萝卜素含量比值的步骤具体为测量该m个所述样点的所述叶片叶绿素类胡萝卜素含量比值。
11.如权利要求10所述的基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演方法,其特征在于,在所述步骤(1)中,所述m为36,所述n个特征波段为2151个特征波段,所述2151个特征波段为从350nm波段至2500nm波段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011249318.5A CN112362812A (zh) | 2020-11-10 | 2020-11-10 | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011249318.5A CN112362812A (zh) | 2020-11-10 | 2020-11-10 | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112362812A true CN112362812A (zh) | 2021-02-12 |
Family
ID=74509535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011249318.5A Pending CN112362812A (zh) | 2020-11-10 | 2020-11-10 | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112362812A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945524A (zh) * | 2021-09-30 | 2022-01-18 | 中国地质大学(武汉) | 基于植被指数的类胡萝卜素与叶绿素比值反演方法及系统 |
CN117664955A (zh) * | 2024-02-01 | 2024-03-08 | 中国林业科学研究院热带林业研究所 | 一种基于苗期叶表型快速鉴别奇楠沉香和普通沉香的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102519886A (zh) * | 2011-12-19 | 2012-06-27 | 浙江大学 | 一种检测农作物叶片中叶绿素a和类胡萝卜素含量的方法 |
CN107024439A (zh) * | 2017-03-23 | 2017-08-08 | 西北农林科技大学 | 一种水稻不同生育期叶绿素含量高光谱估测方法 |
CN109030378A (zh) * | 2018-06-04 | 2018-12-18 | 沈阳农业大学 | 基于pso-elm的粳稻冠层叶绿素含量反演模型方法 |
CN111027523A (zh) * | 2019-12-31 | 2020-04-17 | 塔里木大学 | 一种棉花冠层类胡萝卜素含量的卫星遥感监测方法 |
CN111044516A (zh) * | 2019-12-26 | 2020-04-21 | 沈阳农业大学 | 一种水稻叶绿素含量遥感估测方法 |
CN111781147A (zh) * | 2020-07-06 | 2020-10-16 | 淮阴师范学院 | 基于两波段的内陆湖泊水体叶绿素a浓度遥感反演模型和方法 |
CN111855590A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片淀粉积累量遥感反演模型和方法 |
CN111855589A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片氮积累量遥感反演模型和方法 |
CN111855593A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片淀粉含量遥感反演模型和方法 |
-
2020
- 2020-11-10 CN CN202011249318.5A patent/CN112362812A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102519886A (zh) * | 2011-12-19 | 2012-06-27 | 浙江大学 | 一种检测农作物叶片中叶绿素a和类胡萝卜素含量的方法 |
CN107024439A (zh) * | 2017-03-23 | 2017-08-08 | 西北农林科技大学 | 一种水稻不同生育期叶绿素含量高光谱估测方法 |
CN109030378A (zh) * | 2018-06-04 | 2018-12-18 | 沈阳农业大学 | 基于pso-elm的粳稻冠层叶绿素含量反演模型方法 |
CN111044516A (zh) * | 2019-12-26 | 2020-04-21 | 沈阳农业大学 | 一种水稻叶绿素含量遥感估测方法 |
CN111027523A (zh) * | 2019-12-31 | 2020-04-17 | 塔里木大学 | 一种棉花冠层类胡萝卜素含量的卫星遥感监测方法 |
CN111781147A (zh) * | 2020-07-06 | 2020-10-16 | 淮阴师范学院 | 基于两波段的内陆湖泊水体叶绿素a浓度遥感反演模型和方法 |
CN111855590A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片淀粉积累量遥感反演模型和方法 |
CN111855589A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片氮积累量遥感反演模型和方法 |
CN111855593A (zh) * | 2020-08-04 | 2020-10-30 | 淮阴师范学院 | 水稻叶片淀粉含量遥感反演模型和方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945524A (zh) * | 2021-09-30 | 2022-01-18 | 中国地质大学(武汉) | 基于植被指数的类胡萝卜素与叶绿素比值反演方法及系统 |
CN117664955A (zh) * | 2024-02-01 | 2024-03-08 | 中国林业科学研究院热带林业研究所 | 一种基于苗期叶表型快速鉴别奇楠沉香和普通沉香的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102788752B (zh) | 基于光谱技术的农作物内部信息无损检测装置及方法 | |
CN111855589A (zh) | 水稻叶片氮积累量遥感反演模型和方法 | |
CN111855591A (zh) | 水稻地上部碳氮比遥感反演模型和方法 | |
CN107796764B (zh) | 一种基于三波段植被指数的小麦叶面积指数估算模型的构建方法 | |
CN111855590A (zh) | 水稻叶片淀粉积累量遥感反演模型和方法 | |
CN110567892B (zh) | 一种基于临界氮浓度的夏玉米氮素高光谱预测方法 | |
CN107505271B (zh) | 基于氮素组分辐射传输模型的植株氮素估算方法和系统 | |
CN111829965A (zh) | 水稻地上部淀粉积累量遥感反演模型和方法 | |
CN113268923B (zh) | 一种基于模拟多光谱的夏玉米产量估算方法 | |
CN109187398A (zh) | 一种小麦植株氮含量的高光谱测定方法 | |
CN103868880A (zh) | 基于光谱双峰指数的小麦叶片氮含量监测方法及其监测模型的构建方法 | |
CN111855593A (zh) | 水稻叶片淀粉含量遥感反演模型和方法 | |
CN112270131A (zh) | 基于ard回归算法的水稻叶面积指数遥感反演模型和方法 | |
CN110082309B (zh) | 冬小麦冠层spad值综合光谱监测模型建立方法 | |
CN104778349B (zh) | 一种用于水稻表土氮肥施用等级评定方法 | |
CN111044516B (zh) | 一种水稻叶绿素含量遥感估测方法 | |
CN111855592A (zh) | 水稻单位面积内地上部干物重遥感反演模型和方法 | |
CN112362812A (zh) | 基于Lars算法的水稻叶片叶绿素类胡萝卜素含量比值遥感反演模型和方法 | |
CN112816618A (zh) | 一种氮高效小麦品种的筛选方法 | |
CN114460015A (zh) | 一种用于预测寒地水稻冠层叶绿素含量的方法及模型 | |
CN113049499A (zh) | 水体总氮浓度间接遥感反演方法、存储介质和终端设备 | |
CN112525831A (zh) | 基于CatBoost回归算法的水稻地上部蛋白氮积累量遥感反演模型和方法 | |
CN112630167A (zh) | 基于支持向量机回归算法的水稻叶片碳氮比遥感反演模型和方法 | |
CN112632847A (zh) | 基于XGBoost回归算法的水稻叶片淀粉含量遥感反演模型和方法 | |
CN112270133A (zh) | 基于Huber回归算法的水稻叶片叶绿素含量比值遥感反演模型和方法 |
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 |