一种基于哈尔变换的四进制彩色数字图像盲水印方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110415154A(43)申请公布日2019.11.05(21)申请号201910714537.7(22)申请日2019.08.04(71)申请人鲁东大学地址264025山东省烟台市芝罘区红旗中路186号(72)发明人苏庆堂刘得成袁子涵徐效美王环英张雪婷(51)Int.Cl.G06T1/00(2006.01)权利要求书3页说明书8页附图2页(54)发明名称一种基于哈尔变换的四进制彩色数字图像盲水印方法(57)摘要本发明结合变换域数字水印算法鲁棒性强的优点，公开了一种基于哈尔变换的四进制彩色数字图像盲水印方法。本发明根据哈尔变换具有计算复杂度低的特点，对图像块进行哈尔变换得到变换后矩阵中的最大能量系数，并将最大能量系数量化到多个区间来完成彩色数字水印的嵌入与盲提取。本发明能将彩色图像数字水印嵌入到彩色宿主图像中，不但具有较强的鲁棒性、较好的不可见性和较高安全性的特点，而且具有大容量的特点，适用于高安全性的大容量彩色数字图像版权保护的场合。CN110415154ACN110415154A权利要求书1/3页1.一种基于哈尔变换的四进制彩色数字图像盲水印方法，包含彩色图像数字水印的嵌入过程和提取过程，其水印嵌入的具体过程描述如下：Step1：彩色宿主图像的预处理：将大小为M×M的24位彩色宿主图像H作分层处理，得到红、绿、蓝三个分层宿主图像Hp，并将各层中的像素划分为m×m的非重叠图像块，p=1,2,3分别表示红、绿、蓝三层；Step2：彩色水印图像的预处理：对大小为N×N的24位彩色水印图像W作分层处理，得到红、绿、蓝三个分层水印图像，随之对每层进行基于私钥Kp的仿射变换得到置乱后的分层水印图像Wp，p=1,2,3分别表示红、绿、蓝三层；然后，将各分层水印图像Wp中的每个十进制像素值转换为一个4位的四进制序列，并依次拼接为水印字符串watp，p=1,2,3分别表示红、绿、蓝三层；Step3：构建m阶的哈尔变换矩阵harm，其中，m为宿主图像块的分块大小；Step4：选择图像块：利用选块矩阵selection从分层宿主图像Hp中选择图像块haarblock，其中，选块矩阵selection由MATLAB内置函数randperm伪随机生成，p=1,2,3分别表示红、绿、蓝三层；Step5：求得图像块的最大能量系数：根据公式（1），利用哈尔变换矩阵harm对图像块haarblock进行哈尔变换，得到变换后的矩阵HR，其第一行第一列元素HR(1,1)即为图像块的最大能量系数Emax；（1）-1其中，m为宿主图像块的分块大小，harm为m阶的哈尔变换矩阵，(harm)为矩阵harm的逆矩阵；Step6：计算最大能量系数的上下边界：按照先后顺序从水印字符串中选取水印位w，利用公式（2）、（3），求得最大能量系数的上边界Eupper和下边界Elower；（2）（3）其中，Emax是图像块haarblock的最大能量系数，mod(.)是求余函数，T是量化步长，为量化系数，，；*Step7：计算最优边界值：利用公式（4）求得最大能量系数的最优边界值Emax，并用其替*换原最大能量系数Emax，得到新的矩阵HR；2CN110415154A权利要求书2/3页（4）其中，Eupper和Elower分别为最大能量系数的上下边界，abs(.)是取绝对值函数；Step8：求得含水印图像块：利用公式（5）进行逆哈尔变换，得到含水印图像块haarblock*；（5）-1其中，m为宿主图像块分块的大小，harm为m阶的哈尔变换矩阵，(harm)为矩阵harm的逆矩阵；Step9：获得含水印宿主图像：重复执行本过程的Step4-Step8，直到所有的水印位**都被嵌入为止；最后重组三层含水印的分层宿主图像Hp得到含水印彩色宿主图像H；其水印提取的具体过程描述如下：Step1：彩色含水印宿主图像的预处理：将含水印彩色宿主图像H*作分层处理，得到红、*绿、蓝三层含水印的分层宿主图像Hp，并将各分层宿主图像中的像素划分为m×m的非重叠图像块，p=1,2,3分别表示红、绿、蓝三层；Step2：构建m阶的哈尔变换矩阵harm，其中，m为宿主图像的分块大小；*Step3：选取含水印图像块：利用选块矩阵selection从分层宿主图像Hp中选择含水印图像块haarblock*，其中，选块矩阵selection由MATLAB内置函数randperm伪随机生成，p=1,2,3分别表示红、绿、蓝三层；Step4：求得含水印块的最大能量系数：根据公式（6），利用哈尔变换矩阵harm对含水印图像块haarblock*进行哈尔变换，得到含水印的矩阵HR*，其第一行第一列元素HR*(1,1)即**为图像块haarblo