CN110874617B - 一种冬小麦叶片氮含量估算模型的建立方法 - Google Patents
一种冬小麦叶片氮含量估算模型的建立方法 Download PDFInfo
- Publication number
- CN110874617B CN110874617B CN201911169596.7A CN201911169596A CN110874617B CN 110874617 B CN110874617 B CN 110874617B CN 201911169596 A CN201911169596 A CN 201911169596A CN 110874617 B CN110874617 B CN 110874617B
- Authority
- CN
- China
- Prior art keywords
- wheat
- model
- nitrogen content
- color
- color space
- 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.)
- Active
Links
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 177
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 88
- 241000209140 Triticum Species 0.000 title claims abstract description 79
- 235000021307 Triticum Nutrition 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000001419 dependent effect Effects 0.000 claims abstract description 9
- 238000003064 k means clustering Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 14
- 238000013507 mapping Methods 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 239000004220 glutamic acid Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims 1
- 235000016709 nutrition Nutrition 0.000 abstract description 16
- 230000035764 nutrition Effects 0.000 abstract description 14
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000007637 random forest analysis Methods 0.000 description 42
- 238000012549 training Methods 0.000 description 29
- 238000004422 calculation algorithm Methods 0.000 description 28
- 238000012360 testing method Methods 0.000 description 21
- 230000006870 function Effects 0.000 description 15
- 241000196324 Embryophyta Species 0.000 description 12
- 238000002790 cross-validation Methods 0.000 description 12
- 238000011156 evaluation Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 229930002875 chlorophyll Natural products 0.000 description 8
- 235000019804 chlorophyll Nutrition 0.000 description 8
- 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 description 8
- 238000012417 linear regression Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000010276 construction Methods 0.000 description 5
- 230000004720 fertilization Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000003066 decision tree Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 240000007594 Oryza sativa Species 0.000 description 3
- 235000007164 Oryza sativa Nutrition 0.000 description 3
- 240000008042 Zea mays Species 0.000 description 3
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 238000007667 floating Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000618 nitrogen fertilizer Substances 0.000 description 3
- 235000015097 nutrients Nutrition 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 235000009566 rice Nutrition 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000004611 spectroscopical analysis Methods 0.000 description 3
- 229920000742 Cotton Polymers 0.000 description 2
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 2
- 230000009418 agronomic effect Effects 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 235000005822 corn Nutrition 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 240000006162 Chenopodium quinoa Species 0.000 description 1
- 244000300264 Spinacia oleracea Species 0.000 description 1
- 235000009337 Spinacia oleracea Nutrition 0.000 description 1
- 235000016383 Zea mays subsp huehuetenangensis Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000003608 fece Anatomy 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000010871 livestock manure Substances 0.000 description 1
- 235000009973 maize Nutrition 0.000 description 1
- 201000009240 nasopharyngitis Diseases 0.000 description 1
- 230000004962 physiological condition Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- 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)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种冬小麦叶片氮含量估算模型的建立方法,包括:采集小麦冠层图像,测定小麦氮素营养参数;利用冠层图像H分量的K均值聚类方法分割小麦冠层图像,提取分割处理后冠层图像的非0像素值的R、G、B、H、S、V、L*、a*、b*分量,分别计算其平均像素值作为基础颜色分量,构建融合三个常用颜色空间的多颜色空间;以图像在多颜色空间下的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,选用径向基核函数,建立ε‑SVR模型。本发明建立的多颜色空间下的支持向量回归模型有着更好的泛化性能,具有较高的预测精度,可对小麦叶片氮含量进行快速无损检测。
Description
技术领域
本发明涉及一种冬小麦叶片氮含量估算模型的建立方法。
背景技术
小麦是重要的粮食作物,在全球范围内广泛分布。我国作为小麦的主产国之一,小麦的生产效率就显得尤为重要。提高小麦的生产效率,可以解决我国的粮食储备安全问题,进而实现我国的粮食战略目标。氮素是作物生长的重要养分,影响作物产量,但过量施氮会污染农田生态坏境。对作物及时准确做出营养诊断,并根据诊断结果,制定施肥管理措施,有助于提高肥料利用率,改善农田生态环境,对农业的可持续发展有着重要的意义。
目前,作物氮营养无损诊断方法为多光谱及遥感方法。基于作物冠层的多光谱反射率的航空以及卫星遥感技术适用于监测大尺度下的作物长势和氮素营养状态,且受时空分辨率的限制,导致操作难度大,成本高。基于光谱学原理的地面遥感技术中,便携式叶绿素仪(Soil and Plant Analyzer Development,SPAD)虽然快速无损测定叶绿素相对含量,但对作物的氮饱和状态不敏感,且测定结果受到品种、生长期等因素较大的影响。植物冠层光谱仪GreenSeeker、ASD等获取的冠层光谱数据是包含背景因素的混合光谱,而背景因素影响作物冠层光谱的准确获取。作物的叶片会根据自身不同的营养状态,呈现差异性的颜色和形态,以此反映自身生理状况。有研究表明,基于连续投影算法提取光谱数据的波长位置,大部分位于可见光区域,光谱数据对叶片氮含量(Leaf Nnitrogen Concentration,LNC)最有研究意义的前两个波段均在可见光区域。由此可见,可见光区间对LNC较为敏感。
相关研究成果显示,不同作物,与之营养参数相关的图像参数各不相同。李红军[1]和白金顺[2]的研究发现,与拔节期小麦和灌浆期玉米氮素营养相关性较好的是作物冠层的红光标准值(Normalized Redness Intensity,NRI)。Karcher等试验发现,HSV空间下深绿指数(Dark Ggreen Color Index,DGCI)与草地氮素水平、玉米叶片含氮量存在线性相关性。Wang[3]等分析发现HSV颜色空间的H分量和颜色空间的分量与水稻叶片氮含量存在显著相关性。王方永[4]等研究表明HSI空间的S分量以及颜色空间的分量分别与棉花的氮素浓度、SPAD值的相关关系明显。Rigon[5]分析了RGB、HSB和空间下各颜色分量估算SPAD的线性模型,认为NRI、S以及分量与SPAD值线性关系最强。
可见,作物生长状况受品种、生育期等多个因素的影响,表征作物营养信息的图像参数随之变化。而任何一种颜色分量的变化都会引起作物冠层叶色不同程度的改变,冠层叶色又与氮含量密切相关。上述研究或利用图像单一颜色分量与作物营养指标拟合线性回归方程,未考虑剩余颜色分量对营养检测的影响;或研究单一颜色空间下构建的颜色指数与营养指标的关系,未对比其他颜色空间下构建的营养评价指标。Su[6]和Riccardi[7]综合考虑RGB空间下的各颜色分量,基于最小二乘法建立了估算微藻和藜麦叶绿素含量的多元线性回归(Multivariate Linear Regression,MLR)模型,均取得了较好的效果。在此基础上,Agarwal研究了多个颜色空间下估算菠菜SPAD的线性模型,指出空间下的MLR算法模型优于其他颜色空间下的。为提高估算精度,梁亮等比较了线性回归方法和支持向量回归(Support Vector Regression,SVR)算法遥感反演小麦冠层叶绿素的效果,发现支持向量回归显著提高了反演精度。王丽爱等比较了遥感估算小麦叶片SPAD值的3种机器学习算法模型效果,提出较支持向量回归和人工神经网络,随机森林(Random Forest,RF)算法的学习和回归预测能力在小麦多个生育期估算效果最优。
大量研究基于高光谱,利用机器学习算法进行农学参数估算,缺少对不同算法泛化性能的比较。另外,基于冠层图像可见光颜色分量来估算农学参数的研究较少。
参考文献:
[1]李红军,张立周,陈曦鸣,张玉铭,程一松,胡春胜.应用数字图像进行小麦氮素营养诊断中图像分析方法的研究[J].中国生态农业学报,2011,19(01):155-159.
[2]白金顺,曹卫东,熊静,曾闹华,志水胜好,芮玉奎.应用数码相机进行绿肥翻压后春玉米氮素营养诊断和产量预测[J].光谱学与光谱分析,2013,33(12):3334-3338.
[3]Estimating rice chlorophyll content and leaf nitrogenconcentration with a digital still color camera under natural light[J].PlantMethods,2014,10(1):36.
[4]王方永,王克如,李少昆等.利用数码相机和成像光谱仪估测棉花叶片叶绿素和氮素含量[J].作物学报,2010,36(11):1981-1989.
[5]Rigon J P G,Capuani S,Fernandes D M,et al.A novel method for theestimation of soybean chlorophyll content using a smartphone and imageanalysis[J].Photosynthetica,2016,54(4):559-566.
[6]Su C H,Fu C C,Chang Y C,et al.Simultaneous estimation ofchlorophyll a and lipid contents in microalgae by three-color analysis[J].Biotechnology&Bioengineering,2010,99(4):1034-1039.
[7]Riccardi M,Mele G,Pulvento C,et al.Non-destructive evaluation ofchlorophyll content in quinoa and amaranth leaves by simple and multipleregression analysis of RGB image components[J].Photosynthesis Research,2014,120(3):263-272.
发明内容
本发明的目的是针对现有技术中存在的不足,提供一种冬小麦叶片氮含量估算模型的建立方法,发明人通过分析不同品种、不同种植方案(不同种植密度水平、不同施氮水平)下各基础图像参数对叶片含氮量的表征能力,在各颜色空间内通过各算法构建模型,并通过衡量指标分析比较算法模型的拟合能力与泛化性能。最终确定在3个常用颜色空间(RGB、HSV、L*a*b*空间)融合的多颜色空间下,借助径向基核函数构建ε-SVR模型,具有较高的预测精度,有着更好的泛化性能。基于本发明冬小麦叶片氮含量估算模型可对小麦叶片氮含量进行快速无损检测。
本发明的目的是通过以下技术方案实现的:
一种建立冬小麦叶片氮含量估算模型的方法,包括如下步骤:
步骤(1)、取样:采集小麦冠层图像,测定小麦叶片氮含量;
步骤(2)、利用冠层图像H分量的K均值聚类方法分割小麦冠层图像,提取分割处理后冠层图像的非0像素值的R、G、B、H、S、V、L*、a*、b*分量,分别计算其平均像素值作为基础颜色分量,构建融合三个常用颜色空间的多颜色空间,多颜色空间的基础颜色分量为3种颜色空间下的9种基础颜色分量;
步骤(3)、以图像在多颜色空间下的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,选用径向基核函数,建立ε-SVR模型。
步骤(1)中,小麦冠层图像的获取方法:在小麦拔节期,采用单反相机奥林巴斯E-620,相机设置为自动白平衡;在晴朗天气下,距小麦冠层1m高度、与地面90°拍摄采样。
叶片氮含量的测定方法:小麦冠层图像拍摄当天选择20株长势相近植株进行破坏性取样,采用凯氏定氮法测定小麦的叶片氮含量LNC。
样本点采自不同年份、不同品种、不同种植密度、不同施氮水平。
步骤(2)中,三个常用颜色空间为HSV颜色空间、L*a*b*颜色空间、RGB颜色空间。
步骤(3)中,对于训练样本D={(x1,y1),(x2,y2),…,(xi,yi),…,(xm,ym)},xi是由图像在多颜色空间下的基础颜色分量组成的向量,yi是小麦叶片氮含量;发明人希望能够学习到一个f(x)使得其与y尽可能的接近:f(x)=ωTx+b;
其中,ω和b是待求解的模型参数。在这个模型中,只有当f(x)与y之间完全相同时,损失才为零。
假设能容忍模型输出f(x)与真实输出y最多有ε的偏差,也就是当模型输出和真实输出差别的绝对值大于ε时,才会计算损失。
那么,支持向量回归问题可转化为下式:
其中,C是正则化常数,低的C值使分界面平滑,高的C值通过增加模型自由度以选择更多支持变量来确保所有样品都被正确分类,lε是关于ε的不敏感损失函数,如下式所示:
对ω和b求偏导数,令其为零,可得:
从而得到支持向量回归的对偶问题,如下式:
该过程需要满足KKT条件,即:
可以解得支持向量回归模型为:
其中:
若考虑将样本映射到高维特征空间,则支持向量回归基础式可化为:
f(x)=ωTφ(x)+b;
其中φ(x)表示样本x映射到高维空间后的特征向量。
求解该式会涉及到φ(xi)Tφ(xj),这是样本xi和xj映射到高维空间的内积,由于特征空间的维度较高,通常直接计算该值较为困难。要避开这个困难,通常假设K(xi,xj)满足:
K(xi,xj)=φ(xi)Tφ(xj);
K(xi,xj)即为核函数。本发明使用径向基核函数:
K(xi,xj)=exp(-γ||xi-xj||2),γ>0;
其中,γ是自由参数,定义了单个训练样本的影响大小,可看做被模型选中作为支持向量的样本的影响半径的倒数,值越大影响越小,值越小影响越大。
从而可以得到在高维空间下的ε-SVR模型:
其中:
正则化常数C与核函数参数γ的范围均从0.01到100,以步长10变化;最佳参数以10折交叉验证和网格搜索法确定。
本发明的另一个目的是提供一种冬小麦叶片氮含量估算的方法,包括以下步骤:
步骤(1)、获取待测田块的小麦冠层图像,用基于H分量的K均值聚类方法分割提取小麦冠层图像,计算9种颜色分量R、G、B、H、S、V、L*、a*、b*的平均像素值作为基础颜色分量,将RGB、HSV、L*a*b*这3种常用颜色空间融合为多颜色空间,其基础颜色分量为3种颜色空间下的9种基础分量;
步骤(2)、以多颜色空间下的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,借助径向基核函数,构建ε-SVR模型;
步骤(3)、以多颜色空间的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,将步骤(1)中获得的小麦冠层图像在多颜色空间下的9种颜色基础分量作为自变量代入到ε-SVR模型中,得到预测的小麦叶片氮含量。
本发明法中所有参数的定义及计算方法同上述一种小麦叶片氮含量估算模型的建立方法。
本发明的有益效果:
1、本发明方法通过对不同品种、不同种植密度水平、不同施氮水平,在三种基础颜色空间与融合的多颜色空间下,多元线性回归、支持向量回归和随机森林三种不同的算法建立的模型得到的小麦叶片氮含量预测值与实际值的拟合程度进行分析,表明在多颜色空间下的支持向量回归算法具有较高的预测精度;
2、采样数码相机进行图像的数据采样,方法简便,设备(相机)费用低,方便实际大田环境下的应用;
3、对作物叶片含氮量进行量化评价,可指导农民进行合理施肥,即保证植物所需的氮肥,又可避免过量施肥,提高产量,避免过量施肥污染环境;
4、估计模型的构建方法可以适用于其他作物,但需要参照该方法对不同作物的图像采样数据进行分割并提取其基础颜色分量,从而构建多颜色空间。
附图说明
图1为小麦冠层图像分割效果图;其中,V1N0D2是生选6号、不施氮、行距20cm(3.0×106株/hm2)的试验小区冠层图像;V1N0D1是生选6号、不施氮、行距40cm(1.5×106株/hm2)的试验小区冠层图像。
具体实施方式
1实验与数据
1.1实验设计
2012-2014年,于国家信息农业中心如皋(北纬32°~32°30′、东经120°20′~120°50′)试验基地稻麦轮种田,采用随机裂区设计,12个处理方案,重复3次,共36个试验小区。小区面积35m2(7m×5m),条播,总面积约1080m2。2个供选品种:生选6号(V1)、扬麦18号(V2);2个密度处理:D1(行距40cm,1.5×106株/hm2)、D2(行距20cm,3.0×106株/hm2),株距不变;3个氮肥水平:不施氮(N0)、纯氮150kg·hm-2(N1)、纯氮300kg·hm-2(N2),其中氮肥50%作为基肥施于播种时,50%在拔节期追施(3月中旬)。
1.2数据采集与处理
1.2.1冠层图像采集
拔节期是小麦开花、结实积累养分的重要阶段,此时小麦快速生长,是监测氮素的重要时期。试验取样日期选择了2013年:3月14日(拔节中期),3月22日(拔节中期),4月1日(拔节后期);2014年:3月8日(拔节中期),3月31日(拔节后期),4月15日(拔节后期)。采用单反相机奥林巴斯E-620,相机设置为自动白平衡。在晴朗天气下,距小麦冠层1m高度、与地面90°拍摄采样。每区每次取样1次,36个小区取样36张小麦冠层图像。
1.2.2叶片氮含量的测定
小麦氮含量化学测定与冠层图像采集同期进行。拍摄当天,每个小区选择20株长势相近植株进行破坏性取样,采用凯氏定氮法测定小麦的叶片氮含量LNC。
1.2.3冠层图像处理
大田环境下背景复杂,表现为植株相互遮挡有阴影、光斑,土壤、石砾、杂草造成干扰。这些因素会造成误分割现象。为提高分割精度,利用冠层图像H分量的K均值聚类方法(黄芬,于琪,姚霞等.小麦冠层图像H分量的K均值聚类分割[J].计算机工程与应用,2014,(3):129-134.)分割小麦冠层图像。小麦冠层图像分割效果见图1,H分量的K均值聚类方法将相互遮挡有光斑的叶片与土壤较好分割开,同时叶片边缘细节得以较好地处理。分别计算提取分割处理后冠层图像的非0像素值的R、G、B、H、S、V、L*、a*、b*分量的平均像素值作为基础颜色分量;构建融合三个常用颜色空间(HSV颜色空间、L*a*b*颜色空间、RGB颜色空间)的多颜色空间。
1.3模型的构建
为评估模型的泛化性能,使用嵌套交叉验证:内层交叉验证用来选取模型参数,外层交叉验证用来估计模型泛化性能。嵌套交叉验证中,外层循环遍历数据分为训练集和测试集的所有划分,这样可以减少算法结果对于数据划分的依赖。对于外层循环的每种划分都运行一次网格搜索,每种划分下内层搜索得到的可能是不同的最佳参数。最后,对于外层循环的每种划分,利用最佳参数的模型得到测试集上的评价指标。嵌套交叉验证得到的评价指标可以衡量在网格搜索找到的最佳参数下模型的泛化能力。
1.3.1模型构建方法
采用3种不同的算法建立模型:多元线性回归、支持向量回归和随机森林。模型的自变量为颜色空间的基础颜色分量,因变量为叶片氮含量。
多元线性回归是描述多个自变量与一个因变量的线性模型。其模型参数通过最小二乘法来求取,使得该回归直线到样本的残差平方和最小。
对于训练样本D={(x1,y1),(x2,y2),…,(xi,yi),…,(xm,ym)},yi∈R,传统的回归模型,希望能够学习到一个f(x)使得其与y尽可能的接近:
f(x)=ωTx+b;
其中ω和b是待求解的模型参数。在这个模型中,只有当f(x)与y之间完全相同时,损失才为零。
与传统回归模型不同,支持向量回归假设能容忍模型输出f(x)与真实输出y最多有ε的偏差,也就是当模型输出和真实输出差别的绝对值大于ε时,才会计算损失。
那么,支持向量回归问题可转化为下式:
其中,C是正则化常数,低的C值使分界面平滑,高的C值通过增加模型自由度以选择更多支持变量来确保所有样品都被正确分类,lε是关于ε的不敏感损失函数,如下式所示:
对ω和b求偏导数,令其为零,可得:
从而得到支持向量回归的对偶问题,如下式:
该过程需要满足KKT条件,即:
可以解得支持向量回归模型为:
其中:
若考虑将样本映射到高维特征空间,则支持向量回归基础式可化为:
f(x)=ωTφ(x)+b;
其中,φ(x)表示样本x映射到高维空间后的特征向量。
求解改式会涉及到φ(xi)Tφ(xj),这是样本xi和xj映射到高维空间的内积,由于特征空间的维度较高,通常直接计算该值较为困难。要避开这个困难,通常假设K(xi,xj)满足:
K(xi,xj)=φ(xi)Tφ(xj);
K(xi,xj)即为核函数。本发明使用径向基核函数:
K(xi,xj)=exp(-γ||xi-xj||2),γ>0;
其中,γ是自由参数,定义了单个训练样本的影响大小,可看做被模型选中作为支持向量的样本的影响半径的倒数,值越大影响越小,值越小影响越大。
从而可以得到在高维空间下的ε-SVR模型式:
其中:
随机森林算法是决策树的集成,通过自助法从训练集中采样出n个子集来构建n棵树,每棵树的***节点从含有d个属性的集合中随机选取一个包含k个属性的子集,然后从这个子集中选择一个最优属性用于划分;训练好的n棵回归树,每棵树单独预测,对n个预测值的简单平均作为最终结果;自助采样使得剩下未被抽中的样本组成袋外数据,用于对泛化性能的袋外估计。
1.3.2模型评价标准
选取以下指标衡量算法模型得到的LNC预测值和小麦LNC实际值的拟合程度:
1)决定系数R2:相关系数的平方值,表明预测值解释实际值变差的程度,计算公式为:
2)均方根误差RMSE(Root Mean Square Error):检验预测值和实测值的符合精度,计算公式为:
2结果与分析
2.1颜色分量与氮营养指标相关性
表1为2012-2014年实验实测的小麦叶片氮含量与基础颜色分量之间的相关性分析结果。表1显示,除基础颜色分量a*,其余基础颜色分量与叶片氮含量均显著相关。每个颜色空间中的分量,两个与LNC呈负相关的颜色分量,一个呈正相关。按相关程度来看,颜色分量b*与LNC相关性最高,相关系数绝对值为0.69,B、H、S分量绝对值也都大于0.45,达到显著水平。
表1.颜色分量与叶片氮含量的相关系数
注:“**”表示该颜色变量与叶片氮含量在0.1水平上显著相关。
2.2叶片氮含量估算模型的构建与验证
以冠层图像基础颜色分量为自变量,小麦LNC为因变量,在3个常用颜色空间下构建估算模型。
支持向量回归(SVR)模型的主要参数设置如下:核函数类型为径向基核函数;正则化常数C与核函数参数γ的范围均从0.01到100,以步长10变化;最佳参数以10折交叉验证和网格搜索法确定。
随机森林(RF)模型的回归树数目候选值为{100,200,400,800},用于划分回归树节点的k个属性,其候选值为为避免过度拟合,给随机森林设置两个正则化超参数:回归树的最大深度max_depth和树的叶节点必须有的最小样本数量min_samples_leaf,候选值分别为{4,8,12}、{4,6,8,10}。算法泛化性能的评估采用10×10交叉验证,即外层将数据按照训练集与测试集的方式分为10折,用交叉验证遍历了所有数据的10种划分;对每一种划分采用10折交叉验证寻找最佳参数。
表2显示,L*a*b*颜色空间下,MLR、SVR、RF训练集的R2分别为0.6803、0.8554和0.8748,较RGB颜色空间下,分别提高10.03%,6.58%和4.55%;MLR、SVR、RF训练集的RMSE分别为0.5006,0.3312和0.3131,较RGB空间下分别降低8.48%,15.49%和12.47%。HSV颜色空间下,MLR、SVR、RF训练集的R2分别为0.7108、0.8773和0.8914,较RGB颜色空间下,分别提高14.96%、9.31%和6.54%;RMSE分别为0.4761、0.3100和0.2916,较RGB空间下分别降低12.96%、20.90%和18.48%。结果表明,HSV颜色空间下,各算法训练集的R2最高,RMSE最低;RGB空间下,各算法的R2最低,RMSE最高。
由此可知,HSV颜色空间里各算法模型的学***的小麦LNC状况。高氮水平下小麦叶片叶绿素含量较高,导致导致植物冠层可见光波段的反射特征“绿峰”和“红边”的位置受到影响,表现为绿峰向蓝光方向偏移,波长变短,红边向长波方向移动。这导致了相关性较强的长波段红光反射率降低,相关性弱的短波段蓝光反射率升高。当小麦缺氮时,由于氮素的易运转特性,冠层老叶片的氮素会转向新叶片中,表现为植株下部的叶片先退绿黄化,逐渐扩展到上部叶片。以上因素可能造成RGB空间下,叶色信息的获取存在偏差。HSV和L*a*b*颜色空间通过对RGB基础颜色分量的非线性转换,一定程度上矫正了RGB颜色空间下叶色偏差。同时,RGB空间中包含亮度信息,而大田环境里,光照条件复杂,冠层图像受到土壤、杂草等干扰影响,HSV颜色空间下与小麦冠层叶色相关的基础颜色分量H、S与亮度V分离,L*a*b*颜色空间中色度a*、b*与亮度L*分离,因而减弱了大田复杂背景对冠层反射光谱的影响。
由表2可知,3个颜色空间下RF,训练集上的R2均大于SVR,RMSE均小于SVR,可见RF学习能力优于SVR。分析认为,RF是基于决策树的集成学习模型,该模型不对训练数据做出假设,且该模型参数很多,扩大了模型的自由度,易于紧密拟合训练数据。测试集上,3个颜色空间下SVR的R2均大于RF,RMSE均小于RF,表明SVR算法的泛化能力优于RF,且RF算法有过拟合的情况。分析认为RF虽然在构建每棵树的***节点时引入了随机属性来避免过拟合,但是输入变量小于等于3的情况造成RF会生成一些差异度较小的树,从而降低了算法的泛化性能。SVR是基于最小化结构风险的算法模型,可以得到全局最优解。
结合表2中评价指标的标准差可知,MLR模型在3个颜色空间下,测试集上R2和RMSE浮动区间均包含训练集上的区间,可能出现测试集上评价指标优于训练集。分析认为,缺氮和过氮可能导致小麦生理状态异常,影响了小麦LNC测量的准确性,形成了噪声数据,而基于最小二乘的MLR算法模型对噪声数据比较敏感。SVR模型在RGB空间下,测试集上的评价指标的浮动区间也包含训练集上的区间,说明SVR也受到了噪声数据的干扰,反观RF在3个颜色空间下均没有受到噪声数据的影响,说明其有很好的抗噪能力,可能是RF自助法抽样建立决策树和从输入变量集合中随机抽取变量来***树的节点,这两个随机性的引入使得RF对噪声不敏感。
表2.主要颜色空间下叶片氮含量估算模型交叉验证结果
衡量模型泛化误差的两个方面分别是偏差与方差,偏差用于描述模型的拟合能力,方差用于描述稳定性。训练误差主要是由偏差造成的,而方差通常是由于模型的复杂度相对训练集过高导致,通常体现在测试误差相对训练误差的增量上。MLR测试误差相对训练误差的增量较小,而训练误差较大,造成测试误差相对训练误差的增量与训练误差的比值较小,3个颜色空间下分别为1.97%、0.88%、2.00%;SVR和RF的测试误差相对训练误差的增量较大,而训练误差较小,致使两者的比值较大,分别为17.53%、28.52%、30.65%和47.11%、49.07%、47.24%。由此可知,MLR模型方差的比重小于模型的偏差,即MLR的拟合能力不足,泛化误差主要由偏差造成,而RF模型方差的比重接近于模型的偏差,考虑数据本身噪声也会造成偏差,可认为RF模型的方差主导了泛化误差。SVR模型方差的比重小于模型的偏差,且两者的比值小于RF,但SVR模型的测试集上的评价指标均优于RF,说明RF模型过拟合,SVR模型的正则化常数C对模型的正则化效果优于RF的正则化超参数。
多颜色空间为RGB、HSV、L*a*b*这3种常用颜色空间的融合,其基础颜色分量为3种颜色空间下的9种基础分量。由表3可知,多颜色空间下,MLR、SVR、RF训练集上的R2分别为0.7807、0.8813和0.9152,较HSV颜色空间下分别提高9.83%、0.46%和2.67%;RMSE分别为0.4146、0.3027和0.2578,较HSV颜色空间下分别降低12.92%、2.35%和11.59%。分析认为,品种、光照强度、生育期以及施肥等因素的变化,作物的生长状况和营养状况也会随着改变,导致冠层图像指标对作物的营养状况表达能力的差异,多颜色空间的基础颜色分量包含的信息多于单一颜色空间。
RF测试集上R2为0.7930,高于HSV颜色空间下的0.7372;RMSE为0.3899,低于HSV颜色空间下的0.4347。显然多颜色空间下,RF泛化性能得到提高,其原因可能是RF有多个输入属性值,生成了差异度较大的决策树,使得其泛化性能得到提高。SVR在测试集上的R2略低于RF,且RMSE高于RF,但测试误差相对训练误差的增量为0.092小于RF的0.1321,说明SVR方差小于RF,即SVR比RF要稳定。结合评价指标的方差可知,MLR测试集上R2和RMSE浮动区间均包含训练集上的区间,可能出现测试集上评价指标优于训练集,说明MLR受到噪声数据的干扰。
表3.多颜色空间下叶片氮含量估算模型交叉验证结果
以上研究结果显示,在HSV颜色空间中,各算法模型的估算效果均优于另外两种颜色空间中的效果。融合了RGB、HSV和L*a*b*3种颜色空间的多颜色空间中,各算法模型估算效果优于HSV颜色空间中的效果。
在融合3个颜色空间的多颜色空间下,多元线性回归模型拟合能力不强,偏差主导了泛化误差;随机森林模型拟合能力最强,方差主导了泛化误差;支持向量回归模型正则化常数对模型正则化程度超过随机森林的超参数对模型自身的正则化程度,使得支持向量回归模型有着更好的泛化性能。所以,多颜色空间下的SVR算法模型具有较高的预测精度,可对小麦LNC进行快速无损检测。
Claims (6)
1.一种建立冬小麦叶片氮含量估算模型的方法,其特征在于包括如下步骤:
步骤(1)、取样:采集小麦冠层图像,测定小麦叶片氮含量;
步骤(2)、利用冠层图像H分量的K均值聚类方法分割提取小麦冠层图像,分别计算提取分割处理后冠层图像的非0像素值的R、G、B、H、S、V、L*、a*、b*分量的平均像素值作为基础颜色分量,构建融合三个颜色空间的多颜色空间;
步骤(3)、以图像在多颜色空间下的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,选用径向基核函数,建立ε-SVR模型;
样本D={(x1,y1),(x2,y2),…,(xi,yi),…,(xm,ym)},其中,xi是由图像在多颜色空间下的基础颜色分量组成的向量,yi是小麦叶片氮含量,选用ε-SVR模型:
其中:
K(xi,xj)为径向基核函数;
K(xi,xj)=exp(-γ||xi-xj||2),γ>0(2);
其中,xi、xj为样本D中任意自变量,γ是自由参数,控制映射函数的空间分布复杂度。
2.根据权利要求1所述的冬小麦叶片氮含量估算模型的建立方法,其特征在于步骤(1)中,小麦冠层图像的获取方法:在小麦拔节期,采用单反相机奥林巴斯E-620,相机设置为自动白平衡;在晴朗天气下,距小麦冠层1m高度、与地面90°拍摄采样。
3.根据权利要求1所述的冬小麦叶片氮含量估算模型的建立方法,其特征在于步骤(1)中,叶片氮含量的测定方法:小麦冠层图像拍摄当天对小麦植株进行破坏性取样,采用凯氏定氮法测定小麦的叶片氮含量。
4.根据权利要求1所述的冬小麦叶片氮含量估算模型的建立方法,其特征在于步骤(1)中,样本点采自不同年份、不同品种、不同种植密度、不同施氮水平。
5.根据权利要求1所述的冬小麦叶片氮含量估算模型的建立方法,其特征在于步骤(2)中,三个颜色空间为HSV颜色空间、L*a*b*颜色空间、RGB颜色空间。
6.一种冬小麦叶片氮含量估算的方法,其特征在于包括以下步骤:
步骤(1)、获取待测田块的小麦冠层图像,用基于H分量的K均值聚类方法分割提取小麦冠层图像,计算9种颜色分量R、G、B、H、S、V、L*、a*、b*的平均像素值作为基础特征参数,将3种常用颜色空间RGB、HSV、L*a*b*融合为多颜色空间,多颜色空间的基础颜色分量为3种颜色空间下的9种基础分量;
步骤(2)、以多颜色空间下的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,借助径向基核函数,构建ε-SVR模型;
样本D={(x1,y1),(x2,y2),…,(xi,yi),…,(xm,ym)},其中,xi是由图像在多颜色空间下的基础颜色分量组成的向量,yi是小麦叶片氮含量,选用ε-SVR模型:
其中:
K(xi,xj)为径向基核函数;
K(xi,xj)=exp(-γ||xi-xj||2),γ>0 (2);
其中,xi、xj为样本D中任意自变量,γ是自由参数,控制映射函数的空间分布复杂度;
步骤(3)、以多颜色空间的基础颜色分量作为模型自变量,以叶片氮含量作为模型因变量,将步骤(1)中获得的小麦冠层图像在多颜色空间下的9种颜色基础分量作为自变量代入到ε-SVR模型中,预测得到小麦叶片氮含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911169596.7A CN110874617B (zh) | 2019-11-26 | 2019-11-26 | 一种冬小麦叶片氮含量估算模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911169596.7A CN110874617B (zh) | 2019-11-26 | 2019-11-26 | 一种冬小麦叶片氮含量估算模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110874617A CN110874617A (zh) | 2020-03-10 |
CN110874617B true CN110874617B (zh) | 2023-05-05 |
Family
ID=69718197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911169596.7A Active CN110874617B (zh) | 2019-11-26 | 2019-11-26 | 一种冬小麦叶片氮含量估算模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110874617B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084462B (zh) * | 2020-08-27 | 2023-06-23 | 江苏农林职业技术学院 | 一种基于数字图像的作物氮素营养指标估算模型评价方法 |
CN112329856B (zh) * | 2020-11-06 | 2024-05-14 | 神农智慧农业研究院南京有限公司 | 一种基于光谱与rgb-d融合图像的小麦氮积累量预测方法 |
CN112418073B (zh) * | 2020-11-19 | 2023-10-03 | 安徽农业大学 | 一种基于无人机图像融合特征的小麦植株氮含量估测方法 |
CN112613338B (zh) * | 2020-11-19 | 2023-11-24 | 南京农业大学 | 基于rgb图像融合特征的小麦叶层氮含量估测方法 |
CN113284098B (zh) * | 2021-05-11 | 2024-01-19 | 中国农业大学 | 作物氮含量缺乏分级判别的方法、装置及电子设备 |
CN114119536B (zh) * | 2021-11-22 | 2023-07-04 | 河南省农业科学院农业经济与信息研究所 | 基于云平台技术的作物氮含量快速监测与诊断方法 |
CN115344997B (zh) * | 2022-07-11 | 2024-05-31 | 中国水利水电科学研究院 | 夏玉米植株叶片-冠层-像元尺度氮浓度协同预测方法 |
CN116267454B (zh) * | 2023-02-02 | 2023-11-24 | 安徽农业大学 | 一种筛选不结球白菜耐寒品种的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424637A (zh) * | 2008-12-04 | 2009-05-06 | 浙江大学 | 一种油菜叶片氮素含量遥感估算模型方法 |
US8391565B2 (en) * | 2010-05-24 | 2013-03-05 | Board Of Trustees Of The University Of Arkansas | System and method of determining nitrogen levels from a digital image |
CN102788752B (zh) * | 2012-08-08 | 2015-02-04 | 江苏大学 | 基于光谱技术的农作物内部信息无损检测装置及方法 |
CN105675821B (zh) * | 2016-02-21 | 2018-11-02 | 南京农业大学 | 一种作物氮素营养无损诊断的图像评价指标的建立方法 |
CN109948692B (zh) * | 2019-03-16 | 2020-12-15 | 四川大学 | 基于多颜色空间卷积神经网络及随机森林的计算机生成图片检测方法 |
-
2019
- 2019-11-26 CN CN201911169596.7A patent/CN110874617B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110874617A (zh) | 2020-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110874617B (zh) | 一种冬小麦叶片氮含量估算模型的建立方法 | |
Qiao et al. | UAV-based chlorophyll content estimation by evaluating vegetation index responses under different crop coverages | |
Song et al. | Wavelength selection and spectral discrimination for paddy rice, with laboratory measurements of hyperspectral leaf reflectance | |
Zhou et al. | Diagnosis of winter-wheat water stress based on UAV-borne multispectral image texture and vegetation indices | |
CN109187441B (zh) | 基于冠层光谱信息的夏玉米含氮量监测模型的构建方法 | |
da Silva Junior et al. | Soybean varieties discrimination using non-imaging hyperspectral sensor | |
Banerjee et al. | High-throughput phenotyping using digital and hyperspectral imaging-derived biomarkers for genotypic nitrogen response | |
CN102788752B (zh) | 基于光谱技术的农作物内部信息无损检测装置及方法 | |
CN105675821B (zh) | 一种作物氮素营养无损诊断的图像评价指标的建立方法 | |
Mao et al. | Nondestructive measurement of total nitrogen in lettuce by integrating spectroscopy and computer vision | |
Zhang et al. | Estimating wheat yield by integrating the WheatGrow and PROSAIL models | |
Wen et al. | Estimating leaf nitrogen concentration considering unsynchronized maize growth stages with canopy hyperspectral technique | |
Singh et al. | Assessment of growth, leaf N concentration and chlorophyll content of sweet sorghum using canopy reflectance | |
CN112287886B (zh) | 基于高光谱图像融合图谱特征的小麦植株氮含量估测方法 | |
CN110189793A (zh) | 基于高光谱的小麦氮肥生理利用率估测模型构建及不同氮效率小麦品种分类 | |
Huang et al. | Assimilating remotely sensed information with the wheatgrow model based on the ensemble square root filter for improving regional wheat yield forecasts | |
Chen et al. | Model for estimation of total nitrogen content in sandalwood leaves based on nonlinear mixed effects and dummy variables using multispectral images | |
Bei et al. | Leaf chlorophyll content retrieval of wheat by simulated RapidEye, Sentinel-2 and EnMAP data | |
Brambilla et al. | Application of a low-cost RGB sensor to detect basil (Ocimum basilicum L.) nutritional status at pilot scale level | |
Wang et al. | Estimation of rice aboveground biomass by combining canopy spectral reflectance and unmanned aerial vehicle-based red green blue imagery data | |
Shu et al. | Assessing maize lodging severity using multitemporal UAV-based digital images | |
Fan et al. | Using an optimized texture index to monitor the nitrogen content of potato plants over multiple growth stages | |
Cheng et al. | Hyperspectral and imagery integrated analysis for vegetable seed vigor detection | |
Su et al. | Grain yield prediction using multi-temporal UAV-based multispectral vegetation indices and endmember abundance in rice | |
Wang et al. | Intelligent detection of hard seeds of snap bean based on hyperspectral imaging |
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 |