激光散斑图数值模拟与仿真系统

资 源 简 介

该项目是一个专门用于光测力学、数字图像相关法（DIC）和统计光学研究的专业散斑图合成工具。 其核心功能是利用数值计算方法在计算机上精确模拟物理散斑的形成过程，能够生成符合特定统计规律的随机散斑图案。 系统实现了基于复高斯随机场模型的散斑生成算法，通过在频率域控制随机相位的分布并通过快速傅里叶逆变换（IFFT）回到空间域，从而精确控制散斑的平均颗粒直径、对比度及亮度分布。 此外，系统支持生成主动散斑，即通过对特定形状（如圆形、卵形）的子斑点进行随机排布来模拟喷漆散斑的效果。 为了方便实验仿真，该工具还集成了数字图像旋转、平移、拉伸等几何变换功能，可以根据预设的位移场生成变形后的散斑图对。 应用场景涵盖了DIC位移测量精度验证、光学系统传递函数评估、散斑干涉计量模拟以及深度学习散斑图像数据集的构建。 该系统还提供质量评价模块，能够自动分析生成的散斑图像的平均灰度梯度、子区熵值及自相关函数，确保生成的模拟图像满足实验力学的严苛要求。

详 情 说 明

基于MATLAB的激光散斑图数值模拟与仿真系统

项目介绍

本系统是一款专为光测力学、数字图像相关法（DIC）以及统计光学研究设计的专业级散斑图像合成与仿真工具。在实验力学领域，高质量的散斑图案是确保位移与应变测量精度的前提。该系统通过数值计算模拟物理散斑的产生机理，不仅能够生成符合统计学规律的激光散斑，还能模拟人工喷漆的主动散斑，并支持复杂的亚像素级图像变形仿真，为DIC算法的验证和光学系统的评估提供标准数据集。

功能特性

1. 统计散斑合成：基于复高斯随机场（CGRF）模型，利用频域滤波技术生成具有精确统计特性的激光散斑图，可控散斑平均直径与对比度。
2. 主动散斑模拟：通过随机排布特定几何形状（如圆形子斑点）的方法，真实还原实验中喷漆散斑的视觉特征，支持密度与粒径波动的自定义。
3. 高精度形变仿真：实现包含平移、旋转及单轴应变在内的组合变形场模拟，利用双三次插值确保变形后的图像具备亚像素级精度。
4. 质量评价体系：集成平均灰度梯度（MIG）、信息熵及自相关函数分析功能，从数字化角度量化散斑图像对DIC计算的适用性。
5. 结果可视化与导出：系统自动生成包含参考图、变形图、位移残差及质量分析指标的可视化报表，并支持高位深TIFF格式图像导出。

使用方法

1. 配置环境：确保开发环境已安装MATLAB及图像处理工具箱（Image Processing Toolbox）。
2. 参数设置：在脚本起始位置修改图像分辨率、散斑直径、覆盖率、背景噪声等物理参数，并根据仿真需求设定预设的位移量、旋转角度和应变系数。
3. 执行仿真：运行主函数，系统将依次执行频域随机场生成、主动散斑合成、坐标映射变换、亚像素插值及噪声叠加。
4. 结果分析：程序运行结束后将弹出对比窗口，实时显示生成的散斑图及三维自相关函数分布图，同时在命令行或图形界面中输出质量评价指标。
5. 提取数据：生成的散斑图像（参考图与变形图）将自动保存至当前工作目录下，便于后续DIC软件处理。

系统要求

1. 软件环境：MATLAB R2018b 或更高版本。
2. 必备工具箱：Image Processing Toolbox（用于梯度计算、信息熵分析及图像调整）。
3. 硬件建议：建议内存 8GB 以上，以支持高分辨率图像的快速傅里叶变换与矩阵运算。

核心功能实现逻辑说明

1. 统计散斑生成算法（CGRF模型）

系统首先在频率域生成均匀分布的随机相位场，模拟相干光的相位随机叠加。随后构造一个高斯低通滤波器，其标准差与预设的散斑直径成反比，通过频域滤波控制散斑的颗粒大小。通过快速傅里叶逆变换（IFFT）将信号转回空间域，取模平方获得强度分布。最后通过归一化调整和对比度增强（imadjust）确保图像灰度分布符合实验需求。

1. 主动散斑生成逻辑

该模块模拟人工制备散斑的过程。系统计算达成指定覆盖率所需的子斑点数量，并在图像空间内随机生成中心坐标。每个子斑点具有随机波动的半径和亮度亮度。通过构建圆形掩模（Mask），将子斑点叠加到背景画布上，并进行亮度截断处理，防止局部过曝光。

1. 形变仿真与亚像素插值

形变模拟采用逆向映射法。系统构建原始图像坐标系后，应用包含平移向量、旋转矩阵和应变系数的仿射变换公式计算变形后的坐标映射关系。核心步骤采用双三次插值（Bicubic Interpolation）在非整数像素点上重塑灰度值，从而消除离散化带来的误差，模拟真实的连续形变过程。此外，系统会在图像合成后加入高斯白噪声，以模拟相机传感器的感光噪声。

1. 图像质量评估指标

平均灰度梯度（MIG）：通过计算图像在x和y方向的一阶导数各向异性强度，反映散斑的边缘尖锐度，是衡量DIC子区匹配稳定性的核心指标。 信息熵（Entropy）：表征图像的信息丰富程度，熵值越高代表散斑的特征越显著。 自相关函数：选取图像中心子区进行二维互相关运算，生成的3D峰值图越尖锐，说明散斑的随机性越好，越不易产生误匹配。

关键算法与实现细节分析

1. 频域滤波器构造：利用频率坐标矩阵 [u, v] 计算高斯核，精确建立了光学截止频率与空间域散斑直径的数学对应关系。
2. 变形场映射方程：变形公式 X_deform = (1 + strainX) * ((X - xc) * cos(theta) - (Y - yc) * sin(theta)) + xc + tx 完整涵盖了线性应变、刚体旋转与刚体平移的耦合，且以图像中心为旋转锚点。
3. 动态对比度控制：通过对CGRF结果进行线性拉伸并应用对比度系数（contrastVal），确保了生成的散斑图在直方图分布上占满有效位深，提高了信噪比。
4. 归一化计算：在主动散斑生成中，采用了峰值归一化处理（speckleActive / max），并辅以饱和截断（>1），有效解决了多斑点叠加导致的灰度溢出问题。