🔄 八步相移基本原理
八步相移技术通过投射八幅具有不同相位的正弦条纹图案，利用相位编码信息精确计算物体表面的三维形状。每步相移的光强信号可以表示为：
I_k(x) = A(x) + B(x)cos[φ(x) + δ_k]
参数说明：
• A(x)：背景光强（直流分量）
• B(x)：调制幅度
• φ(x)：被测物体引起的相位调制
• δ_k：第k步的相移量，k = 1,2,3,...,8
八步相移的相位分布
📊 调制对比度的推导过程
步骤1：计算背景光强A(x)
由于八个相移量均匀分布在2π范围内，余弦函数的平均值为零：
A(x) = (1/8) Σ[k=1 to 8] I_k(x)
详细推导：
展开求和：
Σ[k=1 to 8] I_k = Σ[k=1 to 8] [A + B·cos(φ + δ_k)]
= 8A + B·Σ[k=1 to 8] cos(φ + δ_k)
计算余弦和：
Σ[k=1 to 8] cos(φ + δ_k) = cos(φ)·Σcos(δ_k) - sin(φ)·Σsin(δ_k)
由于δ_k均匀分布，两个求和均为0
最终结果：
A(x) = (1/8)·8A = A
步骤2：使用傅里叶分析计算调制幅度
利用离散傅里叶变换的原理，计算N₁和N₂分量：
N₁ = Σ[k=1 to 8] I_k · sin(δ_k)
N₂ = Σ[k=1 to 8] I_k · cos(δ_k)
N₁分量计算及为什么等于4B·sin(φ)的详细推导：
第一步：N₁的定义和展开
N₁ = Σ[k=1 to 8] I_k · sin(δ_k)
其中 I_k = A + B·cos(φ + δ_k)
展开得到：
N₁ = Σ[k=1 to 8] [A + B·cos(φ + δ_k)] · sin(δ_k)
= A·Σ[k=1 to 8] sin(δ_k) + B·Σ[k=1 to 8] cos(φ + δ_k)·sin(δ_k)
第二步：计算第一项 A·Σsin(δ_k)
δ_k = (k-1)·π/4，k = 1,2,3,...,8
即：δ₁=0, δ₂=π/4, δ₃=π/2, δ₄=3π/4, δ₅=π, δ₆=5π/4, δ₇=3π/2, δ₈=7π/4
Σsin(δ_k) = sin(0) + sin(π/4) + sin(π/2) + sin(3π/4) + sin(π) + sin(5π/4) + sin(3π/2) + sin(7π/4)
= 0 + √2/2 + 1 + √2/2 + 0 + (-√2/2) + (-1) + (-√2/2)
= 0 + √2/2 + 1 + √2/2 + 0 - √2/2 - 1 - √2/2 = 0
因此：A·Σsin(δ_k) = A·0 = 0
第三步：计算第二项 B·Σcos(φ + δ_k)·sin(δ_k)
使用三角恒等式：cos(φ + δ_k) = cos(φ)cos(δ_k) - sin(φ)sin(δ_k)
B·Σcos(φ + δ_k)·sin(δ_k)
= B·Σ[cos(φ)cos(δ_k) - sin(φ)sin(δ_k)]·sin(δ_k)
= B·cos(φ)·Σcos(δ_k)sin(δ_k) - B·sin(φ)·Σsin²(δ_k)
第四步：计算Σcos(δ_k)sin(δ_k)
利用恒等式：cos(x)sin(x) = (1/2)sin(2x)
Σcos(δ_k)sin(δ_k) = (1/2)Σsin(2δ_k)
= (1/2)[sin(0) + sin(π/2) + sin(π) + sin(3π/2) + sin(2π) + sin(5π/2) + sin(3π) + sin(7π/2)]
= (1/2)[0 + 1 + 0 + (-1) + 0 + 1 + 0 + (-1)] = 0
第五步：计算Σsin²(δ_k)
Σsin²(δ_k) = sin²(0) + sin²(π/4) + sin²(π/2) + sin²(3π/4) + sin²(π) + sin²(5π/4) + sin²(3π/2) + sin²(7π/4)
= 0 + (√2/2)² + 1² + (√2/2)² + 0² + (-√2/2)² + (-1)² + (-√2/2)²
= 0 + 1/2 + 1 + 1/2 + 0 + 1/2 + 1 + 1/2
= 4
第六步：得出最终结果
N₁ = B·cos(φ)·0 - B·sin(φ)·4
= 0 - 4B·sin(φ)
= -4B·sin(φ)
注意：通常我们取绝对值或者在后续计算中通过平方消除符号，所以表示为：
|N₁| = 4B·sin(φ)
物理意义解释：
• 系数4来源于8个均匀分布相移角度的sin²值之和
• sin(φ)项反映了被测物体相位对正弦分量的调制
• 这个结果说明N₁包含了调制幅度B和相位φ的信息
• 通过N₁和N₂可以完全重构出B和φ
N₂分量计算（类似推导）：
同理可以推导出：
N₂ = Σ[k=1 to 8] I_k · cos(δ_k) = 4B·cos(φ)
因此：
• N₁ = -4B·sin(φ)（考虑符号）
• N₂ = 4B·cos(φ)
• √(N₁² + N₂²) = 4B
• B = √(N₁² + N₂²) / 4
步骤3：计算调制幅度B(x)及除以4的详细解释
B(x) = (1/4)√(N₁² + N₂²)
为什么要除以4？深度解析：
1. 傅里叶变换理论基础：
在连续傅里叶变换中，一个周期函数f(t) = A + B·cos(ωt + φ)的傅里叶系数为：
• 基频余弦分量系数：B/2
• 基频正弦分量系数：B/2
• 幅度：√[(B/2)² + (B/2)²] = B/√2 × √2 = B
2. 离散傅里叶变换(DFT)修正：
对于N点DFT，第k个频率分量的计算公式为：
X[k] = Σ[n=0 to N-1] x[n] · e^(-j2πkn/N)
对于8点相移数据，实际计算得到的N₁和N₂是：
• N₁ = Σ I_k · sin(δ_k) = 4B · sin(φ)
• N₂ = Σ I_k · cos(δ_k) = 4B · cos(φ)
3. 系数4的来源：
在8步相移中，每个三角函数值的平方和为：
Σ[k=1 to 8] sin²(δ_k) = Σ[k=1 to 8] cos²(δ_k) = 4
因此：
√(N₁² + N₂²) = √[(4B·sin(φ))² + (4B·cos(φ))²]
= √[16B²(sin²(φ) + cos²(φ))]
= √[16B²] = 4B
所以：B = √(N₁² + N₂²) / 4
4. 物理意义：
除以4是为了从DFT的原始输出中提取真实的物理调制幅度，消除由于采样点数和三角函数求和引入的比例因子。
数学验证过程：
步骤1：验证三角函数和的性质
对于δ_k = k·π/4 (k=0,1,2,...,7)：
Σ sin²(δ_k) = sin²(0) + sin²(π/4) + sin²(π/2) + ... + sin²(7π/4)
= 0 + 1/2 + 1 + 1/2 + 0 + 1/2 + 1 + 1/2 = 4
Σ cos²(δ_k) = cos²(0) + cos²(π/4) + cos²(π/2) + ... + cos²(7π/4)
= 1 + 1/2 + 0 + 1/2 + 1 + 1/2 + 0 + 1/2 = 4
步骤2：推导N₁和N₂的表达式
I_k = A + B·cos(φ + δ_k)
N₁ = Σ I_k · sin(δ_k)
= Σ [A + B·cos(φ + δ_k)] · sin(δ_k)
= A·Σ sin(δ_k) + B·Σ cos(φ + δ_k)·sin(δ_k)
= 0 + B·[cos(φ)·Σ cos(δ_k)sin(δ_k) - sin(φ)·Σ sin²(δ_k)]
= B·[0 - sin(φ)·4] = -4B·sin(φ)
类似地：N₂ = 4B·cos(φ)
步骤3：计算最终的调制幅度
√(N₁² + N₂²) = √[(-4B·sin(φ))² + (4B·cos(φ))²]
= √[16B²·sin²(φ) + 16B²·cos²(φ)]
= √[16B²(sin²(φ) + cos²(φ))]
= √[16B²] = 4B
因此：B = √(N₁² + N₂²) / 4
步骤4：计算调制对比度M(x)
M(x) = B(x) / A(x) = √(N₁² + N₂²) / (4A(x))
🔍 数值示例分析
示例数据
假设测得八步相移的光强值：
I₁ = 120, I₂ = 145, I₃ = 160, I₄ = 145
I₅ = 80, I₆ = 55, I₇ = 40, I₈ = 55
完整计算过程：
1. 计算背景光强：
A = (120 + 145 + 160 + 145 + 80 + 55 + 40 + 55) / 8 = 100
2. 计算N₁：
N₁ = (√2/2)(145 + 145 - 55 - 55) + (160 - 40)
= (√2/2) × 180 + 120
= 127.3 + 120 = 247.3
3. 计算N₂：
N₂ = (120 - 80) + (√2/2)(145 - 145 - 55 + 55)
= 40 + 0 = 40
4. 计算调制幅度：
B = (1/4)√(247.3² + 40²) = (1/4)√(64,237) ≈ 63.4
5. 计算调制对比度：
M = B/A = 63.4/100 = 0.634
结果分析：调制对比度为0.634，属于良好质量范围，可以保证较高的测量精度。
⚙️ 实际应用中的考虑因素
1. 噪声处理
高斯噪声：八步相移的多点采样能有效降低随机噪声影响
椒盐噪声：可通过中值滤波预处理
系统噪声：八步算法本身具有噪声抑制特性
2. 系统标定
相移标定：确保每步相移准确为π/4
强度标定：消除投影仪和相机的非线性响应
几何标定：建立像素坐标与世界坐标的关系
3. 质量评估指标
调制对比度阈值：
• M > 0.7：优秀质量，精度可达λ/100
• 0.5 < M ≤ 0.7：良好质量，精度约λ/50
• 0.3 < M ≤ 0.5：可接受质量，需要处理
• M ≤ 0.3：质量较差，建议重新测量