[发明专利]一种用于数字图像错误隐藏的数字水印嵌入和提取方法

权利要求书

1.一种用于数字图像错误隐藏的数字水印嵌入方法，其特征在于具体包括以下步骤：

①-1、在数字图像发送端，假设待传输的原始数字图像为8bit灰度图像，记为F，将F中坐标位置为(i,j)的像素点的像素值记为f(i,j)，其中，0≤f(i,j)≤255，1≤i≤I,1≤j≤J，I表示F的宽度，J表示F的高度，不失一般性，在此假设I=J；

①-2、对F进行预处理，得到预处理后的数字图像，记为F'，将F'中坐标位置为(i,j)的像素点的像素值记为f'(i,j)，其中，0≤f'(i,j)≤255；然后对F'进行归一化处理，得到归一化处理后的数字图像，记为F''，将F''中坐标位置为(i,j)的像素点的像素值记为f''(i,j)，f''(i,j)=f'(i,j)/255，其中，0≤f''(i,j)≤1；

①-3、对F''进行L级二维离散小波分解，得到每级对应的一个第一小波逼近子图和三个方向的第一小波细节子图，将对F''进行第l级二维离散小波分解得到的第一小波逼近子图记为FA_l，将对F''进行第l级二维离散小波分解得到的第r个方向的第一小波细节子图记为FD_l(r)，其中，L≥2，1≤l≤L，r=1,2,3，FA_l的分辨率为(I/2^l)×(J/2^l)，FD_l(r)的分辨率为(I/2^l)×(J/2^l)；

①-4、将第L级对应的第一小波逼近子图FA_L中的所有小波域系数的系数值、第L-1级对应的第一小波逼近子图FA_L-1中的所有小波域系数的系数值及第L-1级对应的第r'个方向的第一小波细节子图FD_L-1(r')中的所有小波域系数的系数值均放大q倍，将放大后的第L级对应的第一小波逼近子图记为FA'_L，将放大后的第L-1级对应的第一小波逼近子图记为FA'_L-1，将放大后的第L-1级对应的第r'个方向的第一小波细节子图记为FD'_L-1(r')，其中，FA_L和FA'_L的分辨率均为(I/2^L)×(J/2^L)，FA_L-1、FA'_L-1和FD_L-1(r')、FD'_L-1(r')的分辨率均为(I/2^L-1)×(J/2^L-1)，r'的值为1或2或3，q的取值为FA'_L-1中的最大系数值不超过255时的放大倍数；

①-5、将FA'_L分解为小数部分第一小波逼近子图FAD'_L和整数部分第一小波逼近子图FAI'_L，将FA'_L-1分解为小数部分第一小波逼近子图FAD'_L-1和整数部分第一小波逼近子图FAI'_L-1，将FD'_L-1(r')分解为小数部分第一小波细节子图FDD'_L-1(r')和整数部分第一小波细节子图FDI'_L-1(r')，其中，FAD'_L和FAI'_L的分辨率均为(I/2^L)×(J/2^L)，FAD'_L-1、FAI'_L-1和FDD'_L-1(r')、FDI'_L-1(r')的分辨率均为(I/2^L-1)×(J/2^L-1)；

①-6、将FAI'_L中的每个小波域系数的系数值转化为8bit二进制值，然后按照逐行扫描方式扫描FAI'_L中的每个小波域系数的8bit二进制值，顺序地排列成一个一维二进制序列，记为FAI'_B，再将FAI'_B作为待嵌入的二值数字水印，其中，FAI'_B的长度为I/2^L×J/2^L×8；

①-7、根据FAI'_B中的每位二进制值，按照逐行扫描方式顺序地对FAD'_L-1中的每个小波域系数的系数值进行二值量化处理，实现二值数字水印的嵌入，将经二值量化处理后得到的第L-1级对应的小数部分第一小波逼近子图记为FAD''_L-1，将FAD''_L-1中坐标位置为(i_L-1,j_L-1)的小波域系数的系数值记为fad''_L-1(i_L-1,j_L-1)，当fai'_b((i_L-1-1)×J/2^L-1+(j_L-1-1))=0时，如果fad'_L-1(i_L-1,j_L-1)≥0，则令fad''_L-1(i_L-1,j_L-1)=0.25，如果fad'_L-1(i_L-1,j_L-1)<0，则令fad''_L-1(i_L-1,j_L-1)=-0.25；当fai'_b((i_L-1-1)×J/2^L-1+(j_L-1-1))=1时，如果fad'_L-1(i_L-1,j_L-1)≥0，则令fad''_L-1(i_L-1,j_L-1)=0.75，如果fad'_L-1(i_L-1,j_L-1)<0，则令fad''_L-1(i_L-1,j_L-1)=-0.75，其中，fai'_b((i_L-1-1)×J/2^L-1+(j_L-1-1))表示FAI'_B中的第(i_L-1-1)×J/2^L-1+(j_L-1-1)位二进制值，fad'_L-1(i_L-1,j_L-1)表示FAD'_L-1中坐标位置为(i_L-1,j_L-1)的小波域系数的系数值，FAD''_L-1的分辨率为(I/2^L-1)×(J/2^L-1)；

同样根据FAI'_B中的每位二进制值，按照逐行扫描方式顺序地对FDD'_L-1(r')中的每个小波域系数的系数值进行二值量化处理，实现二值数字水印的嵌入，将经二值量化处理后得到的第L-1级对应的第r'个方向的小数部分第一小波细节子图记为FDD''_L-1(r')，将FDD''_L-1(r')中坐标位置为(i_L-1,j_L-1)的小波域系数的系数值记为fdd''_L-1(r')(i_L-1,j_L-1)，当fai'_b(I/2^L-1×J/2^L-1+(i_L-1-1)×J/2^L-1+j_L-1)=0时，如果fdd'_L-1(r')(i_L-1,j_L-1)≥0，则令fdd''_L-1(r')(i_L-1,j_L-1)=0.25，如果fdd'_L-1(r')(i_L-1,j_L-1)<0，则令fdd''_L-1(r')(i_L-1,j_L-1)=-0.25；当fai'_b(I/2^L-1×J/2^L-1+(i_L-1-1)×J/2^L-1+j_L-1)=1时，如果fdd'_L-1(r')(i_L-1,j_L-1)≥0，则令fdd''_L-1(r')(i_L-1,j_L-1)=0.75，如果fdd'_L-1(r')(i_L-1,j_L-1)<0，则令fdd''_L-1(r')(i_L-1,j_L-1)=-0.75，其中，fai'_b(I/2^L-1×J/2^L-1+(i_L-1-1)×J/2^L-1+j_L-1)表示FAI'_B中的第I/2^L-1×J/2^L-1+(i_L-1-1)×J/2^L-1+j_L-1位二进制值，fdd'_L-1(r')(i_L-1,j_L-1)表示FDD'_L-1(r')中坐标位置为(i_L-1,j_L-1)的小波域系数的系数值，FDD''_L-1(r')的分辨率为(I/2^L-1)×(J/2^L-1)；

①-8、将FAD''_L-1和FAI'_L-1中的所有小波域系数的系数值均缩小q倍，然后将缩小后的FAD''_L-1和缩小后的FAI'_L-1组合成第L-1级对应的新的第一小波逼近子图，记为FA''_L-1，其中，FA''_L-1的分辨率为(I/2^L-1)×(J/2^L-1)；同样将FDD''_L-1(r')和FDI'_L-1(r')中的所有小波域系数的系数值均缩小q倍，然后将缩小后的FDD''_L-1(r')和缩小后的FDI'_L-1(r')组合成第L-1级对应的第r'个方向的新的第一小波细节子图，记为FD''_L-1(r')，其中，FD''_L-1(r')的分辨率为(I/2^L-1)×(J/2^L-1)；

①-9、当r'=1时，对FA''_L-1、FD''_L-1(r')、第L-1级对应的第2个方向的第一小波细节子图FD_L-1(2)、第L-1级对应的第3个方向的第一小波细节子图FD_L-1(3)，进行L-1级二维离散小波逆变换，得到嵌入有二值数字水印的数字图像，记为F'''，其中，F'''的分辨率为I×J；

当r'=2时，对FA''_L-1、FD''_L-1(r')、第L-1级对应的第1个方向的第一小波细节子图FD_L-1(1)、第L-1级对应的第3个方向的第一小波细节子图FD_L-1(3)，进行L-1级二维离散小波逆变换，得到嵌入有二值数字水印的数字图像，记为F'''，其中，F'''的分辨率为I×J；

当r'=3时，对FA''_L-1、FD''_L-1(r')、第L-1级对应的第1个方向的第一小波细节子图FD_L-1(1)、第L-1级对应的第2个方向的第一小波细节子图FD_L-1(2)，进行L-1级二维离散小波逆变换，得到嵌入有二值数字水印的数字图像，记为F'''，其中，F'''的分辨率为I×J；

①-10、将F'''恢复为8bit灰度图像，记为F''''，然后将F''''分块打包传输。