本发明涉及相移干涉测量领域，特别涉及一种基于概率分布函数的相移误差检测方法。
背景技术：
相移干涉测量法是一种高精度的光学相位测量技术，是光学表面面型检测的重要手段之一。现在通用的实现相移的方法主要有机械式和波长调谐式两种，机械式一般采用推动压电陶瓷(pzt)、旋转偏振器件(波片、偏振片)、移动衍射光栅或者倾斜平板等手段来实现相位调制。其中推动压电陶瓷(pzt)产生相移应用最为广泛。波长调谐式则是通过光源频率的改变来实现相位调制，其光源采用波长可调谐激光器。这两种相移方式都存在一定的相移误差，包括初始相移误差和由于硬件老化等原因导致的相移误差等，相移误差被认为是影响测量精度的主要误差之一。因此，要实现高精度的三维面形测量，就必须进行精确相移量的标定。
技术实现要素：
本发明的目的在于克服现有技术中相移干涉测量法测量相移存在相移误差的问题，提供一种基于概率分布函数的相移误差检测方法。
为了实现上述发明目的，本发明提供了以下技术方案：
一种基于概率分布函数的相移误差检测方法，其特征在于，n步相移干涉测量产生n-1个相移误差，n大于等于3，所述相移误差的检测方法包括：
步骤一、获取实际的截断相位分布，并计算出对应的pdf曲线；
步骤二、将所述pdf曲线与事先建立的模拟曲线数据库进行相关运算，找到所述模拟曲线数据库中与所述pdf曲线相关度最大的模拟曲线，该模拟曲线的相移误差作为检测出的相移误差。
优选的，所述模拟曲线数据库的建立方法：
在设定区间内模拟产生一系列不同相移误差情况下的pdf曲线，并建立模拟曲线数据库。
优选的，所述步骤一包括：
获取n幅相移图片,描述所述相移图片中像素的光强，再根据n步相移干涉测量理想情况下的相移量，得到带有相移误差的截断相位；
根据带有相移误差的截断相位得到对应的所述pdf曲线。
优选的，第n幅所述相移图片中像素的光强in(x,y)：
其中，a(x,y)，b(x,y)分别表示背景光强和调制度，
优选的，将检测n步相移干涉测量的相移误差的过程分解为检测
优选的，所述三步相移干涉测量的相移误差检测方法包括：
s1获取三幅相移图片,描述所述相移图片中像素的光强，再根据三步相移干涉测量理想情况下的相移量，得到有误差的截断相位；
根据有误差的截断相位得到对应的三步相移干涉测量pdf曲线；
s2将所述三步相移干涉测量pdf曲线与事先建立的三步相移干涉测量模拟曲线数据库进行相关运算，找到所述三步相移干涉测量模拟曲线数据库中与所述三步相移干涉测量pdf曲线相关度最大的模拟曲线，该模拟曲线的相移误差作为检测出的相移误差。
优选的，所述三步相移干涉测量模拟曲线数据库的建立：在设定区间内模拟产生一系列不同δδ1和δδ2情况下的pdf曲线f0(m,n)，m、n为δδ1、δδ2的编号。
优选的，所述步骤s1，第n幅相移图片任一像素的光强in(x,y)：
当三步相移为满周期三步相移时，相移量0、δ1＝2π/3、δ2＝4π/3，得到带有相移误差的截断相位
所述pdf曲线f：
其中，m代表采样点数，m＝0,1,2,...,m-1。
优选的，所述步骤s2中，所述相关运算采用以下公式来计算相关系数r(m,n):
其中，f为所述pdf曲线，f0(m,n)为模拟曲线；其中，m,n为相移误差的编号。
优选的，所述相关度最大的模拟曲线，为使相关系数r(m,n)达到最大值的模拟曲线。
与现有技术相比，本发明的有益效果：
本发明的一种基于概率分布函数的相移误差检测方法，获取实际的截断相位分布，并计算出对应的pdf曲线；将所述pdf曲线与事先建立的模拟曲线数据库进行相关运算，找到所述模拟曲线数据库中与所述pdf曲线相关度最大的模拟曲线，所述相关度最大的模拟曲线的相移误差作为检测出的相移误差。利用本方法检测出的相移误差与实际相移误差非常接近，通过缩短建立模拟曲线数据库时的采样间隔，还能进一步减少检测误差。本发明还通过将n步相移测量分解为三步相移测量，再进行相移误差检测，只需要建立有两个相移误差的模拟曲线数据库，降低了建立模拟曲线数据库的难度，减小了模拟曲线数据库的体积。
附图说明：
图1为本实施例截断相位及相位误差。
图2为理想情况及四种相移误差情况下的截断相位(a)和pdf曲线(b)。
图3为本实施例的计算过程流程图示意图。
图4为基于概率分布检测相移误差过程示意图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。
在相移技术中，需要获取n幅相移图片用于计算相位值，其中，第n幅条纹图任一像素的光强可以表示为：
其中，a(x,y)，b(x,y)分别表示背景光强和调制度，
n步相移干涉测量具有n-1个相移误差，当n大于3时，将检测n步相移干涉测量的相移误差的过程分解为检测
本实施例讨论n＝3的情况，所提方法同样适用于n>3的情况。在满周期的三步相移干涉测量中，理想情况下的相移量为：0、δ1＝2π/3、δ2＝4π/3，但在实际测量中往往存在两个相移误差量：δδ1和δδ2。相移误差的存在引入了周期性的相位误差，大大影响了三维面形测量的精度。因而，相移量的准确确定是相移干涉测量的核心，如何快速准确的对相移量进行标定，对保证测量的精度具有重要意义。运用统计的方法标定精确相移量，具有很好的噪声抑制能力，可以实现快速实时的相移量标定，消除周期性相位误差，提高测量精度。公式(1)可以化简为：
in(x,y)＝c0(x,y) c1(x,y)cos(δn) c2(x,y)sin(δn)，
(n＝0,1,2),(2)
公式(2)中，c0(x,y)＝a(x,y)，
方便起见，后文中我们将省略(x,y)，由公式(2)计算i1(x,y)-i0(x,y)和i2(x,y)-i0(x,y)可以得到：
解得c1，c2：
以及
若已知精确相移值，将公式(5)和(6)代入公式(3)则可计算出精确相位值。实际上，在基于相移技术的干涉计量中，机械结构的移动误差是相移误差的来源之一，此时实际相移量可认为是预设相移量加上误差项。即δn＇＝δn δδn
本实施例以满周期三步相移为例，但本算法同样适用于任意相移量的三步相移以及3步以上相移，即相移量为，0，δ1＝2π/3，δ2＝4π/3，此时截断相位计算公式可进一步简化为：
不失一般性，假定两次移动引入了两个相移误差，即δδ1和δδ2，图1中展示了几种相移误差情况下的截断相位和相位误差曲线，图1(a)为理想相位，图1(b)-(e)分别为当δδ1＝π/19、δδ2＝π/17，δδ1＝π/19、δδ2＝-π/17，δδ1＝-π/19、δδ2＝π/17，δδ1＝-π/19、δδ2＝-π/17时获取的截断相位，图1(f)-(i)为相应的相移误差引入的相位误差：(a)理想情况下的截断相位。(b)带有相移误差δδ1＝π/19和δδ2＝π/17的截断相位。(c)带有相移误差δδ1＝π/19和δδ2＝-π/17的截断相位。(d)带有相移误差δδ1＝-π/19和δδ2＝π/17的截断相位。(e)带有相移误差δδ1＝-π/19和δδ2＝-π/17的截断相位。(f)-(i)是(b)-(e)中各相移误差情况下相对应的相位误差。
可以看出，由相移误差导致的相位误差具有周期性。而且相位误差随着相移误差的变化而变化，不仅仅是幅值改变，而且还有整体移动。因此，只能是首先得到精确的相移量才能要消除相位误差。本发明采用统计的方法来搜索精确相移量。
1.截断相位的概率分布函数
概率分布函数(pdf)是一种常见的统计方法。我们将pdf作为截断相位的目标函数，可以将其表示为
其中，m代表采样点数，m＝0,1,2,...,m-1，一定范围内，m越大，结果越精确。本文中取m＝64。其含义为：将由三步相移条纹获取的截断相位
理论上，当测量系统中不存在相移误差的时候，每个相位值出现的概率是一样的，此时pdf曲线为均匀曲线。若存在相移误差，此时相位误差值将改变原精确相位值出现的概率，即，使得pdf曲线发生变化。图2(a)展示了几种相移误差情况下的截断相位，红色代表理想情况，即δδ1＝δδ2＝0，另外四种颜色分别代表四种不同δδ1和δδ2组合的情况。图2(b)为对应相移误差情况下的pdf曲线。显而易见，pdf曲线的幅值以及峰值位置会随着相移误差的变化而变化。当δδ1＝δδ2＝0时，pdf曲线是一条直线，每一个采样点的概率都是1/m。例如m＝64时f(m)＝0.0156。
2.基于概率分布函数的相移误差计算方法
由于pdf曲线与相移误差情况一一对应，基于这一特性，本发明提出了一种计算相移误差的方法。本文选择满周期三步相移，设定相移为0，2π/3，4π/3，但这一算法同样适用于任一相移量的三步相移以及3步以上相移，计算过程如图3所示。因此，本发明中实际相移误差的查找步骤如图4所示：
计算机在[-π/10，π/10]区间内模拟产生一系列不同δδ1和δδ2情况下的pdf曲线f0(m,n)，其中m,n为δδ1和δδ2的编号。图4(a)模拟的是δδ1和δδ2均落在区间[-π/10,π/10]时的pdf曲线，采样间隔为π/200，由于δδ1和δδ2的符号会影响pdf曲线出现峰值的位置，为方便理解，只展示了δδ1>0,|δδ2|<|δδ1|的部分。
获取实际的截断相位分布，并计算出对应的pdf曲线f；图4(b)展示的是待测条纹图的pdf曲线，这里模拟的是δδ1＝π/19，δδ2＝-π/17时的情况,即(0.165，-0.185)。
通过相关运算，找到模拟曲线簇f0中与实际曲线相关度最大的曲线，相关系数计算公式为：
当相关系数r取最大值时，模拟曲线对应的相移误差即为待测曲线的相移误差，也即实际测量系统存在的相移误差。待测曲线与模拟曲线的相关系数结果如图4(c)所示，相关系数最大值落在(0.157,-0.189)的位置，因此检测结果为δδ1＝0.157，δδ2＝-0.189，检测误差为(-0.008，-0.004)，几乎接近于正确值(注：检测误差与采样间隔有关，当模拟误差情况为曲线簇中的某一种情况时，检测误差为0)。
1.一种基于概率分布函数的相移误差检测方法，其特征在于，n步相移干涉测量产生n-1个相移误差，n大于等于3，所述相移误差的检测方法包括：
步骤一、获取实际的截断相位分布，并计算出对应的pdf曲线；
步骤二、将所述pdf曲线与事先建立的模拟曲线数据库进行相关运算，找到所述模拟曲线数据库中与所述pdf曲线相关度最大的模拟曲线，该模拟曲线的相移误差作为检测出的相移误差。
2.根据权利要求1所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述模拟曲线数据库的建立方法：
在设定区间内模拟产生一系列不同相移误差情况下的pdf曲线，并建立模拟曲线数据库。
3.根据权利要求1所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述步骤一包括：
获取n幅相移图片,描述所述相移图片中像素的光强，再根据n步相移干涉测量理想情况下的相移量，得到带有相移误差的截断相位；
根据带有相移误差的截断相位得到对应的所述pdf曲线。
4.根据权利要求1所述的一种基于概率分布函数的相移误差检测方法，其特征在于，第n幅所述相移图片中像素的光强in(x,y)：
其中，a(x,y)，b(x,y)分别表示背景光强和调制度，
5.根据权利要求1-4任一项所述的一种基于概率分布函数的相移误差检测方法，其特征在于，将检测n步相移干涉测量的相移误差的过程分解为检测
6.根据权利要求5所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述三步相移干涉测量的相移误差检测方法包括：
s1获取三幅相移图片,描述所述相移图片中像素的光强，再根据三步相移干涉测量理想情况下的相移量，得到有误差的截断相位；
根据有误差的截断相位得到对应的三步相移干涉测量pdf曲线；
s2将所述三步相移干涉测量pdf曲线与事先建立的三步相移干涉测量模拟曲线数据库进行相关运算，找到所述三步相移干涉测量模拟曲线数据库中与所述三步相移干涉测量pdf曲线相关度最大的模拟曲线，该模拟曲线的相移误差作为检测出的相移误差。
7.根据权利要求6所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述三步相移干涉测量模拟曲线数据库的建立：在设定区间内模拟产生一系列不同δδ1和δδ2情况下的pdf曲线f0(m,n)，m、n为δδ1、δδ2的编号。
8.根据权利要求7所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述步骤s1，第n幅相移图片任一像素的光强in(x,y)：
当三步相移为满周期三步相移时，相移量0、δ1＝2π/3、δ2＝4π/3，得到带有相移误差的截断相位
所述pdf曲线f：
其中，m代表采样点数，m＝0,1,2,...,m-1。
9.根据权利要求7所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述步骤s2中，所述相关运算采用以下公式来计算相关系数r(m,n):
其中，f为所述pdf曲线，f0(m,n)为模拟曲线；其中，m,n为相移误差的编号。
10.根据权利要求7所述的一种基于概率分布函数的相移误差检测方法，其特征在于，所述相关度最大的模拟曲线，为使相关系数r(m,n)达到最大值的模拟曲线。
技术总结