CN103440615A - 一种基于余数不变性的高容量可逆水印方法 - Google Patents

一种基于余数不变性的高容量可逆水印方法 Download PDF

Info

Publication number
CN103440615A
CN103440615A CN2013104084529A CN201310408452A CN103440615A CN 103440615 A CN103440615 A CN 103440615A CN 2013104084529 A CN2013104084529 A CN 2013104084529A CN 201310408452 A CN201310408452 A CN 201310408452A CN 103440615 A CN103440615 A CN 103440615A
Authority
CN
China
Prior art keywords
block
sub
image
pixel
watermark
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013104084529A
Other languages
English (en)
Other versions
CN103440615B (zh
Inventor
翁韶伟
张天聪
潘正祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201310408452.9A priority Critical patent/CN103440615B/zh
Publication of CN103440615A publication Critical patent/CN103440615A/zh
Application granted granted Critical
Publication of CN103440615B publication Critical patent/CN103440615B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Editing Of Facsimile Originals (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于余数不变性的高容量可逆水印方法;包括水印嵌入过程及水印提取和原图像恢复过程；本发明提出的整数变换去掉了附加项，降低了嵌入失真，此整数变换仅由一个预测过程组成，它用块的均值去预测块中的每一个像素，仅让子块中最后一个预测误差前的每一个预测误差携带水印信息，而让最后一个像素仅被预测却不能携带任何水印信息，以保证整数变换的可逆性；利用图像子块的均值除以某一给定的嵌入参数所得到的余数在嵌入前后不发生变换的这个特性，整数变换可以无失真恢复出原始像素和水印信息；也因此，本发明无论是在嵌入容量上，还是在恢复图像质量上，都有了相当程度的提高。

Description

一种基于余数不变性的高容量可逆水印方法
技术领域
本发明涉及多媒体信号处理技术领域，具体地讲，涉及一种基于余数不变性的高容量可逆水印方法。
背景技术
传统的数字水印技术会造成宿主图像的永久性失真。但在一些实际应用中却不允许对宿主图像有一丁点的永久修改，比如医疗、军事和司法等领域。以医学图像为例，任何形式的失真都是不允许的。任何一幅医学图像的获取需要精密仪器的支持和昂贵的医疗费用， 更为重要的是失真可能会造成潜在的误诊。例如，对于一幅 ECG（electrocardiographic）信号图，任何一点信号曲线的异常都有可能被解释为某种病理特征。因此，传统的数字水印技术并不适用于医学图像。一种称为可逆水印的技术引起了越来越多人的研究与关注。可逆数字水印能以无损的方式将相关水印信息嵌入到宿主图像中，能在接收端有效地提取水印并精确恢复原始图像。
如何在保持载体图像视觉质量的情况下显著提高水印嵌入的数据容量，成为近年来图像可逆水印算法研究的主要方向之一。Alattar (A. M. Alattar, “Reversible watermark using the difference expansion of a generalized integer transform,” IEEE Trans. Image Process., vol. 13, no. 8, pp. 1147–1156, 2004.) 提出了一种广义差值扩展的可逆水印方法。Alatter提出的广义可逆整数小波变换可概括如下:
(1)
用来表示某一包含了个像素的图像子块，即
Figure 2013104084529100002DEST_PATH_IMAGE004
，
Figure 2013104084529100002DEST_PATH_IMAGE005
表示携带水印后的图像子块，
Figure 2013104084529100002DEST_PATH_IMAGE006
代表子块
Figure 204769DEST_PATH_IMAGE002
的均值，即
Figure 2013104084529100002DEST_PATH_IMAGE007
。假设
Figure 2013104084529100002DEST_PATH_IMAGE008
，代入到等式(1)可得：
(2)
从等式(2)可以看出，Alattar可分成两个部分：一个预测过程（用去预测
Figure 308860DEST_PATH_IMAGE002
中的每一个像素）和一个附加项
Figure 2013104084529100002DEST_PATH_IMAGE011
。Alattar有它自身的弱点：为了保证图像子块的
Figure 204442DEST_PATH_IMAGE010
在嵌入前后不发生变化，除了对每一个像素引入一个预测过程，它还必须要对每一个像素引入一个附加项，这个附加项的引入必定会在预测失真的基础上进一步增大嵌入失真。由于
Figure 2013104084529100002DEST_PATH_IMAGE012
和
Figure DEST_PATH_IMAGE013
（
Figure 2013104084529100002DEST_PATH_IMAGE014
），则
Figure DEST_PATH_IMAGE015
，即
Figure 2013104084529100002DEST_PATH_IMAGE016
。例如，当
Figure DEST_PATH_IMAGE017
，且
Figure 2013104084529100002DEST_PATH_IMAGE018
（
Figure 848919DEST_PATH_IMAGE014
），附加项的值可达到2。因此，由于这个附加项的引入，Alatter对每个像素额外多引入了最高可达2的嵌入失真，这将直接导致图像视觉质量的严重下降。
为了降低嵌入失真，必须去掉这个附加项，因此我们提出了一种基于余数不变性的可逆水印方法，以此提高嵌入容量和图像的视觉质量。而针对Wang (X. Wang, X. L. Li, B. Yang, and Z. M. Guo, “Efficient generalized integer transform for reversible watermarking,” IEEE Signal Processing Letters, vol. 17(6), pp. 567–570, 2010.)提出一种基于有效整数变换的可逆水印技术，在本发明中，我们也将给出本发明与Alattar和Wang在嵌入容量和视觉质量两个方面的比较结果。 
发明内容
本发明的主要目的在于，针对上述现有技术中的不足，提供一种基于余数不变性的高容量可逆水印方法，本方法可以用于医学图像、军事地图和军用图像的管理和授权使用信息中，以此为依据判别医学图像、军事地图和军用图像的来源，版权，真伪等信息。
本发明解决现有技术问题所采用的技术方案是：包括水印嵌入过程及水印提取和原图像恢复过程；
在水印嵌入过程中，先将原始图像分成互不重叠的多个图像子块，用图像子块的均值去预测每个图像子块中的每一个像素，使每个图像子块中的每一个像素得出一个预测误差；
    对每一个图像子块中最后一个预测误差前的每一个预测误差作差扩展，并对每一个扩展后的像素嵌入水印信息，而对每一个图像子块中的最后一个预测误差仅作差扩展而不对其嵌入水印信息，以保证整数变换的可逆性；
对于任一图像子块，用图像子块的方差来评估像素之间的相关性，按照方差的大小，将图像子块分成相关性高的子块和相关性弱的子块，对相关性高的子块的每一子块中的每个像素嵌入水印信息，而对相关性弱的子块不作任何修改；
水印提取和原图像恢复过程为水印嵌入过程的逆过程，通过利用图像子块的均值除以给定的嵌入参数所得到的余数在嵌入前后没有发生变换这一特性，来提取出原始水印信息和恢复出原始图像。
 下面对以上技术方案作进一步阐述：
 进一步地，所述水印嵌入过程包括以下步骤：
先将一副原始图像分成互不重叠的大小为
Figure DEST_PATH_IMAGE019
的图像子块，用
Figure 725608DEST_PATH_IMAGE002
表示某个
Figure 610388DEST_PATH_IMAGE019
大小的图像子块，即
Figure 2013104084529100002DEST_PATH_IMAGE020
；
用图像子块的均值去预测图像子块中的每一个像素，得到
Figure 942668DEST_PATH_IMAGE019
个预测误差；
再对图像子块前
Figure DEST_PATH_IMAGE021
个预测误差中的每一个预测误差进行
Figure 2013104084529100002DEST_PATH_IMAGE022
次差扩展，并嵌
Figure 203885DEST_PATH_IMAGE022
比特水印到每一个扩展后的像素中，
Figure DEST_PATH_IMAGE023
是某一给定的嵌入参数；而对图像子块最后一个像素
Figure 2013104084529100002DEST_PATH_IMAGE024
，仅作差扩展而不对其嵌入水印信息，以保证整数变换的可逆性；
 最后，对于任一图像子块，用图像子块的方差来评估像素之间的相关性，按照方差的大小，将子块分成高相关性的子块和相关性弱的子块，对于高相关性的子块，嵌
Figure DEST_PATH_IMAGE025
比特水印信息到块中的每一个像素；而对于相关性弱的子块不作任何修改。
本发明的有益效果是：本发明提出的整数变换去掉了附加项，因此，本发明所提出的整数变换本身可降低嵌入失真。此整数变换仅由一个预测过程组成，它用块的均值去预测块中的每一个像素，仅让子块中的前
Figure 997397DEST_PATH_IMAGE021
个预测误差携带水印信息，而让最后一个像素仅被预测却不能携带任何水印信息，因为一旦最后一个像素也携带水印，整数变换将不再具有可逆性。利用图像子块的均值除以某一给定的嵌入参数所得到的余数在嵌入前后不发生变换的这个特性，整数变换可以无失真恢复出原始像素和水印信息。也因此，本发明无论是在嵌入容量上，还是在恢复图像质量的主观评价上，都有了相当程度的提高。
附图说明
图1是本发明基于余数不变性的高容量可逆水印方法与Alattar、Wang方案的性能对比图；
其中：图（a１）是对于
Figure 2013104084529100002DEST_PATH_IMAGE026
大小的‘Lena’图像，本发明与Alattar方案、Wang方案的性能比较；图（b１）是对于
Figure 53078DEST_PATH_IMAGE026
大小的‘Baboon’图像，本发明与Alattar方案、Wang方案的性能比较；图（c１）是对于
Figure 869724DEST_PATH_IMAGE026
大小的‘Plane’图像，本发明与Alattar方案、Wang方案的性能比较；
图2是本发明基于余数不变性的高容量可逆水印方法与Wang方案的性能对比图；
其中：图（a２）是对于
Figure 606736DEST_PATH_IMAGE026
大小的‘Lena’图像，本发明与Wang方案的性能比较；图（b２）是对于
Figure 192438DEST_PATH_IMAGE026
大小的‘Baboo’图像，本发明与Wang方案的性能比较；图（c２）是对于
Figure 622282DEST_PATH_IMAGE026
大小的‘Plane’图像，本发明和Wang方案的性能比较。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
以下将结合附图及具体实施例详细说明本发明的技术方案，以便更清楚、直观地理解本发明的发明实质。
本发明的基于余数不变性的高容量可逆水印方法包括水印嵌入过程及水印提取和原图像恢复过程。
在本实施例的水印嵌入过程中，将一副大小为
Figure DEST_PATH_IMAGE027
的原始图像
Figure 2013104084529100002DEST_PATH_IMAGE028
分割成互不重叠的
Figure 991471DEST_PATH_IMAGE019
大小的子块，则有
Figure DEST_PATH_IMAGE029
个图像子块：
Figure 2013104084529100002DEST_PATH_IMAGE030
。用
Figure 328912DEST_PATH_IMAGE002
表示某个
Figure 34700DEST_PATH_IMAGE019
大小的图像子块，即
Figure 369866DEST_PATH_IMAGE020
。整数变换可看作是一预测过程，先用图像子块的均值去预测子块中的每一个像素，得到
Figure 161105DEST_PATH_IMAGE019
个预测误差，即
Figure DEST_PATH_IMAGE031
。对于前
Figure 302236DEST_PATH_IMAGE021
个预测误差中的每一个预测误差进行
Figure 862530DEST_PATH_IMAGE022
次差扩展，即
Figure DEST_PATH_IMAGE033
Figure 2013104084529100002DEST_PATH_IMAGE034
，并嵌比特水印到每一个扩展后的像素中，即
Figure DEST_PATH_IMAGE035
，其中
Figure 2013104084529100002DEST_PATH_IMAGE036
。例如，若
Figure DEST_PATH_IMAGE037
，则对预测误差进行一次差扩展，嵌1比特水印信息到每一个扩展后的预测误差中；若
Figure 2013104084529100002DEST_PATH_IMAGE038
，则对预测误差进行两次差扩展，嵌2比特水印信息到每一个扩展后的预测误差中。对于最后一个预测误差
Figure 2013104084529100002DEST_PATH_IMAGE039
，它仅被扩展却不能携带任何水印信息，因为一旦最后一个预测误差
Figure 837013DEST_PATH_IMAGE039
也携带水印，整数变换将不再具有可逆性。综上所述，整数变换提出如下：
Figure 2013104084529100002DEST_PATH_IMAGE040
 (3)
其中
Figure DEST_PATH_IMAGE041
表示携带水印后的图像子块，
Figure 844152DEST_PATH_IMAGE006
代表子块
Figure 196636DEST_PATH_IMAGE002
的均值，即
Figure 2013104084529100002DEST_PATH_IMAGE042
，
Figure 201501DEST_PATH_IMAGE023
代表给定的嵌入参数，
Figure DEST_PATH_IMAGE043
代表
Figure 701752DEST_PATH_IMAGE022
比特水印信息。本发明中，取
Figure 184686DEST_PATH_IMAGE038
，则
Figure 657256DEST_PATH_IMAGE043
代表两比特水印信息，即
Figure 2013104084529100002DEST_PATH_IMAGE044
，其中
Figure 833022DEST_PATH_IMAGE014
。在此整数变换中，水印嵌入前后，
Figure DEST_PATH_IMAGE045
保持不变，其中
Figure 2013104084529100002DEST_PATH_IMAGE046
返回
Figure DEST_PATH_IMAGE047
对
Figure 2013104084529100002DEST_PATH_IMAGE048
取模运算的结果，即
Figure DEST_PATH_IMAGE049
。
水印提取和原图像恢复过程即为水印嵌入过程的逆过程，就是整数变换的逆过程，描述如下：
分别计算子块
Figure 959852DEST_PATH_IMAGE002
的前
Figure 246477DEST_PATH_IMAGE021
个含水印像素
Figure 2013104084529100002DEST_PATH_IMAGE050
（
Figure DEST_PATH_IMAGE051
）和
Figure 2013104084529100002DEST_PATH_IMAGE052
的差值，消去相同项，得到等式(4)。这就是为什么
Figure 716959DEST_PATH_IMAGE052
不能携带任何水印信息的原因。
Figure DEST_PATH_IMAGE053
 (4)。
仔细观察等式(4)，可发现
Figure 191802DEST_PATH_IMAGE043
是一个小于等于的正整数（即），而
Figure DEST_PATH_IMAGE055
是
Figure 2013104084529100002DEST_PATH_IMAGE056
的整数倍。签于
Figure DEST_PATH_IMAGE057
有以上特点，可使
Figure 591264DEST_PATH_IMAGE057
（
Figure 843254DEST_PATH_IMAGE051
）对执行取模运算，得到模运算结果，即
Figure 2013104084529100002DEST_PATH_IMAGE058
。因此，通过使
Figure 371504DEST_PATH_IMAGE021
个差值对
Figure 735490DEST_PATH_IMAGE056
进行取模运算就可以正确提取出它们所携带的原始水印信息。正确提取出水印后，让
Figure 158381DEST_PATH_IMAGE050
（
Figure 545500DEST_PATH_IMAGE051
）减去它所携带的水印信息
Figure 977618DEST_PATH_IMAGE043
，可得
Figure DEST_PATH_IMAGE059
（具体参见等式(5)）。由于
Figure 2013104084529100002DEST_PATH_IMAGE060
没有携带任何的水印信息，则
Figure DEST_PATH_IMAGE061
。
Figure 2013104084529100002DEST_PATH_IMAGE062
 (5)
  假设
Figure DEST_PATH_IMAGE063
，其中
Figure 2013104084529100002DEST_PATH_IMAGE064
和
Figure DEST_PATH_IMAGE065
。从这个假设可以看出，使
Figure 120411DEST_PATH_IMAGE010
对
Figure 917466DEST_PATH_IMAGE056
执行取模运算，余数就是
Figure 2013104084529100002DEST_PATH_IMAGE066
，即
Figure DEST_PATH_IMAGE067
，将
Figure 650935DEST_PATH_IMAGE010
代入到等式(5)，可得如下等式：
Figure 2013104084529100002DEST_PATH_IMAGE068
 (6)
用
Figure DEST_PATH_IMAGE069
代替
Figure 2013104084529100002DEST_PATH_IMAGE070
以便进一步简化等式(6)。我们通过等式(7)来证明等式
Figure DEST_PATH_IMAGE071
。由此可知，在解码端，只需计算
Figure 2013104084529100002DEST_PATH_IMAGE072
就可以求得块的均值除以
Figure 22060DEST_PATH_IMAGE056
的余数
Figure 55262DEST_PATH_IMAGE066
。
Figure DEST_PATH_IMAGE073
 (7)
在正确得到
Figure 213711DEST_PATH_IMAGE066
后，用等式(8)来计算
Figure 2013104084529100002DEST_PATH_IMAGE074
（
Figure 987632DEST_PATH_IMAGE051
）和
Figure 852820DEST_PATH_IMAGE066
的差值，并用
Figure DEST_PATH_IMAGE075
（
Figure 53994DEST_PATH_IMAGE051
）表示差值。
(8)
  等式(8)的两边同除以4，可得
Figure DEST_PATH_IMAGE077
（）。将所有
Figure 11772DEST_PATH_IMAGE075
（
Figure 59362DEST_PATH_IMAGE051
）求和，可得：
Figure 2013104084529100002DEST_PATH_IMAGE078
 (9)
对等式(9)的两边同除以
Figure 431438DEST_PATH_IMAGE003
，再下取整，可得：
Figure DEST_PATH_IMAGE079
。在
Figure 2013104084529100002DEST_PATH_IMAGE080
和
Figure 555690DEST_PATH_IMAGE066
都被获得之后，就可恢复出
Figure 671413DEST_PATH_IMAGE010
，即
Figure DEST_PATH_IMAGE081
 (10)
将等式(10)恢复出来
Figure 573510DEST_PATH_IMAGE010
的代入到等式(5)即可恢复出来原始像素子块
Figure 116487DEST_PATH_IMAGE002
。
水印嵌入后，含水印的图像子块
Figure 2013104084529100002DEST_PATH_IMAGE082
中的每一含水印像素必须在
Figure DEST_PATH_IMAGE083
范围之内，否则含水印像素将遭受溢出（上溢
Figure 2013104084529100002DEST_PATH_IMAGE084
或下溢
Figure DEST_PATH_IMAGE085
），像素溢出将直接导致在解码端无法正确恢复出原始图像子块
Figure 533562DEST_PATH_IMAGE002
。为了避免像素溢出，我们定义
Figure 2013104084529100002DEST_PATH_IMAGE086
，其中。
更具体地，水印嵌入过程还包括以下步骤：
　　①像素分类：计算子块的方差
Figure 2013104084529100002DEST_PATH_IMAGE088
，即
Figure DEST_PATH_IMAGE089
。对于某一图像子块，水印嵌入后，如果
Figure 988704DEST_PATH_IMAGE082
中所有像素都落在
Figure 34021DEST_PATH_IMAGE083
区间之内，且，则属于
Figure DEST_PATH_IMAGE091
，即
Figure 2013104084529100002DEST_PATH_IMAGE092
；否则属于，即
Figure 2013104084529100002DEST_PATH_IMAGE094
，其中
Figure DEST_PATH_IMAGE095
是一个给定的门限，用来区分出
Figure 227291DEST_PATH_IMAGE091
和
Figure 377649DEST_PATH_IMAGE093
。
Figure 644682DEST_PATH_IMAGE091
中的每一个图像子块，携带
Figure 905900DEST_PATH_IMAGE025
个水印比特，而
Figure 309199DEST_PATH_IMAGE093
中的每一个像素子块保持不变。
②位置图的形成：用一个一维的位置图， 大小为，来标识这两类，位置图中符号‘1’表示中的像素，‘0’表示
Figure 450035DEST_PATH_IMAGE093
中的像素。用算术编码无损压缩位置图。压缩位置图产生的码流用
Figure DEST_PATH_IMAGE097
表示，假设其长度为
Figure 2013104084529100002DEST_PATH_IMAGE098
。因此待嵌入水印信息来自两个方面：1) 压缩后的码流
Figure 249364DEST_PATH_IMAGE097
；2) 载荷
Figure DEST_PATH_IMAGE099
。因此，嵌入容量（即载荷大小）
Figure 2013104084529100002DEST_PATH_IMAGE100
等于
Figure DEST_PATH_IMAGE101
，其中
Figure 2013104084529100002DEST_PATH_IMAGE102
用来表示集合的长度。
③水印信息嵌入过程：对于任一子块
Figure DEST_PATH_IMAGE103
，如果它属于
Figure 959700DEST_PATH_IMAGE091
， 则按照公式(3)嵌两比特水印信息到预测误差
Figure 451861DEST_PATH_IMAGE034
中；如果属于
Figure 490224DEST_PATH_IMAGE093
 ，则保持不变，即
Figure DEST_PATH_IMAGE105
。
④嵌压缩后的码流：在对前
Figure 2013104084529100002DEST_PATH_IMAGE106
中的子块按照步骤③处理完了后，收集
Figure DEST_PATH_IMAGE107
的
Figure 864278DEST_PATH_IMAGE098
个最不重要位（LSB），并把其最不重要位添加到
Figure 199444DEST_PATH_IMAGE099
的后面，用压缩后的码流
Figure 990683DEST_PATH_IMAGE097
代替空的
Figure 69497DEST_PATH_IMAGE098
个最不重要位（LSB）。全部图像子块被处理完了之后，得到水印图像
Figure DEST_PATH_IMAGE109
。
水印提取过程分为以下步骤：
　　第一步，位置图的恢复：收集所有像素的最不重要位，形成一个比特流。 在此比特流中找到算术编码的结束符， 从开始到结束符的比特流被解压缩得到原始位置图。
第二步，提取出水印和恢复出原始像素：对于某一个，即对于第一类水印像素如果它在位置图中所对应的是`0’，则保持不变，即
Figure 2013104084529100002DEST_PATH_IMAGE110
；如果是`1’，即第二类水印像素，则用整数变换的逆过程恢复出原始图像和提取出水印信息。
参照图１、图２给出的实例与实验结果所示：
图1呈现了对`Lena’、 `Baboon’和`Plane’三幅图像仅做了一次嵌入操作时三种方法所得到的实验结果。在所有的嵌入率中，本发明的PSNR值远远高于Alattar方案。在低嵌入率的时候，本发明所获得的PSNR值和Wang方案获得的PSNR几乎一样。随着嵌入率的逐渐增加，本发明的PSNR值在相同嵌入率下远高于Wang方案获得的PSNR。
图2分别呈现了对`Lena’、 `Baboon’和`Plane’三幅图像做了两次嵌入操作时本发明和Wang方案的实验结果。在两次嵌入中，在低嵌入率的时候，本发明所获得的PSNR值和Wang方案获得的PSNR几乎一样。随着嵌入率的逐渐增加，本发明的PSNR值在相同嵌入率下远高于Wang方案获得的PSNR。
大量的实验结果表明，与以往大部分的同类方案相比，我们所设计的基于余数不变性的高容量可逆水印算法，相对于Alattar的方案——无论是在嵌入容量上，还是在恢复图像质量的主观评价上，都有了相当程度的提高。其原因在于：整数变换没有引入附件项，从而有效地降低了嵌入失真。
以上所述仅为本发明的优选实施例，并非因此限制其专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (2)

1.一种基于余数不变性的高容量可逆水印方法，其特征在于：包括水印嵌入过程及水印提取和原图像恢复过程；
在水印嵌入过程中，先将原始图像分成互不重叠的多个图像子块，用图像子块的均值去预测每个图像子块中的每一个像素，使每个图像子块中的每一个像素得出一个预测误差；
 对每一个图像子块中最后一个预测误差前的每一个预测误差作差扩展，并对每一个扩展后的像素嵌入水印信息，而对每一个图像子块中的最后一个预测误差仅作差扩展而不对其嵌入水印信息，以保证整数变换的可逆性；
对于任一图像子块，用图像子块的方差来评估像素之间的相关性，按照方差的大小，将图像子块分成相关性高的子块和相关性弱的子块，对相关性高的子块的每一子块中的每个像素嵌入水印信息，而对相关性弱的子块不作任何修改；
水印提取和原图像恢复过程为水印嵌入过程的逆过程，通过利用图像子块的均值除以给定的嵌入参数所得到的余数在嵌入前后没有发生变换这一特性，来提取出原始水印信息和恢复出原始图像。
2.根据权利要求1所述的基于余数不变性的高容量可逆水印方法，其特征在于：
所述水印嵌入过程包括以下步骤：
先将一副原始图像分成互不重叠的大小为                                                
Figure 2013104084529100001DEST_PATH_IMAGE001
的图像子块，用表示某个
Figure 185383DEST_PATH_IMAGE001
大小的图像子块，即
Figure 2013104084529100001DEST_PATH_IMAGE003
；
用图像子块的均值去预测图像子块中的每一个像素，得到
Figure 669585DEST_PATH_IMAGE001
个预测误差；
再对图像子块前
Figure 843078DEST_PATH_IMAGE004
个预测误差中的每一个预测误差进行
Figure 2013104084529100001DEST_PATH_IMAGE005
次差扩展，并嵌
Figure 663879DEST_PATH_IMAGE005
比特水印到每一个扩展后的像素中，
Figure 694152DEST_PATH_IMAGE006
是某一给定的嵌入参数；而对图像子块最后一个像素
Figure 2013104084529100001DEST_PATH_IMAGE007
，仅作差扩展而不对其嵌入水印信息，以保证整数变换的可逆性；
最后，对于任一图像子块，用图像子块的方差来评估像素之间的相关性，按照方差的大小，将子块分成高相关性的子块和相关性弱的子块，对于高相关性的子块，嵌
Figure 665650DEST_PATH_IMAGE008
比特水印信息到块中的每一个像素；而对于相关性弱的子块不作任何修改。
CN201310408452.9A 2013-09-10 2013-09-10 一种基于余数不变性的高容量可逆水印方法 Expired - Fee Related CN103440615B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310408452.9A CN103440615B (zh) 2013-09-10 2013-09-10 一种基于余数不变性的高容量可逆水印方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310408452.9A CN103440615B (zh) 2013-09-10 2013-09-10 一种基于余数不变性的高容量可逆水印方法

Publications (2)

Publication Number Publication Date
CN103440615A true CN103440615A (zh) 2013-12-11
CN103440615B CN103440615B (zh) 2017-07-07

Family

ID=49694307

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310408452.9A Expired - Fee Related CN103440615B (zh) 2013-09-10 2013-09-10 一种基于余数不变性的高容量可逆水印方法

Country Status (1)

Country Link
CN (1) CN103440615B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933668A (zh) * 2015-06-11 2015-09-23 广东工业大学 一种基于局部定位和Alattar整数变换的可逆水印的处理方法
CN105741224A (zh) * 2016-01-28 2016-07-06 广东工业大学 基于pvo和自适应块分割的可逆水印算法
CN108898542A (zh) * 2018-07-04 2018-11-27 广东工业大学 一种可逆水印的嵌入和提取方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102044061A (zh) * 2010-12-20 2011-05-04 北京大学 一种可逆水印的嵌入方法及其提取方法
CN103050120A (zh) * 2012-12-28 2013-04-17 暨南大学 大容量数字音频可逆水印处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102044061A (zh) * 2010-12-20 2011-05-04 北京大学 一种可逆水印的嵌入方法及其提取方法
CN103050120A (zh) * 2012-12-28 2013-04-17 暨南大学 大容量数字音频可逆水印处理方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ADNAN M. ALATTAR: ""Reversible Watermark Using the Difference"", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
SHAOWEI WENG: ""Invariability of Mean Value Based Reversible"", 《JOURNAL OF INFORMATION HIDING AND MULTIMEDIA SIGNAL PROCESSING》 *
XIANG WANG: ""Efficient Generalized Integer Transform for"", 《IEEE SIGNAL PROCESSING LETTERS》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933668A (zh) * 2015-06-11 2015-09-23 广东工业大学 一种基于局部定位和Alattar整数变换的可逆水印的处理方法
CN104933668B (zh) * 2015-06-11 2019-03-15 广东工业大学 一种基于局部定位和Alattar整数变换的可逆水印的处理方法
CN105741224A (zh) * 2016-01-28 2016-07-06 广东工业大学 基于pvo和自适应块分割的可逆水印算法
CN105741224B (zh) * 2016-01-28 2019-03-19 广东工业大学 基于pvo和自适应块分割的可逆水印算法
CN108898542A (zh) * 2018-07-04 2018-11-27 广东工业大学 一种可逆水印的嵌入和提取方法

Also Published As

Publication number Publication date
CN103440615B (zh) 2017-07-07

Similar Documents

Publication Publication Date Title
CN105741224B (zh) 基于pvo和自适应块分割的可逆水印算法
CN106023058B (zh) 一种高保真可逆水印嵌入方法
Weng et al. Optimal PPVO-based reversible data hiding
CN102147912B (zh) 一种基于自适应差值扩展的可逆图像水印方法
CN101105857B (zh) 一种基于预测和压扩技术的高容量可逆水印方法
CN104933668B (zh) 一种基于局部定位和Alattar整数变换的可逆水印的处理方法
CN101866478B (zh) 一种数字图像中的水印嵌入和提取方法
CN108009975A (zh) 基于二维直方图修改的jpeg图像可逆信息隐藏方法
CN103268588A (zh) 基于载体图像校验码的加密域无损可逆信息隐藏方法
CN106485640A (zh) 一种基于多层次ipvo的可逆水印计算方法
Yang et al. Huffman-code strategies to improve MFCVQ-based reversible data hiding for VQ indexes
CN110047029B (zh) 一种结合多层差扩展具有对比度增强的可逆信息隐藏方法及装置
CN103440615A (zh) 一种基于余数不变性的高容量可逆水印方法
CN101246585B (zh) 一种在数字图像中嵌入和提取水印的方法与装置
Chu et al. A digital image watermarking method based on labeled bisecting clustering algorithm
KR101026081B1 (ko) 가역 워터마크 삽입, 추출 및 원본 영상 복원 방법
CN103440616B (zh) 一种基于自适应预测模型的高容量可逆水印方法
CN111275602A (zh) 人脸图像安全保护方法、系统及存储介质
CN114170059A (zh) 基于奇异值分解的两阶段图像鲁棒可逆信息隐藏方法
Weng et al. Invariability of Remainder Based Reversible Watermarking.
Yu et al. Breaking the YASS algorithm via pixel and DCT coefficients analysis
Saryazdi et al. A blind DCT domain digital watermarking
Wu et al. Anti-JPEG Compression Steganography Based on the High Tense Region Locating Method.
Ma et al. An improved reversible watermarking scheme
Venugopala et al. An approach for bitstream video watermarking on mobile device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170707

Termination date: 20200910

CF01 Termination of patent right due to non-payment of annual fee