发明名称 一种基于压缩感知的数字图像水印嵌入和提取方法
摘要 一种基于压缩感知的数字图像水印嵌入和提取方法,属于信息隐藏和图像处理技术领域。对二值数字图像水印信息进行稀疏化处理,得到稀疏化的一维水印信息,构造测量矩阵,将测量矩阵作为密钥;用密钥稀疏化的一维水印信息进行测量,实现维水印信息的压缩和加密;对原始载体图像信息进行离散余弦变换并对变换后得到的离散余弦变换域信息进行分块,将处理后的水印信息嵌入到载体图像中得到含水印图像。本发明采用压缩感知方法对水印进行处理,实现水印压缩和加密的双重作用,嵌入量增加,使不可见性和安全性增强,水印总体性能显著提高;应用SB字E矩阵,存储空间小,计算速度快,硬件简单;压缩感知重建过程引入稀疏度,极大提高重建准确性及效果。
申请公布号 CN103700061A 申请公布日期 2014.04.02
申请号 CN201310743836.6 申请日期 2013.12.30
申请人 东北大学 发明人 郎俊;马春雷;张正光
分类号 G06T1/00(2006.01)I;G06T9/00(2006.01)I 主分类号 G06T1/00(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种基于压缩感知的数字图像水印嵌入方法,其特征在于:包括以下步骤:步骤1:对二值数字图像水印信息进行稀疏化处理,得到稀疏化的一维水印信息;读取二值数字图像水印信息W=(w<sub>s,q</sub>)<sub>S×Q</sub>,二值数字图像水印信息W是大小为S×Q的矩阵,其中,S表示矩阵的行,Q表示矩阵的列,w<sub>s,q</sub>为二值数字图像水印信息W的第s行、第q列位置上的元素,w<sub>s,q</sub>∈{0, 1};对二值数字图像水印信息W进行逐列扫描,将其转换成一维水印信息<img file="FDA0000449698280000012.GIF" wi="455" he="72" />一维水印信息C<sup>w</sup>是一维列向量,共有L行,且有L=S×Q;w表示二值数字图像水印信息,是watermark的缩写;一维水印信息C<sup>w</sup>经过离散余弦变换,得到离散余弦变换域信息即稀疏化的一维水印信息<img file="FDA0000449698280000013.GIF" wi="475" he="77" />稀疏化的一维水印信息D<sup>w</sup>是一维列向量,共有L行;设定稀疏化一维水印信息的稀疏度K,0&lt;K<<1,采样数据量为K*L,保留K值;步骤2:构造测量矩阵,将测量矩阵作为密钥,用密钥对步骤1确定的稀疏化的一维水印信息进行测量,实现对一维水印信息的压缩和加密;将P个大小为b×b的哈德马矩阵块H<sub>Block</sub>按照如下公式组成H:<img file="FDA0000449698280000011.GIF" wi="677" he="369" />式中,P和b均为整数,取值满足:P*b=L;对矩阵H的列向量进行随机排序,然后对行向量进行随机排序,抽取排序后矩阵H中的任意L<sub>cs</sub>行组成测量矩阵<img file="FDA0000449698280000014.GIF" wi="308" he="75" />测量矩阵Φ是大小为L<sub>cs</sub>×L的矩阵,<img file="FDA0000449698280000015.GIF" wi="88" he="69" />为测量矩阵Φ的第l<sub>cs</sub>行、第l列位置上的元素,L<sub>cs</sub>满足条件:L<sub>cs</sub><<L,L<sub>cs</sub>∈[K*L,L);cs表示压缩感知,是Compressivesensing的缩写;将上述测量矩阵Φ作为密钥保存;用密钥Φ的每个行向量分别与稀疏化的一维水印信息做内积,得到一维水印信息在密钥上的投影,即测量值<img file="FDA0000449698280000016.GIF" wi="495" he="78" />测量值Y<sup>w</sup>是一维列向量,共L<sub>cs</sub>行;因为一维水印信息的总数据量为L,测量值的总数据量为L<sub>cs</sub>,L<sub>cs</sub><<L,所以测量过程实现了数据压缩,又因为在未知密钥Φ的情况下,压缩过程不可逆,所以压缩过程也是对数据的加密过程;步骤3:读取待进行版权保护的数字图像,将其作为原始载体图像信息,保留一份原始载体图像信息的备份,对原始载体图像信息进行离散余弦变换并对变换后得到的离散余弦变换域信息进行分块,将步骤2处理后的水印信息嵌入到载体图像中得到含水印图像;读取待进行版权保护的数字图像信息,即原始载体图像信息A=(a<sub>m,n</sub>)<sub>M×N</sub>,原始载体图像信息A是大小为M×N的矩阵,其中,M表示行,N表示列,a<sub>m,n</sub>为原始载体图像A的第m行、第n列位置上的元素;对原始载体图像信息A逐列进行扫描,将其转换成一维的原始载体图像<img file="FDA00004496982800000211.GIF" wi="469" he="77" />C<sup>A</sup>是一维列向量,共有L<sub>A</sub>行,L<sub>A</sub>=M×N;对一维的原始载体图像C<sup>A</sup>做离散余弦变换,得到原始载体图像的离散余弦变换域信息<img file="FDA00004496982800000212.GIF" wi="487" he="81" />原始载体图像的离散余弦变换域信息D<sup>A</sup>是一维列向量,共有L<sub>A</sub>行;将原始载体图像的离散余弦变换域信息D<sup>A</sup>的低频部分分割成大小为B×1的矩阵块,<img file="FDA0000449698280000021.GIF" wi="312" he="161" />式中<img file="FDA0000449698280000022.GIF" wi="74" he="83" />为向下取整符号,B是矩阵块的行数,β为伸缩系数,且有β∈(0,0.5],用来调节分块的大小,原始载体图像的离散余弦变换域信息D<sup>A</sup>总共被分割成<img file="FDA0000449698280000023.GIF" wi="131" he="145" />块,<img file="FDA0000449698280000024.GIF" wi="76" he="83" />为向上取整符号;然后将水印信息测量值<img file="FDA00004496982800000213.GIF" wi="470" he="82" />嵌入到原始载体图像的离散余弦变换域信息D<sup>A</sup>前L<sub>cs</sub>个分块的每个分块第一行数据中,此过程可描述为:<img file="FDA0000449698280000025.GIF" wi="1092" he="194" />式中,α∈(0,1)是调节嵌入强度的权值,α越大嵌入强度越大,水印鲁棒性越强,嵌入水印后载体图像质量下降也越大,反之,α越小,水印鲁棒性越弱,嵌入水印后载体图像质量越好;嵌入水印后得到含水印的载体图像离散余弦变换域信息<img file="FDA00004496982800000214.GIF" wi="580" he="95" />为一维列向量,行数为L<sub>A</sub>;含水印的载体图像离散余弦变换域信息<img file="FDA0000449698280000028.GIF" wi="64" he="70" />做离散余弦逆变换,得到一维含水印的图像信息<img file="FDA0000449698280000029.GIF" wi="470" he="92" />一维含水印的图像信息<img file="FDA00004496982800000210.GIF" wi="74" he="73" />是一维列向量,共有L<sub>A</sub>行;将一维含水印的图像信息<img file="FDA0000449698280000031.GIF" wi="66" he="72" />按照逐列扫描的逆过程复原成大小为M×N的含水印图像<img file="FDA00004496982800000318.GIF" wi="422" he="88" />为含水印图像<img file="FDA0000449698280000034.GIF" wi="44" he="71" />的第m行、第n列位置上的元素;水印的嵌入过程结束,步骤1、步骤2和步骤3需要分别保留稀疏化一维水印信息的稀疏度K、密钥Φ、伸缩系数β、权值α以及原始载体图像的备份。
地址 110819 辽宁省沈阳市和平区文化路3号巷11号