一种基于泊松方程的矿物图像对比度调节方法
【技术领域】
[0001] 本发明涉及一种图像增强方法，尤其适用于对矿物品质研宄使用的基于泊松方程 的矿物图像对比度调节方法。
【背景技术】
[0002] 图像处理是研宄矿石品位的一个重要方法，近年来，扫描电子显微镜（SEM)越来 越多的用于矿物图像获取，其通过二次电子、背散射电子成像技术，可以清晰的反映出矿石 表面的立体形貌，较其他观察工具而言有很多的优点：超高的分辨率；景深大，视野大，图 像富有立体感，可直接观察各种式样凹凸不平表面的细微结构；样品制备简单等。矿物颗粒 往往具有凹凸不平、表面纹理细节较为复杂且颗粒大小参差不齐等特点，所以SEM成为目 前观察矿物颗粒的最常用工具。
[0003] 但是由于扫描电子显微镜成像原理及一些人为因素，导致其拍摄的图像有一些常 见的图像缺陷：荷电效应、损伤、边缘效应。具体表现为图像出现亮暗异常，即图像对比度有 很大问题，虽然通过调节扫描电镜工作参数等措施可以减轻对图像对比度的影响，但不能 完全消除，所以很有必要对扫描电子显微镜扫描的图像进行处理，调整图像对比度，从而改 善图像质量，增强图像纹理细节。
【发明内容】
[0004] 针对上述问题，本发明的目的是提供一种方法简单，调节效果好的基于泊松方程 的矿物图像对比度调节方法。
[0005] 为实现上述目的，本发明的基于泊松方程的矿物图像对比度调节方法，其包括以 下步骤：
[0006] a.使用扫描电子显微镜获取对被测矿物的扫描图像，并采用高斯滤波器对获取到 的扫描图像进行去噪处理，从而得到去噪后的矿物图像f ;
[0007] b.使用公式：U = log(f+l)，将去噪后的矿物图像f进行取对数操作，得到对数图 像u ;
[0008] c.利用公式：▽ u (X，y) = (ux，uy) = (u (x+1，y) -U (X，y)，u (X，y+1) -U (X，y))，对对 数图像u进行前向差分操作，得到对数图像u的梯度，式中，u (X，y)为对数图像u的第X行 第y列的像素点灰度值；
[0009] d.设定矿物图像f中过暗区域和过亮区域的判别标准：设定过暗区域像素灰度临 界值为T 1，过亮区域像素灰度临界值为T2，由于8位灰度数字图像中，像素灰度值低于50的 区域或像素灰度值大于200的区域，其图像细节人眼均无法识别，所以选取默认值T 1= 50， 选取默认值T2= 200 ;增加矿物图像f的高斯噪声，由于人眼只能识别矿物图像f的高斯噪 声标准方差在4到10之间，因此选取矿物图像f的梯度倍乘因子α = 10/4 = 2. 5 ;
[0010] e.取矿物图像f中的每个像素点（x，y)，由该像素点及与该像素上下左右邻接 的四个像素点构成像素点（X，y)的一个子集：N (X，y) = {(X，y)，（x+1, y)，（X-1, y)，（X ，y+l)，（x，y-l)}，若该子集中的每个像素点的灰度值都小于像素灰度临界默认值T1，即 {f(l，k) ST1, (l，k) eN(x，y)}，则判断该子集区域为矿物图像过暗区域，若该子集中的每 个像素点的灰度值都大于像素灰度临界默认值T2，即{f(l，k)多T 2,（l，k) eN(x，y)}，则 判断该子集区域为矿物图像过亮区域，其余像素点子集不做处理；
[0011] f.选取过暗区域和过亮区域的子集，利用公式：V = α ▽ U(x, y) = ( a Ux, a Uy)， 计算过暗区域和过亮区域对比度改善之后的梯度ν，分别将过暗区域和过亮区域的对数图 像梯度Vu(X，y)乘以梯度倍乘因子α，使乘以梯度倍乘因子α后的矿物图像过暗区域和 过亮区域的对比度提高，拉伸过暗区域和过亮区域的灰度动态范围，在梯度域提高细节，增 强了矿物图像过暗区域和过亮区域的纹理细节，改善视觉效果，得到调节对比度后的矿物 对数图像的梯度V ;
[0012] g.经过步骤f中梯度倍乘因子α调节后，得到调节对比度后的矿物对数图像的梯 度V为
【主权项】
1. 一种基于泊松方程的矿物图像对比度调节方法，其特征在于包括以下步骤： a. 使用扫描电子显微镜获取对被测矿物的扫描图像，并采用高斯滤波器对获取到的扫 描图像进行去噪处理，从而得到去噪后的矿物图像f ; b. 使用公式：u = log(f+l)，将去噪后的矿物图像f进行取对数操作，得到对数图像u ; C.利用公式：▽ u(x, y) = (ux, uy) = (u(x+l, y)_u(x, y), u(x, y+l)_u(x, y))，对对数图 像u进行前向差分操作，得到对数图像u的梯度，式中，u (x，y)为对数图像u的第x行第y 列的像素点灰度值； d. 设定矿物图像f中过暗区域和过亮区域的判别标准：设定过暗区域像素灰度临界值 为T1，过亮区域像素灰度临界值为T 2，由于8位灰度数字图像中，像素灰度值低于50的区 域或像素灰度值大于200的区域，其图像细节人眼均无法识别，所以选取默认值T 1= 50,选 取默认值T2= 200 ;增加矿物图像f的高斯噪声，由于人眼只能识别矿物图像f的高斯噪声 标准方差在4到10之间，因此选取矿物图像f的梯度倍乘因子α = 10/4 = 2. 5 ; e. 取矿物图像f中的每个像素点（x，y)，由该像素点及与该像素上下左右邻接的 四个像素点构成像素点（X，y)的一个子集：N (X，y) = {(X，y)，（x+1, y)，（X-1, y)，（X， y+1)，（x，y-l)}，若该子集中的每个像素点的灰度值都小于像素灰度临界默认值T1，即 {f(l，k) ST1, (l，k) eN(x，y)}，则判断该子集区域为矿物图像过暗区域，若该子集中的每 个像素点的灰度值都大于像素灰度临界默认值T2，即{f(l，k)多T 2,（l，k) eN(x，y)}，则 判断该子集区域为矿物图像过亮区域，其余像素点子集不做处理； f. 选取过暗区域和过亮区域的子集，利用公式：V= a Vu(x，y) = (aux, auy)，计算 过暗区域和过亮区域对比度改善之后的梯度v，分别将过暗区域和过亮区域的对数图像梯 度Vu( X，y)乘以梯度倍乘因子a，使乘以梯度倍乘因子a后的矿物图像过暗区域和过亮 区域的对比度提高，拉伸过暗区域和过亮区域的灰度动态范围，在梯度域提高细节，增强了 矿物图像过暗区域和过亮区域的纹理细节，改善视觉效果，得到调节对比度后的矿物对数 图像的梯度V ; g. 经过步骤f中梯度倍乘因子a调节后，得到调节对比度后的矿物对数图像的梯度V 为：
式中Vx为梯度V的横坐标，V y 为梯度V的纵坐标； h. 根据纽曼边界条件，将对数图像u的像素定义域进行边缘拓展，即在对数图像u的原 始定义域中添加首行、末行、首列、末列，同时将调节对比度后的矿物对数图像的梯度V的 边缘像素点作镜像对称，得到梯度域拓展的调节对比度后的矿物对数图像的梯度V 1; i. 矿物对数图像的梯度'为： a ▽ u (X，y) = ( a ux, a uy)，过暗和过亮区域 V1= {对应的边缘像素点梯度值，边缘镜像像素点} = (V lx，Vly)，其中Vlx为梯度V满 横坐标，VlyV u (x，y) = (ux, uy),其他区域 为梯度V1的纵坐标； j. 分别计算Vlx，Vly的离散傅里叶变换 k. 利用离散傅里叶变换求解泊松方程Au1= divV1:首先利用公式：
求解得到已调节矿物对数图像各个像素点灰度值的离散傅 里叶变换再利用公式：
求解每个像素点离散傅里叶变换 的"",的离散傅里叶反变换Uljl，得到已调节矿物对数图像U1，式中：j =0,……，J_1 ;1 = 〇,......，L-I ;J为步骤h中拓展定义域后的对数图像u的像素点的行数，L为像素点的的 列数，i为虚数单位，将已调节矿物对数图像从梯度域恢复到空间域，所述Uljl为已调节矿 物对数图像^第j行第1列这个像素点的灰度值；
1. 将已调节矿物对数图像^中通过步骤h拓展出的首行、末行、首列、末列像素点裁剪 掉，得到原始定义域大小的矿物对数图像U2; m.通过公式：/' -1，将调节对比度后的矿物对数图像112中的每个像素点灰度值 取指数操作，得到最终对比度调节后的矿物图像f'。
2. 根据权利要求1所述的基于泊松方程的矿物图像对比度调节方法，其特征在于:所述步 骤j中计算Vlx，V1,离散傅里变换得的公式分别为
式中m = 0,--,J-I ；n = 0,--，L-I ; J为步骤h中梯度域拓展 后V1的像素点的行数，L为像素点的的列数。
【专利摘要】一种基于泊松方程的矿物图像对比度调节方法，适用于对矿物品质研究使用。其步骤如下：对矿物图像进行去噪处理，在梯度域中，对去噪后的矿物图像进行变换，利用离散傅里叶变换求解泊松方程的方法，将调节对比度后图像从梯度域恢复到空间域，完成图像对比度调节。本发明通过上述处理步骤，有效改善了矿物图像的对比度，增强了图像纹理细节，有利于对矿物图像进行后续分析。
【IPC分类】G06T5-20
【公开号】CN104766288
【申请号】CN201510195438
【发明人】卢兆林, 陈燕, 李世银
【申请人】中国矿业大学
【公开日】2015年7月8日
【申请日】2015年4月22日