CN110227968A - 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 - Google Patents
基于小波变换的磁流变去除函数寄生条纹消除方法及装置 Download PDFInfo
- Publication number
- CN110227968A CN110227968A CN201910410490.5A CN201910410490A CN110227968A CN 110227968 A CN110227968 A CN 110227968A CN 201910410490 A CN201910410490 A CN 201910410490A CN 110227968 A CN110227968 A CN 110227968A
- Authority
- CN
- China
- Prior art keywords
- striped
- parasitic
- removal function
- coefficient matrix
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000024241 parasitism Effects 0.000 title claims abstract description 34
- 230000009466 transformation Effects 0.000 title claims abstract description 19
- 230000003071 parasitic effect Effects 0.000 claims abstract description 121
- 230000008030 elimination Effects 0.000 claims abstract description 9
- 238000003379 elimination reaction Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 109
- 238000000354 decomposition reaction Methods 0.000 claims description 27
- 230000000737 periodic effect Effects 0.000 claims description 21
- 101100127891 Caenorhabditis elegans let-4 gene Proteins 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 abstract description 7
- 238000005498 polishing Methods 0.000 abstract description 7
- 238000001228 spectrum Methods 0.000 abstract description 7
- 238000001914 filtration Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000010884 ion-beam technique Methods 0.000 abstract description 3
- 238000007514 turning Methods 0.000 abstract description 3
- 239000000758 substrate Substances 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000005350 fused silica glass Substances 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 229940099259 vaseline Drugs 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B1/00—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes
- B24B1/005—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes using a magnetic polishing agent
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B13/00—Machines or devices designed for grinding or polishing optical surfaces on lenses or surfaces of similar shape on other work; Accessories therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Theoretical Computer Science (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了基于小波变换的磁流变去除函数寄生条纹消除方法及装置,该方法能够准确去除寄生条纹,达到磁流变去除函数的精确提取,适用范围广,对于传统频域带阻滤波无法解决的寄生条纹和原始数据频谱严重混叠的情况,仍然适用,是一种较好的寄生条纹去除方法,可以有效应用于消除磁流变去除函数寄生条纹,实现磁流变去除函数的准确提取。同时也可以较好的用于其它超精密加工方式,如气囊抛光、离子束抛光等的去除函数寄生条纹的消除,便于推广应用。
Description
技术领域
本发明涉及光学加工领域,具体是基于小波变换的磁流变去除函数寄生条纹消除方法及装置。
背景技术
磁流变抛光技术作为一种新型的先进光学制造技术,具有高效率、高精度、高表面质量、亚表面损伤小、表面残余应力小等一系列传统加工方法不可比拟的优点,具有良好的应用前景。磁流变去除函数提取的准确程度是决定磁流变加工精度与效率的一个关键因素。目前磁流变去除函数制作流程通常是首先采用平板小基片在磁流变液机床进行采斑,然后用光学干涉仪检测采斑基片的反射或透射面形,然后将检测数据中的去除函数进行提取,最终生成可用于磁流变抛光的去除函数。
由于采斑所用的小基片是平板型光学元件,采用干涉仪进行面形检测时,极易引入寄生条纹,引入寄生条纹后会降低去除函数的提取精度,造成去除函数体积去除效率变化,影响磁流变加工的确定性,同时引入了寄生条纹的去除函数含有复杂的中高频结构,严重影响对去除函数修形能力的评价,因此为了实现磁流变去除函数准确提取必须要消除去除函数中的寄生条纹。
传统抑制干涉检测中寄生条纹的方法是反射测量时在元件背面均匀涂抹凡士林,难以消除寄生条纹;透射测量时将后标准镜与将待测元件形成一定楔角,但是由于做斑基片口径小,通常为Φ50mm,能够允许的楔角很小,难以抑制寄生条纹;采用短相干干涉测量的方法可以消除寄生条纹,如采用4D PhaseCam2000型干涉仪,但是由于其动态测量时对比度较低,倍率放大导致对比度太小而无法测量,使得基片测量时只能采用低倍率,导致检测面形有效数据点很少,难以准确提取去除函数。
磁流变去除函数通常几何尺寸宽度在2~15mm、高度在3~25mm之间,而采斑基片通常采用小口径的平板基片,用干涉仪进行波前检测时发现寄生条纹频谱与磁流变去除函数的频谱往往有存在明显混叠,采用传统频域带阻滤波去除寄生条纹时,往往会明显改变去除函数的形态,使得最终获得的去除函数失真。
因此,如何消除磁流变去除函数中的寄生条纹,提高磁流变去除函数的提取精度和去除函数的修形能力是同行从业人员亟待解决的问题。
发明内容
鉴于上述问题,本发明的目的是为了解决磁流变去除函数干涉仪检测中引入寄生条纹,降低去除函数提取精度和干扰去除函数修形能力评价的问题。
本发明实施例提供基于小波变换的磁流变去除函数寄生条纹消除方法,包括:
获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
在一个实施例中,还包括:
对处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
在一个实施例中,所述N根据所述空间周期T计算获得,包括:
计算多级小波分解层数N,计算公式如下:
其中,N取满足条件下最小整数,T表示寄生条纹的空间周期,pix表示像素尺寸。
在一个实施例中,对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征,包括:
根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
第二方面,本发明还提供基于小波变换的磁流变去除函数寄生条纹消除装置,包括:
获取模块,用于获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
计算模块,用于提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
分解模块,用于采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
记录模块,用于分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
消除模块,用于对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
在一个实施例中,所述装置还包括:
重构模块,用于对消除模块处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
在一个实施例中,所述计算模块中,具体用于计算多级小波分解层数N,计算公式如下:
其中,N取满足条件下最小整数,T表示寄生条纹的空间周期,pix表示像素尺寸。
在一个实施例中,所述消除模块,包括:
确定子模块,用于根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
消除子模块,用于将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
本发明实施例提供的上述技术方案的有益效果至少包括:
本发明实施例提供的一种基于小波变换的磁流变去除函数寄生条纹消除方法,能够准确去除寄生条纹,达到提取磁流变去除函数的精确提取,适用范围广,对于传统频域带阻滤波无法解决的寄生条纹和原始数据频谱严重混叠的情况,仍然适用,是一种较好的寄生条纹去除方法,可以有效应用于消除磁流变去除函数寄生条纹,实现磁流变去除函数的准确提取。同时也可以较好的用于其它超精密加工方式(如气囊抛光、离子束抛光等)的去除函数寄生条纹的消除,便于推广应用。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步地详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明实施例提供的基于小波变换的寄生条纹消除方法的流程图;
图2为本发明实施例提供的步骤S105的流程图;
图3为本发明实施例1提供的基于小波变换的寄生条纹消除方法的具体流程图;
图4为含有寄生条纹的磁流变去除函数检测数据图像;
图5为采用sym5小波4级分解,得到13个系数矩阵图像;
图6为S2系数矩阵消除寄生条纹流程示意图;
图7为滤除寄生条纹后的去除函数数据;
图8为本发明实施例提供的基于小波变换的寄生条纹消除装置的框图;
图9为本发明实施例提供的预消除模块85的框图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
参照图1所示,本发明实施例提供的一种基于小波变换的磁流变去除函数寄生条纹消除方法包括:步骤S101~S105;
S101、获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
S102、提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
S103、采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
S104、分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
S105、对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
本实施例中,能够准确去除寄生条纹,达到提取磁流变去除函数精确提取的要求,适用范围广,可以有效应用于磁流变去除函数寄生条纹的消除,同时也可以较好的用于其它超精密加工方式(如气囊抛光、离子束抛光等)的去除函数寄生条纹的消除,便于推广应用。
进一步地,参照图1所示,该方法还包括:
S106、对处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
在一个实施例中,参照图2所示,上述步骤S105中,对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征,包括:
S1051、根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
S1052、将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
本实施例中,通过对含有寄生条纹的系数矩阵Sj进行处理,实现了系数矩阵中寄生条纹和去除函数特征的有效分离,保证了消除寄生条纹的同时不改变去除函数自身形态,提高了数据处理精度。
具体的,本发明提供的基于小波变换的磁流变去除函数寄生条纹消除方法,包括如下步骤:
1、获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y),其中数据宽度像素数为m,数据高度像素数为n,像素尺寸为pix mm;
2、提取面形数据矩阵f1(x,y)中不含有去除函数形状但是含有寄生条纹的局部数据,观察寄生条纹正截面轮廓结构,计算寄生条纹的空间周期T;
3、计算多级小波分解层数N,计算公式如下
其中N取满足条件下最接近整数;
4、根据磁流变去除函数包含有大梯度结构面形的特点,采用Symlets小波,Symlets小波通常表示为symN(其中N=2~10),根据去除函数结构特点,选用sym5作为小波基函数。采用sym5小波对面形数据矩阵f1(x,y)进行多级分解,分解级数为N,每一级有三个系数矩阵,命名为hi(水平方向),vi(垂直方向),di(对角线方向),其中i=1~N,最后一级还有低频系数矩阵命名为AN,其中N表示分解层数,共得到3N+1个系数矩阵数据。上述Symlets小波是Inrid Daubechies提出的一种正交紧支撑小波函数,一种建立在多分辨率分析和多采样率滤波器理论基础上的离散序列小波变换,相比于经典的db小波,具有更好的对称性,同时其也具有正交性、紧支撑特征,在细节提取上更为优越。
5、分析系数矩阵hi,vi,di(其中i=1~N)的特征,记录含有寄生条纹的系数矩阵为Sj,其中j=1~K,K为含有寄生条纹的系数矩阵的个数。
6、对含有寄生条纹的系数矩阵Sj(其中j=1~K)进行处理,去除其中的寄生条纹特征,根据寄生条纹特征和去除函数特征数值不连续、存在强度差异的特点,采用基于直方图统计的方法确定寄生条纹特征所在取值范围,然后将系数矩阵中数值在寄生条纹特征所在取值范围内的全部数据置零,具体的实现过程为:
首先,分析系数矩阵Sj(其中j=1~K),搜索去除函数的特征分布所在区域,为Rp(p=1~H,H表示系数矩阵Sj中去除函数的个数),对区域Rp采用直方图统计的方法绘制曲线,横坐标为数据点的数值,纵坐标为数据点的个数;
然后,对于系数矩阵Sj中区域Rp,通常去除函数特征的数据点数值大,数据点少,集中在直方图曲线的中心位置,因此计算直方图中心最高点纵坐标数值下降90%对应的横坐标数值,得到去除函数特征数据点的分布范围[Vmin,Vmax],其中Vmin表示最小值,Vmax表示最大值,[]表示闭区间;
最后,将系数矩阵Sj中区域Rp中数据点数值不属于[Vmin,Vmax]的对应数据点置零;同时对系数矩阵Sj中不属于区域Rp的所有数据点置零,即去除寄生条纹的特征;
7、对处理之后的系数矩阵Sj(其中j=1~K)进行重构,重构时小波基函数为sym5,分解层数为N,最终得到了寄生条纹的面形数据矩阵f2(x,y)。
与现有技术相比,本发明的优点在于:
采用了sym5小波对含有寄生条纹的面形数据矩阵进行多级小波分解,实现了寄生条纹的充分分解,该方法对磁流变去除函数的频谱和寄生条纹的频谱严重混叠时仍然适用,而传统的频域滤波法由于无法分辨严重混叠的频谱失效,因此具有更加宽广的适用范围;
采用了去除函数特征区域匹配和直方图统计的方法,实现了系数矩阵中寄生条纹和去除函数特征的有效分离,保证了消除寄生条纹的同时不改变去除函数自身形态,提高了数据处理精度。
下面通过一个完整的实施例来说明基于小波变换的磁流变去除函数寄生条纹消除方法。
实施例1:
比如磁流变去除函数采用ZYGO静态干涉仪测量,采斑基片为Φ50mm,检测数据通光口径为Φ48mm,材料为熔石英,在干涉检测过程中引入了较为明显的寄生条纹。
通过下述步骤对寄生条纹进行滤除:
采用上述试验平台对寄生条纹进行滤除,流程如图3所示。
1、获取含有寄生条纹的干涉仪采斑测量面形数据矩阵,数据宽度像素数m=916,数据高度像素数n=916,像素尺寸为pix=0.052mm,如图4所示;
2、提取面形数据中不含有去除函数形状但是含有寄生条纹的局部数据,计算寄生条纹的空间周期T=0.48mm;
3、计算多级小波分解层数N,计算公式如下:
故N=4;
4、根据去除函数结构特点,选用sym5作为小波基函数。采用sym5小波对面形数据矩阵f1(x,y)进行多级分解,分解级数为N=4,每一级有三个系数矩阵,命名为hi(水平方向),vi(垂直方向),di(对角线方向),其中i=1~4,最后一级还有低频系数矩阵命名为A4,共得到13个系数矩阵数据,如图5所示;
5、分析系数矩阵hi,vi,di(其中i=1~4)的特征,记录含有寄生条纹的系数矩阵,从图5可以发现,系数矩阵v3和v4含有寄生条纹特征,因此Sj={v3,v4},j=1~2,K=2,K为含有寄生条纹的系数矩阵的个数;
6、对含有寄生条纹的系数矩阵Sj(其中j=1~2)进行处理,去除其中的寄生条纹特征,根据寄生条纹特征和去除函数特征数值不连续、存在强度差异的特点,采用基于直方图统计的方法确定寄生条纹特征所在取值范围,然后将系数矩阵中数值在该范围的数据置零,以S2=v4为例阐述实现过程i~iii:
i分析系数矩阵S2,观察图5,v4系数矩阵中含有4个去除函数,故H=4,搜索去除函数的特征分布所在区域,则Rp(p=1~4),如图6(a)所示,采用直方图统计的方法绘制曲线,横坐标为数据点的数值,纵坐标为数据点的个数,如图6(b)所示;
ii以R1区域为例阐述计算过程,计算直方图中心最高点纵坐标数值下降90%对应的横坐标数值,得到去除函数特征数据点的分布范围[Vmin,Vmax],如图6(b)所示,Vmin=-0.08,Vmax=0.07;
iii对系数矩阵S2中R1区域数据点数值不属于[Vmin,Vmax]的对应数据点置零,如图6(c)所示;同样方法对R2、R3、R4区域进行处理,最后对系数矩阵S2中R1、R2、R3、R4区域之外的全部数据点置零,如图6(d)所示;
7、在对系数矩阵S1、S2进行处理后,进行系数矩阵的小波重构,重构时小波基函数为sym5,分解层数为4,最终得到了寄生条纹的面形数据矩阵f2(x,y),如图7所示。上述完成了磁流变去除函数寄生条纹的消除。
基于同一发明构思,本发明实施例还提供了基于小波变换的磁流变去除函数寄生条纹消除装置,由于该装置所解决问题的原理与前述基于小波变换的磁流变去除函数寄生条纹的消除方法相似,因此该装置的实施可以参见前述方法的实施,重复之处不再赘述。
本发明实施例提供的基于小波变换的磁流变去除函数寄生条纹消除装置,参照图8所示,包括:
获取模块81,用于获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
计算模块82,用于提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
分解模块83,用于采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
记录模块84,用于分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
消除模块85,用于对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
在一个实施例中,所述基于小波变换的磁流变去除函数寄生条纹消除装置,还包括:
重构模块86,用于对消除模块85处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
在一个实施例中,所述计算模块82中,具体用于计算多级小波分解层数N,计算公式如下:
其中,N取满足条件下最小整数,T表示寄生条纹的空间周期,pix表示像素尺寸。
在一个实施例中,所述消除模块85,包括:
确定子模块851,用于根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
置零子模块852,用于将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (8)
1.基于小波变换的磁流变去除函数寄生条纹消除方法,其特征在于,包括:
获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
2.如权利要求1所述的消除方法,其特征在于,还包括:
对处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
3.如权利要求1所述的消除方法,其特征在于,所述N根据所述空间周期T计算获得,包括:
计算多级小波分解层数N,计算公式如下:
其中,N取满足条件下最小整数,T表示寄生条纹的空间周期,pix表示像素尺寸。
4.如权利要求1所述的消除方法,其特征在于,对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征,包括:
根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
5.基于小波变换的磁流变去除函数寄生条纹消除装置,其特征在于,所述装置包括:
获取模块,用于获取含有寄生条纹的干涉仪采斑测量面形数据矩阵f1(x,y);
计算模块,用于提取所述面形数据矩阵f1(x,y)中的局部数据,所述局部数据为:含有寄生条纹,且不含有去除函数形状;计算所述寄生条纹的空间周期T;
分解模块,用于采用sym5小波对所述面形数据矩阵f1(x,y)进行多级分解,分解级数为N,得到3N+1个系数矩阵数据;所述N根据所述空间周期T计算获得;
记录模块,用于分析3N+1个所述系数矩阵数据,记录含有寄生条纹的系数矩阵为Sj;其中j=1~K,K为含有寄生条纹的系数矩阵的个数,且K≤3N+1;
消除模块,用于对含有寄生条纹的系数矩阵Sj进行处理,消除其中的寄生条纹特征。
6.如权利要求5所述的消除装置,其特征在于,所述装置还包括:
重构模块,用于对消除模块处理后的系数矩阵Sj进行重构,重构时小波基函数为sym5,分解层数为N,得到消除寄生条纹的面形数据矩阵f2(x,y)。
7.如权利要求5所述的消除装置,其特征在于,所述计算模块中,具体用于计算多级小波分解层数N,计算公式如下:
其中,N取满足条件下最小整数,T表示寄生条纹的空间周期,pix表示像素尺寸。
8.如权利要求5所述的消除装置,其特征在于,所述消除模块,包括:
确定子模块,用于根据寄生条纹特征和去除函数特征数值不连续及存在的强度差异,采用基于直方图统计,确定寄生条纹特征所在取值范围;
消除子模块,用于将所述系数矩阵Sj中数值在所述取值范围内的数据置零,消除其中的寄生条纹特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410490.5A CN110227968B (zh) | 2019-05-16 | 2019-05-16 | 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410490.5A CN110227968B (zh) | 2019-05-16 | 2019-05-16 | 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110227968A true CN110227968A (zh) | 2019-09-13 |
CN110227968B CN110227968B (zh) | 2020-08-21 |
Family
ID=67860705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910410490.5A Active CN110227968B (zh) | 2019-05-16 | 2019-05-16 | 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110227968B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110929599A (zh) * | 2019-11-08 | 2020-03-27 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
CN111122507A (zh) * | 2019-11-06 | 2020-05-08 | 中国工程物理研究院激光聚变研究中心 | 一种基于离散小波变换的光学平板波前的降噪方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040008904A1 (en) * | 2003-07-10 | 2004-01-15 | Samsung Electronics Co., Ltd. | Method and apparatus for noise reduction using discrete wavelet transform |
US6985632B2 (en) * | 2000-04-17 | 2006-01-10 | Canon Kabushiki Kaisha | Image processing system, image processing apparatus, and image processing method |
CN103447891A (zh) * | 2013-08-26 | 2013-12-18 | 中国科学院光电技术研究所 | 一种磁流变高精度定位装置及磁流变去除函数转换方法 |
CN105243650A (zh) * | 2015-11-16 | 2016-01-13 | 成都神州数码索贝科技有限公司 | 一种基于小波变换和傅里叶变换的条纹噪声去除方法及*** |
CN105761231A (zh) * | 2016-03-21 | 2016-07-13 | 昆明理工大学 | 一种用于去除高分辨率天文图像中条纹噪声的方法 |
CN106556940A (zh) * | 2016-11-10 | 2017-04-05 | 武汉精测电子技术股份有限公司 | 一种tft‑lcd屏自动光学检测中的背景抑制方法 |
CN106826401A (zh) * | 2016-07-25 | 2017-06-13 | 中国科学院长春光学精密机械与物理研究所 | 一种磁流变抛光面形误差收敛控制加工方法 |
CN108596850A (zh) * | 2018-04-23 | 2018-09-28 | 中国科学院光电研究院 | 一种成像光谱仪条带噪声去除方法 |
-
2019
- 2019-05-16 CN CN201910410490.5A patent/CN110227968B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6985632B2 (en) * | 2000-04-17 | 2006-01-10 | Canon Kabushiki Kaisha | Image processing system, image processing apparatus, and image processing method |
US20040008904A1 (en) * | 2003-07-10 | 2004-01-15 | Samsung Electronics Co., Ltd. | Method and apparatus for noise reduction using discrete wavelet transform |
CN103447891A (zh) * | 2013-08-26 | 2013-12-18 | 中国科学院光电技术研究所 | 一种磁流变高精度定位装置及磁流变去除函数转换方法 |
CN105243650A (zh) * | 2015-11-16 | 2016-01-13 | 成都神州数码索贝科技有限公司 | 一种基于小波变换和傅里叶变换的条纹噪声去除方法及*** |
CN105761231A (zh) * | 2016-03-21 | 2016-07-13 | 昆明理工大学 | 一种用于去除高分辨率天文图像中条纹噪声的方法 |
CN106826401A (zh) * | 2016-07-25 | 2017-06-13 | 中国科学院长春光学精密机械与物理研究所 | 一种磁流变抛光面形误差收敛控制加工方法 |
CN106556940A (zh) * | 2016-11-10 | 2017-04-05 | 武汉精测电子技术股份有限公司 | 一种tft‑lcd屏自动光学检测中的背景抑制方法 |
CN108596850A (zh) * | 2018-04-23 | 2018-09-28 | 中国科学院光电研究院 | 一种成像光谱仪条带噪声去除方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111122507A (zh) * | 2019-11-06 | 2020-05-08 | 中国工程物理研究院激光聚变研究中心 | 一种基于离散小波变换的光学平板波前的降噪方法和装置 |
CN110929599A (zh) * | 2019-11-08 | 2020-03-27 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
CN110929599B (zh) * | 2019-11-08 | 2023-05-26 | 宁波大学 | 基于小波分析的岩体结构面轮廓曲线形态分离方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110227968B (zh) | 2020-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Adams et al. | Orientation imaging: the emergence of a new microscopy | |
CN105959684A (zh) | 基于双目融合的立体图像质量评价方法 | |
Herranz et al. | Scale-adaptive filters for the detection/separation of compact sources | |
CN110227968A (zh) | 基于小波变换的磁流变去除函数寄生条纹消除方法及装置 | |
CN105205788B (zh) | 一种针对高通量基因测序图像的去噪方法 | |
CN106556940A (zh) | 一种tft‑lcd屏自动光学检测中的背景抑制方法 | |
CN109242924A (zh) | 一种基于深度学习的核磁共振图像的降采样伪影去除方法 | |
CN108345033A (zh) | 一种微地震信号时频域初至检测方法 | |
CN103854267A (zh) | 一种基于变分与分数阶微分的图像融合与超分辨率实现方法 | |
CN109118453A (zh) | 一种背景抑制的图像处理方法 | |
CN102693522A (zh) | 一种彩色图像区域复制篡改检测方法 | |
CN115272303B (zh) | 基于高斯模糊的纺织物缺陷程度评估方法、装置及*** | |
EP1071042B1 (en) | Method for determining the components of image noise patterns of an imaging device and use of this method in an imaging device | |
CN104680488A (zh) | 一种基于中值滤波的dsm进行星载sar影像正射纠正的方法 | |
CN107388963A (zh) | 将小波分析和低通滤波结合的数字莫尔条纹相位提取方法 | |
CN107179550A (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
CN109100813B (zh) | 一种基于协同滤波消除地面核磁共振数据中尖峰噪声的方法 | |
CN103903235B (zh) | 细胞图像超分辨率重构方法 | |
CN103777237A (zh) | 一种基于空变加权镶边波数域滤波的地表高程平滑方法 | |
CN103020907B (zh) | 基于二维集合经验模态分解的dspi条纹滤波*** | |
Links et al. | Improved positron emission tomography quantification by Fourier-based restoration filtering | |
CN104616310B (zh) | 一种图像质量的评估方法和装置 | |
Selvin | Transformations of maps to investigate clusters of disease | |
CN111091601B (zh) | 一种实时白天室外手机图像的pm2.5指数估计方法 | |
CN108333545A (zh) | 基于高通滤波的磁共振图像重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |