CN105741224B - 基于pvo和自适应块分割的可逆水印算法 - Google Patents
基于pvo和自适应块分割的可逆水印算法 Download PDFInfo
- Publication number
- CN105741224B CN105741224B CN201610059600.4A CN201610059600A CN105741224B CN 105741224 B CN105741224 B CN 105741224B CN 201610059600 A CN201610059600 A CN 201610059600A CN 105741224 B CN105741224 B CN 105741224B
- Authority
- CN
- China
- Prior art keywords
- block
- image block
- pixel
- image
- 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.)
- Expired - Fee Related
Links
- 230000002441 reversible effect Effects 0.000 title claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 16
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims description 29
- 238000000605 extraction Methods 0.000 claims description 23
- 238000011084 recovery Methods 0.000 claims description 15
- 238000012986 modification Methods 0.000 claims description 14
- 230000004048 modification Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 230000001174 ascending effect Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 4
- 238000003780 insertion Methods 0.000 abstract description 16
- 230000037431 insertion Effects 0.000 abstract description 16
- 230000011218 segmentation Effects 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 description 9
- 239000000284 extract Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241001504519 Papio ursinus Species 0.000 description 1
- 241000758706 Piperaceae Species 0.000 description 1
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013506 data mapping Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0021—Image watermarking
- G06T1/005—Robust watermarking, e.g. average attack or collusion attack resistant
- G06T1/0071—Robust watermarking, e.g. average attack or collusion attack resistant using multiple or alternating watermarks
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Editing Of Facsimile Originals (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种基于PVO和自适应块分割的可逆水印算法，属于多媒体信息安全领域，该算法先根据局部纹理特征的复杂程度将图像块粗分成两类：平滑块和纹理块，再按照块内像素的相关程度将平滑块细分成两个子类：低平滑性的块和高平滑性的块。对于任一个低平滑性的块，按照像素值由低到高的顺序对其进行排序得到(pσ(1),pσ(2),…,pσ(n))，用pσ(n‑2)去预测pσ(n‑1)和pσ(n)或用pσ(3)去预测pσ(2)和pσ(1)，一个低平滑性的块至多能携带四比特水印信息。而对于高平滑性的块，将其分成只包含3个像素的子块，即(p1,p2,p3)，先对(p1,p2,p3)进行排序得到(pσ(1),pσ(2),pσ(3))，本发明通过仅修改pσ(3)来达到同时修改pσ(1),pσ(2)的目的，极大降低了嵌入失真。
Description
技术领域
本发明涉及一种多媒体信号处理领域，具体是指数字水印技术中的可逆水印算法。
背景技术
传统的数字水印技术会造成宿主图像的永久性失真。但在一些实际应用中却不允许对宿主图像有一丁点的永久修改，比如医疗、军事和司法等领域。以医学图像为例，任何形式的失真都是不允许的。任何一幅医学图像的获取需要精密仪器的支持和昂贵的医疗费用，更为重要的是失真可能会造成潜在的误诊。例如，对于一幅ECG(electrocardiographic)信号图，任何一点信号曲线的异常都有可能被解释为某种病理特征。因此，传统的数字水印技术并不适用于医学图像。一种称为可逆水印的技术引起了越来越多人的研究与关注。可逆数字水印能以无损的方式将相关水印信息嵌入到宿主图像中，能在接收端有效地提取水印并精确恢复原始图像。
如何在保持载体图像视觉质量的情况下显著提高水印嵌入的数据容量，成为近年来图像可逆水印算法研究的主要方向之一。Feng(F.Peng,X.L.Li,B.Yang,“Improved PVO-based reversible data hiding,”Digital Signal Processing,vol.25,pp.255-265,2014.)提出一种基于PVO和像素间位置关系的可逆水印技术。
在Peng等人的方法中，原始图像被分成大小为n＝r×c的图像块。对于任意一个图像块(p1,…,pn)(n≥3)，其全部像素按像素值的大小做升序排列得到(pσ(1),…,pσ(n))，其中σ:{1,...,n}→{1,...,n}是一对一的数据映射，它表示的是像素在排序之前的位置，注意：如果两个像素的值相等，则原先位置在前的像素排序后仍在前而原先位置在后的像素排序后仍在后，即如果pσ(i)＝pσ(j)且i＜j，则pσ(1)≤…≤pσ(n)，σ(i)＜σ(j)。
在考虑了pσ(n)和pσ(n-1)之间的位置关系后，产生出一个新的差值dmax。
dmax＝pu-pv (1)
其中
必须保证映射σ在嵌入前后维持不变，否则本发明就无法恢复出原始像素值和正确抽取出水印信息。注意：当pσ(n-1)＝pσ(n)时，σ(n-1)必定小于σ(n)。因此，对于dmax＝0且σ(n)＞σ(n-1)，d'max∈{-1,0}；对于dmax＜0，d'max∈(-∞,-2]；对于dmax＝1且σ(n)＜σ(n-1)，d'max∈{1,2}；对于dmax＞1，d'max∈[2,∞)，其中d'max是dmax经过扩展和平移后的值。因此，在水印嵌入之后,d'max互不重叠(见等式(3))，也就是说在解码端，只要计算出d'max的取值就可以无失真恢复出原始差值dmax。
其中b∈{0,1}代表1-比特水印信息。相应地，pσ(n)被修改为p'σ(n)：
在解码端，借助d'max＝p'u-p'v来恢复出原始像素pσ(n)和提取出原始水印信息。
●如果d'max＞0，则可知p'u＞p'v，进而推导出σ(n)＜σ(n-1)，u＝σ(n)和v＝σ(n-1)。
如果d'max∈{1,2}，提取的水印b＝d'max-1，恢复出的原始像素为pσ(n)＝p'u-b。
如果d'max＞2，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(n)＝p'u-1。
●如果d'max≤0，则可知pu≤pv，进而推导出σ(n)＞σ(n-1)，u＝σ(n-1)和v＝σ(n)。
如果d'max∈{0,-1}，提取的水印b＝-d'max，恢复出的原始像素为pσ(n)＝p'v-b。
如果d'max＜-1，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(n)＝p'v-1。
其中(u,v)来自于等式(1)。
事实上，将上述的核心思想直接应用到(pσ(1),pσ(2),pσ(3))中，就能得到最小值修改的嵌入方案，细节如下：
假设两个最小值之间的差值为dmin，其值计算如下：
dmin＝ps-pt (5)
其中s＝min(σ(1),σ(2))，t＝max(σ(1),σ(2))。
同理，最小值pσ(1)被修改为p'σ(1)：
在解码端，借助d'min＝p's-p't来恢复出原始像素pσ(1)和提取出原始水印信息。
●如果d'min＞0，则可知p's＞pt'，进而推导出σ(1)＞σ(2)，s＝σ(2)和t＝σ(1)。
如果d'min∈{1,2}，提取的水印b＝d'min-1，恢复出的原始像素为pσ(1)＝p't+b。
如果d'min＞2，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(1)＝p't+1。
●如果d'min≤0，则可知pu≤pv，进而推导出σ(n)＞σ(n-1)，u＝σ(n-1)和v＝σ(n)。
如果d'min∈{0,-1}，提取的水印b＝-d'min，恢复出的原始像素为pσ(1)＝p's+b。
如果d'min＜-1，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(1)＝p's+1。
其中s＝min(σ(1),σ(2))，t＝max(σ(1),σ(2))。
在Peng等人的算法中，用排序后的第二大像素pσ(n-1)去预测最大像素pσ(n)(或用第二小像素pσ(2)去预测最小像素pσ(1))，相对于最近邻预测方法(用每一个像素的右邻或左邻像素来预测这个像素)，预测性能大大提高。然而，Peng等人的算法认为所有的平滑块都具有相同的平滑性，因此对每一个平滑块都采用了相同的处理方式(即通过仅修改pσ(1)和pσ(n)来实现至多2比特的水印信息嵌入)。事实上，块的平滑性也是有很大差异的，平滑级别高的块理应比平滑级别低的块携带更多的水印比特。基于以上想法，本发明提出了一种基于PVO的自适应块分割算法。具体地说，根据块的局部复杂性将平滑块细分成两类：G1和G2，其中G1包含了局部复杂性高的全部平滑块，G2包含了局部复杂性低的全部平滑块。通过使得不同级别的图像块携带不同数量的水印信息，本发明在巨大提高嵌入容量的同时保持了较高的图像视觉效果。又由于本发明仅嵌水印信息在位于平滑区域的图像块中，所以能实现低嵌入率下的高PSNR值。
发明内容
本发明的目的是提出一种降低了嵌入失真的基于PVO的自适应水印嵌入方案。可以用于医学图像、军事地图和军用图像的管理和授权使用信息中，以此为依据判别医学图像、军事地图和军用图像的来源，版权，真伪等信息。
为实现上述目的，本发明所采用的技术方案是：
一种基于PVO(pixel-value-ordering，像素值排序)和自适应块分割的可逆水印算法，包括水印嵌入过程、水印提取及原始图像恢复过程，其中，水印提取及原始图像恢复过程是水印嵌入过程的逆过程，水印嵌入过程包括以下步骤：
①将一幅载体图像素分成互不重叠的图像块，每一个图像块包含n个像素，图像块的(n+1)个相邻像素组成一个像素集合，并计算该集合的方差；设定方差的阈值，将像素集合的方差与阈值比较，如果方差小于某个所述的阈值，则认为图像块属于平滑区域，否则认为图像块位于纹理区域；
②对于位于纹理区域的图像块，在水印嵌入过程中对它们不做任何修改；
对于位于平滑区域的图像块，分成两个部分：第一部分是水印嵌入后没有遭受像素溢出的图像块，第二部分是水印嵌入后遭受像素溢出的图像块；产生一个位置图来标识这两部分图像块的位置信息，并将无损压缩后的位置图和载荷一起嵌入到载体图像中；
③第一部分图像块的细分：按照局部纹理特征的复杂程度将第一部分的像素块分成两类：低平滑性的图像块和高平滑性的图像块；
④对于任意一个低平滑性的图像块，按像素值的大小做升序排列得到(pσ(1),pσ(2),…,pσ(n))，用块中第三大像素pσ(n-2)分别预测最大两个像素pσ(n-1)和pσ(n)得到两个预测误差(或用第三小像素pσ(3)分别预测最小两个像素pσ(2)和pσ(1)得到两个预测误差),因此，通过对四个预测误差进行修改可使得一个低平滑性的块能至多携带4个水印比特；对于高平滑性的图像块，将其细分成包含3个像素的子块，即(p1,p2,p3)，先对(p1,p2,p3)按像素值的大小做升序排列得到(pσ(1),pσ(2),pσ(3))，再用pσ(3)分别预测pσ(1),pσ(2)得到两个预测误差，当两个预测误差同时需要被平移1位时，通过仅修改pσ(3)来达到同时对pσ(1),pσ(2)做平移1位的修改，待所有属于第一部分的图像块经上述修改之后，得到含水印的图像。
进一步地，所述水印提取和原始图像恢复过程的具体方法是：①按照和嵌入过程相同的方式将含水印的图像分成互不重叠的图像块；②抽取压缩后的位置图，并对其进行无损恢复得到原始位置图；③采用和嵌入过程相反的顺序对块进行水印的抽取，即先嵌的后抽取，后嵌的先抽取，对于某一图像块，如果其所对应的方差大于某一给定阈值，则保持其不变，若其方差小于给定阈值，且位置图中所对应的位置标示的是‘0’，也不做任何处理，对于剩余的图像块，对其进行水印的抽取和原始图像的恢复。
进一步地，所述步骤①的具体内容如下：在水印嵌入过程中，将一副大小为R×C的原始图像I分割成互不重叠的n＝r×c大小的图像块，并以行的方向对每块进行编号，记为：对于其中一个图像块对它按照奇数行从左到右，偶数行从右到左的方式进行扫描,将二维的图像块转换成一维像素序列p＝{p1,…,pn},其全部像素按照值的大小排序后得到(pσ(1),…,pσ(n))；
对于某一图像块若存在(r+c+1)个像素环绕着它，则将这些像素(用p1,c+1,…,pr,c+1,pr+1,c+1,pr+1,1,…,pr+1,c表示)构成一个像素集IENP，并用σ表示像素集IENP的方差，由以下公式可计算出σ的值：
其中μENP表示集合IENP的均值，若方差σ小于某一给定的阈值vTh，则认为Bi,j和环绕它的所有像素之间存在着很强的相关性，即认为Bi,j属于平滑区域；否则认为两者之间的相关性弱，Bi,j属于纹理区域。
进一步地，所述步骤③的具体内容如下：
对于属于平滑区域的图像块，如果则它属于低平滑性的图像块G1；如果则它属于高平滑性的图像块G2。
本发明的有益效果：
由于采用上述的算法，具有以下优点：⑴块的平滑性是有很大差异的，平滑级别高的块理应比平滑级别低的块携带更多的水印比特，相对于平滑级别低的块，平滑级别高的块实现相同的嵌入容量时对图像引入的改动更小，从而实现低嵌入率下的高视觉效果；⑵降低了记载图像块的位置信息的位置图的大小。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明：
图1为512×512大小的‘Lena’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图2为512×512大小的‘Baboon’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图3为512×512大小的‘Barbara’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图4为512×512大小的‘Airplane’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图5为512×512大小的‘Boat’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图6为512×512大小的‘Sailboat’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图7为512×512大小的‘Peppers’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较；
图8为512×512大小的‘Elaine’为测试图像，本发明和Wang、Peng、Li、Sachnev和Hong的方案的性能比较。
具体实施方式
一种基于PVO和自适应块分割的可逆水印算法，包括水印嵌入过程、水印提取及原始图像恢复过程，其中，水印提取及原始图像恢复过程是水印嵌入过程的逆过程，其特征在于，水印嵌入过程包括以下步骤：
①将一幅载体图像素分成互不重叠的图像块，每一个图像块包含n个像素，图像块的(n+1)个相邻像素组成一个像素集合，并计算该集合的方差；设定方差的阈值，将像素集合的方差与阈值比较，如果方差小于某个所述的阈值，则认为图像块属于平滑区域，否则认为图像块位于纹理区域；
具体内容如下：在水印嵌入过程中，将一副大小为R×C的原始图像I分割成互不重叠的n＝r×c大小的图像块，并以行的方向对每块进行编号，记为：对于其中一个图像块对它按照奇数行从左到右，偶数行从右到左的方式进行扫描,将二维的图像块转换成一维像素序列p＝{p1,…,pn},其全部像素按照值的大小排序后得到(pσ(1),…,pσ(n))；
对于某一图像块若存在(r+c+1)个像素环绕着它，则将这些像素(用p1,c+1,…,pr,c+1,pr+1,c+1,pr+1,1,…,pr+1,c表示)构成一个像素集IENP，并用σ表示像素集IENP的方差，由以下公式可计算出σ的值：
其中μENP表示集合IENP的均值，若方差σ小于某一给定的阈值vTh，则认为Bi,j和环绕它的所有像素之间存在着很强的相关性，即认为Bi,j属于平滑区域；否则认为两者之间的相关性弱，Bi,j属于纹理区域。
定义图像块集合D＝{PS∈A:0≤p'i≤255(i∈{σ(1),…σ(n)})}，其中PS表示P按照像素值从大到小排序后的像素序列。A＝{PS＝(pσ(1),…,pσ(n)):0≤pi≤255(i∈{σ(1),…σ(n)})}。
按照局部纹理特征的复杂程度将图像块集合D分成两部分：SP和D-SP，其中SP包含了局部复杂性小于vTh的所有块(平滑块)，即SP＝{PS∈D,σ＜vTh}；D-SP包含了所有纹理块，即D-SP＝{PS∈D,σ≥vTh}。
②对于步骤①中，位于纹理区域的图像块，在水印嵌入过程中对它们不做任何修改；对于位于平滑区域的图像块，分成两个部分：第一部分是水印嵌入后没有遭受像素溢出的图像块，第二部分是水印嵌入后遭受像素溢出的图像块；产生一个位置图来标识这两部分图像块的位置信息，并将无损压缩后的位置图和载荷一起嵌入到载体图像中；
像素溢出是指，水印嵌入后，含水印的像素必须在[0,255]范围之内，否则含水印像素将遭受溢出(上溢(＞255)或下溢(＜0))，像素溢出将直接导致在解码端无法正确恢复出原始图像子块P。
位置图生成：当σ＜vTh时，用一个大小为的一维二值位置图，来区分一个平滑块是否遭受到了像素溢出，具体地说，位置图中符号‘1’表示SP中的图像块，‘0’表示(遭受了像素溢出的图像块)中的图像块。这里有一点需要注意：在解码端，只要判断出σ≥vTh，就可以知道哪些图像块属于D-SP，因此没有必要记载D-SP中图像块的位置信息。用算术编码无损压缩位置图。压缩位置图产生的码流用表示，假设其长度为LS。因此待嵌入水印信息来自两个方面：1)压缩后的码流2)载荷
③第一部分图像块的细分：按照局部纹理特征的复杂程度将第一部分的像素块分成两类：低平滑性的图像块和高平滑性的图像块；对于属于平滑区域的图像块，如果则它属于低平滑性的图像块G1；如果则它属于高平滑性的图像块G2。
④对于任意一个低平滑性的图像块G1，按像素值的大小做升序排列得到(pσ(1),pσ(2),…,pσ(n))，用块中第三大像素pσ(n-2)分别预测最大两个像素pσ(n-1)和pσ(n)得到两个预测误差(或用第三小像素pσ(3)分别预测最小两个像素pσ(2)和pσ(1)得到两个预测误差),因此，通过对四个预测误差进行修改可使得一个低平滑性的块能至多携带4个水印比特；
下面详解具体的水印数据嵌入方案:对于属于G1的任一个平滑块，以它的三个最大像素(pσ(n-2),pσ(n-1),pσ(n))为例，本发明用pσ(n-2)分别预测pσ(n-1)和pσ(n)，经等式dmax＝pu-pv(1)之后得到两个差值d1max和d2max。d1max是考虑了pσ(n-2)和pσ(n-1)之间的位置关系(即σ(n-2)和σ(n-1)之间的前后关系)而得到的pσ(n-2)和pσ(n-1)之间的差值，d2max是考虑了pσ(n-2)和pσ(n)之间的位置关系(即σ(n-2)和σ(n)之间的前后关系)而得到的pσ(n-2)和pσ(n)之间的差值。即使对于pσ(n-2)＝pσ(n-1)＝pσ(n)的块，本发明也可通过分别修改d1max和d2max来实现至多两比特水印信息的嵌入。
如将pσ(n)被修改为p'σ(n)，具体修改过程，不再展开详述，由于该修改采用的是现有技术，可参见等式(4)如下：
只要保证在水印嵌入之后，含水印的像素值(即p'σ(n-2)和p'σ(n-1))比pσ(n-3)大，就可以无失真恢复出原始像素值，其中p'σ(n-2)是pσ(n-2)的含水印像素值，同理，p'σ(n-1)是pσ(n-1)的含水印像素值。
以它的三个最小像素(pσ(1),pσ(2),pσ(3))为例，本发明用pσ(3)分别预测pσ(1)和pσ(2)，经等式(5)之后得到两个差值d1min和d2min。d1min是考虑了pσ(3)和pσ(2)之间的位置关系(即σ(3)和σ(2)之间的前后关系)而得到的pσ(3)和pσ(2)之间的差值，d2min是考虑了pσ(3)和pσ(1)之间的位置关系(即σ(3)和σ(1)之间的前后关系)而得到的pσ(3)和pσ(1)之间的差值。通过这种方式，即使对于pσ(1)＝pσ(2)＝pσ(3)的块，本发明也可通过分别修改d1min和d2min来实现至多两比特水印信息的嵌入。对pσ(1)和pσ(2)的具体嵌入过程可参见背景技术中的等式(6)如下
对于属于G2的平滑块：对于高平滑性的图像块，将其细分成只包含3个像素的子块，即(p1,p2,p3)，先对(p1,p2,p3)进行排序得到(pσ(1),pσ(2),pσ(3))，用pσ(3)分别预测pσ(1),pσ(2)得到两个预测误差，当两个预测误差同时需要被平移1位时，通过仅修改pσ(3)来达到同时对pσ(1),pσ(2)做平移1位的修改，待所有属于第一部分的图像块经上述修改之后，得到含水印的图像。
本发明将大小为r×c的图像块重新细分成大小为1×3的子块，对于任意一个子块(p1,p2,p3)，按照像素值从小到大的顺序排列后得到(pσ(1),pσ(2),pσ(3))，本发明用pσ(3)去分别预测pσ(1)和pσ(2)，经等式dmax＝pu-pv(1)之后得到两个差值d1max和d2max，其中d1max是考虑了pσ(3)和pσ(2)之间的位置关系(即σ(3)和σ(2)之间的前后关系)而得到的pσ(3)和pσ(2)之间的差值，d2max是考虑了pσ(3)和pσ(1)之间的位置关系(即σ(3)和σ(1)之间的前后关系)而得到的pσ(3)和pσ(2)之间的差值。对(pσ(1),pσ(2),pσ(3))的具体修改过程可参考表1。
通过表1可知，本发明的创新点在于把两个最小像素看成了一个整体(即像素对(pσ(1),pσ(2)))来处理，当pσ(1)和pσ(2)同时需要被向左平移一位时(即p'σ(1)＝pσ(1)+w1且w1＝1或p'σ(1)＝pσ(1)+1，同时，p'σ(2)＝pσ(2)+w2且w2＝1或p'σ(2)＝pσ(2)+1)，在这种情况下，本发明通过仅修改pσ(3)来间接地达到同时修改pσ(1)和pσ(2)的目的，这样做的好处在于：在容量维持不变的基础上进一步降低了嵌入过程所引入的失真。
其中w1和w2代表1比特水印信息，即w1∈{0,1}和w2∈{0,1}。对(p'σ(1),p'σ(2),p'σ(3))的水印抽取和原始图像恢复过程可参考表2。
对于前个图像块中的每一个像素按照上述步骤修改完之后，收集P中像素的最不重要位(LSB)，并把收集到的LS个最不重要位添加到载荷的后面。剩余的图像块，对其处理的过程和前个图像块类似。但不同的是其P的最不重要位不需要收集。待全部图像块处理完了之后，得到水印图像IW。
水印提取和原始图像恢复过程：
抽取过程分为以下三个步骤：
(1)按照和嵌入过程相同的方式将含水印的图像分成互不重叠的图像块：按照与嵌入过程相同的顺序，将含水印的图像IW分成大小为n的互不重叠的图像块。并以行的方向对每块进行编号，记为：对于其中一个含水印的图像块对它按照奇数行从左到右,偶数行从右到左的方式进行扫描,将二维的图像块转换成一维像素序列P'＝{p'1,p'2…,p'n}；
(2)抽取压缩后的位置图，并对其进行无损恢复得到原始位置图：收集前个像素序列中像素的最不重要位，形成一个比特流。在此比特流中找到算术编码的结束符，从开始到结束符的比特流被解压缩得到原始位置图。
(3)采用和嵌入过程相反的顺序对块进行水印的抽取，即先嵌的后抽取，后嵌的先抽取，对于某一图像块，如果其所对应的方差大于某一给定阈值，则保持其不变，若其方差小于给定阈值，且位置图中所对应的位置标示的是‘0’，也不做任何处理，对于剩余的图像块，对其进行水印的抽取和原始图像的恢复；
具体方案如下：按照和嵌入过程相反的顺序进行水印信息的抽取和原始图像的恢复，即每一个块对应一个一维像素序列。对于和如果其没有(r+c+1)相邻的像素,则保持不变。否则，对于某一像素序列P'，将其所有相邻像素p1,c+1,…,pr,c+1,pr+1,c+1,pr+1,1,…,pr+1,c构成和嵌入过程中相同的IENP(注意：p1,c+1,…,pr,c+1,pr+1,c+1,pr+1,1,…,pr+1,c必须要先于P'恢复出来，才能保证从P'中正确恢复出原始像素序列P)。IENP经等式(7)，即：得到P'所对应的σ，如果σ大于或等于vTh，则保持不变；如果σ小于vTh，且在位置图中所对应的是‘0’，则保持不变；如果在位置图中所对应的是‘1’且则按照表2对P'进行水印的提取和原始像素的恢复；
如果则按照如下方式进行水印的抽取和原始像素的恢复，在解码端，借助d'max＝p'u-p'v来分别恢复出原始像素pσ(n-1)和pσ(n)，并提取出原始水印信息。
·如果d'max＞0，则可知p'u＞p'v，进而推导出σ(n)＜σ(n-1)，u＝σ(n)和v＝σ(n-1)。
如果d'max∈{1,2}，提取的水印b＝d'max-1，恢复出的原始像素为pσ(n)＝p'u-b。
如果d'max＞2，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(n)＝p'u-1。
·如果d'max≤0，则可知pu≤pv，进而推导出σ(n)＞σ(n-1)，u＝σ(n-1)和v＝σ(n)。
如果d'max∈{0,-1}，提取的水印b＝-d'max，恢复出的原始像素为pσ(n)＝p'v-b。
如果d'max＜-1，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(n)＝p'v-1。
其中(u,v)来自于等式(1)。
在解码端，借助d'min＝p's-p't来分别恢复出原始像素pσ(1)和pσ(2)，并提取出原始水印信息。
·如果d'min＞0，则可知p's＞p't，进而推导出σ(1)＞σ(2)，s＝σ(2)和t＝σ(1)。
如果d'min∈{1,2}，提取的水印b＝d'min-1，恢复出的原始像素为pσ(1)＝p't+b。
如果d'min＞2，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(1)＝p't+1。
·如果d'min≤0，则可知pu≤pv，进而推导出σ(n)＞σ(n-1)，u＝σ(n-1)和v＝σ(n)。
如果d'min∈{0,-1}，提取的水印b＝-d'min，恢复出的原始像素为pσ(1)＝p's+b。
如果d'min＜-1，说明没有隐藏任何水印信息，恢复出的原始像素为pσ(1)＝p's+1。
其中s＝min(σ(1),σ(2))，t＝max(σ(1),σ(2))。
本发明的性能远远优于Wang的方法的性能。虽然Wang的方法也是基于块分割和PVO的方法。但是Wang的方法是先将宿主图像分成4×4大小的图像块，再按照局部纹理复杂程度将图像块分成三类：①纹理块；②一般复杂性的平滑块；③低复杂性的平滑块。对于低复杂性的平滑块，Wang等人将4×4的块进行再次分割分成大小为2×2的子块。对于每一个大小为2×2的子块，通过对最大值和最小值进行等式(4)和(6)的修改使得每一个2×2大小的子块能至多携带两个比特水印信息。然而，Wang等人在实验中并没有给出任何证据来证明4×4的图像块能实现最优的嵌入性能。事实上，图像块的大小不应该只限定于4×4，应该为任意大小。而本发明的优势在于可以将宿主图像分成大小为r×c的块。对于属于G2的图像子块，本发明将其分成大小为1×3的子块。即使是r＝4，c＝4，本发明仍能将其分成大小为1×3的子块，只不过4×4的块只能分成5个大小为1×3的子块。
从图1-8也可以看出，对于八幅测试图像，本发明的性能远优于其它五种方法(即Wang的方法、Peng的方法，Li的方法，Sachnev的方法以及Hong的方法)。对于‘Lena’，‘Barbara’和‘Airplane’图像，当嵌入容量较低时(例如，对于‘Lena’图像，嵌入容量低于19000比特)，Wang的方法的性能无法优于Peng的方法的性能。而本发明的优势就在于，块的局部纹理特征越复杂，它所能携带的数据容量就越低。也就是说，块越平滑，它所能携带的水印信息越多。因此，甚至在低嵌入比特率的时候，本发明的PSNR值远远优Wang的方法和Peng方法的PSNR值。
大量的实验结果表明，与以往大部分的同类方案相比，本发明所设计的基于PVO和自适应块分割的可逆水印算法——无论是在嵌入容量上，还是在恢复图像质量的主观评价上，都有了相当程度的提高。其原因在于：对于属于G2的图像子块，采用表1所示的像素修改方法对其进行水印嵌入，在保证嵌入容量的同时进一步降低了嵌入失真；对于属于G1的图像子块，通过同时修改两个最大像素和两个最小像素来达到一个大小为r×c的图像块能至多携带4比特水印信息。
以上所述是本发明的优选实施方式而已，当然不能以此来限定本发明之权利范围，应当指出，对于本技术领域的普通技术人员来说，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的保护范围。
Claims (4)
1.一种基于PVO和自适应块分割的可逆水印算法，包括水印嵌入过程、水印提取及原始图像恢复过程，其中，水印提取及原始图像恢复过程是水印嵌入过程的逆过程，其特征在于，水印嵌入过程包括以下步骤：
①将一幅载体图像素分成互不重叠的r×c大小的图像块，对于每一个图像块，将环绕该图像块的右邻r+1个像素和下邻c个像素构成一个像素集合，并计算该集合的方差σ；设定阈值vTh，将像素集合的方差与阈值比较，如果方差σ小于阈值vTh，则认为图像块属于平滑区域，否则认为图像块位于纹理区域；
②对于位于纹理区域的图像块，在水印嵌入过程中对它们不做任何修改；
对于位于平滑区域的图像块，分成两个部分：第一部分是水印嵌入后没有遭受像素溢出的图像块，第二部分是水印嵌入后遭受像素溢出的图像块；产生一个位置图来标识这两部分图像块的位置信息，并将无损压缩后的位置图和载荷一起嵌入到载体图像中；
③第一部分图像块的细分：按照局部纹理特征的复杂程度将第一部分的像素块分成两类：低平滑性的图像块和高平滑性的图像块；
④对于任意一个低平滑性的图像块，按像素值的大小做升序排列得到(pσ(1),pσ(2),…,pσ(n))，用块中第三大像素pσ(n-2)分别预测最大两个像素pσ(n-1)和pσ(n)得到两个预测误差，同时，用第三小像素pσ(3)分别预测最小两个像素pσ(2)和pσ(1)得到两个预测误差,因此，通过对四个预测误差进行修改可使得一个低平滑性的块能至多携带4个水印比特；对于高平滑性的图像块，将其细分成包含3个像素的子块，即(p1,p2,p3)，先对(p1,p2,p3)按像素值的大小做升序排列得到(pσ(1),pσ(2),pσ(3))，再用pσ(3)分别预测pσ(1),pσ(2)得到两个预测误差，当两个预测误差同时需要被平移1位时，通过仅修改pσ(3)来达到同时对pσ(1),pσ(2)做平移1位的修改，待所有属于第一部分的图像块经上述修改之后，得到含水印的图像。
2.根据权利要求1所述的基于PVO和自适应块分割的可逆水印算法，其特征在于：所述水印提取和原始图像恢复过程的具体方法是：①按照和嵌入过程相同的方式将含水印的图像分成互不重叠的图像块；②抽取压缩后的位置图，并对其进行无损恢复得到原始位置图；③采用和嵌入过程相反的顺序对图像块进行水印的抽取，即先嵌的后抽取，后嵌的先抽取，对于某一图像块，如果其所对应的方差大于某一给定阈值，则保持其不变，若其方差小于给定阈值，且位置图中所对应的位置标示的是‘0’，也不做任何处理，对于剩余的图像块，对其进行水印的抽取和原始图像的恢复。
3.根据权利要求1所述的基于PVO和自适应块分割的可逆水印算法，其特征在于：
所述步骤①的具体内容如下：在水印嵌入过程中，将一副大小为R×C的原始图像I分割成互不重叠的r×c大小的图像块，并以行的方向i和列的方向j对每块进行编号，记为：对于其中一个图像块对它按照奇数行从左到右，偶数行从右到左的方式进行扫描,将二维的图像块转换成一维像素序列p＝{p1,…,pr×c},其全部像素按照值的大小排序后得到(pσ(1),…,pσ(r×c))；
对于某一图像块若该图像块有r+1个右邻像素p1,c+1,…,pr,c+1,pr+1,c+1和c个下邻像素pr+1,1,…,pr+1,c，则将这些邻域像素构成一个像素集IENP，并用σ表示像素集IENP的方差，由以下公式可计算出σ的值：
其中μENP表示集合IENP的均值，若方差σ小于某一给定的阈值vTh，则认为Bi,j和环绕它的所有像素之间存在着很强的相关性，即认为Bi,j属于平滑区域；否则认为两者之间的相关性弱，Bi,j属于纹理区域。
4.根据权利要求1所述的基于PVO和自适应块分割的可逆水印算法，其特征在于：所述步骤③的具体内容如下：
对于属于平滑区域的图像块，如果则它属于低平滑性的图像块G1；如果则它属于高平滑性的图像块G2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610059600.4A CN105741224B (zh) | 2016-01-28 | 2016-01-28 | 基于pvo和自适应块分割的可逆水印算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610059600.4A CN105741224B (zh) | 2016-01-28 | 2016-01-28 | 基于pvo和自适应块分割的可逆水印算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105741224A CN105741224A (zh) | 2016-07-06 |
| CN105741224B true CN105741224B (zh) | 2019-03-19 |
Family
ID=56246959
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610059600.4A Expired - Fee Related CN105741224B (zh) | 2016-01-28 | 2016-01-28 | 基于pvo和自适应块分割的可逆水印算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105741224B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106485640A (zh) * | 2016-08-25 | 2017-03-08 | 广东工业大学 | 一种基于多层次ipvo的可逆水印计算方法 |
| CN106778280A (zh) * | 2016-11-02 | 2017-05-31 | 北京知道未来信息技术有限公司 | 一种填充式远程漏洞PoC编写方法及漏洞检测方法 |
| CN107241597B (zh) * | 2017-05-05 | 2019-10-15 | 中国人民武装警察部队工程大学 | 一种结合四叉树自适应编码的可逆信息隐藏方法 |
| CN108022197B (zh) * | 2017-11-21 | 2021-03-12 | 安徽大学 | 一种基于多纹理度划分的医学图像可逆信息隐藏方法 |
| CN108109101A (zh) * | 2017-12-20 | 2018-06-01 | 广东工业大学 | 一种基于多媒体数据的信息隐藏方法、装置及设备 |
| CN108898542B (zh) * | 2018-07-04 | 2023-07-18 | 广东工业大学 | 一种可逆水印的嵌入方法 |
| CN109522429B (zh) * | 2018-10-18 | 2022-02-01 | 百度在线网络技术（北京）有限公司 | 用于生成信息的方法和装置 |
| CN109389546B (zh) * | 2018-10-30 | 2023-07-25 | 广东工业大学 | 多分块可逆水印嵌入方法及装置、水印提取方法及装置 |
| CN109859090B (zh) * | 2018-12-21 | 2023-07-25 | 广东工业大学 | 基于人类视觉系统的可逆水印方法及装置 |
| CN109741233B (zh) * | 2018-12-29 | 2024-01-23 | 广东工业大学 | 一种可逆水印的嵌入和提取方法 |
| CN110910299B (zh) * | 2019-11-29 | 2021-01-05 | 山东科技大学 | 一种基于整数小波变换的自适应可逆信息隐藏方法 |
| CN111179145B (zh) * | 2019-12-10 | 2023-04-25 | 中国科学院深圳先进技术研究院 | 一种双层可逆水印方法 |
| CN113099068B (zh) * | 2021-04-09 | 2022-05-10 | 华南师范大学 | 基于双图像的大容量像素值排序可逆信息隐藏方法 |
| CN113393364A (zh) * | 2021-07-07 | 2021-09-14 | 福建工程学院 | 一种基于像素值排序的可逆信息隐藏方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004003841A2 (en) * | 2002-06-27 | 2004-01-08 | University Of Bristol | Image attack characterisation by double watermarking |
| CN101105857A (zh) * | 2007-07-20 | 2008-01-16 | 北京交通大学 | 一种基于预测和压扩技术的高容量可逆水印方法 |
| JP2008283479A (ja) * | 2007-05-10 | 2008-11-20 | Fujifilm Software Co Ltd | 可逆的透かしデータ挿入方法、プログラム、装置および画像処理システム |
| CN103440616A (zh) * | 2013-09-10 | 2013-12-11 | 广东工业大学 | 一种基于自适应预测模型的高容量可逆水印方法 |
| CN103440615A (zh) * | 2013-09-10 | 2013-12-11 | 广东工业大学 | 一种基于余数不变性的高容量可逆水印方法 |
| CN104751402A (zh) * | 2015-04-16 | 2015-07-01 | 西安电子科技大学 | 一种鲁棒可逆水印嵌入的方法 |
| CN104933668A (zh) * | 2015-06-11 | 2015-09-23 | 广东工业大学 | 一种基于局部定位和Alattar整数变换的可逆水印的处理方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7319754B2 (en) * | 2002-06-28 | 2008-01-15 | Stmicroelectronics S.A. | Insertion of binary messages in video pictures |
-
2016
- 2016-01-28 CN CN201610059600.4A patent/CN105741224B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004003841A2 (en) * | 2002-06-27 | 2004-01-08 | University Of Bristol | Image attack characterisation by double watermarking |
| JP2008283479A (ja) * | 2007-05-10 | 2008-11-20 | Fujifilm Software Co Ltd | 可逆的透かしデータ挿入方法、プログラム、装置および画像処理システム |
| CN101105857A (zh) * | 2007-07-20 | 2008-01-16 | 北京交通大学 | 一种基于预测和压扩技术的高容量可逆水印方法 |
| CN103440616A (zh) * | 2013-09-10 | 2013-12-11 | 广东工业大学 | 一种基于自适应预测模型的高容量可逆水印方法 |
| CN103440615A (zh) * | 2013-09-10 | 2013-12-11 | 广东工业大学 | 一种基于余数不变性的高容量可逆水印方法 |
| CN104751402A (zh) * | 2015-04-16 | 2015-07-01 | 西安电子科技大学 | 一种鲁棒可逆水印嵌入的方法 |
| CN104933668A (zh) * | 2015-06-11 | 2015-09-23 | 广东工业大学 | 一种基于局部定位和Alattar整数变换的可逆水印的处理方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于直方图平移的大容量可逆水印算法研究;李可;《中国优秀硕士学位论文全文数据库（信息科技辑）》;20140115;第2014卷(第01期);第1-43页 |
| 数字图像的高容量可逆水印的研究;翁韶伟;《中国博士学位论文全文数据库（信息科技辑）》;20091015;第2009卷(第10期);第1-101页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105741224A (zh) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105741224B (zh) | 基于pvo和自适应块分割的可逆水印算法 | |
| CN106023058B (zh) | 一种高保真可逆水印嵌入方法 | |
| Qin et al. | Visible watermark removal scheme based on reversible data hiding and image inpainting | |
| Thabit et al. | Medical image authentication using SLT and IWT schemes | |
| JP5373909B2 (ja) | 圧縮ビデオにおける構文要素のｃａｂａｃ／ａｖｃ準拠の透かし入れ | |
| Weng et al. | Reversible data hiding based on flexible block-partition and adaptive block-modification strategy | |
| CN104933668B (zh) | 一种基于局部定位和Alattar整数变换的可逆水印的处理方法 | |
| Weng et al. | Optimal PPVO-based reversible data hiding | |
| CN108416724B (zh) | 一种基于多预测器的可逆信息隐藏方法 | |
| CN109741233B (zh) | 一种可逆水印的嵌入和提取方法 | |
| CN102044061A (zh) | 一种可逆水印的嵌入方法及其提取方法 | |
| CN106485640A (zh) | 一种基于多层次ipvo的可逆水印计算方法 | |
| CN101105857A (zh) | 一种基于预测和压扩技术的高容量可逆水印方法 | |
| CN109859090B (zh) | 基于人类视觉系统的可逆水印方法及装置 | |
| CN109523453A (zh) | 基于菱形预测与图像像素排序的可逆信息隐藏及提取方法 | |
| Weng et al. | Pairwise IPVO-based reversible data hiding | |
| CN104156909B (zh) | 基于图像内容的鲁棒可逆水印嵌入与检测方法 | |
| Thilagavathi et al. | A survey of reversible watermarking techniques, application and attacks | |
| CN108510425B (zh) | 基于ippvo和优化mhm的可逆水印方法 | |
| Su et al. | Reversible data hiding using the dynamic block-partition strategy and pixel-value-ordering | |
| CN110533569B (zh) | 基于二次差值扩展的水印处理方法 | |
| CN111127288B (zh) | 可逆图像水印处理方法、装置及计算机可读存储介质 | |
| CN103440616A (zh) | 一种基于自适应预测模型的高容量可逆水印方法 | |
| VijayKumar et al. | Histogram shifting based reversible data hiding with controlled contrast enhancement | |
| Arsalan et al. | Intelligent threshold selection for reversible watermarking of medical images |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190319 |