混合混沌虚拟光学图像加密方法.pdf

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910181636.3 (22)申请日 2019.03.11 (71)申请人 齐齐哈尔大学 地址 161006 黑龙江省齐齐哈尔市文化大 街42号 (72)发明人 郭媛魏连锁许鑫敬世伟 (74)专利代理机构 哈尔滨龙科专利代理有限公 司 23206 代理人 高媛 (51)Int.Cl. H04N 1/44(2006.01) H04N 1/32(2006.01) H04L 9/00(2006.01) (54)发明名称 一种混合混沌虚拟光学图像加密方法 (57)摘要 本发明公开

2、了一种混合混沌虚拟光学图像 加密方法， 所述方法包括如下步骤： 步骤一、 假设 明文图像的大小为MN， 从左往右依次读取图像 各点的像素值， 得到其像素矩阵F(x,y)； 步骤二、 对明文图像F(x,y)进行相位编码， 得到Q(x,y)； 步骤二、 Q(x,y)经过随机相位模板RM1调制， 得到 F1(x,y)； 步骤三、 F1(x,y)经傅里叶变换后经过 随机相位模板RM2调制， 得到F2(， )； 步骤四、 F2(， )经过傅里叶逆变换， 得到密文图像U1； 步骤五、 对U1(x,y)进行像素值替代的二次加密， 得到密文图 像U2。 本发明采 用共参数广义 Fibonacci混沌系统和耦合

3、帐篷时空混沌系统分 别构造均匀性好的随机模板， 使两个混沌系统具 有联动效应， 弥补了Logistic混沌系统构造随机 模板元素分布不均匀的缺点， 增强对密文的隐藏 性及密钥的敏感性。 权利要求书2页 说明书8页 附图5页 CN 109981921 A 2019.07.05 CN 109981921 A 1.一种混合混沌虚拟光学图像加密方法， 其特征在于所述方法包括如下步骤： 步骤一、 假设明文图像的大小为MN， 从左往右依次读取图像各点的像素值， 得到其像 素矩阵F(x,y)， 其中： M、 N分别代表明文图像平面矩阵的行数和列数， (x,y)代表明文图像像 素点在空域的位置坐标； 步骤二、

4、 对明文图像F(x,y)进行相位编码， 得到相位编码后的图像Q(x,y)： Q(x,y)expiF(x,y)； 式中， i代表虚数单位； 步骤二、 采用广义Fibonacci混沌系统构造随机相位模板RM1， Q(x,y)经过随机相位模 板RM1调制， 得到经过随机模板RM1调制后的图像F1(x,y)： F1(x,y)Q(x,y)RM1expiF(x,y)expi2 a(, ,x1, 1)； 式中， RM1expi2 a(, ,x1, 1)； 、 、 1、 x1分别为构造随机相位模板RM1的广义 Fibonacci混沌系统的可调参数、 耗散参数、 参数、 初值； 步骤三、 采用耦合帐篷时空混沌系

5、统构造随机相位模板RM2， F1(x,y)经傅里叶变换后经 过随机相位模板RM2调制， 得到经过傅里叶变换和随机模板RM2调制后的图像F2( ， )： 式中， RM2expi2 b( , ,x1, 1)； 、 分别为构造随机相位模板RM2的耦合帐篷时空 混沌系统的耦合系数、 帐篷映射参数； ( ， )代表像素点在频域的位置坐标； 步骤四、 F2( ， )经过傅里叶逆变换， 得到密文图像U1： 式中， FT代表傅里叶变换； FT-1代表傅里叶逆变换； 步骤五、 采用广义Fibonacci混沌系统构造随机相位模板RM3， 对U1(x,y)进行像素值替 代的二次加密， 得到密文图像U2： U2(x,

6、y)newampexp(inewpha)； 式中， newamp表示新的振幅； newpha表示新的相位； RM3c(, ,x2, 2)； x2x0+|mean (U1)|mod1；2、 x0分别为构造随机相位模板RM3的广义Fibonacci混沌系统的参数、 初值。 2.根据权利要求1所述的混合混沌虚拟光学图像加密方法， 其特征在于所述newamp的 计算公式如下： 式中， amp代表取振幅操作。 3.根据权利要求1所述的混合混沌虚拟光学图像加密方法， 其特征在于所述newamp的 计算公式如下： 权利要求书 1/2 页 2 CN 109981921 A 2 式中， 为权重， pha代表取相

7、位操作。 4.根据权利要求1、 2或3所述的混合混沌虚拟光学图像加密方法， 其特征在于所述 (3.74,4)。 5.根据权利要求1、 2或3所述的混合混沌虚拟光学图像加密方法， 其特征在于所述 3.5。 6.根据权利要求1所述的混合混沌虚拟光学图像加密方法， 其特征在于所述x0(0,1)。 7.根据权利要求1、 2或3所述的混合混沌虚拟光学图像加密方法， 其特征在于所述 (0,1)。 8.根据权利要求1、 2或3所述的混合混沌虚拟光学图像加密方法， 其特征在于所述 (0,1)。 权利要求书 2/2 页 3 CN 109981921 A 3 一种混合混沌虚拟光学图像加密方法 技术领域 0001

8、本发明涉及一种虚拟光学图像加密方法， 具体涉及一种混合混沌虚拟光学图像加 密方法。 背景技术 0002 图像作为一种载体， 包括大量信息， 其安全性非常重要。 虚拟光学图像加密 (Virtual optical image encryption,VOIE)既具有光学方法的高维空间、 多自由度等优 点， 又具备计算机信息安全技术的优点， 成为广大学者研究的热门课题。 1995年Refregier 和Javidi提出4f加密系统， 最具有代表性。 该系统利用两块随机模板作为密钥， 增大密钥空 间。 2007年Peng指出该系统对于已知明文， 选择明文和唯密文攻击的脆弱性。 随后该方法又 被推广到基

9、于分数傅里叶变换， 基于菲涅耳变换， 基于计算全息等加密系统。 以上研究中随 机模板虽然增大了密钥空间， 但同时也导致密钥体积过大， 其分发、 传递困难。 由于混沌系 统的伪随机性与VOIE系统中随机模板的要求一致， 许多学者提出基于混沌系统的虚拟光学 图像加密方法， 利用混沌系统的初值和参数替代随机模板作为密钥， 减小密钥体积， 又保留 光学加密快速并行计算的优点。 2013年Song等人提出基于混沌映射和4f的图像加密算法研 究， 提出了基于混沌的应用于灰度图像和彩色图像加密的新型算法。 2014年Zhu等人提出基 于混沌的改进双随机相位编码图像加密算法， 基于菲涅耳域的双随机相位编码提出

10、一种改 进算法。 混沌模型中非线性函数采用实数域Logistic映射， 具有不均匀分布特性， 存在安全 隐患。 发明内容 0003 针对现有技术存在的上述问题， 本发明采用共参数的两种混沌系统生成可置换传 统4f系统随机相位模板的新模板， 提出了一种混合混沌虚拟光学图像加密方法。 本发明采 用共参数广义Fibonacci混沌系统和耦合帐篷时空混沌系统分别构造均匀性好的随机模 板， 使两个混沌系统具有联动效应， 弥补了Logistic混沌系统构造随机模板元素分布不均 匀的缺点， 增强对密文的隐藏性及密钥的敏感性。 将明文图像进行相位编码， 克服系统解密 不需要空域密钥的缺陷， 进行虚拟光学图像加

11、密。 利用广义Fibonacci混沌系统构造第三块 随机模板， 进行振幅和相位像素值替代的二次加密， 更好的抵御统计分析， 并将第一次密文 像素值与混沌系统初值相联系， 提高密文对明文的敏感性和抵御选择明文攻击的能力。 0004 本发明的目的是通过以下技术方案实现的： 0005 一种混合混沌虚拟光学图像加密方法， 包括如下步骤： 0006 步骤一、 假设明文图像的大小为MN， 从左往右依次读取图像各点的像素值， 得到 其像素矩阵F(x,y)， 其中： M、 N分别代表明文图像平面矩阵的行数和列数， (x,y)代表明文图 像像素点在空域的位置坐标； 0007 步骤二、 对明文图像F(x,y)进行

12、相位编码， 得到相位编码后的图像 Q(x,y)： 0008 Q(x,y)expiF(x,y)； 说明书 1/8 页 4 CN 109981921 A 4 0009 式中， i代表虚数单位； 0010 步骤二、 采用广义Fibonacci混沌系统构造随机相位模板RM1， Q(x,y)经过随机相 位模板RM1调制， 得到经过随机模板RM1调制后的图像F1(x,y)： 0011 F1(x,y)Q(x,y)RM1expiF(x,y)expi2 a(, ,x1, 1)； 0012 式中， RM1expi2 a(, ,x1, 1)； Key1(, ,x1, 1)， Key1为密钥； 、 、 1、 x1分别

13、为构造随机相位模板RM1的广义Fibonacci混沌系统的可调参数、 耗散参数、 参数、 初 值； 0013 步骤三、 采用耦合帐篷时空混沌系统构造随机相位模板RM2， F1(x,y)经傅里叶变 换后经过随机相位模板RM2调制， 得到经过傅里叶变换和随机模板RM2调制后的图像： 0014 0015 式中， RM2expi2 b( , ,x1, 1)； Key2( , ,x1, 1)， Key2为密钥； 、 分别为 构造随机相位模板RM2的耦合帐篷时空混沌系统的耦合系数、 帐篷映射参数； ( ， )代表像 素点在频域的位置坐标； 0016 步骤四、 F2( ， )经过傅里叶逆变换， 得到密文图像

14、U1： 0017 0018 式中， FT代表傅里叶变换； FT-1代表傅里叶逆变换； 0019 步骤五、 采用广义Fibonacci混沌系统构造随机相位模板RM3， 对U1(x,y) 进行像 素值替代的二次加密， 得到密文图像U2： 0020 U2(x,y)newampexp(inewpha)； 0021 0022 0023 式中， newamp表示新的振幅； newpha表示新的相位； RM3c(, ,x2, 2)； Key3 (, ,x2, 2)； x2x0+|mean(U1)|mod1， x0(0,1)； 为权重； Key3为密钥；2、 x0分别为构 造随机相位模板RM3的广义Fibon

15、acci混沌系统的参数、 初值； Key1、 Key2、 Key3为主要密钥， Key4， Key4为辅助密钥； amp代表取振幅操作； pha代表取相位操作； x2代表广义 Fibonacci混沌系统的初值x0， 与一次加密密文U1像素平均值的绝对值， 二者相加后对1取 模得到的值。 0024 相比于现有技术， 本发明具有如下优点： 0025 1、 本发明对明文图像进行相位编码， 构建可产生均匀非相关随机序列的广义 Fibonacci混沌系统， 克服原4f系统对第一块相位模板不敏感的缺陷及Logistic混沌系统 说明书 2/8 页 5 CN 109981921 A 5 构造随机模板元素分布

16、不均匀的缺点， 提高密钥传输效率及系统对密钥的敏感性。 0026 2、 本发明对一次加密得到的复值图像进行二次加密， 通过提取其振幅及相位进行 像素值替代操作， 解决复值图像像素值不能按位 “异或” 进行替代的问题， 使密文图像分布 更均匀。 0027 3、 本发明在二次加密中将产生第三块随机模板的混沌初值与一次密文联系， 增强 系统对明文的敏感性， 有效抵御选择明文攻击。 0028 4、 仿真实验表明， 本发明的方法有效增加了密钥空间和密钥敏感性， 提高了加密 系统加密效率和安全性。 附图说明 0029 图1为Logistic混沌系统和广义Fibonacci混沌系统序列分布图， (a) Lo

17、gistic混 沌系统， (b)广义Fibonacci混沌系统； 0030 图2为系统加密原理； 0031 图3为加密过程框图； 0032 图4为算法的加解密图像， (a)明文图像， (b)密文图像， (c)解密图像； 0033 图5为灰度直方图， (a)明文图像， (b)密文图像U1， (c)密文图像U2； 0034 图6为lena图的垂直方向相邻像素分布图， (a)明文图像， (b)密文图像U1， (c)密文 图像U2； 0035 图7为密钥偏差时的解密图像， (a)10-15， (b)10-16， (c) 10-18， (d) 10-19， (e) 210-15， (f) 210-16。

18、 具体实施方式 0036 下面结合附图对本发明的技术方案作进一步的说明， 但并不局限于此， 凡是对本 发明技术方案进行修改或者等同替换， 而不脱离本发明技术方案的精神和范围， 均应涵盖 在本发明的保护范围中。 0037 本发明提供了一种混合混沌虚拟光学图像加密方法， 所述方法包括如下步骤： 0038 一、 构造混沌系统 0039 (1)广义Fibonacci混沌系统 0040 本发明构造一个可以产生均匀非相关随机序列的广义Fibonacci混沌系统， 其函 数模型如下： 0041 Fi(AFi-1+BFi-2+CFi-3)mod M (1)； 0042 Fi、 Fi-1、 Fi-2、 Fi-3

19、分别代表函数模型在第i、 i-1、 i-2、 i-3个数时的状态值； i为指示 量， i1,2,3； A、 B、 C分别为三个变量； mod表示取模运算； M表示模数。 0043 式中， A、 B、 C分别为选取的量子Logistic混沌系统所产生变量xi、 yi和zi的序列组 合， 增加了序列的随机性， 降低相关性； 量子Logistic混沌系统表达式为： 0044 说明书 3/8 页 6 CN 109981921 A 6 0045 式中， 为可调参数， (3.74,4)； 为耗散参数， 3.5； xi、 yi和zi为系统第i次 迭代的状态值； xi+1、 yi+1和zi+1为系统第i+1次

20、迭代的状态值； i为指示量， i1,2,3； 分别是xi和zi的复共轭。 0046 经广义变参数三阶Fibonacci的函数模型处理后， 再对Logistic混沌系统xi+1 1xi(1-xi)进行扰动得到广义Fibonacci混沌系统： 0047 CFQL(F(Q(, )+L(x1, 1)mod1 (3)； 0048 式中， Q(, )表示量子Logistic混沌系统， 、 分别是量子Logistic 混沌系统 的可调参数和耗散参数； F(Q(, )表示Fibonacci系统与量子 Logistic混沌级联； L(x1, 1)表示Logistic混沌系统， x1、 1分别是Logistic

21、混沌系统的初值和参数。 mod表示取模运 算。 不断迭代式(3)产生均匀非相关混沌序列。 为比较Logistic混沌系统与广义Fibonacci 混沌系统的序列分布情况， 迭代500次的序列如图1所示。 从图1可以看出， 广义Fibonacci 混沌系统比Logistic混沌系统分布更加均匀， 弥补了Logistic混沌系统序列分布不均匀的 缺点。 0049 (2)耦合帐篷时空混沌系统 0050 时空混沌是一种在时间和空间方向上都具有混沌行为的非线性动力系统， 该二维 系统有更好的混沌特性。 耦合映像格子CML的模型的定义如下： 0051 0052式中， 为耦合系数；分别表示第i-1、 i、

22、i+1个格子在n 时刻的状 态， i1,2,， M,M为明文图像平面矩阵的行数， n1， 2， ， N,N为明文图像平面矩阵的列 数；分别表示帐篷混沌系统， 为系统第i-1、 i、 i+1个格子在n时刻的 局部状态演化方程。 0053 采用帐篷混沌映射作为CML测局部状态演化方程如下： 0054 0055 式中， 为帐篷映射参数， (0,1)， 当 0.5时得到的是标准帐篷映射； xi-1、 xi 分别表示第i-1， i个格子的状态， i1,2,， M,M为明文图像平面矩阵的行数； f(xi)表示帐 篷混沌系统， 为系统第i个格子的局部状态演化方程； 耦合系数 (0,1)， 当 取值小于0.0

23、1 时， 耦合帐篷时空混沌系统可以产生均匀分布的混沌序列。 将Logistic混沌系统产生的序 列xi作为系统的驱动序列 0056 二、 算法实现 0057 (1)加密原理 0058 本发明针对双随机相位编码存在的问题加以改进， 进行虚拟光学加密。 系统加密 原理如图2所示。 0059 采用广义Fibonacci混沌系统构造随机相位模板RM1， 其中： RM1expi2 a(, , x1, 1)， Key1(, ,x1, 1)。 采用耦合帐篷时空混沌系统构造随机相位模板RM2， 其中： 说明书 4/8 页 7 CN 109981921 A 7 RM2expi2 b( , ,x1, 1)， Ke

24、y2( , ,x1, 1)。 采用广义Fibonacci混沌系统构造随机 模板RM3， 其中： RM3c(, ,x2, 2)， Key3(, ,x2, 2)， x2x0+|mean(U1)|mod1， x0 (0,1)， U1为一次加密的密文。 a、 b、 c分别代表0至1的随机矩阵。 并在像素值替代过程中引 入权重 。 Key1、 Key2、 Key3分别代表三个混沌系统密钥， Key4 。 本发明加密过程如图3所 示。 0060 (2)加密步骤 0061 假设明文图像的大小为MN(M、 N分别代表明文图像平面矩阵的行数和列数)， 从 左往右依次读取图像各点的像素值， 得到其像素矩阵F(x,

25、y)。 加密算法的具体步骤如下： 0062 (a)对明文图像F(x,y)进行相位编码， 得到： 0063 Q(x,y)expiF(x,y) (6)； 0064 式中， Q(x,y)代表相位编码后的图像。 0065 (b)Q(x,y)经过随机模板RM1调制， 得到： 0066 F1(x,y)Q(x,y)RM1expiF(x,y)expi2 a(, ,x1, 1) (7)； 0067 式中， F1(x,y)代表经过随机模板RM1调制后的图像。 0068 (c)F1(x,y)经傅里叶变换后经过随机模板RM2调制， 得到： 0069 0070 式中， F2( ， )代表经过傅里叶变换和随机模板RM2调

26、制后的图像。 0071 (d)F2( ， )经过傅里叶逆变换， 得到密文图像U1： 0072 0073 式中， FT代表傅里叶变换； FT-1代表傅里叶逆变换； 0074 (e)对U1(x,y)进行像素值替代的二次加密， 得到密文图像U2： 0075 U2(x,y)newampexp(inewpha) (10)； 0076 式中， newamp表示新的振幅， 由一次加密的振幅信息与随机模板RM3异或运算得 到； newpha表示新的相位， 由一次加密的相位信息与随机模板RM3 相加平移运算得到， 并引 入权重因子 ， 调节密文图像U1的相位信息和随机模板RM3所占比例的大小， 进一步提高系 统

27、对密钥的敏感性； amp代表取振幅操作； pha代表取相位操作。 0077 0078 说明书 5/8 页 8 CN 109981921 A 8 0079 三、 实验仿真及算法性能分析 0080 选取经典lena(256256)灰度图像作为实验测试仿真对象。 图像的加密算法在 matlab2014a环境下进行实验仿真。 实验用到的密钥数据如下： 量子Logistic混沌系统的可 调参数3.99， 耗散参数 6； Logistic混沌系统的初值x10.4， 参数 13.99； 帐篷映 射的耦合系数 0.001， 参数 0.5； Logistic混沌系统初值x00.39， 参数 24； 权重因 子

28、0.2。 仿真结果如图4所示。 由图4可知， 密文图像已经看不出明文图像的任何轮廓， 类 似于白噪声， 说明加密的效果良好。 解密图像和明文图像在视觉上几乎看不出差别， 说明解 密效果良好。 0081 (1)灰度直方图统计特性分析 0082 直方图是反映图像灰度值分布统计特性的一个重要形式化指标。 从图5 中可明显 看出， 明文图像的灰度值分布不均匀， 呈现出明显的分布统计规律。 密文图像U1的直方图比 明文图像的直方图更加平坦。 二次加密后密文图像U2的灰度值分布非常均匀， 很好地隐藏 了明文图像的灰度信息， 留给密码分析空间很小， 从而有效地抵御了统计分析攻击。 0083 (2)相邻像素点

29、相关性分析 0084 一个好的加密方法能够显著的破坏相邻像素的相关性。 从图6中可明显看出， 明文 图像相邻两个像素点通常具有强相关性， 密文图像U2比密文图像U1的像素点更加均匀， 扩 散性更好， 有效的破坏了明文图像的相邻像素相关性。 计算明文图像和密文图像在水平、 垂 直、 对角方向的相关系数， 结果如表1。 0085 表1相邻像素点相关性计算结果表 0086 0087 从表1中可看出， 明文图像的相关系数都接近1， 相关性很高。 密文图像的相关系数 趋近于0， 说明本发明的加密算法破坏了原始图像的统计特性。 0088 (3)敏感性分析 0089 a、 密钥敏感性分析 0090 密钥敏感

30、性是指当密钥发生微小改变时， 图像发生显著的变化。 图7 (a)和(b)表 明， RM1的密钥偏差10-15时不能正确解密， 在偏差10-16时能恢复原图像； 图7 (c)和(d)表明RM2的密钥偏差 10-18时不能解密图像， 偏差 10-19时能正确解密图 像； 图7(e)和(f)表明， RM3的密钥偏差 210-15时不能正确解密， 而偏差 210-16时能 恢复原图像； 可以认为、 和 1的灵敏度分别为10-16、 10-19、 10-16。 0091 b、 明文敏感性分析 0092 明文敏感性是指当明文图像发生微小的变化， 引起密文图像发生的巨大改变。 用 像素数改变率(NPCR)和

31、归一化平均改变强度(UACI)计算加密算法对明文的敏感性。 设密文 说明书 6/8 页 9 CN 109981921 A 9 图像位于(i,j)点的像素值分别为C1(i,j)和 C2(i,j)。 若C1(i,j)C2(i,j)， D(i,j)0； 若 C1(i,j)C2(i,j)， D(i,j)1。 NPCR和 UACI定义如下： 0093 0094 0095 式中， i、 j分别代表密文图像平面矩阵中像素点的行数和列数。 0096 表2给出了明文敏感性对比数据。 0097 表2明文敏感性对比 0098 0099 0100 NPCR和UACI理想期望值分别为： NPCR99.6094， UAC

32、I33.4635。 从表2中可 看出， VOI和CBVOI算法都对明文不敏感。 本发明的算法较文献相比， 当明文发生微小变化 时， NPCR和UACI更接近于理想期望值， 说明本发明的算法对明文非常敏感。 0101 (4)密钥空间分析 0102 密钥空间是衡量加密系统抵御穷举攻击的重要指标。 本发明以量子 Logistic混 沌系统的可调参数、 耗散参数 ， Logistic的参数 1、 初值x1， 耦合帐篷时空混沌的耦合系 数 、 参数 ， Logistic混沌系统的参数 2、 初值x0、 权重 作为系统加密密钥。 本发明每个加密 密钥采用的数据都为双精度类型， 保留15位有效数字， 可得密

33、钥空间为101510151015 10151015101510151015101510135， 从安全的角度， 密钥空间21001030就能满足 较高的安全级别， 因此对于穷举攻击， 本发明算法的密钥空间是安全的。 0103 (5)信息熵 0104 信息熵是度量信息有序性的一个重要手段， 一个系统越是混乱信息熵就越高。 对 于8-bit的图像来说， 其计算方法如下： 0105 0106 式中， P(Zi)是Zi出现的概率。 8-bit图像的信息熵的最大值为8， 本发明算法的信息 熵是7.9958， 接近信息熵的最大值， 密文像素值分布均匀， 加密效果好。 0107 四、 结论 0108 本发明

34、将光学信息安全技术和混沌密码学有效结合， 提出一种新的混合混沌虚拟 光学图像加密方法。 采用混合混沌系统生成所需要的三块加密模板， 构造出可产生均匀非 说明书 7/8 页 10 CN 109981921 A 10 相关混沌序列的广义Fibonacci混沌系统， 具有随机性强、 相关性低的特点， 增强对密文的 隐藏性， 用于图像加密安全性更高； 对明文图像进行相位编码， 解决系统解密时不需空域密 钥的问题； 构造第三块随机模板进行振幅和相位像素值替代的二次加密， 使得密文的像素 值分布更加均匀， 有效抵御统计攻击。 并且第三块随机模板的混沌系统初值与一次加密密 文联系， 增强明文敏感性， 提高抗选择明文攻击的能力。 实验和仿真结果表明， 本发明的算 法密钥空间大， 敏感性强， 具有高安全性和实用性。 说明书 8/8 页 11 CN 109981921 A 11 图1 说明书附图 1/5 页 12 CN 109981921 A 12 图2 图3 图4 说明书附图 2/5 页 13 CN 109981921 A 13 图5 说明书附图 3/5 页 14 CN 109981921 A 14 图6 说明书附图 4/5 页 15 CN 109981921 A 15 图7 说明书附图 5/5 页 16 CN 109981921 A 16