一种基于红外图像的禽蛋品质鉴别方法与流程

本发明属于图像处理领域，具体涉及一种利用图像处理技术鉴别禽蛋品质的方法。

背景技术：

禽蛋产业不仅是我国农业的支柱产业，也是关系国计民生和社会稳定的重要产业。目前国内外已有的禽蛋品质鉴别方法有：通过人的感官直接鉴别的方法，即通过肉眼观测禽蛋表面是否清洁完整；还有采用物理指标衡量禽蛋品质的方法，例如测量禽蛋尺寸及密度并与国际标准值比较；针对表面品质的无损鉴别方法，即通过采用多个ccd摄像头拍摄禽蛋，并对采集到的数据进行综合处理；还有通过近红外光谱投射或红外光漫反射来观测禽蛋结构与内部化学成分的变化。

然而这些方法都或多或少存在着缺点。首先，通过人的感官指标来对禽蛋进行直接鉴别缺少科学严谨性，采用物理衡量指标方法的理论支撑不够充分，此外密度检测属于有损检测，不利于禽蛋的销售。对禽蛋表面的无损检测方法，只着眼于蛋壳表面情况包括裂纹、污斑，无法检测其内部是否发生变质，无法应用到生产线中。同时，红外光谱分析是利用不同物质对不同频率的红外光的不同吸收性来鉴别物质种类的，利用近红外光谱分析禽蛋品质也存在不能检测到禽蛋内部是否变质的问题。近红外主要表征的是禽蛋壳表面的性质，其内部蛋白质、脂质等有机物特性很难通过不透明的蛋壳反映到外部，故不能很好地判断禽蛋内部有机物的情况，无法有效判断禽蛋变质与否。同时，近红外光谱分析技术存在灵敏度低、抗干扰性差等缺点。

技术实现要素：

针对上述现有技术的不足，本发明提供一种基于红外图像的禽蛋品质鉴别方法，用于解决现有技术中的禽蛋鉴别方法不能准确、高效的对禽蛋内部品质进行无损鉴别的技术问题，突破了通过观测禽蛋表面性质无法深入内部品质的局限性，能够通过检测禽蛋表面性质快速准确的检测出禽蛋内部是否变质，能对禽蛋品质进行等级区分，对禽蛋不会产生损害和诱发变质。

为了解决上述技术问题，本发明采用了如下的技术方案：一种基于红外图像的禽蛋品质鉴别方法，包括以下步骤：

步骤1：利用红外热像仪采集禽蛋的红外图像并将红外图像输入给计算机；

步骤2：对红外图像进行预处理，然后进行二值化，得到二值图像；

步骤3：对二值图像的白色区域进行连通性检测，得到若干个白色连通区域，并根据禽蛋外形特征对白色连通区域进行约束，去除掉不符合要求的连通区域，得到禽蛋的位置分布图；

步骤4：对位置分布图进行降噪，得到表征禽蛋区域的连通区域图像，每一个连通区域对应一个禽蛋，从1到n对各连通区域进行编号；

步骤5：根据红外图像和连通区域图像，计算禽蛋分割图；

步骤6：根据禽蛋分割图计算每个禽蛋区域的灰度均值μi、均方差σi以及面积si，构成特征向量xi＝[μi,σi,si]^t，i＝{1,2,...n}；

步骤7：将特征向量xi带入采用训练样本预先训练好的判决模型，计算出模型值ni；

步骤8：将计算出的模型值ni与判决阈值进行比较，从而完成对禽蛋品质的鉴别。

优选的，训练样本中包含品质合格禽蛋、品质差禽蛋以及品质更差禽蛋；判决阈值包括nt、nr，其中，nt为训练样本中品质差禽蛋按照判决模型所计算出的模型值；nr为训练样本中品质更差禽蛋按照判决模型所计算出的模型值；当ni＞nt，则为品质良好禽蛋；当nt≥ni≥nr，则为品质较差禽蛋；当ni＜nr，则为变质禽蛋。

优选的，判决模型按照如下步骤训练：

步骤201：收集已知品质的禽蛋作为训练样本，令禽蛋的样本特征向量为变量x，按照所述步骤1～6计算出训练样本中各个禽蛋的样本特征向量；其中，训练样本中包含品质合格禽蛋、品质差禽蛋以及品质更差禽蛋；

步骤202：采用品质合格禽蛋训练一个三维空间的高斯模型，并计算模型参数以及∑，为样本特征向量均值，∑为样本特征向量协方差，

步骤203：将模型参数以及∑代入由品质合格禽蛋训练的高斯模型，得到判决模型如下：

步骤204：将品质差禽蛋以及品质更差禽蛋的特征向量代入判决模型，计算出若干模型值，选取品质差禽蛋模型值的平均值作为判决阈值nt，选取品质更差禽蛋模型值的平均值作为判决阈值nr。

优选的，步骤2中图像预处理按照如下步骤进行：首先，对红外图像采用3×3模板进行均值滤波降噪；然后对红外图像中各像素点的灰度值进行线性变换，将灰度值变换到0～255范围内，保证亮度和对比度的统一。

优选的，步骤3中按如下步骤去掉不符合要求的连通区域：

步骤301：根据各个连通区域的像素点的数量计算各个连通区域的面积s′；根据拍摄的实际情况设定参数c1与c2；其中，c1为红外图像中最小禽蛋区域的像素点的数量，c2为红外图像中最大禽蛋区域的像素点的数量；

步骤302：计算各个连通区域的外接矩形，设外接矩形的宽度为w，高度为h；

步骤303：将外接矩形的宽度w与参数c1与c2进行比较，去掉w＜c1和w＞c2的连通区域；

步骤304：根据禽蛋形状近似于椭圆的特点，去掉w/h＞1.5和h/w＞1.5的连通区域；

步骤305：根据禽蛋形状近似于椭圆的特点，去掉s′/(w*h)＜0.7的连通区域。

与现有技术相比，本发明的具有以下有益效果：

1、由于禽蛋内部因微生物侵入或滋生会产生生命代谢而释放热量，因此变质禽蛋表面的温度会比未变质禽蛋高，那么通过红外图像就能间接检测到禽蛋品质的变化情况：红外图像的灰度值反映了禽蛋表面的温度分布情况。这种间接检测方法，突破了观测禽蛋外表面性质无法深入内部品质检测的局限性，同时使得评价标准更加统一客观，符合工业上标准化的要求。

2、本发明的检测方法采用非接触式的检测，而且是采用被动式红外成像，不会对禽蛋造成物理损伤或者产生诱发其变质的不利因素。

3、本发明的检测方法不仅采用温度作为鉴别禽蛋品质的因素，还同时将禽蛋的大小作为品质因素，使得禽蛋品质的分级更加客观、科学。

4、本发明方便易行，大大提升了禽蛋品质检测速度，适合工业化生产场景下的应用。

5、本发明通过灰度值线性变换、滤波降噪以及根据禽蛋外形的约束去除不符合要求的连通区域，大大提高了检测的精度。

附图说明

图1是具体实施方式的流程图；

图2是禽蛋的红外图像；

图3是禽蛋的连通区域图像；

图4是禽蛋分割图；

图5禽蛋的特征向量的三维坐标图。

具体实施方式

下面结合附图和优选实施方式对本发明作进一步的详细说明。

首先对本发明的原理进行说明：禽蛋内部因为微生物侵入或滋生会有生命代谢，释放热量，造成禽蛋外表面温度升高；红外热像仪能够捕捉到不同品质禽蛋热辐射信号，得到禽蛋的红外图像，红外图像的灰度值反映了禽蛋表面温度的分布情况。对禽蛋的红外图像进行处理，实现图像中的禽蛋分割，对得到禽蛋图像信息进行分析计算，获取温度分布及面积特征，然后根据样本训练结构对禽蛋品质进行鉴别，得到禽蛋品质等级的客观评价。

如图1所示，一种基于红外图像的禽蛋品质鉴别方法，包括以下步骤：

步骤1：利用红外热像仪采集禽蛋的红外图像并将红外图像输入给计算机；采用红外热像仪，对禽蛋进行拍摄获得图像，将图像信号通过热像仪的数字接口传输到主控计算机，在主控计算机上编写软件实现网络协议、镜头对焦、数据收发及图像解码，获得清晰的红外图像，禽蛋的红外图像如图2所示。本具体实施方式中采用采用flir公司的photon160x120红外热像仪拍摄禽蛋获得禽蛋的红外图像，为了后续处理方便，拍摄时使各禽蛋留有一定间隔，使得在图像中容易区分。

步骤2：对红外图像进行预处理，然后进行二值化，得到二值图像；

步骤3：对二值图像的白色区域(灰度为255)进行连通性检测，如果两个白色的像素点在横向、纵向和斜向相邻，则认为它们连通，属于同一个连通区域，从而得到若干个白色连通区域。由于白色连通区域包括禽蛋区域，也包括干扰区域，需要根据禽蛋外形特征对白色连通区域进行约束，去除掉不符合要求的连通区域，得到禽蛋的位置分布图；

步骤4：对位置分布图进行降噪，得到表征禽蛋区域的连通区域图像，如图3所示，每一个连通区域对应一个禽蛋，从1到n对各连通区域进行编号；

步骤5：根据红外图像和连通区域图像，计算禽蛋分割图，图4为禽蛋分割图；

步骤6：根据禽蛋分割图计算每个禽蛋区域的灰度均值μi、均方差σi以及面积si，构成特征向量xi＝[μi,σi,si]^t，i＝{1,2,...n}；禽蛋的特征向量参数如下表所示：

禽蛋的特征向量xi在三维空间中可表示为一个点，如图5所示。

步骤7：将特征向量xi带入采用训练样本预先训练好的判决模型，计算出模型值ni；

步骤8：将计算出的模型值ni与阈值进行比较，从而完成对禽蛋品质的鉴别。

本具体实施方式中，训练样本中包含品质合格禽蛋、品质差禽蛋以及品质更差禽蛋；阈值包括nt、nr，其中，nt为训练样本中品质差禽蛋按照判决模型所计算出的模型值；nr为训练样本中品质更差禽蛋按照判决模型所计算出的模型值；当ni＞nt，则为品质良好禽蛋；当nt≥ni≥nr，则为品质较差禽蛋；当ni＜nr，则为变质禽蛋。

本具体实施方式中，判决模型按照如下步骤训练：

步骤201：收集已知品质的禽蛋作为训练样本，样本训练是为了得得模型参数和判决阈值令禽蛋的样本特征向量为变量x，按照所述步骤1～6计算出训练样本中各个禽蛋的样本特征向量；其中，训练样本中包含品质合格禽蛋、品质差禽蛋以及品质更差禽蛋；

步骤202：采用品质合格禽蛋训练一个三维空间的高斯模型，并计算模型参数以及∑，为样本特征向量均值，∑为样本特征向量协方差，

步骤203：将模型参数以及∑代入由品质合格禽蛋训练的高斯模型，得到判决模型如下：

步骤204：将品质差禽蛋以及品质更差禽蛋的特征向量代入判决模型，计算出若干模型值，选取品质差禽蛋模型值的平均值作为判决阈值nt，选取品质更差禽蛋模型值的平均值作为判决阈值nr，本具体实施方式中nt＝0.25，nr＝0.1。

本具体实施方式中的判决模型是在夏季、晴天、常温下进行训练的；若是季节、气候、环境温度发生较大变化，则需要采用训练样本重新对判决模型进行训练。

本具体实施方式中，步骤2中图像预处理按照如下步骤进行：首先，对红外图像采用3×3模板进行均值滤波降噪，3×3模板是指横向、纵向分别有3个像素点的像素区域，均值滤波能够有效减弱红外图像中的噪声干扰；然后对红外图像中各像素点的灰度值进行线性变换，将灰度值变换到0～255范围内，保证亮度和对比度的统一。这样，能够避免不同批次的禽蛋因为拍摄光线等原因造成拍摄出的红外图像的亮度和对比度存在较大差异，而引起的灰度值差异较大，从而使得基于红外图像进行的禽蛋品质鉴定能够具有统一的标准，使得禽蛋品质鉴定更加规范化、标准化。

本具体实施方式中，红外图像进行线性变换的方法如下：

找到红外图像中灰度值最大的像素点以及灰度值最小的像素点，获取红外图像的最大灰度值hmax和最小灰度值hmin；令线性变换前的灰度值为h，变换后的灰度值为g，构造如下方程组以求解线性变换参数k和b：

根据求解到的线性变换参数k和b，得到线性变换方程g＝kh+b，按照该线性变换方程对红外图像中各像素点的灰度值进行线性变换。

本具体实施方式中，步骤3中按如下步骤去掉不符合要求的连通区域：

步骤301：根据各个连通区域的像素点的数量计算各个连通区域的面积s′；根据拍摄的实际情况设定参数c1与c2；其中，c1为红外图像中最小禽蛋区域的像素点的数量，c2为红外图像中最大禽蛋区域的像素点的数量；本具体实施方式中，c1＝1000，c2＝3000；

步骤302：计算各个连通区域的外接矩形，设外接矩形的宽度为w，高度为h；外接矩形的宽度w是指在宽度方向上的像素点的数量，外接矩形的高度h是指在高度方向上的像素点的数量。

步骤303：将外接矩形的宽度w与参数c1与c2进行比较，去掉w＜c1和w＞c2的连通区域；这样能够去除掉红外图像中的长条形连通区域。

步骤304：根据禽蛋形状近似于椭圆的特点，去掉w/h＞1.5和h/w＞1.5的连通区域；

步骤305：根据禽蛋形状近似于椭圆的特点，去掉s′/(w*h)＜0.7的连通区域。

本具体实施方式中，步骤4中采用5×5模板对位置分布图进行中值滤波，以去除图像中孤立的噪声点。中值滤波可以做到既去除噪声又能较好保护图像的细节，从而获得较满意的复原效果，而且运算方便快速，特别适合点、线及尖顶细节很少的禽蛋图像。

本具体实施方式中，步骤2中采用最大类间方差法计算二值化阈值；将灰度值大于等于二值化阈值的像素点赋值为255，将灰度值小于二值化阈值的像素点赋值为0，从而得到二值图像。采用最大类间方差法能够在对背景和目标进行分割时，降低错分概率，从而保证后续步骤能够准确地在二值图像上分割出禽蛋区域和背景区域。