食品中异物的检测方法技术

本发明专利技术涉及一种食品中异物的检测方法，其包括：采集被检食品的Ｘ射线图像；进行对比度增强处理；计算行灰度分布曲线中的峰的数量；获取结构宽度；获得模拟被检食品的无异物Ｘ射线图像；进行图像相减运算步骤，获得食品中异物目标图像。本发明专利技术采用数学形态学的运算来模拟正常食品的图像，然后再通过图像减影方式进行分割，从而可以快速地识别食品Ｘ射线图像中的异物目标，处理方法简单，可应用于在线检测。

【技术实现步骤摘要】

本专利技术涉及一种，尤其是一种通过对食品的X射线图像进行运算处理，提取混入食品中异物图像的检测方法。
技术介绍
目前，对于包装食品中混入异物的检测方法通常是利用X射线采集该食品的X射线图像，然后再通过人为干预，采用直接阈值或自适应阈值分割法，提取食品X射线图像中异物目标图像。由于食品种类较多，而且可能混入的异物的种类、在食品中的位置形态具有不确定性，这就导致了传统提取异物目标图像的分割方法具有较大的难度，尤其是在阈值的预先设定和检测过程中的图像提取都需要加入过多的人为因素，需要检测人员具有较为丰富的实际经验，容易出现误分割目标的现象，因此，传统的通过提取异物图像对不仅缺乏准确性，而且不适合连续在线检测，无法满足现代化食品生产企业的作业要求。
技术实现思路
本专利技术所要解决的技术问题在于针对上述包装食品中异物目标图像的提取方法所存在的不足提供一种，该方法运用数学形态学对食品的X射线图像进行运算处理，能够方便的提取异物目标图像。解决上述技术问题所采用的技术方案为一种，它由以下步骤组成步骤1采集被检食品的X射线图像；步骤2对所采集的X射线图像进行对比度增强处理，获得增强图像；步骤3自动计算增强图像每一行的灰度分布曲线中的峰的数量；所述自动计算依照下述方法进行的逐个选取任意像素点，并计算与该像素点相邻的两个像素点的灰度值，当所述相邻的两个像素点的灰度值同时小于或大于上述任意选取的像素点的灰度值时，该选取的像素点计为一个峰；步骤4利用所得到的峰的数量自动获取结构宽度；所述结构宽度的自动获取方法为当峰的数量小于3时，自动将所述结构宽度设为3；当峰的数量大于3时，首先计算所有峰处像素的元素宽度；然后将最大的元素宽度作为所述的结构宽度；步骤5利用所述结构宽度进行形态学运算，获得模拟被检食品的无异物X射线图像；步骤6将步骤5获得的模拟图像与步骤2中处理后的图像进行图像相减运算，获得食品中异物目标图像，完成食品中异物检测。由以上技术方案可以看出，本专利技术通过采用数学形态学的运算来模拟正常食品--即内部不存在异物的食品的图像，然后再通过图像减影方式进行分割，从而获得感兴趣的异物目标图像。本专利技术克服了传统的食品异物检测方法中存在的图像处理、提取困难以及误分割目标的不足，并将取得了良好的效果。本专利技术的有益效果在于，可以快速地识别食品X射线图像中的异物目标，处理方法简单，可应用于在线检测，也可应用于其他图像处理领域的感兴趣目标识别和提取。以下，通过具体实施方式并结合附图对本专利技术做进一步的详细说明。附图说明图1为本专利技术的流程示意图；图2为本专利技术所涉及一个具体实施例中一个图像行的灰度曲线。具体实施例方式当利用本专利技术所提供的方法对可能混有金属异物的食品进行检测时，其过程如图1所示的流程，其过程如下步骤1采集被检食品的X射线图像。将待检测包装食品以一定的匀速运动通过检测器，设定检测器分辨率为0.5mm×0.5mm，当待检测包装食品处于检测器的正上方0.5mm宽度范围内，其通过部分被扫描成像。此时，X射线扫描的积分时间为Tint，检测器像素分辨率宽度为Spixel，扫描次数位Ns，传送带运动速度位Vs，各参数之间存在的关系式为Tint=Spixel×Ns/Vs]]>在这一采集过程中，采用逐行多次扫描后取平均值的方法获得被检食品的降噪声X射线图像。步骤2对采集的X射线图像进行对比度增强处理，获得增强图像。X射线图像的对比度增强可以依据以下算式进行比度增强处理g(x,y)=Gb-GalnFb-lnFa+Ga]]>其中g(x，y)为变换后的灰度；f(x，y)为变换前的灰度值，Gb，Ga分别为变换后的灰度最大值和最小值，Fb，Fa为变换前灰度最大值和最小值。步骤3自动计算增强图像每一行的灰度分布曲线中的峰的数量。由于金属异物对X射线的阻挡较大，造成图像行的灰度分布呈现不均匀，其灰度分布呈曲线状。在曲线中存在多个峰。峰的计算方式为逐个选取任意像素点，并计算与该像素点相邻的两个像素点的灰度值，当相邻的两个像素点的灰度值同时小于或大于上述任意选取的像素点的灰度值时，该选取的像素点计为一个峰。图2所示为一个X射线图像中某一行的灰度分布曲线，从图中可以看出共有4个峰，每个峰所在像素点的相邻两个像素点的灰度值都同时小于该像素点的灰度值，因此，该图像行峰的数量为4。步骤4根据得到的峰的数量自动获取结构宽度。结构宽度的自动获取方法为当峰的数量小于3时，自动将所述结构宽度设为3；当峰的数量大于3时，依照以下算式首先计算所有峰处像素的元素宽度Width(i，k)＝|m-k|，if|Peak(i，m)-f(i，k)|＜ε其中，f(i，k)是距峰处像素Peak(i，m)处最近的像素点；ε是一个趋近于0的数字； Width(i，k)即为该峰像素的元素宽度。在完成计算后，将最大的元素宽度作为所述的结构宽度。步骤5利用所述结构宽度进行形态学运算，获得模拟被检食品的无异物X射线图像。步骤6对步骤4获得的模拟图像与步骤2中处理后的图像进行图像相减运算，并按照事先设定的阈值对相减运算后获得的图像进行分割而获得精确的食品中异物目标图像，完成食品中异物检测。由于X射线对不同异物的穿透能力不同，因此，图像的行灰度分布曲线中的峰可以呈现波峰或波谷。在图线中表现为曲线的上凸或者曲线的下凹。在利用形态学运算进行处理时，结构开运算能将灰度曲线上凸的部分去除保留凹的部分，闭运算正好与开运算相反。为消除波峰或波谷而采用的形态学运算包括结构闭运算或结构开运算。结构闭运算由数学形态学先膨胀后腐蚀运算得到，其运算式为膨胀运算依式(x)＝max(f(i-v)，...f(i)，...f(f+v))进行；腐蚀运算依式(x)＝min(f(i-v)，...f(i)，…f(i+v))进行；结构开运算为(x)。上述各运算式中，v是结构宽度，i为运算中心像素点位置。食品中混入的金属通常在灰度曲线中表现为凹部分，故使用结构开运算保留上凸部分去除凹部分从而模拟出正常图像。如果食品中混入塑料之类密度远小于食品密度的异物，则使用结构闭运算来模拟不含异物的食品图像。上述图像处理的运算过程完全是通过计算机进行自动运算处理，在开始检测时，只需要人工根据需要检测的异物项目预先设定计算参数以及阈值，就可以实现整个检测过程的全部自动化。最后所应说明的是以上实施例仅用以说明而非限制本专利技术的技术方案，尽管参照上述实施例对本专利技术进行了详细说明，本领域的普通技术人员应当理解依然可以对本专利技术进行修改或者等同替换，而不脱离本专利技术的精神和范围的任何修改或局部替换，其均应涵盖在本专利技术的权利要求范围当中。权利要求1.一种，其特征在于包括以下步骤步骤1采集被检食品的X射线图像；步骤2对所采集的X射线图像进行对比度增强处理，获得增强图像；步骤3自动计算增强图像每一行的灰度分布曲线中的峰的数量；所述自动计算依下述方法进行逐个选取任意像素点，并计算与该像素点相邻的两个像素点的灰度值，当所述相邻的两个像素点的灰度值同时小于或大于上述任意选取的像素点的灰度值时，该选取的像素点计为一个峰；步骤4利用自动计算得到的所述峰的数量自动获取结构宽度；所述结构宽度的自动获取方法为当峰的数量小于3时，自动将所述结构宽度设为3；当峰的数量大于3时本文档来自技高网...

【技术保护点】
一种食品中异物的检测方法，其特征在于包括以下步骤：步骤１：采集被检食品的Ｘ射线图像；步骤２：对所采集的Ｘ射线图像进行对比度增强处理，获得增强图像；步骤３：自动计算增强图像每一行的灰度分布曲线中的峰的数量；所述自动计算 依下述方法进行：逐个选取任意像素点，并计算与该像素点相邻的两个像素点的灰度值，当所述相邻的两个像素点的灰度值同时小于或大于上述任意选取的像素点的灰度值时，该选取的像素点计为一个峰；步骤４：利用自动计算得到的所述峰的数量自动获取结构宽 度；所述结构宽度的自动获取方法为：当峰的数量小于３时，自动将所述结构宽度设为３；当峰的数量大于３时，依照以下算式首先计算所有峰处像素的元素宽度：Ｗｉｄｔｈ（ｉ，ｋ）＝｜ｍ－ｋ｜，ｉｆ｜Ｐｅａｋ（ｉ，ｍ）－ｆ（ｉ，ｋ）｜＜ε其 中，ｆ（ｉ，ｋ）是距峰处像素Ｐｅａｋ（ｉ，ｍ）处最近的像素点，ε是一个趋近于０的数字，Ｗｉｄｔｈ（ｉ，ｋ）即为该峰像素的元素宽度；然后将最大的元素宽度作为所述的结构宽度；步骤５：利用所述结构宽度进行形态学运算，获得模拟被检食品的无异 物Ｘ射线图像；步骤６：将步骤５获得的模拟图像与步骤２中处理后的图像进行图像相减运算，获得食品中异物目标图像，完成食品中异物检测。...

【技术特征摘要】

【专利技术属性】
技术研发人员：韩东海，温朝晖，
申请(专利权)人：中国农业大学，
类型：发明
国别省市：11[中国|北京]