一种基于滑动窗口算法的纹理特征度量方法技术

一种基于滑动窗口算法的纹理特征度量方法，在泡泡浮选领域，本发明专利技术提出一种基于滑动窗口算法的纹理特征度量方法，此方法基于现场设置的数字图像采集系统提取泡泡图像，提出了颗粒区域的概念，并运用滑动窗口算法对泡泡表面的颗粒区域进行准确提取，将所提取的颗粒区域按中心点位置进行分级，继而对颗粒区域之间的密集程度进行度量，定义了纹理特征颗粒密集度，用以反映整幅图像的纹理特征。有效弥补了传统的纹理特征提取方法没有考虑泡泡表面颗粒的缺陷，从而可以更准确的判断工况并有效指导加药。

A texture feature measurement method based on sliding window algorithm

【技术实现步骤摘要】
一种基于滑动窗口算法的纹理特征度量方法
本专利技术属于泡沫浮选
，具体涉及一种锌浮选过程中的纹理特征度量方法。
技术介绍
泡沫浮选是国内外广泛应用的一种选矿方法，该方法能依据矿物表面亲水性与疏水性的差异，有效地将目标矿物分离出来。泡沫浮选过程将目标矿物与其共生的脉石研磨成合适大小的颗粒然后送入浮选槽中，通过添加药剂调整不同矿物颗粒表面性质同时在浮选过程中不断地搅拌和鼓风，使矿浆中形成大量具有不同尺寸、形态、纹理等特征信息的气泡，使有用矿物颗粒粘附在气泡表面，气泡携带矿物颗粒上升至浮选槽表面形成泡泡层，脉石矿物留在矿浆中，从而实现矿物分选。由于浮选工艺流程长、内部机理不明确、影响因素多、涉及变量多且非线性严重、工艺指标不能在线检测等原因，一直以来，浮选过程主要依靠人工肉眼观察浮选槽表面泡泡状态来完成现场操作，这种生产方式主观性强，难以实现浮选泡泡状态的客观评价与认知，造成浮选生产指标波动频繁、矿物原料流失严重、药剂消耗量大、资源回收率低等情况发生，特别在高品位矿物资源日益匮乏的今天，浮选矿源成分复杂、矿物品味低，浮选生产人工操作更加难以有效进行。将机器视觉应用到浮选过程中，利用数字图像处理技术对浮选泡泡图像进行分析，能实现对泡泡状态的客观描述，再进一步寻找与分析泡泡特征参数与工艺指标的关系，从而推进了浮选过程的生产自动化。浮选泡泡随着浮选状态的不同而呈现出特殊的纹理状态，泡泡图像的纹理是泡泡表面粗糙度、对比度和黏性的综合体现，它与加药量、通气量等浮选生产操作变量及精矿品位、尾矿含量等浮选生产指标密切相关，当前的泡泡纹理信息提取方法主要是提取局部特征，存在提取的精度不够且纹理提取过程中没有考虑到泡泡表面颗粒等问题，难以准确反映工况，而实际上，泡泡表面常常附着一些矿石或者杂志小颗粒，造成泡泡表面粗糙不平，而这些小颗粒的出现的数量与分布密集度与锌精矿品位有密切关系，针对之前的研究没有考虑到泡泡表面颗粒的问题，提出一种基于滑动窗口算法的纹理特征度量方法，此方法基于现场设置的数字图像采集系统提取泡泡图像，继而对泡泡表面的颗粒区域进行准确提取，并对颗粒区域之间的密集度进行量化度量，定义了新的纹理特征颗粒密集度，用以反映整幅图像的纹理特征，有效避免了传统的纹理特征提取方法的局限性，从而更准确的判断工况并有效指导加药。
技术实现思路
本专利技术的目的是提供一种基于滑动窗口算法的纹理特征度量方法，在浮选生产中，浮选泡泡表面纹理是反映矿石品位的重要视觉信息，与浮选工况密切相关，直接反映泡泡层的矿化程度。本专利技术针对目前存在的浮选泡泡图像的纹理特征提取没有考虑泡泡表面颗粒的问题，运用滑动窗口算法定义了纹理特征度量方法颗粒密集度。该方法首先对泡泡图像进行分割，然后对泡泡表面颗粒区域进行提取，并对颗粒密集度如何反映工况进行了分析说明，表明它能更准确地对矿物品位进行调控并指导浮选生产。本专利技术采用的技术方案步骤如下：步骤一：利用浮选现场图像采集系统收集锌浮选的泡泡视频并将泡泡视频转换为连续图像，对采集到的锌浮选图像数据进行数据预处理,如下：1)剔除超出正常变化阈值的错误数据；2)剔除不完整的数据；步骤二：将泡泡图像由RGB彩色图像转化为灰度图像，得到图像的灰度矩阵A：egf表示的是灰度图像中的每一个像素点对应的灰度值，其中，g∈N，f∈N，N∈(400,800)；步骤三：泡泡图像中，泡泡常规的形态是表面光滑，高亮点位于单个泡泡凸起曲面的顶端，高亮点区域呈现最小灰度值，而以高亮点为中心，灰度值向下逐渐递增，在到达泡泡边界时达到最大值；而实际上，泡泡常常附着一些矿石或者杂志小颗粒，造成泡泡表面粗糙不平，而这些小颗粒的出现的数量与分布密集度与加药量和锌精矿品位有关系；首先，对泡泡进行分割，采用分水岭的方法将泡泡分割得到h个单个的泡泡，并储存各个泡泡的灰度矩阵，得到单个泡泡的灰度矩阵集合B＝{b1,b2,b3,...,bλ,...,bh}，bλ是第λ个泡泡灰度数值矩阵，h是分割后泡泡总的个数，在分割后的单个泡泡图像集合当中，筛选泡泡尺寸大于1200像素值的泡泡，即为感兴趣区域，对于筛选后的泡泡，用单个泡泡的灰度均值代替此泡泡高亮点部分的灰度值，得到待检测泡泡灰度矩阵集合C＝{c1,c2,c3,...,cε,...,cK}，K是符合泡泡尺寸要求的泡泡个数。步骤四：颗粒区域的检测：对于筛选后得到的的单个泡泡，分别对每一个泡泡进行以下操作：S1：取31×31像素点大小的子块作为窗口，沿着泡泡边界由上到下由左往右进行滑动，步长设定为5；S2：采用滑动窗口算法将单个泡泡图像分割为31×31的矩形子块区域，并对分割后的矩形子块区域进行判定，筛选出颗粒区域，判定的方法具体为：(1)将子块区域的灰度值变化用折线显示，横轴表示列数，纵轴表示灰度值，每一条折线表示灰度矩阵中对应行的灰度值变化，则共有31条折线，所有的折线同时显示在一幅图上；(2)在所有矩形子块区域的灰度图像上标记灰度值为205的灰度线，然后进行以下判断：Q1:判断是否有T(T≥3)条列数差低于2的折线从下向上穿过205灰度线，若无则表示此图像不包含颗粒区域，若有则标记最左侧的穿越点为E1，进而执行Q2判断；Q2:计算这T条折线穿过205灰度线时的斜率k1,k2,...,kη,...,kt，若满足条件kη≥3.6，则执行Q3判断，否则排除此图像；Q3:在该矩形子块区域的灰度图像上，判断这T条折线是否又自上而下穿过205灰度线，若无排除此图像，若有则标记穿回时最右侧穿越点为E2，进而执行Q4判断；Q4:计算E＝E2-E1,若E≤6.8则此矩形子块区域内包含完整的颗粒，定义满足以上条件的矩形子块区域为颗粒区域。步骤五：定义颗粒密集度：所有的颗粒区域均为31×31像素点大小的子块，定义此子块的几何中心即子块灰度矩阵中坐标为(16,16)的像素点为颗粒区域中心点，此中心点在颗粒区域所对应泡泡的灰度矩阵中的位置用笛卡尔坐标表示为：即为颗粒区域中心点位置。重复步骤四搜索出待检测泡泡集合中单个泡泡灰度数值矩阵中的所有颗粒区域，对于泡泡cε,统计颗粒区域的数量Lε，然后确定出各个颗粒区域的中心点坐标以不同颗粒区域中心点之间的直线距离定义一级、二级邻域区域以及其他邻域区域,对颗粒区域Pu,u∈a，定义如下：其中v∈a且v≠u；记颗粒区域Pu的一级邻域区域个数为qu,二级邻域区域个数为su，其他邻域区域个数为wu，则有qu+su+wu＝Lε-1；一级邻域区域个数、二级邻域区域个数以及其他邻域区域个数所占权重分别设定为0.6、0.3、0.1，则对于泡泡cε，定义纹理特征量颗粒密集度Zε，表达式如下：定义整幅图像的颗粒密集度G，如下式所示：步骤六：通过颗粒密集度判定工况：当处于状态①时，泡泡表面纹理较细，药剂过量，泡泡中承载的矿物粒子超过了泡泡的承载量而使泡泡大量破碎，药剂浪费比较严重而且精矿品位低；...

【技术保护点】
1.一种基于滑动窗口算法的纹理特征度量方法,其特征在于，包括以下步骤：/n步骤一：利用浮选现场图像采集系统收集历史时刻的锌浮选的泡泡视频并将泡泡视频转换为多帧的连续图像，对采集到的锌浮选图像数据进行数据预处理；/n步骤二：将泡泡图像由RGB彩色图像转化为灰度图像，得到图像的灰度矩阵A/n

【技术特征摘要】
20190722 CN 20191066174551.一种基于滑动窗口算法的纹理特征度量方法,其特征在于，包括以下步骤：
步骤一：利用浮选现场图像采集系统收集历史时刻的锌浮选的泡泡视频并将泡泡视频转换为多帧的连续图像，对采集到的锌浮选图像数据进行数据预处理；
步骤二：将泡泡图像由RGB彩色图像转化为灰度图像，得到图像的灰度矩阵A



egf表示的是灰度图像中的每一个像素点对应的灰度值，其中，g∈N，f∈N；
步骤三：对泡泡进行分割，采用分水岭的方法将泡泡分割得到h个单个的泡泡，并储存各个泡泡的灰度矩阵，得到单个泡泡的灰度矩阵集合B＝{b1,b2,b3,...,bλ,...,bh}，bλ是第λ个泡泡灰度数值矩阵，在分割后的单个泡泡图像集合当中，筛选出尺寸大于1200像素值的泡泡,记作C＝{c1,c2,c3,...,cε,...,cK}，K是符合泡泡尺寸要求的单个泡泡的个数；
步骤四：颗粒区域的检测
S1：对于分割筛选后的泡泡，取r×r像素点大小的子块作为窗口，沿着泡沫边界由上到下由左往右进行滑动，步长设定为5；
S2：采用滑动窗口的方法将单个泡沫图像分割为r×r的矩形子块区域，并对分割后的矩形子块区域进行判定筛选出颗粒区域；
步骤五：重复步骤四搜索出待检测泡泡集合中单个泡泡灰度数值矩阵中的所有颗粒区域，对于泡泡cε,统计颗粒区域的数量Lε，然后确定出各个颗粒区域的中心点坐标以不同颗粒区域中心点之间的直线距离定义一级、二级邻域区域以及其他邻域区域,对颗粒区域Pu,u∈a，定义如下：



其中v∈a且v≠u；
记颗粒区域Pu的一级邻域区域个数为qu,二级邻域区域个数为su，其他邻域区域个数为wu，则有qu+su+wu＝Lε-1；一级邻域区域个数、二级邻域区域个数以及其他邻域区域个数所占权重分别设定为0.6、0.3、0.1，则对于泡泡cε，定义纹理特征量颗粒密集度Zε，表达式如下：



定义整幅图像的颗粒密集度G，如下式所示：



步骤六：通过颗粒密集度判定工况：



当处于状态①时，泡泡表面纹理较细...

【专利技术属性】
技术研发人员：唐朝晖，罗金，张国勇，李耀国，范影，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南;43