一种基于频率域滤波的干涉图背景去除方法
【专利摘要】本发明公开了一种基于频率域滤波的干涉图背景去除方法，包括以下步骤：首先，对偏振干涉成像光谱仪获得的干涉图进行二维傅里叶变换并平移获得干涉图的频谱图；然后，将干涉图的频谱图与根据该偏振干涉成像光谱仪设计的陷波滤波器对应点相乘；最后，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹。与现有方法相比，本发明具有快速、准确及稳定的特点，尤其适用于复杂背景去除，并且对干涉图噪声具有一定的压制作用，从而提高了复原光谱的信噪比。
【专利说明】
一种基于频率域滤波的干涉图背景去除方法
技术领域
[0001]本发明涉及图像处理技术领域，特别涉及一种干涉图背景去除方法。
【背景技术】
[0002]图像和光谱是描述目标特性的重要信息，对遥感探测具有重要的意义。偏振干涉成像光谱仪能够同时获取目标的图像信息和光谱信息，在环境监测、大气探测、目标识别及生物医学领域有着重要的应用。在偏振干涉成像光谱仪获取的原始数据中，干涉条纹和背景叠加在一起，称之为干涉图。含有背景项的干涉图不能直接用于光谱复原，因为背景项的存在会导致复原光谱的信噪比降低，分辨率降低，谱线位置精度下降甚至出现错误。因此干涉图去背景是偏振干涉成像光谱仪光谱复原过程中的重要环节，干涉条纹切趾、相位修正等步骤均需要在背景去除之后进行。而且，干涉图背景去除的质量对复原光谱的质量有着直接影响。目前工程中常用的干涉图背景去除方法有一阶差分方法和最小二乘拟合方法。一阶差分方法具有原理简单、运算量小的优点，但是对干涉图噪声敏感，容易受到干涉图空间频率变化的影响，对复杂背景的去除效果不好;最小二乘拟合方法由于多项式阶数限制，不能够对复杂背景进行精确的拟合，并且在高阶情况下运算量很大，效率较低。
[0003]除此之外，2011年，穆廷魁等人提出了一种通过光学辅助去除干涉图背景的方法(T Mu,C Zhang ,ff Ren,X Jian,^Static dual-channel polarizat1n imagingspectrometer for simultaneous acquisit1n of inphase and antiphaseinterference images ,"Measurement Science&Technology22.10:105302-105309(8)(2011).)，可以避免复杂的计算机运算，具有集成度高、灵敏度高及稳定性好的优点。但是其对偏振片的位置精度要求很高，并且不易应用于其他类型的偏振干涉成像光谱仪，实用性不好。2013年，张淳民等人提出了一种基于经验模态分解(EMD)的干涉图背景去除方法(CZhang,ff Ren，T Mu，L Fu，C Jia，"Empirical mode decomposit1n based backgroundremoval and de-noising in polarizat1n interference imaging spectrometer，〃Optics Express 21.21:2592-2605(2013))。这种方法无需光谱信息的先验知识，具有良好的自适应性，能够较为彻底地去除背景。但是其算法复杂，运算量很大，耗时较长，限制了其在工程实际中的应用。
【发明内容】
[0004]本发明的目的在于提供一种基于频率域滤波的干涉图背景去除方法，以克服上述现有技术存在的缺点。本发明方法能够根据偏振干涉成像光谱仪的相关参数，实现其相应陷波滤波器的设计，使用该滤波器对干涉图进行频率域滤波，从而获得去除背景的干涉条纹，解决了常规干涉图背景去除方法的精确度低、运算量大及易受噪声影响的问题，适用于工程实际，能够快速准确地将干涉条纹从二维干涉图中直接提取出来。
[0005]为了实现上述目的，本发明采用如下技术方案:
[0006]—种基于频率域滤波的干涉图背景去除方法，包括以下步骤:首先，对偏振干涉成像光谱仪获得的干涉图进行二维傅里叶变换并平移获得干涉图的频谱图；然后，将干涉图的频谱图与根据该偏振干涉成像光谱仪设计的陷波滤波器对应点相乘;最后，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹。
[0007]进一步的，获取的干涉图的尺寸为M像素XN像素，干涉条纹的最大宽度为Wmax像素，最小宽度为Wmin像素;所述陷波滤波器的尺寸与复色光干涉图的尺寸相同；陷波滤波器对称轴左侧的通带起始点的位置为(M/2，N(Wmin-l)/2Wmin)，终止点的位置为(M/2，N(Wmax-l)/2Wmax)，以起始点到终止点连线上的每个点为圆心，作半径为R的圆，所有圆的闭合包络线内部的范围划定为通带，对称轴右侧的通带根据左侧通带按照轴对称方式设置;陷波滤波器上两个通带以外的区域为阻带。
[0008]进一步的，R值为 M/100 ?M/50。
[0009]进一步的，通带中心线值为I，通带边缘值为0，通带中心线到边缘之间的值随着其距中心线的距离D设置为0.5*(1-COS(2jt(R-D)/2R));阻带的值为O。
[0010]进一步的，所述干涉图为复色光干涉图或单色光干涉图。
[0011]相对于现有技术，本发明具有以下有益效果:本发明一种基于频率域滤波的干涉图背景去除方法，实现了快速准确地去除干涉图背景。首先，在运算速度方面远快于6阶以上最小二乘拟合方法以及经验模态分解方法，能够适用于工程实际;其次，在准确度方面不受噪声及复杂背景的影响，提取出的干涉条纹质量高于上述方法，背景项去除彻底，并对噪声有一定的压制作用，从而提高了复原光谱的信噪比；最后，陷波滤波器能够根据不同偏振干涉成像光谱仪的参数进行设计，具有良好的可移植性与实用性，滤波器在同一台偏振干涉成像光谱仪中是确定的，在使用过程中不需更改或修正，具有良好的稳定性。
【附图说明】
[0012]图1为本发明中陷波滤波器的结构示意图；
[0013]图2为本发明中陷波滤波器左侧通带结构示意图；
[0014]图3为本发明中干涉图复杂背景去除过程图；
[0015]图4为本发明干涉图背景去除效果图；
[0016]图5为最小二乘拟合法背景去除效果图；
[0017]图6为经验模态分解法背景去除效果图；
[0018]图7为使用本发明去背景后单色光谱复原效果图；
[0019]图8为使用最小二乘拟合法去背景后单色光谱复原效果图；
[0020]图9为使用经验模态分解法去背景后单色光谱复原效果图；
[0021]其中，I通带、2阻带、21为通带起始点、22为中心线、23为边缘、24为通带终止点、31为复色光干涉图、32为复色光干涉图的频谱图、33为陷波滤波器、34为去背景后的干涉条纹。
【具体实施方式】
[0022]本发明一种基于频率域滤波的干涉图背景去除方法的基本原理为:对于一台特定的偏振干涉成像光谱仪，其具有特定的波长范围，获得的干涉图是一个尺寸为M像素XN像素的二维图像;根据傅里叶变换光谱学原理，单色光干涉条纹的宽度和波长成正比，可以由仪器参数确定。对干涉图进行二维傅里叶变换，可以得到其频谱图，频谱图尺寸与干涉图一致，反映了干涉图的空间频率特征;使用相应的陷波滤波器，就可以从频谱图中选出对应于干涉图中干涉条纹的频谱，从而实现背景的去除。
[0023]为了达到上述目的，首先应依照以下方法设计相应的陷波滤波器:其尺寸与干涉图频谱图的尺寸一致，具有特定的通带I及阻带2;
[0024]通带I分为关于滤波器竖直中心轴线对称的两部分，设滤波器的尺寸为M像素XN像素，干涉条纹的最大宽度为Wmax像素，最小宽度为Wmin像素;滤波器对称轴(垂直于y方向，SPN像素一侧)左侧的通带I起始点21的位置为(M/2，N(Wmin-l)/2Wmin)，终止点24的位置为(M/2，N(Wmax-l)/2Wmax)，以起始点到终止点连线上的每个点为圆心，作半径为R的圆，R值可取为1/100?1/50之间的数值，所有圆的闭合包络线内部的范围划定为通带，对称轴右侧的通带根据左侧通带按照轴对称方式设置。
[0025]通带I中心线22值为I，通带边缘23值为0，通带中心线22到边缘23之间的值随着其距中心线的距离D设置为0.5*(1-cos(2jt(R-D)/2R));阻带2的值设置为O。
[0026]本发明一种基于频率域滤波的干涉图背景去除方法，包含以下步骤:首先，对偏振干涉成像光谱仪获得的复色光干涉图31进行二维傅里叶变换并平移获得复色光干涉图的频谱图32;然后，将复色光干涉图的频谱图32与根据该偏振干涉成像光谱仪设计的陷波滤波器33对应点相乘;最后，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹34。
[0027]下面以两个实例对本发明做进一步详细描述:
[0028]实施例1:
[0029]实验用偏振干涉成像光谱仪所获取的干涉图尺寸为512像素X512像素，探测波长范围为480nm-960nm，干涉条纹的最大宽度为2像素，最小宽度为I像素。则陷波滤波器尺寸为512像素X 512像素，对称轴左侧通带起始点的位置为(256，0)，终止点的位置为(256，128)。以起始点到终止点连线上的每个点为圆心，作半径为8像素的圆，所有圆的闭合包络线内部的范围划定为通带，陷波滤波器对称轴右侧的通带根据左侧通带按照轴对称方式设置。通带中心线值为I，通带边缘值为0，通带中心线到边缘之间的值随着其距中心线的距离D设置为0.5*(l-COS(23i(R-D)/2R));阻带的值设置为0，至此，陷波滤波器设置完成。
[0030]对实验用偏振干涉成像光谱仪所获取的复色光干涉图进行去背景处理。如图3所示:首先，对干涉图进行二维傅里叶变换并平移获得干涉图频谱图；然后，将干涉图频谱图与所述的陷波滤波器对应点相乘;最后，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹。在此实例中，将本发明与最小二乘拟合法去背景和经验模态分解法去背景比较，由相同配置MATLAB运行环境下的程序运行时间可以推算出，本发明的运算速度约为最小二乘拟合法的60倍，经验模态分解方法的10万倍；同时，由图4、图5及图6实验结果证明，在复杂背景去除的效果方面，本发明也是最优秀的。
[0031]实施例2:
[0032]实验用偏振干涉成像光谱仪所获取的干涉图尺寸为512像素X512像素，探测波长范围为480nm-960nm，干涉条纹的最大宽度为2像素，最小宽度为I像素。则陷波滤波器尺寸为512像素X 512像素，对称轴左侧通带起始点的位置为(256，0)，终止点的位置为(256，128)。以起始点到终止点连线上的每个点为圆心，作半径为8像素的圆，所有圆的闭合包络线内部的范围划定为通带，陷波滤波器对称轴右侧的通带根据左侧通带按照轴对称方式设置。通带中心线值为I，通带边缘值为O，通带中心线到边缘之间的值随着其距中心线的距离D设置为0.5*(l-COS(23i(R-D)/2R));阻带的值设置为O，至此，陷波滤波器设置完成。
[0033]对实验用偏振干涉成像光谱仪所获取的单色光干涉图进行去背景处理。首先，对干涉图进行二维傅里叶变换并平移获得干涉图频谱图；然后，将干涉图频谱图与所述的陷波滤波器对应点相乘;接下来，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹；最后，对获取的去除背景的干涉条纹进行复原获得光谱；由图7、图8及图9结果对比表明，本发明对复原光谱的噪声相比于其他方法(最小二乘拟合法、经验模态分解法)具有压制作用，从而提高了复原光谱信噪比。
【主权项】
1.一种基于频率域滤波的干涉图背景去除方法，其特征在于，包括以下步骤: 首先，对偏振干涉成像光谱仪获得的干涉图进行二维傅里叶变换并平移获得干涉图的频谱图； 然后，将干涉图的频谱图与根据该偏振干涉成像光谱仪设计的陷波滤波器对应点相乘； 最后，将相乘结果进行二维傅里叶逆变换获得去除背景的干涉条纹。2.根据权利要求1所述的一种基于频率域滤波的干涉图背景去除方法，其特征在于，获取的干涉图的尺寸为M像素X N像素，干涉条纹的最大宽度为Wmax像素，最小宽度为Wmin像素；所述陷波滤波器的尺寸与干涉图的尺寸相同；陷波滤波器对称轴左侧的通带起始点的位置为(1/2 4(1?-1)/2^，终止点的位置为(1/2，1^_-1)/21_)，以起始点到终止点连线上的每个点为圆心，作半径为R的圆，所有圆的闭合包络线内部的范围划定为通带，对称轴右侧的通带根据左侧通带按照轴对称方式设置；陷波滤波器上两个通带以外的区域为阻带。3.根据权利要求2所述的一种基于频率域滤波的干涉图背景去除方法，其特征在于，R值为 M/100?M/50。4.根据权利要求2所述的一种基于频率域滤波的干涉图背景去除方法，其特征在于，通带中心线值为I，通带边缘值为O，通带中心线到边缘之间的值随着其距中心线的距离D设置为0.5*(1-cos(2jt(R-D)/2R));阻带的值为O。5.根据权利要求1所述的一种基于频率域滤波的干涉图背景去除方法，其特征在于，所述干涉图为复色光干涉图或单色光干涉图。
【文档编号】G06T5/00GK106056542SQ201610338696
【公开日】2016年10月26日
【申请日】2016年5月19日
【发明人】张淳民, 颜廷昱, 李祺伟
【申请人】西安交通大学