相位偏折算法2.5d成像系统 相位偏折的算法实现代码，基于8张原图或者4张，计算出形状图，镜面反射分量，漫反射分量等 py代码和cpp代码都实现了，联系请备注需要哪一种，默认发cpp的

基于对提交的C++源代码文件的详细分析，以下是忠实于原始代码实现的功能说明：

程序功能说明：基于四步相移法的相位偏折测量处理

程序概述

本程序是一个完整的相位偏折测量处理流水线，通过处理8张输入图像（分为X和Y两个方向各4张），计算出相位、梯度、反射分量等多种光学测量结果。

核心处理流程

1. 双方向独立处理

程序分别处理X方向和Y方向的图像序列：

• X方向：使用图像5.jpg、6.jpg、7.jpg、8.jpg
• Y方向：使用图像1.jpg、2.jpg、3.jpg、4.jpg

2. 四步相移法计算包裹相位

每个方向的处理遵循相同的数学原理：

// 计算正弦和余弦分量
SinMat = SinMat + image1 * std::sin(1 * PI / 2.0);
CosMat = CosMat + image1 * std::cos(1 * PI / 2.0);
// ... 对4张图像重复此过程

// 计算包裹相位
float pahval = std::atan2(cosval, sinval);

3. 一维相位解包裹

• X方向：沿图像行方向进行解包裹
• Y方向：沿图像列方向进行解包裹

解包裹逻辑基于相邻像素的相位跳变检测：

if (previousPixlVal - val > 2) PhaIndex++;
if (previousPixlVal - val < -2) PhaIndex--;

4. 多维度结果计算

每个方向生成5种输出结果：

• 包裹相位 (PhaMat)：原始相位计算结果
• 解包裹相位 (DecodePhaMat)：连续相位信息
• 梯度图 (grad_x)：使用Sobel算子计算的相位梯度
• 镜面反射分量 (SRC)：基于图像差分的反射特性计算
  cpp
float val = cv::sqrt(val31 val31 + val42 val42) / 2;

• 平均图像 (AverageImage)：四张输入图像的平均值，表征漫反射分量

5. 结果融合与后处理

在main函数中完成：

• 将X和Y方向的解包裹相位、梯度图、镜面反射分量进行平均融合
• 对融合后的相位和梯度图应用9×9均值滤波平滑处理
• 输出最终的完整结果文件

输出文件体系

程序生成完整的中间和最终结果：

• 方向特定结果：pha.tiff, D_pha.tiff, Shap.tiff, SRC.tiff, AverageImage.tiff
• 融合最终结果：fullDpha.tiff, fullShape.tiff, full_SRC.tiff

技术特点

1. 硬编码路径：所有输入输出路径在代码中固定
2. 单精度浮点：使用CV_32F类型保证计算精度
3. 尺寸验证：在融合前检查矩阵尺寸和类型一致性
4. 控制台反馈：通过std::cout输出处理状态信息

该程序实现了一个完整的相位偏折测量分析流程，从原始图像采集到最终的三维特性表征，为光学测量和表面形貌分析提供了基础数据支撑。

相位偏折算法2.5d成像系统 相位偏折的算法实现代码，基于8张原图或者4张，计算出形状图，镜面反射分量，漫反射分量等 py代码和cpp代码都实现了，联系请备注需要哪一种，默认发cpp的