医药生产异物在线检查措施

1引言

我国有巨大的医药市场，为规范医药生产，国家药典对医药的质量检测，特别是可见异物检测进行了严格规定。可见异物来源…主要有以下几个方面：1)由于容器即安瓿瓶本身质量问题造成的在灭菌后出现的玻屑、小波点和白点等；2)生产过程工作人员不按照规定穿戴工作服等带入的头发丝等；3)生产过程中由于传输管道使用易脱落材质而产生的脱落物和黑点；4)在灌装过程中有未经净化的空气进入，而空气中可能带有大量的尘埃、纤维和色点等。我国的医药企业目前多数采用人工灯检方法进行检测，容易产生人眼疲劳有损人体健康，并且不能保证产品质量。目前目标检测与跟踪研究多集中在智能视频监控、交通监视和军事制导等领域，需要获取百余帧甚至上百帧图像由卡尔曼滤波器学习预测。如文献[2]中所应用的扩展卡尔曼方法(ukf)，文献[3]中的基于多智能体结构的应用处理，以及文献[4]中所应用的模型参考自适应方法等都是针对百余帧图像不断进行学习。文献[5]虽然是针对药液中可见异物的检测，但其使用的差分与形态学滤波相结合的方法，会遗漏许多小的点目标如碳屑等。文献[6]巾结合杂质位置和瓶壁厚度判断是否有杂质存在不能适用丁当溶液rfl古有较多数目的异物时的情况。田此如何设计出合适有效的检测与跟踪方法是灯检机。o研发中的一个关键问题。针对现有技术所存在的上述问题．本文采用概率加权阈值分割与差分相结合的检测方法，井应用夤尔曼滤渡嚣进行预测，从而达到快速跟踪，得到良好的检测和跟踪效果。实验结果表明，该方法可以有效的检铡出异物并进行跟踪以识别。

2检测对象

本文所针对的系统是全自动安瓿灯检机，检测对象主要是安瓿瓶装注射剂药渡中的可见异物。本检测系统针对不同类型的杂质共设计两种成像方法”，针对黑色异物”…．利用其不反光特性．在其背部给光；而针对玻屑、色点、白块等白色异物进行底部给光．形成的图像对比度大，提高检测效率。成像系统如图l所示。

3异物检测与跟踪

3．1运动目标检测

由于机械手的毛刺、环境巾的灰尘、瓶壁上的文字等，使得获取的图像上含有大量噪声干扰，如图2所示。本系统采用概率阈值分翱”“与帧间差分相结合的方法在去除噪声的同时达到检测的效果，同时利用卡尔曼滤波预测下一次目标可能位置，井在预测位置附近搜索，缩小搜索区域，减少卡尔曼滤波器的学习时问，提高检测速度。1)对初步去噪后的当前帧图像与第一帧图像^进行首次差分处理得到一次差分圈像A‘，即：M(x，，)=(x，y)一^(#，y)(1)式巾：，，y表示当前像素点的坐标，，(z，y)表示当前像素点的灰度值。2)对差分后的图像进行分割．进一步取出噪声。①求取图像的灰度概率分布”，井储存。即：t(n)=(∑，(z，y))／n(2)H：；}．．其概率为：P‰1；^(Ⅱ)／("x^)(3)式中：k(n)表示像素值为n的像素点个数(O≤n≤255)，”、^分别表示罔像的宽和高，p蜘，表示像素值为n的像隶点所出现的概率．设初始阕值I=l。②根据厨直t分割图像为两部分A(每个像素点的灰度值都大于E)和B(每个像素点的灰度值都小于‘)。则区域A、B的概率分别为：f叽2上Pm}J…(4)等5‘％2上肌一)区域A、B的平均灰度为：．IoJ地。i刍印“…(5)阻s2去．；．砷m，整幅图像的平均灰度为：p=qM+gp口(6)两个区域的总体方差为：，=吼(p．一p)2+如(脚一p)2(7)④按照最大类问方差””的准则，从。一255改变t，重复步骤@计算类间方差，使方差最大的I值即为区域分割的门限。④利用步骤③所求取的固值进行对步骤(1)中差分后的图像进行分剖。3)对分割后的图像进行第二次差分，进一步去除残留的静态干扰。二次差分后的阿像如图3所示。

3．2运动目标跟踪

检测目标后，引人卡尔曼滤波器进行预测，缩小搜索区域．加快搜索过程．提高跟踪效率．快速识别同标和干扰点，提高检测效率。

321建立系统模型假设一个线性离散系统1。“，其状态方程和观测方程表示为：一24置一-+札-z+q-(8)【z^=^r墨+口l如图4所示。其中，随机信号”H和％．分别表示过程激励噪声和观测噪声。假设它们为相互独立并且正态分布的白色噪声，即：，(")～Ⅳ(0，Q)伊(")一Ⅳ(0，R)(9)在检测中．选择在连续图像中所检测到的运动目标的位置和速度为被估计量，即J；卡尔曼滤渡的目的是要预测下一次目标可能出现的位置，因此输出z=x，从而日=，(4阶单位矩阵)。由于没有外界辅A．因此口=0；因此．针对本系统所建立的模型为：I置24x川+”¨(10)【Z‘=XI+”‘式中：工l=[h，y‘，“¨吖I]，机，y¨似I，叫‘分别表示所检侧到运动异物的质心坐标和在此坐标方向的速度；z。=[以，几]表示当前目标的位置；由于连续两帧的时间间隔只有lO眦，因此假设运动目标在两帧之问的运动轨迹为匀速直线运动，从而可得：rlOf0]A：lo1o7l(11)0Ol0l0000l。式中：r表示前后两帧的时间间隔，在本系统中时间为10如B，因此r=00l。在本系统中。Q、届、P均设置为对角矩阵，口=0lJ洱=0O¨；P(0IO)=100f(J为4阶单位矩阵)。

322卡尔曼预测过程1)获取初始状态信息，初始化预测的变量。经过3l异物检测之后，在所得到的检测结果图像上搜索连通区域，在初始状态下，将所搜索到的连通域作为初始检测目标。记录各自的坐标；。，‰、宽高信息％，^。和面积蛄(表示第0次预测的第•个目标)。由于第一次检测不能得到速度信息，因此设第一次的速度为l像素／秒．即似。、”h各自为l。2)卡尔曼滤波器预测过程①运用前一时刻的得到的状态置，(k—l=0)或状态估计茸．．(^一l>0)和协方差矩阵P¨来预测当前时刻的状态估计置和协方差矩阵只，即时间更新，其方程为：f墨=A邑一-：曰u川(12)L巧=AP‘一lAl+Q②运用时间更新方程在得到氍的情况下，根据预测的协方差矩阵^和观测噪声协方差矩阵也一。计算卡尔曼滤波器的增益甄：墨=P；日：(日。P：研+R。)‘1(13)③根据预测的状态估计墨和实际观测值z。修正系统的状态估计瓦并计算相应的协方差矩阵^即测量更新，其方程为：五=噩+瓦(zt一风墨’(14)L^=(，一甄日。)竹③经过步骤(2)得到预测位置互后，在后一帧图像上由鼍所决定的区域进行搜索，减小搜索范围，提高检测效率。

3．3目标识别

在经过卡尔曼滤波器预测之后，在后一副图像中的预测位置附近进行搜索，其流程如图5所示。在识别过程中的一个主要参数是搜索区域的确定，本文中采取自适应的搜索区域确定方法。在得出卡尔曼滤波器的预测位置后，计算预测位置与原位置的欧式距离噍(后表示第七个可疑目标)。经过分析和实验验证发现，溶液中的异物或干扰存在4种情况：1)对于检测之后遗留的瓶壁上的干扰以及溶液中的气泡噪声u孓埔】，由于在连续图像上，运动距离较小，而且其形状固定，因此面积变化较小，即有喀和S氓，都较小(其中s业，表示第f—l帧图像与第i帧图像中的同一个目标即第后个目标的面积变化幅度)。2)对于面积较大的异物，如玻屑、纤维等由于在拍摄过程中，异物仍在不停的运动，使得在图像上的面积会发生较大的变化。3)对于面积较小的异物，如碳渣等由于质量较重，会有较大的移动距离。因此，综合以上分析，为保证检测精度，不遗漏异物目标，确定搜索半径r龌=m缸(3&，1．5dI)。在确定搜索区域后，根据搜索半径的大小即可确定所跟踪目标是否为异物。分析如表l所示。

4实验结果

根据上述提出的检测疗法，对含有玻屑、纤维等异物的2lnL安瓿瓶装药范进行检测与跟踪识别．在运用差分法检测之后，在用卡尔曼滤波器所预测的范围内进行搜索，并进行识别判断。实验结果如图6、图7所示．表示连续的两幅图片巾当前目标位置和预测目标位置的对比，阿6巾表示所跟踪到的当前位置，*表示对下一次位置的估计，在图7中的+表示对当前位置的估计位置，•表示在预测位置附近进行搜索得到的目标。从图中可以看出在图3中的瓶壁干扰点已经被差分和阈值分割去除，而异物信息仍然得到了保留。图8为连续两次预测中，卡尔曼滤波器的预测目标位置与所检测到的目标实际位置在坐标的绝对数值上的对比。从图中可以看出，基本上预测位置与实际位置相吻合。

5结论

本文根据安瓿溶液中异物检测的实际情况，结合实验提出的视觉检测与跟踪方法经过试验验证既满足检测精度要求也满足了检测速度要求。本文巧妙地利用差分与概率阈值分割结合，即将检测与去除噪声的预处理相结合，压缩了检测过程，提高检测速度。另外本文中利用卡尔曼滤波器进行预测，确定和缩小搜索区域，从而在提高检测精度的同时保证检测速度。