一种基于分区域偏振调控的高反光表面三维重建方法

本发明属于结构光三维测量，尤其涉及一种基于分区域偏振调控的高反光表面三维重建方法，适用于具有高反光表面的金属、陶瓷、涂层等工件的高精度三维形貌重建。

背景技术：

1、结构光三维测量技术因其非接触、高精度、全场快速获取等优势，已成为先进制造、质量检测与逆向工程等领域的关键技术之一。其基本原理是通过投影仪向被测表面投射编码条纹图案，相机同步采集变形条纹图像，经相位解算与系统标定获得被测表面的三维形貌。然而，当被测表面存在金属量块、抛光工件、光滑涂层等高反光区域时，镜面反射导致局部图像饱和或条纹对比度丧失，严重影响相位解算的准确性与完整性，成为制约结构光三维测量精度与适用性的瓶颈问题。

2、针对高反光表面测量难题，现有技术主要分为三类：多曝光高动态范围成像、自适应投影强度调节和偏振抑制。多曝光方法通过采集多组不同曝光时间的图像进行融合，但采集时间长，难以满足在线测量需求；自适应投影方法通过逐像素调整投影强度，但受限于投影仪的动态范围和像素级独立控制的复杂性，实际实施困难。偏振抑制方法利用镜面反射光的部分偏振特性，通过在光路中加入偏振器件并调节偏振角以抑制镜面反射分量，具有硬件结构简单、无需多次采集的优势，是当前研究热点。

3、然而，现有偏振方法大多采用单一全局偏振角，未考虑镜面反射区与漫反射区在最佳偏振角上的本质差异——抑制镜面反射需靠近消光方向，而增强漫反射信号需靠近透光方向。若将偏振角固定于消光方向，虽能抑制高反光区饱和，但会导致低反射区域条纹信号衰减、信噪比不足；若固定于透光方向，则高反光区仍出现过曝。单一全局偏振角无法同时满足两类区域的优化需求，导致抑制效果与条纹质量难以兼顾。

4、为此，本发明提出一种分区域自适应偏振调控方法，根据表面反射特性的空间差异，自动划分高反光区与低反射区，并为各区域独立计算最佳偏振角，以满足高反光区域与低反射区域各自的最佳成像需求。

技术实现思路

1、为了克服现有技术的不足，本发明的目的在于提供一种基于分区域偏振调控的高反光表面三维重建方法，通过物理驱动的区域分割与逆向寻优策略，在保证低反射区条纹信噪比的同时有效抑制高反光区镜面反射，实现高反光表面全区域高完整度、高精度的三维重建。

2、本发明采用如下技术方案：一种基于分区域偏振调控的高反光表面三维重建方法，包括以下步骤：

3、步骤s1：偏振结构光图像采集。搭建单目偏振结构光成像系统，在投影仪前端加装起偏器，在相机前端加装检偏器。投影仪向被测表面投射正弦条纹图案，相机依次采集0°、45°、90°、135°四个偏振角度下的结构光条纹图像，图像灰度值分别用、、、表示。

4、步骤s2：偏振特征计算与鲁棒性增强。根据四个偏振角度下的图像灰度值，解算stokes参量，各参量按以下方式计算：

5、

6、式中，表示总光强；是光的水平与垂直偏振分量之间的强度差；是光的+45°和-45°偏振分量之间的强度差。

7、进而求解各像素的平均光强，偏振调制幅度，偏振相位角,线偏振度。为抑制低反射区域散粒噪声的干扰，构建鲁棒线偏振度，用代表图像像素的横、纵坐标，公式如下：

8、

9、仅对满足r>1且a>10的有效像素计算线偏振度，无效像素置零。同时对计算结果进行值域约束，将最大值截断为1，以规避噪声引发的超物理范围异常值。

10、步骤s3：物体掩膜提取。对四个偏振角度图像逐像素取最大值，生成最大光强图，真实反映场景有效亮度分布。对最大光强图进行全局阈值分割，并执行形态学开运算与闭运算，去除孤立噪点、填充孔洞，获得粗掩膜。在粗掩膜基础上，采用otsu自适应阈值分割结合形态学细化，生成精细工件掩膜，严格限定后续聚类计算域，排除背景干扰。

11、步骤s4：基于fcm的物理驱动区域分割。在精细工件掩膜约束下，以该掩膜内平均光强a的最大值对平均光强a做归一化，得到归一化光强，即每一像素的归一化光强等于其平均光强除以精细工件掩膜内平均光强的最大值，联合鲁棒线偏振度构成二维特征空间。在该特征空间内，采用模糊c均值聚类将像素自适应划分为高反光区和低反射区两类。根据聚类中心平均光强的大小差异，自动标记类别归属——平均光强较高的聚类中心对应高反光区，较低的对应低反射区。最终将精细工件掩膜外的背景有效区域统一归入低反射类别，输出高反光区掩膜与低反射区掩膜，完成全域精准分区。

12、步骤s5：低反射区透光最大化优化。低反射区域以漫反射为主，反射光强较低，首要目标是最大化透射光强以提升条纹信噪比。根据马吕斯定律，出射光强，当偏振角时达到最大值。为抑制单像素噪声，对低反射区掩膜内满足信噪比要求的有效像素，分别统计stokes分量的区域算术平均值，计算低反射区等效偏振相位角，则低反射区最佳偏振角，并通过模运算将其映射至区间。

13、步骤s6：高反光区调制度响应模型与基准角确定。高反光区镜面反射强烈，需在抑制镜面反射与保证条纹可解调之间寻求平衡。首先，分别在四个偏振角下采集一组四步相移条纹图像，用分别代表这四步相移条纹图像的像素灰度值。四步相移的包裹相位、背景光强、条纹调制幅度、条纹调制度。在高反光区域内取空间平均，得到四个平均调制度值。由于光强随偏振角按马吕斯定律呈余弦变化，条纹调制度亦以180°为周期变化，据此建立调制度-偏振角的余弦响应模型：

14、

15、式中，为调制度响应曲线的直流分量（基线值），为调制度响应曲线的调制幅度，为调制度响应曲线的相位偏移；参数由stokes方式求解。

16、同时，利用高反光区平均stokes分量计算，及其透光方向和消光方向。引入低反射区与高反光区的中值光强比：

17、

18、式中，为取中值运算；分别为低反射区、高反光区的像素平均光强，截断至(0,1]区间。在透光方向与消光方向之间按线性插值，定义基准偏振角：

19、

20、式中，表示对180°取模，将角度映射至区间；当高反光区相对越亮（越小），基准角越靠近消光方向；反之则靠近透光方向。

21、步骤s7：自适应窗口搜索与最佳偏振角确定。将透光方向、消光方向、基准角以及调制度响应曲线统一映射至连续空间。以为中心、向两侧各扩展（可根据实际测量装置调整）定义搜索窗口，并截断于透光端与消光端之间。在窗口内以0.1°为步长，利用拟合的调制度响应模型计算各角度对应的调制度预测值。统计窗口内调制度的极大值和极小值，定义可用性阈值：

22、

23、式中，k为阈值系数；作为优选，阈值系数k取值为0.10～0.20，进一步优选为0.15，该阈值仅剔除窗口内底部15%的极低调制度区域。搜索遵循"消光优先"原则：将窗口内角度按从消光端到透光端的方向排序，依次检查是否满足，选取首个满足条件的角度作为高反光区最佳偏振角。若整个窗口内无角度满足条件，则退化选取窗口内调制度最大的角度。该退化机制确保算法在任何极端情况下均能返回可用偏振角。

24、步骤s8：分区采集与点云融合。分别在高反光区最佳偏振角和低反射区最佳偏振角下，采集完整的三频四步相移条纹图像序列。经相位解算、多频外差展开得到绝对相位，再结合偏振结构光成像系统的标定参数重建为两片三维点云。由于相机与投影仪固定，两片点云天然处于同一坐标系。以偏振结构光成像系统的单像素空间分辨率（例如0.1mm）为容差，对两片点云进行体素格栅合并——将整个空间划分为边长为容差的正方体体素，对同一体素内的多个点进行平均或保留距中心最近的点，最终输出覆盖完整、点云密度均匀的被测表面三维形貌。

25、进一步地，步骤s2中鲁棒线偏振度的信噪比阈值可根据实际相机噪声水平自适应调整。

26、进一步地，步骤s4中模糊c均值聚类的类别数固定为2，模糊加权指数取2，迭代终止条件为最大迭代次数100次或隶属度变化小于1e-5。

27、进一步地，步骤s7中搜索窗口的半宽度应使所述高反光区最佳偏振角与消光方向之间留有裕量，优选为10°，步长可取0.1°~0.5°；阈值系数k取值为0.10～0.20，进一步优选为0.15。

28、进一步地，步骤s8中体素格栅合并的容差可设置为偏振结构光成像系统标定平均重投影误差对应的空间分辨率。

29、相比现有技术，本发明的有益效果在于：

30、1.本发明的一种基于分区域偏振调控的高反光表面三维重建方法，其摒弃了传统单一全局偏振角策略，通过fcm聚类自动划分高反光区与低反射区，并分别以透光最大化和消光优先为目标独立计算最佳偏振角，使两类区域均工作于最佳状态，克服了全局偏振角无法兼顾两类区域需求的缺陷。

31、2.本发明的一种基于分区域偏振调控的高反光表面三维重建方法，采用归一化光强与鲁棒线偏振度构建二维特征空间，两类区域在该空间中呈对角分布，具有天然可分性。鲁棒偏振度特征不仅提高了金属高反光场景的分割可靠性，还为方法向非金属高反光物体的泛化提供了物理基础。

32、3.本发明的一种基于分区域偏振调控的高反光表面三维重建方法，针对高反光区，建立了条纹调制度对偏振角的余弦响应模型，并设计了基于窗口动态范围的自适应阈值，在消光优先原则下确保所选偏振角对应的调制度不低于可靠解调水平，实现了镜面反射抑制与条纹质量之间的动态平衡，避免过度消光导致相位解算失败。

33、4.本发明的一种基于分区域偏振调控的高反光表面三维重建方法，可自动完成区域分割、基准角计算、窗口搜索和偏振角输出，无需人工干预；当搜索窗口内无角度满足调制度约束时，退化机制自动选取调制度最大的角度，保证算法在极端条件下仍能返回可用偏振角。

34、5.本发明的一种基于分区域偏振调控的高反光表面三维重建方法，其点云融合策略简单有效，两片点云天然同坐标系，仅需简单体素格栅合并即可获得完整形貌，无需额外配准，计算效率高。本发明在金属量块阶梯、叶片等典型高反光工件上均能有效抑制镜面反射，能够实现形貌完整、无明显空洞的三维重建，适用于复杂曲面及精细纹理工件的三维测量。