1.一种基于泊松方程的矿物图像对比度调节方法,其特征在于包括以下步骤:a.使用扫描电子显微镜获取对被测矿物的扫描图像,并采用高斯滤波器对获取到的扫描图像进行去噪处理,从而得到去噪后的矿物图像f;
b.使用公式:u=log(f+1),将去噪后的矿物图像f进行取对数操作,得到对数图像u;
c.利用公式: 对对数图像u进
行前向差分操作,得到对数图像u的梯度,式中,u(x,y)为对数图像u的第x行第y列的像素点灰度值;
d.设定矿物图像f中过暗区域和过亮区域的判别标准:设定过暗区域像素灰度临界值为T1,过亮区域像素灰度临界值为T2,由于8位灰度数字图像中,像素灰度值低于50的区域或像素灰度值大于200的区域,其图像细节人眼均无法识别,所以选取默认值T1=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-
1)},若该子集中的每个像素点的灰度值都小于像素灰度临界默认值T1,即{f(l,k)≤T1,(l,k)∈N(x,y)},则判断该子集区域为矿物图像过暗区域,若该子集中的每个像素点的灰度值都大于像素灰度临界默认值T2,即{f(l,k)≥T2,(l,k)∈N(x,y)},则判断该子集区域为矿物图像过亮区域,其余像素点子集不做处理;
f.选取过暗区域和过亮区域的子集,利用公式: 计算过暗区
域和过亮区域对比度改善之后的梯度v,分别将过暗区域和过亮区域的对数图像梯度乘以梯度倍乘因子α,使乘以梯度倍乘因子α后的矿物图像过暗区域和过亮区域的对比度提高,拉伸过暗区域和过亮区域的灰度动态范围,在梯度域提高细节,增强了矿物图像过暗区域和过亮区域的纹理细节,改善视觉效果,得到调节对比度后的矿物对数图像的梯度V;
g.经过步骤f中梯度倍乘因子α调节后,得到调节对比度后的矿物对数图像的梯度V为:式中Vx为梯度V的横坐标,Vy为梯度
V的纵坐标;
h.根据纽曼边界条件,将对数图像u的像素定义域进行边缘拓展,即在对数图像u的原始定义域中添加首行、末行、首列、末列,同时将调节对比度后的矿物对数图像的梯度V的边缘像素点作镜像对称,得到梯度域拓展的调节对比度后的矿物对数图像的梯度V1;
i.矿物对数图像的梯度V1为:
其中V1x为梯度V1的横坐标,V1y为梯度V1的纵坐标;
j.分别计算V1x,V1y的离散傅里叶变换
k.利用离散傅里叶变换求解泊松方程Δu1=divV1:首先利用公式:求解得到已调节矿物对数图像各个像素点灰度值的离散傅里叶变换 再利用公式:求解每个像素点离散傅里叶变换 的离散傅里叶反变换u1jl,
得到已调节矿物对数图像u1,式中:j=0,......,J-1;l=0,......,L-1;J为步骤h中拓展定义域后的对数图像u的像素点的行数,L为像素点的的列数,i为虚数单位,将已调节矿物对数图像从梯度域恢复到空间域,所述u1jl为已调节矿物对数图像u1第j行第l列这个像素点的灰度值;
l.将已调节矿物对数图像u1中通过步骤h拓展出的首行、末行、首列、末列像素点裁剪掉,得到原始定义域大小的矿物对数图像u2;
m.通过公式: 将调节对比度后的矿物对数图像u2中的每个像素点灰度值取指数操作,得到最终对比度调节后的矿物图像f'。
2.根据权利要求1所述的基于泊松方程的矿物图像对比度调节方法,其特征在于:所述步骤j中计算V1x,V1y的离散傅里变换得 的公式分别为:式中m=0,....,J-1;n=0,....,L-1;J为步骤h中梯度域拓展后V1的像素点的行数,L为像素点的的列数。