照度不均匀电子散斑干涉条纹图像增强方法与流程

本发明属于光信息处理技术和图像处理技术领域，涉及一种基于retinex理论的照度不均匀的电子散斑干涉条纹图像增强方法。

背景技术：

电子散斑干涉(electronicspecklepatterninterferometry,简称espi)测量技术，是一种现代光学测量技术，它具有全场非接触、高精度和高灵敏度、不避光等优点。该技术在机械、土木、水利、电器、航空航天、兵器工业及生物医学等领域的检测中具有非常重要的地位。电子散斑干涉技术利用光的干涉原理对被测物体进行精密测量，当被测的工件产生微小的位移时，在其表面上所形成的散斑图形也会发生移动，利用光学照相的方法将散斑图形记录下来，对变形前后采集的两幅散斑图像做绝对值相减、相加或相乘的运算，即可得到对应变形信息的散斑条纹图，因此，准确地从干涉条纹图像中提取相位是获得物体变形位移信息的关键。

电子散斑干涉技术虽然有很多优良特性，可是随着现代科学技术的飞速发展，从工程及应用科学方面提出的测试问题越来越多，产生的条纹图也越来越复杂，尽管目前光条纹图像信息提取方法取得了很大进展，但是在实验获取的电子散斑条纹图像中有时会出现照度不均匀的现象，现有的方法在对照度不均匀的条纹图进行处理时无法获得的完整的条纹信息，进而无法得到完整的相位信息，这是光学测试领域的一个难点。因此，找到一种能够有效解决电子散斑条纹图像的照度不均匀的处理方法，对于后续处理显得尤为重要。

近些年基于人类视觉的图像增强技术成为一种研究趋势，其中，retinex理论是通过模拟人眼成像过程而建立起来的一种颜色恒常知觉计算理论。retinex这个词是由视网膜(retina)和大脑皮层(cortex)两个词组合构成，其理论是land和mccann提出的，它的基本思想是物体的颜色由物体对长波、中波和短波光线的反射能力决定的，并不取决于该点进入人眼的绝对光线，并且物体的色彩不受光照非均匀性的影响，具有一致性。retinex可以分为两种形式：基于路径的retinex以及基于中心环绕retinex。常用的retinex算法是基于中心环绕的单尺度retinex算法(single-scaleretinex,ssr)和多尺度retinex算法(multi-scaleretinex,msr)。retinex算法基本原理是将给定图像分解成两个不同的图像，反射图像和入射图像，反射图像就是原始没有受到干扰的图像。该方法是要在给定图像中，降低入射图像的影响，尽量保留物体本质的反射属性，从而得到图像最本质的样子。因此可以将retinex算法引入到电子散斑干涉条纹图像中，针对照度不均匀这一现象，结合retinex理论对图像光线的处理，可以去除光线对图像的影响，得到照度均匀的结果。

技术实现要素：

为克服现有技术的不足，本发明旨在实现对照度不均匀的电子散斑干涉图像的对比度增强，图像滤波后的条纹能够尽量保持完整，使得后续提取骨架线时提取的更为精准。本发明采用的技术方案是，利用retinex理论，选取基于中心环绕的多尺度retinex算法msr(multi-scaleretinex)，对低质量且照度不均匀的电子散斑干涉条纹图像进行处理，调节尺度参数，实现图像增强。

照度不均匀电子散斑干涉条纹图像增强方法，具体步骤如下：

步骤1：获得低质量、照度不均匀电子散斑干涉条纹图像f；

步骤2：构建retinex理论模型；

步骤3：针对每幅图调节msr算法的各个尺度参数，使得到的增强图像i1的照度均匀，效果最佳。

步骤2：构建retinex理论模型，具体步骤如下：

步骤2-1：首先，根据retinex理论，利用取对数的方法将照射分量和反射分量分离，即：

logs(x，y)＝logr(x，y)+logl(x，y)

其中，s代表我们所看到的图像，l是入射光图像，r是物体的反射性质。retinex就是通过图像s来得到物体的反射性质r，也就是设法去除或者降低入射光l的影响从而得到物体原本该有的样子；

步骤2-2：通过计算图像中像素点与周围区域中像素的加权平均来对图像中照度变化做估计，并将其去除，最后只保留图像中物体的反射属性，用高斯模板对原图像作卷积，即相当于对原图像作低通滤波，得到低通滤波后的图像d(x，y)，f(x，y)表示高斯滤波函数：

d(x，y)＝s(x，y)*f(x，y)

步骤2-3：在对数域中，用原图像减去低通滤波后的图像，得到高频增强的图像g(x，y)：

g(x，y)＝logs(x，y)-logd(x，y)

步骤2-4：对g(x，y)取反对数，得到增强后的图像r(x，y)：

r(x，y)＝expg(x，y)

步骤3：根据msr尺度的要求，分别用大、中、小三种尺度对图像分别进行步骤2-2，步骤2-3及步骤2-4，并且对每个尺度下的结果进行加权求和。

还包括效果比较步骤：

步骤1：分别利用直方图均衡化、常微分方程增强对输入的图像进行增强，得到增强图像i2、i3；

步骤5：将得到的所有增强图像与原图像进行滤波，得到滤波后的结果并进行对比。

本发明的特点及有益效果是：

本发明提出的基于retinex理论的照度不均匀的电子散斑干涉条纹图像增强方法，相比于传统的电子散斑干涉条纹图像增强方法，增强效果更好，几乎去除了照度不均匀给图像带来的问题，尽可能多的保持了条纹的边缘，得到的条纹更加完整，更利于后续的骨架线提取和相位展开。

附图说明：

图1是本发明基于retinex理论的照度不均匀电子散斑干涉条纹图像增强方法的方法流程图；

图2(a)(b)是低密度电子散斑干涉条纹图像及其滤波结果，(c)(d)是对图2(a)的直方图均衡化的增强图像及其滤波结果，(e)(f)是对图2(a)的常微分方程增强图像及其滤波结果，(g)(h)是对图2(a)的msr增强图像及其滤波结果；

图3(a)(b)是高密度电子散斑干涉条纹图像及其滤波结果，(c)(d)是对图3(a)的直方图均衡化的增强图像及其滤波结果，(e)(f)是对图3(a)的常微分方程增强图像及其滤波结果，(g)(h)是对图3(a)的msr增强图像及其滤波结果；

图4(a)(b)是变密度电子散斑干涉条纹图像及其滤波结果，(c)(d)是对图4(a)的直方图均衡化的增强图像及其滤波结果，(e)(f)是对图4(a)的常微分方程增强图像及其滤波结果，(g)(h)是对图4(a)的msr增强图像及其滤波结果；

具体实施方式

接下来应用本发明提出的基于retinex理论的照度不均匀电子散斑干涉条纹图像增强方法进行增强。具体步骤如下：

步骤1：进行实验获得照度不均匀的实验图像。根据电子散斑干涉条纹图像的密度特征进行实验，实验获得低密度和高密度的照度不均匀电子散斑干涉条纹图像，变密度的照度不均匀电子散斑条纹图4(a)是经过电脑模拟得到的。

步骤2：构建retinex理论模型。具体步骤如下：

步骤2-1：首先，根据retinex理论，利用取对数的方法将照射分量和反射分量分离，即：

logs(x，y)＝logr(x，y)+logl(x，y)

其中，s代表我们所看到的图像，l是入射光图像，r是物体的反射性质。retinex就是通过图像s来得到物体的反射性质r，也就是设法去除(或者降低)入射光l的影响从而得到物体原本该有的样子。

步骤2-2：最终的反射图像会被假设地估计为空间平滑图像(其物理解释就是通过计算图像中像素点与周围区域中像素的加权平均来对图像中照度变化做估计，并将其去除，最后只保留图像中物体的反射属性)，所以可以用高斯模板对原图像作卷积，即相当于对原图像作低通滤波，得到低通滤波后的图像d(x，y)，f(x，y)表示高斯滤波函数：

d(x，y)＝s(x，y)*f(x，y)

步骤2-3：在对数域中，用原图像减去低通滤波后的图像，得到高频增强的图像g(x，y)：

g(x，y)＝logs(x，y)-logd(x，y)

步骤2-4：对g(x，y)取反对数，得到增强后的图像r(x，y)：

r(x，y)＝expg(x，y)

步骤3：根据msr尺度的要求，尺度的大小选择最经典的三尺度，分别用大、中、小三种尺度对图像分别进行步骤2-2，步骤2-3及步骤2-4，并且对每个尺度下的结果进行加权求和。

步骤4：为了突出本发明所提出的方法的优越性，本次实验还选取了两个传统的图像增强方法，分别为直方图均衡化增强方法和常微分方程图像增强方法，对实验图像进行增强。

步骤5：将得到的所有增强图像与原图像进行滤波，得到滤波后的结果。

该算法的流程图见附图1。

为了验证方法的有效性，给出实验结果。

图2(a)为低密度电子散斑干涉条纹图像，图2(b)是图2(a)滤波之后的结果。同样的，图2(c)(e)(g)，分别为通过直方图均衡化、常微分方程及msr算法增强后的图像，而图2(d)(f)(h)分别为上述三幅图像的滤波结果。图3(a)为高密度的电子散斑干涉条纹图像，图3(b)是图3(a)滤波之后的结果。同样的，图3(c)(e)(g)，分别为通过直方图均衡化、常微分方程及msr算法增强后的图像，而图3(d)(f)(h)分别为上述三幅图像的滤波结果。图4(a)为变密度的电子散斑干涉条纹图像，图4(b)是图4(a)滤波之后的结果。同样的，图4(c)(e)(g)，分别为通过直方图均衡化、常微分方程及msr算法增强后的图像，而图4(d)(f)(h)分别为上述三幅图像的滤波结果。

无论是从增强图像的对比还是滤波之后的图像对比，都可以看出本发明提出的图像增强方法的能够很好的对各种密度特征的照度不均匀的电子散斑干涉条纹图像进行增强，并且滤波之后的结果也有优越性，几乎消除了照度不均匀的现象。

尽管上面结合图示对本发明进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨的情况下，还可以作出很多变形，这些均属于本发明的保护之内。

本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。