CN105894436A - 一种基于Gibbs抽样的图像隐写方法 - Google Patents
一种基于Gibbs抽样的图像隐写方法 Download PDFInfo
- Publication number
- CN105894436A CN105894436A CN201610200934.9A CN201610200934A CN105894436A CN 105894436 A CN105894436 A CN 105894436A CN 201610200934 A CN201610200934 A CN 201610200934A CN 105894436 A CN105894436 A CN 105894436A
- Authority
- CN
- China
- Prior art keywords
- image
- steganography
- distortion
- steganographic
- array
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005070 sampling Methods 0.000 title claims abstract description 16
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 230000011218 segmentation Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000003709 image segmentation Methods 0.000 claims description 2
- 238000005192 partition Methods 0.000 claims description 2
- 238000003491 array Methods 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000000682 scanning probe acoustic microscopy Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000007704 transition Effects 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
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
本发明公开了一种基于Gibbs抽样的图像隐写方法。首先,基于领域聚类的思想对隐写图像进行阵列分割,并对分割子阵列进行Gibbs抽样和STC编码;其次,充分利用子阵图像的方向信息(包括水平、垂直以及对角),采用加权求和的方式给出了隐写失真函数,维持了图像高阶统计特征。并利用Gibbs抽样理论实现信息的嵌入,从而在最小化嵌入失真的同时维持了隐写图像的特征模型。在Fisher准则的框架内优化函数参量,有效避免了图像特征集合过渡优化问题。实验显示,本文方法在降低图像模型畸变的同时,有效提升了隐写图像的分类误差,增强了隐写的安全性。
Description
技术领域
本发明涉及图像隐写方法,具体地,涉及一种基于Gibbs抽样的图像隐写方法。
背景技术
图像隐写算法的主要目的就是在以图像为载体,传输秘密信息,要求隐写过程中,畸变低、效率高、检测误差低。随着图像大数据信息技术的发展,图像隐写技术的需求越来越高,逐渐成为现代图像处理的研究热点。目前,针对图像隐写算法的研究主要集中在隐写畸变效应和载体图像模型保持两个方面。其中,隐写畸变效应的研究目的是降低隐写对载体图像造成的畸变干扰,STC(Syndrome-trellis code)编码技术实现了隐写与畸变效应的编码规则问题,但是,如何建立有效的失真度量函数,是该方面研究的重点。现有研究种失真函数的选择多数假设图像干扰为加性高斯白噪声,严格限制了失真函数优化与图像写隐分析的同步性,导致隐写嵌入过程中的畸变影响无法度量和优化。为克服这一缺陷,在充分考虑邻帧图像关联性的基础上,提出了Gibbs采样失真度量方法,初步探索了隐写过程中图像畸变特性度量及同步优化问题。但是,该方法只考虑相邻帧图像的关联问题,精度和范围都无法满足现代图像写隐技术的需求。针对载体图像模型保持的研究,目前多数都是针对图像某单一特征进行保持和优化,而忽略其他特征信息,虽然计算简单,但是安全性较差。其中比较具有代表性的是基于局部复杂度的自适应隐写方法,较好地保持了图像一阶直方图特征模型,但是图像高阶特征畸变严重。
发明内容
本发明的目的在于,针对上述问题,提出一种基于Gibbs抽样的图像隐写方法。该方法首先基于领域聚类的思想对隐写图像进行阵列分割,并对分割子阵列进行Gibbs抽样和STC(Syndrome-trellis code)编码;接着,充分利用子阵图像的方向信息(包括水平、垂直以及对角),采用加权求和的方式给出了隐写失真函数,维持了图像高阶统计特征。并利用Gibbs抽样理论实现信息的嵌入,从而在最小化嵌入失真的同时维持了隐写图像的特征模型。推导了失真函数与子阵列的统计相关性,在Fisher准则的框架内优化函数参量,有效避免了图像特征集合过渡优化问题。仿真结果显示,本文方法在降低图像模型畸变的同时,提升了嵌入的安全性。
为实现上述目的,本发明采用的技术方案是:
包括以下步骤:
步骤1、采用STC编码方法初始化隐写图像像素阵列;
步骤2、对初始化阵列进行分割,利用子阵图像的方向信息,包括
水平、垂直以及对角,采用加权求和的方式构造隐写失真函数;
步骤3、采用Gibbs抽样方法将需隐写信息嵌入阵列;
步骤4、利用上述计算结果更新阵列。
本发明的技术方案具有以下有益效果:
建立了图像隐写特征的集合模型,实现了隐写图像的多特征交叉实现,维持了较好的安全性、抗检测能力,并维持了图像高阶统计特征,从而保持较好的隐写安全性的同时,增强了模型的保持能力,优化了模型的畸变失真。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为具体实施例算法流程示意图
图2a为单幅隐写图像的输入效果,图2b为单幅隐写图像的输出效果。
图3a为水平方向的离散度加权和,3b为垂直方向的离散度加权和,3c为对角方向的离散度加权和。
图4为NIP与ITP特征比较图。
图5为SPAM与NIP特征比较图。
图6为SPAM与本发明特征集合比较图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
步骤1、采用STC编码方法初始化隐写图像像素阵列;
假定隐写目标图像x的特征空间参量I,则隐写目标图像的集合可以表示为x=(x1,x2,...,xn)∈In,嵌入信息后隐写图像构成的集合表示为Y,单个图像可以表示为Y,其取值范围满足Y的概率分布表示为在给定隐写秘密容量m时,隐写图像的最优分布需要满足满映射嵌入的条件下,失真度最小,表示为
隐写图像的分布特征进行STC(Syndrome-trellis code)编码后为
式(2)、(3)中,λ的取值主要取决于式(1)的约束条件。
其中为最小失真度,为Y的概率分布,D(y)为隐写图像的元素向量,s.t.H为约束函数,m为给定隐写秘密容量。
步骤2、对初始化阵列进行分割,利用子阵图像的方向信息,包括水平、垂直以及对角,采用加权求和的方式构造隐写失真函数;
为构造满足条件的失真函数,需要明确隐写图像的相邻关系,在给定相邻关系判定准则的基础上进行隐写图像的阵列分割。假设载体图像集合为κ={1,2,...,n},利用η={η(i)|i∈κ}表示划分邻域图像集合,如果η(i)满足以下三个条件:
则可以将η(i)中对应的图像元素表示为该图像的相邻元素,进而基于该准则进行图像子阵分割,在分辨率为n1×n2图像x中,可以将划分的子阵表示为
其中κst为s行,t列图像元素,分割过程中,为便于图像邻域的表示,给出了图像基团的定义。在隐写图像分割子集满足情况下,可以将c定义为同类子集组成的基团,在元素基团集合中,只有最大值的基团才能描述图像高阶统计特性,因此,在进行基团建模过程中,只需要采用极值基团建立相应的失真函数,将定义失真函数的基团定义为C。
隐写优化问题等效为图像在加性高斯噪声扰动情况下的失真最小求解问题。xi、xi为基团内两个元素,有
水平方向:
隐写信息以后,对载体图像的畸变表示为
其中,y为写入秘密后的图像,为写隐在水平方向的变化参量,D(y)为参量矩阵。同理,可以获取垂直方向和对角方向的畸变参量为 其中Hc(y)主要是反应相应方向基团变化的大小,则写隐图像全局失真函数可以表示为
C=C1∪C2∪C3∪C4 (8)
C=C1∪C2∪C3∪C4 (10)
ω度量失真畸变的取值为
基于(4)式可以将图像元素局部失真函数表示为
C1、C2、C3、C4为各个方向上失真函数的基团,ω为度量失真畸变的取值,ρi(y′iy~i)为图像元素局部失真函数。
Di,j的变化范围在正负半轴成对称状态,设置-T≤Di,j≤T,则在满足条件的区域内,可以将度量失真函数的参量θ表示为
θ=(θ-T,θ-t+1,...,θT-1,θT) (13)
权值ω计算为
步骤3、将吉布斯抽样信息嵌入原矩阵;
在图像集合采集服从分布的一个样本,通过吉布斯抽样,得出隐写图像集合Y上定义式(15)
其中,为yi转化为y′i的条件概率,为吉布斯抽样后隐写图像集合。在吉布斯抽样扫描中,针对访问对象σ进行全局扫描,计算
其中Iσ(i)为图像σ(i)的特征空间参量,为全局扫描后隐写图像集合。
步骤4、利用上述计算结果更新阵列。
进一步的,可在Fisher准则的框架内优化函数参量:
具体为采用具有较强抗检测性能的Fisher线性分类器对参量集合θ进行优化。
优化参量样本集合为U∈Rm×N,假设隐写图像通过阵列分割为两类子阵列,阵列1的元素数为N1,阵列2的元素个数为N2,整体优化隐写元素个数为N1+N2=N。两个子阵参数优化训练构成的子集分别为U1和U2。通过投影变换可以将m维矢量uj表示为标量vj,即vj=δTuj,uj∈U。Fisher分类器的准则函数可以表示为
为矢量矩阵通过投影变换到一维空间样本集合的均值,相应的映射样本离散度为
通过δ投影变换,利用拉格朗日变换可以获取极值点
δ*=(Sδ)-1(M1-M2) (19)
为隐写元素投影样本的均值向量。为计算的离散度总和。在给定元素集合的条件下,可以通过Fisher准则获取极值为
为确保隐写算法的安全性,需要展现较差的分类效果,失真函数的优化目标是为了取得较小的JFmax值。图1为含优化的本发明算法流程示意图。
为对本发明方法的优越性进行分析,以下针对理论分析进行了两类仿真实验。第一类是单幅隐写图像的离散度加权仿真;第二类是多幅图像隐写分类误差对比仿真。
离散值加权和越大,越集中,图像隐写的畸变影响就越小。实验结果如图2a、2b和图3a-3c所示。其中图2a、2b为单幅隐写图像的输入输出效果。图3a-3c为不同方向的离散度加权和,可以看出,单幅隐写图像的特征元素为45000*45000,通过图3a-3c的离散加权分布可以看出,本文方法的图像特征元素加权在三个方向上都保持了较高的集中度,鲁棒性较强。
针对多幅图像隐写分类误差计算,定义分类误差PE的计算表达式为
根据前面的分析可以知道,为了保持失真度小,需要PE具有较稳定的隐写特性;为了保持安全性,需要PE的值尽可能的大。针对仿真抽取的5000维特征,分别采用SPAM特征、NIP特征和ITF特征进行不同隐写像素的分类误差仿真分析,具体的结果如图4、图5、图6所示。从图4中可以看出,NIP特征优于ITF特征,能够保持较好的安全性,具有很好的模型保持特性,在隐写像素超过15以后,分类误差趋于稳定;从图5中可以看出,SPAM特征具有较优秀的隐写安全性,在像素超过15以后,分类误差趋于稳定;从图6中可以看出,本文方法明显由于SPAM特征,能够快速的趋于稳定,在像素超过4以后,就保持了平稳的收敛状态,并且维持在较高的误差范围,保持隐写安全性的同时,较好地保持了模型的失真畸变。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于Gibbs抽样的图像隐写方法,其特征在于,包括以下步骤:
步骤1、采用STC编码方法初始化隐写图像像素阵列;
步骤2、对初始化阵列进行分割,利用分割后子阵图像的方向信息,
采用加权求和的方式构造隐写失真函数,所述方向信息,包括水平、垂直以及对角;
步骤3、采用Gibbs抽样方法将需隐写信息嵌入阵列;
步骤4、利用上述计算结果更新阵列。
2.权利要求1所述的图像隐写方法,其特征在于,步骤2构造隐写失真函数过程为:
明确隐写图像的相邻关系,在相邻关系判定准则的基础上进行隐写图像的阵列分割:
假设载体图像集合为κ={1,2,...,n},利用η={η(i)|i∈κ}表示划分邻域图像集合,如果η(i)满足以下三个条件:
则可以将η(i)中对应的图像元素表示为该图像的相邻元素,进而基于该准则进行图像子阵分割;
将图像分割子集c定义为同类子集组成的基团,进行基团建模,采用极值基团建立相应的失真函数,将定义失真函数的基团定义为C;
隐写优化问题等效为图像在加性高斯噪声扰动情况下的失真最小求解问题;xi、xj为基团内两个元素,有水平方向:
隐写信息以后,对载体图像的畸变表示为
其中,y为写入秘密后的图像,为写隐在水平方向的变化参量,D(y)为参量矩阵;同理,可以获取垂直方向和对角方向的畸变参量为其中Hc(y)是表示相应方向基团变化的大小,则写隐图像全局失真函数可以表示为
C=C1∪C2∪C3∪C4
C=C1∪C2∪C3∪C4
图像元素局部失真函数为
其中,ω为度量失真畸变的取值,ρi(y′iy~i)为图像元素局部失真函数,C1、C2、C3为三个方向上失真函数的基团,Di,j的变化范围在正负半轴成对称状态,设置-T≤Di,j≤T,则在满足条件的区域内,可以将度量失真函数的参量θ表示为
θ=(θ-T,θ-t+1,...,θT-1,θT)。
3.权利要求2所述的图像隐写方法,其特征在于,步骤3采用Gibbs抽样方法将需隐写信息嵌入阵列过程为:
在图像集合采集服从分布的一个样本,通过吉布斯抽样,得出隐写图像集合Y上定义式
其中,为yi转化为y′i的条件概率,为吉布斯抽样后隐写图像集合;在吉布斯抽样扫描中,针对访问对象σ进行全局扫描,计算y′σ(i)∈Iσ(i),
其中Iσ(i)为图像σ(i)的特征空间参量,为全局扫描后隐写图像集合。
4.权利要求2所述的图像隐写方法,其特征在于,隐写失真函数通过多次优化获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610200934.9A CN105894436A (zh) | 2016-03-31 | 2016-03-31 | 一种基于Gibbs抽样的图像隐写方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610200934.9A CN105894436A (zh) | 2016-03-31 | 2016-03-31 | 一种基于Gibbs抽样的图像隐写方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105894436A true CN105894436A (zh) | 2016-08-24 |
Family
ID=57011846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610200934.9A Pending CN105894436A (zh) | 2016-03-31 | 2016-03-31 | 一种基于Gibbs抽样的图像隐写方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105894436A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107979711A (zh) * | 2017-11-28 | 2018-05-01 | 上海大学 | 基于最优化失真修改阻止隐蔽传输的方法 |
CN109644196A (zh) * | 2016-08-31 | 2019-04-16 | 迈克菲有限责任公司 | 消息保护 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006050364A2 (en) * | 2004-10-28 | 2006-05-11 | University Of Utah Research Foundation | Image denoising with information-theoretic, adaptive filtering |
CN103034853A (zh) * | 2013-01-08 | 2013-04-10 | 武汉大学 | 一种jpeg图像通用隐写分析方法 |
-
2016
- 2016-03-31 CN CN201610200934.9A patent/CN105894436A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006050364A2 (en) * | 2004-10-28 | 2006-05-11 | University Of Utah Research Foundation | Image denoising with information-theoretic, adaptive filtering |
CN103034853A (zh) * | 2013-01-08 | 2013-04-10 | 武汉大学 | 一种jpeg图像通用隐写分析方法 |
Non-Patent Citations (3)
Title |
---|
TOMÁS FILLER等: "Gibbs Construction in Steganography", 《IEEE TRANSACTIONS ON INFORMATION FORENSICS AND SECURITY》 * |
罗永有: "基于图像块间像素特征的隐写算法设计与实现", 《广西名族大学学报(自然科学版)》 * |
高瞻瞻等: "融合最小化失真与模型保持的隐写算法", 《中国图象图形学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109644196A (zh) * | 2016-08-31 | 2019-04-16 | 迈克菲有限责任公司 | 消息保护 |
CN109644196B (zh) * | 2016-08-31 | 2022-06-03 | 迈克菲有限责任公司 | 消息保护 |
CN107979711A (zh) * | 2017-11-28 | 2018-05-01 | 上海大学 | 基于最优化失真修改阻止隐蔽传输的方法 |
CN107979711B (zh) * | 2017-11-28 | 2019-07-23 | 上海大学 | 基于最优化失真修改阻止隐蔽传输的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zeng et al. | Large-scale JPEG image steganalysis using hybrid deep-learning framework | |
US9449395B2 (en) | Methods and systems for image matting and foreground estimation based on hierarchical graphs | |
CN106204415A (zh) | 一种新型的图像配准方法 | |
CN110830043B (zh) | 一种基于混合加权全变分和非局部低秩的图像压缩感知重构方法 | |
CN111160400A (zh) | 一种基于修正边界攻击的对抗攻击方法 | |
CN101866490A (zh) | 基于差分免疫克隆聚类的图像分割方法 | |
CN111652240A (zh) | 一种基于cnn的图像局部特征检测与描述方法 | |
CN115311502A (zh) | 基于多尺度双流架构的遥感图像小样本场景分类方法 | |
Xiang et al. | Embedding new data points for manifold learning via coordinate propagation | |
Li et al. | Constructing immunized stego-image for secure steganography via artificial immune system | |
Yang et al. | A Novel Universal Steganalysis Algorithm Based on the IQM and the SRM. | |
CN105894436A (zh) | 一种基于Gibbs抽样的图像隐写方法 | |
Tang et al. | Reinforcement learning of non-additive joint steganographic embedding costs with attention mechanism | |
CN112966748B (zh) | 基于边缘感知双分支fcn的极化sar图像分类方法 | |
Zheng et al. | Extracting non-negative basis images using pixel dispersion penalty | |
CN112529050A (zh) | 一种广义弹性网络非平行支持向量机分类算法 | |
CN105894437A (zh) | 一种基于马尔科夫链的隐写算法 | |
Bacciu et al. | Compositional generative mapping of structured data | |
Calderhead et al. | Sparse approximate manifolds for differential geometric mcmc | |
Meng et al. | Passage method for nonlinear dimensionality reduction of data on multi-cluster manifolds | |
CN113344049B (zh) | 一种基于Sinkhorn距离的盲高光谱解混模型的构建方法 | |
CN114155443B (zh) | 一种基于多感受野图注意力网络的高光谱图像分类方法 | |
CN113610125B (zh) | 一种基于编码器-解码器图神经网络的高光谱分类方法 | |
Wang et al. | Moiré pattern removal with a generative adversarial network | |
Li et al. | Geometric structure based intelligent collaborative compressive sensing for image reconstruction by l0 minimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160824 |