本文主要是介绍数字水印 | 图像标准化论文:Digital Watermarking Robust to Geometric Distortions(二),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
目录
- C 变换参数的确定
- D 水印的影响
- E 可替代的标准化过程
🤖原文: Digital Watermarking Robust to Geometric Distortions
🤖前言: 这是一篇 2005 年的 SCI 一区 + CCF-A,但是网上关于它的讲解貌似挺少的。文中提出了两种数字水印方案,但是我只关注第一种方案中的图像标准化技术。由于本人很菜,因此可能存在翻译或者理解的错误,请各位指正!
C 变换参数的确定
在本节中,我们展示了如何确定与变换相关的参数,使它们达到各自的标准化目标。
矩阵 A x = ( 1 β 0 1 ) \mathbf{A}_x=\begin{pmatrix} 1 & \beta \\ 0 & 1 \end{pmatrix} Ax=(10β1)
回顾前文公式:
μ p q ′ = ∑ i = 0 p ∑ j = 0 q ( p i ) T ( q j ) a 11 i ⋅ a 12 p − i ⋅ a 21 j ⋅ a 22 q − j ⋅ μ i + j , p + q − i − j \mu'_{pq} = \sum_{i=0}^{p}\sum_{j=0}^{q} \begin{pmatrix} p \\ i \end{pmatrix}^T \begin{pmatrix} q \\ j \end{pmatrix} a^i_{11}\cdot a^{p-i}_{12}\cdot a^{j}_{21}\cdot a^{q-j}_{22}\cdot \mu_{i+j,p+q-i-j} μpq′=i=0∑pj=0∑q(pi)T(qj)a11i⋅a12p−i⋅a21j⋅a22q−j⋅μi+j,p+q−i−j
个人理解:这里提到前文公式,是为了告诉读者 μ 30 ( 2 ) \mu^{(2)}_{30} μ30(2) 是怎么求出来的。应该就是把 p , q p,q p,q 和矩阵 A x \mathbf{A}_x Ax 中的参数代入上式,从而得到 μ 30 ( 2 ) \mu^{(2)}_{30} μ30(2)。可是我代入进去的结果不对啊?这里的转置是我自己加的,不加求不了矩阵乘法啊!
我们得到:
μ 30 ( 2 ) = μ 30 ( 1 ) + 3 β μ 21 ( 1 ) + 3 β 2 μ 12 ( 1 ) + β 3 μ 03 ( 1 ) \mu^{(2)}_{30}=\mu^{(1)}_{30}+3\beta\mu^{(1)}_{21}+3\beta^2\mu^{(1)}_{12}+\beta^3\mu^{(1)}_{03} μ30(2)=μ30(1)+3βμ21(1)+3β2μ12(1)+β3μ03(1)
其中, μ p q ( 1 ) \mu^{(1)}_{pq} μpq(1) 是图像 f 1 ( x , y ) f_1(x,y) f1(x,y) 的中心矩。
令 μ 30 ( 2 ) = 0 \mu^{(2)}_{30}=0 μ30(2)=0,我们得到:
μ 30 ( 1 ) + 3 β μ 21 ( 1 ) + 3 β 2 μ 12 ( 1 ) + β 3 μ 03 ( 1 ) = 0 \mu^{(1)}_{30}+3\beta\mu^{(1)}_{21}+3\beta^2\mu^{(1)}_{12}+\beta^3\mu^{(1)}_{03}=0 μ30(1)+3βμ21(1)+3β2μ12(1)+β3μ03(1)=0
参数 β \beta β 就是通过这个式子得到的。
注意到上式在 μ 03 ( 1 ) ≠ 0 \mu^{(1)}_{03}\neq 0 μ03(1)=0 的情况下最多可以有三个根,这对于大多数自然图像来说是普遍成立的。特别地,我们可能有以下两种情况:
- 三个根中一个是实根,另外两个是复根;
- 三个根都是实根。
对于第一种情况,我们简单地取 β \beta β 为实根;对于第二种情况,我们取 β \beta β 为三个实根的中位数。参见附录,这样的选择保证了得到的标准化图像的唯一性。
当然,在一些非常不寻常的条件下,根的个数可能会发生变化。例如,当上式涉及到的所有矩都为 0 0 0 时,它将有无穷多个解。这可以发生在图像是旋转对称的时候,比如圆盘或者圆环。我们参考文献 [16] 和 [17] 来更详细地介绍一般的标准化过程。
矩阵 A y = ( 1 0 γ 1 ) \mathbf{A}_y=\begin{pmatrix} 1 & 0 \\ \gamma & 1 \end{pmatrix} Ay=(1γ01)
回顾前文公式:
μ p q ′ = ∑ i = 0 p ∑ j = 0 q ( p i ) T ( q j ) a 11 i ⋅ a 12 p − i ⋅ a 21 j ⋅ a 22 q − j ⋅ μ i + j , p + q − i − j \mu'_{pq} = \sum_{i=0}^{p}\sum_{j=0}^{q} \begin{pmatrix} p \\ i \end{pmatrix}^T \begin{pmatrix} q \\ j \end{pmatrix} a^i_{11}\cdot a^{p-i}_{12}\cdot a^{j}_{21}\cdot a^{q-j}_{22}\cdot \mu_{i+j,p+q-i-j} μpq′=i=0∑pj=0∑q(pi)T(qj)a11i⋅a12p−i⋅a21j⋅a22q−j⋅μi+j,p+q−i−j
我们得到:
μ 11 ( 3 ) = γ μ 20 ( 2 ) + μ 11 ( 2 ) \mu^{(3)}_{11}=\gamma\mu^{(2)}_{20}+\mu^{(2)}_{11} μ11(3)=γμ20(2)+μ11(2)
令 μ 11 ( 3 ) = 0 \mu^{(3)}_{11}=0 μ11(3)=0,我们得到:
γ = − μ 11 ( 2 ) μ 20 ( 2 ) \gamma=-\frac{\mu^{(2)}_{11}}{\mu^{(2)}_{20}} γ=−μ20(2)μ11(2)
因此,参数 γ \gamma γ 具有唯一的解。
矩阵 A s = ( α 0 0 δ ) \mathbf{A}_s=\begin{pmatrix} \alpha & 0 \\ 0 & \delta \end{pmatrix} As=(α00δ)
缩放参数 α \alpha α 和 δ \delta δ 的大小通过在水平和垂直方向上将图像 f 3 ( x , y ) f_3(x,y) f3(x,y) 缩放到规定的标准尺寸来确定。它们的符号都是确定的,这样 μ 50 ( 4 ) \mu^{(4)}_{50} μ50(4) 和 μ 05 ( 4 ) \mu^{(4)}_{05} μ05(4) 都是正的,可以通过水平翻转或垂直翻转来改变。
简而言之,参数 α \alpha α 和 δ \delta δ 的大小是根据规定的标准尺寸确定的,标准尺寸的大小决定了参数 α \alpha α 和 δ \delta δ 的取值。由于缩放倍数是非负的,因此参数 α \alpha α 和 δ \delta δ 的符号一定为正(?)
D 水印的影响
值得注意的是,对于水印嵌入,标准化是对原始图像进行的;对于水印提取,标准化是对含水印图像进行的。因此,重要的是设计水印信号,使其对标准化图像的影响最小。
令 w ( x , y ) w(x,y) w(x,y) 表示添加到原始图像 f ( x , y ) f(x,y) f(x,y) 中的水印信号。令 m p q ( w ) m^{(w)}_{pq} mpq(w) 表示 w ( x , y ) w(x,y) w(x,y) 的原点矩。根据前文公式:
d 1 = m 10 m 00 , d 2 = m 01 m 00 d_1=\frac{m_{10}}{m_{00}},\ d_2=\frac{m_{01}}{m_{00}} d1=m00m10, d2=m00m01
可以取 m 10 ( w ) = m 01 ( w ) = 0 m^{(w)}_{10}=m^{(w)}_{01}=0 m10(w)=m01(w)=0,使得 w ( x , y ) w(x,y) w(x,y) 对标准化过程的中心步骤没有影响。
此外,我们希望对于 p + q p+q p+q 等于 2 2 2 和 3 3 3 有 m p q ( w ) = 0 m^{(w)}_{pq}=0 mpq(w)=0,使得水印不影响其余的标准化变换。这里假设 w ( x , y ) w(x,y) w(x,y) 和 f ( x , y ) f(x,y) f(x,y) 是统计独立的,因此它们的二阶和三阶中心矩是可加的。
正如后面将要讨论的那样,水印是由零均值高斯或均匀源产生的 C D M A \mathsf{CDMA} CDMA 信号,它被添加到图像的中频 D C T \mathsf{DCT} DCT 系数中。从我们的数值例子可以看出,这样的水印几乎满足上述所描述的所有期望的性质,并且对标准化图像几乎没有影响。
E 可替代的标准化过程
上述标准化过程由一系列初等仿射变换(即剪切和缩放操作)组成。我们指出,其他变换过程也可以用类似的方式构造,以实现标准化图像中的放射变换不变性。比如以下过程:
A = ( c o s ϕ s i n ϕ − s i n ϕ c o s ϕ ) ( α 0 0 δ ) ( 1 β 0 1 ) \mathbf{A}=\begin{pmatrix} cos\phi & sin\phi \\ -sin\phi & cos\phi \end{pmatrix} \begin{pmatrix} \alpha & 0 \\ 0 & \delta \end{pmatrix} \begin{pmatrix} 1 & \beta \\ 0 & 1 \end{pmatrix} A=(cosϕ−sinϕsinϕcosϕ)(α00δ)(10β1)
其中包括:
- 1)在 x x x 方向上的剪切;
- 2 )在 x x x 方向和 y y y 方向上的缩放;
- 3)旋转 ϕ \phi ϕ 角度;
上述过程中的参数可以通过执行每个步骤的一组预定义的矩来确定。感兴趣的读者可以参考文献 [15]。
这篇关于数字水印 | 图像标准化论文:Digital Watermarking Robust to Geometric Distortions(二)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!