CN112862656A - 一种抗打印扫描拍摄的数字图像水印方法 - Google Patents
一种抗打印扫描拍摄的数字图像水印方法 Download PDFInfo
- Publication number
- CN112862656A CN112862656A CN202110271615.8A CN202110271615A CN112862656A CN 112862656 A CN112862656 A CN 112862656A CN 202110271615 A CN202110271615 A CN 202110271615A CN 112862656 A CN112862656 A CN 112862656A
- Authority
- CN
- China
- Prior art keywords
- image
- watermark
- zernike
- embedding
- moment
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 128
- 238000000605 extraction Methods 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000013139 quantization Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 9
- 230000000877 morphologic effect Effects 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 230000021615 conjugation Effects 0.000 claims description 5
- 230000006835 compression Effects 0.000 abstract description 3
- 238000007906 compression Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000002411 adverse Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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
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
本发明提出一种抗打印扫描拍摄的数字图像水印方法,解决了如何抵抗几何变换、打印扫描及拍摄对数字图像水印的不利影响,兼顾水印提取与水印图像主观质量的问题,数字图像水印嵌入过程基于原始图像的Zernike矩,通过量化嵌入的方法进行水印嵌入,根据水印嵌入前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像,提升了水印图像的主观质量,基于Zernike矩变换,设计了一种匹配的数字图像水印提取算法,其中包含图像自动定位方法,提取水印过程可以抵抗旋转、缩放、噪声、压缩等带来的影响,从而保证了水印信息的有效提取。
Description
技术领域
本发明涉及数字图像认证的技术领域,更具体地,涉及一种抗打印扫描拍摄的数字图像水印方法。
背景技术
数字水印是将特定的信息嵌入数字信号中,数字信号可能是音频、图像或视频等,数字水印可分为浮现式和隐藏式两种,前者是可被看见的水印,其所包含的信息可在观看图片或视频时同时被看见;后者是以数字资料的方式加入音频、图片或视频中,一般情况下无法被看见,本发明属于隐藏式水印,其重要应用之一是保护著作权,期望能借此避免或阻止数字媒体未经授权的复制和拷贝。
数字水印图像在使用过程中可能遭受到几何变换或者打印扫描/拍摄等攻击,这些攻击能够破坏水印检测的同步性。微小的几何变换攻击就可能使得水印检测失败,打印扫描攻击涉及到数模转换和模数转换过程,其中包含多次不均匀采样和量化,经过打印扫描视觉效果和原图相似,但图像质量已被破坏,出现像素失真和几何失真的现象;经过打印拍摄,不仅存在像素失真,还会出现倾斜拍摄导致的非仿射畸变。因此,设计出能够稳健地抵抗这些攻击的水印方法十分有必要。
有些人提出了许多稳健的数字水印方法来保护数字水印。
2007年,Y.Xin,S.Liao及M.Pawlak等人在Pattern Recognition上发表“Geometrically robust image watermarking on a circular domain”的文章(vol.40,no.1,pp.3740-3752),利用Zernike矩的不变性实现了对几何变换和有损压缩的鲁棒性;近年来,一种可以恢复载体图像的可逆水印算法也得到广泛研究,2020年,R.Hu及S.Xiang等人在IEEE Transactions on Image Processing发表了“Cover-Lossless Robust ImageWatermarking Against Geometric Deformations”的文章(vol.30,pp.318-331),提出在低阶归一化Zernike幅值中嵌入鲁棒水印,并将鲁棒水印嵌入带来的失真作为补偿信息进行可逆隐藏,实现载体无损恢复,2020年6月5日,中国发明专利(CN111242831A)中公开了一种基于Zernike矩抗几何攻击的可逆鲁棒水印方法,首先计算原始图像的Zernike矩,采用量化水印方法进行水印嵌入,处理Zernike正反变换过程中存在的失真信息并判断图像是否受到攻击;图像受到攻击时,计算受到攻击后带有水印信息图像的Zernike矩,进行量化水印提取得到水印信息,该专利中提出的方法考虑了利用Zernike矩帮助图像水印过程抵抗旋转、缩放、噪声、压缩等带来的影响,在未受到攻击时可以提取水印和恢复原始图像,在受到攻击时可以有效提取出水印。
发明内容
为解决如何抵抗几何变换、打印扫描及拍摄对数字图像水印的不利影响,兼顾水印提取与水印图像主观质量的问题,本发明提出一种抗打印扫描拍摄的数字图像水印方法,提升水印图像的主观质量,保证水印信息比特流的有效提取。
为了达到上述技术效果,本发明的技术方案如下:
一种抗打印扫描拍摄的数字图像水印方法,所述方法包括数字图像水印嵌入过程及数字图像水印提取过程,所述数字图像水印嵌入过程首先确定待嵌入水印的原始图像,基于原始图像的Zernike矩,通过量化嵌入的方法进行水印嵌入,根据水印嵌入前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像;所述数字图像水印提取过程配合数字图像水印嵌入过程,将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,基于Zernike矩变换,进行水印提取。
优选地,所述数字图像水印嵌入过程具体包括:
S1.确定原始图像的种类,根据原始图像的种类选择水印嵌入通道;
S2.确定原始图像Zernike矩的最大阶数N,计算单通道图像的n阶m重Zernike矩Anm;
S3.筛选Zernike矩,将筛选后的Zernike矩组成可嵌入水印的矩集合P,固定Anm,生成矩集合P中每个Zernike矩Anm的重构图像;
S4.计算重构图像的区域指标,根据指标大小确定待嵌入水印信息的Zernike矩位置;
S5.对Zernike矩进行归一化,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,嵌入水印,生成嵌入水印后的归一化Zernike矩幅值,并计算嵌入水印后的Zernike矩;
S6.根据嵌入水印前后的Zernike矩重构图像,计算两幅重构图像的差值并进行选择性区域截断,将优化后的差值图像叠加到原始图像上生成水印图像。
优选地,步骤S1所述确定原始图像的种类,根据原始图像的种类选择水印嵌入通道的过程包括:
确认原始图像是彩色图像还是灰度图像,若是彩色图像,则将原始图像变换到YCbCr色度空间,选择Y通道进行水印嵌入;若是灰度图像,则直接在原始图像上进行水印嵌入。
优选地,步骤S2所述n阶m重Zernike矩Anm的计算公式为:
其中,x,y表示原始图像的像素坐标,n为非负整数,表示阶数;m为整数,表示重复度;n和m满足:0≤n≤N,n-|m|为偶数,n≥|m|;Vnm(x,y)表示Zernike多项式,是单位圆内一类完备的正交基,表示为:
其中,Rnm(ρ)为径向多项式,定义为:
在此,在计算数字图像f(x,y)的Zernike矩时,取图像的中心为原点,把所有像素坐标映射到单位圆内,即x2+y2≤1,单位圆外的点不参加计算。至此,求得f(x,y)从0到N的各阶Zernike矩。
优选地,步骤S3所述筛选Zernike矩的过程包括:
根据Zernike矩的共轭性和对称性,去除重复的矩;
根据部分Zernike矩的正交性偏离问题,去除重复度为4的倍数的矩;
筛选后的Zernike矩组成可嵌入水印的矩集合P表示为:
P={Anm,n≤N,m≥0,m≠4i};
x,y均表示图像的像素坐标;
步骤S4所述区域指标为每个重构图像内切圆区域与圆环区域的像素指标,计算的具体过程包括:
计算区域内像素值的均值平方与方差之和Gq:
计算矩集合P中所有Zernike矩在圆环区域的值的平均值保留的Zernike矩,构成集合P2;将集合P2中Zernike矩根据其在内切圆区域的值从小到大排序,选取对应最小的L个Zernike矩作为嵌入水印的位置,筛选后的L个Zernike矩组成Zernike矩向量L表示待嵌入水印信息的比特长度,考虑了对重构图像主观质量的影响。
优选地,步骤S5中所述对Zernike矩进行归一化的过程为:
利用Zernike零阶矩A00对Apq中的Zernike矩进行归一化,公式为:
利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值的过程为:
设待嵌入的水印信息序列表示为w={wi,i=1,L,36,wi∈{0,1}},对于水印信息序列中的每个水印比特,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,生成嵌入水印后的归一化Zernike矩幅值的公式为:
优选地,步骤S6所述的根据嵌入水印前后的Zernike矩重构图像的过程为:
利用嵌入水印前的Zernike矩重构得到图像Ilow,表达式为:
对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像的过程包括:
计算嵌入水印前后重构图像的差值图像,表达式为:
其中,Irw表示差值图像;设差值图像Irw圆环区域q2的每个像素点表示为(xi,yi),若某一个像素点在原始图像中位于平滑区域,则对其在差值图像Irw中的像素值vi进行截断优化,优化公式为:
其中,I表示原始图像。
在此,若原始图像为灰度图,则IW为最终含水印的灰度图;若原始图像为彩色图像,则用含水印的单通道图像IW代替原始图像YCbCr色度空间中的Y通道并变换到RGB色度空间,得到最终含水印的彩色图像。
优选地,数字图像水印提取过程具体包括:
SA.将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,对待测图像进行预处理得到待测对象;
SB.采用形态学方法对待测对象中目标水印图像区域的顶点进行自动定位,再经过透视变换进行图像矫正,并裁剪出目标水印图像,进而计算图像边缘的行列像素均值,剔除无关背景白边;
SC.对待测图像进行与水印嵌入过程相同的Zernike矩计算,并确定嵌入水印信息的Zernike矩位置,组成水印提取阶段筛选后的Zernike矩向量;
SD.对Zernike矩进行归一化操作,利用对应的量化提取方法对Zernike矩幅值进行水印提取,得到水印比特信息。
在此,所述数字图像水印提取过程是与数字图像水印嵌入过程配合的,在步骤SC中确定嵌入水印信息的Zernike矩位置的做法与前述水印嵌入过程相同,通过自动定位图像顶点、透视矫正、白边剔除等步骤可以实现目标水印区域的精准裁剪,进而准确提取有意义的水印信息比特流。
优选地,步骤SA所述对待测图像进行预处理的过程包括:
确定水印图像的种类,所述种类包括灰度图像及彩色图像;当水印图像为灰度图像时,将经扫描的水印图像的灰度通道作为待测对象,当水印图像为彩色图像时,将经扫描的水印图像YCbCr空间的Y通道作为待测对象;
将经打印拍摄的水印图像的YCbCr空间的Y通道作为待测对象;
若水印图像经几何攻击,则不进行步骤SB所述的过程;
利用Zernike零阶矩A′00对A′pq中的Zernike矩进行归一化,公式为:
其中,T=10κ,κ∈Ν*,与水印嵌入过程取值相同;
利用对应的量化提取方法对Zernike矩幅值进行水印提取,得到水印比特信息的过程满足:
根据每个归一化后的Zernike矩的幅值提取水印比特信息,公式为:
其中,L表示Zernike矩位置的个数。
与现有技术相比,本发明技术方案的有益效果是:
本发明提出一种抗打印扫描拍摄的数字图像水印方法,包括数字图像水印嵌入过程及数字图像水印提取过程,数字图像水印嵌入过程基于原始图像的Zernike矩,通过量化嵌入的方法进行水印嵌入,根据水印嵌入前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像,提升了水印图像的主观质量,基于Zernike矩变换,通过数字图像水印提取过程配合数字图像水印嵌入过程,由于Zernike矩具有旋转不变性、对噪声不敏感的特点,且归一化后的Zernike矩具有缩放不变性,提取水印过程可以抵抗旋转、缩放、噪声、压缩等带来的影响,从而保证了水印信息的有效提取。
附图说明
图1表示本发明实施例中提出的抗打印扫描拍摄的数字水印图像方法的整体示意流程框图;
图2表示本发明实施例中提出的Zernike矩重构图像的像素值分析区域中图像内切圆区域示意图;
图3表示本发明实施例中提出的Zernike矩重构图像的像素值分析区域中图像内圆环区域示意图;
图4表示本发明实施例中提出的实验图像的示意图;
图5表示本发明实施例中提出的水印图像经打印扫描后的示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好地说明本实施例,附图某些部位会有省略、放大或缩小,并不代表实际尺寸;
对于本领域技术人员来说,附图中某些公知内容说明可能省略是可以理解的。
附图中描述位置关系的用于仅用于示例性说明,不能理解为对本专利的限制;
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示的抗打印扫描拍摄的数字水印图像方法的整体示意流程框图,参见图1,所述方法包括数字图像水印嵌入过程及数字图像水印提取过程,所述数字图像水印嵌入过程首先确定待嵌入水印的原始图像,基于原始图像的Zernike矩,通过量化嵌入的方法进行水印嵌入,根据水印嵌入前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像;所述数字图像水印提取过程配合数字图像水印嵌入过程,将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,基于Zernike矩变换,进行水印提取。
其中,数字图像水印嵌入过程具体包括:
S1.确定原始图像的种类,根据原始图像的种类选择水印嵌入通道;
步骤S1所述确定原始图像的种类,根据原始图像的种类选择水印嵌入通道的过程包括:
确认原始图像是彩色图像还是灰度图像,若是彩色图像,则将原始图像变换到YCbCr色度空间,选择Y通道进行水印嵌入;若是灰度图像,则直接在原始图像上进行水印嵌入。
S2.确定原始图像Zernike矩的最大阶数N,计算单通道图像的n阶m重Zernike矩Anm;
n阶m重Zernike矩Anm的计算公式为:
其中,x,y表示原始图像的像素坐标,n为非负整数,表示阶数;m为整数,表示重复度;n和m满足:0≤n≤N,n-|m|为偶数,n≥|m|;Vnm(x,y)表示Zernike多项式,是单位圆内一类完备的正交基,表示为:
其中,Rnm(ρ)为径向多项式,定义为:
在计算数字图像f(x,y)的Zernike矩时,取图像的中心为原点,把所有像素坐标映射到单位圆内,即x2+y2≤1,单位圆外的点不参加计算。至此,求得f(x,y)从0到N的各阶Zernike矩。
S3.筛选Zernike矩,将筛选后的Zernike矩组成可嵌入水印的矩集合P,固定Anm,生成矩集合P中每个Zernike矩Anm的重构图像;筛选Zernike矩的过程包括:
根据Zernike矩的共轭性和对称性,去除重复的矩;
根据部分Zernike矩的正交性偏离问题,去除重复度为4的倍数的矩;
筛选后的Zernike矩组成可嵌入水印的矩集合P表示为:
P={Anm,n≤N,m≥0,m≠4i};
x,y均表示图像的像素坐标。
S4.计算重构图像的区域指标,根据指标大小确定待嵌入水印信息的Zernike矩位置;
步骤S4所述区域指标为每个重构图像内切圆区域与圆环区域的像素指标,图2表示Zernike矩重构图像的像素值分析区域中图像内切圆区域,图3表示Zernike矩重构图像的像素值分析区域中图像内圆环区域,计算每个重构图像内切圆区域与圆环区域的像素指标,具体过程包括:
计算区域内像素值的均值平方与方差之和Gq:
计算矩集合P中所有Zernike矩在圆环区域的值的平均值保留的Zernike矩,构成集合P2;将集合P2中Zernike矩根据其在内切圆区域的值从小到大排序,选取对应最小的L个Zernike矩作为嵌入水印的位置,筛选后的L个Zernike矩组成Zernike矩向量L表示待嵌入水印信息的比特长度,考虑了对重构图像主观质量的影响。
S5.对Zernike矩进行归一化,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,嵌入水印,生成嵌入水印后的归一化Zernike矩幅值,并计算嵌入水印后的Zernike矩;
利用Zernike零阶矩A00对Apq中的Zernike矩进行归一化,公式为:
利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值的过程为:
设待嵌入的水印信息序列表示为w={wi,i=1,L,36,wi∈{0,1}},对于水印信息序列中的每个水印比特,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,生成嵌入水印后的归一化Zernike矩幅值的公式为:
S6.根据嵌入水印前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断优化的方式,叠加到原始图像上生成水印图像。
根据嵌入水印前后的Zernike矩重构图像的过程为:
利用嵌入水印前的Zernike矩重构得到图像Ilow,表达式为:
对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像的过程包括:
计算嵌入水印前后重构图像的差值图像,表达式为:
其中,Irw表示差值图像;设差值图像Irw圆环区域q2的每个像素点表示为(xi,yi),若某一个像素点在原始图像中位于平滑区域,则对其在差值图像Irw中的像素值vi进行截断优化,优化公式为:
其中,I表示原始图像。若原始图像为灰度图,则IW为最终含水印的灰度图;若原始图像为彩色图像,则用含水印的单通道图像IW代替原始图像YCbCr色度空间中的Y通道并变换到RGB色度空间,得到最终含水印的彩色图像。
在本实施例中,数字图像水印提取过程具体包括:
SA.将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,对待测图像进行预处理得到待测对象;
SB.采用形态学方法对待测对象中目标水印图像区域的顶点进行自动定位,再经过透视变换进行图像矫正,并裁剪出目标水印图像,进而计算图像边缘的行列像素均值,剔除无关背景白边;
具体的,对待测对象应用高斯平滑滤波器后使用Sobel算子进行边缘检测;对边缘点图像进行形态学闭运算平滑轮廓、填充缝隙,进行形态学开运算删除小面积连通分量,再跟踪图像中连通对象的外边界,选择最大连通区域作为目标轮廓;计算包含该目标轮廓的最小凸多边形,得到其顶点集合,根据两点确定一条直线,得到顶点集合对应的线段集合;剔除线段集合中相邻冗余线段,确定水印图像边缘线段,进一步得到图像的顶点,使用透视变换进行图像矫正并裁剪出变换后的水印图像Itrans。
剔除水印图像Itrans边缘的白边。对图像Itrans边缘的行、列分别进行单行、单列像素均值计算,根据像素阈值剔除无关背景白边,确定图像边缘,得到最终的待测目标图像Iextract。
SC.对待测图像进行与水印嵌入过程相同的Zernike矩计算,并确定嵌入水印信息的Zernike矩位置,组成水印提取阶段筛选后的Zernike矩向量;
SD.对Zernike矩进行归一化操作,利用对应的量化提取方法对Zernike矩幅值进行水印提取,得到水印比特信息。
在此,所述数字图像水印提取过程是与数字图像水印嵌入过程配合的,在步骤SC中确定嵌入水印信息的Zernike矩位置的做法与前述水印嵌入过程相同,,通过自动定位图像顶点、透视矫正、白边剔除等步骤可以实现目标水印区域的精准裁剪,进而准确提取有意义的水印信息比特流。
在本实施例中,步骤SA所述对待测图像进行预处理的过程包括:
确定水印图像的种类,所述种类包括灰度图像及彩色图像;当水印图像为灰度图像时,将经扫描的水印图像的灰度通道作为待测对象,当水印图像为彩色图像时,将经扫描的水印图像YCbCr空间的Y通道作为待测对象;
将经打印拍摄的水印图像的YCbCr空间的Y通道作为待测对象;
若水印图像经几何攻击,则不进行步骤SB所述的过程;
利用Zernike零阶矩A′00对A′pq中的Zernike矩进行归一化,公式为:
其中,T=10κ,κ∈Ν*,与水印嵌入过程取值相同;
利用对应的量化提取方法对Zernike矩幅值进行水印提取,得到水印比特信息的过程满足:
根据每个归一化后的Zernike矩的幅值提取水印比特信息,公式为:
其中,L表示Zernike矩位置的个数。
下面结合具体实验对象进一步说明本发明所提方法的实施过程,如图4所示,左侧表示原始图像Lena,是512×512的灰度图像,嵌入36比特的水印信息,具体做法为:
一、确定图像Zernike矩的最大阶数N=31,对大小为512×512的Lena原始图像f(x,y),计算其n阶m重Zernike矩:
其中,n为非负整数,称为阶数;m为整数,称为重复度。n和m满足:
0≤n≤N,n-|m|为偶数,n≥|m|。
其中Zernike多项式Vnm(x,y)为:
Vnm(x,y)=Vnm(ρ,θ)=Rnm(ρ)ejmθ
径向多项式Rnm(ρ)为:
在计算数字图像f(x,y)的Zernike矩时,取图像的中心为原点,把所有像素坐标映射到单位圆内,即x2+y2≤1,单位圆外的点不参加计算。至此,求得f(x,y)从0到N的各阶Zernike矩。
二、确定待嵌入水印信息的Zernike矩位置。根据Zernike矩的共轭性和对称性,去除重复的矩;根据重复度为4×i(i=0,1,2,L)的Zernike矩存在正交性偏离的问题,去除此类不准确的矩,不作为嵌入水印的位置。经上述筛选后的矩组成可嵌入水印的矩集合P,即P={Anm,n≤N,m≥0,m≠4i}。
令集合P中Anm值相同:Anm=10+10i,计算集合P中每个Zernike矩Anm的重构图像:
计算每个重构图像内切圆区域与圆环区域(如图2所示)的像素指标——分析区域内像素值的均值平方与方差之和Gq:
计算集合P中所有Zernike矩在圆环区域的值的平均值保留的Zernike矩,构成集合P2;将集合P2中Zernike矩根据其在内切圆区域的值从小到大排序,考虑对重构图像主观质量的影响,选取对应最小的36个Zernike矩作为嵌入水印的位置,筛选后的36个Zernike矩组成Zernike矩向量
三.对选定的Zernike矩进行归一化操作后使用量化嵌入方法对其幅值进行水印嵌入,具体步骤如下:
使用Zernike零阶矩A00对Apq中的Zernike矩进行归一化:
其中,T是归一化常数,令T=1000。
准备待嵌入的水印信息序列,令其为w={wi,i=1,L,36,wi∈{0,1}},如“101000010010110011111010000100101100”。对于每个水印比特,采用量化水印方法修改对应位置的Zernike矩幅值来实现嵌入水印的目的,最终生成嵌入水印后的归一化幅值
计算嵌入水印后的Zernike矩,表达式:
四.通过嵌入水印前后的Zernike矩重构空域图像,使用选择性区域截断差值图像的方式达到优化水印图像主观质量的目的,并生成最终水印图像IW。具体步骤如下:
利用未修改的Zernike矩重构得到图像Ilow:
如图5所示,将水印图像经过打印扫描处理,作为待测图像,进行水印抗打印扫描攻击检测。
在打印扫描场景下,水印图像是灰度图,因此将打印图像的灰度通道作为待测对象。
在打印扫描场景下,自动定位水印图像区域,进行图像矫正。对待测对象应用高斯平滑滤波器后使用Sobel算子进行边缘检测;对边缘点图像进行形态学闭运算平滑轮廓、填充缝隙,进行形态学开运算删除小面积连通分量,再跟踪图像中连通对象的外边界,选择最大连通区域作为目标轮廓;计算包含该目标轮廓的最小凸多边形,得到其顶点集合,根据两点确定一条直线,得到顶点集合对应的线段集合;剔除线段集合中相邻冗余线段,确定水印图像边缘线段,进一步得到图像的顶点,使用透视变换进行图像矫正并裁剪出变换后的水印图像Itrans。剔除水印图像Itrans边缘的白边。对Itrans图像边缘的行、列分别进行单行、单列像素均值计算,根据像素阈值剔除无关背景白边,确定图像边缘,得到最终的待测目标图像Iextract。
对选定的Zernike矩进行归一化操作后使用对应的量化提取方法对其幅值进行水印提取,具体步骤如下:
使用Zernike零阶矩A′00对A′pq中的Zernike矩进行归一化:
其中,T=1000,与嵌入过程相同取值;
根据每个归一化后的Zernike矩幅值提取水印比特:
其中,S=32,与嵌入过程相同取值。
附图中描述位置关系的用于仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种抗打印扫描拍摄的数字图像水印方法,其特征在于,所述方法包括数字图像水印嵌入过程及数字图像水印提取过程,所述数字图像水印嵌入过程首先确定待嵌入水印的原始图像,基于原始图像的Zernike矩,通过量化嵌入的方法进行水印嵌入,根据水印嵌入前后的Zernike矩重构图像,对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像;所述数字图像水印提取过程配合数字图像水印嵌入过程,将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,基于Zernike矩变换,进行水印提取。
2.根据权利要求1所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,所述数字图像水印嵌入过程具体包括:
S1.确定原始图像的种类,根据原始图像的种类选择水印嵌入通道;
S2.确定原始图像Zernike矩的最大阶数N,计算单通道图像的n阶m重Zernike矩Anm;
S3.筛选Zernike矩,将筛选后的Zernike矩组成可嵌入水印的矩集合P,固定Anm,生成矩集合P中每个Zernike矩Anm的重构图像;
S4.计算重构图像的区域指标,根据指标大小确定待嵌入水印信息的Zernike矩位置;
S5.对Zernike矩进行归一化,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,嵌入水印,生成嵌入水印后的归一化Zernike矩幅值,并计算嵌入水印后的Zernike矩;
S6.根据嵌入水印前后的Zernike矩重构图像,计算两幅重构图像的差值并进行选择性区域截断,将优化后的差值图像叠加到原始图像上生成水印图像。
3.根据权利要求2所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,步骤S1所述确定原始图像的种类,根据原始图像的种类选择水印嵌入通道的过程包括:
确认原始图像是彩色图像还是灰度图像,若是彩色图像,则将原始图像变换到YCbCr色度空间,选择Y通道进行水印嵌入;若是灰度图像,则直接在原始图像上进行水印嵌入。
5.根据权利要求4所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,步骤S3所述筛选Zernike矩的过程包括:
根据Zernike矩的共轭性和对称性,去除重复的矩;
根据部分Zernike矩的正交性偏离问题,去除重复度为4的倍数的矩;
筛选后的Zernike矩组成可嵌入水印的矩集合P表示为:
P={Anm,n≤N,m≥0,m≠4i};
x,y均表示图像的像素坐标;
步骤S4所述区域指标为每个重构图像内切圆区域与圆环区域的像素指标,计算的具体过程包括:
计算区域内像素值的均值平方与方差之和Gq:
6.根据权利要求5所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,步骤S5中所述对Zernike矩进行归一化的过程为:
利用Zernike零阶矩A00对Apq中的Zernike矩进行归一化,公式为:
利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值的过程为:
设待嵌入的水印信息序列表示为w={wi,i=1,L,36,wi∈{0,1}},对于水印信息序列中的每个水印比特,利用量化嵌入方法修改对应Zernike矩位置的Zernike矩幅值,生成嵌入水印后的归一化Zernike矩幅值的公式为:
7.根据权利要求6所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,步骤S6所述的根据嵌入水印前后的Zernike矩重构图像的过程为:
利用嵌入水印前的Zernike矩重构得到图像Ilow,表达式为:
对重构图像的差值采用选择性区域截断的方式,叠加到原始图像上生成水印图像的过程包括:
计算嵌入水印前后重构图像的差值图像,表达式为:
其中,Irw表示差值图像;设差值图像Irw圆环区域q2的每个像素点表示为(xi,yi),若某一个像素点在原始图像中位于平滑区域,则对其在差值图像Irw中的像素值vi进行截断优化,优化公式为:
其中,I表示原始图像。
8.根据权利要求7所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,数字图像水印提取过程具体包括:
SA.将数字图像水印嵌入过程生成的水印图像经打印扫描/拍摄后作为待测图像,对待测图像进行预处理得到待测对象;
SB.采用形态学方法对待测对象中目标水印图像区域的顶点进行自动定位,再经过透视变换进行图像矫正,并裁剪出目标水印图像,进而计算图像边缘的行列像素均值,剔除无关背景白边;
SC.对待测图像进行与水印嵌入过程相同的Zernike矩计算,并确定嵌入水印信息的Zernike矩位置,组成水印提取阶段筛选后的Zernike矩向量;
SD.对Zernike矩进行归一化操作,利用对应的量化提取方法对Zernike矩幅值进行水印提取,得到水印比特信息。
9.根据权利要求8所述的抗打印扫描拍摄的数字图像水印方法,其特征在于,步骤SA所述对待测图像进行预处理的过程包括:
确定水印图像的种类,所述种类包括灰度图像及彩色图像;当水印图像为灰度图像时,将经扫描的水印图像的灰度通道作为待测对象,当水印图像为彩色图像时,将经扫描的水印图像YCbCr空间的Y通道作为待测对象;
将经打印拍摄的水印图像的YCbCr空间的Y通道作为待测对象;
若水印图像经几何攻击,则不进行步骤SB所述的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110271615.8A CN112862656B (zh) | 2021-03-12 | 2021-03-12 | 一种抗打印扫描拍摄的数字图像水印方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110271615.8A CN112862656B (zh) | 2021-03-12 | 2021-03-12 | 一种抗打印扫描拍摄的数字图像水印方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112862656A true CN112862656A (zh) | 2021-05-28 |
CN112862656B CN112862656B (zh) | 2024-02-02 |
Family
ID=75994382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110271615.8A Active CN112862656B (zh) | 2021-03-12 | 2021-03-12 | 一种抗打印扫描拍摄的数字图像水印方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112862656B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023000991A1 (zh) * | 2021-07-20 | 2023-01-26 | 北京沃东天骏信息技术有限公司 | 针对水印载体的水印处理方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903075A (zh) * | 2012-10-15 | 2013-01-30 | 西安电子科技大学 | 基于图像特征点全局校正的鲁棒水印方法 |
WO2016114950A1 (en) * | 2015-01-15 | 2016-07-21 | Gopro, Inc. | Watermarking digital images to increase bit dept |
CN110097490A (zh) * | 2019-04-26 | 2019-08-06 | 上海理工大学 | 基于IWT-Schur的感兴趣区域可逆水印嵌入及提取方法 |
CN111242831A (zh) * | 2020-01-20 | 2020-06-05 | 暨南大学 | 一种基于Zernike矩抗几何攻击的可逆鲁棒水印方法 |
-
2021
- 2021-03-12 CN CN202110271615.8A patent/CN112862656B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903075A (zh) * | 2012-10-15 | 2013-01-30 | 西安电子科技大学 | 基于图像特征点全局校正的鲁棒水印方法 |
WO2016114950A1 (en) * | 2015-01-15 | 2016-07-21 | Gopro, Inc. | Watermarking digital images to increase bit dept |
CN110097490A (zh) * | 2019-04-26 | 2019-08-06 | 上海理工大学 | 基于IWT-Schur的感兴趣区域可逆水印嵌入及提取方法 |
CN111242831A (zh) * | 2020-01-20 | 2020-06-05 | 暨南大学 | 一种基于Zernike矩抗几何攻击的可逆鲁棒水印方法 |
Non-Patent Citations (2)
Title |
---|
叶天语 等: "抗亮度和对比度调整的盲鲁棒量化水印算法", 光子学报, vol. 40, no. 10, pages 1577 - 1584 * |
沃焱 等: "基于修改的Zernike矩的抗几何攻击的数字水印方法", 计算机科学, vol. 36, no. 5, pages 247 - 250 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023000991A1 (zh) * | 2021-07-20 | 2023-01-26 | 北京沃东天骏信息技术有限公司 | 针对水印载体的水印处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112862656B (zh) | 2024-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5566811B2 (ja) | プリントアンドスキャン文書の画像評価のためのボケ除去および監視適応的スレッショルディング | |
Qureshi et al. | A bibliography of pixel-based blind image forgery detection techniques | |
Sencar et al. | Overview of state-of-the-art in digital image forensics | |
US7068812B2 (en) | Decoding hidden data from imagery | |
EP1485863B1 (en) | Authenticating printed objects using digital watermarks associated with multidimensional quality metrics | |
US9349237B2 (en) | Method of authenticating a printed document | |
US7965894B2 (en) | Method for detecting alterations in printed document using image comparison analyses | |
EP2014082A1 (en) | Generating a bitonal image from a scanned colour image | |
JP4362538B2 (ja) | 画像処理装置、画像形成装置、画像送信装置、画像読取装置、画像処理システム、画像処理方法、画像処理プログラムおよびその記録媒体 | |
JP5534411B2 (ja) | 画像処理装置 | |
Gupta et al. | A study on source device attribution using still images | |
Ebrahimnejad et al. | Adaptive Removal of high-density salt-and-pepper noise (ARSPN) for robust ROI detection used in watermarking of MRI images of the brain | |
JP4595014B2 (ja) | 電子透かし埋め込み装置および検出装置 | |
CN112862656A (zh) | 一种抗打印扫描拍摄的数字图像水印方法 | |
CN109829890B (zh) | 一种jpeg图像载体安全性评价方法 | |
CN116757904A (zh) | 基于数字图像水印的信息隐藏和验证方法及装置 | |
JP6006675B2 (ja) | マーカ検出装置、マーカ検出方法、及びプログラム | |
JP2006107034A (ja) | 画像解析方法及び装置、画像処理方法及びシステム、これらの動作プログラム | |
JP4310031B2 (ja) | 画像処理装置、画像処理方法及び記憶媒体 | |
CN117635411B (zh) | 基于混合域分解技术的数字水印处理方法 | |
Kumar et al. | A review of image features extraction techniques and their applications in image forensic | |
CN111325650A (zh) | 一种可逆水印方法 | |
CN111861845B (zh) | 一种基于阈值分割和直方图均衡的可逆水印方法 | |
Singh et al. | Digital Image Forensic Techniques with Attacks: A Review | |
Wang et al. | Improving scanned binary image watermarking based on additive model and sampling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |