改进的基于小波分解的旋滤波算法制造技术

本发明专利技术公开了一种改进的基于小波分解的旋滤波算法，包括：先用小波进行分解，得到含有携带大部分图像信息的低频图像，然后对低频图像进行旋滤波，定义低频图像的条纹方向，求各个方向均值，计算该点各个方向上灰度值与该方向线的均值的差，取绝对值，累加求和。找出累加和最小的方向及切线的方向，在该方向上进行中值或均值滤波。本发明专利技术的有益结果，是在时间效率优于原始旋滤波的同时，提高了降噪效果，为下一步进行干涉图像的相位提取和相位解包裹提供了良好的前提。

【技术实现步骤摘要】

本专利技术涉及基于计算机的电子散斑干涉测量
，尤其涉及散斑干涉条纹图的降噪技术。
技术介绍
散斑在最初的研究中一直是被作为影响图像的噪声来进行去除的，直到1966年，Ennos发现它存在可以测量的强度和准确的相位，确定了散斑的应用价值。散斑不仅是一种噪声，而且是测量物体的信息载体。人们一般将用电子硬件处理的电子散斑干涉技术(ESPI)和数字软件处理的数字散斑干涉技术(DSPI)统称为电子散斑干涉技术。其基本原理是让激光散斑承载被测物体变化的信息，根据被测物体在激光照射后产生的散斑场的相关条纹测量光束前后的相位变化。它利用计算机软件处理得到的记录物体表面特征的变化信息，将散斑干涉条纹图在计算机上显示，并通过各种运算和处理把这些条纹转化为人们熟悉的物理量。用于抑制图像噪声的方法有空域的均值、中值和维纳滤波等以及变换域的傅立叶变换和小波变换等。传统的滤波技术在去除条纹图的噪声部分同时降低了条纹图的对比度使图像变得模糊，而且有时会破坏条纹的结构，不利于下一步提取相位进行测量的研究。现有技术中提出了在整幅图像中而不是像旋滤波在局部的窗口中寻找相似像素，滤波效果更加有效，但是该方法在全图中寻找的时间相对来说比较长，算法速度较慢。现有技术中提出了偏微分方程方法在散斑干涉条纹图滤波中的应用，将待处理图像看作偏微分算子，利用偏微分方程把初始图像变形，从而将偏微分方程与图像联系起来，取得了不错的结果，但是，该方法需要较多的参数，所需处理时间较长。目前国内外对于散斑条纹图去噪的研究分析，并没有一种方法能够适应所有的散斑干涉条纹图。由于条纹图的高噪声和条纹结构的不同，怎样在去噪的同时保持条纹图细节信息，减少图像的变形，提高条纹图的对比度，提高计算效率。
为我们下一步提取相位和相位的解包裹提供一个良好的条件。
技术实现思路
本专利技术为了解决上述问题，提出了一种改进的基于小波分解的旋滤波算法，首先对即将处理的图像进行小波分解，然后对低频图像进行旋滤波，在提高算法效率的同时保持了图像的降噪效果。为实现上述目的，本专利技术采用下述技术方案：一种改进的基于小波分解的旋滤波算法，包括以下步骤：(1)采集待处理的含有噪声的图像，用Haar小波进行小波分解，得到图像分解后的低频散斑干涉条纹图；(2)定义得到的低频散斑干涉条纹图中，当前点为中心的设定窗口内以等角度间隔的条纹方向滤波线；(3)计算每一条纹方向滤波线上各点的灰度平均值；(4)根据条纹方向滤波线上各点的灰度值与该条纹方向滤波线上各点的灰度平均值，计算每一条纹方向滤波线的灰度变化；(5)选取灰度变化最小的方向滤波线作为条纹的切线方向；(6)在条纹的切线方向进行中值或者均值滤波；(7)对低频散斑干涉条纹图中的每一个点重复步骤(2)-(6)的过程，得到降噪以后的散斑干涉条纹图。进一步地，所述步骤(1)中，采用Haar小波进行单尺度的二维离散小波变换，得到一幅低频图和三幅不同方向上的高频图。进一步地，所述步骤(2)中，在图像当前点(i,j)为中心的n×n像素内，定义间隔相等的方向滤波直线；所述条纹方向滤波线取8条、16条或者32条。进一步地，所述步骤(4)中，每一条纹方向滤波线的灰度变化具体为：条纹方向滤波线上每一点的灰度值与该条纹方向滤波线上各点的灰度平均值之差的绝对值的累加和。进一步地，所述步骤(6)中，中值滤波是将所述条纹切线方向上各点的灰度值进行排序，用中间点的值取代当前点的灰度值。进一步地，所述步骤(6)中，均值滤波是将所述条纹切线方向上各点的灰度平均值取代当前点灰度值。本专利技术的有益效果是：通过对真实散斑干涉条纹图和模拟散斑干涉条纹图进行降噪处理，本专利技术新提出的改进的算法能提高降噪效率的同时保持散斑干涉条纹图的降噪效果。适于实际应用，为散斑干涉条纹图下一步的提取相位和相位的解包裹提供良好的前提条件。附图说明图1(a)-图1(b)为图像一层小波分解和图像两层小波分解；图2(a)-图2(f)为旋滤波方向滤波线示意图；图3(a)-图3(c)为模拟生成的无噪声条纹图像、添加噪声以后的条纹图、真实的散斑干涉条纹图；图4(a)-图4(c)为上述图3(a)-图3(c)三幅图的剖面图；图5(a)-图5(d)为模拟噪声图像旋滤波降噪后图像及剖面图、真实散斑干涉条纹图像旋滤波降噪后图像及剖面图；图6(a)-图6(d)分别为模拟噪声图像新改进的旋滤波降噪后图像及其剖面图、真实散斑干涉条纹图像新改进旋滤波降噪后图像及其剖面图；图7为旋滤波与本文改进旋滤波的时间对比示意图；图8为本专利技术的整体流程。具体实施方式下面结合附图与实施例对本专利技术作进一步说明。首先模拟生成一幅条纹图，考虑到散斑干涉条纹图中含有较多的乘性噪声和少量的加性噪声，在生成的模拟条纹图添加一定的乘性散斑噪声和少量的高斯加性噪声，使之与真实的散斑干涉条纹图相近。选取两幅真实散斑干涉条纹图为下面实验做准备。小波分解我们采用的是Haar小波。在众多的小波函数中，Haar小波与其他的正交函数相比，它具有构造简单，对应的滤波器具有线性相位性、计算方便等特点。Haar小波的定义如下： Ψ ( Ω ) = j 4 Ω sin 2 ( Ω a ) e - j Ω 2 / 2 ; ]]>图1(a)-图1(b)为图像一层小波分解和图像两层小波分解，图像经过一层分解以后被分成一个低频子带和三个高频子带，第二层将低频子带分成四个子带。图2(a)-图2(f)为旋滤波及其方向滤波线示意图。图3(a)-图3(c)为模拟生成的无噪声条纹图像、添加噪声以后的模拟散斑条纹图、真实的散斑干涉条纹图，仿真模拟图是为了更好的对比试验效果。图4(a)-图4(c)为上述图3(a)-图3(c)三幅图的剖面图，剖面图可以明显的看出降噪效果和条纹结构的保持，提取剖面图可以更好的对比降噪效果。图5(a)-图5(d)为模拟噪声图像旋滤波降噪后图像及剖面图、真实散斑干涉条纹图像旋滤波降噪后图像及剖面图，可以看出旋滤波降噪以后图像的剖面图5(b)和5(d)明显优于未降噪图像的剖面图4(b)和4(c)。图6(a)-图6(d)分别为模拟噪声图像新改进的旋滤波降噪后图像及其剖面图、真实散斑干涉条纹图像新改进的旋滤波降噪后图像及其剖面图，新改进旋滤波降噪后图像剖面图6(b)和6(d)效果优于旋滤波降噪以后图像的剖面图5(b)和5(d)，降噪效果提高。图7选取3组大小不同的图像进行旋滤波与改进旋滤波降噪时间比较图，从图7中可以看出新改进的旋滤波降噪时间变短，而且图像越大，新改进的旋滤波效率提高越明显。如图8所示，本专利技术所述改进的降噪算法的具体过程如下：(1)输入一幅含有噪声的图像，用Haar小波进行小波分解，得到一幅低频图和三幅不同方向上的高频图。(2)定义得到的低频散斑干涉条本文档来自技高网...

【技术保护点】
一种改进的基于小波分解的旋滤波算法，其特征是，包括以下步骤：(1)采集待处理的含有噪声的图像，用Haar小波进行小波分解，得到图像分解后的低频散斑干涉条纹图；(2)定义得到的低频散斑干涉条纹图中，当前点为中心的设定窗口内以等角度间隔的条纹方向滤波线；(3)计算每一条纹方向滤波线上各点的灰度平均值；(4)根据条纹方向滤波线上各点的灰度值与该条纹方向滤波线上各点的灰度平均值，计算每一条纹方向滤波线的灰度变化；(5)选取灰度变化最小的方向滤波线作为条纹的切线方向；(6)在条纹的切线方向进行中值或者均值滤波；(7)对低频散斑干涉条纹图中的每一个点重复步骤(2)‑(6)的过程，得到降噪以后的散斑干涉条纹图。

【技术特征摘要】
1.一种改进的基于小波分解的旋滤波算法，其特征是，包括以下步骤：(1)采集待处理的含有噪声的图像，用Haar小波进行小波分解，得到图像分解后的低频散斑干涉条纹图；(2)定义得到的低频散斑干涉条纹图中，当前点为中心的设定窗口内以等角度间隔的条纹方向滤波线；(3)计算每一条纹方向滤波线上各点的灰度平均值；(4)根据条纹方向滤波线上各点的灰度值与该条纹方向滤波线上各点的灰度平均值，计算每一条纹方向滤波线的灰度变化；(5)选取灰度变化最小的方向滤波线作为条纹的切线方向；(6)在条纹的切线方向进行中值或者均值滤波；(7)对低频散斑干涉条纹图中的每一个点重复步骤(2)-(6)的过程，得到降噪以后的散斑干涉条纹图。2.如权利要求1所述的一种改进的基于小波分解的旋滤波算法，其特征是，所述步骤(1)中，采用Haar小波进行单尺度的二维离散小波变换，得到一幅低频图和...

【专利技术属性】
技术研发人员：辛化梅，王喜连，
申请(专利权)人：山东师范大学，
类型：发明
国别省市：山东;37