基于二元T分布的数字水印检测方法

基于二元t分布的数字水印检测方法
技术领域
1.本发明涉及多媒体信息安全技术领域，具体为基于二元t分布的数字水印检测方法。


背景技术：

2.多媒体技术的发展，不仅使人们获取资源的途径更加多样化，同时也带来了众多安全隐患，例如：网站的非法入侵、信息的非法传播、版权安全等问题。因此，多媒体版权保护和内容认证问题亟待解决。图像水印技术为多媒体安全领域提供了有效的技术手段。
3.数字图像水印技术包括两部分，分别是水印嵌入和水印提取。水印嵌入是通过嵌入策略将水印嵌入到图像中，既保留了图像的使用价值，又将版权信息隐藏到图像中。水印提取则是利用技术手段从含水印图像中提取水印信息，从而证明该图像的所有权归属。数字图像水印方案的优劣性需要从不可感知性、水印容量和鲁棒性三个方面进行综合评价，从信息学角度讲，三者具有固有的矛盾关系。因此，保持三者平衡是一项具有挑战性的工作。近年来，随着数字图像水印技术的飞速发展，众多具有创新性的水印方案陆续被提出，但实现不可感知性、水印容量、鲁棒性三者之间的平衡问题仍是数字图像水印技术领域的一大难题。
4.为实现不可感知性、水印容量、鲁棒性之间的良好平衡，数字图像水印技术领域通常先利用变换方法对图像进行处理，再使用乘性嵌入策略嵌入水印。这是因为变换域对常见信号攻击类型具有鲁棒性，而乘性嵌入方法的优点在于，水印的嵌入强度能随着图像信号的强弱改变，二者的互相结合可以在一定程度上保持鲁棒性与不可感知性的平衡。基于统计模型的数字图像水印方案的核心是利用统计分布有效地提取水印，进而确定图像所属权。然而，近年来对图像水印技术的深入研究发现，传统的基于统计模型的水印方案存在一些不足之处：首先，仅使用变换域对图像信号进行处理，其鲁棒性还远远不足；其次，在当前大多数水印方案中，图像信号间的相关特性未被充分利用；最后，当前大多数水印方案利用传统最大似然估计方法估计模型参数，其参数精度相对较低，估计时间相对较长。


技术实现要素：

5.本部分的目的在于概述本发明的实施方式的一些方面以及简要介绍一些较佳实施方式。在本部分以及本技术的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊，而这种简化或省略不能用于限制本发明的范围。
6.为解决上述技术问题，根据本发明的一个方面，本发明提供了如下技术方案：
7.基于二元t分布的数字水印检测方法，
8.约定：i和i'分别表示原始图像信号和含水印图像信号，n1×
n1表示图像尺寸；dj[x,y]表示非下采样剪切波变换(nsst)域系数，j为子带方向；w
l
(l＝1,2,...,l)代表水印序列，l是水印长度；q
l
为奇异值系数矩阵；bk为高熵块，k＝1,2,...,l；xi和yi分别表示原始
图像和含水印图像的nsst域奇异值系数；α表示水印的嵌入强度；p表示维数，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设；
[0009]
所述基于二元t分布的数字水印检测方法包括水印嵌入和水印提取过程，
[0010]
所述水印嵌入流程如下：
[0011]
s1：初始设置：通过初始化设置得到原始图像信号i；
[0012]
s2：水印嵌入：
[0013]
s2.1：采用nsst对尺寸大小为n1×
n1的原始图像i进行分解，设置参数为[2 3]，得到一个低频子带以及第一尺度8个方向高频子带和第二尺度4个方向高频子带；
[0014]
s2.2：计算第二尺度能量最大高频子带：
[0015][0016]dj
[x,y]为nsst域系数，j表示子带方向；
[0017]
s2.3：将第二尺度能量最大的高频子带划分成不重叠的大小为8
×
8的图像子块，按照熵值大小降序排列，选择前l个高熵块b1,b2,...,b
l
进行奇异值分解，得到l个奇异值系数集合q
l
(l＝1,2,...,l)；
[0018]
s2.4：通过乘性水印嵌入策略修改奇异值集合q
l
(l＝1,2,...,l)中的奇异值系数，得到含水印的奇异值集合；
[0019]
s2.5：进行svd重构和逆nsst变换获取含水印图像i'；
[0020]
所述水印提取流程如下所示：
[0021]
s3：基于二元t分布建立模型：
[0022]
s3.1：采用nsst对尺寸为n1×
n1的含水印图像i'进行分解，设置参数为[23]，得到一个低频子带以及第一尺度8个方向和第二尺度4个方向高频子带；
[0023]
s3.2：分别选取第一尺度和第二尺度高频子带中能量最大的高频方向子带，将这两个高频子带分别划分成不重叠的大小为8
×
8的图像子块，找到嵌入水印时记录的前l个高熵块b
l
(l＝1,2,...,l)进行奇异值分解，得到奇异值系数集合来构造最大似然检测方法实现水印提取；
[0024]
s3.3：对含水印的奇异值系数集合利用二元t分布进行统计建模，假设p维随机向量服从多元t分布，则多元t分布的联合概率密度函数为：
[0025][0026]
其中：
[0027][0028]
其中，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；
[0029]
s3.4：通过对数矩估计方法估计参数。首先假设向量x＝(x1,x2)
t
服从二元t分布，
求样本x的r阶中心绝对阶矩，然后对mr求导，接着令r＝0，最后综合估计表达式计算参数值；
[0030]
s3.5：将参数估计值保存，为构造检测方法实现水印提取做准备；
[0031]
s4：设计最大似然解码器：
[0032]
s4.1：根据统计理论，水印提取过程被认为是一个二值假设检验问题，设h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设，表示为：
[0033][0034]
其中，xi和yi分别为原始图像和含水印图像的nsst域奇异值系数，α表示嵌入强度参数；
[0035]
s4.2：设计最大似然解码器：
[0036][0037]
s4.3：根据基于二元t分布的ml解码器，进行具体水印信息提取：
[0038][0039]
s4.4：根据最大似然检测方法，从含水印图像中提取水印序列。
[0040]
作为本发明所述的基于二元t分布的数字水印检测方法的一种优选方案，其中，所述步骤2.4中乘性水印嵌入策略如下所示：
[0041][0042]
其中，xi和yi分别表示原始图像和含水印图像的nsst域奇异值系数，α表示水印的嵌入强度。
[0043]
与现有技术相比，本发明具有的有益效果是：第一，采用参数为[2 3]的nsst和svd相结合构造建模对象，这种方式一定程度上提高了水印方案的不可感知性和鲁棒性；
[0044]
第二，引入二元t分布描述nsst域奇异值系数，一方面，可以准确地模拟出尖峰重尾的特征，另一方面，可以有效捕捉尺度间相关性；
[0045]
第三，使用对数矩估计方法估计二元t分布模型参数，该方法降低了算法的复杂度，计算简单，可提供更精确的估计值。
附图说明
[0046]
为了更清楚地说明本发明实施方式的技术方案，下面将结合附图和详细实施方式对本发明进行详细说明，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。其中：
[0047]
图1为本发明实施例中nsst域奇异值系数尺度间相关性分析图；
[0048]
图2为本发明实施例中二元t分布(实线)拟合nsst域奇异值双变量经验pdf(虚线)的等值线图；
[0049]
图3为本发明实施例中原始图像与含1024位水印图像；
[0050]
图4为本发明实施例中含水印图像及所提取的水印图；
[0051]
图5为本发明实施例中含1024位水印图像受到各种攻击的误码率变化图；
[0052]
图6为本发明水印嵌入流程图；
[0053]
图7为本发明水印提取流程图。
具体实施方式
[0054]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。
[0055]
其次，本发明结合示意图进行详细描述，在详述本发明实施方式时，为便于说明，表示器件结构的剖面图会不依一般比例作局部放大，而且所述示意图只是示例，其在此不应限制本发明保护的范围。此外，在实际制作中应包含长度、宽度及深度的三维空间尺寸。
[0056]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明的实施方式作进一步地详细描述。
[0057]
本发明方法操作过程为：水印嵌入、基于二元t分布建立模型及设计ml解码器。
[0058]
请参阅图6-7，基于二元t分布的数字水印检测方法，
[0059]
约定：i和i'分别表示原始图像信号和含水印图像信号，n1×
n1表示图像尺寸；dj[x,y]表示非下采样剪切波变换(nsst)域系数，j为子带方向；w
l
(l＝1,2,...,l)代表水印序列，l是水印长度；q
l
为奇异值系数矩阵；bk为高熵块，k＝1,2,...,l；xi和yi分别表示原始图像和含水印图像的nsst域奇异值系数；α表示水印的嵌入强度；p表示维数，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设；
[0060]
所述基于二元t分布的数字水印检测方法包括水印嵌入和水印提取过程，
[0061]
所述水印嵌入流程如下：
[0062]
s1：初始设置：通过初始化设置得到原始图像信号i；
[0063]
s2：水印嵌入：
[0064]
s2.1：采用nsst对尺寸大小为n1×
n1的原始图像i进行分解，设置参数为[2 3]，得到一个低频子带以及第一尺度8个方向高频子带和第二尺度4个方向高频子带；
[0065]
s2.2：计算第二尺度能量最大高频子带：
[0066]
[0067]dj
[x,y]为nsst域系数，j表示子带方向；
[0068]
s2.3：将第二尺度能量最大的高频子带划分成不重叠的大小为8
×
8的图像子块，按照熵值大小降序排列，选择前l个高熵块b1,b2,...,b
l
进行奇异值分解，得到l个奇异值系数集合q
l
(l＝1,2,...,l)；
[0069]
s2.4：通过乘性水印嵌入策略修改奇异值集合q
l
(l＝1,2,...,l)中的奇异值系数，得到含水印的奇异值集合；
[0070]
s2.5：进行svd重构和逆nsst变换获取含水印图像i'；
[0071]
所述水印提取流程如下所示：
[0072]
s3：基于二元t分布建立模型：
[0073]
s3.1：采用nsst对尺寸为n1×
n1的含水印图像i'进行分解，设置参数为[23]，得到一个低频子带以及第一尺度8个方向和第二尺度4个方向高频子带；
[0074]
s3.2：分别选取第一尺度和第二尺度高频子带中能量最大的高频方向子带，将这两个高频子带分别划分成不重叠的大小为8
×
8的图像子块，找到嵌入水印时记录的前l个高熵块b
l
(l＝1,2,...,l)进行奇异值分解，得到奇异值系数集合来构造最大似然检测方法实现水印提取；
[0075]
s3.3：对含水印的奇异值系数集合利用二元t分布进行统计建模，假设p维随机向量服从多元t分布，则多元t分布的联合概率密度函数为：
[0076][0077]
其中：
[0078][0079]
其中，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；
[0080]
s3.4：通过对数矩估计方法估计参数。首先假设向量x＝(x1,x2)
t
服从二元t分布，求样本x的r阶中心绝对阶矩，然后对mr求导，接着令r＝0，最后综合估计表达式计算参数值；
[0081]
s3.5：将参数估计值保存，为构造检测方法实现水印提取做准备；
[0082]
s4：设计最大似然解码器：
[0083]
s4.1：根据统计理论，水印提取过程被认为是一个二值假设检验问题，设h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设，表示为：
[0084][0085]
其中，xi和yi分别为原始图像和含水印图像的nsst域奇异值系数，α表示嵌入强度参数；
[0086]
s4.2：设计最大似然解码器：
[0087][0088]
s4.3：根据基于二元t分布的ml解码器，进行具体水印信息提取：
[0089][0090]
s4.4：根据最大似然检测方法，从含水印图像中提取水印序列。
[0091]
其中，所述步骤2.4中乘性水印嵌入策略如下所示：
[0092][0093]
其中，xi和yi分别表示原始图像和含水印图像的nsst域奇异值系数，α表示水印的嵌入强度。
[0094]
实施例1
[0095]
实验测试和参数设置：
[0096]
实验环境：matlab r2011a，nsst参数设置为[2 3]。
[0097]
图1为nsst域奇异值分解尺度间相关性分析图；
[0098]
图2为二元t分布(实线)拟合nsst域奇异值双变量经验pdf(虚线)的等值线图；
[0099]
图3为原始图像与含1024位水印图像；
[0100]
图4为含水印图像及所提取的水印图；
[0101]
图5为各种攻击的误码率变化图。
[0102]
实验结果表明，本发明的方法能够利用二元t分布捕捉nsst域奇异值的尺度间相关性，实现精准建模，在保证不可感知性的同时提高了鲁棒性。
[0103]
虽然在上文中已经参考实施方式对本发明进行了描述，然而在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，本发明所披露的实施方式中的各项特征均可通过任意方式相互结合起来使用，在本说明书中未对这些组合的情况进行穷举性的描述仅仅是出于省略篇幅和节约资源的考虑。因此，本发明并不局限于文中公开的特定实施方式，而是包括落入权利要求的范围内的所有技术方案。

技术特征：
1.基于二元t分布的数字水印检测方法，其特征在于：约定：i和i'分别表示原始图像信号和含水印图像信号，n1×
n1表示图像尺寸；d
j
[x,y]表示非下采样剪切波变换(nsst)域系数，j为子带方向；w
l
(l＝1,2,...,l)代表水印序列，l是水印长度；q
l
为奇异值系数矩阵；b
k
为高熵块，k＝1,2,...,l；x
i
和y
i
分别表示原始图像和含水印图像的nsst域奇异值系数；α表示水印的嵌入强度；p表示维数，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设；所述基于二元t分布的数字水印检测方法包括水印嵌入和水印提取过程，所述水印嵌入流程如下：s1：初始设置：通过初始化设置得到原始图像信号i；s2：水印嵌入：s2.1：采用nsst对尺寸大小为n1×
n1的原始图像i进行分解，设置参数为[2 3]，得到一个低频子带以及第一尺度8个方向高频子带和第二尺度4个方向高频子带；s2.2：计算第二尺度能量最大高频子带：d
j
[x,y]为nsst域系数，j表示子带方向；s2.3：将第二尺度能量最大的高频子带划分成不重叠的大小为8
×
8的图像子块，按照熵值大小降序排列，选择前l个高熵块b1,b2,...,b
l
进行奇异值分解，得到l个奇异值系数集合q
l
(l＝1,2,...,l)；s2.4：通过乘性水印嵌入策略修改奇异值集合q
l
(l＝1,2,...,l)中的奇异值系数，得到含水印的奇异值集合；s2.5：进行svd重构和逆nsst变换获取含水印图像i'；所述水印提取流程如下所示：s3：基于二元t分布建立模型：s3.1：采用nsst对尺寸为n1×
n1的含水印图像i'进行分解，设置参数为[23]，得到一个低频子带以及第一尺度8个方向和第二尺度4个方向高频子带；s3.2：分别选取第一尺度和第二尺度高频子带中能量最大的高频方向子带，将这两个高频子带分别划分成不重叠的大小为8
×
8的图像子块，找到嵌入水印时记录的前l个高熵块b
l
(l＝1,2,...,l)进行奇异值分解，得到奇异值系数集合来构造最大似然检测方法实现水印提取；s3.3：对含水印的奇异值系数集合利用二元t分布进行统计建模，假设p维随机向量服从多元t分布，则多元t分布的联合概率密度函数为：其中：
其中，ν表示形状参数，μ表示均值向量，σ表示协方差矩阵；s3.4：通过对数矩估计方法估计参数。首先假设向量x＝(x1,x2)
t
服从二元t分布，求样本x的r阶中心绝对阶矩，然后对m
r
求导，接着令r＝0，最后综合估计表达式计算参数值；s3.5：将参数估计值保存，为构造检测方法实现水印提取做准备；s4：设计最大似然解码器：s4.1：根据统计理论，水印提取过程被认为是一个二值假设检验问题，设h1和h0分别为在nsst域奇异值中嵌入“1”和
“‑
1”的假设，表示为：其中，x
i
和y
i
分别为原始图像和含水印图像的nsst域奇异值系数，α表示嵌入强度参数；s4.2：设计最大似然解码器：s4.3：根据基于二元t分布的ml解码器，进行具体水印信息提取：s4.4：根据最大似然检测方法，从含水印图像中提取水印序列。2.根据权利要求1所述的基于二元t分布的数字水印检测方法，其特征在于，所述步骤2.4中乘性水印嵌入策略如下所示：其中，x
i
和y
i
分别表示原始图像和含水印图像的nsst域奇异值系数，α表示水印的嵌入强度。

技术总结
本发明公开了基于二元T分布的数字水印检测方法,首先，利用NSST进行参数为[2 3]的分解，将尺度二能量最大高频子带划分成不重叠的大小为8


技术研发人员：牛盼盼 代钰 王向阳
受保护的技术使用者：辽宁师范大学
技术研发日：2023.04.26
技术公布日：2023/7/12