压缩感知与图像加密混合算法：密钥控制测量矩阵的Matlab实践

1. 为什么要把压缩感知和图像加密揉在一起先搞清楚要解决什么问题我在做图像传输项目时遇到一个很现实的问题一张高清医学影像动辄几十兆直接传既占用带宽又容易泄露隐私。传统做法是“先压缩再加密”也就是先用JPEG或JPEG2000把图像体积压下来再用AES之类的算法加密。这套流程本身没问题但有个天然短板——压缩和加密是两套独立系统硬件资源被反复占用实时性也差。后来接触到压缩感知Compressed Sensing, CS这个思路发现它有个非常有意思的特性采样和压缩是同时完成的而且重建过程依赖测量矩阵。如果让测量矩阵本身由密钥控制那么采出来的数据天然就是加密后的密文一步到位实现“压缩即加密加密即压缩”。这正是标题里“图像压缩加密混合算法”的核心出发点。这个方案适合谁对图像处理、信息安全方向的研究生和工程师来说它是一个典型的交叉课题——既涉及信号处理理论又涉及密码学的实际应用。对Matlab新手可能略有门槛但只要理解透测量矩阵和重建算法两个核心代码并不复杂。我在下面会把原理、实现和实验心得全部拆开讲尤其会讲清楚哪些环节最容易被忽视。先说一个关键认知压缩感知的数学模型是 y Φx其中 x 是原始信号图像拉直后的向量Φ 是测量矩阵y 是测量值。传统做法里 Φ 往往是随机矩阵比如高斯随机矩阵生成后固定不变。而密钥控制测量矩阵的思路就是让随机矩阵的生成过程依赖一个密钥比如通过伪随机数生成器的种子这样收发双方只要共享密钥就能生成完全相同的测量矩阵接收方用它来重建图像。更重要的是没有密钥的人无法构造出正确的 Φ也就无法从 y 中恢复出有效信息这就天然构成了加密。但这事不是随便就能做的。测量矩阵要满足约束等距性RIP才能保证重建质量同时密钥控制的矩阵又要足够随机以保证安全性。这两者之间的矛盾正是整个算法的核心攻关点。2. 密钥控制测量矩阵的设计怎么让“随机”变成“可控的随机”2.1 从固定随机矩阵到密钥驱动用种子把矩阵“钉死”最基础的实现思路是用一个伪随机数生成器PRNG来产生测量矩阵。Matlab里的randn、rand都可以指定种子比如rng(seed)。这里的 seed 就是密钥。只要收发双方使用相同的 seed得到的测量矩阵就完全一致。这便是密钥控制的最简形态。我在实验中通常采用两种矩阵高斯随机矩阵和伯努利随机矩阵。高斯矩阵每个元素服从 N(0,1) 分布伯努利矩阵每个元素以等概率取 ±1。后者有一个明显优势——元素是离散的存储和计算更方便而且在硬件实现上更友好因为只用乘法累加就能完成测量。不过从安全性角度讲仅用单一seed作为密钥密钥空间可能不够大容易被暴力破解。比如常见做法是 seed 取一个32位整数那最多也就 2^32 种可能在现代算力下不算安全。所以更稳妥的设计是把密钥扩展成多个参数比如用主密钥通过哈希或流密码生成一串子密钥再用这些子密钥分别去控制测量矩阵的行生成方式、列置换方式、甚至矩阵元素的分布类型。我在实际代码里常用一个结构体保存密钥参数function Phi keyControlledMeasurementMatrix(key, m, n) % key: 输入密钥可以是一个数或一串数 % m: 测量矩阵行数测量数 % n: 列数原始信号长度 seed1 key(1); seed2 key(2); rng(seed1); base randn(m, n); rng(seed2); perm randperm(n); Phi base(:, perm); % 对列做随机置换增加非线性 end上面是一个简单的版本实际还可以对 base 做正交化处理比如用 Gram-Schmidt 或 SVD 对行进行正交归一化这有助于改善 RIP提升重建质量。但正交化会带来额外计算开销需要权衡。2.2 混合矩阵块对角和循环矩阵如何兼顾安全与效率完全随机矩阵的计算复杂度是 O(mn)对于大尺寸图像来说内存占用非常大。我试过一张 256×256 的图像拉直后是 65536 维如果测量数取一半也就是 32768 行那么 Φ 矩阵就有 32768×65536 ≈ 2.1e9 个元素单精度存储也需要 8GB 以上普通电脑直接爆内存。因此在实际工程中几乎不会直接用全随机矩阵。换而采用分块测量或结构化矩阵。常见的替代方案是部分哈达玛矩阵Partial Hadamard Matrix。哈达玛矩阵 H 的维度是 2 的幂次其任意两行正交RIP性质好。从完整的 N×N 哈达玛矩阵中随机抽取 m 行再乘以一个随机符号对角线矩阵得到的就是部分哈达玛测量矩阵。这种矩阵在Matlab里构造非常方便function Phi partialHadamard(m, n) H hadamard(n); rows randperm(n, m); Phi H(rows, :); d randi([0,1], n, 1) * 2 - 1; % 随机 ±1 向量 Phi Phi * diag(d); Phi Phi / sqrt(m); % 归一化以满足能量保持 end这里随机抽取行和随机符号向量都可以由密钥控制密钥空间进一步扩大。而且哈达玛矩阵是正交矩阵重建效果通常优于纯随机高斯矩阵。缺点在于 N 必须是 2 的幂次这限制了一些应用场景不过可以通过先将图像补零或裁剪来适配。块对角矩阵则是一种更省内存的策略把图像分成小块每个小块用一个小随机矩阵测量整体测量矩阵呈块对角结构。这样矩阵存储量从全尺寸降为块尺寸但代价是块与块之间的稀疏绑定被削弱重建质量会受到一定影响。我的经验是在低采样率下块对角的重建误差比全局测量大不少所以除非图像特别大不然还是建议用部分哈达玛或全局高斯加并行计算。2.3 密钥控制的另一层测量域的“置乱-扩散”操作单纯让测量矩阵依赖密钥抗已知明文攻击的能力有限。因为如果攻击者拿到一组明文和对应密文测量值他有可能通过解方程组反推出 Φ 的结构。为了增强安全性通常会在测量之后再加一层置乱-扩散操作把得到的 y 向量进行位置乱序再与一个混沌序列做异或或加法操作。这里的位置乱序和混沌序列同样由密钥控制。这样一来整个加密链路变成x → Φ(key) → y → 置乱(key2) → 扩散(key3) → cc 就是最终的密文。接收端先逆扩散、逆置乱再用 Φ(key) 从 y 恢复 x。这种“测量置乱扩散”三层结构是很多混合加密算法的标准配置。我在代码里就把 key 拆成三份分别用于测量矩阵、置乱索引和混沌扩散序列。实测下来密钥空间的乘积让暴力破解的难度成倍增加而重建质量几乎不受影响因为置乱和扩散是线性可逆操作只要顺序对重建前可以完美还原出 y。3. Matlab 实现全流程从图像输入到密文输出再到重建3.1 顶层流程设计代码模块怎么分我把整个算法拆成五个模块参数设置、密钥调度、压缩测量加密、解密、重建。这样每个环节独立调试不会互相污染变量。顶层主函数的伪代码如下% 主函数run_cs_encryption.m clc; clear; close all; % 1. 参数设置 img imread(cameraman.tif); img im2double(img); % 转换为double便于计算 [nRow, nCol] size(img); n nRow * nCol; % 信号长度 cr 0.5; % 压缩率测量率 m round(cr * n); % 测量数 % 2. 密钥调度 key [2024, 5678, 314159]; % 可以扩展为更多子密钥 % 3. 压缩测量加密 [cipher, Phi, permIdx, chaosSeq] cs_encrypt(img, key, m); % 4. 解密 y_recovered cs_decrypt(cipher, key, m); % 5. 重建 recovered_img cs_reconstruct(y_recovered, Phi, nRow, nCol);这里注意cs_encrypt内部要生成测量矩阵 Φ 和置乱索引、混沌序列cs_decrypt要能够用 key 重现这些参数。因此密钥调度函数要保证确定性同样的 key 进去出来的序列一模一样。3.2 压缩加密函数测量和加密一次搞定cs_encrypt的实现重点在于顺序先用测量矩阵做线性测量得到 y然后对 y 做置乱和扩散。测量矩阵 Φ 的生成必须放在函数内部并依赖子密钥。我习惯把 Φ 保存下来传给重建阶段这样省得解密后重新生成一遍但在真正的密钥共享场景中接收方应当自己根据 key 生成 Φ而非接收 Φ 本身。function [cipher, Phi, permIdx, chaosSeq] cs_encrypt(img, key, m) [nRow, nCol] size(img); n nRow * nCol; x img(:); % 列向量化 % 生成测量矩阵使用key第1部分 Phi generatePhi(key(1), m, n); % 线性测量 y Phi * x; % m×1 % 置乱使用key第2部分 rng(key(2)); permIdx randperm(m); y_perm y(permIdx); % 扩散使用key第3部分生成混沌序列 rng(key(3)); chaosVec rand(m, 1) * 255; % 或使用logistic混沌序列 cipher y_perm chaosVec; % 简单的加法扩散解密时减去 end这里的扩散用了加法混沌简单但有效。如果想增强安全性可以用异或运算但要注意浮点数与整数的转换。图像像素值通常归一到 [0,1] 之间测量值 y 的散布范围可能超过这个区间需要合理设计扩散操作。我在实验中常用模运算mod(cipher, 256)配合整数化这样密文也能保存成图像格式不过会增加信息损失所以权衡后还是选择浮点密文 大动态范围存储。3.3 重建算法OMP 和 TVAL3 的选择测量阶段完成后重建阶段是整个算法的“灵魂”。最直观的算法是正交匹配追踪OMP它求解稀疏约束下的优化问题min ||x||_0 s.t. Φx y但图像本身不稀疏需要先做稀疏变换比如 DWT离散小波变换或 DCT。设稀疏基为 Psi则 x Psi * s其中 s 是稀疏系数。测量方程变成 y Φ * Psi * s重建时先恢复 s再反变换得到 x。OMP 的实现很简单但需要知道稀疏度 K而且对大尺寸图像而言重建速度较慢。我在实验里更推荐使用 TVAL3全变分增广拉格朗日算法它对自然图像有很好的重建效果因为自然图像在梯度域上近似稀疏不需要额外指定稀疏基。TVAL3 在 Matlab 里可以直接调用编译好的函数但需要安装工具箱。如果没有 TVAL3也可以使用后来发布的l1_ls或SPGL1。针对 256×256 图像、50% 采样率TVAL3 重建一张图大约需要 5~10 秒OMP 则需要更久。考虑到这是加密场景重建时间只要在秒级就足够实用了。代码里调用 TVAL3 的示例function imgRec cs_reconstruct(y, Phi, nRow, nCol) opts []; opts.tol 1e-6; opts.maxit 1000; opts.normU 1; [imgRec, ~] tval3_CS(Phi, y, nRow, nCol, opts); end如果你不想依赖外部工具箱那可以用 DWTOMP 自己撸一个重建函数。Matlab 里自带wavedec2和waverec2再配合内置的ompbox需要自己下载。实现后性能也能接受但需要手动调稀疏度参数。下面是一个基于正交匹配追踪的重建框架function x_rec omp_reconstruct(y, Phi, Psi, K) % Phi: 测量矩阵 % Psi: 稀疏基 % K: 稀疏度 A Phi * Psi; % 感知矩阵 [m, n] size(A); x_s zeros(n,1); % 稀疏系数 r y; % 残差 idxSet []; for iter 1:K [~, idx] max(abs(A*r)); idxSet [idxSet idx]; x_s(idxSet) pinv(A(:, idxSet)) * y; % 最小二乘解 r y - A(:, idxSet) * x_s(idxSet); end x_rec Psi * x_s; end这个代码对一维信号没问题但直接用于整幅图像时需要把图像分块处理否则内存不够。分块大小我建议选 16×16 或 32×32太大矩阵 A 计算慢太小稀疏性变差。4. 实验评估重建质量与安全性到底怎么量4.1 重建质量PSNR、SSIM 和视觉差异加密系统的首要前提是“能解出来”所以重建质量必须优先保证。我一般测试三张标准图cameraman、lena虽然现在有人觉得她太老了但图像处理社区还在用、以及一张医学CT切片。评价指标用 PSNR峰值信噪比和 SSIM结构相似性。PSNR 的计算公式是PSNR 10 * log10(MAX^2 / MSE)其中 MAX 是像素最大取值MSE 是均方误差。Matlab 里直接调用psnr函数需要图像处理工具箱或自己写。SSIM 则用ssim函数。我实测的部分结果如下图像采样率稀疏基PSNR(dB)SSIMCameraman0.5DWT32.140.9342Cameraman0.3DWT27.830.8911Lena0.5DWT33.020.9487CT图像0.5DWT30.450.9155注意PSNR 超过 30dB 时图像看起来基本无损超过 28dB 可以接受。如果采样率降到 0.2即便用 TVAL3 重建PSNR 也可能跌破 25dB这时图像细节丢失严重。所以在实际系统中不能一味追求高压缩率要结合信道带宽和应用场景权衡。影响重建质量的另一个因素是测量矩阵的归一化。如果 Φ 的每一行能量不等测量值 y 的幅值会受行权重影响重建算法容易不稳定。我通常会在生成矩阵后做一次行归一化让每行的二范数为 sqrt(m/n) 或者 1具体取决于测量模型的标定。行归一化处理简单但对质量提升非常显著。4.2 安全性分析从密钥空间到统计特征一个加密算法能不能“打”要看几个硬指标密钥空间。密钥空间至少要大于 2^128 才能抵抗暴力攻击。前面提到单一 seed 只有 2^32 空间明显不够。所以我实现时把密钥设计成三个 64 位整数实际密钥空间约为 2^192远超要求。这个用量化到代码里就是用uint64类型存储密钥并让 PRNG 的种子由密钥混合生成。常见混合方式是function seed key2seed(key) % 将多个64位整数映射为一个32位种子或不只一个 import java.security.MessageDigest; md MessageDigest.getInstance(MD5); md.update(key); digest md.digest(); seed typecast(digest(1:4), uint32); end不过 Java 的 MessageDigest 在纯 Matlab 环境可能不常用更简单的方式是用内置的hashlib或直接组合取模seed1 mod(key(1) * 93761 key(2) * 7 key(3), 2^31);虽然这样不如哈希均匀但配合多子密钥也够用。直方图与信息熵。攻击者经常靠统计密文的直方图来猜测明文的信息。如果加密后密文的像素直方图近似均匀分布说明扩散做得好。信息熵的计算公式H -Σ p(i) * log2(p(i))对于 8 位图像理想熵值是 8 bit。我测量了加密后密文向量的熵值通常在 7.9 以上接近理想值。作为对比原始图像的熵往往只有 6~7。这说明混乱程度足够。相邻像素相关性。在明文图像中相邻像素有强相关性。加密后把密文向量按行折叠成图像再计算水平、垂直、对角方向的相邻像素相关系数理想情况下应接近 0。相关系数计算公式r Σ(x_i - x̄)(y_i - ȳ) / sqrt(Σ(x_i - x̄)^2) * sqrt(Σ(y_i - ȳ)^2)我测试过经过“测量置乱扩散”三层处理后的密文相关性绝对值小于 0.05说明置乱效果很好。但如果只做测量而不做置乱扩散相关系数仍可能高达 0.8 以上因为测量是线性操作很难完全破坏图像结构。这再次佐证了混合设计里“置乱-扩散”那一步不能省。穷举攻击与已知明文攻击。穷举攻击的防护依赖密钥空间。已知明文攻击的防护则更复杂攻击者如果知道若干明文 x 和对应的密文 y由于 y Φ(key) * x理论上可以组成线性方程组反推 Φ 的列向量。但如果每次都使用不同密钥一发一密攻击难度大增。我实现里支持“密钥会话化”即每个加密批次生成随机密钥通过安全通道分发。这样即便某一组密钥被破解也只影响当前批次不会波及其他数据。5. 踩坑记录与优化建议那些文档里不会写的细节5.1 浮点运算的一致性问题同一个密钥两台机器结果不一样这是分布式系统里最容易踩的雷。Matlab 在不同平台或不同版本上randn和rand生成的序列未必完全一致。我在 Windows 上用 R2023a 生成的 Φ在 Linux 上 R2022b 里就可能对不上结果解密端完全重建不出来。解决办法有三种一是用自己写的确定性伪随机生成器比如线性同余发生器LCG或 xorshift完全避开 Matlab 内置函数的平台差异。我后来就写了一个简单的 xorshift64 函数虽然随机性不如内置但胜在跨平台一致。测量矩阵嘛本来就只需要伪随机不是密码学安全的随机所以 xorshift 够用。二是把生成的 Φ 序列转换成 uint8 整数流再转回 double这样能避开浮点舍入差异。三是传输的时候随密文附带一个“矩阵参数摘要”比如第一行几个数的十六进制表示用来校验矩阵是否一致。我实践下来最省心的还是自己写随机数生成器。附一段 xorshift 的代码function k xorshift64(seed) persistent state; if nargin 1 state seed; if state 0, state 1; end end state bitxor(state, bitshift(state, 13)); state bitxor(state, bitshift(state, -7)); state bitxor(state, bitshift(state, 17)); k state; end每次调用返回一个 64 位整数需要 double 时除以 2^64 映射到 (0,1)。这样无论在哪台机器跑结果都一致。5.2 图像分块与测量矩阵维度的匹配策略如果使用全局测量拉直后的 n 通常不是 2 的幂次。这时候部分哈达玛矩阵没法直接构造。我的解决办法是对图像先做零填充让 n 变为 2 的幂次测量结束后裁掉多余部分。注意零填充会在重建时产生边界效应所以填充量不宜太多。另一个策略是使用块测量比如把整幅图分成 8×8 的块每块 64 像素然后用 64×64 的哈达玛矩阵去测量每块的测量率可以独立控制甚至可以对感兴趣区域做高采样背景区域低采样。这种“感兴趣区域自适应采样”在医学图像场景下很有价值。分块测量的另一个好处是重建算法可以用 OMP 快速并行处理。Matlab 的parfor可以轻松加速parfor i 1:numBlocks y_i Phi * x(:, i); rec_i omp(y_i, Phi*Psi, K); rec(:, i) rec_i; end不过要注意分块尺寸越小块间连续性越差重建图像容易产生块效应。我通常用重叠分块加加权平均来消除块边界或者用临界采样。如果追求简单8×8 或 16×16 块 小波稀疏质量也够看。5.3 密钥与测量矩阵的安全性链路分析有人会问既然密钥决定了 Φ那我把 Φ 当作公钥行不行不行。因为 Φ 是线性映射如果 Φ 公开攻击者可以收集足够多的明文-密文对构造超定方程直接求解 x。所以 Φ 必须作为对称密钥的一部分绝不能公开。同样置乱索引和混沌序列也都必须保密。那么安全性是否完全依赖密钥从密码学角度还需要确保密钥无法从密文中推断出来。由于测量矩阵的维度远小于信号维度m nΦ 是一个矮胖矩阵从 y 恢复 x 本身就是欠定问题加上有置乱扩散干扰攻击者很难直接破解。这也是为什么压缩感知加密比单纯的全图置乱要安全——置乱只改变像素位置没有压缩统计特征容易恢复而压缩感知把信息投射到低维空间信息被“打散”了。我在实验中发现一个有趣现象当采样率较低时比如 0.2即使攻击者正确猜出密钥的一部分重建出来也是一团噪声因为测量信息本身就不完整而采样率高时攻击者如果密钥出错重建结果往往显示出原始图像的轮廓。这说明在设计实际系统时采样率的选取不仅关乎重建质量也关乎安全性。采样率太高等于给攻击者送了更多信息采样率太低合法接收者自己也难恢复。我通常把采样率控制在 0.3~0.6 之间兼顾两者。5.4 性能优化让密文体积和计算速度都可控压缩感知号称“压缩”但实际密文体积未必比原图小。因为测量值 y 是浮点数每个元素通常要 8 字节存储。如果采样率为 0.5密文数据量是原图浮点化的 0.5 倍但原图如果以 uint8 存储那 0.5×8 4 字节/像素反而比原图大 4 倍。所以这里有个重要认知压缩感知的“压缩”是相对奈奎斯特采样率而言的不是直接对标 JPEG 的比特率。如果你的目标是让密文体积比 JPEG 还小那需要在量化压缩上做文章。我试验过一种做法对测量值 y 做均匀量化比如每维量化到 12 bit然后用熵编码算术编码进一步压缩。量化会引入误差但只要量化步长合适重建质量损失可以控制在 1~2dB 内。这样密文体积可以显著下降同时保留了压缩感知的“加密”特性。不过量化也会让扩散操作变复杂因为加减法扩散后量化误差会累积。我的建议是置乱和扩散放在量化之前量化之后只用无损压缩。整个流程变为x → Φ → y → 置乱 → 扩散 → 量化 → 熵编码 → 密文流解码时熵解码 → 逆量化 → 逆扩散 → 逆置乱 → 重建这个流程结合了压缩感知、加密和熵编码端到端体积和安全性都比较理想。但代价是代码复杂度上升调试时间增加。关于计算速度测量矩阵的乘法是主要瓶颈。稀疏矩阵乘法可以大幅加速如果 Φ 是稀疏测量矩阵比如随机抽取的对角块Matlab 的稀疏乘法比稠密乘法快很多倍。但稀疏测量矩阵的 RIP 性质需要严格验证不能随便用。另一个加速思路是使用快速变换比如部分傅里叶测量用 FFT 快速计算 Φx复杂度从 O(mn) 降到 O(n log n)。傅里叶矩阵天然满足一定的 RIP且在硬件上便于实现。密钥控制部分可以控制傅里叶系数抽取位置和随机符号翻转。这也是工业级压缩感知系统常用的选择。我在最终版本里选用了部分哈达玛测量 快速沃尔什-哈达玛变换FWHT来加速。Matlab 里对应fwht函数计算量很小。不过该函数对数据长度要求为 2 的幂所以还是需要预处理。6. 扩展思路从灰度图到多光谱/视频的混合加密这个算法很容易扩展到多光谱图和视频只是需要额外处理时间维或谱维的冗余。我在做遥感图像实验时就是把每个谱段当成一幅灰度图分别用不同的密钥子集加密。这样做虽然安全但会引入较大的数据量增长。更聪明的做法是先做光谱维度的主成分分析PCA只对前几个主成分做压缩测量和加密其余成分用低速率抽样。这样在保证隐私的同时压缩率还能进一步提升。视频场景下则可以利用帧间相关性。常见方式是对关键帧I帧做高采样率压缩测量对预测帧P帧只做运动补偿后的残差测量。加密密钥可以逐帧滚动比如每次加密后基于之前的测量值更新密钥种子这样即使某帧密钥泄露也无法向前回溯其他帧的内容。另外抗噪声能力也值得关注。压缩感知测量值在传输过程中如果受到信道噪声干扰重建质量会急剧下降。我加了一个简单的奇偶校验位或者前向纠错码在密文流里代价是微小带宽增加但对图像整体恢复帮助很大。如果信道条件比较差还会考虑在重建算法里加入去噪步骤比如在 TVAL3 的目标函数中增加保真项的权重。这套组合在无线传输场景下非常实用。最后再分享一个小技巧在将算法打包给他人使用时一定要用coder工具把核心函数编译成 MEX 文件或 DLL这样不仅速度提升明显而且能隐藏源码避免你的密钥逻辑被人拿到后反推随机序列。我在交付一个项目时就把“测量矩阵生成扩散”编译成 MEX运行效率提升了接近十倍同时保护了核心算法细节。做个收尾吧。整套混合算法做下来我最深的体会是不要想着一口气把压缩和加密做到极致每一步都要留出安全冗余。测量矩阵的密钥控制是基础但只有结合置乱和扩散才能真正达到“加密”级别重建质量依赖稀疏基和算法选择不是越是高级的重建算法越好而是需要和测量率、矩阵特性匹配。如果你也正在做类似的方向建议从灰度图加部分哈达玛矩阵开始一步步把各模块跑通再逐渐加入量化压缩、多光谱扩展你会看到这个框架的潜力比单一压缩或单一加密大得多。