留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于轨道角动量全息和频移的大容量光学信息加密技术

郑伟 张迪 原昊 于娜娜 席思星 王桂林 马帅 王晓雷 郞利影

郑伟, 张迪, 原昊, 于娜娜, 席思星, 王桂林, 马帅, 王晓雷, 郞利影. 基于轨道角动量全息和频移的大容量光学信息加密技术[J]. 红外与激光工程, 2023, 52(7): 20230313. doi: 10.3788/IRLA20230313
引用本文: 郑伟, 张迪, 原昊, 于娜娜, 席思星, 王桂林, 马帅, 王晓雷, 郞利影. 基于轨道角动量全息和频移的大容量光学信息加密技术[J]. 红外与激光工程, 2023, 52(7): 20230313. doi: 10.3788/IRLA20230313
Zheng Wei, Zhang Di, Yuan Hao, Yu Nana, Xi Sixing, Wang Guilin, Ma Shuai, Wang Xiaolei, Lang Liying. High capacity optical information encryption technology based on OAM holography and frequency shift[J]. Infrared and Laser Engineering, 2023, 52(7): 20230313. doi: 10.3788/IRLA20230313
Citation: Zheng Wei, Zhang Di, Yuan Hao, Yu Nana, Xi Sixing, Wang Guilin, Ma Shuai, Wang Xiaolei, Lang Liying. High capacity optical information encryption technology based on OAM holography and frequency shift[J]. Infrared and Laser Engineering, 2023, 52(7): 20230313. doi: 10.3788/IRLA20230313

基于轨道角动量全息和频移的大容量光学信息加密技术

doi: 10.3788/IRLA20230313
基金项目: 国家自然科学基金项目(11904073,61875093);河北省高等学校科学技术研究项目(BJ2020028);邯郸市科学技术研究与发展计划项目(19422031008-3);山西省科技创新项目(2019L0877,2022L476);来晋科研专项项目(QZX-2018012);运城学院博士科研启动项目(YQ-2017015);山西省研究生教育教学改革课题(2022YJJG263);运城市科技计划项目(YCKJ-202268);运城学院学科建设项目
详细信息
    作者简介:

    郑伟,男,讲师,博士,主要从事光学信息处理和太赫兹功能器件方面的研究

    通讯作者: 席思星,男,副教授,博士,主要从事光场调控和光学信息安全方面的研究。; 王晓雷,男,教授,博士,主要从事光场调控和太赫兹技术方面的研究。
  • 中图分类号: O438.1; TP391

High capacity optical information encryption technology based on OAM holography and frequency shift

Funds: Natural Science Foundation of China (11904073, 61875093); Science and Technology Project of Hebei Education Department (BJ2020028); Handan Science and Technology Research and Development Program (19422031008-3); Shanxi Provincial Science and Technology Innovation Project (2019L0877, 2022L476); Research Special Project of Shanxi (QZX-2018012); The doctoral research initiation project of Yuncheng College (YQ-2017015); Shanxi Province Graduate Education and Teaching Reform Project (2022YJJG263); Yuncheng Science and Technology Plan Project (YCKJ-202268); Discipline Construction Project of Yuncheng College
  • 摘要: 提出了一种基于轨道角动量全息(Orbital Angular Momentum, OAM)和频移的大容量光学信息加密方法。该方法实现了对多个图像信息的并行加密。首先,对多幅原始图像进行采样,采样阵列的采样间隔取决于具有不同拓扑荷数的螺旋相位的空间频率。然后,多个采样图像信息经过随机相位调制、傅里叶变换和频移相位调制后相干叠加构成轨道角动量保留全息图。最后,将不同拓扑荷的螺旋相位分别编码到轨道角动量保留全息图中,得到轨道角动量选择全息图,进行相干叠加后构成最终的单个加密全息图。解密时,轨道角动量复合选择全息图被加载到空间光调制器上,用包含特定拓扑荷数的涡旋光束照射,并经过傅里叶变换获得多个解密信息。该加密系统具有极高的加密灵活性和极大的加密容量,不仅可以在同一拓扑荷下,设计不同的频移因子来并行加密一组多个图像信息,还可以利用不同拓扑荷对多组图像信息进行加密。该方法将涡旋光束的模式设定为一个新的光学密钥,极大地提高了光学加密系统的安全性。此外,该光学加密方法中,待加密图像信息的尺寸不受空间光调制器的像元数量限制,极大地提高了光学实现信息加密的可行性和有效性。仿真实验结果表明该方法具有较高的安全性、抗噪性和抗剪切能力。
  • 图  1  原始图像信息

    Figure  1.  Original image information

    图  2  螺旋相位板空间频率及采样阵列的示意图

    Figure  2.  Schematic diagram of spatial frequency and sampling array of spiral phase plates

    图  3  不同拓扑荷下采样常数的示意图

    Figure  3.  Sampling constants diagram of under different topological charges

    图  4  OAM保留全息图获取过程的示意图

    Figure  4.  Schematic diagram of OAM-preserved hologram acquisition process

    图  5  OAM复合选择全息图获取过程的示意图

    Figure  5.  Schematic diagram of OAM composite selective hologram acquisition process

    图  6  解密光路示意图

    Figure  6.  Schematic diagram of the decryption optical path

    图  7  当$m{\text{ = 1,\;2}}$的OAM光束照射时的解密结果

    Figure  7.  Decryption result when $m{\text{ = 1,\;2}}$ OAM beam is irradiated

    图  8  当$m{\text{ = 1}}$的OAM光束照射时的解密结果

    Figure  8.  Decryption result when $m{\text{ = 1}}$ OAM beam is irradiated

    图  9  当$m{\text{ = 2}}$的OAM光束照射时的解密结果

    Figure  9.  Decryption result when $m{\text{ = 2}}$ OAM beam is irradiated

    图  10  当$m \ne {\text{1,\;2}}$的OAM光束照射时的解密结果

    Figure  10.  Decryption result when $m \ne {\text{1,\;2}}$ OAM beam is irradiated

    图  11  抗剪切性分析结果

    Figure  11.  Shear resistance simulation results

    图  12  抗噪性分析结果

    Figure  12.  Noise immunity simulation results

    表  1  原始图像信息和密文图像信息的相关性系数

    Table  1.   Correlation coefficient of original image and ciphertext image

    Image
    Horizontal correlation

    Vertical correlation

    Diagonal correlation
    Frog0.99770.99790.9963
    Tree0.97150.97110.9504
    “optics”0.86610.93830.8415
    “像”0.90730.91060.8139
    Flower0.98810.99910.9977
    “光”0.98840.97910.9716
    Panda0.99760.99090.9841
    Chessboard0.98810.98610.9743
    Encrypted text−0.00130.01640.0072
    下载: 导出CSV

    表  2  明文敏感性分析

    Table  2.   Plaintext sensitivity analysis

    Index change when pixel
    increases by 0.1
    Index changes when
    two pixels are exchanged
    NPCRUACINPCRUACI
    99.85%34.75%99.84%34.58%
    下载: 导出CSV

    表  3  原始图像和解密图像的相似性

    Table  3.   Similarity of original image and decrypted image

    Similarity indexFrogTree“optics”“像”Flower“光”PandaChessboard
    SSIM0.97030.96050.96450.95090.98290.97570.96990.9742
    CC0.99430.99320.99100.99240.99400.99220.99540.9936
    下载: 导出CSV
  • [1] 付时尧,黄磊,吕燕来等.光束轨道角动量谱的测量技术研究进展(特邀)[J].红外与激光工程,2021,50(09):9-23. doi:  10.3788/IRLA20210145

    Fu Shiyao, Huang Lei, Lv Yanlai, et al. Advances on the measurement of orbital angular momentum spectra for laser beams (Invited) [J]. Infrared and Laser Engineering, 2021, 50(9): 20210145. (in Chinese) doi:  10.3788/IRLA20210145
    [2] Xiao Qiang, Ma Qian, Yan Tao, et al. Orbital-angular-momentum-encrypted holography based on coding information metasurface [J]. Advanced Optical Materials, 2021, 9(11): 2002155. doi:  10.1002/ADOM.202002155
    [3] Li Fajing, Ding Hao, Nie Shouping, et al. Multiple-image encryption using phase jump gradient factors -based OAM multiplexing holography [J]. Optics and Lasers in Engineering, 2023, 160: 107303. doi:  10.1016/j.optlaseng.2022.107303
    [4] Wang Jian, Yang Jengyuan, Fazal Irfan M, et al. Terabit free-space data transmission employing orbital angular momentum multiplexing [J]. Nature Photonics, 2012, 6(7): 488-496. doi:  10.1038/nphoton.2012.138
    [5] Bozinovic Nenad, Yue Yang, Ren Yongxiong, et al. Terabit-scale orbital angular momentum mode division multiplexing in fibers [J]. Science, 2013, 340(6140): 1545-1548. doi:  10.1126/science.1237861
    [6] Yan Yan, Xie Guodong, Lavery J, et al. High-capacity millimetre-wave communications with orbital angular momentum multiplexing [J]. Nature Communications, 2014, 5(1): 4876. doi:  10.1038/ncomms5876
    [7] 南久航, 韩一平. 双路多进制涡旋光通信[J]. 光学学报, 2021,41(12):1206901. doi:  10.3788/AOS202141.1206001

    Nan Jiuhang, Han Yiping. Dual-channel multiband vortex optical communication [J]. Acta Optica Sinica, 2021, 41(12): 1206001. (in Chinese) doi:  10.3788/AOS202141.1206001
    [8] 苏志锟, 王发强, 路轶群, 等. 基于光子轨道角动量的密码通信方案研究[J]. 物理学报, 2008, 57(5): 3016-3021. doi:  10.3321/j.issn:1000-3290.2008.05.064

    Su Zhikun, Wang Faqiang, Lu Yiqun, et al. Study on quantum cryptography using orbital angular momentum states of photons [J]. Acta Physica Sinica, 2018, 57(5): 3016-3021. (in Chinese) doi:  10.3321/j.issn:1000-3290.2008.05.064
    [9] Zhou Hongqiang, Wang Yongtian, Li Xin, et al. A deep learning approach for trustworthy high-fidelity computational holographic orbital angular momentum communication [J]. Applied Physics Letters, 2021, 119(4): 044104. doi:  10.1063/5.0051132
    [10] Shi Zijian, Wan Zhensong, Zhan Ziyu, et al. Super-resolution orbital angular momentum holography [J]. Nature Communications, 2023, 14(1): 1869. doi:  10.1038/s41467-023-37594-7
    [11] Xie Xiangsheng, Chen Yongzhu, Yang Ken, et al. Harnessing the point-spread function for high-resolution far-field optical microscopy [J]. Physical Review Letters, 2014, 113(26): 263901. doi:  10.1103/PhysRevLett.113.263901
    [12] Kozawa Yuichi, Matsunaga Daichi, Sato Shunichi. Superresolution imaging via superoscillation focusing of a radially polarized beam [J]. Optica, 2018, 5(2): 86-92. doi:  10.1364/OPTICA.5.000086
    [13] Lavery Martin P J, Speirits F C, Barnett S M, et al. Detection of a spinning object using light’s orbital angular momentum [J]. Science, 2013, 341(6145): 537-540. doi:  10.1126/science.1239936
    [14] Kravets V G, Schedin F, Jalil R, et al. Singular phase nano-optics in plasmonic metamaterials for label-free single-molecule detection [J]. Nature Materials, 2013, 12(4): 304-309. doi:  10.1038/nmat3537
    [15] Li Ye, Yu Lin, Zhang Yixin. Influence of anisotropic turbulence on the orbital angular momentum modes of Hermite-Gaussian vortex beam in the ocean [J]. Optics Express, 2017, 25(11): 12203-12215. doi:  10.1364/OE.25.012203
    [16] Fang Xinyuan, Ren Haoran, Gu Min. Orbital angular momentum holography for high-security encryption [J]. Nature Photonics, 2020, 14(2): 102-108. doi:  10.1038/s41566-019-0560-x
    [17] Ahmed Hammad, Kim Hongyoon, Zhang Yuebian, et al. Optical metasurfaces for generating and manipulating optical vortex beams [J]. Nanophotonics, 2022, 11(5): 941-956. doi:  10.1515/nanoph-2021-0746
    [18] Ren Haoran, Fang Xinyuan, Jang Jaehyuck, et al. Complex-amplitude metasurface-based orbital angular momentum holography in momentum space [J]. Nature Nanotechnology, 2020, 15(11): 948-955. doi:  10.1038/s41565-020-0768-4
    [19] Ren Haoran, Fang Xinyuan, Jang Jaehyuck, et al. 3D-printed complex-amplitude metasurface for orbital angular momentum holography[C]//Meta 2021, 2021: 1-2.
    [20] Ren Haoran, Maier Stefan A. Complex amplitude metasurfaces for orbital angular momentum multiplexing holography[C]//2020 Conference on Lasers and Electro-Optics Pacific Rim (CLEO-PR), Sydney, NSW, Australia, 2020: 1-2.
    [21] Zhou Hongqiang, Sain Basudeb, Wang Yongtian, et al. Polarization-encrypted orbital angular momentum multiplexed metasurface holography [J]. ACS Nano, 2020, 14(5): 5553-5559. doi:  10.1021/acsnano.9b09814
    [22] Fang Xinyuan, Yang Haocheng, Yao Wenzhe, et al. High-dimensional orbital angular momentum multiplexing nonlinear holography [J]. Advanced Photonics, 2021, 3(1): 015001.
  • [1] 赵丽娟, 吴雨静, 徐志钮.  高阶轨道角动量传输光纤设计及传输特性研究(内封底文章) . 红外与激光工程, 2024, 53(3): 20240007-1-20240007-13. doi: 10.3788/IRLA20240007
    [2] 李阳, 范晨晨, 郝修路, 马小雅, 姚天甫, 许将明, 曾祥龙, 周朴.  高功率涡旋拉曼光纤激光器 . 红外与激光工程, 2023, 52(6): 20230292-1-20230292-6. doi: 10.3788/IRLA20230292
    [3] 闫丽凝, 吕宏, 贺蕊, 庞令阳, 孔妍.  矢量涡旋光束经光阑-透镜系统的OAM与偏振特性 . 红外与激光工程, 2022, 51(12): 20220250-1-20220250-9. doi: 10.3788/IRLA20220250
    [4] 任元, 吴昊, 王琛, 刘政良, 刘通, 熊振宇.  量子涡旋陀螺若干关键参数的仿真计算(特邀) . 红外与激光工程, 2022, 51(4): 20220004-1-20220004-10. doi: 10.3788/IRLA20220004
    [5] 郑伟, 席思星, 王桂林, 李永宏, 姜其畅.  结合计算全息和频移的JTC系统多图像光学加密方法 . 红外与激光工程, 2022, 51(5): 20220175-1-20220175-7. doi: 10.3788/IRLA20220175
    [6] 曹康, 姜成昊, 朱精果, 杜鹃, 乔治, 叶征宇.  激光多普勒移频特性研究 . 红外与激光工程, 2021, 50(11): 20210116-1-20210116-8. doi: 10.3788/IRLA20210116
    [7] 付时尧, 黄磊, 吕燕来, 高春清.  光束轨道角动量谱的测量技术研究进展(特邀) . 红外与激光工程, 2021, 50(9): 20210145-1-20210145-15. doi: 10.3788/IRLA20210145
    [8] 于策, 王天枢, 张莹, 林鹏, 郑崇辉, 马万卓.  大气湍流信道中OAM光束与高斯光束传输性能的实验研究 . 红外与激光工程, 2021, 50(8): 20200400-1-20200400-10. doi: 10.3788/IRLA20200400
    [9] 王花花, 马健, 杨智新, 杜浩然, 高禄, 张泽.  热光非定域螺旋相衬成像(特邀) . 红外与激光工程, 2021, 50(9): 20210450-1-20210450-7. doi: 10.3788/IRLA20210450
    [10] 吕浩然, 白毅华, 叶紫微, 董淼, 杨元杰.  利用超表面的涡旋光束产生进展(特邀) . 红外与激光工程, 2021, 50(9): 20210283-1-20210283-16. doi: 10.3788/IRLA20210283
    [11] 汪莎, 张志成, 邓国亮, 周寿桓.  直接输出的超短脉冲轨道角动量涡旋光产生技术研究进展(特邀) . 红外与激光工程, 2020, 49(12): 20201061-1-20201061-11. doi: 10.3788/IRLA20201061
    [12] 李子乐, 周舟, 梁聪玲, 郑国兴.  融合计算全息术与纳米印刷术的多功能超表面研究进展 . 红外与激光工程, 2020, 49(9): 20201036-1-20201036-15. doi: 10.3788/IRLA20201036
    [13] 邓子岚, 涂清安, 李向平.  多维度超表面及其在信息加密防伪上的应用 . 红外与激光工程, 2020, 49(9): 20201034-1-20201034-16. doi: 10.3788/IRLA20201034
    [14] 白秀丽, 陈鹤鸣, 张凌菲.  轨道角动量模传输的圆环形光子晶体光纤 . 红外与激光工程, 2019, 48(2): 222002-0222002(8). doi: 10.3788/IRLA201948.0222002
    [15] 柯熙政, 宁川, 王姣.  大气湍流下轨道角动量复用态串扰分析 . 红外与激光工程, 2018, 47(11): 1122002-1122002(7). doi: 10.3788/IRLA201847.1122002
    [16] 黄应清, 苏健, 陈祎贝, 闫兴鹏, 蒋晓瑜.  全息打印技术研究进展 . 红外与激光工程, 2018, 47(4): 406008-0406008(13). doi: 10.3788/IRLA201847.0406008
    [17] 柯熙政, 薛璞.  轨道角动量叠加态的产生及其检验 . 红外与激光工程, 2018, 47(4): 417007-0417007(6). doi: 10.3788/IRLA201847.0417007
    [18] 步敬, 张莉超, 豆秀婕, 杨勇, 张聿全, 闵长俊.  任意拓扑荷光学旋涡的产生及应用 . 红外与激光工程, 2017, 46(6): 634001-0634001(5). doi: 10.3788/IRLA201746.0634001
    [19] 柯熙政, 郭新龙.  大气斜程传输中高阶贝塞尔高斯光束轨道角动量的研究 . 红外与激光工程, 2015, 44(12): 3744-3749.
    [20] 黄民双, 黄军芬.  布里渊光子晶体光纤环形移频器 . 红外与激光工程, 2012, 41(8): 2125-2129.
  • 加载中
图(12) / 表(3)
计量
  • 文章访问数:  131
  • HTML全文浏览量:  29
  • PDF下载量:  40
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-21
  • 修回日期:  2023-06-20
  • 刊出日期:  2023-07-25

基于轨道角动量全息和频移的大容量光学信息加密技术

doi: 10.3788/IRLA20230313
    作者简介:

    郑伟,男,讲师,博士,主要从事光学信息处理和太赫兹功能器件方面的研究

    通讯作者: 席思星,男,副教授,博士,主要从事光场调控和光学信息安全方面的研究。; 王晓雷,男,教授,博士,主要从事光场调控和太赫兹技术方面的研究。
基金项目:  国家自然科学基金项目(11904073,61875093);河北省高等学校科学技术研究项目(BJ2020028);邯郸市科学技术研究与发展计划项目(19422031008-3);山西省科技创新项目(2019L0877,2022L476);来晋科研专项项目(QZX-2018012);运城学院博士科研启动项目(YQ-2017015);山西省研究生教育教学改革课题(2022YJJG263);运城市科技计划项目(YCKJ-202268);运城学院学科建设项目
  • 中图分类号: O438.1; TP391

摘要: 提出了一种基于轨道角动量全息(Orbital Angular Momentum, OAM)和频移的大容量光学信息加密方法。该方法实现了对多个图像信息的并行加密。首先,对多幅原始图像进行采样,采样阵列的采样间隔取决于具有不同拓扑荷数的螺旋相位的空间频率。然后,多个采样图像信息经过随机相位调制、傅里叶变换和频移相位调制后相干叠加构成轨道角动量保留全息图。最后,将不同拓扑荷的螺旋相位分别编码到轨道角动量保留全息图中,得到轨道角动量选择全息图,进行相干叠加后构成最终的单个加密全息图。解密时,轨道角动量复合选择全息图被加载到空间光调制器上,用包含特定拓扑荷数的涡旋光束照射,并经过傅里叶变换获得多个解密信息。该加密系统具有极高的加密灵活性和极大的加密容量,不仅可以在同一拓扑荷下,设计不同的频移因子来并行加密一组多个图像信息,还可以利用不同拓扑荷对多组图像信息进行加密。该方法将涡旋光束的模式设定为一个新的光学密钥,极大地提高了光学加密系统的安全性。此外,该光学加密方法中,待加密图像信息的尺寸不受空间光调制器的像元数量限制,极大地提高了光学实现信息加密的可行性和有效性。仿真实验结果表明该方法具有较高的安全性、抗噪性和抗剪切能力。

English Abstract

    • 涡旋光束是一种具有独特波前结构的光束。涡旋光束携带的轨道角动量(Orbital Angular Momentum, OAM)与光的振幅、相位、偏振和波长等参量相互独立,并且模式不同的OAM之间具有正交性,故而可以无限叠加[1]。因此,轨道角动量在光学加密[2]、光学通信[3-7]、全息成像[8]、光操纵[9]、超分辨成像[10-12]、高精度测量[13-15]等领域具备巨大的应用潜力,而将OAM应用于图像加密领域有望为高安全、超大容量加密系统的研究开辟出一条道路。

      虽然OAM是一个独立的光学参量,但由于布拉格衍射公式中缺乏拓扑荷数的选择性,使得传统的数字全息图并不能显示OAM的选择性,即全息图在OAM光束照明的情况下,所有的OAM模式都以相同的方式衍射,这就限制了以OAM作为密钥的光学加密系统的发展。为了解决这个问题,2019年12月,Fang等人[16]提出一种保留OAM选择特性全息图的获取方法。该方法通过编码一个OAM特有的采样阵列将OAM特性保存到全息图像上。2020年1月,Ren等人[17-20]在此基础上提出了一种基于超表面编码的轨道角动量加密全息成像方法,该方法通过使用超高维OAM复用的大尺度复振幅超表面演示了多达200个正交图像通道的无透镜重建和全息视频显示,该方法验证了OAM大容量加密系统的可行性。为了提高系统的安全性,同年5月,Zhou等人[21]提出了偏振加密轨道角动量多路复用超表面全息术。该方法使用双折射超表面并通过轨道角动量复用和不同偏振通道进行全息加密,只有利用正确的拓扑荷和特定的偏振态才能重建OAM选择全息图的信息。2021年,该课题组[22]又提出了高维轨道角动量复用非线性全息术,成功地解决了非线性光学系统中非线性谐波场与输入场耦合的问题。

      为了进一步提高轨道角动量加密系统的容量和突破空间光调制器的像元数量限制,文中提出了一种基于轨道角动量全息术和频移的大容量光学信息加密新方法。该方法在利用不同拓扑荷进行加密的基础上,结合频域复用技术实现了在不同拓扑荷下对多组不同信息的加密。多组待加密图像信息经过空间频率采样、随机相位调制、傅里叶变换和频移相位调制后相干叠加得到多个轨道角动量保存全息图,再将不同拓扑荷的螺旋相位编码到这些轨道角动量保存全息图上,从而得到多个轨道角动量选择全息图,将多个不同拓扑荷对应的轨道角动量选择全息图相干叠加构成一幅轨道角动量复用选择全息图。多个解密信息只有在拓扑荷全部正确时才能通过解密系统获取。该系统将多个信息加密成单幅图像,并设置轨道角动量的拓扑荷为密钥,极大地提高了光学信息加密系统的安全性和加密容量。但由于实验设备的影响,轨道角动量所能提供的容量是有限的,因此结合频域复用技术能够极大地提升加密系统的容量。通过仿真实验验证了该大容量光学信息加密方案的有效性和可行性。实验结果表明该方法具有很高的安全性、抗噪性和抗剪切能力。此外,该方法在一定程度上可以突破SLM像元数量的限制,实现对超越SLM像元数量的大尺寸图像信息的加密。

    • 该加密方法可以对多幅不同尺寸和类型的图像信息实现并行加密。为了验证该方法对多幅图像信息的加密性能,文中选取的八幅待加密图像信息如图1所示,分别为512×512个像素的“蛙”和“熊猫”,512×256个像素的“树”和“棋盘”,256×512个像素的“optics”和“光”,256×256个像素的“像”字和“花”,将八幅图像信息分为左右两组进行加密实验。

      图  1  原始图像信息

      Figure 1.  Original image information

      首先,对原始图像信息进行采样,得到对应的采样图像信息。其中,采样阵列的采样常数d由不同拓扑荷数的螺旋相位板的空间频率来确定。在圆柱坐标下,螺旋相位板的分布可以表示为:

      $$ \psi (r,\theta ) = {\text{circ}}\left(\frac{r}{R}\right)\exp ({{i}}m\theta ) $$ (1)

      式中:$R$为螺旋相位板的半径;$r$和$\theta $分别为极坐标系中的半径和方位角。螺旋相位板的空间频谱的振幅模量为:

      $$ \left| {E(\rho ,\varphi )} \right| = \left| {{\text{FT\{ }}\psi (r,\theta )\} } \right| = \left| {\frac{{{{( - 1)}^{m + 1}}k}}{f}\exp ({{i}}m\varphi ) \int\limits_0^R {r{J_m}\left(\frac{k}{f}r\rho \right){\text{d}}r} } \right| $$ (2)

      式中:$k = {{2{\pi}} \mathord{\left/ {\vphantom {{2{\pi}} \lambda}} \right. } {\lambda}}$和$f$分别为入射光的波数和傅里叶透镜的焦距;FT表示傅里叶变换;$\;\rho $和$\varphi $分别为全息图平面上的半径和方位角;${J_m}$表示第一类贝塞尔函数。由于振幅模量具有径向对称性,当方位角为零时,其分布为:

      $$ \left| {E(\rho ,0)} \right| = \left| {\frac{{{{( - 1)}^{m + 1}}k}}{f}\int\limits_0^R {r{J_m}\left(\frac{k}{f}r\rho \right){\text{d}}r} } \right| $$ (3)

      故在近轴极限下,不同拓扑荷数的螺旋相位板的空间频率可以用基于傅里叶变换的圆环状强度分布来表示,如图2所示。其中,采样常数d与不同拓扑荷下螺旋相位板的空间频率相对应[16],其被定义为最大振幅模量30%时所对应的像素宽度。在计算过程中,激光的波长为$ \lambda $,SLM的单个像素尺寸为$ h $。因此通过以下公式计算得到全息图的最大数值孔径:

      图  2  螺旋相位板空间频率及采样阵列的示意图

      Figure 2.  Schematic diagram of spatial frequency and sampling array of spiral phase plates

      $$ {\text{sin}}\theta {\text{ = }}\frac{\lambda }{{2h}} $$ (4)

      为了将全息图所携带的全部信息进行光学传输,可计算重建全息图像的傅里叶透镜的最小的数值孔径取值,然后通过计算最大模振幅的30%所对应的像素宽度就可以得出采样常数d的取值。拓扑荷$m = 1,2, 3,4$时采样常数d大小如图3所示。

      图  3  不同拓扑荷下采样常数的示意图

      Figure 3.  Sampling constants diagram of under different topological charges

      然后,利用傅里叶变换的频移特性将8个采样图像信息合成为两幅保留全息图,其获取过程示意图如图4所示。采样图像信息分别经过随机相位调制、逆傅里叶变换和频移相位调制后相干叠加构成两个OAM保留全息图$ {F_1}(u,v) $和$ {F_2}(u,v) $。

      图  4  OAM保留全息图获取过程的示意图

      Figure 4.  Schematic diagram of OAM-preserved hologram acquisition process

      $$ {g_i}(x,y){\text{ = }}{f_i}(x,y)\exp [{{i}}2{\text{π }} \times {\text{rand}}({m_i},{n_i})] $$ (5)
      $$ {G_i}(u,v) = {\text{IFT[}}{g_i}(x,y)] $$ (6)
      $$ {F_1}(u,v) = \sum\limits_{i = 1}^4 {{G_i}(u,v)\exp [{{i}}2{\text{π }}({a_i}u + {b_i}v)} ] $$ (7)
      $$ {F_2}(u,v) = \sum\limits_{i = 5}^8 {{G_i}(u,v)\exp [{{i}}2{\text{π }}({a_i}u + {b_i}v)} ] $$ (8)

      式中:$ {f_i}(x,y) $表示采样图像信息;$\exp [{{i}}2{\text{π }} \times {\text{rand}} ({m_i},{n_i})]$为随机相位因子;$ \text{IFT} [\cdot ] $表示逆傅里叶变换;$\exp [{i}2{\text{π }} \cdot ({a_i}u + {b_i}v)]$为频移相位因子,其中$ ({a_i},{b_i}) $表征各个图像的位移系数,具体数值与坐标轴选取以及成像平面大小有关。文中,分别为“蛙”$ ({a_{\text{1}}},{b_{\text{1}}}){\text{ = }}(1{\text{/}}3,1{\text{/}}6) $,“树”$ ({a_{\text{2}}},{b_{\text{2}}}){\text{ = }}(1{\text{/}}12,{\text{1/6}}) $,“optics”$ ({a_{\text{3}}},{b_{\text{3}}}){\text{ = }}({\text{1/3}}, - {\text{1/3}}) $,“像”$ ({a_{\text{4}}},{b_{\text{4}}}){\text{ = }}(1{\text{/}}12, - {\text{1/3}}) $,“花”$ ({a_5},{b_5}){\text{ = }}( - {\text{1/12}},{\text{1/3}}) $,“光”$ ({a_6},{b_6}){\text{ = }}( - 1{\text{/}}3,{\text{1/3}}) $,“熊猫”$ ({a_7},{b_7}){\text{ = }}( - 1{\text{/}}6, - {\text{1/6}}) $,“棋盘”$ ({a_8},{b_8}){\text{ = }}( - {\text{5/12}}, - {\text{1/6}}) $。

      最后,为了在保留全息图中实现OAM选择特性,将拓扑荷$m = - 1$和$m = - 2$的螺旋相位板分别编码到两个OAM保留全息图上,这样就得到了两个OAM选择全息图,将其相干叠加构成最终的OAM复合选择全息图,即密文图像,该过程如图5所示。

      图  5  OAM复合选择全息图获取过程的示意图

      Figure 5.  Schematic diagram of OAM composite selective hologram acquisition process

      图4图5可知,相比于其他的多图像加密算法,该方法通过两次加密过程可以将多个不同尺寸和类型的待加密图像信息加密为单个全息图,并且密文图像中未显示任何原始信息。该加密系统具有极高的加密灵活性和极大的加密容量,不仅可以在同一拓扑荷下,设计不同的频移因子来加密一组多个图像信息,还可以利用不同拓扑荷对多组图像信息进行同时加密。该方法将OAM光束的特定拓扑荷设定为解密密钥,使得加密系统的安全性得到极大的提高。

    • 为了验证加密结果的有效性和鲁棒性,本节对加密结果的相关性进行了分析。图像的相关性可以表征图像相邻位置像素值的相关程度,加密图像相关性越低,加密方法的安全性越高。本节分别对原始图像信息和密文图像信息在水平、垂直和对角3个方向上的相关性系数进行了计算,计算结果如表1所示。计算公式如下:

      $$ {{\rm{cov}}} (x,y) = \frac{1}{N} \times \sum\limits_{i = 1}^N {[{x_i} - E(x)} ][{y_i} - E(y)] $$ (9)
      $$ CC = \frac{{{{\rm{cov}}} (x,y)}}{{\sqrt {D(x)D(y)} }} $$ (10)

      表1的计算结果可见,密文图像信息各个方向上的相关性都很低,文中的加密方法有效地降低了高度相关的原始图像信息的统计特性,能够有效抵抗基于像素相关性的统计攻击,具有较高的安全性。

      表 1  原始图像信息和密文图像信息的相关性系数

      Table 1.  Correlation coefficient of original image and ciphertext image

      Image
      Horizontal correlation

      Vertical correlation

      Diagonal correlation
      Frog0.99770.99790.9963
      Tree0.97150.97110.9504
      “optics”0.86610.93830.8415
      “像”0.90730.91060.8139
      Flower0.98810.99910.9977
      “光”0.98840.97910.9716
      Panda0.99760.99090.9841
      Chessboard0.98810.98610.9743
      Encrypted text−0.00130.01640.0072
    • 明文敏感性是一种评价加密系统质量的常用指标,明文敏感性越好,加密系统的安全性越好。可以通过计算某一幅明文图像发生微小变化(某一像素值增大0.1和某两个像素值交换)时密文图像的像素值改变率(NPCR)和归一化平均变化强度(UACI),以此来验证文中加密方法的明文敏感性,计算结果如表2所示。计算公式如下:

      $$ NPCR = \frac{{\displaystyle\sum\limits_i {\displaystyle\sum\limits_j {z(i,j)} } }}{{M \times N}} $$ (11)
      $$ UACI = \frac{1}{{M \times N}}\left[\sum\limits_i {\sum\limits_j {\left| {I(i,j) - {I'}(i,j)} \right|} } \right] $$ (12)

      式中:当$I(i,j){\text{ = }}{I'}(i,j)$时,$z(i,j){\text{ = 0}}$;当$I(i,j) \ne {I'}(i,j)$时,$z(i,j){\text{ = 1}}$。

      表 2  明文敏感性分析

      Table 2.  Plaintext sensitivity analysis

      Index change when pixel
      increases by 0.1
      Index changes when
      two pixels are exchanged
      NPCRUACINPCRUACI
      99.85%34.75%99.84%34.58%

      表2可知,在该加密系统中当某一明文图像信息发生上述两种的微小变化时密文像素值改变率均在99%以上,其归一化平均改变强度均在30%以上,这说明该系统具有较好的明文敏感性。

    • 多个图像信息的解密可以通过使用特定拓扑荷数的OAM光束和解密光路来实现。文中解密光路设置如图6所示,不同拓扑荷数的OAM光束由SLM1产生,SLM2上加载复合全息图,傅里叶透镜将频域中叠加的各个图像按照对应的频移因子进行分离,最后各个解密图像由CCD在焦平面上接收。

      图  6  解密光路示意图

      Figure 6.  Schematic diagram of the decryption optical path

      由于密文全息图存在OAM选择性,因此只有在包含$m = 1$和$m = 2$模式的OAM光束照射时才能解密出所有图像信息。通过使用拓扑荷为$m = 1$和$m = 2$的OAM光束以及解密光路进行多个图像信息的解密,其解密结果如图7所示。

      图  7  当$m{\text{ = 1,\;2}}$的OAM光束照射时的解密结果

      Figure 7.  Decryption result when $m{\text{ = 1,\;2}}$ OAM beam is irradiated

      将解密结果和原始图像进行对比,采用结构相似性指数SSIM和相关性系数CC来描述两图像间的相似度,SSIM和CC的计算结果如表3所示。结构相似性 (Structural SIMilarity, SSIM),衡量两张图像的相似度,通过SSIM指标可以衡量样本和预测图像之间的亮度、对比度和结构,SSIM指标更加关注图像生成的感知面的图像质量。

      表 3  原始图像和解密图像的相似性

      Table 3.  Similarity of original image and decrypted image

      Similarity indexFrogTree“optics”“像”Flower“光”PandaChessboard
      SSIM0.97030.96050.96450.95090.98290.97570.96990.9742
      CC0.99430.99320.99100.99240.99400.99220.99540.9936
      $$ I(x,y) = \frac{{2{\mu _x}{\mu _y} + {c_1}}}{{{\mu _x}^2 + {\mu ^2}_y + {c_1}}} $$ (13)
      $$ c(x,y) = \frac{{2{\sigma _x}{\sigma _y} + {c_2}}}{{{\sigma _x}^2 + {\sigma ^2}_y + {c_2}}} $$ (14)
      $$ s(x,y) = \frac{{{\sigma _{xy}} + {c_3}}}{{{\sigma _x}{\sigma _y} + {c_3}}} $$ (15)

      公式(13)表示的是图像亮度的相似程度,其中$\; {\mu _x} $以及$ \;{\mu _y} $分别表示图像的$ x $和$ y $的均值;公式(14)表示的所示图像的对比度相似性,其中$ {\sigma _x} $和$ {\sigma _y} $分别表示的是图像$ x $和图像$ y $的方差;而公式(15)表示图像的结构相似度,其中$ {\sigma _{xy}} $表示图像$ x $和图像$ y $的协方差,$ {c_1} $、$ {c_2} $、$ {c_3} $均为常数。将上述的公式联立,可以得到公式(16),SSIM的取值在0~1之间,SSIM的值越大,说明两种图像更相似。

      $$ {\rm{SSIM}}(x,y) = \frac{{2({\mu _x}{\mu _y} + {c_1})(2{\sigma _{xy}} + {c_2})}}{{({\mu ^2}_x{\mu ^2}_y + {c_1})({\sigma ^2}_x + {\sigma ^2}_y + {c_2})}} $$ (16)
    • 设置3个对照试验来对密钥的有效性进行验证,分别是仅使用拓扑荷为$m = 1$的OAM光束照射、仅使用拓扑荷为$m = 2$的OAM光束照射以及使用拓扑荷不为$m = 1,2$的OAM光束照射,将实验结果与图7进行对比分析来对密钥的有效性进行验证。其结果如下:

      当仅使用拓扑荷为$m = 1$的OAM光束照射时只能解密出图像“蛙”、“树”、“optics”、“像”,如图8所示;当仅使用拓扑荷为$m = 2$的OAM光束照射时只能解密出图像“花”、“光”、“熊猫”、“棋盘”,如图9所示;当使用不包含$m = 1$和$m = 2$任一模式的OAM光束照射时不能解密出任何图像,解密结果如图10所示。

      图  8  当$m{\text{ = 1}}$的OAM光束照射时的解密结果

      Figure 8.  Decryption result when $m{\text{ = 1}}$ OAM beam is irradiated

      图  9  当$m{\text{ = 2}}$的OAM光束照射时的解密结果

      Figure 9.  Decryption result when $m{\text{ = 2}}$ OAM beam is irradiated

      图  10  当$m \ne {\text{1,\;2}}$的OAM光束照射时的解密结果

      Figure 10.  Decryption result when $m \ne {\text{1,\;2}}$ OAM beam is irradiated

      图7~10可知,只有在使用正确的密钥时才可以解密出对应的多个图像,当密钥错误时,其无法获取任何原始图像信息。因此,该加密方法中以拓扑荷数$m$作为密钥具有很好的有效性,从而使得该加密系统具有很高的安全性。

    • 当密文图像受到不同面积的剪切后使用正确密钥进行解密,通过解密结果质量的优劣来对该加密系统的抗剪切性进行分析。从密文图像的边缘进行圆环剪切,图11为解密图像信息与原始信息的平均结构相似性指数SSIM随密文丢失像素面积比S%变化的关系曲线图,并分别给出了密文丢失像素面积比为40%、55%和70%时的解密结果。

      图  11  抗剪切性分析结果

      Figure 11.  Shear resistance simulation results

      图11可以看出,在该加密方法中,当密文丢失像素面积比小于50%时,其仍可以通过正确的密钥和解密光路解密出多个尺寸和类型不同的图像,说明该系统具有较好的抗剪切攻击能力。其中由于二值图的像素值较为简单,其抗剪切能力较灰度图像更强。

      通过对密文图像添加不同密度的椒盐噪声并使用正确密钥进行解密,观察分析解密结果质量来对该加密系统的抗噪性进行验证。图12为解密图像与原始图像的平均结构相似性指数SSIM随密文图像中加入的椒盐噪声密度变化的关系曲线,并分别给出了添加椒盐噪声密度为10%、30%和50%时的解密结果。

      图  12  抗噪性分析结果

      Figure 12.  Noise immunity simulation results

      图12可以看出,在该加密方法中,当密文图像加入椒盐噪声的密度小于30%时,其可以通过正确的密钥和解密光路解密出多个质量较好的图像,说明该系统具有较好的抗噪攻击能力。其中由于二值图的像素值较为简单,其抗噪能力略强于灰度图像。同时,对该系统进行了多种暴力攻击测试,如:滤波攻击、高斯噪声攻击等均进行了攻击测试,结果表明该系统对于抵抗这些暴力攻击均有着不错的表现。

    • 为改善加密系统容量,文中提出了一种基于轨道角动量全息和频移的大容量光学信息加密方法。多个待加密图像信息经过两次加密调制后构成一幅光学全息图,解密时由特定拓扑荷数的涡旋光束通过解密系统获取多个解密图像信息。该方法在保留OAM不同拓扑荷正交性的基础上,结合傅里叶变换的频移特性实现了在不同拓扑荷下对多组不同图像信息的加密。该系统将多个图像信息加密成单个全息图,并设置轨道角动量的拓扑荷作为密钥,极大地提高了光学图像加密系统的安全性和加密容量。此外,待加密图像信息的尺寸不受空间光调制器的像元数量限制,极大地提高了光学实现信息加密的可行性和有效性。

参考文献 (22)

目录

    /

    返回文章
    返回