图像加密Multiple images encryption based on 3D scrambling and hyper-chaotic system

论文信息：

---

abstract

利用混沌理论和elliptic curve Elgamal cryptosystem的概念来生成密文图像和共享密钥。使用多个平面图像生成3D图像，并经历用于生成密码数据的置换和替换阶段。实验结果和分析表明，该算法具有密钥空间大、密码数据统计特性好、抗攻击能力强等特点。

---

1.Introduction

创新点：

1. 通过使用多幅图像形成3D立方体结构来加密多幅图像。
2. 为了防止已知明文攻击的算法，将SHA-512应用于输入图像以生成散列值，该散列值又被用来生成使用椭圆曲线密码算法计算的共享密钥。
3. 为了在一定百分比的遮挡或噪声攻击下提供更好的图像检索，3D立方体沿所有轴进行排列和旋转，以生成置乱的3D立方体。
4. 为了减少计算混沌序列的迭代次数，采用了基变换运算。

2. Proposed multiple image encryption scheme

2.1. Creating 3D cube

多个2D平面图像被分割成128×128，并通过堆叠在彼此之上而结合形成立方体。

2.2. Key exchange

SHA-512从3D图像立方体中生成512位散列值h。利用ECC（Elliptic Curve Cryptography）中点乘的概念，用G生成一个点hG(hG_x,hG_y)，其中G是有限域上椭圆曲线的生成元。Ep: y2≡ x³+ax+b mod p
在通信方之间秘密共享为密文hG’，计算为hg’={kg,(hG+kP_receiver)}其中k是512位的随机整数,n_receiver为接收端的私钥，P_receiver为n_receiver与G的点乘（具体计算可见下方）。
点加：P1（x1，y1），P2（x2，y2）

点加倍：P1（x1，y1）
P3（x3，y3）=P1（x1，y1）+P1（x1，y1）

点减：P1（x1，y1），P2（x2，y2）
P3（x3，y3）=P1（x1，y1）+P2（x2，-y2）
点乘：用点加和点倍增混合计算kP1.

2.3. Generating chaotic sequence using hyper-chaotic system

超混沌系统

a = 35; b =3; c= 12; d = 7 ； e ∈（0.085,0:798）

2.4. Encryption algorithm

加密过程如下：

1. 将n个图像转变成128×128×128的立方体（具体组成所需n的大小如表1）；

2. SHA-512 算法从上述立方体中生成哈希值h,利用ECC中的点乘生成一个点hG(hG_x,hG_y);

3. 将上述坐标转换成512位二进制，并分成4部分：{(b1;;b128)、(b129;;b256)、(b257;;b384)、(b385;;b512)}，再将4部分转变成整数，即{(p;q;r;s)}。

4. 计算超混沌系统的初始参数（w₀，x₀，y₀，z₀）：w₀ =p/2¹²⁸; x₀ =q/2¹²⁸; y₀ =r/2¹²⁸; z₀ =s/2¹²⁸;

5. 由hG_x产生的（w₀，x₀，y₀，z₀）作为超级混沌系统初始输入，另e取初试输入的平均值。每隔e个距离取值，形成128的混沌序列（x_i，y_i，z_i）。对（x_i，y_i，z_i）排序得到（sx_i，sy_i，sz_i）,对比两者位置变化得到置换表（P_x,P_y,P_z）。

6. 置乱：取Pz的值，首先沿z平面对3D图像立方体进行切片置换。将所获得的3D图像立方体沿x平面旋转90°，并使用Px的值沿z平面进行切片置换。类似地，将所获得的3D图像立方体沿y平面旋转90°，并通过沿z平面取Py的值来执行切片排列，从而得到置乱图像S_IMG。

7. 扩散：由hG_y产生的（w₀，x₀，y₀，z₀）作为超级混沌系统初始输入，得到序列（w_i，x_i，y_i，z_i），进而得到取值范围（0_{255）的序列。此处（w}i_，xi_，yi_，zi~）的取值间隔为e/100,
   w_byte=(|w_i|×10¹⁵)₂₅₆[[2;;7]];
   x_byte=(|x_i|×10¹⁵)₂₅₆[[2;;7]];
   y_byte=(|y_i|×10¹⁵)₂₅₆[[2;;7]];
   z_byte=(|z_i|×10¹⁵)₂₅₆[[2;;7]];
   (以256为基数执行基数转换操作，基数转换后，取位置2到7的值，因为第一位的值多在0~50，导致数值大小不够平均，故去掉。）
   举例 ：514181197950012={1;211;165;35;167;108;60}₂₅₆
   
   再将上述w_byte，x_byte，y_byte，z~byte进行组合得到
   P₁=S(w_byte,x_byte)
   P₂=S(w_byte,y_byte)
   P₃=S(w_byte,z_byte)
   P₄=S(x_byte,y_byte)
   P₅=S(x_byte,z_byte)
   P₆=S(y_byte,z_byte)
   P₇=S(P₁,P₆)
   S混洗操作，如果w_byte中的项目数小于x_byte，则在w_n处终止。如果w_byte中的项目数大于x_byte，则终止发生在w_n，从x_n取值。如果项目的数量在w_byte和x_byte中都相等，则在x_n处终止。

8. 应用所有的P_i生成混沌立方体C，大小128×128×128，有三个通道。

9. 对S_IMG和C的每个像素进行XORed异或操作得到加密的C_IMG。

2.5. Decryption algorithm

解密过程如下：

1. 得到加密的C_IMG。
2. 从hG’中计算hG。
3. hG_x和hG_y通过超混沌系统生成两个混沌序列。
4. 用hG_y根据7.8步中的算法，计算出C。
5. C_IMG与C异或的到S_IMG。
6. 利用hG_x，得到步骤5中的（P_x,P_y,P_z）。
7. 由（P_x,P_y,P_z）得到逆置换表（RP_x,RP_y,RP_z）。
8. 用RP_y沿着z平面对加扰的3D图像立方体（C）进行解扰。生成的3D图像立方体沿y轴旋转90度。再使用RP_x沿z平面解扰，并沿x轴旋转90。最后，利用RP_z对z平面进行解扰，得到纯3D图像立方体.
9. 通过从3D图像立方体中取出2D切片并重新组装，可以从解密的3D图像立方体中获得n个平面图像。

3. Simulation and analysis of the proposed scheme

3.1. Correlation analysis
3.2. Histogram analysis
3.3. Key space and key sensitivity analysis
3.4. Complexity of the proposed algorithm
3.5. Stability of the algorithm
3.6. Differential attack
3.7. Entropy
3.8. PSNR and SSIM
3.9. Occlusion and salt pepper noise attack
3.10. Execution Speed
3.11. Randomness test

4. Conclusion

---

关于加密中的扩散操作笔者还未弄清楚具体实现！！！！