发明名称 基于QR分解的数字图像水印嵌入方法和数字图像认证方法
摘要 本发明的基于QR分解的数字图像水印嵌入方法和数字图像认证方法相互配合进行认证,所述的认证方法不需要原始的图像,只需要将数据S,T1,T2的取值与所述嵌入方法中的相同即可,且在图像的QR分解中具有一条重要的性质,即R矩阵第一行的元素的绝对值大于其它各行对应的元素,因此在R矩阵的第一行元素中嵌入水印信息,对图像的视觉质量不会产生大的影响,经试验证明R矩阵这种性质的普遍存在性以及该认证方法的有效性、可靠性,既不会引起图像视觉质量的下降,又能够对图像的恶意篡改做到精确认证和定位;本发明技术方案不需要额外的数据管理;认证快捷迅速,适合数字图像这种大数据量的计算处理;可以判断篡改类型;可以精确地进行篡改定位。
申请公布号 CN102156954B 申请公布日期 2012.08.22
申请号 CN201110079913.3 申请日期 2011.03.31
申请人 河南省臻嘉科技有限公司 发明人 汪萍;刘粉林;巩道福
分类号 G06T1/00(2006.01)I 主分类号 G06T1/00(2006.01)I
代理机构 郑州联科专利事务所(普通合伙) 41104 代理人 刘建芳;马柯柯
主权项 1.一种基于QR分解的数字图像水印嵌入方法,其特征在于:具体步骤如下:(1)、对大小为M×N的图像I分块,分解成大小均为n×n的图像块阵列,表示为:<img file="FSB00000804960500011.GIF" wi="1078" he="448" />M、N分别为图像I的行像素数、列像素数,n为分解得到的各图像块的行、列像素数,设B<sub>i,j</sub>为分解后得到的图像块阵列中位于第i行、第j列的图像块,1≤i≤t,1≤j≤r,<img file="FSB00000804960500012.GIF" wi="179" he="109" /><img file="FSB00000804960500013.GIF" wi="153" he="109" />2≤n≤min(M,N),min(M,N)为取M和N中最小的一个,M,N分别是n的整数倍,M,N若不是n的整数倍则添0补足;(2)、对每一个图像块B<sub>i,j</sub>(i=1,2,…,t,j=1,2,…,r)执行如下步骤:1A)、对图像块B<sub>i,j</sub>进行QR分解,表示为B<sub>i,j</sub>=Q<sub>i,j</sub>×R<sub>i,j</sub>,其中Q<sub>i,j</sub>、R<sub>i,j</sub>分别为分解得到的大小为n×n的Q矩阵和R矩阵,QR分解的过程为:Q<sub>i,j</sub>是一个具有标准正交向量的n×n矩阵,Q<sub>i,j</sub>的列由B<sub>i,j</sub>中的列通过格拉姆一施密特正交化处理得到,R<sub>i,j</sub>是一个n×n的上三角矩阵,设B<sub>i,j</sub>和Q<sub>i,j</sub>分别为<img file="FSB00000804960500014.GIF" wi="399" he="99" /><img file="FSB00000804960500015.GIF" wi="424" he="73" />其中<img file="FSB00000804960500016.GIF" wi="145" he="64" />分别为B<sub>i,j</sub>和Q<sub>i,j</sub>的列向量,那么矩阵R<sub>i,j</sub>可由下式计算得到:<img file="FSB00000804960500017.GIF" wi="1523" he="380" />其中&lt;·,·&gt;表示向量的内积;1B)、由<img file="FSB00000804960500021.GIF" wi="400" he="66" />生成认证水印W<sub>i,j</sub>={w<sub>1</sub>,w<sub>2</sub>,…,w<sub>n</sub>},其中<img file="FSB00000804960500022.GIF" wi="374" he="63" />(l=1,2,…,n),其中C(·)为水印生成函数,具体的水印生成过程为:令<img file="FSB00000804960500023.GIF" wi="394" he="82" />首先对<img file="FSB00000804960500024.GIF" wi="36" he="51" />使用MD5方法,得到128位的二进制数据<img file="FSB00000804960500025.GIF" wi="511" he="54" />然后将<img file="FSB00000804960500026.GIF" wi="41" he="54" />按位进行异或,得到1位的水印信息w<sub>l</sub>,即<maths num="0001"><![CDATA[<math><mrow><msub><mi>w</mi><mi>l</mi></msub><mo>=</mo><msubsup><mover><mi>w</mi><mo>^</mo></mover><mi>l</mi><mn>1</mn></msubsup><mo>&CirclePlus;</mo><msubsup><mover><mi>w</mi><mo>^</mo></mover><mi>l</mi><mn>2</mn></msubsup><mo>&CirclePlus;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&CirclePlus;</mo><msubsup><mover><mi>w</mi><mo>^</mo></mover><mi>l</mi><mi>k</mi></msubsup><mo>&CirclePlus;</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>&CirclePlus;</mo><msubsup><mover><mi>w</mi><mo>^</mo></mover><mi>l</mi><mn>128</mn></msubsup><mo>,</mo></mrow></math>]]></maths>其中<img file="FSB00000804960500028.GIF" wi="37" he="39" />为异或运算,<img file="FSB00000804960500029.GIF" wi="196" he="63" />为<img file="FSB000008049605000210.GIF" wi="41" he="54" />的第k位数据,k=1,2,…,128;1C)、将水印W<sub>i,j</sub>通过量化嵌入的水印嵌入方法嵌入在R<sub>i,j</sub>矩阵的第一行的n个元素r<sub>1,1</sub>,r<sub>1,2</sub>,…,r<sub>1,n</sub>中,表示为:{r′<sub>1,1</sub>,r′<sub>1,2</sub>,…,r′<sub>1,n</sub>}=E({r<sub>1,1</sub>,r<sub>1,2</sub>,…,r<sub>1,n</sub>},{w<sub>1</sub>,w<sub>2</sub>,…,w<sub>n</sub>}),其中W<sub>i,j</sub>={w<sub>1</sub>,w<sub>2</sub>,…,w<sub>n</sub>}为生成的水印,{r<sub>1,1</sub>,r<sub>1,2</sub>,…,r<sub>1,n</sub>}为嵌入水印前的原始的R<sub>i,j</sub>矩阵中的第一行的元素,{r′<sub>1,1</sub>,r′<sub>1,2</sub>,…,r′<sub>1,n</sub>}为嵌入水印后的R<sub>i,j</sub>矩阵的第一行的元素,其中:<img file="FSB000008049605000211.GIF" wi="801" he="156" />(l=1,2,…,n)其中mod为模运算,S为用户设定的量化因子,S的值越大,则水印抵抗噪声的能力越强,但对图像视觉质量的影响也越大,S的值越小,则水印抵抗噪声的能力越弱,但对图像视觉质量的影响也越小,<img file="FSB000008049605000212.GIF" wi="229" he="109" /><img file="FSB000008049605000213.GIF" wi="205" he="110" />为阈值;如此,得到含有水印的R矩阵R′<sub>i,j</sub>:<img file="FSB000008049605000214.GIF" wi="656" he="300" />1D)、进行逆QR分解变换,得到含水印的图像块B′<sub>i,j</sub>,表示为B′<sub>i,j</sub>=Q<sub>i,j</sub>×R′<sub>i,j</sub>,其中R′<sub>i,j</sub>为含水印的R矩阵;(3)、将各个含水印的图像块按照步骤(1)中分解时的顺序组合成含水印的图像I′,组合方法为:<img file="FSB00000804960500031.GIF" wi="1366" he="325" />
地址 450001 河南省郑州市文化路85号A时代广场2110室