本文基于离散量子漫步算法和像素分块算法,提出了一种既安全又有效的量子彩色图像水印方案。离散量子漫步,作为随机行走的量子版本,可以被视作量子态与位置概率分布间的一种非线性映射关系。其独特的数学性质,也作为混沌行为的印记,使得它成为CNOT门控制位的理想选择。在水印嵌入过程中,采取了像素分块策略,通过精心比较各像素块内像素的嵌入能力,从而优化了嵌入水印后的图像视觉质量。该方案不仅增强了水印方案的安全性,同时也保证了图像的观感不受影响。 Based on the discrete quantum walk algorithm and pixel blocking algorithm, this paper proposes a safe and effective quantum color image watermarking scheme. Discrete quantum walk, as the quantum version of random walk, can be regarded as a nonlinear mapping relationship between quantum states and position probability distributions. Its unique mathematical properties, also serving as an imprint of chaotic behavior, make it an ideal choice for the control bit of the CNOT gate. In the watermark embedding process, a pixel blocking strategy is adopted, carefully comparing the embedding capabilities of pixels within each pixel block, thereby optimizing the visual quality of the image after embedding the watermark. This scheme not only enhances the security of the watermarking scheme but also ensures that the image perception is not affected.
本文基于离散量子漫步算法和像素分块算法,提出了一种既安全又有效的量子彩色图像水印方案。离散量子漫步,作为随机行走的量子版本,可以被视作量子态与位置概率分布间的一种非线性映射关系。其独特的数学性质,也作为混沌行为的印记,使得它成为CNOT门控制位的理想选择。在水印嵌入过程中,采取了像素分块策略,通过精心比较各像素块内像素的嵌入能力,从而优化了嵌入水印后的图像视觉质量。该方案不仅增强了水印方案的安全性,同时也保证了图像的观感不受影响。
量子图像处理,量子水印,量子漫步
Qianhui Wang, Zhenzhen Li, Zichen Li
School of Information Engineering, Beijing Institute of Graphic Communication, Beijing
Received: Apr. 26th, 2024; accepted: May 24th, 2024; published: May 31st, 2024
Based on the discrete quantum walk algorithm and pixel blocking algorithm, this paper proposes a safe and effective quantum color image watermarking scheme. Discrete quantum walk, as the quantum version of random walk, can be regarded as a nonlinear mapping relationship between quantum states and position probability distributions. Its unique mathematical properties, also serving as an imprint of chaotic behavior, make it an ideal choice for the control bit of the CNOT gate. In the watermark embedding process, a pixel blocking strategy is adopted, carefully comparing the embedding capabilities of pixels within each pixel block, thereby optimizing the visual quality of the image after embedding the watermark. This scheme not only enhances the security of the watermarking scheme but also ensures that the image perception is not affected.
Keywords:Quantum Image Processing, Quantum Watermarking, Quantum Walks
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着信息时代的飞速发展,人们虽然能享受前所未有的便捷与高效,但网络安全面临的风险也日益增加。网络上的信息内容时常面临被非法窥探和被泄露的威胁,这不仅涉及个人隐私,更可能触及企业机密,乃至国家安全的红线,其潜在损失难以估量。数字图像 [
文献 [
2016年,Jianzhi Sang等人基于国际上流行的量子图像表示方法NEQR [
| I 〉 = 1 2 n ∑ i = 0 2 2 n − 1 | c i 〉 ⊗ | i 〉 (1)
其中, | c i 〉 表示相应像素点的颜色信息,用二进制序列表示,用其定义如下:
| c i 〉 = | R q − 1 ⋯ R 0 G q − 1 ⋯ G 0 B q − 1 ⋯ B 0 〉 (2)
R j , G j , B j ∈ { 0 , 1 } , j = 0 , 1 , ⋯ , q − 1 (3)
| R q − 1 ⋯ R 0 〉 表示RGB图像的像素中R通道的颜色值, | G q − 1 ⋯ G 0 〉 表示G通道的颜色值, | B q − 1 ⋯ B 0 〉 表示B通道的颜色值。
| i 〉 表示相应像素点的位置信息,它表示为
| i 〉 = | y 〉 | x 〉 = | y n − 1 y n − 2 ⋯ y 0 〉 | x n − 1 x n − 2 ⋯ x 0 〉 (4)
x k , y k ∈ { 0 , 1 } , k = 0 , 1 , ⋯ , n − 1 (5)
量子漫步是经典随机游走地量子版本,在经典物理世界中,并不存在真正的随机。虽然其假设粒子或个体的每一步移动都是独立的随机事件,但在实际应用中,粒子或个体的运动可能会受到各种因素的影响,如环境因素、个体间的相互作用等。像“抛硬币”或计算机的随机数都是伪随机,即它们的结果在某种程度上是可预测的。为实现真正的随机游走,需要一个不可预测的随机变量。而量子力学通过其独特的态叠加性质,为我们提供了一种产生真正随机数的途径,即将量子力学引入到经典随机游走中,从而得到量子漫步这一量子计算的通用模型 [
离散量子漫步主要由以下几个部分构成:
位置空间(Position Space):在离散量子漫步中,位置空间是由一系列离散的点或节点构成的,每个点代表漫步者可能的位置。这通常可以表示为一个希尔伯特空间,其基态对应于不同的位置。而条件位移算子 S 的结构取决于量子漫步正在行走的图。例如,等式2-1给出了在不受限制的线上进行的离散量子漫步的条件位移算子:
S ^ = ∑ x ( | x + 1 , 0 〉 〈 x , 0 | + | x − 1 , 1 〉 〈 x , 1 | ) (6)
硬币空间(Coin Space):除了位置空间外,离散量子漫步还引入了一个额外的硬币空间。这个空间用于描述漫步者下一步可能的移动方向。在最简单的模型中,硬币可以是一个二态系统,表示向左或向右移动。二维硬币算子 C 的一般矩阵表示由等式7给出:
C = ( ρ 1 − ρ e i ω 1 − ρ e i ω − ρ e i ( ω + ϕ ) ) , 0 ≤ ρ ≤ 1 , 0 ≤ ω , ϕ ≤ π (7)
漫步者的态(Walker's State):漫步者的态是位置空间和硬币空间的张量积。这意味着漫步者的状态同时包括了它的位置信息和移动方向的信息。
演化算子(Evolution Operator):离散量子漫步的演化是通过一个么正算子来实现的。如等式8所示。
U ^ = S ^ ( C ^ ⊗ I ^ ) (8)
其中, U ^ 代表演化算子; C ^ 代表硬币算子,负责更新硬币的状态,即决定下一步的移动方向; S ^ 代表条件位移算子,根据硬币的状态来更新漫步者的位置。
测量(Measurement):在量子漫步的某个时间点,我们可以对漫步者的状态进行测量,以确定它的确切位置或其他相关信息。
综上所述,离散量子漫步系统是通过结合位置空间、硬币空间、漫步者的态、演化算子和测量来构建的。
本文在方案 [
(1) 通过在具有N个节点的圆上运行r步量子漫步,以平方等式的概率振幅来计算概率分布 P 。
r步离散量子漫步的表达形式如下所示:
| ψ 〉 t n = U ^ | ψ 〉 t 0 = S ^ ( C ^ ⊗ I ^ ) | ψ 〉 t 0 (9)
其中, C ^ 和 S ^ 分别是硬币算子和位移算子。
对于具有N个节点的圆,等式10和11是更合适的硬币算子和条件位移算子:
C = ( cos θ sin θ sin θ − cos θ ) (10)
S ^ = { | 2 , 0 〉 〈 1 , 0 | + | N , 1 〉 〈 1 , 1 | when x = 1 | 1 , 0 〉 〈 N , 0 | + | N − 1 , 1 〉 〈 N , 1 | when x = N | x + 1 , 0 〉 〈 x , 0 | + | x − 1 , 1 〉 〈 x , 1 | when x ≠ 1 , N (11)
使用等式12中给出三个硬币算子 C ^ 0 、 C ^ 1 、 C ^ 2 分别构造三种不同的演化算子 U ^ 0 、 U ^ 1 、 U ^ 2 ,并使用这三个演化算子在圆上进行一维单粒子量子漫步而产生 P 。
C 0 = ( cos θ 1 sin θ 1 sin θ 1 − cos θ 1 ) , C 1 = ( cos θ 2 sin θ 2 sin θ 2 − cos θ 2 ) , C 2 = ( cos θ 3 sin θ 3 sin θ 3 − cos θ 3 ) (12)
其中, θ 1 , θ 2 , θ 3 ∈ { 0 , 2 π } 。
U ^ 0 = S ^ ( I ^ ⊗ C ^ 0 ) , U ^ 1 = S ^ ( I ^ ⊗ C ^ 1 ) , U ^ 2 = S ^ ( I ^ ⊗ C ^ 2 ) (13)
演化算子 U ^ 0 、 U ^ 1 、 U ^ 2 按照以下步骤进行选择:生成一个二进制序列m。并设置好量子漫步的步数r。然后,量子漫步的每一步的演化算子 U ^ i j ( j ∈ { 1 , ⋯ , r } )将根据等式14中给出的规则来选择。
U ^ i j = { U ^ 0 , if m j = 0 ( i .e .,the j thbitof m isequalto 0 ) U ^ 1 , if m j = 1 ( i .e .,the j thbitof m isequalto1 ) U ^ 2 , ifthe Q W stepnumber j isgreaterthelengthof m (14)
量子漫步运行后,在位置x处找到粒子的概率可以表述如下:
P ( x , r ) = ∑ c ∈ { 0 , 1 } | 〈 x , c | ( U ^ ) r | ψ 〉 0 | 2 (15)
(2) 将 P 的元素分成 n 个集合 P i ,其中每个 P i 由 n 个数组成(即 N = n × n , i ∈ { 1 , ⋯ , n } )。此外,产生一个临时阶 n 矩阵 T ,其中 T i 是集合 P i 。换言之,就是将概率分布 P 调整大小并重塑为 n × n 矩阵。
(3) 通过等式16,将矩阵 T 变换成一个新的临时矩阵 S ,它的条目是0到255之间的整数,即我们将每个条目 t i , j ∈ [ 0 , 1 ] 映射成一个整数 s i , j ∈ { 0 , ⋯ , 255 } 。
s i , j = f l o o r ( ( t i , j × 10 8 ) mod 256 ) (16)
为了计算上述等式,我们使用了MATLAB函数 f i x ,即,
s i , j = f i x ( ( t i , j × 10 8 ) mod 256 ) (17)
最终,将矩阵 S 通过NCQI模型,产生矩阵 | K 〉 。
将量子图像 | K 〉 作为控制量子比特,对水印图像 | I 〉 运用CNOT门生成置乱后的水印图像,用 | X 〉 表示。
具体过程如下。
| R 7 x ⋯ R 0 x G 7 x ⋯ G 0 x B 7 x ⋯ B 0 x 〉 ⊗ | i 〉 = C ^ n o t ⊗ 24 | R 7 k ⋯ R 0 k G 7 k ⋯ G 0 k B 7 k ⋯ B 0 k 〉 | R 7 c ⋯ R 0 c G 7 c ⋯ G 0 c B 7 c ⋯ B 0 c 〉 ⊗ | i 〉 = | C n o t ( R 7 k , R 7 c ) ⋯ C n o t ( R 0 k , R 0 c ) ⋯ C n o t ( G 0 k , G 0 c ) ⋯ C n o t ( B 0 k , B 0 c ) 〉 ⊗ | i 〉 (18)
| X 〉 表示如下。
| X 〉 = 1 2 n ∑ i = 0 2 2 n − 1 | R 7 x ⋯ R 0 x G 7 x ⋯ G 0 x B 7 x ⋯ B 0 x 〉 ⊗ | i 〉 (19)
在嵌入水印之前,将载体图像分割成多个大小为 2 × 2 的像素块,以一个大小为 2 n × 2 n 的载体图像为例,可以被分割成 2 n − 1 × 2 n − 1 个像素块,每个像素块包含4个紧密相连的像素。为了比较每个像素块中4个像素的嵌入能力,本文通过双量子比特比较器对每个像素的颜色通道的最后三位量子比特进行测量,这一步旨在计算其中0的数量,因为0的数量多少直接反映了该像素承载额外信息的能力。完成测量后,根据每个像素最后三位中0的数量,从高到低进行排序,以此生成序列K。图1给出了大小为 4 × 4 的图像分块示例。
在载体图像被分割成多个像素块之后,按照从前往后的顺序,逐一将水印图像的像素信息嵌入到载体图像的像素块中。具体来说,水印图像的第一个像素点会被嵌入到 B 0 像素块,水印图像的第二个像素点会被嵌入到 B 1 像素块,以此类推。以水印图像中每个颜色通道值为8为例,本文会将8位的二进制序列从左到右分成4部分:分别是3位量子比特,2位量子比特,2位量子比特和1位量子比特。这样的划分是为了后续更精准的进行嵌入操作。
同时在嵌入过程中,会结合序列K,对于一个像素块中嵌入能力最强的像素点,选择其对应颜色通道值的最后3位,利用LSB算法进行水印信息的嵌入。而对于嵌入能力最弱的像素点,只修改其最后1位。至于嵌入能力位于中间的两个像素点,则分别修改它们的最后2位。图2呈现了一个大小为 1 × 1 的水印图像的R通道的嵌入示例,展示了水印信息如何被嵌入到载体图像中。
通过本文所提出的嵌入方法,能够在确保水印嵌入容量的前提下,实现最佳的视觉质量,从而在保障信息完整性的同时,提供近乎无损的图像观赏体验。
图1. (a)未分块前的载体图像,(b) 载体图像分块示例
图2. (a) 水印图像,(b) 载体图像
在水印的提取过程中,需要用到载体图像 | C 〉 和嵌入水印后的载体图像 | C W 〉 和序列K。具体提取过程如下:
将经典的载体图像转换为NCQI模型下的量子图像 | C 〉 。
在载体图像 | C 〉 和嵌入水印后的载体图像 | C W 〉 之间执行异或运算,得到量子图像 | X ′ 〉 。
结合序列K和 | X ′ 〉 ,得到量子图像 | X 〉 。
提取出的水印图像恢复成原图像的操作就是置乱操作的逆操作,需要用到提取出的量子图像 | X 〉 和控制图像 | K 〉 。得到置乱前的水印图像 | I 〉 的过程可以描述为:
| I 〉 = C − n o t ( | K 〉 , | X 〉 ) 。 (20)
在本章中,对上述水印方案在经典计算机上进行了仿真实验,所有实验均在MATLAB R2014a进行实现,实验用到大小为 256 × 256 的图像 Lena、Pepper、Sugar和Tree。
图3. (a)~(d)为原始图像,(e)~(f)为置乱后的图像,(i)~(l)为恢复后的图像
图3给出了图像置乱打的仿真结果,下文进行了直方图分析,评估置乱的效果。
图像的直方图是一种统计图表,它直观地展示了图像中像素值的分布情况。在灰度图像中,直方图特别揭示了从0到255的色调值的分布情况。在彩色图像中,可以针对每个颜色通道分别计算直方图,或者首先将图像转换为灰度图像再计算直方图。直方图分析对于评估图像置乱的效果非常重要,因为它能够清晰地反映出图像中各个像素值的出现频率。好的置乱方法应该能够确保不同置乱图像的直方图具有高度的一致性,以此来抵御潜在的统计攻击。图4给出了原始图像Lena及置乱后图像的对比的直方图,图5、图6和图7则分别给出了sugar、tree和pepper置乱后图像的直方图。而这些图像在经过3.2节中提出置乱方法处理后,其直方图呈现出高度的一致性。这一结果充分证明了本章所提的量子图像置乱方法能够有效抵御直方图分析攻击,从而确保了图像数据的安全性。
图4. Lena图像的直方图
图5. 水印图像sugar置乱后直方图
图6. 水印图像tree置乱后直方图
图7. 水印图像pepper置乱后直方图
图8. 嵌入实验结果:(a) 水印图像,(b) 载体图像,(c) 嵌入水印后的载体图像
在仿真实验中,本章以Lena、Pepper、Sugar和Tree四个图像作为载体,进行了水印嵌入的测试。视觉质量是评估水印算法的重要指标,而峰值信噪比(RGB-PSNR)经常用于评估水印图像的视觉质量 [
P S N R R G B = ( P S N R R + P S N R G + P S N R B ) / 3 , (21)
P S N R X = 10 log 10 M A X 2 M S E = 20 log 10 M A X I M S E , X ∈ { R , G , B } (22)
图8以Lena图为例,清晰地展示了嵌入水印前后的Lena图像。通过对比,可以观察到,载体图像在视觉上并未出现显著的变化,这充分证明了本章所提出的嵌入方案在保持图像视觉质量方面具有较高的性能。表1详细列出了这四个载体图像在嵌入水印后的峰值信噪比(PSNR)值。这些数值数据提供了客观的嵌入效果评价,从而再次印证了本章所提出的水印嵌入方案不仅具有良好的视觉效果,同时也具备了高度的技术可行性。
载体图像 | 水印图像 | RGB-PSNR |
---|---|---|
Lena | Pepper | 56.4591 |
Lena | Sugar | 56.5702 |
Lena | Tree | 56.7287 |
Pepper | Lena | 57.1023 |
Pepper | Sugar | 56.9806 |
Pepper | Tree | 57.2542 |
Sugar | Lena | 56.9905 |
Sugar | Pepper | 57.0100 |
Sugar | Tree | 57.2019 |
Tree | Lena | 56.7842 |
Tree | Pepper | 56.6948 |
Tree | Sugar | 56.7963 |
表1. 仿真实验的实验结果
本文深入探讨了一种创新的彩色图像量子图像水印方案,该方案以NCQI模型为基础,并涉及量子漫步算法和像素分块算法,在第3节详细阐述了水印的置乱、嵌入以及提取过程。为了验证方案的实用性和性能,本文在第4节中,通过MATLAB平台,进行了严谨的仿真实验。实验结果充分证明了在量子计算环境下,本方案不仅切实可行,而且效果显著。特别值得一提的是,本方案在图像处理过程中保持了出色的视觉质量,这一优势使得本方案在实际应用中具有更广泛的前景和潜力。
北京市教委科研计划项目(KM202010015009);北京印刷学院网络空间安全(数字版权保护技术)培育 学科建设项目(21090123010);北京印刷学院青年卓越项目(Ea202411)。
王骞卉,李祯祯,李子臣. 基于量子漫步和像素分块算法的量子彩色图像水印方案Quantum Color Image Watermarking Scheme Based on Quantum Walk and Pixel Blocking Algorithm[J]. 计算机科学与应用, 2024, 14(05): 194-205. https://doi.org/10.12677/csa.2024.145128
https://doi.org/10.1109/79.879337
https://doi.org/10.1063/1.1495073
https://doi.org/10.1109/CSA.2008.24
https://doi.org/10.1007/s10773-012-1087-9
https://doi.org/10.1007/s10773-013-1687-z
https://doi.org/10.1016/j.ins.2011.09.028
https://doi.org/10.1007/s11128-012-0423-6
https://doi.org/10.1007/s11128-013-0561-5
https://doi.org/10.1007/s11128-013-0629-2
https://doi.org/10.1007/s11128-014-0783-1
https://doi.org/10.1007/s11128-016-1260-9
https://doi.org/10.1007/s11128-017-1577-z
https://doi.org/10.1007/s11128-018-2138-9
https://doi.org/10.1016/j.physa.2022.128248
https://doi.org/10.1007/s11128-023-03926-1
https://doi.org/10.1007/s11128-013-0567-z
https://doi.org/10.1007/s11128-016-1463-0
https://doi.org/10.1145/380752.380757
https://doi.org/10.1103/PhysRevA.70.022314
https://doi.org/10.1038/srep07784
https://doi.org/10.1007/s11128-019-2386-3
https://doi.org/10.1016/j.physa.2019.123869
https://doi.org/10.1109/AQTR.2010.5520671