[发明专利]一种基于扩展二维经验小波变换的工程表面多尺度滤波方法

说明书

本发明涉及图像处理技术领域，针对现有滤波方法中拟合误差、末端效应和模态混叠的问题，公开了一种基于扩展二维经验小波变换的工程表面多尺度滤波方法，所述方法包括下述步骤：（1）生成灰度图像；（2）获得去全局趋势的频谱；（3）自适应频谱分割；（4）二维经验小波分解；（5）经验模态函数合并；（6）表面成分重构。本发明基于自适应尺度空间表示技术和Riesz变换，扩展了原始的二维经验小波变换，建立了每个提取的固有模态函数与其物理尺寸之间的关系，实现工程表面多尺度表征和滤波，从而指示制造过程和零件功能的不规则性，为后续建立表面成分和加工工艺系统、零件性能的映射关系等工作的开展提供了基础。

技术领域

本发明涉及图像处理技术领域，尤其是涉及一种基于扩展二维经验小波变换的工程表面多尺度滤波方法。

背景技术

表面形貌是制造过程中零件与其预期功能行为之间的重要联系。为了满足更严格的公差要求和更高的性能标准，工程师往往需要零件表面质量满足指定的尺寸、形状和峰值等各种表面功能属性。工程表面的表面形貌存在多尺度性，由三个波长范围的成分组成，其中较长的波长成分是形状误差，而较短的波长成分是粗糙度，波纹度是中等波长成分，介于形状误差和粗糙度之间。不同尺度的表面成分具有不同的来源，并且以不同的方式影响零件的功能，如振动、刀具磨损、颤振和密封。因此，实现表面的多尺度滤波对工程表面形貌分析具有重要意义。

滤波技术是将表面分离成形状误差、波纹度和粗糙度的过程，也一直是表面形貌分析方面的研究热点。2RC滤波器和高斯滤波器是最早提出的滤波技术，且高斯滤波器仍被视为当前的标准滤波器。作为非线性和非平稳信号处理工具，小波变换和经验模态分解逐渐用于工程表面分析。与传统的滤波方法不同，它们具有多尺度几何结构，可以将给定的表面形貌划分为一系列频率分量，每个分量都与其尺度相匹配，分离的频率分量的数量多于三个，能更好地反映更精细的表面信息。已有技术中，Jerome Gilles在论文“EmpiricalWavelet Transform”(IEEE Transactions on Signal Processing.2013；61:3999-4010)中提出经验小波变换 (EWT)方法，它继承了小波变换和经验模态分解的滤波特性，不仅克服了小波滤波在小波基上的非自适应性，而且填补了经验模态分解的理论基础，避免了经验模态分解滤波的模态混叠，但EWT缺乏对表面各尺度成分的表征。

近年来出现的高清晰测量技术综合了接触式测量技术测量范围大和非接触式测量技术采样密度高的优点，能够为零件加工表面预测复现更为清晰的表面形貌。专利号为CN103544694A名称为“零件表面形貌三维高密度点云数据转化为灰度图像的方法”的中国专利，提出了把测量点云转化为灰度图像的方法，为后续表面多尺度滤波提供了新的数据基础。

发明内容

为了克服传统滤波方法中拟合误差、末端效应和模态混叠的缺点，弥补EWT无法表面成分表征的不足，本发明提出了一种基于扩展二维经验小波变换(Extended Bi-dimensional Empirical Wavelet Transform，EBEWT)的工程表面多尺度滤波方法。本发明基于自适应尺度空间表示技术和Riesz变换，扩展了原始的二维经验小波变换，建立了每个提取的固有模态函数与其物理尺寸之间的关系，实现工程表面各尺度成分的表征和滤波，从而指示制造过程和零件功能的不规则性。

为了实现上述目的，本发明采用以下技术方案：

一种基于扩展二维经验小波变换的工程表面多尺度滤波方法，所述方法包括下述步骤：

(1)生成灰度图像：对测量后的零件表面形貌三维高密度点云数据进行处理，生成能够反映表面形貌高度值的灰度图像；

(2)获得去全局趋势的频谱：通过伪极快速傅里叶变换和多项式拟合得到；

(3)自适应频谱分割：通过尺度空间表示技术生成频谱的尺度空间平面，并找出全部有意义的模态；

(4)二维经验小波分解：建立二维小波滤波器组，得到经验模态函数；

(5)经验模态函数合并：将波长符合特定表面成分的经验模态函数进行合并；