CN105913056B - 一种基于工业相机特征提取的润后水渍烟叶识别方法 - Google Patents
一种基于工业相机特征提取的润后水渍烟叶识别方法 Download PDFInfo
- Publication number
- CN105913056B CN105913056B CN201610195473.0A CN201610195473A CN105913056B CN 105913056 B CN105913056 B CN 105913056B CN 201610195473 A CN201610195473 A CN 201610195473A CN 105913056 B CN105913056 B CN 105913056B
- Authority
- CN
- China
- Prior art keywords
- tobacco leaf
- matrix
- tobacco leaves
- original
- water
- 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
- 241000208125 Nicotiana Species 0.000 title claims abstract description 204
- 235000002637 Nicotiana tabacum Nutrition 0.000 title claims abstract description 204
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000000605 extraction Methods 0.000 title claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title abstract description 17
- 239000013598 vector Substances 0.000 claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 25
- 230000002159 abnormal effect Effects 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 69
- 238000007781 pre-processing Methods 0.000 claims description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000007787 solid Substances 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 abstract 1
- 238000002834 transmittance Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000010025 steaming Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000003763 carbonization Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000012569 chemometric method Methods 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/10—Image acquisition
- G06V10/12—Details of acquisition arrangements; Constructional details thereof
- G06V10/14—Optical characteristics of the device performing the acquisition or on the illumination arrangements
- G06V10/141—Control of illumination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/56—Extraction of image or video features relating to colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/68—Food, e.g. fruit or vegetables
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
- Manufacture Of Tobacco Products (AREA)
Abstract
本发明公开了一种基于工业相机特征提取的润后水渍烟叶识别方法,所述方法包括如下步骤:1)采集在线烟叶样本的所有原始烟叶图像,将烟叶样本人工分类为正常烟叶、水渍烟叶、近正常烟叶和近水渍烟叶,并获取所有原始烟叶图像的信息;2)对所述原始烟叶图像的信息进行处理获得所有原始烟叶图像的特征向量;3)对所述特征向量进行SNV预处理以消除异常数值信息;4)对预处理后的特征采取有监督模式的二维投影处理获得二维投影向量。快速的解决了一润后水渍烟叶的快速识别以及比例的计算,为通过工艺的控制来减少水渍烟叶的比例打下了坚实的基础,积累了原始的水渍烟叶识别数据。
Description
技术领域
本发明涉及一种图像特征信息提取及识别领域,具体涉及一种基于工业相机特征提取的润后水渍烟叶识别方法。
背景技术
在润叶工序受高温高湿的影响,烟叶容易产生烟叶的水渍烟及蒸片,在复烤环节很难把水渍烟及蒸片调控到合适的水分,造成水渍片水分偏大,在片烟装箱过程中,水渍片温度及水分都较高,造成局部烟叶水分及温度过高,在保管过程中容易产生霉变及碳化,造成损失。水渍烟叶的识别一般根据烟叶的叶面和叶背是一种颜色基调,对着阳光通透性较强。但是这种人工的操作一方面不能在生产线上完成对水渍烟叶的事实识别以及比例的测算,另外一方面,现场没有光源,使得根据烟叶自身的透光性来判别生产线的烟叶是否为水渍烟叶不具备可行性;随着近年来在烟叶外观分析领域工业相机的逐步应用,颜色特征提取以及图像处理技术的发展,对水渍烟叶在烟叶的透光性的研究,以及图像的特征抽提上来完成对水渍烟叶的识别成为可能。
图像由像素与颜色空间构成的,一方面在实验室离线分析的状态来识别烟叶是否透光,可以根据实验室开灯与不开灯的状态下烟叶自身的色彩空间的变化来完成,另一方面在实际的生产线可以根据烟叶颜色的自身特征来识别烟叶是否为水渍烟,根据所含有的特征区域占整个区域的比例来完成现场的水渍烟叶的统计监控分析。
而对于如何准确并简单的从工业相机中提取的水渍烟叶的特征信息,来对烟叶是否为水渍烟叶进行识别,结合图像的拼接,并根据该组特征信息的面积占比来对水渍烟叶一段时间内的实际比例进行测算,现有技术中并没有关于此种的技术方案记载。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种基于工业相机特征提取的润后水渍烟识别方法,以克服现有技术中的缺陷。
为了实现上述目的或者其他目的,本发明是通过以下技术方案实现的。
一种基于工业相机特征提取的润后水渍烟叶识别方法,其特征在于,所述识别方法包括如下步骤:
采集在线烟叶样本的所有原始烟叶图像,将烟叶样本人工分类为正常烟叶、水渍烟叶、近正常烟叶和近水渍烟叶,并获取所有原始烟叶图像的信息;
对所述原始烟叶图像的信息进行处理获得所有原始烟叶图像的特征向量;
对所述特征向量进行SNV预处理以消除异常数值信息;
对预处理后的特征采取有监督模式的二维投影处理获得二维投影向量;
将所述二维投影向量按照同种分类进行同列拼接,不同种分类之间上下拼接的原则形成新图像信息矩阵;
根据新图像信息矩阵和阈值识别新图像信息矩阵中包含水渍烟叶特征的连通区域面积;
计算包含水渍烟叶特征的连通区域面积占新图像信息矩阵的总面积的比例即为水渍烟的比例。
可选地,步骤1)中,采用工业相机采集在线烟叶样本的所有原始图像。
可选地,步骤1)中,将烟叶样本人工分类的依据为透光率。
本发明中通过实验室相机获取开灯和不开灯条件下水渍烟叶和正常烟叶的图像,并对图像进行信号分析,计算图像的透光率。以此先确认了水渍烟叶和正常烟叶在透光率上有很大差别,可以通过透光率将水渍烟叶和正常烟叶区分开来。
更为优选地,人工分类的标准为:水渍烟叶是指烟叶的叶面和叶背是一种颜色,对着阳光看通透性较强;近水渍烟叶是指具备水渍烟叶透光性的特征但是局部偏黄且透光不明显;近正常烟叶是指具备正常烟叶的特征但是透光性略强的;正常烟叶是指正常经过一润加工后的烟叶。
可选地,步骤2)中,获得所有原始烟叶图像的特征的方法包括如下步骤:
21)将所有原始烟叶图像的信息进行灰度处理;
22)通过PCA分解获得原始烟叶图像的特征向量即为所有原始烟叶图像的特征。
更为可选地,步骤21)中灰度处理获得的为二维灰度矩阵。
更为可选地,步骤22)中,PCA分解是指主成分分解,是指将二维灰度矩阵进行分解形成得分矩阵,选取得分矩阵的第一列作为特征向量,将各特征向量合并形成特征向量矩阵。
更优选地,具体地处理步骤为:
对每张原始烟叶图像由RGB空间进行灰度处理,得到二维灰度矩阵G;
对每个分类的每张图像的灰度矩阵Gi采取如下的方式进行特征抽提:
对Gi进行特征矩阵的分解获得:Gi=Pi*Si+Ei;
其中Pi为载荷矩阵,Si为得分矩阵;Ei是指分解的残差矩阵;
取得分矩阵Si的第一列形成列向量作为该图像的特征向量,记为S(:,1);
分别把各分类中的每个图像的特征向量合并成特征向量矩阵TS=[Si(:,1)]。
可选地,步骤3)中,SNV预处理为标准正态变量变换预处理。
更为优选地,SNV预处理是以特征向量矩阵中的每列数据为单位,将各列中各点值减去此列的平均值,然后除以此列的标准方差获得。
可选地,步骤4)中,有监督模式的二维投影是指线性判别LDA或偏最小二乘的线性判别。上述所述投影实质上就是从不同的角度去观察数据,寻找最能充分挖掘数据特征的作为最优投影方向。
可选地,步骤6)中,所述阈值为二维投影向量的阈值。
更为优选地,步骤6)中,所述阈值为1.14。
优选地,通过工业相机获得在生产线上一润工序后的近正常烟叶、近水渍烟叶、正常烟叶和水渍烟叶的图像;将图像由RGB空间转换成基于颜色模型的一维颜色,并获得颜色值,通过颜色值确认一维颜色值可以将水渍烟叶与其他烟叶区分开来。进一步地,将图像信息进行处理以获得近正常烟叶、近水渍烟叶、正常烟叶和水渍烟叶的特征向量;对特征向量进行预处理消除异常的数值信息。进一步地,对预处理之后的特征向量,采用有监督模式的二维投影处理获得二维投影向量。进一步地,把二维投影向量从上到下进行拼接,形成新图像信息;根据新图像信息计算水渍烟叶占样本烟叶中的比例。
本发明提供一种快速从图像信息中提取有效的水渍烟叶的特征信息,为水渍烟叶的分类以及识别指导提供依据。
本发明的积极效果:
快速的解决了一润后水渍烟叶的快速识别以及比例的计算,为通过工艺的控制来减少水渍烟叶的比例打下了坚实的基础,积累了原始的水渍烟叶识别数据;
本专利是图像特征抽提的经典范例,首次把化学计量学的方法与图像处理技术融合在一块,为运用化学计量学技术来推动图像处理技术的发展在图像特征抽提以及图像相似度的比对方面做了较为新颖的尝试。
附图说明
图1为正常烟叶的透光率T1,T1均值为29.34%,且T1<40%。
图2为水渍烟叶的透光率T2,T2均值为90%,且T2>85%。
图3为近正常烟叶的图像及编号。
图4为近水渍烟叶的图像及编号。
图5为正常烟叶的图像及编号。
图6为水渍烟叶的图像及编号。
图7为不同烟叶的类型对应的颜色值。
图8为不同种类的烟叶的特征图谱。
图9为识别结果。
图10为IS图像。
图11为水渍烟特征面积图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效。
在进一步描述本发明具体实施方式之前,应理解,本发明的保护范围不局限于下述特定的具体实施方案;还应当理解,本发明实施例中使用的术语是为了描述特定的具体实施方案,而不是为了限制本发明的保护范围。下列实施例中未注明具体条件的试验方法,通常按照常规条件,或者按照各制造商所建议的条件。
当实施例给出数值范围时,应理解,除非本发明另有说明,每个数值范围的两个端点以及两个端点之间任何一个数值均可选用。除非另外定义,本发明中使用的所有技术和科学术语与本技术领域技术人员通常理解的意义相同。除实施例中使用的具体方法、设备、材料外,根据本技术领域的技术人员对现有技术的掌握及本发明的记载,还可以使用与本发明实施例中所述的方法、设备、材料相似或等同的现有技术的任何方法、设备和材料来实现本发明。
在烟草行业中打叶复烤工艺流程主要包括以下步骤:真空回潮→一润→电子称→二润→打叶风分等步骤。本发明实施例中所述的水渍烟叶是指在打叶复烤工艺中一润处取得的烟叶。
在润叶工序受高温高湿的影响,烟叶容易产生水渍烟及蒸片,在复烤环节很难把水渍烟及蒸片调控到合适的水分,造成水渍片水分偏大,在片烟装箱过程中,水渍片温度及水分都较高,造成局部烟叶水分及温度过高,在保管过程中容易产生霉变及碳化,造成损失。
本发明实施例中识别方法识别的对象均为一润工序的烟叶。但不排除将此种识别方法及原理使用在其他工序中产生的水渍烟叶的识别和计算中。
实施例1
使用A复烤厂一润后烟叶作为样本烟叶,利用实验室相机分别获取开灯或不开灯条件下的正常烟叶和水渍烟叶的图像;对图像进行信号分析,用矩阵来表示数字图像,根据矩阵计算透光率。
本实施例中实验室相机的像素为2848×4288;但不限于使用此。
透光率T的计算公式为:
T为透光率;
R1为不开灯条件下图像对应的灰度矩阵;
R2为开灯条件下图像对应的灰度矩阵。
average具体指对按照上式计算获得的矩阵求平均值;
获得T为向量。
本实施例中,获得的透光率如图1和图2所示。其中图1为正常烟叶的透光率T1,T1均值为29.34%,且T1<40%。图2为水渍烟叶的透光率T2,T2均值为90%,且T2>85%。
由实施例可以看出仅通过透光率就可以将正常烟叶和水渍烟叶区分开来。
实施例2
本实施例为本申请中基于工业相机特征提取的润后水渍烟叶识别方法。
(1)在A复烤厂一润工序后用工业相机采集获取所有原始烟叶图像,将烟叶样本人工分类为正常烟叶、水渍烟叶、近正常烟叶和近水渍烟叶,并获取所有原始烟叶图像的信息。
人工分类的标准为:水渍烟叶是指烟叶的叶面和叶背是一种颜色,对着阳光看通透性较强;近水渍烟叶是指具备水渍烟叶透光性的特征但是局部偏黄且透光不明显;近正常烟叶是指具备正常烟叶的特征但是透光性略强的;正常烟叶是指正常经过一润加工后的烟叶。
进一步地,分别对四种分类的烟叶图像进行编号。其中:近正常烟叶的图像及编号分别为1、2、3、4、5、6、7、8、9、10、11、12、13和14,见图3;近水渍烟叶的图像及编号分别为1、2、3、4、5、6、7、8、9和10,见图4;正常烟叶的图像及编号为1、2、3、 4、5、6、7、8、9、10、11、12、13和14,见图5;水渍烟叶的图像及编号为1、2、3、4、 5、6、7、8、9、10、11、12和13,见图6。
(2)对所述原始烟叶图像的信息进行处理获得所有原始烟叶图像的特征,具体为:
21)将所有原始烟叶图像的信息进行灰度处理;具体地,灰度处理获得的为二维灰度矩阵;
22)通过PCA分解获得原始烟叶图像的特征向量即为所有原始烟叶图像的特征。
PCA分解是指主成分分解,是指将二维灰度矩阵进行分解形成得分矩阵,选取得分矩阵的第一列作为特征向量,将各特征向量合并形成特征向量矩阵。
PCA分解具体的处理步骤为:
对每张原始烟叶图像由RGB空间进行灰度处理,得到二维灰度矩阵G;
对每个分类的每张图像的灰度矩阵Gi采取如下的方式进行特征抽提:
对Gi进行特征矩阵的分解获得:Gi=Pi*Si+Ei;
其中Pi为载荷矩阵,Si为得分矩阵;Ei是指分解的残差矩阵;
取得分矩阵Si的第一列形成列向量作为该图像的特征向量,记为S(:,1);
分别把各分类中的每个图像的特征向量合并成特征向量矩阵TS=[Si(:,1)]。
根据TS作图如图8所示。因为TS的每一行是根据S(:,1)组成的,每一行的信息代表着每一幅图像降维后的信息,对矩阵TS的每一行画折线图就可以得到图8。
另外,将上述获取的原始烟叶图像由RGB空间转换成基于颜色模型的一维颜色,以进行烟叶颜色的特征抽提,获得颜色值。本发明中将烟叶图像由RGB空间转换成基于颜色模型的一维颜色的方法参见申请号为2015101338058的专利。
各样本的具体颜色值结果如下表表1所示:
表1
近正常烟叶 | 近水渍烟 | 正常烟叶 | 水渍烟叶 | |
1 | 1.451 | 1.311 | 0.968 | 0.752 |
2 | 1.557 | 1.263 | 1.632 | 1.135 |
3 | 1.299 | 1.006 | 1.447 | 1.137 |
4 | 1.536 | 1.136 | 1.523 | 1.034 |
5 | 1.554 | 1.182 | 1.536 | 0.164 |
6 | 1.427 | 1.429 | 1.594 | 0.247 |
7 | 1.509 | 1.128 | 1.5 | 1.534 |
8 | 1.418 | 1.089 | 1.53 | 0.539 |
9 | 1.47 | 1.443 | 1.46 | 0.497 |
10 | 1.414 | 1.432 | 1.374 | 0.549 |
11 | 1.556 | 1.53 | 0.935 | |
12 | 1.418 | 1.426 | 0.833 | |
13 | 1.53 | 1.199 | 1.285 | |
14 | 1.563 | 1.634 |
从表1中可以看出:水渍烟叶中的颜色值大部分偏小,正常烟叶中的颜色值大部分偏大,根据一维颜色值就可以对水渍烟叶与其他烟叶区分开来。如图7所示。图7是根据表1的数据按照近正常烟叶、近水渍烟、正常烟叶和水渍烟叶中各颜色值点依次直接划曲线得到的。
(3)对所述特征进行SNV预处理以消除异常数值信息;
所述SNV为标准正态变量变换,所述SNV预处理为:以特征向量矩阵中的每列数据单位,将各列中各点值减去此列的平均值,然后除以此列的标准方差。
SNV预处理是以特征向量矩阵中的每列数据为单位,将各列中各点值减去此列的平均值,然后除以此列的标准方差获得。
其具体的计算方式为:
其中,
m代表TS的行数,即总烟草样本个数;
TSi,k是指特征向量矩阵TS的第i行,第k列的元素;
TSi是TSi,k相应行的均值。
TSi,SNV是预处理之后的特征向量矩阵即为预处理后的特征。
(4)对预处理后的特征采取有监督模式的二维投影处理获得二维投影向量;具体为:
有监督模式的二维投影是指线性判别LDA或偏最小二乘的线性判别。上述所述投影实质上就是从不同的角度去观察数据,寻找最能充分挖掘数据特征的作为最优投影方向。
所谓投影实质上就是从不同的角度去观察数据,寻找最能充分挖掘数据特征的作为最优投影方向。
可在单位超球面中随机抽取若干个初始投影方向a(a1,a2,Λ,an),b(b1,b2,Λ,bn),计算其投影指标的大小,根据指标选大的原则,最后确定最大指标对应的解为最优投影方向。
若(a1,a2,Λ,an)为n维的单位向量,n为TS的列数,则样本i在一维线性空间的投影特征值zi的表达为
其中,aj为预处理后的图像的第一维投影向量的第j个值;
TSsnv(i,j)为通过PCA降维获取的特征矩阵,经过SNV预处理的得到的数值矩阵的第i行第j列的数值;
bj为预处理后的图像的第二维投影向量的第j个值;
zi是一行两列的数据,i=1:m,m是指样本的个数,对每个zi画散点图得到图9。
二维投影向量的获得原则是:在类间的距离与类内的距离的比值最大。
类间的距离R计算方式为:
ri,k=||zi-zk||(i,k=1,2,L,m) (4)
zk为二维投影值的第k个样本的投影值,此处的zk
为一个二维值;
ri,k为二维投影值的第k个样本与第i个样本的距离;
R表示的是所有TSsnv矩阵投影后投影值的距离和;
类内的距离S的计算方式为:
S(a)代表的是TSsnv矩阵投影后投影值的标准偏差
是指zi的均值。
S=s(a);
Q(zi)=R/S
Q(zi)代表的是,类间的距离与类内的距离比值。
选择目标函数Q(zi)取的最大值Max Q(zj)时对应的投影方向即为最优投影方向。
在最优投影方向上对应的二维投影向量即为采取有监督模式的二维投影处理获得二维投影向量。
(5)将所述二维投影向量按照同种分类进行同列拼接,不同种分类之间上下拼接的原则形成新图像信息矩阵。
把所有的样本的图像信息矩阵从上到下进行拼接;形成一个新图像矩阵信息,即为一个新的图像IS,具体结果如图10所示。图10是不同样本矩阵拼接在一块的图像,其目的在于对每一幅图像的局部能够计算,从图10可以看出,不同烟叶类型的样本,在图像上有明显的差异,水渍烟叶的颜色偏深。
(6)根据新图像信息矩阵和阈值识别新图像信息矩阵中包含水渍烟叶特征的连通区域面积。所述阈值为二维投影向量的阈值。对新图像计算每张图像包含水渍烟叶特征C1< PC(阈值)区域的面积总和;或者先在新图像上,按照像素点为单元,依据新图像的投影向量z1,z2计算新图像像素点,每块区域的类别属性,计算水渍烟叶类别区域的面积比例。
(7)计算包含水渍烟叶特征的连通区域面积占新图像信息矩阵的总面积的比例即为水渍烟的比例
具体地,本实施例中,根据新图像矩阵信息及新的图像IS,并基于图像的颜色特征,计算所有小于阈值1.14的数值信息的个数占整个新图像矩阵信息中数值个数的比值,即为水渍烟叶占样本烟叶中的比例。本实施例中为:3702495/12533760=29.54%。
对拼接的图像,整体进行颜色特征与图像特征S(i,1)的判断,根据经验选取阈值为p (此处的p=1.14);在整个拼接的图像上寻找颜色值与图像特征值s(i,1)均小于阈值p的区域,并计算每一区域的面积(此处的面积是指的像素和),画出在颜色值与图像特征值s(i, 1)均小于阈值p的区域的面积分布得到图11。
进一步地,根据上述技术方案原理及公开实施例,如果所有的烟叶是在不同的时间取的,那么这段时间内的拼接图像的特征区域面积比即为该段时间内水渍烟叶的比例。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (10)
1.一种基于工业相机特征提取的润后水渍烟叶识别方法,其特征在于:所述方法包括如下步骤:
采集在线烟叶样本的所有原始烟叶图像,将烟叶样本人工分类为正常烟叶、水渍烟叶、近正常烟叶和近水渍烟叶,并获取所有原始烟叶图像的信息;
对所述原始烟叶图像的信息进行处理获得所有原始烟叶图像的特征向量;
对所述特征向量进行SNV预处理以消除异常数值信息;
对预处理后的特征采取有监督模式的二维投影处理获得二维投影向量;
将所述二维投影向量按照同种分类进行同列拼接,不同种分类之间上下拼接的原则形成新图像信息矩阵;
根据新图像信息矩阵和阈值识别新图像信息矩阵中包含水渍烟叶特征的连通区域面积;
计算包含水渍烟叶特征的连通区域面积占新图像信息矩阵的总面积的比例即为水渍烟的比例。
2.如权利要求1所述方法,其特征在于:采用工业相机采集在线烟叶样本的所有原始图像。
3.如权利要求1所述方法,其特征在于:获得所有原始烟叶图像的特征的方法包括如下步骤:
21)将所有原始烟叶图像的信息进行灰度处理;
22)通过PCA分解获得原始烟叶图像的特征向量即为所有原始烟叶图像的特征。
4.如权利要求1所述方法,其特征在于:获得所有原始烟叶图像的特征的方法包括如下步骤:
对每张原始烟叶图像由RGB空间进行灰度处理,得到二维灰度矩阵G;
对每个分类的每张图像的灰度矩阵Gi采取如下的方式进行特征抽提:
对Gi进行特征矩阵的分解获得:Gi=Pi*Si+Ei;
其中Pi为载荷矩阵,Si为得分矩阵;Ei是指分解的残差矩阵;
取得分矩阵Si的第一列形成列向量作为该图像的特征向量,记为S(:,1);
分别把各分类中的每个图像的特征向量合并成特征向量矩阵TS=[Si(:,1)]。
5.如权利要求4所述方法,其特征在于:SNV预处理为标准正态变量变换预处理,是以特征向量矩阵中的每列数据为单位,将各列中各点值减去此列的平均值,然后除以此列的标准方差获得。
6.如权利要求1所述方法,其特征在于:有监督模式的二维投影是指线性判别LDA或偏最小二乘的线性判别。
7.如权利要求6所述方法,其特征在于:二维投影的原则是类间的距离与类内的距离的比值最大。
8.如权利要求7所述方法,其特征在于:类间的距离R的计算方式为:
ri,k=||zi-zk||(i,k=1,2,L,m)
zi为第i个样本在二维投影空间的投影特征值;
zk为第k个样本在二维投影空间的投影特征值;
m为所有烟叶样本的数量;
ri,k为在二维投影空间,第i个样本与第k个样本的投影值距离;
R表示的是所有TSsnv矩阵投影后投影值的距离和;
所述TSsnv为通过PCA降维获取的特征矩阵,经过SNV预处理得到的数据矩阵。
9.如权利要求7所述方法,其特征在于:类内的距离S的计算方式为:
S(a)代表的是TSsnv矩阵投影后投影值的标准偏差
是指zi的均值;
n是指在二维投影空间投影样本的数量;此处的n=m;
S=s(a)。
10.如权利要求1所述方法,其特征在于:所述阈值为二维投影向量的阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610195473.0A CN105913056B (zh) | 2016-03-29 | 2016-03-29 | 一种基于工业相机特征提取的润后水渍烟叶识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610195473.0A CN105913056B (zh) | 2016-03-29 | 2016-03-29 | 一种基于工业相机特征提取的润后水渍烟叶识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105913056A CN105913056A (zh) | 2016-08-31 |
CN105913056B true CN105913056B (zh) | 2019-03-29 |
Family
ID=56745379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610195473.0A Active CN105913056B (zh) | 2016-03-29 | 2016-03-29 | 一种基于工业相机特征提取的润后水渍烟叶识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105913056B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268895A (zh) * | 2018-01-12 | 2018-07-10 | 上海烟草集团有限责任公司 | 基于机器视觉的烟叶部位识别方法、电子设备及存储介质 |
CN111947927B (zh) * | 2020-07-16 | 2021-09-07 | 南京航空航天大学 | 一种基于色度理论的滚动轴承故障检测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200931323A (en) * | 2008-01-10 | 2009-07-16 | Int Currency Tech | License identification device |
CN104268834B (zh) * | 2014-09-19 | 2017-09-08 | 贵州大学 | 一种基于数字图像处理技术的标准等级烟叶图像数据库的合成方法 |
CN104331702A (zh) * | 2014-11-03 | 2015-02-04 | 黄辉 | 基于图像识别的茶叶鲜叶分叉数量识别方法 |
CN104713835B (zh) * | 2015-03-25 | 2017-09-29 | 上海烟草集团有限责任公司 | 一种烟叶颜色在线数值化识别方法 |
CN104914105A (zh) * | 2015-06-09 | 2015-09-16 | 中国农业科学院茶叶研究所 | 一种基于图像识别技术的茶叶等级鉴定方法 |
-
2016
- 2016-03-29 CN CN201610195473.0A patent/CN105913056B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105913056A (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110569837B (zh) | 优化损伤检测结果的方法及装置 | |
Weyler et al. | Joint plant instance detection and leaf count estimation for in-field plant phenotyping | |
Barbedo | A novel algorithm for semi-automatic segmentation of plant leaf disease symptoms using digital image processing | |
CN110705655A (zh) | 一种基于光谱和机器视觉耦合的烟叶分类方法 | |
CN110472489B (zh) | 一种人脸美丽等级预测方法、装置及存储介质 | |
CN104463199A (zh) | 基于多特征及分割再校正的矿岩块度分类方法 | |
CN117115147A (zh) | 一种基于机器视觉的纺织品检测方法及系统 | |
CN106228150A (zh) | 基于视频图像的烟雾检测方法 | |
CN110261329A (zh) | 一种基于全谱段高光谱遥感数据的矿物识别方法 | |
CN105023027A (zh) | 基于多次反馈机制的现场鞋底痕迹花纹图像检索方法 | |
CN104713835A (zh) | 一种烟叶颜色在线数值化识别方法 | |
CN105913056B (zh) | 一种基于工业相机特征提取的润后水渍烟叶识别方法 | |
CN114723994A (zh) | 一种基于双分类器对抗增强网络的高光谱图像分类方法 | |
CN109035223A (zh) | 一种用于卫星遥感图像可用度的智能评估方法 | |
CN111178405A (zh) | 一种融合多种神经网络的相似物体识别方法 | |
Calixto et al. | Development of a computer vision approach as a useful tool to assist producers in harvesting yellow melon in northeastern Brazil | |
Su | [Retracted] Data Research on Tobacco Leaf Image Collection Based on Computer Vision Sensor | |
CN112528726B (zh) | 一种基于光谱成像和深度学习的棉蚜虫害监测方法及系统 | |
Pei et al. | State recognition and temperature rise time prediction of tobacco curing using multi-sensor data-fusion method based on feature impact factor | |
Li et al. | Image processing for crop/weed discrimination in fields with high weed pressure | |
Buyukarikan et al. | Classification of physiological disorders in apples using deep convolutional neural network under different lighting conditions | |
CN114549494A (zh) | 一种温室生产环境下草莓白粉病快速检测方法 | |
CN113283378A (zh) | 一种基于梯形区域归一化像素差异特征的猪脸检测方法 | |
Joycy et al. | Notice of Violation of IEEE Publication Principles: Survey on automatic segmentation of relevant textures in agricultural images | |
CN111626409B (zh) | 一种图像质量检测的数据生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Xu Qimin Inventor after: Dou Jiayu Inventor after: Wang Jizhou Inventor after: Wei Jun Inventor after: Chen Qing Inventor after: Zhang Jun Inventor before: Xu Qimin Inventor before: Dou Jiayu Inventor before: Wang Jizhou Inventor before: Chen Qing Inventor before: Zhang Jun |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |