[发明专利]基于迭代非线性双随机相位编码的图像加密方法

权利要求书

1.一种迭代非线性双随机相位编码的图像加密方法，其特征是按如下步骤进行：

(1)加密：

(i)f(x，y)代表待加密的原始图像，R₁(x，y)和R₂(u，v)是计算机随机生成的两个相位板，分别可以具体表示成exp[2πα(x，y)]和exp[2πβ(u，v)]，α(x，y)和β(u，v)代表两个在区间[0，1]上具有均匀概率分布并且统计无关的随机矩阵，其中(x，y)和(u，v)分别表示输入平面和菲涅耳衍射输出平面的坐标，运用迭代相位恢复算法进行加密时，假定在第n-1次(n＝1，2，3…)迭代过程中已经得到两个相位板P_n(x，y)和P′_n(u，v)，当n＝1时，特别规定P_n(x，y)＝R₁(x，y)、P′_n(u，v)＝R₂(u，v)，即两个计算机随机生成的相位板被用作第一次迭代运算过程使用的两个密钥；

(ii)第n次迭代过程中，首先对相位板P_n(x，y)进行一次波长为λ，距离为z₁的菲涅耳变换，接着对变换后得到的复振幅进行取振幅和取相位操作，分别得到振幅分布g_n(u，v)和相位分布r_n(u，v)，即：

gn(u,v)=PT{FrTz1,λ[Pn(x,y)]}---(1)]]>

rn(u,v)=PR{FrTz1,λ[Pn(x,y)]}---(2)]]>

其中PT{}代表取振幅运算，即除去复振幅的相位信息，PR{}代表取相位运算，即除去复振幅的振幅信息，FrT[]代表菲涅耳变换，以某一函数U₀(x，y)为例，在波长为λ的平面光波的照射下，传播方向上距离为z处的菲涅耳衍射分布U(u，v)数学上可以表示为：

U(u,v)=FrTz,λ[U0(x,y)]=exp(jkz)jλz∫∫-∞∞U0(x,y)exp{jπλz[(u-x)2+(v-y)2]}dxdy---(3)]]>

其中k是波数，大小为式(1)的逆变换可以表示为：

U₀(x，y)＝IFrT_z，λ[U(u，v)]            (4)

其中IFrT[]代表逆菲涅耳变换；

(iii)对g_n(u，v)和P′_n(u，v)的乘积作一次波长为λ，距离为z₂的菲涅耳变换，对变换后的结果进行取振幅和取相位操作，分别得到振幅图像f′_n(x，y)和相位分布r′_n(x，y)，即：

fn′(x,y)=PT{FrTz2,λ[gn(u,v)Pn′(u,v)]}---(5)]]>

rn′(x,y)=PR{FrTz2,λ[gn(u,v)Pn′(u,v)]}---(6)]]>

(iv)如果n不大于预先设定的某一整数N，则对r′_n(x，y)和f(x，y)的乘积作一次逆菲涅耳变换，对变换后得到的复振幅分布进行取振幅和取相位操作后分别得到振幅分布g′_n+1(u，v)和相位分布P′_n+1(u，v)，接着对g′_n+1(u，v)和r_n(u，v)的乘积作逆菲涅耳变换，对变换后的结果进行取相位操作，得到相位分布P_n+1(x，y)，计算公式分别如下：

gn+1′(u,v)=PT{IFrTz2,λ[f(x,y)rn′(x,y)]}---(7)]]>

Pn+1′(u,v)=PR{IFrTz2,λ[f(x,y)rn′(x,y)]}---(8)]]>

Pn+1(x,y)=PR{IFrTz1,λ[gn+1′(u,v)rn(u,v)]}---(9)]]>

其中P′_n+1(u，v)、P_n+1(x，y)将用于下一轮迭代运算；

(v)重复步骤(ii)-(iv)，当迭代次数n达到N时，分别由式(8)、式(9)得到两个相位密钥P′_N+1(u，v)、P_N+1(x，y)，利用P′_N+1(u，v)、P_N+1(x，y)根据式(1)和式(5)得到振幅图像f′_N+1(x，y)，即

fN+1′(x,y)=PT{FrTz2,λ{PT{FrTz1,λ[Pn+1(x,y)]}Pn+1′(u,v)}}---(10)]]>

两个用于光学解密的相位板P₁(x，y)、P₂(u，v)分别为

P₁(x，y)＝P_N+1(x，y)              (11)

P2(u,v)=PN+1′(u,v)rN+1*(u,v)---(12)]]>

其中rN+1(u,v)=PR{FrTz1,λ[PN+1(x,y)]},]]>*表示相位共轭，由式(10)、式(11)和式(12)得到迭代结束后输出的振幅图像：

f′(x,y)=fN+1′(x,y)=PT{FrTz2,λ{FrTz1,λ[P1(x,y)]P2(u,v)}}---(13)]]>

(2)解密：

(i)对加密过程中得到的P₁(x，y)作一次波长为λ，距离为z₁的菲涅耳变换，变换后的结果为

h1(u,v)=FrTz1,λ[P1(x,y)]---(14)]]>

(ii)h₁(u，v)与加密过程中得到的另一相位板P₂(u，v)相乘后作一次波长为λ，距离为z₂的菲涅耳变换，对变换后的结果进行取振幅操作，得到最终的解密结果D(x，y)，即

D(x,y)=PT{FrTz2,λ[h1(u,v)P2(u,v)]}---(15)]]>

综合以上过程，解密结果可以表述为：

D(x,y)=PT{FrTz2,λ{FrTz1,λ[P1(x,y)]P2(u,v)}}=f′(x,y)---(16)]]>

因此，迭代过程中使用的波长λ、传播距离z₁和z₂以及加密过程生成的两个相位板P₁(x，y)、P₂(u，v)都是光学解密过程中所必需的密钥；解密得到的图像D(x，y)与迭代加密过程中计算得到的振幅图像f′(x，y)相同。