CN102495068B - 大米垩白米在线检测方法 - Google Patents
大米垩白米在线检测方法 Download PDFInfo
- Publication number
- CN102495068B CN102495068B CN201110403152.2A CN201110403152A CN102495068B CN 102495068 B CN102495068 B CN 102495068B CN 201110403152 A CN201110403152 A CN 201110403152A CN 102495068 B CN102495068 B CN 102495068B
- Authority
- CN
- China
- Prior art keywords
- rice
- chalk
- theta
- image data
- chalkiness
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明实施例提供了一种大米垩白米在线检测方法和***,用于实时、准确检测大米垩白率和垩白度,获取大米产线生产信息。本发明实施例方法包括:周期性的采集数目相同的米粒作为检测样本;采用摄像装置对所述检测样品进行样本图像采集,形成采集图像数据;将所述采集图像数据进行灰度图像处理,形成灰度图像数据;将所述灰度图像数据分出正常区域、反光区域和垩白区域;通过期望最大值算法对垩白像素点进行进一步检测;根据垩白像素点检测出垩白米的垩白率和垩白度。
Description
技术领域
本发明涉及图像数据处理领域,尤其涉及一种大米垩白米在线检测方法和***。
背景技术
我国是世界上最大的稻米生产国。而稻米品质问题已经成为制约我国大米生产、销售和出口的瓶颈。大米的垩白度是大米品质的一项十分重要的指标,也是影响大米市场价格的重要因素。优质稻米国家标准GB/T17981-1999中定义:“垩白”是指米粒胚乳中的白色不透明部分,根据垩白在大米胚乳中所处的位置分为腹白、心白和背白。在优质稻米检测指标中,垩白指标有两个:垩白度和垩白粒率,国家规定的检测方法如下:
垩白粒率:从优质稻米精米试样中随机取整精米数100粒,拣出有垩白的米粒,按下式求出垩白粒率:
垩白度:在检测垩白粒率时拣出的垩白米中,随机取10粒,不足10粒按实有数取,将米粒平放,正视观察,逐粒目测垩白面积占整个籽粒投影面积的百分率,求出平均值,垩白度按下式计算:
垩白度(%)=垩白粒率×垩白大小
垩白米蒸煮以后,饭粒上会出现裂纹,米饭蓬松中空,食味较差,商品性差,市场价格低,因此有必要在产线生产时对垩白米进行检测。
发明内容
本发明实施例提供了一种大米垩白米在线检测方法和***,用于实时、准确检测大米垩白率和垩白度,获取大米产线生产信息。
本发明一实施例提供的一种大垩白米在线检测方法,用于在线检测大米垩白率和垩白度,判断大米生产状态,包括以下步骤:
周期性的采集数目相同的米粒作为检测样本;
采用摄像装置对所述检测样品进行样本图像采集,形成采集图像数据;
将所述采集图像数据进行灰度图像处理,形成灰度图像数据;
将所述灰度图像数据分出正常区域、反光区域和垩白区域;
通过期望最大值算法对垩白像素点进行进一步检测,
根据垩白像素点检测出垩白米的垩白粒和垩白度。
进一步地,通过期望最大值算法对垩白像素点进行进一步的准确检测,具体为:
A、初始设定图像的像素值均匀分布,取米粒区域的灰度和色度的最小值作为正常区域均值,平均值作为反光区域均值,最大值作为垩白区域均值,假设其出现的初始概率值分别为0.4、0.3、0.3,均方差为128;
B、利用EM算法对参数模型进行迭代计算,
E步骤:计算完整数据的对数似然函数的期望,记为:
上式中,p(yi|xi,θ(t))表示第t次迭代之后的后验概率, 表示第t次迭代时第i类像素出现的概率,θ(t)表示第t次迭代时的参数集,xi表示样本数据(即灰度值),yi表示隐含的参数(即样本属于哪一类);
C、M步骤:将E步骤的p(yi|xi,θ(t))和θ(t)代入下列迭代公式进行计算:
D、用M步骤中得到的参数作为θ(t+1)代替θ(t),不断交替计算E步骤和M步骤,对于给定的ε>0,如果|θ(t)-θ(t+1)|<ε,则停止迭代。此时可以得到此4个高斯分布的一组参数θ(*);
E、利用参数集θ(*),计算每个像素分属于此3类的概率,即若p(yi=j|xi,θ(*))=max{p(yi=j|xi,θ(*))},j=1,2,3,则将xi归为j类,从而准确检测到垩白米粒的像素点。
优选地,采集图像数据的像素为800*600以上,帧率25fps以上。
本发明又一实施例提供了一种大米垩白米在线检测***,用于在线检测大米垩白率和垩白度,判断大米生产状态,包括:
采样装置,周期性的对在线大米进行采样收集,形成检测样本;
图像信号采集装置,对所述检测样本进行摄像以及图像数据采集,形成采集图像数据;
图像识别与处理装置,用于接收所述采集图像数据,对所述采集图像数据进行灰度化处理,形成灰度图像数据;将所述灰度图像数据分出正常区域、反光区域和垩白区域;通过期望最大值算法对垩白像素点进行进一步的准确检测,根据垩白像素点检测出垩白米的垩白粒和垩白度;
客户端装置,用以实时监测大米垩白率和垩白度,获取产线运行状态信息。
优选地,采样装置安装在大米生产线上,周期性对大米进行采样。
优选地,周期性对大米进行采样的时间为每20秒采样一次。
优选地,采样装置为筛状带有凹坑的采样器。
优选地,图像信号采集装置包括摄像装置和数据采集卡,
优选地,摄像装置具有800×600以上像素,采集检测样品图像数据;
优选地,数据采集卡包括多路数据输入选择模块、模拟/数字转换模块和通信接口模块,多路数据输入选择模块接收图像数据,并通过模拟/数字转换模块将模拟的图像数据转换为数字图像数据,再通过通信接口模块发送给图像识别与处理装置进行图像处理。
从以上技术方案可以看出,可以准确实时的对大米垩白米粒的垩白率和垩白度进行检测,有效控制大米生产质量。
附图说明
图1是本发明实施例一大米垩白米在线检测方法流程图;
图2是本发明实施例二大米垩白米在线检测***结构示意图。
具体实施方式
本发明实施例提供了一种大米碎米粒在线检测方法和***,以下将结合具体实施例详细描述本发明实现过程。
一种大米碎米粒在线检测方法,用于在线检测大米碎米粒,判断大米生产状态,流程图如图1所示:
S101,周期性的采集数目相同的米粒作为检测样本;
本实施例中采用每分钟3次即每20秒一次的周期对样本进行检测,每次采样的数目是100粒,本领域的技术人员可知,可以根据实际产线生产的需求控制检测样本的时间和样本的大小。
S102,采用摄像装置对检测样品进行样本图像采集,形成采集图像数据;
本实施例中摄像装置为24位真彩高速摄像装置,有效像素在800×600以上,每像素24bit,帧率为25fps以上,可保证采集图像数据的准确性,保证检测效果。
S103,将采集图像数据进行灰度图像处理,形成灰度图像数据;
将采集的图像数据通过接口接入数据处理装置进行数据处理,先将24位彩色图像转换为8位灰度图像,以便数据处理装置进行处理。
S104,将所述灰度图像数据分出正常区域、反光区域和垩白区域;
S105,通过期望最大值算法对垩白像素点进行进一步检测;
S106,根据垩白像素点检测出垩白米的垩白粒和垩白度。
通过期望最大值算法对垩白像素点进行进一步的准确检测,具体为:A、初始设定图像的像素值均匀分布,取米粒区域的灰度和色度的最小值作为正常区域均值,平均值作为反光区域均值,最大值作为垩白区域均值,假设其出现的初始概率值分别为0.4、0.3、0.3,均方差为128;
B、利用EM算法对参数模型进行迭代计算,
E步骤:计算完整数据的对数似然函数的期望,记为:
上式中,p(yi|xi,θ(t))表示第t次迭代之后的后验概率, 表示第t次迭代时第i类像素出现的概率,θ(t)表示第t次迭代时的参数集,xi表示样本数据(即灰度值),yi表示隐含的参数(即样本属于哪一类);
C、M步骤:将E步骤的p(yi|xi,θ(t))和θ(t)代入下列迭代公式进行计算:
D、用M步骤中得到的参数作为θ(t+1)代替θ(t),不断交替计算E步骤和M步骤,对于给定的ε>0,如果|θ(t)-θ(t+1)|<ε,则停止迭代。此时可以得到此4个高斯分布的一组参数θ(*);
E、利用参数集θ(*),计算每个像素分属于此3类的概率,即若p(yi=j|xi,θ(*))=max{p(yi=j|xi,θ(*))},j=1,2,3,则将xi归为j类,从而准确检测到垩白米粒的像素点。
一种大米碎米粒在线检测***,结构示意图如图2所示,包括:
采样装置21,周期性的对在线大米进行采样收集,形成检测样本;
本实施例用采用推进机构将大米样品从生产线推入直径60-100mm的采样筛进行采样,每隔20秒对大米进行一次采样,采样时采用筛状带有凹坑的采样器,且采样器每次固定采样大米100粒,采样器倾斜放置,凹坑中的大米很容易被产险上流动的大米带走,方便下一次采样。
图像信号采集装置22,对所述检测样本进行摄像以及图像数据采集,形成采集图像数据;
图像信号采集装置22包括摄像装置和数据采集卡,摄像装置具有800*600以上像素,每像素24bit,帧率为25fps以上,可保证采集图像数据的准确性,保证检测效果,采集检测样品图像数据;数据采集卡包括多路数据输入选择模块、模拟/数字转换模块和通信接口模块,多路数据输入选择模块接收图像 数据,并通过模拟/数字转换模块将模拟的图像数据转换为数字图像数据,再通过通信接口模块发送给所述图像识别与处理装置进行图像处理。
图像识别与处理装置23,用于接收所述采集图像数据,对所述采集图像数据进行灰度化,形成灰度图像数据;将所述灰度图像数据分出正常区域、反光区域和垩白区域;通过期望最大值算法对垩白像素点进行进一步的准确检测,根据垩白像素点检测出垩白米的垩白粒和垩白度;;
将采集的图像数据通过接口接入数据处理装置进行数据处理,先将24位彩色图像转换为8位灰度图像,以便数据处理装置进行处理。
客户端装置24,用以实时监测碎米粒数,获取产线运行状态信息。
客户端装置24接收图像识别与处理装置通过网络接口发送过来的碎米率百分比,实时检测碎米粒是否超标,如果超过产线定义的碎米率,则可将产线停止,重新进行前序生产工序的调解,以保证大米的生产质量。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上对本发明所提供的一种大米垩白米在线检测方法和***进行了详细介绍,对于本领域的一般技术人员,依据本发明实施例的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (1)
1.一种大米垩白米在线检测方法,用于在线检测大米垩白米粒,判断大米生产状态,其特征在于,包括以下步骤:
周期性的采集数目相同的米粒作为检测样本;
采用摄像装置对所述检测样品进行样本图像采集,形成采集图像数据;
将所述采集图像数据进行灰度图像处理,形成灰度图像数据;
将所述灰度图像数据分出正常区域、反光区域和垩白区域;
通过期望最大值算法对垩白像素点进行进一步检测;
根据垩白像素点检测出垩白米的垩白粒和垩白度;
其中,通过期望最大值算法对垩白像素点进行进一步的准确检测,具体为:
A、初始设定图像的像素值均匀分布,取米粒区域的灰度和色度的最小值作为正常区域均值,平均值作为反光区域均值,最大值作为垩白区域均值,假设其出现的初始概率值分别为0.4、0.3、0.3,均方差为128;
B、利用EM算法对参数模型进行迭代计算,
E步骤:计算完整数据的对数似然函数的期望,记为:
上式中,p(yi|xi,θ(t))表示第t次迭代之后的后验概率,表示第t次迭代时第i类像素出现的概率,θ(t)表示第t次迭代时的参数集,xi表示样本数据(即灰度值),yi表示隐含的参数(即样本属于哪一类);
C、M步骤:将E步骤的p(yi|xi,θ(t))和θ(t)代入下列迭代公式进行计算:
D、用M步骤中得到的参数作为θ(t+1)代替θ(t),不断交替计算E步骤和M步骤,对于给定的ε>0,如果|θ(t)-θ(t+1)|<ε,则停止迭代。此时可以得到此4个高斯分布的一组参数θ(*);
E、利用参数集θ(*),计算每个像素分属于此3类的概率,即若p(yi=j|xi,θ(*))=max{p(yi=j|xi,θ(*))},j=1,2,3,则将xi归为j类,从而准确检测到垩白米粒的像素点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110403152.2A CN102495068B (zh) | 2011-12-01 | 2011-12-01 | 大米垩白米在线检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110403152.2A CN102495068B (zh) | 2011-12-01 | 2011-12-01 | 大米垩白米在线检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102495068A CN102495068A (zh) | 2012-06-13 |
CN102495068B true CN102495068B (zh) | 2014-01-08 |
Family
ID=46186909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110403152.2A Expired - Fee Related CN102495068B (zh) | 2011-12-01 | 2011-12-01 | 大米垩白米在线检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102495068B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103411929B (zh) * | 2013-05-10 | 2015-09-02 | 中国水稻研究所 | 一种稻米透明度检测方法 |
CN103674816B (zh) * | 2013-12-25 | 2016-05-11 | 浙江工商大学 | 垩白米粒检测方法及*** |
CN105319223B (zh) * | 2015-12-07 | 2018-03-20 | 浙江工商大学 | 一种垩白米的检测方法和检测*** |
CN106872471A (zh) * | 2017-03-27 | 2017-06-20 | 中国水稻研究所 | 稻米垩白检测方法 |
CN111435541A (zh) * | 2019-01-14 | 2020-07-21 | 珠海格力电器股份有限公司 | 获取米粒的垩白度的方法、装置和烹饪器具 |
CN109781730B (zh) * | 2019-01-28 | 2020-06-09 | 江南大学 | 一种快速识别垩白米的方法 |
CN110009609B (zh) * | 2019-03-26 | 2021-03-30 | 江南大学 | 一种快速检测黄粒米的方法 |
CN110954551B (zh) * | 2019-12-05 | 2023-02-10 | 浙江五芳斋实业股份有限公司 | 一种用于青团品质的在线检测方法 |
CN113776993A (zh) * | 2021-07-28 | 2021-12-10 | 深圳市麦稻智联科技有限公司 | 一种大米在线检测***、检测方法及检测设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1845949B1 (en) * | 2005-02-03 | 2008-05-21 | Nycomed Pharma AS | Melt granulation of a composition containing a calcium-containing compound |
CN101551376A (zh) * | 2008-03-31 | 2009-10-07 | 国家粮食储备局成都粮食储藏科学研究所 | 一种米类加工精度的检验方法 |
CN202471601U (zh) * | 2011-12-06 | 2012-10-03 | 李同强 | 大米垩白米在线检测*** |
-
2011
- 2011-12-01 CN CN201110403152.2A patent/CN102495068B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1845949B1 (en) * | 2005-02-03 | 2008-05-21 | Nycomed Pharma AS | Melt granulation of a composition containing a calcium-containing compound |
CN101551376A (zh) * | 2008-03-31 | 2009-10-07 | 国家粮食储备局成都粮食储藏科学研究所 | 一种米类加工精度的检验方法 |
CN202471601U (zh) * | 2011-12-06 | 2012-10-03 | 李同强 | 大米垩白米在线检测*** |
Non-Patent Citations (4)
Title |
---|
《基于灰度图像的大米垩白检测算法研究》;孙翠霞等;《广西工学院学报》;CNKI;20100630(第1期);37,38页 * |
基于高光谱图像技术的大米垩白检测;黎静等;《南方农机》;20090220(第01期);27-29 * |
孙翠霞等.《基于灰度图像的大米垩白检测算法研究》.《广西工学院学报》.CNKI,2010,(第1期),37-38页. |
黎静等.基于高光谱图像技术的大米垩白检测.《南方农机》.2009,(第01期),27-29. |
Also Published As
Publication number | Publication date |
---|---|
CN102495068A (zh) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102495068B (zh) | 大米垩白米在线检测方法 | |
CN202421024U (zh) | 大米碎米率在线检测*** | |
CN108509976A (zh) | 动物的识别装置和方法 | |
CN102494977B (zh) | 大米碎米率在线检测方法和*** | |
CN113408328B (zh) | 基于毫米波雷达的手势分割与识别算法 | |
CN105608446A (zh) | 一种视频流异常事件的检测方法及装置 | |
CN106874929B (zh) | 一种基于深度学习的珍珠分类方法 | |
CN109410190B (zh) | 基于高分辨率遥感卫星影像的杆塔倒断检测模型训练方法 | |
CN109377487B (zh) | 一种基于深度学习分割的水果表面缺陷检测方法 | |
CN109376589A (zh) | 基于卷积核筛选ssd网络的rov形变目标与小目标识别方法 | |
CN106772379B (zh) | 基于傅里叶变换的双参数cfar舰船检测方法 | |
CN108960142B (zh) | 基于全局特征损失函数的行人再识别方法 | |
Sun et al. | Wheat head counting in the wild by an augmented feature pyramid networks-based convolutional neural network | |
CN109086794A (zh) | 一种基于t-lda主题模型的驾驶行为模式识方法 | |
CN103927516A (zh) | 基于数字图像处理的海水珍珠鉴别*** | |
CN110110610A (zh) | 一种用于短视频的事件检测方法 | |
CN113096178A (zh) | 猪只体重的估测方法、装置、设备和存储介质 | |
CN202471601U (zh) | 大米垩白米在线检测*** | |
CN111476307A (zh) | 一种基于深度领域适应的锂电池表面缺陷检测方法 | |
CN106407966A (zh) | 一种应用于考勤的人脸识别方法 | |
CN115496891A (zh) | 一种小麦倒伏程度分级方法和装置 | |
CN109919215A (zh) | 基于聚类算法改进特征金字塔网络的目标检测方法 | |
CN108634925A (zh) | 一种基于手势识别的视力测试*** | |
CN108717522A (zh) | 一种基于深度学习和相关滤波的人体目标跟踪方法 | |
CN102521874B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140108 Termination date: 20141201 |
|
EXPY | Termination of patent right or utility model |