多频外差结构光三维测量与绝对相位恢复系统

资 源 简 介

本项目旨在开发一套完整的结构光三维测量算法库，核心聚焦于利用多频外差技术进行高精度的相位解包裹。系统主要功能涵盖从条纹生成、图像采集、相位计算到三维重建的全过程。首先，程序能够生成多组不同频率（如70、64、59等）的四步或N步相移正弦条纹图样。其次，利用相移算法处理相机采集到的变形条纹图像，计算出高精度的截断相位（包裹相位）。核心模块通过多频外差原理，将不同高频条纹的相位信息进行叠加差拍，合成等效波长极大的相位图，从而在不依赖空间邻域信息的情况下，将[-pi, pi]区间的截断相位展开为连续的绝对相位，有效解决相位模糊问题，特别适合测量表面具有不连续台阶或复杂几何形状的物体。此外，该项目还包含图像预处理（如高斯滤波、掩膜处理）、相位误差补偿以及基于标定参数的相位-深度映射功能，最终输出高精度的三维点云模型，适用于工业零部件检测、逆向工程及生物特征识别等场景。

详 情 说 明

基于多频外差原理的高精度绝对相位恢复与三维测量系统

项目介绍

本项目实现了一套完整的结构光三维测量仿真系统。该系统基于Matlab开发，核心旨在演示如何利用多频外差（Multi-frequency Heterodyne）技术解决相位解包裹（Phase Unwrapping）难题，从而实现从高频截断相位到连续绝对相位的恢复，并最终重建出被测物体的三维点云模型。

程序完整模拟了结构光测量的全流程，包括虚拟三维场景生成、条纹投影与图像采集模拟、四步相移法相位计算、多频外差相位解包裹以及相位-深度映射。该代码特别适用于研究和理解高精度三维光学测量中的核心算法逻辑。

功能特性

• 全流程仿真：无需硬件设备，通过数学模型生成虚拟物体（类高斯曲面）并模拟投影机和相机的成像过程。
• 多频条纹生成：支持生成互质频率组合（代码中默认设置为70、64、59）的正弦条纹，用于实现高精度的外差解相。
• 模拟图像采集：模拟了环境光、物体反射率、条纹调制以及随机噪声对采集图像的影响。
• 四步相移算法：采用标准的N步（默认为4步）相移算法，点对点计算高精度的截断相位和调制度。
• 多频外差解包裹：核心算法模块，利用三个频率的拍频效应，级联恢复绝对相位，有效解决相位模糊问题。
• 三维重建与可视化：基于简化的线性标定模型将绝对相位转换为深度信息，并提供包括2D相位图、截面线对比及3D点云在内的多维度可视化结果。

算法实现细节与逻辑

本项目在主程序中依次执行以下六个主要步骤，逻辑严格对应代码实现：

1. 系统参数初始化与场景生成

程序首先定义了图像分辨率（800x600）和核心测量参数，包括三个互质的条纹频率（$f_1=70, f_2=64, f_3=59$）和相移步数（4步）。 同时，利用混合高斯函数生成了一个具有复杂几何形状的虚拟物体高度图（GroundTruthZ），并生成了圆形掩膜（Mask），用于模拟物体在视场中的有限范围。

2. 条纹投影与图像采集模拟

利用数学模型模拟投影过程。对于每一个频率，程序生成4幅相移条纹图。 实现逻辑：结合了基础载波相位（由频率决定）和物体高度引起的相位调制（由scan_sensitivity系数决定）。为了模拟真实环境，图像生成过程中叠加了：

• 环境背景光（A）
• 表面反射调制度（B）
• 高斯随机噪声
• 背景区域的低对比度处理

3. N步相移相位计算

对采集到的图样应用并行的四步相移算法。 算法细节：通过计算四幅图像强度的代数运算（atan2函数），获取每个像素点的截断相位（Wrapped Phase），范围限制在 $(-pi, pi]$。同时计算调制度（Modulation），并据此生成有效像素掩膜，过滤掉背景噪声区域（调制度阈值设为0.05）。

4. 多频外差绝对相位恢复（核心）

这是系统的核心算法模块，用于将截断相位展开为连续相位。 级联逻辑：

1. 一级外差：计算 $f_1$ 与 $f_2$、$f_2$ 与 $f_3$ 之间的拍频相位（$phi_{12}, phi_{23}$）及对应的拍频频率（$f_{12}, f_{23}$）。通过相位减法函数 subPhase 确保差值相位被正确归一化。
2. 二级外差：计算 $f_{12}$ 与 $f_{23}$ 之间的二级拍频相位（$phi_{123}$）及频率（$f_{123}$）。理论上 $f_{123}$ 极小，对应波长覆盖全视场，因此 $phi_{123}$ 可视为无截断的绝对相位基准。
3. 相位展开：采用由低频向高频逐级引导的策略。

* 利用 $phi_{123}$ 展开 $phi_{12}$。 * 利用展开后的 $phi_{12}$ 展开主频率 $phi_1$（最高频率），从而获得具备高精度的绝对相位图。

5. 相位转深度（三维重建）

将恢复的绝对相位转换为空间几何坐标。 实现逻辑：

1. 计算参考平面相位（RefPhase），即不存在物体时的载波相位。
2. 计算相位差：物体绝对相位减去参考相位。
3. 深度映射：应用简化的线性模型 $Z = K_{sys} times DeltaPhi$（其中 $K_{sys}=2.5$），将相位差直接映射为高度数据。
4. 点云合成：结合像素坐标网格（X, Y）与计算出的 Z 值，生成最终的三维点云数据。

6. 结果可视化

程序最后调用可视化子函数，在一个窗口中分块展示：

• 模拟采集的原始条纹图像。
• 主频率下的截断相位图（锯齿状）。
• 中心行截面线对比图（直观展示截断相位与展开后相位的关系）。
• 最终合成的绝对相位伪彩图。
• 经过渲染的三维点云网格模型（应用了光照效果）。

关键函数说明

• generateVirtualObject：构建虚拟被测物，利用多个高斯分布函数的叠加模拟平滑但起伏变化的表面。
• simulateProjectedPattern：物理成像仿真器，负责将高度信息编码进相移条纹中，并引入噪声模型。
• phaseShiftAlgorithm：标准相移解相器，输入图像序列，输出截断相位和调制度图。
• multiFrequencyHeterodyne：多频外差逻辑的具体实现者，负责执行相位差拍和级联解包裹运算。
• subPhase：辅助函数，处理相位减法中的 $2pi$ 跳变问题，保证拍频相位落在对应的线性区间内。
• unwrapDualFreq：辅助函数，实现基于频率比（Ratio）的相位展开公式，利用低频无模糊相位辅助高频相位确定级次。

系统要求

• Matlab R2016a 或更高版本（代码使用基础矩阵运算，无特殊工具箱依赖，但在较低版本中需确保 imagesc 或 colormap 等绘图函数的兼容性）。
• 足够的内存以存储多组高分辨率图像矩阵（建议 4GB 以上）。

使用方法

1. 确保Matlab的工作路径已设置为包含本脚本的文件夹。
2. 直接运行 main 函数。
3. 控制台将输出各阶段的处理进度及耗时信息。
4. 运行结束后，会弹出一个综合图形窗口显示仿真与重建结果。