基于计算全息和θ调制的彩色图像加密方法*

席思星 于娜娜 王晓雷 朱巧芬 董昭 王微 刘秀红 王华英

1) (河北工程大学数理科学与工程学院,邯郸 056038)

2) (南开大学现代光学研究所,天津 300350)

1 引 言

在双随机相位图像加密技术[1]提出后,不断出现了许多改进的加密系统,如分数傅里叶变换加密系统[2]、菲涅耳变换加密系统[3]、混沌和置乱加密系统[4],这些系统在加密时除了两个随机相位密钥以外,还有附加密钥,增强了系统的安全性和抵御攻击的能力.但这些加密系统都采用单色光照射,解密的图像都是二值图像或者灰度图像.近年来,Shi等[5]首次提出将叠层衍射成像技术应用于光学图像加密,他们在4f系统中用光学探针扫描原始灰度图像并经过随机相位密钥调制,最后记录加密图像的振幅完成加密; 解密通过叠层衍射成像的迭代运算完成.随后他们[6]又提出了一种基于可视密码和体全息光栅技术的光学图像加密方法,该方法将二值图像加密到可视密码中并将其伪装成二维码,最后通过体全息光栅曝光记录完成加密; 解密通过体全息光栅的光学再现完成.然而,彩色图像在图像细节、逼真度和色彩表现等方面是灰度图像所不能比拟的,因此彩色图像加密有巨大的研究价值和应用空间,彩色图像加密也逐渐成为国内外学者和研究人员的研究热点[7−10].当前的彩色图像加密方案主要基于空域彩色图像加密和变换域彩色图像加密,前者主要利用混沌系统直接置乱图像的像素(值)和位置,后者是对变换域系数进行加密处理.例如肖迪和谢沂均[11]采用彩色图像JPEG压缩编码提出了一种彩色图像加密方法.秦怡和郑长波[12]提出了基于傅里叶变换和双随机相位编码的彩色图像加密方案.Yuan等[13]利用棋盘光栅和衍射成像实现了彩色图像的单通道光学加密.

目前彩色图像的加密方法主要基于彩色三基色原理和彩色图像的编码方法[14],其在加密过程中把彩色图像分解为红、绿、蓝三个通道的子图像,每个子图像釆用灰度图像的加密方法.每个通道的参数和随机相位掩模都可作为密钥,只要有一条通道上的解密密钥是未知的,就无法获得正确的明文信息[15].这类方法提高了加密的安全性,但是其加密系统相对比较复杂,实现成本比较高,不利于实际应用.因此,杨晓苹等[16,17]利用彩色图像空间转化提出了基于双相位编码的单通道彩色图像加密方法; Zhou等[18]提出单通道彩色图像加密方法,将彩色图像R,G,B三个分量编码成一个灰度图像,再对灰度图像进行加密,降低了加密装置的要求和成本.本文在双随机相位图像加密和计算全息的基础上,提出基于θ调制原理的彩色图像加密技术,即对彩色图像的 R,G,B 三个分量分别进行加密并在空域叠加,其中θ调制的作用是使R,G,B三个分量在空域重叠而在频域分离,从而可以在解密过程中恢复出每个色彩分量.由于该方法只需要一个菲涅耳衍射和θ调制的加密和解密装置,因此降低了对实验装置和操作性的要求.最终的加密结果为一幅二元实值的计算全息图(computer generated hologram,CGH),不仅完全隐藏了彩色图像的信息,而且易于存储和传输.除了传统的双随机相位密钥外,该方法中的菲涅耳衍射距离和θ调制的滤波器等都是关键的密钥,并且计算全息加密图像具有高的抗剪切和抗噪能力.同时,本文采取了相位截断非对称加密系统,并进行低频滤波处理的方法,使得安全性被进一步提高[19,20].因此,本文提出的彩色图像加密技术在信息安全领域具有重要的应用价值.

2 加密过程

本文提出的彩色图像的加密过程分为三基色分量θ调制加密和计算全息编码两个步骤,第一个步骤在图1所示的菲涅耳衍射和θ调制的4f系统中完成.

在图1所示的加密系统中,待加密彩色图像放置于物平面Σ1,紧贴待加密图像的是随机相位板p1(x,y)(random phase mask,RPM1).红色、绿色和蓝色激光器同时照射物平面Σ1加载待加密彩色图像信息,然后经过RPM1的调制和距离为z0的菲涅耳衍射,到达菲涅耳衍射物平面Σ2,衍射过程中依次使用红、绿、蓝三个单色滤波光片分别选出三基色分量.每个分量经过平面Σ2上的透射式振幅型正弦光栅和傅里叶频谱面Σ3的滤波器调制.每加密一个分量,正弦光栅和滤波器同步逆时针旋转 60°,–1级频谱经过傅里叶变换到达像面Σ4被二维图像探测器记录.

本文所选取的待加密彩色图像如图2所示,其图像分布可表达为o(x,y),像素数为 128×128 的三维矩阵,其三基色分量分别为oR(x,y),oG(x,y),oB(x,y),因此可以得到

菲涅耳衍射距离选取z0=0.3m,加密过程可表示为:

其中p1(x,y) 为随机相位,FrTzo[•,633nm] 表示 衍射距离为z0、波长为633 nm 的菲涅耳衍射,Hi(x,y)为光栅条纹与x轴夹角为i°的光栅函数,其中i分别为 0°,60°和 120°,F1{•} 表示与H0(x,y) 对应的滤波过程,其余同理.

图1 彩色图像的菲涅尔衍射和 θ 调制系统Fig.1.Fresnel diffraction and θ modulation system for color image.

图2 待加密彩色图像Fig.2.Color image to be encrypted.

图3 彩色图像加密的计算全息编码流程图Fig.3.Flow chart of color image encryption by computer generated hologram.

第二个步骤为计算全息编码过程,编码流程如图3所示.

为了提高加密系统的安全性,将振幅o′(x,y) 保留而将相位信息r(x,y)去除,以破坏加密系统的线性关系.振幅o′(x,y) 经 过第二个随机相位p2(x,y) 的调制,获得加密复振幅o′(x,y)p2(x,y),其强度分布如图4(a)所示.最后对所获得的加密复振幅o′(x,y)p2(x,y)进行计算全息编码得到二元实值计算全息加密图,在此过程中,采用的计算全息编码单元为 9 pixel× 9 pixel,所以最终得到加密全息图总像素数为 128×128×9×9,如图4(b)所示,完成加密.

传统彩色图像加密方法的加密结果为R,G,B三幅图或单幅彩色图[7−12],而图4(b)显示的加密结果为一幅灰度二元实值计算全息图[21],该图完全隐藏了原始彩色图像的灰度和色彩信息,在传输和存储过程中更具有一般性和迷惑性.

图4 (a)加密结果图; (b)加密计算全息图Fig.4.(a) The encrypted image; (b) the encrypted computer generated hologram.

3 解 密

彩色图像解密为加密的逆过程,在图5所示的系统中完成.该系统包含两个4f系统和一次菲涅耳衍射过程,光源为与加密过程相同的红色、绿色和蓝色激光.在第一个4f系统的入射面放置一个纯相位型空间光调制器(SLM),用以加载加密计算全息图(图4); 在该4f系统的出射面放置第二个SLM,用以加载计算全息相位密钥.其中相位密钥表达式为

图5中由第一个4f系统连接的两个SLM,其作用相当于将计算全息相位密钥紧贴计算全息加密图放置,在该4f系统的出射面上的光场可表示为

图5 彩色图像解密系统Fig.5.Decryption system of color image.

在图5中的第二个4f系统中,频谱面上放置带有三色滤光片的滤波器.由于在加密过程中对三基色分别进行了不同角度的θ调制,因此携带颜色信息的–1级衍射光将分别经过滤波器的不同位置,并通过相应的滤光片从而恢复出色彩值.在此之后,经过距离z0的菲涅耳衍射,此时所恢复的三基色光上述三基色分量叠加得到解密彩色图像,由电荷耦合器接收.在正确的解密密钥、正确的滤波器和正确的菲涅耳衍射距离下,解密结果如图6(a)所示,解密图像为原始彩色图像的共轭图像,经过反转可获得原始彩色图像,如图6(b)所示,可见较好地恢复了彩色图像的灰度和色彩信息,完成解密.

图6 解密彩色图像Fig.6.Decrypted color image.

4 安全性分析和测试

4.1 密钥安全性分析

加密之后的图像通过公共通信信道传输后,存在信息失真的多种可能性.因此,为进一步说明和验证本文所提出方法的可行性和有效性,引入相关系数(correlation coefficien,CC)和图像逼真度 (image fidelity,IF)来评价解密结果的质量,分别定义如下:

本文提出的光学图像加密方法中,除双随机相位密钥以外,菲涅耳衍射距离和滤波器都可作为附加密钥.图7(a)为随机相位密钥p2′错误、其他密钥都正确时的解密结果图,可见解密结果图完全丢失了原始彩色图像的信息,解密结果图与原始图像的相关系数为 C C=0.007,图像逼真度 I F=0.011.图7(b)为菲涅耳衍射距离z0=0.35m 、其他密钥都正确时的解密结果图,解密结果图与原始图像的相关系数为 C C=0.021,图像逼真度 I F=0.043,可见菲涅耳衍射距离作为图像加密密钥也获得了很好的加密效果.

图7 密钥错误时的解密结果 (a)相位密钥 p 2′ 错误;(b)菲涅耳衍射距离错误(z0=0.35m)Fig.7.Decrypted results with wrong keys: (a) Wrong key p2′; (b) wrong key Fresnel diffraction distance with z0=0.35m.

解密时,设置滤波器分别只恢复彩色图像的R,G 和 B 三个分量,如图8(a)—(c)所示,图8(a)—(c)是相应的原始三基色分量,它们对应的相关系数分别为 C C=0.876,C C=0.891,C C=0.903,图像逼真度分别为 I F=0.881,I F=0.900,I F=0.911.从图8可以看出,恢复出来的图像丢失了部分细节信息.造成该现象的原因有两个: 一是由于第一个随机相位密钥的散射作用使得部分高频分量在截取傅里叶频谱时被舍去; 二是加密结果的计算全息量化编码引入了一定误差.经验证,当降低第一个随机相位密钥的动态取值范围并增加计算全息量化编码的像素点数时,解密图像的质量会出现很大提升.因此,原始图像的灰度和色彩信息可以被恢复出来,其恢复质量在一定程度上取得了令人满意的效果.

图8 解密彩色图像三基色分量 (a)红色分量; (b) 绿色分量; (c)蓝色分量; (a'),(b'),(c')原始图像三基色分量Fig.8.Tricolor components of color image: (a) Red component; (b) green component; (c) blue component; (a'),(b'),(c') red,green and blue components of original color image.

图9 不同程度的剪切攻击和剪切攻击后的恢复图像Fig.9.Shear attack of different levels and image restoration after shear attack.

图10 (a)—(d)分别加入密度 0.1,0.3,0.5 和 1.0 的椒盐噪声后的加密图像; (a')—(d')相应的解密结果图Fig.10.(a)－(d) Encrypted images with salt and pepper noises with density of 0.1,0.3,0.5 1.0; (a')－(d') corresponding decryption results.

4.2 抗剪切攻击测试

本文还针对加密图像做了抗剪切和抗噪声模拟实验,图9给出了对加密图像做不同程度剪切后的实验结果和经剪切攻击后恢复得到的解密图像.

表1为针对加密图像做剪切攻击实验的数据,重构的解密图像仍然具有较好的视觉效果.

表1 剪切攻击处理Table 1.Treatment of shear attack.

4.3 抗噪声攻击测试

为了测试该图像加密系统的抗噪声性能,将不同密度的椒盐噪声加入到加密图像中.图10是密度分别为0.1,0.3,0.5和1.0的噪声攻击下的加密图像及其解密结果图.

图10表明,当攻击噪声的密度增大时,加密图像变化明显,但对解密结果影响不大,明文图像的全部信息能够被解密出来,获得了高质量的解密结果图.这说明本方法具有很高的抗噪性能,该加密系统能够抵御噪声的攻击.

将本方法与文献[13,16,18]进行了抗噪性能对比,结果如表2所列.从表2中可以看出,本算法对噪声攻击具有较好的鲁棒性,在噪声密度较小时与文献[13,16,18]相当; 当噪声度较大时,本文方法具有明显优势.因此,该加密系统是安全可行的.

表2 噪声攻击处理Table 2.Treatment of noise attack.

5 结 论

利用θ调制的空间频谱复用和计算全息技术提出了一种新型的彩色图像光学加密方法.加密过程中,首先利用红、绿、蓝三色激光器和三色滤光片将彩色图像分为三个信道,每个信道引入不同方向的透射式振幅型正弦光栅,使得三基色分量在空间频谱面上分离; 然后对三基色分量进行叠加和相位截断非对称处理; 最后进行罗曼型计算全息编码,从而将彩色图像加密成二元实值的灰度CGH.该CGH加密图完全隐藏了原始彩色图像的灰度和色彩信息,在传输和存储过程中更具有一般性和迷惑性.解密过程中,CGH加密图和CGH密钥在4f和菲涅耳衍射系统中完成,计算机模拟结果证明了该方法的有效性和安全性.此外,对该彩色图像加密方法的抗噪声攻击和抗剪切攻击性能进行了研究,并与文献中的方法进行了对比,发现该方法对噪声攻击和剪切攻击具有较好的鲁棒性; 而且当攻击噪声密度较大时,本文方法的优势更为明显.因此,该方法具有高的安全性、抗剪切和抗噪能力,这在信息安全领域具有重要的应用价值.