基于图像技术的粉体混合均匀性检测方法与流程

本发明涉及粉体混合均匀性检测领域，具体涉及一种基于图像技术的粉体混合均匀性检测方法，为各行业粉体混合分析及终点判断提供了方法。

背景技术：

粉体混合是一个复杂的随机过程，混合均匀性的检测方法一直是困扰着人们的棘手问题。随着时代的发展，依赖人工经验判定混合均匀性的日子已经过去。用科学、数量的形式来判定混合均匀性，就是粉体混合的定量分析。要做到定量分析，必然要有采样、预处理、检测、统计分析(数据处理)等过程，从而得到单一或数个的量值来表达粉体的混合均匀性。

混合均匀性是指在外力的作用下，各种粉体相互混合，使之在任何容积里每种组分粉体微粒均匀分布的程度。混合过程中粉体的运动十分复杂，如何正确评价粉体的混合均匀性一直是研究的热点。粉体混合是工业生产中常见的操作，混合效果的好坏直接决定着产品的质量。目前在化工、制药、建筑等领域都有非常广泛的应用，是药品、混凝土等产品生产最关键的操作工序之一。

技术实现要素：

本发明所要解决的技术问题是针对现有技术的不足提供一种基于图像技术的粉体混合均匀性检测方法。

本发明的技术方案如下：

一种基于图像技术的粉体混合均匀性检测方法，包括以下步骤：

s1：粉体样本采集；

s2：粉体图像采集；

s3：图像预处理：将采集得到的粉体图像进行图像预处理，依次包括图像灰度化、图像降采样、高斯滤波和图像锐化；图像灰度化将粉体图像变换为灰度矩阵；图像降采样将图像分辨率降低以减少程序运行时间；高斯滤波用于过滤图像采集过程中可能出现的误差；图像锐化用于增大不同粉体原料之间的色差，提高分析结果的准确性；

s4：灰度直方图提取特征参数：对预处理后的粉体图像，采用灰度直方图准确提取图像特征参数，该特征参数反映了粉体图像中每一灰度级与其出现频率间的统计关系；

s5：检测图像灰度级分布的相似性：根据s4得到的结果，首先根据特征参数计算不同采样点粉体图像之间的变权重欧氏距离，并通过每组图像之间的距离表征图像间灰度级分布的相似性；

s6：根据实际粉体特性和加工要求设定阈值并判断混合均匀性：根据不同粉体特性，设定合适阈值，通过图像间灰度级分布的相似性判断粉体混合均匀性，并进行混合终点判断。

所述的基于图像技术的粉体混合均匀性检测方法，所述的图像灰度化：采用平均值法实现灰度化：求出每个像素点的r、g、b三个分量的平均值，然后将这个平均值赋予给对应灰度图像该像素点的灰度值：

f＝(r+g+b)/3

式中，f为灰度图像每个像素点的灰度值；r、g、b为彩色图像中每个像素点颜色的的三个分量。

所述的基于图像技术的粉体混合均匀性检测方法，所述的高斯滤波：根据选择的固定窗口大小及窗口内任意像素与中心像素点的距离，利用高斯函数实现系数权值的分配，即高斯滤波可表示为：

式中，wx,y表示中心像素(x,y)的m×m(m为奇数)大小的领域；wd为空间距离相似度权重因子；

所述的图像降采样：图像降采样即是采样点数减少；对于一幅n*m的图像来说，如果降采样系数为k，则即是在原图中，每行每列每隔k个点取一个点组成一幅图像。

所述的基于图像技术的粉体混合均匀性检测方法，所述的图像锐化：采用线性锐化滤波，通过灰度变换函数压缩低端扩展高端，提高粉体颗粒色差的对比度。

所述的基于图像技术的粉体混合均匀性检测方法，所述的改进的变权重欧式距离的算法为：

s51：提取图像特征参数：通过灰度直方图算法提取图像特征参数，该特征参数反映了粉体图像中每一灰度级与其出现频率间的统计关系；

若图像具有l级灰度,则大小为m×n的灰度图像f(x,y)的灰度直方图hist[0,…,l-1]可用如下计算获得：

(1)初始化hist[k]＝0,1,…l-1；

(2)对所有的x＝0,1,…m-1和y＝0,1,…n-1计算：hist[f(x,y)]＝hist[f(x,y)]+1；

(3)归一化得到灰度直方图的向量：h[f(x,y)]＝hist[f(x,y)/(m×n)]；

s52：中心图像特征参数：根据上步提取的多图像特征参数，采用均值法计算得中心图像的特征参数：

式中，i为每组检测中不同采样点粉体图像的个数，取i＝36；hi[f(x,y)]为第i幅粉体图像的特征参数向量；

s53：归一化：计算中心图像特征参数各灰度级像素点个数及占比；

s54：根据占比确定权重：中心图像特征参数各灰度级所占比例即为该灰度级的权重系数；

s55：各图像与中心图像带权重的欧式距离：

式中，j为灰度级数；ε(i)为基于多图像特征参数的平均特征参数各灰度级像素点所占的比例；

s56：距离均值：

式中，s(h1,h2,…,h36)为多张图像距离均值，即为图像间灰度级分布的相似性；s的值越靠近1，表示多张图像灰度级分布越相似；越靠近0，表示多张图像灰度级分布越不相似。

所述的基于图像技术的粉体混合均匀性检测方法，选择数字显微镜采集粉体图像，将待观察的粉体样本置于物镜正下方，调整至合适的光源，调节焦距找到所要观察的样品，切换到高倍模式，通过调节样品位置，找到所要观察的视场。选好视场后，放大到需要的聚焦倍数，使粉体颗粒图像达到最佳状态，然后拍摄图像。旋转转盘，便于拍摄不同位置的粉体图像，每旋转10度采集一张图片，每组样本采集36张图像，构建数据库，实时保存粉体图像。

采用上述方案，本发明利用图像技术，通过分析图像灰度级分布的相似性，实现混合均匀性在线分析，并判断混合终点。具有实时性强、适用性广并可提高经济效益的优点。

附图说明

图1为基于图像处理的粉体混合均匀性检测方法的流程图；

图2为图像预处理流程图；

图3为改进的变权重欧式距离算法流程图。

具体实施方式

以下结合具体实施例，对本发明进行详细说明。

如图1所示，为基于图像处理的粉体混合均匀性检测方法的流程图，包括：

s1：粉体样本采集。

s2：粉体图像采集：选择数字显微镜采集粉体图像。将采样装置固定在载物台上，将待观察的粉体样品置于物镜正下方，调整至合适的光源，调节焦距找到所要观察的样品，切换到高倍模式，通过调节样品位置，找到所要观察的视场。选好视场后，放大到需要的聚焦倍数，使粉体颗粒图像达到最佳状态，然后拍摄图像。拍摄不同位置的粉体图像，每组样品拍摄36张照片，构建图像数据库，实时保存粉体图像。

s3：图像预处理：为了解决图像采集中可能出现的误差和图像处理时数据量过大的问题，根据采集得到的粉体图像，将采集得到的粉体图像进行图像预处理。其中包括图像灰度化、图像降采样、高斯滤波和图像锐化。

首先是图像灰度化将粉体图像变换为灰度矩阵；然后进行图像降采样，将图像分辨率降低以减少程序运行时间；高斯滤波是过滤图像采集过程中可能出现的误差；图像锐化是增大不同粉体原料之间的色差，提高分析结果的准确性。

s4：灰度直方图提取特征参数：对预处理后的粉体图像，采用灰度直方图准确提取图像特征参数，该特征参数反映了粉体图像中每一灰度级与其出现频率间的统计关系。该方法运算简单,具有平移、旋转和缩放不变性。

s5：检测图像灰度级分布的相似性：根据s4得到的结果，首先根据特征参数计算不同采样点粉体图像之间的变权重欧氏距离，并通过每组36张图像之间的距离表征图像灰度级分布相似性。为了验证该距离的准确性，计算不同灰度级像素个数的标准差,归一化处理后对不同等级像素标准差求平均值，通过该均值来表征图像间灰度级分布的相似性。对比两种数据处理方法，得到最终不同采样点的粉体图像灰度级分布相似性结果。

s6：根据实际粉体特性和加工要求设定阈值并判断混合均匀性：根据不同粉体特性，设定合适阈值将图像灰度级分布相似性检测结果与粉体混合均匀性相联系。通过不同采样点粉体图像灰度级分布相似性判断粉体是否混合均匀，并进行混合终止判断。若混合均匀，则停止混合；若没有混合均匀，则间隔五分钟后重新取样进行分析判断。

s7：为了满足操作简单、实时显示的粉体混合均匀性检测的要求，在上位机建立方便、直观、友好的人机界面，分为主界面、趋势、存储记录和参数设置四个画面，实时显示、记录、存储相关数据、变化趋势以及均匀性判断结果。

主界面用来显示数据库图像调用、图像灰度级分布相似性结果和混合均匀性判断；趋势画面是由当前时刻的图像灰度级分布相似性检测的特征参数和本次混合过程中过去时刻的图像灰度级分布相似性检测的特征参数以图形点连成的曲线图，通过趋势显示画面可以直观明了的掌握粉体混合均匀性的变化规律；存储记录画面是用来对历史混合过程的图像灰度级分布相似性检测和混合均匀性判断的参数进行存储；参数设置画面是根据实际的粉体原料，具有特定权限的人员正确登录后可以对判断混合均匀性的阈值进行设定。

如图2所示，为图像预处理流程框图。通过该方法，可以解决图像采集中可能出现的误差和图像处理时数据量过大的问题，提高实时性；并且增大不同粉体原料之间的色差，提高准确性。

s1：图像灰度化:采用平均值法实现灰度化：求出每个像素点的r、g、b三个分量的平均值，然后将这个平均值赋予给对应灰度图像该像素点的灰度值：

f＝(r+g+b)/3

式中，f为灰度图像每个像素点的灰度值；r、g、b为彩色图像中每个像素点颜色的三个分量。

s2：高斯滤波：根据选择的固定窗口大小及窗口内任意像素与中心像素点的距离，利用高斯函数实现系数权值的分配，即高斯滤波可表示为：

式中，wx,y表示中心像素(x,y)的m×m(m为奇数)大小的领域；wd为空间距离相似度权重因子。

s3：图像降采样：图像降采样即是采样点数减少。对于一幅n*m的图像来说，如果降采样系数为k，则即是在原图中，每行每列每隔k个点取一个点组成一幅图像。

s4：图像锐化：由于大多数粉体图像其灰度分布集中在较窄的区间,引起图像细节不够清晰。经过锐化可使图像的灰度距离拉开或使灰度分布均匀,从而增大反差,图像细节清晰,易于寻找其特征值。图像锐化是补偿图像的轮廓，增强图像的边缘及灰度跳变的部分，使图像变得清晰。本专利采用线性锐化滤波，通过灰度变换函数压缩低端扩展高端，提高粉体颗粒色差的对比度。

如图3所示，为改进的变权重欧式距离算法流程图。通过该方法，可解决多张图像灰度级分布相似性分析的问题。本专利通过度量多幅图像与中心图像特征参数的变权重欧氏距离来计算图像灰度级分布的相似性。

s1：提取图像特征参数：通过灰度直方图算法提取图像特征参数，该特征参数反映了粉体图像中每一灰度级与其出现频率间的统计关系。

若图像具有l(通常l为256,即8位灰度级)级灰度,则大小为m×n的灰度图像f(x,y)的灰度直方图hist[0,…,l-1]可用如下计算获得:

(1)初始化hist[k]＝0,1,…l-1；

(2)对所有的x＝0,1,…m-1和y＝0,1,…n-1计算：hist[f(x,y)]＝hist[f(x,y)]+1；

(3)归一化得到灰度直方图的向量：h[f(x,y)]＝hist[f(x,y)/(m×n)]。

s2：中心图像特征参数：根据上步提取的多图像特征参数，采用均值法计算得中心图像的特征参数。

式中，i为每组检测中不同采样点粉体图像的个数，本专利取i＝36；hi[f(x,y)]为第i幅粉体图像的特征参数向量。

s3：归一化：计算中心图像特征参数各灰度级像素点个数及占比。

s4：根据占比确定权重：中心图像特征参数各灰度级所占比例即为该灰度级的权重系数。

s5：各图像与中心图像带权重的欧式距离：

式中，j为灰度级数；ε(i)为基于多图像特征参数的平均特征参数各灰度级像素点所占的比例。

s6：距离均值：

式中，s(h1,h2,…,h36)为多张图像距离均值，即为图像灰度级分布的相似性；s的值越靠近1，表示多张图像灰度级分布越相似；越靠近0，表示多张图像灰度级分布越不相似。

上述方法首次将图像灰度级分布相似性用于粉体均匀性的检测，并取得良好的效果。使用图像之间的灰度级分布相似性检测粉体的混合均匀性，可以实现实时在线的监测。在工业领域有良好的应用前景。

应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。