发明名称 面向数字作品分级权益管理的免疫图像数字水印方法
摘要 本发明属于信息技术领域,涉及一种面向数字作品分级权益管理的免疫图像数字水印方法,包括水印嵌入及分级质量锁控处理、分级权益管理与保护两大部分,其中,水印嵌入及分级质量锁控处理由水印嵌入和分级质量锁控过程构成,水印嵌入包括水印预处理、版权水印嵌入和权益水印嵌入及其质量锁控,水印嵌入后再对各层级图像质量进行分级锁控;分级权益管理与保护由密钥分配与权益分配构成。本发明能够实现分级使用权益和版权权益保护的同时,进行相关权益的有效分离,使其彼此没有重叠,有效阻止依据相邻像素相似性的攻击,并进而使其成为授权方式下的一种权益管理方式,具有高度安全性,为数字社区内安全共享数字图像作品提供了强有力的技术支持。
申请公布号 CN106408503A 申请公布日期 2017.02.15
申请号 CN201611069097.7 申请日期 2016.11.29
申请人 海南大学 发明人 姚孝明;王豪
分类号 G06T1/00(2006.01)I 主分类号 G06T1/00(2006.01)I
代理机构 海口翔翔专利事务有限公司 46001 代理人 莫臻
主权项 一种面向数字作品分级权益管理的免疫图像数字水印方法,其特征在于,包括水印嵌入及分级质量锁控处理、分级权益管理与保护两大部分,其具体过程如下:第一部分:水印嵌入及分级质量锁控处理该部分由水印嵌入和分级质量锁控两个关键过程构成,一、水印嵌入水印嵌入过程包括水印预处理、版权水印嵌入和权益水印嵌入及其质量锁控;1)、水印预处理:水印包括版权水印和权益水印两种,分别标记为W<sub>b</sub>、W<sub>1</sub>,二者均为含有版权特征或权益特征且像素值取0与1两种数值的二值图像;采用现有常规方法对其进行置乱或加密处理,并采用相应置乱方法将上述水印信息转化为密文,相应的参数作为密钥保存以备后用,将其二值形式按以下公式由0,1变换为‑1,1,即:对于水印矩阵第i行第j列元素W(i,j);<img file="FDA0001164818770000011.GIF" wi="598" he="174" />2)、版权水印嵌入:给定原始图像I,版权水印W<sub>b</sub>,嵌入深度alpha,质量分级参数为n,具体水印嵌入过程如下:(ⅰ)n级小波分解:根据现有小波分解理论,选择某种小波基做n级小波分解,其中第n级低尺度近似分量、低尺度近似分量分别与水平方向、垂直方向和对角方向的高频分量分别记为A<sub>n</sub>,H<sub>n</sub>,V<sub>n</sub>及D<sub>n</sub>;原始图像I转化为尺度不等的系数矩阵集合{A<sub>n</sub>,H<sub>n</sub>,V<sub>n</sub>,Dn,H<sub>n‑1</sub>,V<sub>n‑1</sub>,D<sub>n‑1</sub>,…,H<sub>i</sub>,V<sub>i</sub>,D<sub>i</sub>,…,H<sub>1</sub>,V<sub>1</sub>,D<sub>1</sub>},由此可得其中A<sub>n</sub>;(ⅱ)版权水印嵌入:对A<sub>n</sub>采用现有非盲水印的方法按以下公式嵌入版权水印W<sub>b</sub>,得含水印分量矩阵<img file="dest_path_image001.GIF" wi="36" he="26" />即:<img file="dest_path_image002.GIF" wi="192" he="26" />其中,alpha是版权水印嵌入深度;(ⅲ)n级小波重构:根据现有小波重构理论,由<img file="dest_path_image003.GIF" wi="37" he="26" />H<sub>n</sub>,V<sub>n</sub>,D<sub>n</sub>重构得<img file="dest_path_image004.GIF" wi="49" he="26" />由<img file="dest_path_image005.GIF" wi="49" he="26" />H<sub>n‑1</sub>,V<sub>n‑1</sub>,D<sub>n‑1</sub>可以重构<img file="dest_path_image006.GIF" wi="49" he="26" />…;由A<sub>1</sub>’,H<sub>1</sub>,V<sub>1</sub>,D<sub>1</sub>可以重构得到含有版权水印的图像<img file="dest_path_image007.GIF" wi="22" he="22" />3)、权益水印嵌入及其质量锁控:给定含版权水印图像<img file="dest_path_image008.GIF" wi="22" he="22" />权益水印W<sub>1</sub>,嵌入深度α<sub>1</sub>(取值范围为(0,1]),密钥种子key1,具体嵌入及其质量锁控过程如下:(ⅰ)伪随机矩阵的准备:依据密钥种子key1采用既有伪随机序列生成算法得到行列大小与A<sub>n‑1</sub>相等的伪随机矩阵β<sub>1</sub>;分别计算权益水印W<sub>1</sub>的行列奇偶校验和Ch与Cl,并将其尾部填充0使其长度大小与key1相同,将三者逐位做XOR运算的结果作为新密钥种子得到另一行列大小与A<sub>n‑1</sub>相等的伪随机矩阵random;(ⅱ)权益水印嵌入:将A<sub>n‑1</sub>视为IA区域,按以下公式嵌入其权益水印W<sub>1</sub>得嵌入水印后的<img file="dest_path_image009.GIF" wi="49" he="22" /><img file="FDA00011648187700000210.GIF" wi="1430" he="175" />(ⅲ)图像质量锁控:对于权益水印所在的子图A<sub>n‑1</sub>,其质量锁控方法有些 特殊,分以下三个步骤处理:步骤1:计算补偿矩阵:SRC<sub>n‑1</sub>(i,j)=A<sub>n‑1</sub>(i,j)‑A'<sub>n‑1</sub>(i,j)其中,A<sub>n‑1</sub>与A'<sub>n‑1</sub>分别是嵌入权益水印前后的第n‑1级小波分解低尺度分量子图矩阵;步骤2:将权益水印W<sub>1</sub>嵌入到SRC<sub>n‑1</sub>内,得到:SRC'<sub>n‑1</sub>(i,j)=SRC<sub>n‑1</sub>(i,j)‑β<sub>1</sub>(i,j)×W<sub>1</sub>(i,j)其中,W<sub>1</sub>为权益水印;β<sub>1</sub>为上述嵌入权益水印时采用的伪随机矩阵;步骤3:依据加密密钥K(ID)对上述矩阵加密,得到:SSRC<sub>n‑1</sub>=EN(SRC'<sub>n‑1</sub>,K(ID))其中:EN()为现有典型对称加密算法,,密钥K(ID)由系统根据用户注册ID作为种子自动生成;(ⅳ)数据标准化:选用min‑max标准化方法对上述得到的A'<sub>n‑1</sub>和SSRC<sub>n‑1</sub>做标准化处理,具体也分以下三个步骤来完成:步骤一:计算:<img file="dest_path_image010.GIF" wi="351" he="26" /><img file="dest_path_image011.GIF" wi="415" he="26" /><img file="dest_path_image012.GIF" wi="294" he="26" />其中:min()为取最小值函数,max()为取最大值函数,P<sub>1</sub>为结构化系统数据变量;步骤二:计算标准幅值范围:根据理论分析及实验,按下式计算:对于小波第k级子图分量,可取A<sub>mxk</sub>=255×2<sup>k</sup>,A<sub>mnk</sub>=S<sub>mxk</sub>=40×2<sup>k</sup>,Smnk=‑40×2<sup>k</sup>。步骤三:根据下列标准化变换式,可得:<img file="FDA0001164818770000041.GIF" wi="1621" he="310" />(ⅴ)小波重构复原:将上述计算得到的标准化结果SSRC<sup>*</sup><sub>n‑1</sub>作为<img file="dest_path_image013.GIF" wi="41" he="22" />替换原D<sub>n‑1</sub>,再由A<sup>*</sup><sub>n‑1</sub>,H<sub>n‑1</sub>,V<sub>n‑1</sub>,<img file="dest_path_image014.GIF" wi="41" he="26" />一起通过小波变换重构得到<img file="dest_path_image015.GIF" wi="49" he="26" />将标准化过程得到的结构化系统变量P1最后作为权益密钥之一记录到数据库中保存;二、分级质量锁控水印嵌入后,需要对各层级图像质量进行分级锁控,具体按如下四个步骤处理:步骤1:计算补偿矩阵:将D<sub>n‑k</sub>分量矩阵视为UIA区域,则第k级补偿矩阵数据SRC<sub>n‑k</sub>由以下公式获得:SRC<sub>n‑k</sub>(i,j)=A<sub>n‑k</sub>(i,j)‑A'<sub>n‑k</sub>(i,j)其中:A<sub>n‑k</sub>为未嵌入权益水印的原始第n‑k级小波分解低尺度子图分量,<img file="dest_path_image016.GIF" wi="49" he="26" />为由第n‑k+1级小波分解嵌入水印后各个分量即<img file="dest_path_image017.GIF" wi="64" he="26" />H<sub>n‑k+1</sub>,V<sub>n‑k+1</sub>,<img file="dest_path_image018.GIF" wi="60" he="26" />的重构值;步骤2:加密补偿矩阵:依据加密密钥K(ID)对上述矩阵加密,得到:SSRC<sub>n‑k</sub>=EN(SRC,K(ID))其中:EN()为现有典型对称加密算法;密钥K(ID)由系统根据用户注册ID作为种子自动生成;步骤3:标准化处理:(ⅰ)计算:<img file="dest_path_image019.GIF" wi="355" he="26" /><img file="dest_path_image020.GIF" wi="404" he="26" /><img file="dest_path_image021.GIF" wi="291" he="26" />其中:min()为取最小值函数,max()为取最大值函数,Pk为结构化系统数据变量;(ⅱ)对于小波第k级子图分量,可取A<sub>mxk</sub>=255×2<sup>k</sup>,A<sub>mnk</sub>=S<sub>mxk</sub>=40×2<sup>k</sup>,S<sub>mnk</sub>=‑40×2<sup>k</sup>;(ⅲ)对<img file="dest_path_image022.GIF" wi="45" he="26" />和<img file="dest_path_image023.GIF" wi="79" he="26" />按下式进行数据标准化处理,可得:<img file="FDA0001164818770000056.GIF" wi="1622" he="310" />步骤4:小波重构:将SSRC<sup>*</sup><sub>n‑k</sub>作为<img file="dest_path_image024.GIF" wi="45" he="26" />替换原D<sub>n‑k</sub>,再由A<sup>*</sup><sub>n‑k</sub>,H<sub>n‑k</sub>,V<sub>n‑k</sub>,<img file="dest_path_image025.GIF" wi="49" he="26" />一起通过小波变换重构<img file="dest_path_image026.GIF" wi="64" he="26" />重复以上步骤1~4,直到所有各级权益质量分级锁控处理完毕;如此,最终由A*<sub>1</sub>,H<sub>1</sub>,V<sub>1</sub>,D*<sub>1</sub>可以重构得到公开免疫图像<img file="dest_path_image027.GIF" wi="22" he="22" />将标准化过程得到的结构化系统变量P<sub>k</sub>(k=1,2,…,n‑1)最后作为权益密钥之一记录到数据库中保存。;第二部分:分级权益管理与保护分级权益管理与保护由以下密钥分配与权益分配两个环节构成;一、密钥分配根据不同权益需求,将相应密钥分为相对独立的四类:A、一类权益密钥:{SRC加解密密钥K(ID)、幅值范围P1、水印置乱密钥}B、二类权益密钥:{SRC加解密密钥K(ID)、幅值范围P1、伪随机矩阵β<sub>1</sub>种子密钥key1}C、三类权益密钥:{SRC加解密密钥K(ID)、幅值范围Pk}D、四类权益密钥:{版权水印密钥}二、权益分配依据上述四类密钥可以合理分配相应的权益,具体实现方法如下:①针对第一类密钥所有者,有权查看水印信息,由以下五个步骤完成:步骤1:将公开免疫图像<img file="dest_path_image028.GIF" wi="15" he="22" />进行n级小波分解,得第n‑1级子图分量A<sub>n‑1</sub>,H<sub>n‑1</sub>,V<sub>n‑1</sub>,D<sub>n‑1</sub>,从而有如下计算公式:<img file="FDA0001164818770000062.GIF" wi="526" he="183" />步骤2:据密钥P<sub>1</sub>,按下式进行反标准化处理得到SSRC<sub>n‑1</sub>:<img file="FDA0001164818770000063.GIF" wi="1564" he="151" />其中min()为取最小值函数,max()为取最大值函数;<img file="dest_path_image029.GIF" wi="120" he="26" />即为由密钥P<sub>1</sub>中提取的原嵌入水印后SSRC幅值范围;步骤3:根据SRC加解密密钥K(ID)可得补偿矩阵<img file="dest_path_image030.GIF" wi="64" he="22" />SRC'<sub>n‑1</sub>=DE(SSRC<sub>n‑1</sub>,K(ID))其中:DE()为对称密码算法;步骤4:按下式计算权益水印密文数据:<img file="FDA0001164818770000071.GIF" wi="694" he="158" />步骤5:依据水印置乱密钥,做逆变换处理得到权益水印明文;②第二类密钥所有者,有权获得水印密文以及较高质量的近似图像<img file="dest_path_image031.GIF" wi="52" he="26" />具体由以下五个步骤完成:步骤1:重复第一类密钥拥有者的步骤1~4,获得权益水印密文<img file="dest_path_image032.GIF" wi="45" he="26" />步骤2:由密钥P<sub>1</sub>提取嵌入水印后A<sub>n‑1</sub>的幅值范围,依下式进行反标准化处理得到<img file="dest_path_image033.GIF" wi="52" he="26" /><img file="FDA0001164818770000075.GIF" wi="1318" he="150" />其中:min()为取最小值函数,max()为取最大值函数,A<sub>mn1</sub>,A<sub>mx1</sub>即为由密钥P<sub>1</sub>提取的嵌入水印后A<sub>n‑1</sub>的幅值范围;步骤3:由伪随机矩阵β<sub>1</sub>密钥key1,根据以下公式计算α<sub>1</sub>近似值,具体过程如下:1)相邻像素对:在<img file="dest_path_image034.GIF" wi="41" he="26" />中依<img file="dest_path_image035.GIF" wi="37" he="26" />取值相反即分别为1和‑1的条件按照先水平后垂直的扫描次序寻找相邻像素对,其行列坐标分别记作(i,j)和(ii,jj);2)若记abs()为求绝对值的函数,则可按下列公式估算α'值:<img file="FDA0001164818770000078.GIF" wi="758" he="157" /><img file="FDA0001164818770000079.GIF" wi="862" he="127" /><img file="FDA00011648187700000710.GIF" wi="598" he="135" />3)估算α'<sub>1</sub>值:重复步骤1)、2),直到遍历<img file="dest_path_image036.GIF" wi="22" he="22" />区域所有相邻像素对,对所有符合要求的α'求平均值即α<sub>1</sub>的近似值,记作α'<sub>1</sub>;步骤4:根据以下公式得到恢复近似图像A<sub>n‑1</sub>:<img file="FDA0001164818770000082.GIF" wi="934" he="158" />步骤5:将含有补偿矩阵的D<sub>n‑1</sub>区域数据置零,由更新后的A<sub>n‑1</sub>,H<sub>n‑1</sub>,V<sub>n‑1</sub>,D<sub>n‑1</sub>重构A<sub>n‑2</sub>,得到近似恢复图;③针对第三类密钥所有者,有权获得相应层级的高质量恢复图像<img file="dest_path_image037.GIF" wi="64" he="26" />具体步骤如下:(a).根据公开图像小波分解得到的A<sup>*</sup><sub>n‑k</sub>,按照下式计算得到<img file="dest_path_image038.GIF" wi="52" he="26" /><img file="FDA0001164818770000085.GIF" wi="1373" he="151" />其中min()为取最小值函数,max()为取最大值函数;<img file="dest_path_image039.GIF" wi="120" he="26" />即为由密钥P<sub>k</sub>提取的原A<sub>n‑k</sub>嵌入水印后的幅值范围;(b).根据公开图像小波分解得到的SSRC<sup>*</sup><sub>n‑k</sub>,按照下式计算得到SSRC<sub>n‑k</sub>:<img file="FDA0001164818770000087.GIF" wi="1613" he="151" />其中min()为取最小值函数,max()为取最大值函数;<img file="dest_path_image040.GIF" wi="117" he="26" />即为由密钥P<sub>k</sub>中提取的原SSRC幅值范围;(c).根据SRC加解密密钥K(ID)可得补偿矩阵SRC<sub>n‑k</sub>:SRC<sub>n‑k</sub>=DE(SSRC<sub>n‑k</sub>,K(ID))其中:DE()为前述对称密码算法如AES、DES等的解密算法;(d).根据补偿矩阵恢复图像:A<sub>n‑k</sub>(i,j)=SRC<sub>n‑k</sub>(i,j)+A'<sub>n‑k</sub>(i,j)(e).将含有补偿矩阵的D<sub>n‑k</sub>区域数据置零,由更新后的A<sub>n‑k</sub>,H<sub>n‑k</sub>,V<sub>n‑k</sub>,D<sub>n‑k</sub>可以重构得到第k+1级精确恢复图像A<sub>n‑k‑1</sub>;④第四类密钥拥有者,有权检测或查看版权水印,具体步骤如下:1)分别将公开图像I’和原图像I做n级小波分解,得到<img file="dest_path_image041.GIF" wi="34" he="26" />与A<sub>n</sub>;2)依据现有非盲鲁棒水印算法中其版权水印嵌入的公式及其相应的解密算法求出版权水印的估计值<img file="dest_path_image042.GIF" wi="41" he="26" />3)将W<sub>b</sub>和<img file="dest_path_image043.GIF" wi="37" he="26" />按图像方式显示并加以对比,并由此证实版权权益所属;或按现有水印检测方式来确认。
地址 570228 海南省海口市美兰区海甸岛人民大道58号