基于EMD修改方向的高效图像隐写算法仿真实现

资 源 简 介

本算法实现了一种高效的数字图像隐写方法，其核心理念是利用n个载体像素的修改方向来承载一个(2n+1)进制的秘密数字。 该项目的基本原理是将载体图像划分为长度为n的非重叠像素组，并定义一个特定的提取函数，通过对像素灰度值进行加权求和并取模(2n+1)运算，建立像素组到秘密信息的映射关系。 在嵌入过程中，算法首先将二进制秘密信息转换为(2n+1)进制序列。通过比较当前像素组的提取函数计算结果与待嵌入的秘密数字之间的差值，确定对该组内哪一个像素进行+1或-1的操作，或者保持像素值不变。

详 情 说 明

基于修改方向（EMD）的高效图像隐写算法仿真实现

1. 项目介绍

本算法实现了一种经典的数字图像隐写方法——基于修改方向（Exploiting Modification Direction, EMD）的高效隐写算法。该算法的核心理念是利用 $n$ 个载体像素构成的像素组，通过至多修改其中一个像素的值（加1或减1），来承载一个 $(2n+1)$ 进制的秘密数字。这种方法在保证较高嵌入效率的同时，极大地缩减了对图像像素的改动幅度，从而实现了极高的视觉不可感知性。

2. 功能特性

• 高效嵌入能力：每 $n$ 个像素最多只需修改一个像素的值，即可嵌入一个 $(2n+1)$ 进制数据。
• 边界溢出控制：自动识别并处理 0（极黑）和 255（极白）像素，通过预处理将像素值限制在 1-254 范围内，彻底防止加减操作导致的灰度溢出，确保信息提取的 100% 准确率。
• 灵活的进制转换：内置二进制流到多进制（2n+1 进制）的转换逻辑，支持通过调整 $n$ 值平衡容量与隐蔽性。
• 多维评估指标：全自动计算并输出均方误差（MSE）、峰值信噪比（PSNR）、每像素比特数（bpp）以及误码率（BER）。
• 结果可视化：直观展示原始图像、载密图像以及像素差异图，并自动弹出验证成功提示框。

3. 使用方法

1. 启动 MATLAB 软件。
2. 将包含主功能的脚本文件放置在 MATLAB 当前工作路径下。
3. 在命令行窗口输入该脚本定义的函数名并按回车。
4. 程序将自动读取图像、嵌入随机秘密信息、执行提取过程并生成详尽的性能报告及可视化图像。

4. 系统要求

• 软件环境：MATLAB R2016a 或更高版本。
• 依赖工具箱：Image Processing Toolbox（图像处理工具箱）。
• 硬件建议：标准的桌面或笔记本电脑环境即可流畅运行。

5. 核心逻辑与实现过程

载体准备与预处理 程序首先尝试加载标准图像（如 cameraman.tif），若读取失败则自动生成灰度渐变图作为补充。 核心预处理步骤是将图像转换为双精度浮点型，并对处于灰度边缘的值进行强制调整：

• 所有 0 值修改为 1。
• 所有 255 值修改为 254。

这一步骤通过极小的代价保证了隐写后的像素值不会超出 [0, 255] 的合法区间。

信息编码转换 秘密信息生成为二进制比特流。根据 $n$ 的取值计算出 $(2n+1)$ 进制基数。 为了简化转换并避免大型矩阵运算导致的精度失效，程序采用分组转换法：将二进制流按照固定的位长度（由进制基数的对数向下取整决定）切割，并转换为对应的十进制值，作为待嵌入的符号序列。

EMD 嵌入逻辑 图像被划分为不重叠的像素组。对于每一个组：

1. 计算提取函数 $f(g_1, g_2, dots, g_n) = (sum_{i=1}^{n} g_i times i) mod (2n+1)$。
2. 计算待嵌入数字 $d$ 与当前函数值 $f$ 的差值 $s = (d - f) mod (2n+1)$。
3. 执行如下修改准则：

* 若 $s = 0$，则不修改任何像素。 * 若 $1 le s le n$，则将该组内第 $s$ 个像素的值加 1。 * 若 $s > n$，则将该组内第 $(2n+1-s)$ 个像素的值减 1。

信息提取逻辑 在提取阶段，接收者只需将载密图像同样划分为长度为 $n$ 的像素组，通过相同的加权求和并取模公式，即可直接计算出嵌入的 $(2n+1)$ 进制数字，无需原始载体参考。

6. 算法关键细节分析

提取函数的设计 程序中实现的提取函数 $f$ 利用了系数向量 $[1, 2, dots, n]$。这种线性加权取模的方法确保了对组内任何一个像素进行 $pm 1$ 的修改，都能在结果 $f$ 上产生唯一的偏移，从而覆盖 $[0, 2n]$ 的所有可能状态。

修改方向的唯一性 通过对 $s$ 取值的判断，算法能够精准定位需要修改的像素下标。这种映射机制保证了对于每一个像素组，最多只需要变动一个像素的一级灰度值，这种极简的修改方案是 EMD 算法高 PSNR 的基础。

7. 性能评估体系

程序在执行完毕后会自动输出以下技术指标：

• 容量（Capacity/bpp）：衡量单位像素承载的信息位数，受 $n$ 值影响。
• 均方误差（MSE）：计算原始图像与载密图像之间的平方误差均值。
• 峰值信噪比（PSNR）：衡量图像失真度的关键指标。由于 EMD 算法每次仅修改 $pm 1$ 灰度级，该值通常远高于 50dB。
• 误码率（BER）：通过比对原始二进制流与提取出的比特流，验证算法的鲁棒性与正确性。