基于LWT算法和AES加密算法的水印嵌入和提取方法与流程

本发明涉及水印技术领域，更具体的涉及基于lwt算法和aes加密算法的水印嵌入和提取方法。

背景技术：

信息技术发展的趋势增强了数据传输安全性的需要。通信信道的数据传输为了数据的安全目的必须对数据进行加密。因特网的出现已经改变了人们的生活方式，在因特网上人们可以交互他们的多媒体文本和数据，如：语音、数据和图像等。一个突出的问题就是版权保护和认证问题。近些年数字图像水印技术已经提出来用于多媒体文本和数据的版权保护。

现有的水印嵌入方法是基于dwt的水印嵌入方法，但是dwt精确度低、计算速度慢，且没有对水印图像进行加密，从而使得水印的安全性低。

综上，现有的水印嵌入方法存在精确度低、计算速度慢以及安全性低的问题。

技术实现要素：

本发明实施例提供基于lwt算法和aes加密算法的水印嵌入和提取方法，用以解决现有技术中精确度低、计算速度慢以及安全性低的问题。

本发明实施例提供基于lwt算法和aes加密算法的水印嵌入和提取方法，包括：水印嵌入方法和水印提取方法；

所述水印嵌入方法包括：

利用提升小波变换lwt将原始图像分解为ll子带、lh子带、hl子带和hh子带；

将ll子带划分成多个l×l大小的模块k，并对每一个模块k进行奇异值分解，获得模块k的奇异值λi；

将水印图像进行aes加密，并将加密后的图像d进行奇异值分解，获得加密后的图像d的奇异值λwi；

根据如下公式修正模块k的奇异值；

其中，α为修正因子，为模块k修正后的奇异值；

基于所述模块k修正后的奇异值构成的奇异值矩阵对模块k进行重构，得到嵌入水印的模块，所有嵌入水印的模块构成嵌入水印的ll子带；

将嵌入水印的ll子带、lh子带、hl子带和hh子带做lwt逆变换，获得嵌入水印图像。

所述水印提取方法包括：

利用提升小波变换lwt将嵌入水印图像进行分解，获得提升小波变换的ll'子带；

将ll'子带划分成多个l×l大小的模块k^d，并对每一个模块k^d进行奇异值分解，获得模块k^d的奇异值；

按照如下公式，确定aes加密水印图像的奇异值；

其中，α为修正因子，为aes加密水印图像的奇异值，为模块k^d的奇异值，λi为模块k的奇异值；

基于aes加密水印图像的奇异值构成的奇异值矩阵，对aes加密水印图像进行重构，得到重构的aes加密水印图像；

将重构的aes加密水印图像进行aes解密，获得原始水印图像。

较佳地，将水印图像分解成8×8个子块，对每一个子块采用aes加密算法进行加密时，采用128bit的密钥，并且对每一个子块均进行10轮的数据加密。

较佳地，每一轮的数据加密具体包括：字节代换、行位移、列混合和轮密钥加。

较佳地，所述提升小波变换，具体包括：

分裂步骤：将输入信号x(n)分解为偶点和奇点样值序列xe(n)和xo(n)；

预测步骤：利用偶点序列xe(n)的预测值p[xe(n)]来预测奇点序列xo(n)，将差值d(n)＝xo(n)-p[xe(n)]定义为信号x(n)的高频分量，其中，p[·]为预测算子；

更新步骤：利用高频分量d(n)来更新偶点序列xe(n)，近似低频分量c(n)＝xe(n)+u[d(n)]，其中u[·]为更新算子。

较佳地，所述奇异值分解变换，具体包括：

尺寸大小为m×n像素的图像的奇异值分解表达式可以展成：c＝usv′，其中，u是一个正交矩阵u′u＝i，s＝diag(λi)为奇异值λi,i＝1,2...,r的对角矩阵，并且奇异值以降序排列，v也为正交矩阵v′v＝i；u矩阵的列向量u为左边奇异值向量，v矩阵的列向量v为右边奇异值向量，则图像的奇异值分解表达式也可以通过下式表示：

本发明实施例中提供的基于lwt算法和aes加密算法的水印嵌入和提取方法，由于本发明水印嵌入方法是基于提升小波变换(lwt)和高级加密标准(aes)数据加密技术，lwt精确度高，计算速度快，运算效率高，是一种整数到整数的变换；因此，提高了水印的运算速率，且水印图像采用aes数据加密，数据更安全，为水印技术又多了一层保护，提高了水印的安全性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的基于lwt算法和aes加密算法的水印嵌入的流程示意图；

图2为本发明实施例提供的提升小波变换的分解框图；

图3为本发明实施例提供的aes加密算法流程图；

图4为本发明实施例提供的基于lwt算法和aes加密算法的水印提取方法的流程示意图；

图5(1)为本发明实施例提供的原始图像；

图5(2)为本发明实施例提供的嵌入水印图像；

图6(1)为本发明实施例提供的水印图像；

图6(2)为本发明实施例提供的aes加密图像；

图6(3)为本发明实施例提供的aes解密图像；

图7(1)为本发明实施例提供受到高斯噪声干扰时提取的水印图像；

图7(2)为本发明实施例提供的受到椒盐噪声干扰时提取的水印图像；

图7(3)为本发明实施例提供的受到斑纹噪声干扰时提取的水印图像。

具体实施方式

下面结合附图，对本发明的一个具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。

图1示例性的示出了本发明实施例提供的基于lwt算法和aes加密算法的水印嵌入方法的流程示意图，该基于lwt算法和aes加密算法的水印嵌入方法：

步骤1、利用提升小波变换lwt将原始图像分解为ll子带、lh子带、hl子带和hh子带。

其中，离散小波变换已经成功地应用于图像压缩和图像处理领域中。由于传统小波系数的浮点类型，使得在图像重构时引入舍入误差失去了图像的精确性。提升小波变换在直接分析整数域问题时简化可逆性问题，精确度高。提升小波方案采用了原位运算，运算速度快，运算效率高，并且所需存储空间少。该算法分为如下三步：分裂、预测和更新。提升小波变换的分解框图如图2所示,重构过程为分解的逆过程。

另外：该提升小波变换，具体包括：

分裂步骤：将输入信号x(n)分解为偶点和奇点样值序列xe(n)和xo(n)。

预测步骤：利用偶点序列xe(n)的预测值p[xe(n)]来预测奇点序列xo(n)，将差值d(n)＝xo(n)-p[xe(n)]定义为信号x(n)的高频分量，其中，p[·]为预测算子。

更新步骤：利用高频分量d(n)来更新偶点序列xe(n)，近似低频分量c(n)＝xe(n)+u[d(n)]，其中u[·]为更新算子。

步骤2、将ll子带划分成多个l×l大小的模块k，并对每一个模块k进行奇异值分解，获得模块k的奇异值λi。

另外，该奇异值分解，具体包括：

尺寸大小为m×n像素的图像的奇异值分解表达式可以展成：c＝usv′，其中，u是一个正交矩阵u′u＝i，s＝diag(λi)为奇异值λi,i＝1,2...,r的对角矩阵，并且奇异值以降序排列，v也为正交矩阵v′v＝i；u矩阵的列向量u为左边奇异值向量，v矩阵的列向量v为右边奇异值向量，则图像的奇异值分解表达式也可以通过下式(1)表示：

具体地，对每一个模块k进行奇异值分解，获得模块k的奇异值λi，按照如下公式(2)进行：

k＝usv′

s＝diag(λi)(2)

其中，u，v均是根据模块k求得。

步骤3、将水印图像进行aes加密，并将加密后的图像d进行奇异值分解，获得加密后的图像d的奇异值λwi。

其中，将加密后的图像d进行奇异值分解，获得加密后的图像d的奇异值λwi，按照如下公式(3)进行：

d＝uwswvw′，

其中，sw＝diag(λwi)，uw，vw均是根据加密后的图像d求得。

具体地，aes和des算法的安全性非常好。aes算法是对称的，在加密和解密过程中公用一个密钥，而des算法是非对称的；aes算法中选择长度128bit、192bit和256bit的密钥，aes算法的密钥比des更难破解。初始密钥经过密钥扩展形成轮密钥(roundkey),每一轮数据包括四个操作：字节代换、行位移、列混合和轮密钥加。图3为aes加密算法流程图(第k轮算法)。.

步骤4、根据如下公式(4)修正模块k的奇异值。

其中，α为修正因子，为模块k修正后的奇异值。

步骤5、基于该模块k修正后的奇异值构成的奇异值矩阵对模块k进行重构，得到嵌入水印的模块，所有嵌入水印的模块构成嵌入水印的ll子带。

其中，基于该修正后模块k的奇异值构成的奇异值矩阵对模块k进行重构，得到嵌入水印的模块时，具体按照如下公式(5)进行重构。

其中，为嵌入水印的模块，为奇异值构成的奇异值矩阵，

步骤6、将嵌入水印的ll子带、lh子带、hl子带和hh子带做lwt逆变换，获得嵌入水印图像。

可选地，本发明实施例，将水印图像分解成8×8个大小相同的子块，对每一个子块采用aes加密算法进行加密时，采用128bit的密钥，并且对每一个子块均进行10轮的数据加密。

图4示例性的示出了本发明实施例提供的基于lwt算法和aes加密算法的水印提取方法的流程示意图，该基于lwt算法和aes加密算法的水印提取方法包括：

步骤s1、利用提升小波变换lwt将嵌入水印图像进行分解，获得提升小波变换的ll'子带。

步骤s2、将ll'子带划分成多个l×l大小的模块k^d，并对每一个模块k^d进行奇异值分解，获得模块k^d的奇异值。

需要说明的是，对模块k^d进行奇异值分解，获得模块k^d的奇异值的方法与获得第一奇异值的方法类似，本发明实施例不再详细阐述。

步骤s3、按照如下公式，确定aes加密水印图像的奇异值；

其中，α为修正因子，为aes加密水印图像的奇异值，为模块k^d的奇异值，λi为原始图像矩阵ll子带模块，也即模块k的奇异值。

步骤s4、基于aes加密水印图像的奇异值构成的奇异值矩阵，对aes加密水印图像进行重构，得到重构的aes加密水印图像。

其中，基于aes加密水印图像的奇异值构成的奇异值矩阵，对aes加密水印图像进行重构，得到重构的aes加密水印图像时，按照如下公式(7)进行重构。

其中，为重构的aes加密水印图像，uw和vw为嵌入过程中所保留的密钥。

步骤s5、将重构的aes加密水印图像进行aes解密，获得原始水印图像。

本发明实施例提供的基于lwt算法和aes加密算法的水印嵌入和提取方法的实验仿真及其性能：

仿真环境为matlab7.0.选择数字图像处理标准测试图库中的标准图像，尺寸大小为256×256像素，如图5(1)所示，水印图像尺寸大小为32×32像素，如图6(1)所示。首先将水印图像分解成8×8个子块，每子块大小为4×4像素。选择128bit初始密钥进行密钥扩展生成轮密钥，每个子块数据经过10轮的字节代换、行位移、列混合和轮密钥加形成灰度水印图像的aes加密图像，如图6(2)所示。再将原始图像进行lwt变换分解为四个子带，取出ll子带，大小为128×128像素.提取ll子带和aes水印图像的主值，按照水印嵌入方法嵌入水印。嵌入水印图像如图5(2)所示。水印提取按照相逆的算法进行。图6(3)为提取的水印aes解密图像。

随着网路技术的发展，网络数据加密的要求日益提高，aes的应用首先体现在网络信息安全领域中。灰度水印图像经过aes数据加密后比eds数据加密更难破解，aes算法抗旁道攻击以及其它攻击性能非常强。因此该水印的安全性非常好。

其中，水印的性能由其性能指标峰值信噪比psnr和相关系数r衡量。峰值信噪比其中，maxi为原始图像中的最大值像素，mse为原始图像和嵌入水印图像的均方误差。相关系数其中，w为原始水印图像，为所提取的水印图像，为原始水印图像的均值，为所提取水印图像的均值。原始图像如图5(1)所示，嵌入水印图像如图5(2)所示，水印图像如图6(1)所示，其加密图像如图6(2)所示，提取的水印图像如图6(3)所示，此时psnr值非常大，均方误差非常小，相关系数r＝1，由此可以看出这足以保证嵌入水印图像的透明性。当遭受攻击时，psnr值和相关系数r要相应地减小。

图7(1)为遭受高斯噪声攻击时提取的水印图像，噪声均值为0，方差为0.0001，此时psnr＝39.95db,相关系数r＝0.6572。图7(2)为遭受椒盐噪声攻击时提取的水印图像，椒盐密度为0.0001，此时psnr＝46.38db,相关系数r＝0.7366。图7(3)为遭受斑纹噪声攻击时提取的水印图像，噪声方差为0.0001，此时psnr＝45.44db,相关系数r＝0.7017.由此可以看出该水印的鲁棒性。

综上该，本发明利用提升小波理论运算效率高、精确度高的优点和灰度水印图像aes加密数据更安全的特点，提出了基于提升小波变换的aes数据加密水印技术，以提高水印的鲁棒性和安全性。

本发明提供基于lwt算法和aes加密算法的水印嵌入和提取方法，由于本发明水印嵌入方法，是基于提升小波变换(lwt)和高级加密标准(aes)数据加密技术，lwt精确度高，计算速度快，运算效率高，是一种整数到整数的变换；因此，提高了水印的运算速率，且水印图像采用aes数据加密，数据更安全，为水印技术又多了一层保护，提高了水印的安全性。

尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。