基于相位恢复算法和干涉原理的数字图像加密方法
【专利摘要】一种基于相位恢复算法和干涉原理的数字图像加密方法。按如下两大步骤进行：一是加密：利用一块随机相位模板作为加密密钥，将加密密钥与待加密图像分别作为相位恢复算法中输入平面与输出平面的约束条件，通过迭代运算得到加密结果，即输入面上不同于加密密钥的另一个相位模板和傅立叶频域上的相位模板；二是解密：解密密钥与加密过程中生成的输入面上的相位模板叠加后进行傅立叶变换，变换结果经傅立叶频域上的相位模板调制后作逆傅立叶变换，对变换后的结果取振幅得到的图像即为解密结果。本发明提出的图像加密方法具有迭代算法收敛速度快、解密简单快速和安全性高的优点。
【专利说明】基于相位恢复算法和干溃原理的数字图像加密方法 【技术领域】
[0001] 【【背景技术】】本发明设及一种基于相位恢复算法和干设原理的数字图像加密方法， 属于信息安全技术领域。 【【背景技术】】
[0002] 近年来，基于光电信息处理的图像加密技术越来越受到人们的关注。其中，基于4f 系统的光学双随机相位编码技术的应用较为广泛。但双随机相位编码技术在加密过程中需 要同时记录信息的振幅和相位，在图像解密时光学效率不高。相位恢复算法是一种通过可 测量的光场强度确定光场相位分布的方法，属于解决"逆问题"的重要技术，其设计思想在 光学信息安全中也就有着广泛应用。1996年，Johnson和化asher利用相位恢复算法将图像 信息加密到两块相位模板中，解密则在经典的光学双随机相位编码系统中完成。基于光学 干设原理的图像加密方法的研究则是最近几年光学图像安全处理研究领域出现的新热点。 此类方法具有加密运算简单，光学解密方便的特点。只需运用光强探测器，如CCD就可W在 解密系统的输出面上方便地得到解密结果。但是直接利用干设方法得到的相位模板存在信 息泄露问题，当作为加密结果的其中任意一块相位模板放置在解密装置中时，在输出面上 就能看到原图像的轮廓。
[0003] 许多基于相位恢复算法或干设原理的加密方法都是采用数值计算的方式将图像 加密成几块相位分布不同的相位模板，解密时采用光学手段。光学解密所需硬件成本高，灵 活性低，对相位模板的空间排列精度要求非常高，因此在实际应用上目前仍受到很大限制。 基于光学信息安全处理技术的加密方法与发展程度已经非常高的数字技术相结合，能有效 推动光学信息安全技术的实际应用。 【
【发明内容】
】
[0004] 本发明要解决的技术问题是提供基于相位恢复算法和干设原理的数字图像加密 方法。
[0005] 解决上述技术问题采用如下技术措施:基于相位恢复算法和干设原理的数字图像 加密方法按如下步骤进行：
[0006] (1)加密：
[0007] 。)^义，7)代表待加密图像，3^，7)、义；^(乂，>〇:和^：^1^〇，1/)是计算机生成的 Ξ个随机相位模板，它们分别可W具体表示成e邱[;r0(x,y) ]、exp[;ri(x,y)]和e邱[r2(u， 乂）]，其中^，7)和(11，乂)分别表示空间域和傅立叶频域的坐标^〇^，7)^1^，7)和〇(11^) 代表Ξ个在区间[0,2π]上具有均匀概率分布并且统计无关的随机矩阵，他们的尺寸大小与 ^义，7)相同，运用迭代傅立叶变换进行加密时，3(^，7)作为加密密钥使用，足1^(^乂，^)和 则仅用于初次迭代运算，第j次（j = l，2,3···)迭代运算能够生成第j + 1次迭代 运算所需要的两个相位模板
[000引（ii)当进行第j次迭代运算时，首先对f(x，y)和哀:的乘积作傅立叶变 换，得到波函数g//,W(w，i/)，即
[0009]
(1)
[0010]其中FT代表傅立叶变换，接着，的W(W，V)的振幅和反户（W，F)的乘积作一次 逆傅立叶变换，得到新的波函数(；χ，，即
[0011]
(25
[0012]其中PT代表取振幅运算，即除去复振幅的相位信息，IFT代表逆傅立叶变换；
[OOU] 和加密密钥R(x，y)计算得到一个相位模板PiW(x，y)，计算公 式为：
[0014]
(3)
[001引其中α是一个大于零但数值很小的常数，相位模板PiW(X，y)与加密密钥R(X，y)相 加后作一次傅立叶变换，得到一个新的波函数於，即
[0019] 其中PR代表取相位运算，即除去复振幅的振幅信息；
[0020] (iv)利用y兮'>O，^0的相位信息和由式（l)得到的(6l("如v)计算得到新的相位 模板巧八(//，|/)，计算公式为：
[0021]
(6)
[00剖其中*表示复共辆运算，对巧V)与的./)(//，I/)的乘积作一次逆傅立叶变换 后得到波函数的J')(x，y)，即
[0023]
(7)
[0024] 对变换后得到的复振幅进行取相位和取振幅操作，分别得到用于下次迭代的相位 模板必1(川，片，切和振幅分布//W〇，_y)，即
[0027]由此，在第j次迭代运算过程中分别由式(5)、式(8)得到了第j+1次迭代运算过程 所需的两个相位模板《·''+"〇，!/)和巧州>仁，少，），分别由式(3)、式(6)得到的两块相位 模板巧和巧八(//，1/)则作为本次（即第j次}迭代运算图像fU'y}的两个加密结 果；
[00巧](V)重复步骤（i)-(iv)，直到定义的均方差(Mean Square E：;r;ro;r，MSE)达到设计的 精度值或者循环次数达到设置的最大迭代次数为止，均方差定义为：
[0029]
(10)
[0030] 其中M，N表示图像的尺寸，f(m，n)和巧;)分别表示两幅振幅图像在像素点 (m，n)的值，通过MSE可W反映出所进行的迭代运算的收敛性，当迭代次数总共完成J次时迭 代运算终止，根据式(3)、式(6)分别得到两块相位模板PiW(x，y)和，它们就是 原始图像f(x，y)的加密结果，分别重新记作Pi(x，y)和Ρ2(μ，ν);
[0031] (2)解密：
[0032] (i)R(x，y)作为解密密钥与相位模板Pi(x，y)相加后作一次傅立叶变换，变换后得 到波函数恥(μ, V)，即
[0033] φ<ι(μ,ν)=ΡΤ[Ρι(χ,γ)+Κ(χ,γ)] (11)
[0034] (?)Φ<ι(μ，ν)与另一相位模板Ρ2(μ，ν)相乘后作一次逆傅立叶变换，对变换后得到 的复振幅进行取振幅操作得到解密结果，即
[003引 fd(x，y)=PT{IFT[ilM(y，v)P2(y，v)]} (12)
[0036] 由式(4)、（7)、（9)、（11)和（12)可^得到：
[0037]
(13)
[0038] 也就是说，解密得到的图像fd(x，y)就是加密过程第J次迭代运算得到的振幅图像 fd'\x,y) 〇
[0039] 本发明的有益效果在于：（1)本发明结合相位恢复算法与光学干设原理，加密最终 生成两块相位模板，不需要生成振幅板，并且克服了大部分基于干设原理的加密方法存在 的相位模板信息泄露问题；（2)加密过程和解密过程均采用数字方式，不需要进行相位的光 学全息记录；（3)迭代算法收敛速度快，运算过程中引入一块随机相位模板作为加密密钥 (同时也是解密密钥），增加了系统的安全性。 【【附图说明】】
[0040] 图1为加密过程流程图。
[0041] 图2为解密过程流程图。
[0042] 图3(a)待加密图像^义，7);(6)加密密钥3山7)的相位分布图；（。化1">托3/)的 相位分布图；（d)公f C//，V：)的相位分布图。
[0043] 图4迭代30次后得到的加密结果和解密图像。（a)相位模板Pi; (b)相位模板P2;(c) 解密图像fd。
[0044] 图5原图f(x，y)和第j次迭代运算后得到的图像/;/W(JC，J；)之间的MSE值与迭代运 算次数j之间的关系图；
[0045] 图6相位模板信息泄露问题测试结果。（a)密钥R和相位模板Pi用于解密时得到的 结果；（b)密钥R和相位模板P2用于解密时得到的解密图像；（c)Pi和P2用于解密时得到的解 密图像;（d)P2单独用于解密时得到的解密图像.
[0046] 图7(a)解密密钥R发生错误时得到的解密图像；（b)相位模板Pi发生错误时得到的 解密图像;（C)相位模板P2发生错误时得到的解密图像。 【【具体实施方式】】
[0047] 本发明所述方法的【具体实施方式】如下：
[004引（1)图像的加密过程(如图1所示)分如下几个步骤：
[0049] (i)f (x，y)代表待加密图像，R(x，y)、
是计算机生成的 Ξ个随机相位模板，它们分别可W具体表示成e邱[;r0(x,y) ]、exp[;ri(x,y)]和e邱[r2(u， 乂）]，其中^，7)和(11^)分别表示空间域和傅立叶频域的坐标^〇^，7)、〇^，7)和^(11^) 代表Ξ个在区间[0,2π]上具有均匀概率分布并且统计无关的随机矩阵，他们的尺寸大小与 ^义，7)相同，运用迭代傅立叶变换进行加密时，3^，7)作为加密密钥使用，必1^(；_^，>〇和 则仅用于初次迭代运算，第j次（j = l，2,3···)迭代运算能够生成第j + 1次迭代 运算所需要的两个相位模板
[0化0] (ii)当进行第j次迭代运算时，首先对f(x，y)和足的乘积作傅立叶变 换，得到波函蠻
其中FT代表傅 立叶变换，接着，的八批1/)的振幅和足V)的乘积作一次逆傅立叶变换，得到新的 波函i
，其中PT代表取 振幅运算，即除去复振幅的相位信息，IFT代表逆傅立叶变换；
[0051] 和加密密钥R(x，y)计算得到一个相位模板PiW(x，y)，计算公式
小的常数，相位模板PiW(x，y)与加密密钥R(x，y)相加后作一次傅立叶变换，得到一个新的 波函I
里 为： 其中α是一个大于零但数值很 ， 行取相位操作，得到相位模巧
其中PR 代表取相位运算，即除去复振幅的振幅信息；
[0052] (iv)利用的相位信息和由式（1)得到的Φι^(μ，ν)计算得到新的相位 模板巧^从，·^^ ，计算公式为：
， 其中*表示复共辆运算，对
的乘积作一次逆傅立叶变换后得到 波函数〇，_)0,印
，对变换后得到的复 振幅进行取相位和取振幅操作，分别得到用于下次迭代的相位模板和振幅分 布
由此，在第j次迭代运算过程中分别由式(5)、式(8)得到了第j+l次迭代运算过程所需的两 个相位模板
，分别由式(3)、式(6)得到的两块相位模板PiW (x，y)和巧W〇，F)则作为本次(即第j次)迭代运算图像f(x，y)的两个加密结果；
[0053] (V)重复步骤（i)-(iv)，直到定义的均方差(Mean Square E：;r;ro;r，MSE)达到设计的 精度值或者循环次数达到设置的最大迭代次数为止，均方差定义为：
，其中M，N表示图像的尺寸，f (m，n)和/dW(^W，巧^分别表示两幅振幅图像在像素点(m，n)的值，通过MSE可W反映出所进 行的迭代运算的收敛性，当迭代次数总共完成J次时迭代运算终止，根据式(3)、式(6)分别 得到两块相位模板PiW(x，y)和
，它们就是原始图像f(x，y)的加密结果，分别 重新记作Pi(x，y)和Ρ2(μ，ν);
[0054] (2)图像的解密过程(如图2所示）：
[0055] (i)R(x，y)作为解密密钥与相位模板Pi(x，y)相加后作一次傅立叶变换，变换后得 到波函数恥（4^)，即恥（4卢）=尸1'[口1^，7)+1?(义，7)];
[0056] (ii)恥(μ, V)与另一相位模板Ρ2(μ，ν)相乘后作一次逆傅立叶变换，对变换后得到 的复振幅进行取振幅操作得到解密结果，即fd(x，y)=PT{IFT[恥(4^化2(4^)]}，由式(4)、 (7)、（9)、（11)和（12)可^得到：
也就是说，解密得到的图像fd (X，y)就是加密过程第J次迭代运算得到的振幅图像/?0, V)。
[0057] 下面结合实施例和附图对本发明的内容进行进一步的解释。
[005引首先，选择大小为200X200的灰度图"Lena"作为待加密图像，归一化后分别如图3 (a)所示，加密密钥R(x，y)的相位分布如图3(b)所示，它在数值上是一个区间在[0,2π]上具 有均匀概率分布且大小为200X200的随机矩阵。计算机生成的用于初次迭代运算的两个随 机相位
分别如图3(c)和图3(d)所示。
[0059]根据加密流程图图1进行图像加密运算，实验中设定最大迭代次数J = 30,式(3)中 的常数α大小设置为1 X 10-7,经30次迭代后加密运算终止，此时图像的信息已经被加密到如 图AW所示的巧3〇>(X，y)和图4(ω所示的巧3。)(//，1/)两块相位模板中，运两块相位模板 分别就是加密结果Pi(x，y)和Ρ2(μ，ν)。迭代过程中得到的图像
如图4(c)所示， 由式(13)可知，它也是相位模板Pi(x，y)和Ρ2(μ，ν)的解密结果fd(x，y)。可W看出，解密图像 质量极高，无法用肉眼分辨出它与原图之间的区别。
[0060] 图5给出了解密图像的均方差函数随迭代次数的变化情况。当加密算法终止运算 时即30次迭代完成后，均方差函数已经接近ΙΟΛ运说明本发明提出的加密算法的收敛速度 很快。
[0061] 下面考察本发明的安全性。由式（11)和式（12)可知，如果只使用相位模板Pi进行 解密，即式（11)的R和式（12)中的P2分别用数值0和1代替，得到的解密结果显然是一个元素 大小全为1的矩阵。当使用密钥R和相位模板Pi进行解密时，即式(12)中P2用数值1代替，得到 的解密结果如图6(a)所示；当密钥R和相位模板P2用于解密时，即式（11)中Pi用数值0代替， 得到的解密结果如图6(b)所示；当只使用相位模板Pi和P2用于解密时，即式(11)中R用数值0 代替，得到的解密结果如图6(c)所示;而只使用相位模板P2进行解密，即式(12)中恥用数值 1代替，得到的结果则如图6(d)所示。上述分析和实验结果表明，密钥R、相位模板Pi和相位 模板P2是实现图像解密的Ξ个必备要素，缺一不可，相位模板不存在信息泄露问题。
[0062] 进一步测试系统的安全性。当用于解密的密钥R、相位模板Pi和Ρ2Ξ者中的其中任 意一个发生错误时，如图7所示，得到的Ξ种解密结果均为噪声分布图，运些实验结果反映 出本发明提出的加密方法具有很高的安全性。
【主权项】
1. 一种基于相位恢复算法和干涉原理的数字图像加密方法，其特征是按如下步骤进 行： (1)加密： (i) f (x，y)代表待加密图像，R(x，y)和是计算机生成的三个 随机相位模板，它们分别可以具体表示成61口[1'()(1，7)]、61口[1'1(1，7)]和61口|^2(11，￥)]，其 中（x，y)和(u，v)分别表示空间域和傅立叶频域的坐标，r〇(x，y)、ri(x，y)和r2(u，v)代表三 个在区间[〇,2π]上具有均匀概率分布并且统计无关的随机矩阵，他们的尺寸大小与f( x，y) 相同，运用迭代傅立叶变换进行加密时，R(x，y)作为加密密钥使用，和 则仅用于初次迭代运算，第j次（j = l，2,3···)迭代运算能够生成第j+Ι次迭代 运算所需要的两个相位模板JC[(V+1>和i：f+1); (ii) 当进行第j次迭代运算时，首先对f(x，y)和尺的乘积作傅立叶变换，得 到波函数ν/?Ακμ)，即 ψ\η{ιι,ν) = Ft[/(x, v)</!(x,>〇J Cl) 其中FT代表傅立叶变换，接着，p 〇, 1/)的振幅和ATf 〇,幻的乘积作一次逆傅 立叶变换，得到新的波函数vi73 即 ψ^{χ.γ) = (2) 其中PT代表取振幅运算，即除去复振幅的相位信息，IFT代表逆傅立叶变换； (iii) (X, 和加密密钥R(x，y)计算得到一个相位模板Piw(X，y)，计算公式为：其中α是一个大于零但数值很小的常数，相位模板与加密密钥R(x，y)相加后 作一次傅立叶变换，得到一个新的波函数^/^7)(//,1/)，即 y/, η(μ, ν) = FT[/f n (x, y) + R(x, y)\ ⑷ 对以力(//, v)进行取相位操作，得到相位模板(从i/)，即 0 川)(，，,v) = PR[〇，v)] (5) 其中PR代表取相位运算，即除去复振幅的振幅信息； (iv) 利用以1//,0的相位信息和由式⑴得到的ihu)(y，v)计算得到新的相位模板 计算公式为： Ρ^ η {μ, ν) - {PR^7) (μ,ν)]} Ρ^-[ψΙ η(μ, ν}\ (6) 其中*表示复共辄运算，对/^ω〇,ν)与的乘积作一次逆傅立叶变换后得 至ij波函数,即 ψ{/'(χ,χ) = \Fl^Ah{u,v)[in{u,v)\ (?) 对变换后得到的复振幅进行取相位和取振幅操作，分别得到用于下次迭代的相位模板 < y+1)(工，J)和振幅分布/)力(·υ)，即 Kli+')(x,y) = PR\i/^)(x,y)\ (8) ./:;.7.) (X，少)=PT[W) (.y)] (9； 由此，在第j次迭代运算过程中分别由式(5)、式(8)得到了第j+l次迭代运算过程所需 的两个相位模板J^+1)Ca,0和尺y+U〇,：^)，分别由式⑶、式(6)得到的两块相位模板 ?1(巧^7)和/^(//，1〇则作为本次(即第」次)迭代运算图像以1，7)的两个加密结果； (v)重复步骤（i)-(iv)，直到定义的均方差(Mean Square Error，MSE)达到设计的精度 值或者循环次数达到设置的最大迭代次数为止，均方差定义为：其中M，N表示图像的尺寸，f(m，n)和分别表示两幅振幅图像在像素点(m，n) 的值，通过MSE可以反映出所进行的迭代运算的收敛性，当迭代次数总共完成J次时迭代运 算终止，根据式(3)、式(6)分别得到两块相位模板PPUj)和0，它们就是原始 图像f(x，y)的加密结果，分别重新记作Pi(x，y)和Ρ2(μ，ν); (2)解密： (i)R(x，y)作为解密密钥与相位模板PKxj)相加后作一次傅立叶变换，变换后得到波 函数Φ<?(μ，ν)，即 itd(y,v)=FT[Pi(x,y)+R(x,y)] (11) 与另一相位模板Ρ2(μ，ν)相乘后作一次逆傅立叶变换，对变换后得到的复 振幅进行取振幅操作得到解密结果，ΒΡ fd(x，y)=PT{IFT[ilM(y，v)P2(y，v)]} (12) 由式(4)、（7)、（9)、（11)和（12)可以得到： v)=/j,y)(x, v) (13) 也就是说，解密得到的图像fd(X，y)就是加密过程第J次迭代运算得到的振幅图像
【文档编号】G06T1/00GK106067156SQ201610211568
【公开日】2016年11月2日
【申请日】2016年4月6日 公开号201610211568.7, CN 106067156 A, CN 106067156A, CN 201610211568, CN-A-106067156, CN106067156 A, CN106067156A, CN201610211568, CN201610211568.7
【发明人】汪小刚, 陈翼翔, 徐一清
【申请人】浙江农林大学