MATLAB图像正则化修复工具箱

资 源 简 介

该软件包是一套专门为科学研究和工程应用设计的数字图像修复平台，其核心聚焦于通过正则化技术解决图像处理中的病态逆问题。工具箱完整实现了包括Tikhonov正则化、全变分（Total Variation, TV）正则化、L1范数稀疏正则化在内的多种主流算法模型。其主要功能涵盖了图像去噪、图像去模糊、图像补全以及超分辨率重建。在实现过程中，该软件包采用了高效的数学优化策略，如交替方向乘子法（ADMM）、分裂Bregman方法和带约束的最小二乘优化。用户可以通过该软件包在保持图像边缘锐度与去除随机噪声之间取得平衡

详 情 说 明

MATLAB图像处理正则化工具箱

项目介绍

本工具箱是一个专门用于数字图像修复的科研与工程实验平台，核心致力于解决图像处理中的病态逆问题（Ill-posed Inverse Problems）。通过集成多种经典的正则化理论与现代数学优化算法，该工具箱能够有效处理图像在获取与传输过程中产生的模糊和噪声。软件包通过模块化的设计，实现了从退化模型构造、正则化求解到多维度定量评价的完整实验流程，特别强调在滤除噪声的同时保留图像的关键边缘特征。

---

主要功能特性

• 多模型集成：实现了Tikhonov正则化（L2范数）、全变分（Total Variation, TV）正则化以及基于L1范数的稀疏重构模型。
• 高效优化算法：内置了二阶闭式频率域解法、交替方向乘子法（ADMM）以及分裂Bregman方法（Split Bregman），确证了复杂非平滑目标函数的高速收敛。
• 全流程仿真：涵盖了从原始图像归一化、高斯模糊核构造、加性高斯白噪声注入到最终图像重构的全套仿真环境。
• 多维评价体系：提供峰值信噪比（PSNR）、结构相似性（SSIM）及均方误差（MSE）三种定量评估指标。
• 动态可视化：实时生成算法收敛曲线对比图以及修复效果对比阵列，便于直观评估算法性能。

---

系统要求

• 运行环境：MATLAB R2016b 或更高版本。
• 依赖工具箱：主要算法基于底层数学运算实现，对工具箱依赖度低，推荐安装 Image Processing Toolbox 以获得最佳图像读取与显示体验。

---

主程序实现逻辑说明

工具箱的主执行流程严格遵循以下科学实验步骤：

1. 实验数据准备：

加载测试图像（如Cameraman）并执行浮点归一化。通过构造 15x15 的高斯点扩散函数（PSF）对图像进行模糊处理，并添加标准差为 0.02 的高斯随机噪声，模拟真实的成像退化过程 $y = Kx + n$。

1. Tikhonov 正则化实现：

利用频率域闭式解法。通过拉普拉斯算子构造正则化项，在傅里叶领域利用最小二乘原理计算修复图像。该方法执行速度最快，适用于对计算效率要求极高的场景。

1. TV-ADMM 算法实现：

采用交替方向乘子法求解全变分最小化问题。将原问题拆解为 x-子问题（图像平滑更新）、z-子问题（利用软阈值算子进行梯度收缩）以及乘子更新（u-step）。该算法在抑制噪声的同时能显著保持图像边缘的锐度。

1. Split Bregman 迭代实现：

针对L1范数优化问题，采用分裂Bregman迭代策略。通过引入辅助变量和Bregman变量，将复杂的非线性约束优化转化为一系列子问题的迭代求解，具有良好的收敛性和数值稳定性。

1. 性能分析与可视化：

自动计算各个算法产出图像的最佳性能指标，并绘制目标函数迭代收敛曲线，直观展示不同正则化策略在残差下降速度上的差异。

---

核心算法与关键函数分析

1. 频率域图像恢复逻辑 程序充分利用了循环卷积在频率域变为点乘的特性。通过 psf2otf 函数将空域的点扩散函数转换为光学传递函数，从而在 tikhonov_restoration 中实现了一步到位的反卷积计算，极大降低了计算复杂度。

2. 梯度算子与转置算子 为了支持变分问题的求解，工具箱实现了自定义的梯度计算函数 get_grad（计算水平与垂直梯度）及其转置算子 xt_grad。这两个算子配合使用，能够构造出用于ADMM和Bregman方法中线性方程组的右端项（RHS）。

3. 软阈值收缩算子 (Soft-Thresholding) 在TV-ADMM和Split Bregman的实现中，核心非线性步骤是梯度向量的收缩操作。程序通过计算梯度的模长并应用收缩公式，实现了对图像小跳变（噪声）的抑制和对大跳变（边缘）的保留。

4. 结构相似性 (SSIM) 的底层实现 工具箱内置了简化的SSIM计算逻辑，综合考虑了图像的亮度均值、对比度方差以及结构协方差，相比单一的PSNR指标，更能反映人类视觉系统的感知质量。

---

使用方法

1. 将软件包所有文件置于MATLAB当前工作路径下。
2. 在命令行窗口输入主程序命令并回车。
3. 程序将自动依次执行三种正则化方案，并在控制台输出包含PSNR、SSIM和运行时间的性能评估总结表。
4. 系统将弹出两个图形窗口：一个用于展示不同算法修复后的视觉效果对比，另一个用于展示优化算法的收敛历程。