一种基于指数柯西瑞利建模的数字水印检测方法技术

本发明专利技术公开一种基于指数柯西瑞利建模的数字水印检测方法，首先对载体图像执行UDWT变换，选出目标子带，而后执行PCET，确定嵌入位置，通过构造非线性嵌入函数，将水印以乘性嵌入方式嵌入到载体图像中，再经过矩重构以及UDWT逆变换得到含水印图像；然后通过分析UDWT

【技术实现步骤摘要】
一种基于指数柯西瑞利建模的数字水印检测方法


[0001]本专利技术涉及数字图像版权保护
，具体为一种基于指数柯西瑞利建模的数字水印检测方法。

技术介绍

[0002]近年来互联网的飞速发展，使人们足不出户就可以在网络上搜索到任何想要获取的信息。但是互联网在带来便利的同时，也带来了很多安全隐患，诸如多媒体作品的版权问题等。为此数字水印技术横空出世。数字水印技术是指利用数字隐藏原理将版权信息以水印信号的形式、按不可察觉的方式嵌入到所保护的载体中，这种操作方式的优势是不会损害原作品的质量，即人的肉眼无法察觉嵌入的水印信号，只能依靠特定解码器或检测器来提取识别。水印可以应用于数字图像、视频以及音频等诸多方向。该技术在保护作者的著作权、追踪数字作品的非法使用等方面发挥了重要作用。
[0003]不可感知性、鲁棒性和水印容量是衡量数字水印算法性能优劣的三个重要指标。其中，水印的不可感知性是指嵌入水印后的载体图像与原图像在感知上没有差别；鲁棒性是指水印算法的抗攻击能力，一般来说，水印算法的抵抗攻击能力越强，表明鲁棒性越好；水印容量指的是可以隐藏在数字图像中的水印信息总量，从理论上来说，水印容量越大越好，但实际应用中，过大的水印容量会降低含水印图像的视觉质量。同时提高鲁棒性、不可感知性和水印容量一直是数字水印算法领域面临的挑战。
[0004]近年来，国内外研究团队对数字图像水印做了大量的工作，开发出了一系列极具研究价值的统计模型水印算法。但依旧存在一些不足之处：1)大多数的数字水印算法仅在变换域中嵌入水印，然而变换域系数不稳定、鲁棒性不足，导致水印算法的性能不理想；2)现有的水印嵌入方法多采用线性嵌入，而样本系数通常是非线性的，因此现有的线性嵌入方法并不完美；3)水印方案的性能依赖于对图像系数的有效建模，提升检测精度的关键在于更精确地描述图像信息特征，确立精确的统计模型，现有的基于传统分布函数的模型对图像存在的阴影、边缘模糊等微小细节的描述还不够精准。

技术实现思路

[0005]本部分的目的在于概述本专利技术的实施方式的一些方面以及简要介绍一些较佳实施方式。在本部分以及本申请的说明书摘要和专利技术名称中可能会做些简化或省略以避免使本部分、说明书摘要和专利技术名称的目的模糊，而这种简化或省略不能用于限制本专利技术的范围。
[0006]为解决上述技术问题，根据本专利技术的一个方面，本专利技术提供了如下技术方案：
[0007]一种基于指数柯西瑞利建模的数字水印检测方法，
[0008]约定：I表示载体图像；N
b
表示目标子带分割的总块数；L表示N
b
块中能量排序最高的块数；n表示阶数；m表示重复度；B
l
(l＝1,2,...,L)表示第l个块中UDWT
‑
PCET幅值系数的集合；B
l
'表示含水印图像的第l个块中UDWT
‑
PCET幅值系数的集合；w
l
表示由{
‑
1，1}生成的
伪随机序列；I'表示含水印图像；λ表示尺度参数；β表示形状参数；x
i
表示原始UDWT
‑
PCET幅值系数；y
i
表示含水印图像的UDWT
‑
PCET幅值系数；
[0009]所述基于指数柯西瑞利建模的数字水印检测方法包括水印嵌入和水印提取：
[0010]所述水印嵌入步骤如下：
[0011]S1：初始设置：输入载体图像I并初始化设置；
[0012]S2：水印嵌入：
[0013]S2.1：对载体图像I进行二级非抽样离散小波变换(UDWT)，尺度二将获得三个高频子带，分别计算三个高频子能量，选择能量最大子带作为目标子带，能量计算公式如下：
[0014][0015]其中，X
j
[a,b]为第j个方向子带系数；
[0016]S2.2：将目标子带划分为N
b
个不重叠且等大小的块，计算熵值并降序排列，对N
b
个块中熵值最高的前L个块进行n阶极坐标复指数变换(PCET)，以一块为例，选取[
‑
m,m]&(m≠0)&[1,n]区域为目标区域，计算得到L个目标PCET幅值系数集合B
l
(l＝1
,
2,...,L)用于水印信息的嵌入；
[0017]S2.3：利用非线性嵌入强度函数修改每块UDWT
‑
PCET幅值系数以隐藏水印信息，水印嵌入表达式如下：
[0018][0019]其中f1(x)和f0(x)分别是嵌入水印为1和
‑
1时，对应的嵌入强度函数，f1(x)和f0(x)具体表示如下：
[0020][0021][0022]a1,b1,a2,b2为需要结合仿真实验确定的参数；
[0023]S2.4：极坐标复指数矩重构，再执行逆UDWT，从而获得水印图像I'；
[0024]所述水印提取步骤如下：
[0025]S3：基于指数柯西瑞利分布的UDWT
‑
PCET幅值系数统计建模：
[0026]S3.1：与上述嵌入过程相同，对载体图像进行UDWT二级分解，选取尺度二能量最大子带作为目标子带，将目标子带划分为N
b
个不重叠且等大小的块，并进行n阶PCET，在每块中选取与嵌入位置相同位置上的幅值系数作为参数估计的训练样本；
[0027]S3.2：利用指数柯西瑞利分布对UDWT
‑
PCET幅值系数进行统计建模，假设随机变量x服从指数柯西瑞利分布，则其表达式如下：
[0028][0029]S3.3：通过基于退火遗传算法的ML优化参数估计方法，求解模型中的参数，用于解
码器的构造；
[0030]S4：构造最大似然解码器：
[0031]S4.1：根据ML决策准则推导出基于指数柯西瑞利分布的最大似然解码器，首先，根据二元假设，H1和H0分别表示UDWT
‑
PCET幅值系数嵌入“1”和嵌入
“‑
1”的两种情况：
[0032][0033]S4.2：为提取水印w
l
，根据ML决策准则，取对数后得到的对数似然比为：
[0034][0035]其中，i∈B
′
l
；
[0036]S4.3：计算嵌入函数的反函数，g1(y)、g0(y)分别表示的和反函数；
[0037]S4.4：将统计模型f
Y
(y
i
|w
l
＝
±
1)引入L
l
(y)进一步计算，其中f
X
(x)为指数柯西瑞利的概率密度函数。经计算，L
l
(y)表示为：
[0038][0039]其中G1(y
i
)、G0(y<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于指数柯西瑞利建模的数字水印检测方法，其特征在于，约定：I表示载体图像；N
b
表示目标子带分割的总块数；L表示N
b
块中能量排序最高的块数；n表示阶数；m表示重复度；B
l
(l＝1,2,...,L)表示第l个块中UDWT
‑
PCET幅值系数的集合；B
l
'表示含水印图像的第l个块中UDWT
‑
PCET幅值系数的集合；w
l
表示由{
‑
1，1}生成的伪随机序列；I'表示含水印图像；λ表示尺度参数；β表示形状参数；x
i
表示原始UDWT
‑
PCET幅值系数；y
i
表示含水印图像的UDWT
‑
PCET幅值系数；所述基于指数柯西瑞利建模的数字水印检测方法包括水印嵌入和水印提取：所述水印嵌入步骤如下：S1：初始设置：输入载体图像I并初始化设置；S2：水印嵌入：S2.1：对载体图像I进行二级非抽样离散小波变换(UDWT)，尺度二将获得三个高频子带，分别计算三个高频子能量，选择能量最大子带作为目标子带，能量计算公式如下：其中，X
j
[a,b]为第j个方向子带系数；S2.2：将目标子带划分为N
b
个不重叠且等大小的块，计算熵值并降序排列，对N
b
个块中熵值最高的前L个块进行n阶极坐标复指数变换(PCET)，以一块为例，选取[
‑
m,m]&amp;(m≠0)&amp;[1,n]区域为目标区域，计算得到L个目标PCET幅值系数集合B
l
(l＝1,2,...,L)用于水印信息的嵌入；S2.3：利用非线性嵌入强度函数修改每块UDWT
‑
PCET幅值系数以隐藏水印信息，水印嵌入表达式如下：其中f1(x)和f0(x)分别是嵌入水印为1和
‑
1时，对应的嵌入强度函数，f1(x)和f0(x)具体表示如下：表示如下：a1,b1,a2,b2为需要结合仿真实验确定的参数；S2.4：极坐标复指数矩重...

【专利技术属性】
技术研发人员：王向阳，马润桐，牛盼盼，
申请(专利权)人：辽宁师范大学，
类型：发明
国别省市：