CN106595522A - 一种光栅投影三维测量系统的误差校正方法 - Google Patents

Abstract

本发明公开了一种光栅投影三维测量系统的误差校正方法，其目的在于当测量表面反射率低的物体时，利用图像调制度中承载的物体反射率信息，对原始图像进行像素调整，从而提高相位质量。其实现步骤为：首先用相机采集受物体影响的变形条纹；然后用采集到的图像计算调制度系数I"并将其进行归一化；设定分割阈值T，将所有I"＜T的像素点分类，每一类中所有像素对应的I"值相近；对于每一类中所有像素，可以获取若干组N个灰度值，将这些组灰度值进行求均值，可获取一组灰度值曲线；最后用该组平均灰度值来替换该类中所有像素对应的N个灰度值，完成图像的校正工作。校正之后的图像可以用于计算相位，相位信息可以最终转化为待测物体的三维信息。

Description

一种光栅投影三维测量系统的误差校正方法

技术领域

本发明涉及一种光栅投影三维测量系统误差校正方法，属于计算机视觉中三维重构的领域。

背景技术

基于光栅投影的三维测量技术FPP(fringe projection profilometry)由于其精度高，速度快，受环境光影响较小等优点，近年来受到了广泛的研究和应用。作为一种基于主动光投影的三维测量方法，FPP也有相应的局限性。其中最为明显的一个是测量系统中所采用的工业相机在实际测量过程中会产生多种噪声，该噪声会降低采集到的图片的信噪比(signal-to-noise ratio)，影响光栅图像的质量，从而影响求解的相位质量以及最终的三维重构精度。该现象尤其在测量低反射率物体时较为严重。

对于表面纹理复杂，尤其是具有较暗纹理的物体，FPP系统采集到的图片信噪比较小，从而在暗部纹理图像部分求解得的相位质量较差。针对这一问题，目前绝大部分解决方法或是通过提高相机光圈大小，快门速度和相机增益，或是整体提高投影仪投影光强，从而使得拍摄到的图片中，暗部纹理部分获得充足曝光值。该类方法可以有效提高暗部细节的相位信息，但其缺点为增加光圈大小，曝光时间和相机增益会造成部分图像饱和，即图像灰度值达到255(对于图像格式为8位的工业相机)，尤其是当物体表面纹理非常暗的时候，该类方法会造成严重的图像饱和。同时，该类方法往往需要对一个物体进行若干次测量，然后将不同的测量结果融合成为一个最终的结果。整个测量过程操作较为繁琐，且相机参数的调节难以量化。

造成暗部表面纹理对光栅图像质量有影响的主要原因有两个：相机的采样效应和随机噪声。相机的采样效应对光栅图像质量的影响可以通过使用更高位数的相机或者采用上述多重曝光和调整投影仪亮度的方法进行改善；随机噪声则需要额外设计算法进行补偿，这也是本专利的主要内容。

发明内容

发明目的：光栅投影三维测量系统在测量具有较暗纹理的物体时，会出现较为明显的相位误差。针对此问题，本发明提供一种基于对所采集的光栅条纹图像进行调制度分析并直接对图像进行校正，从而校正最终的相位误差的方法。该方法无需借助除测量系统本身以外任的硬件，也无需对同一物体进行多次测量，仅通过分析光栅条纹图的调制度信息，实现低反射率部分的相位误差补偿。

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

一种光栅投影三维测量系统的误差校正方法，包括：

(1)用投影仪向待测物体表面投射N幅相移正弦光栅图像并用相机采集；

(2)根据采集得到的N幅相移条纹图求解调制度系数；

(3)基于调制度系数确定物体表面发射率较低的部分，对应的像素点为需要处理的像素点；

(4)将所有要处理的像素点根据调制度系数进行分类，同一类的像素点的调制度系数相近；

(5)对于每一类中若干个像素，每个像素对应一组N个灰度值，将这些像素的若干组N个灰度值对应位置进行求均值操作，然后用均值替换原像素的灰度值，实现图像校正；

(6)基于校正后的图像计算主值相位并最终求解物体的三维信息。

所述步骤(2)中根据采集到的N幅相移条纹图I_i,i＝1,2,..,N，求解调制度系数I″的公式为：

其中δ_i是每一步的相移量。

作为优选，所述步骤(2)中求解得到调制度系数后根据如下公式归一化至区间[0,1]：

I_z″＝(I″-I″_min)/(I″_max-I″_min)

其中I_z″为归一化之后的调制度系数，I″_max，I″_min分别表示求取的调制度系数I″的最大和最小值。

所述步骤(3)中通过设定一个分割阈值，将归一化后的调制度系数低于设定阈值的像素点确定为需要处理的像素点；作为优选，所述分割阈值取值范围为0.25～0.4。

所述步骤(4)中通过设定一个最小阈值，将每两个像素点对应的调制度系数差值小于设定阈值的像素点划分在同一类中；作为优选，所述最小阈值取值范围为0.008～0.012。

所述步骤(6)中主值相位的计算公式为：

其中，I_n为校正后的第n幅图像的灰度值。

有益效果：本发明提出一种光栅投影三维测量系统测量低反射率物体时的误差校正方法，相比现有技术，具有以下有益效果：

本发明针对传统光栅投影三维测量系统在测量具有较暗纹理的物体时容易产生相位误差的问题，提出了基于条纹调制度分析的相位误差校正算法。在光栅投影三维测量中，物体表面反射率相同的点在成像之后对应像素的调制度也相同，这些相同的像素点对应的图像噪声的方差也一样。所以在实际测量中，可以将调制度较低的像素点进行分类，每一类中的若干个像素点对应的调制度系数近似相等(两两之间的差值小于0.008～0.012)。将这些像素点看作是对物体上同一反射率的点进行的多次采样，则可以通过对这些像素点对应的若干组，每组N个灰度值进行类似均值滤波处理。处理后的条纹图即可用来获取校正后的相位。不同于采用多重曝光和改变投影仪亮度的方法，本发明算法实现过程简便，无需对同一物体进行多次测量，用数学方法即可有效减小随机噪声对相位的影响，明显改进暗部纹理区域对应的相位质量，从而提高暗部纹理对应的三维重构精度。

附图说明

图1是发明的整个过程的流程图。

图2是光栅投影三维测量系统框架图。

图3是具有丰富纹理的待测物体示意图。

图4是暗部纹理处某像素对应的N个灰度值曲线图。

图5是M类中某一类若干个像素对应的多组N个灰度值曲线示意图。

图6是图4中数据的均值曲线图。

图7是用本专利方法校正之前的相位误差结果图。

图8是用本专利方法校正之后的相位误差结果图。

图9是应用本专利校正方法之前相位误差的直方图。

图10是应用本专利校正方法之后相位误差的直方图。

具体实施方式

下面结合附图和具体实施例，进一步阐明本发明。在Windows操作系统下选用MATLAB作为编程工具，对计算机生成的正弦光栅以及CCD相机采集到的光栅图像进行处理。该实例采用具有黑色纹理的白色平面作为被测物体，证实本专利提出的误差校正方法的有效性。应理解这些实例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

本发明实施例公开的一种光栅投影三维测量系统测量低反射率物体时的误差校正方法，首先用投影仪向待测物体表面投射N幅相移正弦光栅图像并用相机采集；然后根据采集得到的N幅相移条纹图求解调制度系数，基于调制度系数确定物体表面发射率较低的部分，对应的像素点为需要处理的像素点；接着将所有要处理的像素点根据调制度系数进行分类，每一类中若干像素对应的若干组N个灰度值相应位置进行求均值，然后用均值替换原像素灰度值，实现图像校正；最后基于校正后的图像计算主值相位并最终求解物体的三维信息。算法主要流程如图1所示。测量系统结构框图如图2所示，I^s为计算机生成用于测量的标准相移条纹图；I^p为投影仪投射出的光栅图；I^o和I^c分别为经物体反射之后的条纹图和最终相机采集到的条纹图。在整个条纹的采集过程中，采集到的条纹图质量决定了相位质量，也影响最终的三维重构精度。对于具有较暗纹理的待测物体，计算得到的相位会收到明显影响。本专利通过对采集到的条纹进行调制度分析，应用统计学思想对条纹灰度值进行处理，有效减少了暗部纹理处的相位误差，提高系统的重构精度。

本发明实施例的具体实施过程包括以下步骤：

步骤1:根据主动光投影三维测量系统中的硬件三角关系固定投影仪和摄像机，将表面纹理复杂的待测物体放置在合适的位置。使用投影仪在物体上投射所需的N幅标准相移正弦光栅图像I，条纹灰度值设置为：

其中，I(i,j)为光栅图像I在第i行第j列处的灰度值，p为光栅条纹周期，为光栅的相移量。为了简便，在下文描述中省略掉像素坐标(i,j)。

步骤2:将摄像机相关参数：光圈大小，快门速度和感光度进行合理设置，使得采集回来的图像不会出现图像饱和(即最亮区域灰度值小于255)。在此相机参数下对N幅条纹图进行采集。图3为采集到的一幅相移光栅条纹图。采集到的条纹灰度值为：

I_n＝I′+I″cos[φ+2πn/N]

其中，n＝1,2,...,N，I_n为采集到的第n幅图像的灰度值，I′为条纹光强的背景值，I″为调制强度，φ为待求的主值相位分布。I′和I″与条纹图I_n拥有相同的分辨率。

步骤3:对于步骤2中采集得到的条纹图，求解条纹图的调制度系数I″，并将该系数进行归一化处理。I″反映所拍摄到的物体表面每个像素的反射率信息，与步骤2中采集到的图片中的像素是一一对应的关系。条纹调制度系数的计算方法：

步骤3.1:对于采集到的N幅相移条纹图I_i,i＝1,2,..,N，根据下式计算条纹调制度系数：

其中δ_i是每一步的相移量。

所述步骤3中条纹调制度系数的归一化方法如下：

步骤3.2:对于计算得到的调制度系数I″，根据下式将其归一化至区间0～1。

I_z″＝(I″-I″_min)/(I″_max-I″_min)

其中I_z″为归一化之后的调制度系数，I″_max，I″_min分别表示求取的调制度系数I″的最大和最小值。

步骤4:在区间0.25～0.4内设置一个阈值T,如T＝0.3，根据该阈值将步骤2中归一化之后的调制度系数I″进行阈值分割，对应调制度系数低于该阈值的像素点可以当作物体表面低反射率暗部纹理部分，也是本专利要处理的部分；高于该阈值的像素点则是对应物体表面反射率较高的区域，本专利中对此区域不做处理。阈值T的具体选择应根据不同的测量场景：对于不同的场景，分析其调制度系数I″，未归一化之前的调制度系数最大值I″_max和最小值I″_min差值越大，则阈值T应设置为越小，即越接近0.25；该差值越小，则T应越接近0.4。对于本专利所要处理的每一个像素，在N幅相移光栅条纹图中可以获取N个灰度值，则该N个灰度值的分布应呈正弦性。

图4是图3中A点对应的N个灰度值的曲线图，A点是物体上暗部纹理处反射率为0.1的像素点，容易看到由于相机的采样效应和随机噪声的影响，其正弦性很差。

步骤5:将步骤3中需要处理的所有像素进一步分为M类，分类规则是每一类中的像素点两两之间调制度系数I″相差不超过设定的阈值，阈值取值可在0.008～0.012范围内选取，如0.01。则每一类中的若干个像素点可以近似看作是从物体上同一反射率的纹理部分采集得到的。

图5是图3中与A点处于同一类的10个像素的灰度值曲线，可以看到该10个像素对应的灰度值曲线由于随机噪声和采样效应的影响，正弦性均较差。

步骤6:对于M类中每一类的若干像素，可以看作是对同一反射率的点的若干次采样。若某一类中有Q个像素，则表示第q个像素在第n步相移图中的灰度值，q＝1,2,...,Q，n＝1,2,...,N。对于每一步n相移，进行如下平均操作：

其中，为求均值操作后得到的平均灰度值，即校正之后的图像灰度值。该平均操作可以有效减小随机噪声对N个灰度值分布的正弦性影响。用平均之后的N个灰度值替换原始该类中的所有像素对应的N个灰度值，则完成了该类间的灰度值校正工作。图5中的10个像素对应的灰度值曲线图取平均之后如图6所示。对比图4与图6，能发现图6的正弦性要明显好于图4校正前的正弦性。对所有的M类重复这一校正操作，则完成了原始图像的校正工作。校正之后的图像可以用于获取更准确的相位。相位计算公式为：

用38步相移测量同一物体20次，将20次解得的相位的均值作为标准相位。将本专利的条纹校正方法用于4步相移，校正前后的条纹图求得的相位值与标准相位的差如图7和图8所示。可以发现经过本专利方法的校正，相位质量得到明显提升，平均相位误差降至校正前的1/5。

图9和图10为上述误差补偿前后的相位误差直方图，随机噪声引起的相位误差是零均值高斯分布。可以看出采用本专利误差校正方法之后，相位误差的方差明显减小。

步骤7：对主值相位进行展开得到绝对相位，根据经典光栅投影的相位到高度的转换公式，最终求得测量物体的三维信息。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.一种光栅投影三维测量系统的误差校正方法，其特征在于，包括：

(1)用投影仪向待测物体表面投射N幅相移正弦光栅图像并用相机采集；

(2)根据采集得到的N幅相移条纹图求解调制度系数；

(3)基于调制度系数确定物体表面发射率较低的部分，对应的像素点为需要处理的像素点；

(4)将所有要处理的像素点根据调制度系数进行分类，同一类的像素点的调制度系数相近；

(5)对于每一类中若干个像素，每个像素对应一组N个灰度值，将这些像素的若干组N个灰度值对应位置进行求均值操作，然后用均值替换原像素的灰度值，实现图像校正；

(6)基于校正后的图像计算主值相位并最终求解物体的三维信息。

2.根据权利要求1所述的一种光栅投影三维测量系统的误差校正方法，其特征在于，所述步骤(2)中根据采集到的N幅相移条纹图I_i,i＝1,2,..,N，求解调制度系数I″的公式为：

$I^{\′\′}=\frac{2}{N}\sqrt{({\left({\mathrm{\ΣI}}_i{\mathrm{cos\δ}}_i\right)}^2+{\left({\mathrm{\ΣI}}_i{\mathrm{sin\δ}}_i\right)}^2)}$

其中δ_i是每一步的相移量。

3.根据权利要求1所述的一种光栅投影三维测量系统的误差校正方法，其特征在于，所述步骤(2)中求解得到调制度系数后根据如下公式归一化至区间[0,1]：

I″_z＝(I″-I″_min)/(I″_max-I″_min)

其中I″_z为归一化之后的调制度系数，I″_max，I″_min分别表示求取的调制度系数I″的最大和最小值。

4.根据权利要求3所述的一种光栅投影三维测量系统的误差校正方法，其特征在于，所述步骤(3)中通过设定一个分割阈值，将归一化后的调制度系数低于设定阈值的像素点确定为需要处理的像素点；所述分割阈值取值范围为0.25～0.4。

5.根据权利要求3所述的一种光栅投影三维测量系统的误差校正方法，其特征在于，所述步骤(4)中通过设定一个最小阈值，将每两个像素点对应的调制度系数差值小于设定阈值的像素点划分在同一类中；所述最小阈值取值范围为0.008～0.012。

6.根据权利要求1所述的一种光栅投影三维测量系统的误差校正方法，其特征在于，所述步骤(6)中主值相位的计算公式为：

$\mathrm{\φ}=arctan\[\frac{\underset{n=1}{\overset{N}{\Σ}}{I}_n^{r}sin(2\mathrm{\π}n/N)}{\underset{n=1}{\overset{N}{\Σ}}{I}_n^r\cos (2\mathrm{\π}n/N)}\]$

其中，为校正后的第n幅图像的灰度值。