一种光学相干层析变形场降噪方法与流程

本发明涉及光学技术领域，尤其涉及一种光学相干层析变形场降噪方法。

背景技术：

光学相干层析(OCT)是一种利用低相干光干涉原理的层析测量技术，它具有测量速度快和无损检测等优点，可以实现材料内部微米级轮廓的精确测量。在结合了光谱仪结构后发展出的谱域光学相干层析技术(SOCT)，进一步提高了该方法的测量速度与信噪比。近年来，以SOCT技术为基础又发展出了一系列新的功能型OCT方法，例如测量材料内部双折射特性的偏振敏感SOCT，测量材料内部流速场的多普勒SOCT，以及测量材料内部变形场的相位对照SOCT等，并被广泛应用于生物医学和先进材料等多项领域的研究与检测当中。其中，相位对照SOCT是一种利用干涉光谱相位敏感特性的SOCT方法，可以实现材料内部变形场的纳米级层析测量，被先后应用于眼角膜，高聚物等材料在不同载荷下的微变形场动态检测当中。

然而，相位对照SOCT方法是基于对干涉信号进行测量的技术，不可避免地会受到散斑噪声的影响。同时，材料内部散射信号很弱，且随着深度的增加不断减弱。这些都导致了测量结果的信噪比低，往往需要人为的参与，才能辨识出噪声污染下被测样件的截面轮廓，以及内部变形场细节上的变化趋势。

亟待解决上述问题，从而实现对测量结果的降噪，使截面轮廓更加清晰，变形场分布更加细致。

技术实现要素：

本发明提供一种光学相干层析变形场降噪方法，在相位对照谱域光学相干层析技术的基础上结合了模糊逻辑的方法，通过对测量结果提取特征参数，并进行模糊逻辑推理和判断，对样件的截面轮廓以及样件内部的变形场进行降噪，使截面轮廓更加清晰，变形场分布更加细致。

一种光学相干层析变形场降噪方法包括：

计算干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)；

根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄，其中C₁为附近区域的幅值之和、C₂为单一方向的幅值之和、C₃为相邻区域内相位差的连续程度、C₄为相邻区域内相位差拟合的误差；

对特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4；

对模糊化特征参数进行判断，输出O₁和O₂，其中O₁为该点是样件内部一点的可能性，O₂为该点是噪声点的可能性；

计算降噪后的输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)；

根据输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)解调出降噪后的变形场。

优选的，

步骤计算干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)包括：

将成像在CCD相机上的干涉光谱用光强表示为

其中，k表示波数，I₀表示直流分量，I_R表示参考面，I_i表示被测样件的反射光强，M表示被测样件的表面个数，φ_i表示干涉光强的相位，φ_i＝φ_i0+2kz，φ_i0表示初始相位；

将光强信号沿波数k轴进行傅里叶变换得到原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)；

通过公式计算被测样件变形前后的位移，其中深度方向的坐标z＝π·f_k。

优选的，

步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₁的取值范围；

通过公式计算C₁。

优选的，

步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

以(y,z)为中心分别沿y方向、z方向和yz轴之间45度夹角方向为特征参数C₂取值；

通过公式

计算C₂。

优选的，

步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₃的取值范围；

对区域内的相位差进行三阶曲面拟合得到Δφ_f(y,z)；

通过公式C₃(y,z)＝|Δφ_m(y,z)-Δφ_f(y,z)|计算C₃。

优选的，

步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₄的取值范围；

对区域内的相位进行三阶曲面拟合；

通过公式计算C₄。

优选的，

步骤对特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4包括：

通过公式分别计算得到模糊化特征参数C_n1、C_n2、C_n3和C_n4，其中，C_max表示相应特征参数的最大值，C_min表示相应特征参数的最小值；

选取3个模糊状态，分别标记为B、M和S；

绘制模糊集隶属度函数。

优选的，

步骤对模糊化特征参数进行判断，并输出O₁和O₂包括：

参照模糊集隶属度函数进行判断；

若满足条件1，则输出O₁是B；

否则，输出O₁是S；

条件1为参数C₁是B、参数C₂是B、参数C₄是S、参数C₁与C₄均是M和参数C₂与C₄均是M中的任意一种。

优选的，

步骤对模糊化特征参数进行判断，并输出O₁和O₂还包括：

若满足条件2，则输出O₂是B；

否则，输出O₂是S；

条件2为参数C₁是S、参数C₂是S、参数C₃是B、参数C₄是B和参数C₄是S且参数C₃是M中的任意一种。

优选的，

步骤所述计算降噪后的输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)包括：

通过公式

计算降噪后的幅值A_e(y,z)和相位差Δφ_e(y,z)。

从以上技术方案可以看出，本发明实施例具有以下优点：

本发明的光学相干层析变形场降噪方法，首先，计算出干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)。其次，根据原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算提取特征参数C₁、C₂、C₃和C₄，其中，C₁为附近区域的幅值之和，即选取一个区域对其内部的幅值求和；C₂为单一方向的幅值之和，即选取几个方向分别求和，取其中的最大值后再减去中心点自身的幅值；C₃为相邻区域内相位差的连续程度，即以(y,z)为中心取一个区域，对区域内的相位差进行三阶曲面拟合，再取拟合后中心点的相位差与原始相位差之间差的绝对值；C₄为相邻区域内相位差拟合的误差，即以(y,z)为中心选取一个区域，先对区域内的相位进行三阶曲面拟合，再用拟合后的相位差分别与原始相位差作差并求和。再次，对求得的特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4，并对模糊化特征参数进行判断，输出O₁和O₂，其中O₁为该点是样件内部一点的可能性，O₂为该点是噪声点的可能性。最后计算降噪后的幅值A_e(y,z)和相位差Δφ_e(y,z)，根据输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)解调出降噪后的变形场，实现对样件的截面轮廓以及样件内部的变形场进行降噪，使截面轮廓更加清晰，变形场分布更加细致。

附图说明

图1为本发明中一种光学相干层析变形场降噪方法实施例的流程图；

图2为本发明中一种光学相干层析变形场降噪方法实施例中模糊集隶属度函数图；

图3为本发明中一种光学相干层析变形场降噪方法实施例中含有散斑噪声截面轮廓y-z平面的幅值分布图；

图4为图3中y＝1mm处幅值随深度变化的曲线图；

图5为本发明中一种光学相干层析变形场降噪方法实施例中经过本方法降噪后截面轮廓y-z平面的幅值分布图；

图6为图5中y＝1mm处幅值随深度变化的曲线图；

图7为降噪处理之前样件内部的相位差分布图；

图8为使用本方法处理后样件内部的相位差分布图。

具体实施方式

本发明提供一种光学相干层析变形场降噪方法，在相位对照谱域光学相干层析技术的基础上结合了模糊逻辑的方法，通过对测量结果提取特征参数，并进行模糊逻辑推理和判断，对样件的截面轮廓以及样件内部的变形场进行降噪，使截面轮廓更加清晰，变形场分布更加细致。

本发明的一种光学相干层析变形场降噪方法实施例包括：

101、计算干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)；

102、根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄，其中C₁为附近区域的幅值之和、C₂为单一方向的幅值之和、C₃为相邻区域内相位差的连续程度、C₄为相邻区域内相位差拟合的误差；

103、对特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4；

104、对模糊化特征参数进行判断，输出O₁和O₂，其中O₁为该点是样件内部一点的可能性，O₂为该点是噪声点的可能性；

105、计算降噪后的输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)；

106、根据输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)解调出降噪后的变形场。

本实施例中，首先，计算出干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)。其次，根据原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算提取特征参数C₁、C₂、C₃和C₄，其中，C₁为附近区域的幅值之和，即选取一个区域对其内部的幅值求和；C₂为单一方向的幅值之和，即选取几个方向分别求和，取其中的最大值后再减去中心点自身的幅值；C₃为相邻区域内相位差的连续程度，即以(y,z)为中心取一个区域，对区域内的相位差进行三阶曲面拟合，再取拟合后中心点的相位差与原始相位差之间差的绝对值；C₄为相邻区域内相位差拟合的误差，即以(y,z)为中心选取一个区域，先对区域内的相位进行三阶曲面拟合，再用拟合后的相位差分别与原始相位差作差并求和。再次，对求得的特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4，并对模糊化特征参数进行判断，输出O₁和O₂，其中O₁为该点是样件内部一点的可能性，O₂为该点是噪声点的可能性。最后计算降噪后的幅值A_e(y,z)和相位差Δφ_e(y,z)，根据输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)解调出降噪后的变形场，参照如公式实现对样件的截面轮廓以及样件内部的变形场进行降噪，使截面轮廓更加清晰，变形场分布更加细致。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤101、计算干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)包括：

将成像在CCD相机上的干涉光谱用光强表示为

其中，k表示波数，I₀表示直流分量，I_R表示参考面，I_i表示被测样件的反射光强，M表示被测样件的表面个数，φ_i表示干涉光强的相位，φ_i＝φ_i0+2kz，φ_i0表示初始相位；

将光强信号沿波数k轴进行傅里叶变换得到原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)；

通过公式计算被测样件变形前后的位移，其中深度方向的坐标z＝π·f_k。

本实施例中，在不考虑材料自身相互干涉产生干涉信号的前提下，将干涉光谱用光强表示为再将光强信号沿波数k轴进行傅里叶变换得到原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)，其中，傅里叶变换后得到的幅值A(y,z)能够描述y-z平面上被测样件内部的截面轮廓结构，而被测样件变形前后的位移可以用相位差Δφ表示，从而为后续计算提取特征参数做准备。变形场与相位差之间是线性关系，通过公式计算，因此后面对相位差的去噪即是对变形场的去噪。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₁的取值范围；

通过公式计算C₁。

本实施例中，选取适当的特征区域，如以(y,z)为中心取5×5的一个区域，取n＝2，则有即对其内部的幅值求和，计算得出特征参数C₁的取值。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

以(y,z)为中心分别沿y方向、z方向和yz轴之间45度夹角方向为特征参数C₂取值；

通过公式

计算C₂。

本实施例中，以(y,z)为中心沿y方向，z方向，yz轴之间45度夹角的方向分别求和，取其中的最大值后再减去中心点自身的幅值，从而计算出特征参数C₂。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₃的取值范围；

对区域内的相位差进行三阶曲面拟合得到Δφ_f(y,z)；

通过公式C₃(y,z)＝|Δφ_m(y,z)-Δφ_f(y,z)|计算C₃。

本实施例中，以(y,z)为中心取7×7的一个区域为例，对区域内的相位差进行三阶曲面拟合，再取拟合后中心点的相位差与原始相位差之间的差值，再求绝对值，从而求得特征参数C₃。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤根据干涉光谱的原始幅值A_m(y,z)和原始相位差Δφ_m(y,z)计算特征参数C₁、C₂、C₃和C₄包括：

选取以(y,z)为中心的一个区域为特征参数C₄的取值范围；

对区域内的相位进行三阶曲面拟合；

通过公式计算C₄。

本实施例中，如以(y,z)为中心取7×7的一个区域，上式中取a＝3，则有即首先对区域内的相位进行三阶曲面拟合，再用拟合后的相位差分别与原始相位差作差并求和，从而计算特征参数C₄。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤对特征参数C₁、C₂、C₃和C₄进行模糊化处理得到模糊化特征参数C_n1、C_n2、C_n3和C_n4包括：

通过公式分别计算得到模糊化特征参数C_n1、C_n2、C_n3和C_n4，其中，C_max表示相应特征参数的最大值，C_min表示相应特征参数的最小值；

选取3个模糊状态，分别标记为B、M和S；

绘制模糊集隶属度函数。

本实施例中，参考图2，横坐标为特征参数，纵坐标为隶属度。隶属度函数是模糊控制的应用基础，是否正确地构造隶属度函数是能否用好模糊控制的关键之一。隶属度函数的确定过程,本质上说应该是客观的，但每个人对于同一个模糊概念的认识理解又有差异，因此，隶属度函数的确定又带有主观性。用3个模糊状态B、M和S，绘制模糊集隶属度函数，B表示大、M表示中和S表示小。本发明以如下函数为例：

绘制隶属度函数。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤对模糊化特征参数进行判断，并输出O₁和O₂包括：

参照模糊集隶属度函数进行判断；

若满足条件1，则输出O₁是B；

否则，输出O₁是S；

条件1为参数C₁是B、参数C₂是B、参数C₄是S、参数C₁与C₄均是M和参数C₂与C₄均是M中的任意一种。

本实施例中，参照模糊集隶属度函数进行判断，关于输出O₁的规则包括：如果参数C₁是B，则输出O₁是B；如果参数C₂是B，则输出O₁是B；如果参数C₄是S，则输出O₁是B；如果参数C₁是M，且参数C₄是M，则输出O₁是B；如果参数C₂是M，且参数C₄是M，则输出O₁是B；若不满足上述任一规则，则输出O₁是S。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤对模糊化特征参数进行判断，并输出O₁和O₂包括：

若满足条件2，则输出O₂是B；

否则，输出O₂是S；

条件2为参数C₁是S、参数C₂是S、参数C₃是B、参数C₄是和参数C₄是S且参数C₃是M中的任意一种。

本实施例中，参照模糊集隶属度函数进行判断，关于输出O₂的规则包括：如果参数C₁是S，则输出O₂是B；如果参数C₂是S，则输出O₂是B；如果参数C₃是B，则输出O₂是B；如果参数C₄是B，则输出O₂是B；如果参数C₄是S，且参数C₃是M，则输出O₂是B；若不满足上述任一规则，则输出O₂是S。

下面对本发明的一种光学相干层析变形场降噪方法实施例做进一步说明，步骤所述计算降噪后的输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)包括：

通过公式

计算降噪后的幅值A_e(y,z)和相位差Δφ_e(y,z)。

本实施例中，根据前述对O₁和O₂的范围的判断，计算降噪后的输出幅值A_e(y,z)和输出相位差Δφ_e(y,z)，从而对样件的截面轮廓以及样件内部的变形场进行降噪，使截面轮廓更加清晰，变形场分布更加细致。降噪前后的截面轮廓幅值可参阅图3-图6，相位差分布情况可参阅图7-图8，改图显示的是测量材料内部的一个截面，为三维空间中的y-z平面。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。在本申请所提供的几个实施例中，应该理解到，所揭露的装置可以通过其它的方式实现。

以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。