删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于空间角度复用和双随机相位的多图像光学加密方法

本站小编 Free考研考试/2021-12-29

摘要:提出了基于空间角度复用和双随机相位的多图像光学加密新方法. 加密过程中, 首先将原始图像进行随机相位调制和不同距离的菲涅耳衍射; 其次, 将携带调制后图像的参考光与携带随机相位且具有不同立体角的参考光相干叠加, 产生干涉条纹; 最后, 将不同方向的干涉条纹叠加形成复合加密图像. 解密为加密的逆过程, 将复合加密图像置于空间滤波和菲涅耳衍射系统中, 经过不同相位密钥解调和正确距离的菲涅耳衍射完成解密, 得到多幅解密图像. 该方法可以同时对多幅图像进行高效的加密, 计算简单、安全可靠、抗噪声能力强. 利用相关系数评估了该方法的加密效果, 并通过仿真实验验证了该方法的有效性和安全性.
关键词: 多图像加密/
空间角度复用/
数字全息/
干涉加密

English Abstract


--> --> -->
光学图像加密技术已在信息安全领域显示出巨大的应用潜力, 成为****们研究的热点. 1995年, Refregier和Javidi[1]首次提出了安全性高、鲁棒性强的双随机相位光学加密技术, 此后一系列派生的光学图像加密方法被不断提出, 如分数傅里叶变换加密系统[2]、菲涅耳变换加密系统[3]、分数Hartley变换加密系统[4]、混沌[5]和置乱加密系统[6]、相移干涉加密系统[7]、联合变换相关器加密系统[8]、叠层衍射成像加密系统[9]、偏振加密系统[10]、数字全息加密系统[11]和计算全息加密系统[12]等. 但是, 上述加密方法仅针对单个图像, 随着大数据技术的快速发展和信息传输能力的不断增强, 传统的单图像加密传输已经不能满足日益增长的信息需求. 因此, 越来越多的****开始研究多图像加密技术[13-15].
多图像加密技术的要点在于图像在加密过程中如何合成, 合成方法直接影响整个算法的计算效率和最终解密图像的质量. 目前的多图像光学加密技术主要基于复用、数字全息、压缩感知、混沌和特殊光学变换等. 例如, Situ和Zhang[16]引入了波长复用来实现多图像加密. Xu等[17]提出了一种基于随机振幅板和菲涅耳全息图的多图像加密方法. Deepan等[18]将压缩感知技术应用于基于双随机相位密钥的空间复用加密系统, 实现多图像加密. Tang等[19]结合位平面分解和混沌映射算法对多幅图像进行了加密, Kong等[20]使用级联分数傅里叶变换将多个图像叠加成单个图像进行加密. 目前, 由于基于单一技术的加密方法具有局限性, 这些多图像加密方法大多是基于多种技术手段的组合, 尽管提高了加密图像的数量, 但也增加了系统的复杂性. 同时, 数据处理的时间和复杂度也随着加密能力的增加而增加. 此外, 受限于解密实验中随机相位密钥逐像素对齐的要求[21], 这些方法主要集中在数字系统或光电混合系统, 难以光学实现.
针对上述问题, 本文提出一种基于空间角度复用和双随机相位的光学多图像加密方法. 该方法利用数字全息的空间角度复用技术实现多个待加密图像的同时加密; 利用基于干涉原理的双随机相位光学图像加密技术, 将两个随机相位密钥分别置于物光束和参考光束中, 并以干涉条纹的振幅形式保存和传输, 简化了加密系统; 解密过程是正确密钥调制的空间滤波和菲涅耳衍射过程, 解决了传统方法解密实验中随机相位密钥难以逐像素对齐的问题. 此外, 本文给出了可行的实验加密和解密系统, 并通过仿真实验验证了光学加密方法的可行性. 该方法具有存储效率高、计算简单等特点, 在提高信息传输效率和多用户认证方面具有重要的应用前景.
本文提出的多图像加密系统如图1所示, 是包含两个4f成像系统和两个相位型空间光调制器(SLM)的马赫-曾德尔干涉系统. 首先待加密图像由4f系统成像到第一个SLM1上, 被加载到SLM1上的随机相位密钥调制. 调制后的图像经过衍射距离为${z_i}$的菲涅耳衍射过程后与携带加载在SLM2上的第二个随机相位密钥的参考光相干叠加, 形成加密干涉条纹, 并由电荷耦合器(CCD)记录. 然后, 将待加密图像和第一个随机相位去除, 由同一CCD记录第二个随机相位与平行光的干涉条纹. 该干涉条纹可再现恢复第二个随机相位密钥及其共轭. 变换待加密图像时, 改变双随机相位和菲涅耳衍射距离, 并通过旋转CCD和调整光路改变参考光的立体角, 形成间距和方向不同的加密干涉条纹. 最后, 将多个全息图叠加构成复合加密图像.
图 1 多图像光学加密系统(SLM是空间光调制器, $f$是透镜焦距, ${z_i}$是菲涅耳衍射距离, $R(\theta )$表示CCD旋转后与x轴的夹角, ${\alpha _i}$是物光O与参考光R的夹角)
Figure1. Optical setup of multiple-image encryption process. SLM is spatial light modulator, $f$ is focal length, ${z_i}$ is the distance of Fresnel diffraction, $R(\theta )$ is the rotation angle of CCD, ${\alpha _i}$ is angle between object light O and reference light R.

本文选取8幅原始图像, 分别是“A”和“B”2个字母, “Lena”和“baboon”两个256灰度图像, 以及“光”、“学”、“全”和“息”4个汉字, 其表达式分别为$a(x, y)$, $b(x, y)$, $c(x, y)$, $d(x, y)$, $e(x, y)$, $f(x, y)$, $g(x, y)$$h(x, y)$, 如图2所示.
图 2 (a)?(h) 8幅待加密图像
Figure2. (a)?(h) Multiple-image to be encrypted.

在如图1所示的光学系统中, 每个原始图像被平行光垂直照射, 并被第一随机相位${p_{{\rm{1}}i}}$($i$为图像序号)调制, 第一随机相位可以表示为
$\begin{split}& {p_{11}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\;{p_{12}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\\& {p_{13}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\;{p_{14}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\\& {p_{15}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\;{p_{16}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\\& {p_{17}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\;{p_{18}} = {\rm{exp[j{\text{π}}\, rand}}(x,y)],\end{split}$
其中, ${\rm{rand}}(x, y)$表示0到1的随机矩阵. 为了尽可能地恢复原始图像信息提高解密图像质量, 第一个随机相位的动态范围选取为$0\text{—}{\text{π}}$. 由于不同的衍射距离可以提高多图像加密的安全性, 加密过程中的菲涅耳衍射距离分别设置为${z_1} = 0.5$ m, ${z_2} = 0.55$ m, ${z_3} = 0.6$ m, ${z_4} = 0.65$ m, ${z_5} = 0.7$ m, ${z_6} = 0.75$ m, ${z_7} = 0.8$ m和${z_8} = 0.85$ m. 经过随机相位调制和菲涅耳衍射后的光场分布可表示为
$\begin{split}& {a_1}(x,y) = {\rm{FrT}}{z_1}{\rm{[}}a(x,y){p_{11}}(x,y),\;{\lambda _1}],\\& {b_1}(x,y) = {\rm{FrT}}{z_2}{\rm{[}}b(x,y){p_{12}}(x,y),\;{\lambda _2}],\\& {c_1}(x,y) = {\rm{FrT}}{z_3}{\rm{[}}c(x,y){p_{13}}(x,y),\;{\lambda _3}],\\& {d_1}(x,y) = {\rm{FrT}}{z_4}{\rm{[}}d(x,y){p_{14}}(x,y),\;{\lambda _4}],\\& {e_1}(x,y) = {\rm{FrT}}{z_5}{\rm{[}}e(x,y){p_{15}}(x,y),\;{\lambda _5}],\\& {f_1}(x,y) = {\rm{FrT}}{z_6}{\rm{[}}f(x,y){p_{16}}(x,y),\;{\lambda _6}],\\& {g_1}(x,y) = {\rm{FrT}}{z_7}{\rm{[}}g(x,y){p_{17}}(x,y),\;{\lambda _7}],\\& {h_1}(x,y) = {\rm{FrT}}{z_8}{\rm{[}}h(x,y){p_{18}}(x,y),\;{\lambda _8}],\end{split}$
其中, ${\rm{FrT}}{{{z}}_i}[ \cdot, \lambda ]$表示衍射距离为${z_i}$、波长为${\lambda _i}$的菲涅耳衍射.
然后, (2)式表示的光场与具有不同空间角度且携带第二个随机相位的参考光束相干叠加. 如图1的内插图所示, 参考光的空间角用$({\alpha _i}, \;{\theta _i})$表示, 其中, ${\alpha _i}$是参考光束与z轴之间的夹角, $\theta {}_i$ 是CCD逆时针旋转的角度. 当加密更多图像时, 参考光空间角度可相应调整. 本文选取参考光的立体角度$({\alpha _i}, \;{\theta _i})$分别为$(3^\circ, 0^\circ )$, $(3^\circ, 45^\circ )$, $(3^\circ, 90^\circ )$, $(3^\circ, 135^\circ )$, $(6^\circ, 0^\circ )$, $(6^\circ, 45^\circ )$, $(6^\circ, 90^\circ )$$(6^\circ, 135^\circ )$, 相应的参考光表达式为
$\begin{split}& {R_1} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}3^\circ }}{{{\lambda _1}}} + {p_{{{21}}}}\left(x,y\right)\right)\right],\\& {R_2} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}3^\circ }}{{{\lambda _2}}} + {p_{22}}\left(x,y\right)\right)\right],\\& {R_3} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}3^\circ }}{{{\lambda _3}}} + {p_{23}}\left(x,y\right)\right)\right],\\& {R_4} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}3^\circ }}{{{\lambda _4}}} + {p_{24}}\left(x,y\right)\right)\right],\\& {R_5} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}6^\circ }}{{{\lambda _5}}} + {p_{25}}\left(x,y\right)\right)\right],\\& {R_6} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}6^\circ }}{{{\lambda _6}}} + {p_{26}}\left(x,y\right)\right)\right],\\& {R_7} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}6^\circ }}{{{\lambda _7}}} + {p_{27}}\left(x,y\right)\right)\right],\\& {R_8} = {\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}6^\circ }}{{{\lambda _8}}} + {p_{28}}\left(x,y\right)\right)\right],\end{split}$
其中${p_{2 i}}(x, y)$为第二随机相位密钥. 为了简化解密过程, 令8幅待加密图像对应的第二随机相位密钥和激光波长相同, 即
$\begin{split}{p_2}_i(x,y)\, & = {p_2}(x,y) = {\rm{exp}}[2{\rm{j{\text{π}} \,rand(}}x,y)] \\ &= {\rm{exp}}[2{\rm{j{\text{π}} }}M(x,y)],\end{split}$
${\lambda _i} = \lambda = 532\;{\rm{nm}}.$
将8个待加密图像分别与(3)式表示的参考光束干涉. 为了降低带宽和采样点, 提高重建图像的质量, 用均匀场代替原来的偏置分量. 因此, 形成的干涉场分布分别为
$\begin{split}{I_1}\left(x,y\right) =\, & 1 + {a_1}\left(x,y\right)R^ * _1 + a^*_1\left(x,y\right){R_1}\\=\, & 1 + {a_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + a^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_2}\left(x,y\right) =\,& 1 + {b_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + b^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_3}\left(x,y\right) =\, & 1 + {c_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + c^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_{\rm{4}}}\left(x,y\right) =\, & 1 + {d_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + d^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_{\rm{5}}}\left(x,y\right) =\, & 1 + {e_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + e^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_{\rm{6}}}\left(x,y\right) =\, & 1 + {f_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + f^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_{\rm{7}}}\left(x,y\right) =\, & 1 + {g_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + g^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\\{I_{\rm{8}}}\left(x,y\right) =\, & 1 + {h_1}\left(x,y\right){\rm{exp}}\left[ - {\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right] + h^ * _1\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + M\left(x,y\right)\right)\right],\end{split}$
其中*表示图像共轭. 8幅原始图像产生的加密干涉条纹如图3所示, 其中干涉条纹方向与参考光的空间角${\theta _i}$参数相关, 而干涉条纹周期与参考光的空间角${\alpha _i}$参数相关.
图 3 (a)—(h)8幅原始图像对应的加密干涉条纹
Figure3. (a)-(h) Encrypted interference fringes corresponding to 8 original images.

最后将8个加密干涉条纹叠加构成最终的加密图像如图4所示, 其表达式为
图 4 加密结果图
Figure4. Encrypted image.

$\begin{split} I(x,y) =\,& {I_{1}}(x,y) + {I_{2}}(x,y) + {I_{3}}(x,y) + {I_{4}}(x,y) \\&+ {I_{5}}(x,y) + {I_{6}}(x,y) + {I_{7}}(x,y) + {I_{8}}(x,y). \end{split}$
加密结果图是一幅接近随机分布的灰度图像, 相对于传统图像加密方法的加密复值图像结果, 更便于存储和传输. 该图完全隐藏了所有原始图像的信息和加密图像个数, 在存储和传输过程中更具有一般性.
多图像解密为加密的逆过程, 在图5所示的系统中完成. 该系统包含两个4f成像系统和一次菲涅耳衍射过程. 在第一个4f系统的入射面放置加密结果图(图4), 出射面放置一个相位型SLM, 加载随机相位密钥${p_3}(x, y)$.
图 5 多图像解密系统
Figure5. Decryption system of multiple-image.

其中, CCD记录的第二个随机相位与平面波的干涉条纹再现恢复得到${p_2}$或其共轭${p_2}^*$, 可作为解密相位密钥, 首先应用${p_3} = {p_2}^ * $, 相当于将$p{}_3$紧贴加密图像放置. 经过随机相位密钥${p_3}$调制后的光场分布为
$\begin{split}& I\left(x,y\right) \times {p_3}\left(x,y\right)\\= & \left[{I_1}\left(x,y\right) + {I_2}\left(x,y\right) + {I_3}\left(x,y\right) + {I_4}\left(x,y\right) + {I_5}\left(x,y\right) + {I_{\rm{6}}}\left(x,y\right) + {I_7}\left(x,y\right) + {I_8}\left(x,y\right)\right] \times p_2^ * \left(x,y\right) \\ = \;& p_2^ * \left(x,y\right) + {a_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(x\frac{{{\rm{sin}}3^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + a_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}x\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {b_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + b_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x + y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {c_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(y\frac{{{\rm{sin}}3^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + c_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}y\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {d_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + d_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x - y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {{\rm{e}}_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(x\frac{{{\rm{sin}}6^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + e_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}x\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {f_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(\left(x + y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + f_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x + y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {g_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(y\frac{{{\rm{sin}}6^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + g_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}y\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \\& + p_2^ * \left(x,y\right) + {h_1}\left(x,y\right){\rm{exp}}\left[ - {\rm{j}}2{\text{π}}\left(\left(x - y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } + 2M\left(x,y\right)\right)\right] + h_1^ * \left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x - y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda }\right]. \end{split} $
从(8)式可见, 经过解密相位密钥调制后, 每个图像的0级和–1级信息携带随机相位, 经过傅里叶变换后成为弥散整个频谱面的噪声, 而+1级信息经过傅里叶变换后成为原始图像共轭的正确频谱, 频谱图如图6(a)所示.
图 6 (a) $p_2^* $为解密密钥时频谱图, 其中F1—F8为相应的滤波器; (b) ${p_2}$为解密密钥时频谱图
Figure6. (a) Decrypted spectrum with right key $p_2^ * $, where F1-F8 are filters; (b) decrypted spectrum with key ${p_2}$.

图6(a)所示频谱图可见, 每个原始图像的+1级频谱被恢复, 频谱位置受参考光的立体角度$({\alpha _i}, \;{\theta _i})$控制, 0级和–1级频谱被调制为噪声弥散到整个频谱空间. 为了进行对比, 本文将${p_2}$作为解密密钥时的频谱放置于图6(b), 可见此时每个原始图像的–1级频谱被恢复, 其位置受参考光的立体角度$({\alpha _i}, \;{\theta _i})$控制, 0级和+1级频谱被调制为噪声弥散到整个频谱空间. 此时(8)式变为
$ \begin{split} & I\left(x,y\right) * {p_3}\left(x,y\right) \\ = \,&\left[{I_1}\left(x,y\right) \!+\! {I_2}\left(x,y\right) \!+\! {I_3}\left(x,y\right) \!+\! {I_4}\left(x,y\right) \!+\! {I_5}\left(x,y\right) \!+\! {I_6}\left(x,y\right) \!+\! {I_7}\left(x,y\right) \!+\! {I_8}\left(x,y\right)\right] \times {p_2}\left(x,y\right) \\ =\,&{p_2}\left(x,y\right) \!+\! {a_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}x\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \!+\! a _1^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}3^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ & \!+\! {p_2}\left(x,y\right) \!+\! {b_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}\left(x \!+\! y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \!+\! b_1 ^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x \!+\! y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ & \!+\! {p_2}\left(x,y\right) \!+\! {c_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}y\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \!+\! c _1^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}3^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ &\!+\! {p_2}\left(x,y\right) \!+\! {d_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}\left(x \!-\! y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda }\right] \!+\! d_1 ^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x \!-\! y\right)\frac{{{\rm{sin}}3^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ & \!+\! {p_2}\left(x,y\right) \!+\! {e_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}x\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \!+\! e _1^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(x\frac{{{\rm{sin}}6^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ &\!+\! {p_2}\left(x,y\right) \!+\! {f_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}\left(x \!+\! y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \!+\! f _1^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x \!+\! y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ & \!+\! {p_2}\left(x,y\right) \!+\! {g_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}y\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \!+\! g _1^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(y\frac{{{\rm{sin}}6^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right] \\ &\!+\! {p_2}\left(x,y\right) \!+\! {h_1}\left(x,y\right){\rm{exp}}\left[ \!- {\rm{j}}2{\text{π}}\left(x \!-\! y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda }\right] \!+\! h_1 ^ *\left(x,y\right){\rm{exp}}\left[{\rm j}2{\text{π}}\left(\left(x \!-\! y\right)\frac{{{\rm{sin}}6^\circ }}{\lambda } \!+\! 2M\left(x,y\right)\right)\right]. \end{split} $
经过图6(a)所示的相应滤波器滤波, 并经过正确衍射距离${z_i}$的菲涅耳衍射后光场分布为
$ \begin{split} & {a^ * }(x,y)p_{11}^ * (x,y) = {\rm{FrT}}{z_1}[{{\rm F}_1}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_1}[a_1^ * (x,y),\lambda ], \\ & {b^ * }(x,y)p_{12}^ * (x,y) = {\rm{FrT}}{z_2}[{{\rm F}_2}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_2}[b_1^ * (x,y),\lambda ], \\ & {c^ * }(x,y)p_{13}^ * (x,y) = {\rm{FrT}}{z_3}[{{\rm F}_3}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_3}[c_1^ * (x,y),\lambda ], \\ & d^ * (x,y)p_{14}^ * (x,y) = {\rm{FrT}}z_4[{{\rm F}_4}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_4}[d_1^ * (x,y),\lambda ], \\ & e^ * (x,y){p_{15}}^ * (x,y) = {\rm{FrT}}{z_5}[{{\rm F}_5}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_5}[e_1^ * (x,y),\lambda ], \\ & f^ * (x,y)p_{16}^ * (x,y) = {\rm{FrT}}{z_6}[{{\rm F}_6}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_6}[f_1^ * (x,y),\lambda ], \\ & g^ * (x,y)p_{17}^ * (x,y) = {\rm{FrT}}{z_7}[{{\rm F}_7}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_7}[g_1^ * (x,y),\lambda ], \\ & h^ * (x,y)p_{18}^ * (x,y) = {\rm{FrT}}{z_8}[{{\rm F}_8}(I(x,y)p_2^ * (x,y)),\lambda] = {\rm{FrT}}{z_8}[h_1^ * (x,y),\lambda]. \end{split} $
从(10)式可见, 经过滤波和菲涅耳衍射后, 用CCD接收到原始图像的共轭, 完成解密, 获得解密图像如图7所示.
图 7 (a)—(h) 8个图像的正确密钥解密结果
Figure7. (a)-(h) 8 Decrypted images with all right keys.

图7可见, 用正确的密钥和滤波器可解密获得原始图像的共轭图像, 由于每个图像的0级和–1级频谱被调制为噪声弥散到整个频谱空间, 通过滤波器后形成了噪声, 对解密结果有一定的干扰. 但弥散在整个频谱的噪声通过滤波器的占比仍然较小, 对解密结果影响不大, 因此仍获得了令人满意的解密结果.
2
4.1.安全性分析和测试
-->加密之后的图像通过公共通信信道传输后, 存在信息失真的多种可能性. 因此, 为进一步说明和验证本文所提出的多图像加密方法的可行性和有效性, 引入相关系数$CC$来评价解密结果的质量, 定义如下:
$CC \!=\! \frac{{\sum\limits_m {\sum\limits_n {(o(m,n) - \overline o )(o'(m,n) - \overline {o'} )} } }}{{\sqrt {\!\left(\!{\sum\limits_m {\sum\limits_n {{{(o(m,n) \!- \!\overline o )}^2}} } }\!\right)\!\!\left(\!{\sum\limits_m {\sum\limits_n {{{(o'(m,n) \!-\! \overline {o'})}^2}} } }\! \right)} }}, $
其中, $o(x, y)$$o'(x, y)$分别表示原始图像和解密图像, $\overline o $$\overline {o'} $分别表示$o(x, y)$$o'(x, y)$的平均值.
本文提出的多图像光学加密方法中, 除随机相位密钥以外, 菲涅耳衍射距离和激光波长都可作为附加密钥. 首先, 当所有密钥都正确时解密结果(图7)与原始图像(图2)的CC值分别为CC1 = 0.9239, CC2 = 0.9077, CC3 = 0.8432, CC4 = 0.8072, CC5 = 0.9029, CC6 = 0.8883, CC7 = 0.8929, CC8 = 0.9040. 可见, 当所有密钥都正确时通过解密可以很好地获得原始图像的信息. 将解密结果中图7(c)图7(d)与其他二值图像解密结果对比发现, 灰度图像解密结果质量要差一些, 其CC值也相对较低, 可见本文提出的多图像加密方法更适用于二值图像, 因为输入平面的振幅和相位随机性对解密结果影响很大. 图8为随机相位密钥p3错误、其他密钥都正确时的解密结果图, 以原始图“A”(图8(a))和原始图“光”(图8(b))为例.
图 8 (a)随机相位密钥$p_2^ * $错误时原始图“A”解密结果; (a')随机相位密钥$p_2^ * $错误且无滤波器时原始图“A”解密结果; (b)随机相位密钥$p_2^ * $错误时原始图“光”解密结果; (b')随机相位密钥$p_2^ * $错误且无滤波器时原始图“光”解密结果
Figure8. Decrypted results with wrong key $p_2^ * $: (a) For original image “A”; (a') for original image “A” without filter; (b) for original image “光”; (b') for original image “光” without filter.

图8可见, 当随机相位密钥错误时, 有滤波器解密得到图8(a)图8(b), 解密结果图与原始图像的相关系数仅有$CC = 0.0873$$CC = 0.0914$, 无滤波器时解密得到图8(a')图8(b'), 均类似均匀噪声, 没有原始图像的任何信息.
图9为菲涅耳衍射距离波长$\lambda $错误其他密钥正确时的解密结果, 以原始图“B”(图9(a))和原始图“学”(图9(b))为例.
图 9 (a)波长$\lambda = 633$ nm其他密钥正确时原始图“B”的解密图像; (b)原始图“学”的解密图像
Figure9. Decrypted results with wrong key $\lambda = 633$ nm: (a) For original image “B”; (b) for original image “学”.

图9可见, 当光源波长错误$\lambda = 633$ nm时, 解密得到图9(a)图9(b), 解密结果图与原始图像的相关系数仅有$CC = 0.3038$$CC = 0.3391$, 解密结果图像质量极大地下降, 并且对每个图像应用不同波长进行加密时更具有迷惑性, 因此, 光源波长可作为该多图像加密方法的附加密钥.
解密结果与另一个重要参数菲涅耳衍射距离${z_i}$依赖关系如图10所示(以原始图像“息”为例). 当所有密钥正确时, $CC$值在0.9左右, 可以获得良好的解密结果. 虽然本文将随机相位的动态范围从${\text{0—2π}}$压缩到${\text{0—π}}$以提高解密图像的质量, 大大降低了衍射距离的灵敏度, 但$CC$值仍随${z_i}$偏差的增大而迅速减小, 当误差大于${\rm{1}}{\rm{.5}}\;{\rm{cm}}$时, $CC$小于0.2, 无法区分插图中所示的解密图像. 因此, 该多图像光学加密方法对菲涅耳衍射距离${z_i}$是高度敏感的, ${z_i}$可作为附加密钥来提高安全性.
图 10${p_1}$的动态范围为${\text{0—π}}$时原始图“息”的CC随菲涅耳衍射距离${z_8}$的变化
Figure10. The $CC$ dependence on ${z_8}$ when the dynamic range of ${p_1}$ is 0-π.

图9图10可见, 附加密钥波长$\lambda $和菲涅耳衍射距离${z_i}$错误且与正确密钥差别较小时, 虽然$CC$值迅速减小, 但仍然有原始图像轮廓出现. 这是因为将随机相位${p_1}$的动态范围从${\text{0—2π}}$压缩到${\text{0—π}}$后提高了解密图像的质量, 降低了图像加密系统对附加密钥的灵敏度. 因此, 本文将随机相位${p_1}$的动态范围提高到${\text{0—1.5π}}$, 验证解密结果与菲涅耳衍射距离${z_i}$依赖关系, 结果如图11所示.
图 11${p_1}$的动态范围为${\text{0—1.5π}}$时原始图“息”的$CC$随菲涅耳衍射距离${z_8}$的变化
Figure11. The $CC$ dependence on ${z_8}$ when the dynamic range of ${p_1}$ is 0-1.5π.

通过对比图10图11, 当随机相位${p_1}$的动态范围提高时, 正确解密密钥${z_i}$应用时, 解密结果质量降低, 对应的$CC$值由高于0.9下降到0.8以上, 但是$CC$值随密钥${z_i}$误差变化的曲线斜率明显变大, 即图像加密系统对附加密钥${z_i}$的灵敏度大幅提高. 因此, 随机相位${p_1}$的动态范围应合理选取, 既要保证图像加密系统的安全性, 又要兼顾解密图像的质量.
2
4.2.加密系统的复用容量分析
-->对于多图像加密而言, 多图像加密系统的加密容量是一个关键参数. 文献[16]给出了复用容量的概念, 即达到设定的图像解密质量评判阈值时, 该项复用技术可以加密的最大原始图像数量. 为了分析本文提出的多图像加密系统的复用容量, 以同时加密多个二值图像“A”字母为例, 设定所有解密图像的$CC$平均值为0.9, 即$\overline {CC} = 0.9000$为图像解密质量阈值, 其中$\overline \bullet $表示平均值. 同时加密二值图像“A”字母的个数可通过改变参考光立体角度$({\alpha _i}, \;{\theta _i})$${\alpha _i}$${\theta _i}$参数设定, 经过模拟实验和计算获得该值为32, 当阈值设定为$\overline {CC} = 0.8000$时, 得到够同时加密二值图像“A”字母的个数为64. 可见本文提出的多图像光学加密系统具有较高的复用容量.
本文提出了基于空间角度复用和双随机相位的多图像光学加密方法. 该方法利用数字全息的空间角度复用技术将多个图像加密为单个灰度图像, 易于保存和传输; 利用基于干涉原理的双随机相位光学图像加密技术, 将双随机相位分别置于物光和参考光, 降低了加密系统的复杂度; 同时将第二个随机相位板放置于参考光束可克服传统方法解密实验中随机相位密钥难以逐像素对齐的问题, 而且参考光可进行随机相位复用, 进一步提高了加密系统的容量, 同时通过多图像加密系统的复用容量分析, 获得了该系统的加密容量. 因此, 该方法可以同时对多幅图像进行高效的加密, 计算简单、安全可靠、抗噪声能力强, 在信息安全领域具有重要的应用价值.
相关话题/图像 系统 空间 光学 质量

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于超强耦合量子点-纳米机械振子系统的全光学质量传感
    摘要:提出一种复合量子点-纳米机械振子系统,该系统以半导体芯片为基底,量子点嵌入倒置半导体圆锥纳米线的底端,通过光学抽运-探测技术来驱动量子点-纳米机械振子系统,研究该系统中的相干光学特性.通过探测吸收谱给出确定机械振子频率和量子点-纳米机械振子耦合强度的全光学方法.此外,基于该系统理论上提出一种在 ...
    本站小编 Free考研考试 2021-12-29
  • 多体系统中相干资源的一般化理论
    摘要:量子相干理论是一类重要的量子资源理论,其自由操作是各种类型的非相干操作.在单体相干资源理论中,最大相干态是最重要的量子资源态,它被定义为在非相干操作下可以转化为任何其他纯态的量子态.但是,这一情形在多体系统中发生了巨大改变:不仅在有些相干度量下不存在唯一的最大相干态,而且在有些非相干操作下几乎 ...
    本站小编 Free考研考试 2021-12-29
  • 经典场驱动对量子系统生存概率的影响
    摘要:考虑一个受经典场驱动的二能级系统与零温玻色子库相互作用,研究经典场驱动对量子Zeno效应和量子反Zeno效应中量子系统存活概率的影响.结果表明,经典场驱动可以降低量子系统的有效衰减率,即提高量子系统的存活概率.此外,环境的欧姆性对于提高量子系统的存活概率也起着重要作用,设置适当的环境欧姆参数可 ...
    本站小编 Free考研考试 2021-12-29
  • 腔自旋波混合系统的研究进展
    摘要:近年来腔自旋波混合系统引起人们的研究兴趣.基于自旋波体系的优点,有望建立一个以自旋波量子为核心的、实现不同物理系统之间信息传递的平台.本文简要介绍了腔自旋波混合系统的发展进程,阐明自旋波量子与微波腔光子的耦合机制;着重介绍了近期在腔自旋波混合系统中关于非线性和赝厄米性方面的研究进展,其中包括非 ...
    本站小编 Free考研考试 2021-12-29
  • 基于转角样品杆的脉冲强磁场电输运测量系统
    摘要:测量物质在不同外加磁场方向下的电输运性质是近年来研究拓扑量子材料的一种重要实验方法,为探索物质的新奇特性提供了独特的视角和手段.研究表明,在超高强度的脉冲强磁场下,材料的电输运研究可能扩展至量子极限区域,将观察到更加丰富的物理现象.而现有的电输运测量系统中,稳态场下的样品测量杆受限于尺寸和材料 ...
    本站小编 Free考研考试 2021-12-29
  • 基于全保偏光纤结构的主振荡脉冲非线性放大系统
    摘要:提出了基于全保偏光纤结构的主振荡脉冲非线性放大系统,该系统由基于半导体可饱和吸收镜锁模的直线型光纤振荡器、二级放大结构脉冲非线性光纤放大器和具有负色散的单模传导光纤的脉冲压缩器构成.通过此系统获得了中心波长为1560nm,重复频率为200MHz的超短激光脉冲,脉冲半高全宽为44fs,单脉冲能量 ...
    本站小编 Free考研考试 2021-12-29
  • 空间电子辐射环境中绝缘介质电荷沉积特性及陷阱参数研究综述
    摘要:空间电子辐射环境中绝缘介质充放电特性与介质表面电荷交换过程或内部电荷迁移过程密切相关.介质表面/内部电荷运动很大程度上取决于材料的微观特性,空间电荷与陷阱是反映绝缘介质微观特性的重要参数.本文综述了电子辐射环境中绝缘介质内部空间电荷和陷阱的形成、作用机理、测量方法、存在的问题及国内外研究现状. ...
    本站小编 Free考研考试 2021-12-29
  • 石墨烯莫尔超晶格体系的拓扑性质及光学研究进展
    摘要:当两个晶格常数或晶格转角不同的二维材料叠加在一起时会出现周期性的莫尔条纹结构,这种莫尔超晶格形成了一个新的二维周期势,可以大大改变原有体系的物理性质.最近石墨烯与石墨烯、石墨烯与六方氮化硼形成的莫尔超晶格提供了一个非常有趣的体系,在该体系中石墨烯的电子能带结构发生了根本性的改变,在原本的能谷处 ...
    本站小编 Free考研考试 2021-12-29
  • 声子系统中弹性波与热输运的拓扑与非互易现象
    摘要:声子是晶格集体振动模式二次量子化之后的准粒子激元.在声子的框架下,可以对固体中的力学、弹性波以及热现象进行统一描述.随着对固体系统认识的提高,声子成为补充和替代电子、光子的另一种操控固体器件的重要手段.其中,对声子体系中弹性波和热传导的调控在理论和应用上都具有非常巨大的价值.弹性波作为信息载体 ...
    本站小编 Free考研考试 2021-12-29
  • 微光学元件宽带光源照明下无色散边界条件探讨
    摘要:基于宽带光源微光学元件(如集成导光板)在衍射色散方面的设计需求,本文构建了宽带光源微光学元件衍射理论分析模型,探讨分析了衍射光谱的色度规律特性,提出并定义了能准确定量衡量衍射光束色散程度的色散量C,同时明确给出了零色散的边界判据点.通过对研制的矩形位相光栅进行测试分析,所得的光谱色度特性规律与 ...
    本站小编 Free考研考试 2021-12-29