Security Analysis of a Color Image Encryption Algorithm Using a Fractional-Order Chaos

Abstract

Fractional-order chaos has complex dynamic behavior characteristics, so its application in secure communication has attracted much attention. Compared with the design of fractional-order chaos-based cipher, there are fewer researches on security analysis. This paper conducts a comprehensive security analysis of a color image encryption algorithm using a fractional-order hyperchaotic system (CIEA-FOHS). Experimental simulation based on excellent numerical statistical results supported that CIEA-FOHS is cryptographically secure. Yet, from the perspective of cryptanalysis, this paper found that CIEA-FOHS can be broken by a chosen-plaintext attack method owing to its some inherent security defects. Firstly, the diffusion part can be eliminated by choosing some special images with all the same pixel values. Secondly, the permutation-only part can be deciphered by some chosen plain images and the corresponding cipher images. Finally, using the equivalent diffusion and permutation keys obtained in the previous two steps, the original plain image can be recovered from a target cipher image. Theoretical analysis and experimental simulations show that the attack method is both effective and efficient. To enhance the security, some suggestions for improvement are given. The reported results would help the designers of chaotic cryptography pay more attention to the gap of complex chaotic system and secure cryptosystem.

1. Introduction

Nowadays, with the rapid development of optical fiber broadband access network, 5G and other communication technologies, the security of multimedia data, especially digital images, is of particular interest in communication networks [1]. As everyone knows, encryption is an effective means of achieving security enhancements [2]. However, traditional text encryption algorithms such as AES, DES, and IDEA are not suitable for digital images because they featured with strong correlation between adjacent pixels. To deal with the problem, various methodologies are introduced to design different image ciphers. Among them, chaos-based image encryption is the most popular one, because chaos has characteristics of sensitivity to initial values, dense periodic points, and long-term unpredictability of orbits [3,4,5]. In the past two decades, chaotic image encryption technology has been widely discussed and has become a research hotspot [6]. To improve the security performance of chaotic image encryption technology, various chaotic systems with resistance to dynamic degradation are studied, including quantum chaotic map [7], fractional-order chaos [8], non-degenerated hyperchaos [9], economic chaotic map [10], and cascaded chaotic systems [11], etc. However, chaotic cryptography still lacks authoritative metrics, especially in terms of security. Accordingly, many reported chaotic encryption algorithms have been broken [12,13,14,15]. As shown in

Table 1. Some chaos-based ciphers broken by various attack methods.

Ciphers	Broken by	Attack Methods
Fridrich et al. [34] in 1998	Xie et al. [16] in 2017	Chosen-ciphertext attack
Zhao et al. [35] in 2015	Norouzi et al. [36] in 2017	Chosen-plaintext attack
Ye [37] in 2010	Li et al. [18] in 2017	Cipher-only attack
Zhou [38] in 2015	Chen et al. [17] in 2016	Differential cryptanalysis
Song et al. [15] in 2015	Wen et al. [13] in 2019	Chosen-plaintext/cipertext attacks
Shafique et al. [14] in 2018	Wen et al. [12] in 2019	Chosen-plaintext attack

, some previous chaos-based ciphers are vulnerable upon various attack methods, including chosen-ciphertext attack [16], chosen-/known-plaintext attack [12], differential cryptanalysis [17], even cipher-only attack [18]. Therefore, research on security is extremely important and has received much attention [19,20,21,22,23,24,25,26,27,28,29,30,31,32,33].

As described in Ref. [39], fractional-order chaotic systems have higher complexity and more optional key parameters and can be used as a competitive encryption scheme. Correspondingly, image encryption algorithms based on fractional-order chaotic systems have attracted the attention of researchers in recent years [35,40,41,42]. In 2013, Wang et al. [40] introduced a fractional-order chaos into image encryption for the first time, and gave some experiments to verify its performance. Since then, many image encryption schemes based on fractional-order chaotic systems have been proposed [35,41,42]. For example, in 2017, Zhang et al. [41] proposed a color image encryption scheme combing with fractional-order hyperchaotic system and DNA encoding. Yet, cryptanalysts have reported that some fractional-order chaotic image encryption algorithms have some fatal security issues. Exactly, Norouzi et al. [36] pointed out that the image cipher that using an improper fractional-order chaotic system was insecure, which was proposed in [35]. As far as we know, there are still few research studies concerning cryptanalysis on the ciphers based on fractional-order chaotic systems. Moreover, considering that each cryptosystem has its intrinsic characteristics, it is necessary and urgent to perform cryptanalysis on these existing ciphers.

In 2015, a color image encryption algorithm based on a fractional-order hyperchaotic system was proposed [42]. In color image encryption algorithm using a fractional-order hyperchaotic system (CIEA-FOHS), using the pseudo-random sequences generated by the fractional-order hyperchaotic system, RGB-inter permutation, RGB-intra permutation and pixel diffusion are successively performed to get cipher images from plain images. Meanwhile, the relevant pixel correlation, histogram and other experimental analysis are given to verify its security performance. However, from the perspective of cryptanalysis, we found some security defects as follows:

• The existence of an equivalent key. CIEA-FOHS encrypts the image using a pseudo-random sequence generated by fractional-order chaos. However, these sequences are not related to plaintext. Thus, these sequences can be considered as equivalent keys.
• Two-stage permutations can be equivalently simplified to only once. The reason is that the two permutations only change the position of the pixel without changing the value of the pixel.
• The paradigm of the diffusion part is insecure. According to the conclusion of Ref. [43], a class of diffusion encryption using module addition and XOR operations can be cracked with only two special plain images and their corresponding cipher images. Unfortunately, CIEA-FOHS is also the case.

Based on the three points, CIEA-FOHS cannot resist against a chosen-plaintext attack method with the divide-and-conquer strategy. More specifically, under the scenario of chosen-plaintext attack, firstly an equivalent diffusion key is obtained, and then an equivalent permutation key is achieved, and finally the original images can be restored from the encrypted images with the equivalent keys.

2. The Encryption Algorithm under Study

In this section, the fractional-order hyperchaotic system used in Reference [42] is presented, and then the specific steps of CIEA-FOHS are introduced.

2.1. Fractional-Order Hyperchaotic System

The fractional-order hyperchaotic system used in CIEA-FOHS is derived from Ref. [39], given as

$$\left\{\begin{array}{c}{D}_t^{\alpha }x\left(t\right)=-z-w\hfill \\ D_t^{\alpha }y\left(t\right)=2y+z\hfill \\ D_t^{\alpha }z\left(t\right)=14x-14y\hfill \\ D_t^{\alpha }w\left(t\right)=100(x-g\left(w\right))\hfill \end{array}\right.$$

(1)

where $x,y,z,w$ are the four state variables, $g\left(w\right)=w-(\left|w-0.4\right|-\left|w-0.8\right|-\left|w+0.4\right|-\left|w+0.8\right|)$, $D_t^{\alpha }$ is the fractional derivative under the definition of Caputo and $\alpha$ is the derivative order. The attractor of the fractional-order hyperchaotic system is shown in Figure 1.

2.2. Description of CIEA-FOHS

As shown in Figure 2, CIEA-FOHS consists of three main parts: inter-permutation, intra-permutation and pixel diffusion. It is noted that, a two-dimensional image is transformed into an one-dimensional sequence in raster scan order. Specifically, a color plain image $\mathit{I}$ of size $H\times W\times 3$ is converted into three sequences of length $H\times W$ expressed as: $\mathit{IR}$, $\mathit{IG}$, and $\mathit{IB}$, which correspond to the three RGB channels of the image. The main contents are briefly introduced as follows:

• The Secret Key:

The secret keys of CIEA-FOHS include $(t_f,\alpha ,h,x_0,y_0,z_0,w_0)$, where $t_f$ is the fractional derivative defined by Caputo definition, $\alpha$ is the dimension, h is the step size for discretization, and $(x_0,y_0,z_0,w_0)$ are the four initial values of the fractional-order hyperchaotic system defined in Equation (1), respectively. In CIEA-FOHS, these keys are used to generate some chaos-based pseudo-random sequences for encryption [42].

• Initialization:

In Equation (1), by selecting the secret key as the initial values and parameters and iterating L times, one gets four chaos-based pseudo-random sequences ${\left\{x_i\right\}}_{i=1}^L$, ${\left\{y_i\right\}}_{i=1}^L$, ${\left\{z_i\right\}}_{i=1}^L$ and ${\left\{w_i\right\}}_{i=1}^L$, where $L=H\times W$ represents the number of pixels in a single image channel.

• Stage 1. RGB-inter permutation:

The RGB-inter permutation refers to the process of pixel replacement between channels. This stage is implemented by two control vectors ${\left\{sel{E}_i\right\}}_{i=1}^L$ and ${\left\{selLe{n}_i\right\}}_{i=1}^L$, which are given as

$$\left\{\begin{array}{c}sel{E}_i=\left(\right|x_i|\times {10}^{14})mod6\hfill \\ selLe{n}_i=\left(\right|z_i|\times {10}^{14})mod3\hfill \end{array}\right.$$

(2)

where $i=1\sim L$. More specifically, ${\left\{sel{E}_i\right\}}_{i=1}^L$ is used to switch channels, as shown in

Table 2. The stutas of RGB-inter permutation under six rules.

Rule $\mathit{selE}\left(\mathit{i}\right)$	0	1	2	3	4	5
	$R\to R$	$R\to R$	$R\to G$	$R\to B$	$R\to G$	$R\to B$
Permutation status	$G\to G$	$G\to B$	$G\to R$	$G\to R$	$G\to B$	$G\to G$
	$B\to B$	$B\to G$	$B\to B$	$B\to G$	$B\to R$	$B\to R$

, and ${\left\{selLe{n}_i\right\}}_{i=1}^L$ is to control the position and length of the permutation pixel, given as

$$\left\{\begin{array}{c}length=(sum\left(ER(pos:pos+length-1)\right)mod64),ifselLe{n}_i=0\hfill \\ length=(sum\left(EG(pos:pos+length-1)\right)mod64),ifselLe{n}_i=1\hfill \\ length=(sum\left(EB(pos:pos+length-1)\right)mod64),ifselLe{n}_i=2\hfill \end{array}\right.$$

(3)

where $pos$ is the starting position, $length$ is the length of the permautation pixels, and $sum$ is the cumulative function.

• Stage 2. RGB-intra permutation:

Sort ${\left\{y_i\right\}}_{i=1}^L$, ${\left\{z_i\right\}}_{i=1}^L$, and ${\left\{w_i\right\}}_{i=1}^L$ to get three index sequences ${\left\{I{Y}_i\right\}}_{i=1}^L$, ${\left\{I{Z}_i\right\}}_{i=1}^L$, and ${\left\{I{W}_i\right\}}_{i=1}^L$ respectively, and their values range $[1,L]$. Use ${\left\{I{Y}_i\right\}}_{i=1}^L$, ${\left\{I{Z}_i\right\}}_{i=1}^L$, and ${\left\{I{W}_i\right\}}_{i=1}^L$ to permute $\mathit{ER}$, $\mathit{EG}$ and $\mathit{EB}$ respectively, given as $E{R}_i=ER\left(I{Y}_i\right)$, $E{G}_i=EG\left(I{Z}_i\right)$ and $E{B}_i=EB\left(I{W}_i\right)$.

• Stage 3. Pixel diffusion:

Perform pixel diffusion on $\mathit{ER}$, $\mathit{EG}$ and $\mathit{EB}$, and then get three channels of the cipher image $\mathit{C}$. Exactly, the three channels $\mathit{CR}$, $\mathit{CG}$ and $\mathit{CB}$ are defined as

$$\left\{\begin{array}{c}C{R}_i=S{X}_i\oplus ((E{R}_i+S{X}_i)mod256)\oplus C{R}_{i-1}\hfill \\ C{G}_i=S{Y}_i\oplus ((E{G}_i+S{Y}_i)mod256)\oplus C{G}_{i-1}\hfill \\ C{B}_i=S{Z}_i\oplus ((E{B}_i+S{Z}_i)mod256)\oplus C{B}_{i-1}\hfill \end{array}\right.$$

(4)

where $i=1\sim L$, ⊕ is bitwise XOR operation, $mod$ represents modulo operation, and $C{R}_0=S{X}_L$, $C{G}_0=S{Y}_L$, and $C{B}_0=S{Z}_L$. Here, three diffusion sequences $\mathit{SX}$, $\mathit{SY}$ and $\mathit{SZ}$ are generated by $S{X}_i=round\left(x_i\right)\times {10}^{14}$, $S{Y}_i=round\left(y_i\right)\times {10}^{14}$ and $S{Z}_i=round\left(z_i\right)\times {10}^{14}$ respectively, where $round$ is a rounding operation on real numbers.

Decryption is the inverse of encryption and is not described in detail here.

3. Security Analysis of CIEA-FOHS

3.1. Preliminary Analysis of CIEA-FOHS

Referring to the basic assumptions of cryptanalysis, everything about the cryptosystem is public and only the secret key is unknown for attackers [13]. Chosen-plaintext attack is a common and powerful method of cryptanalysis. It assumes that attackers can arbitrarily choose the plaintext that is conducive to deciphering and obtain the corresponding ciphertext [12]. Under the scenario of chosen-plaintext attack, attackers can construct special plain images, such as all black and all white, and obtain the corresponding cipher images to analyze the target cipher.

From the perspective of cryptanalysis, two-stage permutations of CIEA-FOHS can be treated as a global pixel permutation because they only change the pixels’ position without their values. The difference is that the number of pixels performing the permutation is $3HW$ instead of $HW$. Then, the algorithm structure of CIEA-FOHS is actually a classic single-round permutation-diffusion. Moreover, the generation process of all chaos-based pseudo-random sequences is independent of the plain image, which means that these sequences can be regarded as an equivalent key. The reason is that, in the case of a given secret key, these sequences are fixed for encrypting different plain images with the same size. Then, CIEA-FOHS can be equivalently simplified as Figure 3, where $\mathit{PM}$ is an equivalent permutation key and three diffusion sequences $\mathit{SX}$, $\mathit{SY}$ and $\mathit{SZ}$ serve as an equivalent diffusion key.

Based on the above, under the scenario of chosen-plaintext attack and the strategy of divide and conquer, one can get the equivalent keys and then recover the original plain images. Specifically, firstly choose some plain images with same pixel values to cancel the permutation and get the corresponding plain images to obtain the diffusion key; then achieve the permutation key by the method of Reference [12]; finally, recover the images by the equivalent keys.

3.2. Analysis on the Diffusion Part

In this section, based on chosen-plaintext attack, it is assumed that the plaintext image with the same pixel value is selected, and the corresponding ciphertext image is obtained.

• Step 1. Choose the all-zero plain image $\mathit{I}$${}^{\left(0\right)}$ and get the corresponding cipher image ${\mathit{C}}^{\left(0\right)}$ to determine $S{X}_L,S{Y}_L,S{Z}_L$.

The reason for choosing the all-zero image is that the permutation is invalid at this time, and the diffusion can be eliminated to the greatest extent. Then, Equation (4) becomes

$$\left\{\begin{array}{c}C{R}_i^{\left(0\right)}=C{R}_{i-1}^{\left(0\right)}\hfill \\ C{G}_i^{\left(0\right)}=C{G}_{i-1}^{\left(0\right)}\hfill \\ C{B}_i^{\left(0\right)}=C{B}_{i-1}^{\left(0\right)}\hfill \end{array}\right.$$

(5)

when $i=1$, one has $C{R}_1^{\left(0\right)}=C{R}_0$. Since $C{R}_0=S{X}_L$, thus $S{X}_L=C{R}_i^{\left(0\right)}$. Similarly, one further gets $S{Y}_L=C{G}_i^{\left(0\right)}$ and $S{Z}_L=C{B}_i^{\left(0\right)}$.

• Step 2. Choose two special plain images and get the corresponding cipher images to determine $S{X}_i,S{Y}_i,S{Z}_i$ for $i=1\sim L-1$.

Referring to [43,44], the two chosen plaintexts are pure-color images with pixel values of 85 and 170, represented as $\mathit{I}$${}^{\left(85\right)}$ and $\mathit{I}$${}^{\left(170\right)}$, respectively. Because for the combined operation of module addition and bitwise XOR, choosing these two plain images can minimize the number of solutions for $\mathit{SX}$, $\mathit{SY}$, $\mathit{SZ}$. Under the plain image $\mathit{I}$${}^{\left(85\right)}$ and its corresponding cipher image $\mathit{C}$${}^{\left(85\right)}$, one gets

$$\left\{\begin{array}{c}C{R}_i^{\left(85\right)}=S{X}_i\oplus ((85+S{X}_i)mod256)\oplus C{R}_{i-1}^{\left(85\right)}\hfill \\ C{G}_i^{\left(85\right)}=S{Y}_i\oplus ((85+S{Y}_i)mod256)\oplus C{G}_{i-1}^{\left(85\right)}\hfill \\ C{B}_i^{\left(85\right)}=S{Z}_i\oplus ((85+S{Z}_i)mod256)\oplus C{B}_{i-1}^{\left(85\right)}\hfill \end{array}\right.$$

(6)

Similarly, given the plain image $\mathit{I}$${}^{\left(170\right)}$ and its corresponding cipher image ${\mathit{C}}^{\left(170\right)}$, one has

$$\left\{\begin{array}{c}C{R}_i^{\left(170\right)}=S{X}_i\oplus ((170+S{X}_i)mod256)\oplus C{R}_{i-1}^{\left(170\right)}\hfill \\ C{G}_i^{\left(170\right)}=S{Y}_i\oplus ((170+S{Y}_i)mod256)\oplus C{G}_{i-1}^{\left(170\right)}\hfill \\ C{B}_i^{\left(170\right)}=S{Z}_i\oplus ((170+S{Z}_i)mod256)\oplus C{B}_{i-1}^{\left(170\right)}\hfill \end{array}\right.$$

(7)

By performing bitwise on Equations (6) and (7), one further gets

$$\left\{\begin{array}{c}\left(85\dot{+}S{X}_i\right)\oplus \left(170\dot{+}S{X}_i\right)=C{R}_i^{\left(85\right)}\oplus C{R}_{i-1}^{\left(85\right)}\oplus C{R}_i^{\left(170\right)}\oplus C{R}_{i-1}^{\left(170\right)}\hfill \\ \left(85\dot{+}S{Y}_i\right)\oplus \left(170\dot{+}S{Y}_i\right)=C{G}_i^{\left(85\right)}\oplus C{G}_{i-1}^{\left(85\right)}\oplus C{G}_i^{\left(170\right)}\oplus C{G}_{i-1}^{\left(170\right)}\hfill \\ \left(85\dot{+}S{Z}_i\right)\oplus \left(170\dot{+}S{Z}_i\right)=C{B}_i^{\left(85\right)}\oplus C{B}_{i-1}^{\left(85\right)}\oplus C{B}_i^{\left(170\right)}\oplus C{B}_{i-1}^{\left(170\right)}\hfill \end{array}\right.$$

(8)

where $\dot{+}$ is defined as $a\dot{+}b\stackrel{\Delta }{=}mod(a+b,256)$. It is worth pointing out that the reason why 85 and 170 are chosen as the attack images is that their binary are 01010101 and 10101010 respectively. At this time, the number of possible solutions of $S{X}_i,S{Y}_i,S{Z}_i$ is the smallest, which is two. More precisely, the difference between the two solutions is 128. Then, based on Equation (8), we propose Alogrithm 1 to determine $S{X}_i,S{Y}_i,S{Z}_i$, where $i=1\sim L-1$.

Algorithm 1: Determining $S{X}_i,S{Y}_i,S{Z}_i$ for $i=1\sim L-1$

• Step 3. Eliminate the diffusion part by $\mathit{SX}$, $\mathit{SY}$, $\mathit{SZ}$.

Corresponding to Equation (4), the decryption process of diffusion is given as

$$\left\{\begin{array}{c}E{R}_i=(S{X}_i\oplus C{R}_i\oplus C{R}_{i-1}-S{X}_i)mod256\hfill \\ E{G}_i=(S{Y}_i\oplus C{G}_i\oplus C{G}_{i-1}-S{Y}_i)mod256\hfill \\ E{B}_i=(S{Z}_i\oplus C{B}_i\oplus C{B}_{i-1}-S{Z}_i)mod256\hfill \end{array}\right.$$

(9)

Thus, $\mathit{ER}$, $\mathit{EG}$, $\mathit{EB}$ can be restored from $\mathit{CR}$, $\mathit{CG}$, $\mathit{CB}$ with $\mathit{SX}$, $\mathit{SY}$, $\mathit{SZ}$, respectively.

3.3. Analysis on the Permutation Part

Once the diffusion part is broken, CIEA-FOHS degenerates into a permutation-only cipher. Based on existing research, it cannot resist a chosen-plaintext attack. The basic idea of attacking permutation-only is to construct a special plain image with unequal element values, and get the corresponding permuted image. Taking $2\times 2\times 3$ as an example, the process of solving $\mathit{PM}$ is described below. First, a chosen plain image and the corresponding permuted image are given as

$$\mathit{IR}=\left[\begin{array}{cc}0& 1\\ 2& 3\end{array}\right];\mathit{IG}=\left[\begin{array}{cc}4& 5\\ 6& 7\end{array}\right];\mathit{IB}=\left[\begin{array}{cc}8& 9\\ 10& 11\end{array}\right]$$

$$\mathit{ER}=\left[\begin{array}{cc}5& 8\\ 3& 11\end{array}\right];\mathit{EG}=\left[\begin{array}{cc}1& 10\\ 2& 9\end{array}\right];\mathit{EB}=\left[\begin{array}{cc}6& 4\\ 0& 7\end{array}\right]$$

For ease of explanation, a matrix of size $H\times 3W$ is obtained by connecting three channels of size $H\times W$ in a row connection manner. Then, the permutation process can be described by

$$\left[\begin{array}{cccccc}0& 1& 4& 5& 8& 9\\ 2& 3& 6& 7& 10& 11\end{array}\right]\stackrel{\mathit{PM}}{\to }\left[\begin{array}{cccccc}5& 8& 1& 10& 6& 4\\ 3& 11& 2& 9& 0& 7\end{array}\right]$$

where $\mathit{PM}$ is the permutation matrix of size $H\times 3W$. Finally, $\mathit{PM}$ is determined as

$$\mathit{PM}=\left[\begin{array}{cccccc}(2,5)& (1,3)& (1,6)& (1,1)& (1,2)& (2,4)\\ (2,3)& (2,1)& (1,5)& (2,6)& (1,4)& (2,2)\end{array}\right]$$

(10)

Obviously, one can recover $(\mathit{IR},\mathit{IG},\mathit{IB})$ from $(\mathit{ER},\mathit{EG},\mathit{EB})$ with $\mathit{PM}$. However, the situation may be more complicated for large size images. For an 8-bit image, the pixel value range is $[0,255]$. Thus, when $3HW>256$, $\mathit{PM}$ cannot be determined by only one chosen plain image and its corresponding cipher image. Fortunately, this problem has been solved in our latest research [12,13]. The basic idea is to combine multiple chosen plain images in a weighted manner to form a matrix with different elements, and the number of chosen plain images required for attacking permutation is $⌈{log}_{256}\left(3HW\right)⌉$, where $⌈.⌉$ is the rounding up operation.

Based on the above, the steps for attacking permutation are briefly summarized as follows:

• Step 1. Choose some special plain images and get their corresponding cipher images to determine the permutation matrix $\mathit{PM}$;
• Step 2. Use the permutation matrix $\mathit{PM}$ to recover the original images from the permuted images.

3.4. The Proposed Chosen-Plaintext Attack Method

Following the above-mentioned discussion, CIEA-FOHS cannot resist the attack method proposed in this paper. The flowchart of the attack method is shown in Figure 4, and the specific steps based on chosen-plaintext attack are given as: firstly, get an equivalent diffusion key $(\mathit{SX},\mathit{SY},\mathit{SZ})$ by the method in Section 3.2; secondly, achieve the permutation matrix $\mathit{PM}$ by the method in Section 3.3; finally, recover the original images with the equivalent keys.

Moreover, the complexity required for the attack method is discussed here. In terms of data complexity, for color images of size $H\times W\times 3$, the number of chosen plain images required to decipher diffusion and permutation is 3 and $⌈{log}_{256}\left(3HW\right)⌉$, respectively. Hence, the total data complexity required is $\mathit{O}(3+⌈{log}_{256}\left(3HW\right)⌉)$.

4. Experimental Verifications and Discussions

To verify our security analysis, the algorithm steps of CIEA-FOHS strictly follow Ref. [42]. Although Due to the complexity of fractional-order chaos, some parameters may not be completely consistent, but this does not affect the effectiveness of security analysis. We conduct simulation verification on the proposed image cryptosystem based on a PC (personal computer) with MATLAB r2018b. The running PC is installed with Windows 10 64-bit OS (operating system), Intel(R) Core(TM) i5-8265U CPU @ 1.60 GHz and 8 GB memory. We select some typical images listed in

Table 3. The time required for breaking CIEA-FOHS by our proposed attack method (unit: second).

Images	Sizes	Encrytion Time	Attacking Diffusion			Attacking Permutation		Totol Attacking Time
			Step 1	Step 2	Step 3	Step 1	Step 2
Figure 5a	2 × 2 × 3	0.0280	0.1559	0.1811	1.0297	0.0244	2.7151	4.1502
Figure 13b	100 × 100 × 3	0.1539	0.0920	19.6092	1.1407	0.2764	2.7102	24.0427
Figure 13d	300 × 200 × 3	0.3280	0.5092	101.7737	0.7872	0.9055	2.4353	106.8545
Figure 12e	256 × 256 × 3	0.6391	0.6913	120.4768	1.6147	1.9642	3.7725	129.4039
Figure 13f	512 × 512 × 3	3.5386	2.8134	988.3704	1.9930	4.2884	5.0459	1004.4617

for experiments. Among them, the image “Lenna” of size $256\times 256\times 3$ given in Ref. [42] is also included. In Equation (1), we set the experimental secret key parameters for h = 0.001, $\alpha$ = 104, $t_f$ = 100, $x_0$ = 1.002, $y_0$ = 0.949, $z_0$ = 0.997 and $w_0$ = 1.103.

• Case 1. Breaking CIEA-FOHS with an image of size $2\times 2\times 3$:

In order to better illustrate the attack process, we first adopt an extremely simple image with a size of 2 × 2 × 3. A pair of the given target plain and cipher images $\mathit{I}$ and $\mathit{C}$ is shown in Figure 5a,c respectively, and their histograms are shown in Figure 5b,d respectively. Accordingly, the numerical matrices of $\mathit{I}$ and $\mathit{C}$ are:

$$\mathit{IR}=\left[\begin{array}{cc}11& 22\\ 33& 44\end{array}\right];\mathit{IG}=\left[\begin{array}{cc}55& 66\\ 77& 88\end{array}\right];\mathit{IB}=\left[\begin{array}{cc}99& 100\\ 111& 122\end{array}\right]$$

$$\mathit{CR}=\left[\begin{array}{cc}70& 165\\ 103& 145\end{array}\right];\mathit{CG}=\left[\begin{array}{cc}231& 154\\ 118& 28\end{array}\right];\mathit{CB}=\left[\begin{array}{cc}181& 24\\ 171& 165\end{array}\right]$$

Firstly, following Step 1 in Section 3.2, choose the all-zero plain image ${\mathit{I}}^{\left(0\right)}$ shown in Figure 6a and temporarily use the encryption machine of CIEA-FOHS, and then get the corresponding cipher image ${\mathit{C}}^{\left(0\right)}$, as shown in Figure 6c. The all-zero plain image ${\mathit{I}}^{\left(0\right)}$ and the corresponding cipher image ${\mathit{C}}^{\left(0\right)}$ and their histograms are shown in Figure 6b,d, respectively. Similarly, the numerical matrices of ${\mathit{I}}^{\left(0\right)}$ and ${\mathit{C}}^{\left(0\right)}$ are:

$${\mathit{IR}}^{\left(0\right)}=\left[\begin{array}{cc}0& 0\\ 0& 0\end{array}\right];{\mathit{IG}}^{\left(0\right)}=\left[\begin{array}{cc}0& 0\\ 0& 0\end{array}\right];{\mathit{IB}}^{\left(0\right)}=\left[\begin{array}{cc}0& 0\\ 0& 0\end{array}\right]$$

$${\mathit{CR}}^{\left(0\right)}=\left[\begin{array}{cc}77& 77\\ 77& 77\end{array}\right];{\mathit{CG}}^{\left(0\right)}=\left[\begin{array}{cc}174& 174\\ 174& 174\end{array}\right];{\mathit{CB}}^{\left(0\right)}=\left[\begin{array}{cc}109& 109\\ 109& 109\end{array}\right]$$

Then, one has $S{X}_L=77$, $S{Y}_L=174$ and $S{Z}_L=109$ because $S{X}_L=C{R}_0$, $S{Y}_L=C{G}_0$ and $S{Z}_L=C{B}_0$, where $L=2\times 2=4$.

Secondly, based on Step 2 in Section 3.2, choose the two plain images ${\mathit{I}}^{\left(85\right)}$ and ${\mathit{I}}^{\left(170\right)}$, and get the corresponding cipher images, ${\mathit{C}}^{\left(85\right)}$ and ${\mathit{C}}^{\left(170\right)}$, which are shown in Figure 7a–d, respectively. The values of their RGB three channels are:

$${\mathit{IR}}^{\left(85\right)}=\left[\begin{array}{cc}85& 85\\ 85& 85\end{array}\right];{\mathit{IG}}^{\left(85\right)}=\left[\begin{array}{cc}85& 85\\ 85& 85\end{array}\right];{\mathit{IB}}^{\left(85\right)}=\left[\begin{array}{cc}85& 85\\ 85& 85\end{array}\right]$$

$${\mathit{CR}}^{\left(85\right)}=\left[\begin{array}{cc}176& 186\\ 77& 85\end{array}\right];{\mathit{CG}}^{\left(85\right)}=\left[\begin{array}{cc}5& 181\\ 110& 24\end{array}\right];{\mathit{CB}}^{\left(85\right)}=\left[\begin{array}{cc}184& 94\\ 229& 241\end{array}\right]$$

$${\mathit{IR}}^{\left(170\right)}=\left[\begin{array}{cc}170& 170\\ 170& 170\end{array}\right];{\mathit{IG}}^{\left(170\right)}=\left[\begin{array}{cc}170& 170\\ 170& 170\end{array}\right];{\mathit{IB}}^{\left(170\right)}=\left[\begin{array}{cc}170& 170\\ 170& 170\end{array}\right]$$

$${\mathit{CR}}^{\left(170\right)}=\left[\begin{array}{cc}231& 235\\ 177& 81\end{array}\right];{\mathit{CG}}^{\left(170\right)}=\left[\begin{array}{cc}120& 24\\ 174& 238\end{array}\right];{\mathit{CB}}^{\left(170\right)}=\left[\begin{array}{cc}199& 123\\ 45& 1\end{array}\right]$$

Then, combining Algorithm 1, we determine $\mathit{SX}$ $\mathit{SY}$ $\mathit{SZ}$ as

$$\mathit{SX}=\left[\begin{array}{cccc}84& 86& 89& 77\end{array}\right];\mathit{SY}=\left[\begin{array}{cccc}63& 31& 71& 46\end{array}\right];\mathit{SZ}=\left[\begin{array}{cccc}64& 36& 119& 109\end{array}\right]$$

or

$$\mathit{SX}=\left[\begin{array}{cccc}212& 214& 217& 205\end{array}\right];\mathit{SY}=\left[\begin{array}{cccc}191& 159& 199& 174\end{array}\right];\mathit{SZ}=\left[\begin{array}{cccc}192& 164& 247& 237\end{array}\right]$$

Thirdly, by Step 3 in Section 3.2, the corresponding permuted image shown in Figure 8c can be restored from the targeted cipher image Figure 8a with $\mathit{SX}$ $\mathit{SY}$ $\mathit{SZ}$. Fourthly, following Step 1 in Section 3.3, construct some special attack images to obtain the permutation matrix $\mathit{PM}$. For images of size $2\times 2\times 3$, the process of solving $\mathit{PM}$ is exactly the same as Section 3.3. Then, we determine the $\mathit{PM}$ as Equation (10). Fifth, by Step 2 in Section 3.3, recover $(\mathit{IR},\mathit{IG},\mathit{IB})$ from $(\mathit{ER},\mathit{EG},\mathit{EB})$ with $\mathit{PM}$. Thus, the original plain image shown in Figure 8e can be recovered.

• Case 2. Breaking CIEA-FOHS with “Lenna” of size $256\times 256\times 3$:

Firstly, following Step 1 in Section 3.2, choose the all-zero plain image ${\mathit{I}}^{\left(0\right)}$ shown in Figure 9a and temporarily use the encryption machine of CIEA-FOHS, and then get the corresponding cipher image ${\mathit{C}}^{\left(0\right)}$, as shown in Figure 9b, and the corresponding three channel images and their histograms of ${\mathit{C}}^{\left(0\right)}$ are shown in Figure 9c,d, respectively. Exactly, one has $S{X}_L=238$, $S{Y}_L=168$ and $S{Z}_L=91$ owing to $S{X}_L=C{R}_0$, $S{Y}_L=C{G}_0$ and $S{Z}_L=C{B}_0$.

Secondly, based on Step 2 in Section 3.2, choose the two plain images, ${\mathit{I}}^{\left(85\right)}$ and ${\mathit{I}}^{\left(170\right)}$, and get the corresponding cipher images, ${\mathit{C}}^{\left(85\right)}$ and ${\mathit{C}}^{\left(170\right)}$, which are shown in Figure 10a–d, respectively.

Furthermore, one determines $S{X}_i,S{Y}_i,S{Z}_i$ for $i=1\sim L-1$ by Algorithm 1.

Thirdly, by the method in Section 3.3, choose the three plain images (shown in Figure 11a–f) and get the corresponding cipher images (shown in Figure 11g–l), and then use Algorithm 1 again to obtain their corresponding permuted images (shown in Figure 11m–r). Then, we can get $\mathit{PM}$.

Finally, we recover the original image from the cipher image of “Lenna” shown in Figure 12a. First, the permuted image shown in Figure 12c is obtained from the cipher image with ($\mathit{SX}$, $\mathit{SY}$, $\mathit{SZ}$). Then, the plain image is restored by $\mathit{PM}$, which is shown in Figure 12e.

Without loss of generality, we do the experiments based on other images with different sizes. The experimental results are shown in Table 3 and Figure 13. They both verify the effectiveness of our attack method. Besides, it can be seen from Table 3 that the proposed attack is efficient. Taking the image “Lenna” of size $256\times 256\times 3$ as an example, when the encryption time is 0.6391 s, the time needed for the corresponding attack is just 129.4039 s. Even if the image size increases, the time required for the attack is still within an acceptable range. Thus, it verifies that our method is computationally feasible.

Moreover, we verified the data complexity required for the attack. As discussed in Section 3.4, the total data complexity required for breaking CIEA-FOHS is $\mathit{O}(3+⌈{log}_{256}\left(3HW\right)⌉)$. In our experiment with chosen-plaintext attack, the number of attack images required for sizes $2\times 2\times 3$ and $100\times 100\times 3$ are 4 and 5, respectively. And for sizes $300\times 200\times 3$, $256\times 256\times 3$ and $512\times 512\times 3$, the number of attack images required are all 6. Therefore, the experimental verification is consistent with the theoretical calculation.

5. Suggestions for Improvement

On the basis of the above, CIEA-FOHS is insecure against a chosen-plaintext attack method because of its inherent security defects. To enhance the security, some suggestions for improvement are listed below:

• Suggestion 1. Ensuring the substantial security contribution of the fractional-order chaos to the corresponding cipher. The attractor phase diagram of the fractional-order hyperchaotic system is shown in Figure 1, which shows the extremely complex dynamics. Undoubtedly, fractional-order chaos is one of the preferred sources of entropy for encryption. However, due to the negligence of algorithm design, CIEA-FOHS has serious security defects and is attacked.
• Suggestion 2. Security analysis should be implemented from the perspective of cryptography, not limited to numerical statistical verification. As Ref. [45] points out, many encryption algorithms have excellent statistical analysis results, but they are still insecure. In fact, good statistical analysis results are only a necessary and not a sufficient condition for security. Some security flaws are difficult to reflect with numerical statistical results, but they can be clearly revealed by theoretical security analysis. For example, the existence of an equivalent key makes CIEA-FOHS vulnerable to cryptographic attacks. Given the implementation of detailed cryptographic security analysis, these flaws can be avoided, thereby improving security.

6. Conclusions

In this paper, a detailed security analysis of a color image encryption algorithm named CIEA-FOHS using a fractional-order chaos was performed. From the perspective of cryptanalysis, this paper found that CIEA-FOHS can be broken by a chosen-plaintext attack method, owing to its some inherent security defects. Theoretical analysis and experimental simulations show that the attack method is both effective and efficient for attacking CIEA-FOHS. Although the fractional-order chaotic system has complex dynamics, the algorithm defects may cause insecurity. The reported results would help the designers of chaotic cryptography pay more attention to the gap between complex chaotic system and secure cryptosystem.