[0042] 本发明的目的是提供一种融合离散小波变换的空域彩色数字图像盲水印方法，其特征在于通过具体的水印嵌入过程和提取过程来实现的，其水印嵌入过程描述如下：
[0043] 第一步：彩色图像数字水印的预处理：首先，将一幅大小为32×32的彩色图像数字水印W依照红、绿、蓝三基色的顺序分成3个分层水印图像Wi；然后，将每个分层水印图像进行基于Kai的Arnold置乱变换；最后，将每个十进制数表示的像素用8位二进制数表示(例2
如：可将215转换成二进制数11010111)，并依次连接形成长度为8×32＝8192(bits)的分层水印位序列SWi，其中i＝1,2,3,分别表示红、绿、蓝三层；
[0044] 第二步：获取宿主图像的嵌入块：将一幅大小为512×512的原始彩色宿主图像H分成3个分层宿主图像Hi；同时，把每一个分层宿主图像Hi分成大小为2×2的图像块；根据分层水印信息长度8192bits，利用基于密钥Kbi的伪随机置乱算法在分层宿主图像Hi中选择合适位置的图像块，其中8192<＝(512×512)/(2×2),其中i＝1,2,3,分别表示红、绿、蓝三层；
[0045] 第三步：选取一个图像块A，不需要真正的二维离散小波变换，而该变换的低频近似分量中第一个系数值CA11就可以在空域中按照公式(1)直接求出；
[0046]
[0047] 其中，ap，q表示图像块A中第p行第q列元素，1≤p，q≤2；
[0048] 此处，设选取的图像块A为 则CA11＝(211+210+213+209)/2＝421.5000；
[0049] 第四步：从分层水印序列SWi中按先后顺序取出一个嵌入水印信息w，依据该嵌入水印信息及公式(2)、(3)，在空域中计算当将该水印信息量化嵌入到CA11上时，将得到两个极限值CAlow与CAhigh；
[0050] 如果w＝′0′，则
[0051] 如果w＝′1′，则
[0052] 其中，floor(.)是向下取整函数，T是正整数量化步长；此时，设w＝‘0’，T＝40，则根据公式(2)、(3)，得CAlow＝410，CAhigh＝450；
[0053] 第五步：依据公式(4)，将水印w嵌入到CA11，并得到其含水印的结果[0054]
[0055] 其中，abs(.)是绝对值函数， 表示低频近似分量中第一个系数在嵌入水印后的值；此时，CA11＝421.5000，CAlow＝410，CAhigh＝450，所以根据公式(4)[0056] 第六步：依据公式(5)，计算在图像块A的所有元素中嵌入水印后的结果 并用*其替换原有相应元素，得到真正含水印的图像块A ；这样在空域中通过修改该图像块完成了在离散小波变换后修改其低频分量中第一个系数值CA11来嵌入水印信息w的过程；
[0057]
[0058] 其中，ap，q表示嵌入水印前图像块A中第p行第q列的值， 表示嵌入水印后图像块A中第p行第q列的值，1≤p，q≤2；此时，原始图像块 CA11＝421.5000，根据公式(5)在空域中得到真正含水印的图像块
[0059] 第七步：重复执行上述第三步到第六步，直到所有的水印信息都被嵌入完成为止；*
最后，将含水印的红、绿、蓝分层图像重新组合并获得含水印图像H；
[0060] 其水印提取过程描述如下：
[0061] 第一步：将含水印图像H*分成3个分层含水印图像 其中i＝1，2，3，分别表示红、绿、蓝三层，并将每一分层含水印图像 进一步分成2×2的非重叠图像块；
[0062] 第二步：在分层含水印图像 中，利用上述水印嵌入过程中所提到的基于密钥Kbi的伪随机置乱算法选择图像块；
[0063] 第三步：选取一个图像块A*，在空域中利用公式(6)计算其离散小波变换后其低频近似分量中第一个系数值
[0064]
[0065] 其中， 表示含水印图像块A*中第p行第q列的值，1≤p，q≤m，m是图像块的尺寸*大小；此时m＝2，设选取的图像块A 为 则利用公式(6)得其
[0066] 第四步：利用公式(7)，提取图像块A*中所含有的水印w*；
[0067]
[0068] 其中，mod(.)是求余函数，T为正整数量化步长；此时， T＝40，利*用公式(7)提取所含有的水印w＝‘0’；
[0069] 第五步：重复执行第三步、第四步，提取每层二进制水印序列SWi*，然后每8位二进制信息为一组转换成十进制的像素值，其中i＝1，2，3，分别表示红、绿、蓝三层；
[0070] 第六步：对转化后的每层十进制像素执行基于密钥Kai的逆Arnold变换并获得分*层水印Wi，其中i＝1，2，3，分别表示红、绿、蓝三层；
[0071] 第七步：组合分层水印Wi*形成最终的提取水印W*。
[0072] 该方法既具有频域数字水印算法的优点，也具有空域数字水印算法的优点，鲁棒性强，实时件高，水印不可见性好；算法简单快捷，适用于彩色图像作为数字水印的高效版权保护。
[0073] 本发明有效性验证
[0074] 为了证明本发明的有效性，选择如图1(a)、图1(b)所示的两幅大小为512×512的24位标准图像作为宿主图像，并分别用如图2(a)、图2(b)所示的两幅大小为32×32的24位彩色图像作为数字水印进行验证。
[0075] 图3(a)、图3(b)是将图2(a)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像，其结构相似度SSIM值依次是0.97380、0.97125，其峰值信噪比PSNR值依次是41.9139dB、41.9532dB；图4(a)、图4(b)是依次从图3(a)、图3(b)中提取的水印，其归一化互相关系数NC值分别是1.00000、1.00000；图5(a)、图5(b)、图5(c)、图5(d)、图5(e)、图5(f)是将图3(a)所示的含水印图像依次进行JPEG2000压缩(5：1)、椒盐噪声(2％)、中值滤波(3×1)、锐化(1.0)、缩放(400％)、剪切(25％)等攻击后所提取的水印，其归一化互相关系数NC值分别是0.99695、0.99759、0.98473、0.99892、0.99971、0.97580。
[0076] 图6(a)、图6(b)是将图2(b)所示的水印依次嵌入到宿主图像图1(a)、图1(b)后所得到的含水印图像，其结构相似度SSIM值依次是0.97388、0.97227，其峰值信噪比PSNR值依次是41.9038dB、42.0572dB；图7(a)、图7(b)是依次从图6(a)、图6(b)中提取的水印，其归一化互相关系数NC值分别是1.00000、1.00000；图8(a)、图8(b)、图8(c)、图8(d)、图8(e)、图8(f)是将图6(a)所示的含水印图像依次进行JPEG2000压缩(5：1)、椒盐噪声(2％)、中值滤波(3×1)、锐化(1.0)、缩放(400％)、剪切(25％)等攻击后所提取的水印，其归一化互相关系数NC值分别是0.99716、0.99433、0.97838、0.99933、0.99956、0.98880。
[0077] 该算法在平台2.27GHZ CPU,2.00GB RAM,Win 7,MATLAB 7.10.0(R2010a)上进行过近万次运行，其数字水印的平均嵌入时间是0.374421秒，平均提取时间是0.338252秒，总计时间为0.712673秒。
[0078] 综上所述，所嵌入的彩色图像数字水印具有较好的不可见性，满足了水印算法的不可见性要求；同时，从各种受攻击图像中所提取的彩色图像数字水印具有较好的可鉴别性和较高的NC值，说明该方法具有较强的鲁棒性；另外，该算法的平均运行总时间小于1秒，满足了多媒体大数据快速版权保护的需要。
