CN113049530A - 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 - Google Patents
一种基于近红外高光谱的单粒玉米种子水分含量检测方法 Download PDFInfo
- Publication number
- CN113049530A CN113049530A CN202110286296.8A CN202110286296A CN113049530A CN 113049530 A CN113049530 A CN 113049530A CN 202110286296 A CN202110286296 A CN 202110286296A CN 113049530 A CN113049530 A CN 113049530A
- Authority
- CN
- China
- Prior art keywords
- sample
- corn seed
- characteristic
- image data
- moisture content
- 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
- 240000008042 Zea mays Species 0.000 title claims abstract description 72
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 title claims abstract description 72
- 235000002017 Zea mays subsp mays Nutrition 0.000 title claims abstract description 72
- 235000005822 corn Nutrition 0.000 title claims abstract description 72
- 238000001514 detection method Methods 0.000 title claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 230000003595 spectral effect Effects 0.000 claims abstract description 32
- 238000012549 training Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000004927 fusion Effects 0.000 claims abstract description 21
- 238000007637 random forest analysis Methods 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 238000003066 decision tree Methods 0.000 claims description 21
- 239000002245 particle Substances 0.000 claims description 21
- 238000012937 correction Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 9
- 238000000701 chemical imaging Methods 0.000 claims description 6
- 210000001161 mammalian embryo Anatomy 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000002835 absorbance Methods 0.000 claims description 3
- 238000000149 argon plasma sintering Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000513 principal component analysis Methods 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 235000011888 snacks Nutrition 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000010899 nucleation Methods 0.000 abstract description 3
- 238000007781 pre-processing Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000009331 sowing Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000001035 drying Methods 0.000 description 3
- 239000002994 raw material Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000007500 overflow downdraw method Methods 0.000 description 2
- 208000003643 Callosities Diseases 0.000 description 1
- 206010020649 Hyperkeratosis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 235000020985 whole grains Nutrition 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
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Chemical & Material Sciences (AREA)
- Software Systems (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
Abstract
本发明提供一种基于近红外高光谱的单粒玉米种子水分含量检测方法,步骤为:1)同批次采集不同水分含量的N个玉米种子样本,并测量样本的水分含量数据yi以及高光谱图像数据xi;2)对高光谱图像数据xi进行预处理得到预处理光谱图像数据x′i;3)提取预处理光谱图像数据x′i的特征光谱以及高光谱图像数据xi的纹理特征,将提取出的特征光谱和纹理特征拼接得到光谱图像融合矩阵Pi;4)将样本集划分为训练集与测试集,利用随机森林算法构建水分预测模型;5)待测样本预测:根据随机森林模型对待测样本的水分含量进行自适应加权预测。本发明利用近红外高光谱技术,将光谱和图像数据融合,构建单粒玉米种子水分含量检测模型,对于玉米精量播种技术的应用具有重要意义。
Description
技术领域
本发明涉及玉米种子水分检测领域,特别是一种基于近红外高光谱的单粒玉米种子水分含量检测方法。
背景技术
玉米是目前我国种植面积最广、产量第一的农作物,不仅是“饲料之王”,还是重要的工业原料。近年来随着由于美国先锋公司“单粒播”玉米种子的推出,以及欧美等发达国家的玉米机械化单粒精量播种技术的引入,我国传统玉米播种模式发生了深刻的改变。单粒精量播种技术较传统播种模式更有利于机械化操作,省工、省种、高产,但是对每颗种子都提出了更高的检测需求,其中对单粒播种的玉米种子水分要求不能高于14%,否则会导致种子呼吸增大、消耗养分、活力下降。根据我国国标GB4404.1-2008,水分是我国农作物种子质量四大必检项目之一。因此单粒玉米种子水分的快速、高通量、无损检测对于单粒精量播种具有重要的现实意义。
农作物种子检验规程水分测定(GB/T 3543.6-1995)中规定了种子水分测量方法有低恒温烘干法、高温烘干法、高水分预先烘干法,上述传统检测方法虽然发展早,但普遍存在试样破坏性、耗时长、无法单颗测定等问题,无法满足当前种业市场快速发展的需求,因此急需建立一种快速、无损、稳健、高精度的检测方法。
近年来近红外高光谱技术以其快速、无损、绿色分析特点在种子水分检测领域表现出极大的应用潜力[20-25]。我国在2010年发布了玉米水分含量近红外法测定国标GB/T24900-2010[12]用于玉米种子批水分测定,表明近红外光谱技术在玉米种子批水分检测领域具有实际应用可行性,但是目前尚无针对单粒种子水分近红外高光谱检测的方法或标准等。
发明内容
本发明的目的就是提供一种基于近红外高光谱的单粒玉米种子水分含量检测方法,它可以用于单粒玉米种子水分含量的检测。
本发明的目的是通过这样的技术方案实现的,具体步骤如下:
1)数据采集:同批次采集不同水分含量的玉米种子样本,去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本,对N个玉米种子样本进行水分测定,得到水分含量数据yi,利用高光谱成像系统采集玉米种子样本的高光谱图像数据xi;
2)数据处理:采用多元散射校正MSC方法消除步骤1)中高光谱图像数据xi的光散射、基线漂移和光程变化的干扰,得到预处理光谱图像数据x′i;
3)特征提取:采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱,采用灰度共生矩阵提取特征波长下的高光谱图像数据xi的纹理特征,所述纹理特征包括对比度、能量、熵、逆方差、相关性;将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi;
4)构建模型:采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本集划分为训练集与测试集,利用基于集成策略的随机森林算法构建水分预测模型;
5)待测样本预测:根据步骤4)中构建的随机森林模型对待测样本的水分含量进行自适应加权预测。
进一步,步骤1)中数据采集的具体步骤为:
1-1)同批次采集不同水分含量的玉米种子样本,对种子进行初步的筛选去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本;
1-2)利用高光谱成像系统,将N个玉米种子样本排列整齐,胚朝上、胚朝下各扫描一次,采集高光谱图像数据xi,并进行黑白板校正,采用主成分分析确定感兴趣区域;
1-3)对N个玉米种子样本进行水分测定。
进一步,步骤3)特征提取的具体步骤如下:
3-1)采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱;
3-2)提取特征波长下的高光谱图像数据xi的灰度图像,将灰度图像进行直方图均衡化的处理,计算图像特征值的步距和方向,步距选择d=1,将中心像素与其相邻像素点心与其相邻像素点直接做比较运算,灰度共生矩阵的方向选择为为0°、45°、90°、135°四个方向,并求出这四个方向特征值的平均值作为最终特征值的共生矩阵,将共生矩阵归一化后,通过计算对比度、能量、熵、逆方差、相关性参数来表征图像纹理特征;
3-3)将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi。
进一步,步骤4)中构建模型的具体步骤为:
4-1)采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本按照3:1的比例划分为训练集与测试集;
4-2)对于训练集D={(y1,P1),(y2,P2),...,(yi,Pi)},利用Bootstrap方法从训练集中采取随机有放回的抽取N个子训练集,抽取的子训练集个数与随机森林中树模型数量相同;
4-3)对单个子训练集应用决策树的计算方法,即从数据的n维特征波长中无放回的随机抽取mtry维特征波长,即用于二叉树的特征波长个数;
4-4)对于一个特征值序列进行各个点位的二分,计算每一次二分之后的两个点位的残差平方和:
式(1)中,j为特征波长,s为分割点,c1和c2为叶子节点的样本值,R1和R2为分割后的子区域,xi为样本在每个特征波长下的吸光度,yi为样本的水分值;
选择最小的那个特征值作为候选分割点,在针对候选的分割点的选出残差平方和最小的作为二叉树当前层的节点,然后提出到选定的特征,再从两个子区域中继续该操作分别获得新的子节点,继续分割,直至满足预设停止条件;
4-5)重复步骤4-3)至步骤4-4),最后将输入空间划分为M个区域R1 R2…RM,生成N棵决策树,构成随机森林预测模型:
进一步,步骤5)中对待测样本进行预测的具体步骤如下:
5-1)首先根据欧式距离公式计算待测样本与随机森林模型中第i棵决策树样本集中心的距离di:
5-2)根据步骤5-1)中的欧式距离确定在待测样本时每棵决策树的权重;第i棵决策树的权重ai:
式(4)中,i∈N;
5-3)对步骤5-2)的每棵决策树的权重进行加权,输出待预测样本水分的预测值Y:
由于采用了上述技术方案,本发明具有如下的优点:
1、本发明利用近红外高光谱技术,将光谱和图像数据融合,构建基于集成策略的自适应权重的随机森林算法的水分预测模型,对于玉米精量播种技术的应用具有重要意义;2、通过本发明中的制备不同水分含量的样本,增大了样本水分梯度范围,提高了模型的泛化能力;3、光谱预处理减少了非线性因素的影响,提高了信噪比;光谱降维大大降低了数据冗余性,减少模型的输入变量,提高了单粒玉米种子水分检测的效率;4、通过特征层融合方法,将提取出的特征波长作为光谱特征和提取出的图像特征进行融合建模,高光谱的光谱信息可以提供不同品种、不同水分含量的玉米种子的内部差异,高光谱的图像信息可以提供不同品种、不同水分含量的玉米种子的外部属性之间的差异,融合高光谱的光谱信息和图像信息可以更有效的检测玉米种子水分含量;5、本发明基于自适应权重策略,在预测阶段对不同训练样本集赋予不同的权重进行预测,提高了预测精度。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
本发明的附图说明如下。
图1为本发明的方法流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
如图1所示的一种基于近红外高光谱的单粒玉米种子水分含量检测方法,其特征在于,具体步骤如下:
1)数据采集:同批次采集不同水分含量的玉米种子样本,去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本,对N个玉米种子样本进行水分测定,得到水分含量数据yi,利用高光谱成像系统采集玉米种子样本的高光谱图像数据xi,具体步骤为:
1-1)同批次采集不同水分含量的玉米种子样本,对种子进行初步的筛选去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本;
1-2)利用高光谱成像系统,将N个玉米种子样本排列整齐,胚朝上、胚朝下各扫描一次,采集高光谱图像数据xi,并进行黑白板校正,采用主成分分析确定感兴趣区域;
1-3)依据《GB/T 10362-2008.粮油检验玉米水分测定》中“整粒试样”测定方法对样本进行水分测定。
2)数据处理:采用多元散射校正MSC方法消除步骤1)中高光谱图像数据xi的光散射、基线漂移和光程变化的干扰,得到预处理光谱图像数据x′i;
在发明实例中,由于复杂样本光谱信号往往会受到杂散光、噪声、基线漂移等因素的干扰,样本粒径是影响光谱测量的一个重要参数,随着样本粒径的增加,所测光谱的重现性变差,光谱的变动性随粒径的增加呈指数形式增加;每颗玉米籽粒形状、直径都不相同,所测光谱差异性变大,造成测量误差无法消除,从而影响最终的定量分析结果;多元散射校正可以有效消除近红外漫反射光谱中由于样品的镜面反射及不均匀造成的噪声,消除光谱的基线漂移现象及光谱的不重复性,多元散射校正方法是现阶段多波长标定建模常用的一种数据处理方法,经过散射校正后得到的光谱数据可以有效地消除散射的影响,增强了与成分含量相关的光谱吸收信息,因此基于多元散射校正相对其他预处理方法而言,建立的单粒种子水分模型性能较好。
3)特征提取:采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱,采用灰度共生矩阵提取特征波长下的高光谱图像数据xi的纹理特征,所述纹理特征包括对比度、能量、熵、逆方差、相关性;将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi;具体步骤如下:
3-1)采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱;
3-2)提取特征波长下的高光谱图像数据xi的灰度图像,将灰度图像进行直方图均衡化的处理,计算图像特征值的步距和方向,步距选择d=1,将中心像素与其相邻像素点心与其相邻像素点直接做比较运算,灰度共生矩阵的方向选择为为0°、45°、90°、135°四个方向,并求出这四个方向特征值的平均值作为最终特征值的共生矩阵,将共生矩阵归一化后,通过计算对比度、能量、熵、逆方差、相关性参数来表征图像纹理特征;
3-3)将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi。
在本发明实例中,不同水分含量、不同品种的玉米种子的纹理、饱满度等外部属性可能因为内部化学成分含量的微小区别呈现变化,而这些变化仅从个人感官无法发现。通过高光谱的图像信息可以表征一定差异性,在近红外波段下,成分不同对光谱吸收也不同,在某个特定波长下的图像对某个缺陷会有较显著的反映,可以有效地快速检测出玉米种子的外部属性如尺寸、形状、颜色、纹理和缺陷,充分利用高光谱的图像信息;通过特征层融合方法,将提取出的特征波长作为光谱特征和提取出的图像特征进行融合建模,高光谱的光谱信息可以提供不同品种、不同水分含量的玉米种子的内部差异,高光谱的图像信息可以提供不同品种、不同水分含量的玉米种子的外部属性之间的差异,因此融合高光谱的光谱信息和图像信息可以更有效的检测玉米种子水分含量。
4)构建模型:采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本集划分为训练集与测试集,利用基于集成策略的随机森林算法构建水分预测模型,具体步骤为:
4-1)采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本按照3:1的比例划分为训练集与测试集;
4-2)对于训练集D={(y1,P1),(y2,P2),...,(yi,Pi)},利用Bootstrap方法从训练集中采取随机有放回的抽取N个子训练集,抽取的子训练集个数与随机森林中树模型数量相同;
4-3)对单个子训练集应用决策树的计算方法,即从数据的n维特征波长中无放回的随机抽取mtry维特征波长,即用于二叉树的特征波长个数;
4-4)对于一个特征值序列进行各个点位的二分,计算每一次二分之后的两个点位的残差平方和:
式(1)中,j为特征波长,s为分割点,c1和c2为叶子节点的样本值,R1和R2为分割后的子区域,xi为样本在每个特征波长下的吸光度,yi为样本的水分值;
选择最小的那个特征值作为候选分割点,在针对候选的分割点的选出残差平方和最小的作为二叉树当前层的节点,然后提出到选定的特征,再从两个子区域中继续该操作分别获得新的子节点,继续分割,直至满足预设停止条件;
4-5)重复步骤4-3)至步骤4-4),最后将输入空间划分为M个区域R1 R2…RM,生成N棵决策树,构成随机森林预测模型:
5)待测样本预测:根据步骤4)中构建的随机森林模型对待测样本的水分含量进行自适应加权预测,具体步骤如下:
5-1)首先根据欧式距离公式计算待测样本与随机森林模型中第i棵决策树样本集中心的距离di:
5-2)根据步骤5-1)中的欧式距离确定在待测样本时每棵决策树的权重;第i棵决策树的权重ai:
式(4)中,i∈N;
5-3)对步骤5-2)的每棵决策树的权重进行加权,输出待预测样本水分的预测值Y:
在本发明实例中,由于在训练集中类别不平衡的情况下,对某些贡献度较大的样本数据,通常需要有小的训练误差。因此当预测未知样本X的水分Y时,基于自适应权重策略,对不同训练样本集赋予不同的权重进行预测。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (5)
1.一种基于近红外高光谱的单粒玉米种子水分含量检测方法,其特征在于,具体步骤如下:
1)数据采集:同批次采集不同水分含量的玉米种子样本,去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本,对N个玉米种子样本进行水分测定,得到水分含量数据yi,利用高光谱成像系统采集玉米种子样本的高光谱图像数据xi;
2)数据处理:采用多元散射校正MSC方法消除步骤1)中高光谱图像数据xi的光散射、基线漂移和光程变化的干扰,得到预处理光谱图像数据x′i;
3)特征提取:采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱,采用灰度共生矩阵提取特征波长下的高光谱图像数据xi的纹理特征,所述纹理特征包括对比度、能量、熵、逆方差、相关性;将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi;
4)构建模型:采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本集划分为训练集与测试集,利用基于集成策略的随机森林算法构建水分预测模型;
5)待测样本预测:根据步骤4)中构建的随机森林模型对待测样本的水分含量进行自适应加权预测。
2.如权利要求1所述的一种基于近红外高光谱的单粒玉米种子水分含量检测方法,其特征在于,步骤1)中数据采集的具体步骤为:
1-1)同批次采集不同水分含量的玉米种子样本,对种子进行初步的筛选去掉虫蛀孔洞颗粒、霉变颗粒、破损颗粒获得N个玉米种子样本;
1-2)利用高光谱成像系统,将N个玉米种子样本排列整齐,胚朝上、胚朝下各扫描一次,采集高光谱图像数据xi,并进行黑白板校正,采用主成分分析确定感兴趣区域;
1-3)对N个玉米种子样本进行水分测定。
3.如权利要求1所述的一种基于近红外高光谱的单粒玉米种子水分含量检测方法,其特征在于,步骤3)特征提取的具体步骤如下:
3-1)采用去噪自编码器DAE提取步骤2)中预处理光谱图像数据x′i的特征光谱;
3-2)提取特征波长下的高光谱图像数据xi的灰度图像,将灰度图像进行直方图均衡化的处理,计算图像特征值的步距和方向,步距选择d=1,将中心像素与其相邻像素点心与其相邻像素点直接做比较运算,灰度共生矩阵的方向选择为为0°、45°、90°、135°四个方向,并求出这四个方向特征值的平均值作为最终特征值的共生矩阵,将共生矩阵归一化后,通过计算对比度、能量、熵、逆方差、相关性参数来表征图像纹理特征;
3-3)将提取出的特征光谱和纹理特征拼接进行特征层融合,得到光谱图像融合矩阵Pi。
4.如权利要求1所述的一种基于近红外高光谱的单粒玉米种子水分含量检测方法,其特征在于,步骤4)中构建模型的具体步骤为:
4-1)采用SPXY算法及玉米种子样本的光谱图像融合矩阵Pi、水分含量数据yi将玉米种子样本按照3:1的比例划分为训练集与测试集;
4-2)对于训练集D={(y1,P1),(y2,P2),...,(yi,Pi)},利用Bootstrap方法从训练集中采取随机有放回的抽取N个子训练集,抽取的子训练集个数与随机森林中树模型数量相同;
4-3)对单个子训练集应用决策树的计算方法,即从数据的n维特征波长中无放回的随机抽取mtry维特征波长,即用于二叉树的特征波长个数;
4-4)对于一个特征值序列进行各个点位的二分,计算每一次二分之后的两个点位的残差平方和:
式(1)中,j为特征波长,s为分割点,c1和c2为叶子节点的样本值,R1和R2为分割后的子区域,xi为样本在每个特征波长下的吸光度,yi为样本的水分值;
选择最小的那个特征值作为候选分割点,在针对候选的分割点的选出残差平方和最小的作为二叉树当前层的节点,然后提出到选定的特征,再从两个子区域中继续该操作分别获得新的子节点,继续分割,直至满足预设停止条件;
4-5)重复步骤4-3)至步骤4-4),最后将输入空间划分为M个区域R1 R2 … RM,生成N棵决策树,构成随机森林预测模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110286296.8A CN113049530B (zh) | 2021-03-17 | 2021-03-17 | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110286296.8A CN113049530B (zh) | 2021-03-17 | 2021-03-17 | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113049530A true CN113049530A (zh) | 2021-06-29 |
CN113049530B CN113049530B (zh) | 2022-07-26 |
Family
ID=76513030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110286296.8A Active CN113049530B (zh) | 2021-03-17 | 2021-03-17 | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113049530B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210169003A1 (en) * | 2019-12-04 | 2021-06-10 | Monsanto Technology Llc | Combine Harvesters For Use In Harvesting Corn, And Related Methods |
CN113433081A (zh) * | 2021-07-06 | 2021-09-24 | 江南大学 | 基于动态神经网络的果蔬干燥过程中品质检测方法及系统 |
CN113762150A (zh) * | 2021-09-07 | 2021-12-07 | 北京工商大学 | 一种尾矿库特征分析及模型构建方法及系统 |
CN113777104A (zh) * | 2021-08-24 | 2021-12-10 | 北京工商大学 | 一种单粒玉米种子高光谱检测成熟度的方法 |
CN114486799A (zh) * | 2022-04-15 | 2022-05-13 | 广东省农业科学院动物科学研究所 | 一种家禽肠道丁酸梭菌检测方法及系统 |
CN115345239A (zh) * | 2022-08-17 | 2022-11-15 | 无锡迅杰光远科技有限公司 | 样品含量识别方法、装置和存储介质 |
CN116843998A (zh) * | 2023-08-29 | 2023-10-03 | 四川省分析测试服务中心 | 一种光谱样本加权方法及系统 |
CN117074353A (zh) * | 2023-08-18 | 2023-11-17 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
CN118032677B (zh) * | 2024-04-09 | 2024-06-07 | 甘肃农业大学 | 当归冻干产品加工质量测试方法和系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103487411A (zh) * | 2013-09-01 | 2014-01-01 | 西北大学 | 一种随机森林算法结合激光诱导击穿光谱识别钢材牌号的方法 |
CN110108644A (zh) * | 2019-05-15 | 2019-08-09 | 北京师范大学 | 一种基于深度级联森林和高光谱图像的玉米品种鉴别方法 |
CN110163101A (zh) * | 2019-04-17 | 2019-08-23 | 湖南省中医药研究院 | 中药材种子区别及等级快速判别方法 |
WO2019174572A1 (zh) * | 2018-03-12 | 2019-09-19 | 浙江大学 | 一种利用种胚谱图特征波段识别谷物种子冻伤情况的方法 |
CN110728302A (zh) * | 2019-09-10 | 2020-01-24 | 武汉纺织大学 | 一种基于HSV和Lab颜色空间的色纺织物组织识别的方法 |
CN111272668A (zh) * | 2020-01-22 | 2020-06-12 | 中国农业科学院农产品加工研究所 | 小麦品种鉴别模型的构建方法 |
CN111598786A (zh) * | 2019-11-28 | 2020-08-28 | 南京航空航天大学 | 一种基于深度去噪自编码网络的高光谱图像解混方法 |
CN112418653A (zh) * | 2020-11-19 | 2021-02-26 | 重庆邮电大学 | 一种基于机器学习算法的携号转网潜客识别系统及方法 |
CN112446522A (zh) * | 2019-09-02 | 2021-03-05 | 中国林业科学研究院资源信息研究所 | 面对多尺度分割的产草量估算方法、装置及存储介质 |
-
2021
- 2021-03-17 CN CN202110286296.8A patent/CN113049530B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103487411A (zh) * | 2013-09-01 | 2014-01-01 | 西北大学 | 一种随机森林算法结合激光诱导击穿光谱识别钢材牌号的方法 |
WO2019174572A1 (zh) * | 2018-03-12 | 2019-09-19 | 浙江大学 | 一种利用种胚谱图特征波段识别谷物种子冻伤情况的方法 |
CN110163101A (zh) * | 2019-04-17 | 2019-08-23 | 湖南省中医药研究院 | 中药材种子区别及等级快速判别方法 |
CN110108644A (zh) * | 2019-05-15 | 2019-08-09 | 北京师范大学 | 一种基于深度级联森林和高光谱图像的玉米品种鉴别方法 |
CN112446522A (zh) * | 2019-09-02 | 2021-03-05 | 中国林业科学研究院资源信息研究所 | 面对多尺度分割的产草量估算方法、装置及存储介质 |
CN110728302A (zh) * | 2019-09-10 | 2020-01-24 | 武汉纺织大学 | 一种基于HSV和Lab颜色空间的色纺织物组织识别的方法 |
CN111598786A (zh) * | 2019-11-28 | 2020-08-28 | 南京航空航天大学 | 一种基于深度去噪自编码网络的高光谱图像解混方法 |
CN111272668A (zh) * | 2020-01-22 | 2020-06-12 | 中国农业科学院农产品加工研究所 | 小麦品种鉴别模型的构建方法 |
CN112418653A (zh) * | 2020-11-19 | 2021-02-26 | 重庆邮电大学 | 一种基于机器学习算法的携号转网潜客识别系统及方法 |
Non-Patent Citations (1)
Title |
---|
付琼莹等: "高光谱影像的DAE分类", 《测绘科学技术学报》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210169003A1 (en) * | 2019-12-04 | 2021-06-10 | Monsanto Technology Llc | Combine Harvesters For Use In Harvesting Corn, And Related Methods |
US11968925B2 (en) * | 2019-12-04 | 2024-04-30 | Monsanto Technology Llc | Methods for harvesting corn using combine harvesters |
CN113433081A (zh) * | 2021-07-06 | 2021-09-24 | 江南大学 | 基于动态神经网络的果蔬干燥过程中品质检测方法及系统 |
CN113433081B (zh) * | 2021-07-06 | 2022-05-24 | 江南大学 | 基于动态神经网络的果蔬干燥过程中品质检测方法及系统 |
CN113777104A (zh) * | 2021-08-24 | 2021-12-10 | 北京工商大学 | 一种单粒玉米种子高光谱检测成熟度的方法 |
CN113777104B (zh) * | 2021-08-24 | 2023-08-29 | 北京工商大学 | 一种单粒玉米种子高光谱检测成熟度的方法 |
CN113762150B (zh) * | 2021-09-07 | 2024-04-02 | 北京工商大学 | 一种尾矿库特征分析及模型构建方法及系统 |
CN113762150A (zh) * | 2021-09-07 | 2021-12-07 | 北京工商大学 | 一种尾矿库特征分析及模型构建方法及系统 |
CN114486799A (zh) * | 2022-04-15 | 2022-05-13 | 广东省农业科学院动物科学研究所 | 一种家禽肠道丁酸梭菌检测方法及系统 |
CN115345239B (zh) * | 2022-08-17 | 2023-10-13 | 无锡迅杰光远科技有限公司 | 样品含量识别方法、装置和存储介质 |
CN115345239A (zh) * | 2022-08-17 | 2022-11-15 | 无锡迅杰光远科技有限公司 | 样品含量识别方法、装置和存储介质 |
CN117074353A (zh) * | 2023-08-18 | 2023-11-17 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
CN117074353B (zh) * | 2023-08-18 | 2024-05-03 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
CN116843998B (zh) * | 2023-08-29 | 2023-11-14 | 四川省分析测试服务中心 | 一种光谱样本加权方法及系统 |
CN116843998A (zh) * | 2023-08-29 | 2023-10-03 | 四川省分析测试服务中心 | 一种光谱样本加权方法及系统 |
CN118032677B (zh) * | 2024-04-09 | 2024-06-07 | 甘肃农业大学 | 当归冻干产品加工质量测试方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113049530B (zh) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113049530B (zh) | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 | |
Cavallo et al. | Non-destructive and contactless quality evaluation of table grapes by a computer vision system | |
Williams et al. | Maize kernel hardness classification by near infrared (NIR) hyperspectral imaging and multivariate data analysis | |
Yao et al. | Nondestructive detection for egg freshness grade based on hyperspectral imaging technology | |
Shenderey et al. | NIRS detection of moldy core in apples | |
CN109470648B (zh) | 一种单籽粒作物不完善粒快速无损判定方法 | |
CN106706553A (zh) | 一种快速无损测定玉米单籽粒直链淀粉含量的方法 | |
CN108680515B (zh) | 一种单粒水稻直链淀粉定量分析模型构建及其检测方法 | |
CN110132856B (zh) | 小麦赤霉病染病籽粒识别的光谱病害指数构建及应用 | |
Wang et al. | Discrimination of cracked soybean seeds by near-infrared spectroscopy and random forest variable selection | |
Lu et al. | Hyperspectral imaging with machine learning to differentiate cultivars, growth stages, flowers, and leaves of industrial hemp (Cannabis sativa L.) | |
CN116559111A (zh) | 一种基于高光谱成像技术的高粱品种识别方法 | |
Rathna Priya et al. | Characterising corn grain using infrared imaging and spectroscopic techniques: a review | |
Miart et al. | MuSeeQ, a novel supervised image analysis tool for the simultaneous phenotyping of the soluble mucilage and seed morphometric parameters | |
Cheng et al. | Hyperspectral and imagery integrated analysis for vegetable seed vigor detection | |
Liu et al. | Rapid determination of rice protein content using near-infrared spectroscopy coupled with feature wavelength selection | |
CN106706554A (zh) | 一种快速无损测定玉米单穗籽粒直链淀粉含量的方法 | |
CN102519903B (zh) | 一种近红外光谱测定双孢蘑菇白度值的方法 | |
Zhang et al. | A transfer learning method for spectral model of moldy apples from different origins | |
Sun et al. | Nondestructive identification of barley seeds varieties using hyperspectral data from two sides of barley seeds | |
CN110609011A (zh) | 单籽粒玉米种子淀粉含量近红外高光谱检测方法及系统 | |
Xu et al. | Hyperspectral imaging with machine learning for non-destructive classification of Astragalus membranaceus var. mongholicus, Astragalus membranaceus, and similar seeds | |
CN116297236A (zh) | 一种基于高光谱的单粒玉米种子活力鉴别方法及装置 | |
CN114813593A (zh) | 基于高光谱成像技术的糟醅总酸含量的检测方法 | |
CN115661813A (zh) | 基于高光谱成像技术的多粮混杂碎颗粒快速分类检测方法 |
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 |