基于图像处理技术的矿石图像粒度检测算法的制作方法

本发明属于矿山开发中的矿石粒度检测分析技术领域，尤其是涉及一种基于图像处理技术的矿石图像粒度检测算法。

背景技术：

矿石的粒度分布信息是能够直接反映破碎机破碎效果的重要数据指标，高效率、高精度的矿石粒度分布检测是实现矿石破碎优化的基础，通过观察矿石颗粒的粒度，能清晰地了解破碎过程中各级破碎机的工作状况，其参数信息甚至会直接影响矿山开发过程中的矿石产率。目前基于计算机视觉与图像处理技术的矿石粒度信息的检测已经普遍应用在矿山开发过程中，以能实时获取破碎后矿石的粒度信息并反馈到plc控制系统，及时调节破碎机等选矿设备的参数，从而达到提高选矿效率、降低破碎机故障率、减少能耗、最大限度地利用矿产资源的目的，同时也会为矿山的生产带来巨大的经济效益。

但在其应用过程中也出现了有待于优化和解决的诸多问题。由于矿石颗粒通常具有数量多、易粘连、粒度差异大、形状不规则等特点。基于阈值分割的算法在很大程度上依赖于阈值的选择。针对矿石这类易粘连的颗粒，我们通常采用分水岭算法对其进行分割。但是单一的采用分水岭算法又容易引起过分割和欠分割现象的发生。

技术实现要素：

本发明的目的是提出了一种基于图像处理技术的矿石图像粒度检测算法，通过精准的矿石图像分割获得高精度、高准确率的粒度检测，降低破碎机的故障率，增高矿石产率。

本发明的目的是通过下述技术方案来实现的：

本发明的基于图像处理技术的矿石图像粒度检测算法，其特征在于包括如下步骤：

(1)利用工业相机完成对矿石图像的采集；

(2)对采集来的矿石图像进行预处理，依次进行图像灰度化、中值滤波去噪和二值化处理；

(3)对得到的二值化图像应用链码技术进行距离变换和形态学重构；

(4)采取canny算子边缘检测和基于区域生长的分水岭算法相结合的矿石图像分割算法；

(5)利用所得矿石参数及粒度模型对矿石粒度进行分析计算。

所述的步骤(1)和步骤(2)的具体步骤包括：

a.对矿石图像进行灰度化处理，将图像上各像素点的灰度值调节到rgb颜色分量的平均值，满足计算机图像分割处理的需要；

b.采用加权中值滤波去噪，把滤波窗口中灰度最大值用中值替换；

c.采取最大类间方差二值化算法，将两个邻域窗口最优阈值中较小的阈值认定为最佳阈值，以区分开目标矿石和矿石传送带周围的背景区域。

所述的步骤(3)包括把二值图中的前景目标灰度化，应用链码技术准确跟踪矿石颗粒边界；进行形态学重构变换，将距离变换后图像像素点灰度值与重构后的像素值做差，获取局部极大值区域，保护矿石图像的边缘信息。

所述的步骤(4)包括从x，y，45°，135°四个方向对canny算子进行边缘检测，计算像素点梯度幅值，抑制非极大值；结合医学领域中的区域生长算法弥补分水岭算法存在的过分割与欠分割问题。

所述的步骤(5)包括精确快速地提取分割后目标矿石图像的边缘信息，获得矿石颗粒周长l及投影面积s参数，利用kemeny的经验公式计算出目标矿石颗粒的粒径d，再根据椭球的体积公式，建立矿石体积的数学模型，从而求得矿石颗粒的体积参数；最后根据图像参数与实际矿石参数的比例关系，实现对矿石实际大小的测量。

本发明的优点：

本发明提出的基于图像处理技术的矿石图像粒度检测算法，将边缘检测算法与基于区域生长的分水岭算法结合,以此获得高精度的矿石粒度信息。并将参数信息反馈到破碎机的plc控制系统，调整排矿口宽度、进出料速度等一系列可能造成破碎机故障率增加的因素。通过精准的矿石图像分割实现高精度、高准确率的粒度检测，降低破碎机的故障率，增高矿石产率。

附图说明

图1为本发明基于图像处理技术的矿石图像粒度检测功能实现流程图。

具体实施方式

下面结合附图进一步说明本发明的具体实施方式。

如图1所示，本发明的基于图像处理技术的矿石图像粒度检测算法，其特征在于包括如下步骤：

(1)利用工业相机完成对矿石图像的采集；

(2)对采集来的矿石图像进行预处理，依次进行图像灰度化、中值滤波去噪和二值化处理；

(3)对得到的二值化图像应用链码技术进行距离变换和形态学重构；

(4)采取canny算子边缘检测和基于区域生长的分水岭算法相结合的矿石图像分割算法；

(5)利用所得矿石参数及粒度模型对矿石粒度进行分析计算。

所述的步骤(1)和步骤(2)的具体步骤包括：

a.对矿石图像进行灰度化处理，将图像上各像素点的灰度值调节到rgb颜色分量的平均值，满足计算机图像分割处理的需要；

b.采用加权中值滤波去噪，把滤波窗口中灰度最大值用中值替换；

c.采取最大类间方差二值化算法，将两个邻域窗口最优阈值中较小的阈值认定为最佳阈值，以区分开目标矿石和矿石传送带周围的背景区域。

所述的步骤(3)包括把二值图中的前景目标灰度化，应用链码技术准确跟踪矿石颗粒边界；进行形态学重构变换，将距离变换后图像像素点灰度值与重构后的像素值做差，获取局部极大值区域，保护矿石图像的边缘信息。

所述的步骤(4)包括从x，y，45°，135°四个方向对canny算子进行边缘检测，计算像素点梯度幅值，抑制非极大值；结合医学领域中的区域生长算法弥补分水岭算法存在的过分割与欠分割问题。

所述的步骤(5)包括精确快速地提取分割后目标矿石图像的边缘信息，获得矿石颗粒周长l及投影面积s参数，利用kemeny的经验公式计算出目标矿石颗粒的粒径d，再根据椭球的体积公式，建立矿石体积的数学模型，从而求得矿石颗粒的体积参数；最后根据图像参数与实际矿石参数的比例关系，实现对矿石实际大小的测量。

下面将结合说明书附图，针对具体实施方案对本发明做更详细完整的阐释说明。

1、要想在矿石图像分割步骤中得到理想的分割效果，首先要研究如何提高工业摄像机采集来的矿石图像的预处理效果，更好的消除潜在因素影响。

a.灰度化处理就是为了避免在一定的灰度范围中，由于各种不可抗因素而对采集的矿石图像造成无法分辨的结果。所以需要对矿石图像上各像素点的灰度值进行处理，将图像上各像素点的灰度值调节到rgb颜色分量的平均值，突出反映我们所需提取的有用信息，以达到图像增强的目的。

b.随后对得到的灰度图像进行图像滤波，由于我们最终目的是获得高精度的粒度参数，提高参数的准确度就需要我们标定出最接近理想状态的矿石图像轮廓边缘，所以在选取滤波方法这一环节，本发明采取了中值滤波这一方法去噪。中值滤波可以在消除图像噪声的同时，保护矿石颗粒的边缘信息。

中值滤波是一种非线性平滑滤波技术，选取一种二维模板，将滤波窗口n取奇数，这样的目的是为了愈加容易在区域中取得中心点。对输入的矿石图像，以任一像素为中心设置一个确定的邻域a，a的边长为2n+1(n＝0,1,2,…)。将邻域内各像素的强度值按大小顺序排列，将数字图像或序列中的某一点用该点的周围区域中各个点的中值代替，遍历整幅图像，从而消除孤立像素点。但传统的中值滤波算法每次都要将滤波窗口内各像素的强度大小进行排序，虽然可以获得较好的滤波去噪效果，但却耗时耗力，工作效率也相对低下，难以满足高效、快速的矿石粒度检测要求。所以本发明在传统中值滤波的基础上稍加改进，把滤波窗口中灰度最大值用中值替换。

中值滤波在i处信号样本的输出值用下式来定义：

y(i)＝med[x(i-n),...,x(i),...x(i+n)]

(1)

其中x(i-n)，...，x(i)，...x(i+n)为滤波窗口内采集的信号样本，n为正整数。邻域a可以用x(i,j)来表示:

x(i,j)＝maxx(i-1:i+1,j-1:j+1)(2)

在中值滤波图像中，将像素灰度值大于最优阈值的像素点当作是矿石，反之认为是背景。

c.运用二值化算法，基于形态学图像处理，利用结构元对矿石图像进行腐蚀并对图像中的矿石颗粒进行探测。如果结构元探测的区域内包含所有的前景色，那么就用前景色像素点替换中心位置的点，如果结构元探测的区域内不全包含前景色，那么就用背景色取代该点，以此把目标矿石和它周围的非矿石背景准确分离开来。二值化的关键一步是设置一个较为准确的阈值，将小于这个阈值的像素点调整成0，而大于等于这个阈值的像素点调整成255。

运用最大类间方差二值化算法，取两个邻域窗口最优阈值中较小阈值的为最佳阈值。

在此方法中，第一步要得到中心像素p的两个邻域窗口(即和)的最优阈值和两个窗口大小分别为n*n和m*m。

式中，t为灰度级范围取值[0,1]；l为最大灰度级(一般为255)；为最大类间方差。若则认定为最佳阈值。

第二步，要利用最优阈值将中值滤波之后的灰度图f(x,y)二值化，fb(x,y)为二值化的结果。

最终得到的矿石颗粒二值化图像才能称为是完全意义上的黑白图像。由于采矿现场条件恶劣，所需采集的矿石颗粒图像周围环境比较复杂，存在矿石颗粒堆积重叠、光照不均匀、灰尘颗粒密集等不利于粒度检测的情况，但经过改进的最大类间方差算法可以很好地适应该环境，从而得到更好效果的二值图像。

2、随后我们需要对上一环节得到的二值化图像进行距离变换，距离变换的结果不是另一幅二值图像，而是一幅灰度级图像，图像中每个像素的灰度值为该像素与距其最近的背景像素间的距离。本发明中采用链码技术按不同层次轮廓点灰度级递增的方式逐层对目标物体进行遍历，直至完成图像的距离转换。改进后的距离变换算法可以准确跟踪矿石颗粒边界，并有效地完成矿石颗粒这种易粘连物体的分离。与此同时，还弥补了传统形态学距离变换算法多次腐蚀、多次扫描、易耗时的不足。最重要的是提高了图像分割的处理速度和分割精度，具有极高的现实意义。

3、经距离变换后的矿石颗粒图像仍然存在多个峰头，而在去除多余峰头的过程中很有可能引起过分割问题，因此我们还需要对其进行形态学重构。消除矿石图像中剩余的细小图像和那些不包含任何结构因素的像素点。本发明对矿石图像先进行闭运算后腐蚀，然后进行开运算，以清除微小的边界颗粒。对极值点进行调整，使距离变换后的矿石图像中像素点的灰度值与重构变换后的图像灰度值做差，以获取局部极大值区域。更好的保护矿石图像的边缘信息，从而优化和完善矿石图像分割效果。

4、本发明针对现有分割算法的问题，提出了多特征融合的多尺度矿石图像分割算法。目前的分割算法主要是针对图像亮度特征进行分割，但对复杂矿石图像处理难以达到理想的分割效果。本发明中，由于矿料皮带上运送的矿石边界和轮廓并不是特别明显，而且在色彩上来说也是十分接近，这时候若是单采用分水岭分割算法，会造成矿石图像的过分割和欠分割问题，因此本发明中采用基于边缘检测和区域分割相结合的矿石图像分割算法进行优化。

a.本发明选取改进后的canny算子，克服了传统canny算法使用一组固定参数处理不同图像的缺点，缓解了平滑滤波与保持图像边缘细节，连续边缘与假边缘的矛盾。

改进后的canny算子从四个方向进行边缘检测，选择在3*3领域内计算梯度幅值，改善了2*2邻域内由于图像离散特性，易受噪声干扰的情况。梯度幅值计算公式如下：

gx(i,j)＝i(i,j+1)-i(i,j-1)(6)

gy(i,j)＝i(i+1,j)-i(i-1,j)(7)

g45°(i,j)＝i(i-1,j+1)-i(i+1,j-1)(8)

g135°(i,j)＝i(i+1,j+1)-i(i-1,j-1)

(9)

分别计算出x，y，45°，135°四个方向的偏导数gx(i,j)、gy(i,j)、g45°(i,j)及g135°(i,j)，并以此计算x水平方向和y垂直方向的差分方程：

可得矿石图像的像素点梯度幅值为：

改进后的3*3领域中像素的梯度方向θ(i,j)∈(-90°,90°)，将其领域划分为四个象限，随后沿像素的梯度方向判断该图像像素点的梯度幅值在邻域范围内是否为极大值。将灰度变化最大的保留下来，其他的不予保留，从而抑制非极大值。

b.为进一步削弱分水岭算法的过分割情况，利用医学领域中的区域生长算法对分水岭算法进行优化。

区域生长算法要先在矿石颗粒图像内设置种子区域，对种子区域的相邻像素进行估算，来确定他们是否应该被加入到种子区域，如果是，则被添加到该区域，使种子区域生长。而判断像素是否落入生长区间的准则是观察该像素的所有邻域像素点的灰度值是否均包含在预定值范围内。用邻域像素值来代替当前像素值，可以降低小结构被并入生长区域的概率。重复判断像素是否落入生长区间，直到没有新的像素被加入。

步骤一：设置阈值区间

首先要设定两个参数下限阈值tdown和上限阈值tup来构成生长区间，运用迭代算法思想求得阈值信息。

1)迭代法采取逼近的思想，求出图像的最大、最小灰度值，设为ratemax和ratemin，则初始阈值t0为：

2)根据t0，将图像的平均灰度值分为m1和m2两组。

3)利用下式分别求出两组的平均灰度值:

4)求出加权平均后的阈值为：

将步骤2)～4)重复n次，当阈值t不再变换时，结束迭代过程。此时的阈值设定为最佳阈值。阈值的选取会直接影响到种子生长区域的增涨，降低其机动性。

步骤二：定义初始生长区域

根据矿石颗粒图像像素点的灰度值，选取种子点坐标位置(x，y)，在阈值相同和不同时进行分割。选取其邻域半径，来判断相邻像素是否位于分割区域内。

将改进后的canny算子边缘检测和基于区域生长算法的分水岭算法相结合，不但可以精准定位矿石颗粒边缘，还能弥补分水岭算法存在的过分割与欠分割的不足。而且优化后的分割算法可以有效提升矿石分割效果，为之后的高精度粒度检测提供基础。

5、矿石图像粒度检测中所需要检测的参数主要有面积、周长、粒径和体积等。根据分水岭算法分割后的矿石图像，配合8链码技术，能够精确快速的提取目标矿石的边缘信息，根据边界信息就能够求得其周长lc及矿石颗粒的投影面积s参数。通过工业相机采集到的矿石图像中的目标矿石参数是二维数据，但是在计算破碎矿石的粒度信息时要使用的是矿石参数的三维数据，所以在实际统计过程中，我们需要从其二维数据中推测出相应矿石颗粒的三维数据信息。kemeny的经验公式就很好的解决了矿石图像数据信息从二维到三维的转化，并发现矿石颗粒的粒径的大小既非它的最大线性长度也非其最小线性长度，而是与其等效的最佳匹配椭圆的长、短轴a和b有关。

通过研究利用kemeny的经验公式便可以计算出目标矿石颗粒的粒径d为：

再根据椭球的体积公式，建立矿石体积的数学模型，可以求得颗粒的体积为:

通过上述方法得到的数据结果都是以像素为单位的，但是图像中像素的分辨率不是一定的，它会受外界环境的影响，工业相机拍照的高度不同，得到的像素结果就不同。在本发明中，用于图像采集的工业相机距离传输带的高度设置为1.3m。由于矿石图像的单位像素与实际的矿石单位尺寸成正比，可根据图像参数与实际矿石参数的比例关系，实现对矿石实际大小的测量。通过上述一系列过程，可以得到传送带上破碎矿石粒度分布信息。

综上所述，本发明所提及的矿石粒度检测方案，可以计算出高精度、高准确率的矿石粒度信息，从而反馈到破碎机plc控制系统，调节破碎机相关硬件设备参数，降低破碎机的故障率，增加矿石产率。

本发明基于图像处理技术的矿石图像粒度检测算法，将边缘检测算法与基于区域生长的分水岭算法结合,以此获得高精度的矿石粒度信息，并将参数信息反馈到破碎机的plc控制系统，调整排矿口宽度、进出料速度等可能造成破碎机故障率增加的参数，通过精准的矿石图像分割获得高精度、高准确率的粒度检测，降低破碎机的故障率，增高矿石产率。