结构光三维重建（四步相移&三频外差）MATLAB实战：从条纹生成到相位展开

1. 结构光三维重建基础入门

结构光三维重建就像给物体拍"立体照片"的过程。想象一下，你用手电筒斜着照向墙面，看到的亮斑形状会随着墙面凹凸变化——这就是最基础的结构光原理。在实际应用中，我们会用投影仪代替手电筒，用相机记录变形后的光斑图案，再通过计算还原出物体的三维形状。

MATLAB在这个领域特别受欢迎，主要因为三个原因：一是矩阵运算能力强，处理图像数据得心应手；二是可视化功能完善，能实时查看相位图、三维点云等结果；三是代码可读性好，特别适合算法验证阶段。我刚开始接触这个领域时，就是通过MATLAB的demo程序理解了相位计算的核心逻辑。

整套流程可以分解为五个关键步骤：

1. 生成多组正弦条纹图案（就像制作特殊的条形码）
2. 用投影仪将这些图案投射到物体表面
3. 相机拍摄变形后的条纹图像
4. 通过相位计算得到"包裹相位"（相当于加密的高度信息）
5. 最后进行相位展开，得到真正的三维坐标

初学者最容易混淆的是"包裹相位"和"绝对相位"的区别。包裹相位就像钟表指针——虽然显示3点，但可能是15点、27点...因为每12小时就循环一次。而我们的任务就是要解开这个循环，找到真正的"绝对相位"。

2. 正弦条纹生成实战

生成高质量的正弦条纹是整套流程的第一步，也是很多新手容易翻车的地方。先来看一个典型参数设置：

matlab复制width = 1240;  % 图像宽度
height = 1110; % 图像高度
A = 0.5;       % 条纹对比度
B = 0.5;       % 条纹振幅
WaveLengthArr = [70, 64, 59];  % 三个频率的周期(像素)
nStepPS = 4;   % 四步相移

这里有个关键点：WaveLengthArr中的70/64/59不是随便选的数字。根据我的实测经验，这三个数应该互为质数，且满足1/70 - 1/64 ≈ 1/6的关系。这样的频率组合能确保后续三频外差时，解包裹过程更加稳定。

条纹生成的核心函数如下：

matlab复制function [FringeSet] = GetFringeSet(height, width, nStepPS, WaveLengthArr, ObjectHeight,A,B)
    FringeSet = cel