Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 6797386, 12 pages
https://doi.org/10.1155/2018/6797386

Research Article

Secure Image Encryption Based on an Ideal New Nonlinear Discrete Dynamical System

R. M. Lin and T. Y. Ng

School of Mechanical and Aerospace Engineering, Nanyang Technological University, Singapore 639798

Correspondence should be addressed to R. M. Lin; mrmlin@ntu.edu.sg

Received 2 November 2017; Accepted 26 March 2018; Published 3 May 2018

Academic Editor: Bruno G. M. Robert

Copyright © 2018 R. M. Lin and T. Y. Ng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

It is well known that encryption algorithms developed based on Logistic map suffer from limited key space due to the narrow regions of system parameters which can be used, potential risk of security in the presence of numerous periodic windows within the key space, and weakness in known-plain-text attack due to the inherent correlation among the chaotic sequence used for encryption. To overcome these existing problems, this paper presents a secure image encryption algorithm based on a new highly nonlinear discrete dynamical system with ideal chaotic characteristics. Transcendental functions are introduced together with modulo operations which physically represent discontinuous time-varying nonlinearities, leading to extremely complex chaotic behavior that is highly sensitive to system parameters and initial conditions, both of which are considered as the key for the cryptosystem. Extensive numerical experiment results have revealed that the proposed image encryption algorithm offers advantages of unlimited key space and high-level security, since those problematic periodic windows are no longer present within the key space, and it is extremely robust against known-plain-text attack, since the chaotic sequence generated bears no correlation whatsoever due to the folding effect of modulo operation. The algorithm makes truly efficient yet highly secure image encryption based on chaotic systems a reality.

1. Introduction

Chaos is an inherent dynamic behavior associated with many nonlinear systems called chaotic systems [1]. These systems have some very important and unique characteristics such as extreme sensitivities of system responses to system parameters and initial conditions, leading to unpredictable system behavior, and periodic doublings, leading to nonperiodicity and seeming randomness. These complex dynamic properties are exactly what are desired in cryptography [2] since when they are used in data encryption, they serve to achieve the required diffusion and mixing of data to enhance security. These observations have recently led to the development of a number of data encryption algorithms based on chaotic sequence, which have since found useful and practical applications.

Synchronized continuous chaotic systems have been recently employed to achieve secure communications. Most of these methods are based on the principle that message signals can be mixed with a carrier signal from a chaotic system and subsequently sent through the transmission channels. The message is then recovered through chaos synchronization and subtraction of the received masked signals. To ensure successful communications, synchronization of the chaotic carrier signals between sender and receiver becomes paramount. A Lorenz chaotic system was implemented at circuit level and was used to demonstrate that synchronization of such a system can be achieved when used in communication [3]. The possibility of synchronization coupled with chaos to be used to implement secure communication systems which is inherently error-free in contrast with the signal-masking schemes proposed in [3] was discussed in [4]. A general approach for constructing chaotic synchronized dynamical systems used for communication was examined based on a decomposition of a given system into active and passive parts [5]. For enhanced security, additive noise was added to chaotic signals to generate encryption keys [6] and an effective error function attack approach was proposed to critically evaluate the latent security of existing chaos synchronization based encryption schemes [7].

In digital communications, data encryption such as image encryption schemes has been vigorously investigated to meet the ever-growing demand for real-time secure image transmissions over the Internet and through wireless networks. Traditional data encryption methods such as the Data Encryption Standard (DES) suffer from extremely low efficiency when bulk data such as images are involved [8]. To significantly improve efficiency needed while maintaining required security, encryption schemes developed based on discrete chaotic dynamical systems have been recently sought. The earliest attempt in this direction seemed to be the chaotic encryption algorithm proposed by Matthews [9]. Based on the classical one-dimensional Logistic map, which has been known to become chaotic under some system parameters, a number of image encryption methods have been developed and have met some success. By directly using the chaotic sequence generated from a Logistic map, a pixel by pixel image encryption algorithm was designed and was found to be effective [10]. Detailed cryptanalysis of the Logistic map based algorithms was conducted by Li et al. [11] and Arroyo et al. [12], revealing some of the pitfalls and recommending some possible improvement. Logistic map was also employed to construct Block Encryption Ciphers to enhance efficiency [13]. Intertwining Logistic map was investigated by Wang and Xu [14] to overcome some of the shortcomings of the conventional Logistic map.

Higher dimensional chaotic maps have also been investigated in their potential for image encryption applications. A 2D logistic map together with DNA sequence operations was examined by Wang et al. [15] for image encryption applications. A 2D Bakers map was first generalized and then discretized to create a symmetric block encryption scheme [16] and its detailed security analysis was carried out by Solak et al. [17]. A 2D Arnold cat map was examined and found to have desirable characteristics for image encryption [18]. The map was further generalized and extended to three dimensions (3D), which was then effectively employed to encrypt image data [19]. Attempts have been undertaken to couple two different chaotic systems to further enhance security [20, 21]. Cryptanalysis of these encryption schemes based on composition maps was accomplished by Li et al. [22]. Such coupling produces additional security feature similar to those of high dimensional chaotic maps. Further, general spatiotemporal chaotic systems have been proposed to achieve better trade-off between level of security and efficiency [23], as well as chaotic maps with local and global feedback and multilevel perturbations [24]. In addition to applications of chaotic cryptography to image encryptions, further researches on applications of joint image encryption with compression [25] as well as video encryptions [26, 27] have been recently carried out.

Much progress has been made to date on image encryption based on chaotic systems. However, most existing methods suffer from some or all of these problems: (i) limited key space since chaotic behavior only occurs under certain system parameter settings; (ii) presence of periodic behavior within key space which undermines security; and (iii) weakness in any possible known-plain-text attack due to the known correlation among the chaotic sequence used. The data sequences generated based on most chaotic systems within stable strange attractors are chaotic but not completely random, while ideal encryption security requires truly random data sequences. To overcome these existing problems, this paper presents a new highly nonlinear discrete dynamical system with ideal chaotic characteristics for secure image encryption. Transcendental functions are introduced together with modulo operations which physically represent discontinuous time varying nonlinearities, leading to extremely complex chaotic behavior which is highly sensitive to system parameters and initial conditions, both of which can be considered as the keys for the cryptosystem. Extensive numerical experiment results have revealed that proposed image encryption algorithm offers advantages of unlimited key space and high-level security, since those problematic periodic windows are no longer present within the key space, and it is extremely robust against known-plain-text attack, since the chaotic sequence generated bears no correlation whatsoever due to the folding effect of modulo operation. The algorithm is believed to make truly efficient yet highly secure image encryption a reality.

2. A New Highly Nonlinear Discrete Dynamical System

The proposed discrete dynamical system and its application to encryption can be considered as an improved Logistic map [28] into which a new transcendental function and a new floating modulo operation are introduced. Mathematically, the system is simply defined as

Without the modulo operation, the thus defined new map is very similar to the classical Logistic map except that, in this case, a transcendental sine function is used instead of a linear function. As a result, the proposed new system essentially becomes more nonlinear and possesses, as expected, much wider ranges of system parameters in which complex desirable chaotic behavior is observed. Such increases in the complexity of system dynamic chaotic behavior as well as in system parameter domain over which complex chaotic behavior occurs are ideal for effective, efficient, and more secure data encryption applications.

2.1. Improved Performance over Logistic Map

To compare these characteristics, bifurcation diagrams and Lyapunov exponents for both the Logistic map and the new system have been computed with all possible system parameter λ values before divergence occurs and these are shown in Figures 1–4.

Figure 1: Bifurcation diagram of Logistic map.

Figure 2: Lyapunov exponent of Logistic map.

Figure 3: Bifurcation diagram of the new system.

Figure 4: Lyapunov exponent of the new system.

The maximum and minimum values of the proposed new system without the modulo operation for every given value of λ can be determined mathematically. Since , the derivative becomes . Then, by letting the derivative be zero, , the values in the range become and , corresponding to a maximum value of and a minimum value of .

2.2. Dynamic Characteristics of the Proposed Discrete Dynamical System

It can be seen clearly that the new proposed system has a much wider range of system parameters over which the Lyapunov exponents are positive and the magnitudes of these positive Lyapunov exponents are much larger than those observed from the conventional Logistic map. These indicate that the chaotic behavior of the proposed discrete dynamical system can be far more complex than that of the conventional Logistic map. From data encryption point of view, while the new proposed system certainly offers much larger key space and better security, there still exist finite windows of system parameters in which periodic behaviors are present, as in the case of the conventional Logistic map, and these may present potential security risk if the keys chosen happen to be within these windows. For a seemingly continuous system parameter region of chaotic behavior, as shown in Figures 1 and 3, brief periodic windows are yet observed under much refined resolutions as shown in Figures 5 and 6.

(a) 3.974 ≤ λ ≤ 3.996

(b) 3.9845 ≤ λ ≤ 3.9852

Figure 5: Logistic map periodic windows under refined resolutions.

(a) 7.4 ≤ λ ≤ 7.8

(b) 7.6511 ≤ λ ≤ 7.6522

Figure 6: Periodic windows under refined resolutions of the new system.

Many more such periodic windows can perhaps be found if yet better resolutions are employed during numerical simulations. In fact, due to the fractal nature of all chaotic systems, it can be said that once such periodic windows coexist among largely chaotic behavior, there are infinitely many! This is easily seen from Figure 6(b), where periodic windows are interwoven infinitesimally close with chaotic ones with fractal boundaries between them. As a result, one shall obtain very similar behavior no matter how many times the boundaries are magnified as shown in Figure 7, where 100 and 10,000 times are used to magnify the region shown in Figure 6(b).

(a) Magnified by 100x

(b) Magnified by 10000x

Figure 7: Fractal boundaries of periodic and chaotic behavior.

To overcome such difficulties when applied to image encryption, modulo operation is further introduced to possibly eliminate these periodic windows to further enhance security. A modulo operation, from system dynamics point of view, can be considered as discontinuity and time-varying nonlinearities which are the strongest nonlinear physical mechanisms possibly encountered in practice. Introduction of such additional nonlinear mechanism can be expected to lead to improved chaos characteristics desirable for image encryption applications. A modulo 1 operation is used as defined in (1) since each time it folds data back to the desired region of . The bifurcation map and Lyapunov exponents of the complete new system are shown in Figures 8 and 9. Due to the very large system parameter range (λ = 0–10¹⁶) in this case and in different regions, characteristically distinctive dynamic behavior has been observed; the bifurcation map is presented in interval format to better illustrate the key performance of the system.

(a)λ = 0–20 (fixed points and bifurcations and periodic windows)

(b)λ = 20–200 (periodic windows increasingly becoming smaller)

(c)λ = 200–10¹⁰ (ideally chaotic and uniformly random working region)

(d)λ = 10¹⁰–10¹³ (gradual return to periodic behavior)

(e)λ = 10¹³–10¹⁶ (periodic behavior)

Figure 8: Complete bifurcation diagram of the new proposed system.

Figure 9: Lyapunov exponents of the new proposed system.

As shown in Figure 8, for small values of system parameter λ, the diagram is the same as that of Figure 3, since the modulo operation bears no effect due to small values of . As λ increases, bifurcations and periodic windows are observed but the physical widths of these windows become shorter. As λ further increases, the visible periodic windows become narrower and narrower and completely disappear when λ reaches around . Then a very wide range of ideally chaotic and uniformly random working region is obtained. This region is called working region, since only data generated for system parameters in the region will be used for subsequent image encryption. This is further supported by the Lyapunov exponents obtained as shown in Figure 9, in which the Lyapunov exponents increase almost linearly with until λ = 10¹⁰ beyond which the Lyapunov exponent start to fluctuate, indicating occurrence of possible periodic behavior as shown in Figure 8(d), even though the values are largely positive due to the numerical algorithm used. It is worth mentioning that exactly the same numbers of data points are plotted in Figures 8(c) and 8(d), but the latter appears to be containing far less data points due to the overlapping of many of the data points pertaining to periodic or near periodic behavior. Further increasing λ leads to gradual return to more periodic behavior (Figure 8) and when , various periodic responses have been observed with period as high as 28. As a result, very large λ values () are not recommended. In fact, extensive numerical simulations have shown that when , no periodic windows have been observed! This may be because, for these values of λ, each iteration produces a remote data point (far from ) that is basically the previous local value (within ) magnified λ times. When the modulo operation is applied, that data point is then being folded back to the desired interval . As a result, λ becomes essentially the stretching index that measures how far a data point is brought back to through the folding operation of the modulo function which is repeated throughout the iterations. Such repeated applications of folding processes lead to extremely complex dynamic behavior, in which periodic windows are totally removed, an ideal chaotic characteristic for secure data encryption applications. Typical iteration results for some large value of λ within the working region () are shown in Figure 10; an almost uniformly distributed random data sequence is obtained.

(a) Iteration results

(b) Data distribution histogram

Figure 10: Iteration map and data distribution for .

Another important characteristic of the new system which is very important for secure image encryption is that the chaotic data sequences generated for large values of λ within the working region are totally uncorrelated, as shown in Figure 11(b), making the inverse recovery of the sequences given some known values such as known-plain-text in the case of encryption impossible to realize. For Logistic map, however, the chaotic data sequences are perfectly correlated as shown in Figure 11(a); only few known data points are required before the whole map is mathematically recovered!

(a) Logistic map ()

(b) New complete system ()

Figure 11: Data correlations of Logistic map and the new complete system.

As compared with the existing Logistic map, the new complete system has large continuous key space and more desirable chaotic characteristics. Image encryption based on the new complete system will certainly become more secure, since it eliminates the potential security risk due to the existence of periodic windows and removes the known data correlation which significantly undermines security.

3. New Encryption Algorithm

An image encryption algorithm has been developed based on the proposed new nonlinear discrete system. Chaotic data sequence can be generated using the iteration formula shown in (1) for some suitable large value of λ within the working region () and some chosen initial condition , both of which can be considered as the keys of encryption algorithm. The original set of chaotic data sequence are decimal fractions; however, all images are specified as a grid of integer values. A simple transformation is needed to turn the chaotic data sequence into another sequence of integers. The plain image is then encrypted using a simple XOR operation with the chaotic integer sequence generated using the proposed new complete system. The block diagram of the encryption algorithm is shown in Figure 12.

Figure 12: Block diagram of the proposed encryption algorithm.

The complete encryption/decryption algorithm can be summarized in the following steps:

(1)Set encryption keys for the plain image which are the values for λ and .(2)Perform some initial iterations to eliminate any harmful effect of transient dynamics of the system (50 is used in all simulation presented in this paper).(3)Iterate one more time to generate a decimal fraction which is in double precision and choose the first 9 most significant digits to obtain 3 integers of 3 digits each.(4)Three bytes of chaotic sequence are generated by performing modulo 256 operations on the 3 integers and subsequently they are XORed with 3 bytes of plain-image data (grey value or color RGB values) to form 3 bytes of encrypted image data.(5)Repeat steps and until all plain-image data are encrypted and then output the completely encrypted object image data.(6)Send the encrypted object image data through public communication channel.(7)Pass the encryption keys through secure communication channel.(8)End.

The decryption algorithm is similar to the encryption algorithm except that the same given encryption keys should be used and the input image data are those of the ciphered image. It is noted that only the first 9 most significant digits are extracted in a double precision decimal fraction to ensure absolute safety in accuracy across different computing architectures.

4. Experimental Analysis

The proposed new encryption algorithm has been applied to an amount of image data and very encouraging results have been achieved. Experimental analysis of the standard-color Lena plain 512 × 512 image shown in Figure 13(a) has been carried out and the results are presented. Figure 13(b) shows the encrypted cipher image with chosen encryption keys , indicating that the original image has been completely destroyed and randomized. On the other hand, if the cipher image is then decrypted using the same keys, the exact original image is discovered as shown in Figure 13(c). However, if wrong keys are used, unexpected randomized decrypted images are obtained as shown in Figure 13(d), in which only very slightly different keys are used, demonstrating the extreme sensitivity of the process to system parameter λ. Similar sensitivity to initial condition has also been observed throughout the cases examined.

(a) Plain image

(b) Encrypted image

(c) Decrypted image with correct keys

(d) Decrypted image with wrong keys

Figure 13: Image encryption and decryption results.

Statistical analysis of the Lena image and its encrypted counterpart has been examined and their histograms are shown in Figures 14 and 15. The fact that almost uniform random distributions of the encrypted histograms in all 3 color components are achieved demonstrates that the proposed algorithm is capable of covering up all characters of the plain image and has very good performance of balanced 0-1 ratios, zero correlation, and enhanced security.

Figure 14: Histograms of the plain image.

Figure 15: Histograms of encrypted cipher image.

To further examine the statistical correlation between adjacent pixels in plain and ciphered images, coherence values have been computed, which is defined as

where are the RGB values of the th pair of 2 adjacent pixels, is the number of such pairs formed, and and are the mean values of and (). By randomly selecting 10000 () such pairs in the horizontal, vertical, and diagonal directions, the coherence values for the RGB components of the plain and ciphered images are computed and are shown in Table 1. One typical such pair of data in the horizontal direction for the red color components of the plain and ciphered images is plotted as shown in Figure 16. As can be seen, the correlation that exists in the plain image is completely destroyed in its encrypted counterpart when the encryption algorithm is applied. This is further illustrated by the coherence values obtained, indicating that the encryption algorithm satisfies zero correlation.

Table 1: Coherence values of adjacent pixels of plain and ciphered images.

Plain image
Ciphered image

Red
Green
Blue
Red
Green
Blue

Horizontal

Vertical

Diagonal

(a) Plain image (red),

(b) Ciphered image (red),

Figure 16: Correlation of horizontally adjacent pixels of plain and ciphered images.

Computationally, the proposed encryption/decryption algorithm is very fast and efficient. The notebook used for this test is a 2.4 GHz Fujitsu Lifebook with an Intel® Core™ i5 CPU and a 2.99 GB RAM. The total CPU time clocked for the encryption and the subsequent decryption of the Lena color image of size 512 × 512 is 0.40221 seconds, which is much shorter than those reported [10].

5. Security Analysis

The new encryption algorithm developed in this paper based on the new highly nonlinear discrete dynamical system belongs to one-time one-password system that can be considered as secure according to the requirements of information communication system discussed by Shannon [2]. Compared with existing encryption algorithms based on classical Logistic map, the new algorithm has large continuous key space by eliminating the periodic windows that undermine security, and, as shown in Figure 11(b), the chaotic sequences generated are totally uncorrelated, thereby making their reverse recovery totally hopeless! These have essentially resolved all existing security weaknesses pertaining to chaos based image encryption methods and have opened a new era for their ever-increasing practical applications.

Extensive experimental analyses have shown that chaotic data sequences generated by the proposed system are truly random, which have an ideal property of uniform distribution. When such data sequences are employed, trial tests have shown that the proposed encryption algorithm possesses many ideal characteristics of traditional cryptography, such as highly complex system behavior, ideal discontinuous and time-varying nonlinearity due to the introduction of modulo operation, zero correlation, and the irreversibility of the process of generating chaotic sequences. Because of all these desirable properties, the algorithm is believed to be able to resist grey code attack and statistical attack strongly [4].

Brute-force attack is becoming increasingly popular nowadays and the strength of a cryptosystem under such attack is mainly measured by its keys. The proposed algorithm is capable of resisting Brute-force attack, since its keys consist of 2 floating-point numbers . If one uses the first 15 digits of a floating-point number in double precision, then the key will have altogether digits. This represents a total possible number of keys of 10³², which is far greater than 2⁵⁶ of the most commonly used 56-bit DES [29], a publicly acknowledged secure encryption algorithm. Further, the key length of the proposed algorithm can be easily further increased. Instead of (1), the more general system can be written as

As a result, 5 system parameters, λ, , μ (), α, and β (both are short integers), can be used to form the keys. The key in this case will have digits, in which the two short integers contribute roughly about 5 digits, leading to a possible number of keys of , which is believed to be beyond the reach of the most sophisticated attacks! One typical bifurcation diagram for the more general system with μ = 0.75, α = 2, and β = 3 is shown in Figure 17 and, as expected, similar chaotic characteristics are observed.

(a)λ = 0–150 (fixed points and bifurcations and periodic windows)

(b)λ = 150–1000 (periodic windows increasingly becoming smaller)

(c)λ = 1000–10¹⁰ (ideally chaotic and uniformly random working region)

(d)λ = 10¹⁰–10¹³ (gradual return to periodic behavior)

(e)λ = 10¹³–10¹⁶ (periodic behavior)

Figure 17: Bifurcation diagram of the more general system.

Most chaotic encryption algorithms are very vulnerable to known-plain-text attack due to the correlation that exists among the values of the chaotic sequence, as shown in Figure 11(a). In the case of encryption based on the Logistic map, theoretically, only 2 pixels of known-plain-text are needed to compute the hidden keys. This is a fatal weakness associated with most of the existing chaotic encryption algorithms, as far as their practical applications are concerned. The new proposed algorithm, however, has a perfect ability to resist such an attack, since there is no such correlation whatsoever among the chaotic sequence, as shown in Figure 11(b), due to the application of the modulo operation. For the case of possible possession of any known plain text, no matter how much that text may be, it is not possible to use these texts to further increase the knowledge of the hidden plain image or to use them to possibly recover the hidden keys. The keys and the remaining portion of the plain image are not in any sense undermined due to the availability of some known-plain-text, since the encryption process is irreversible unless the keys are known. In addition, only partial digits of the values of the chaotic sequence are used during encryption and this adds further mathematical complexity to that already irreversible encryption process. This enhanced security is one of the most attractive features of the new algorithm.

6. Conclusions

In this paper, a new highly nonlinear discrete dynamical system has been developed and investigated. Chaotic characteristics of the system such as bifurcations, fractal boundaries, and statistical properties of chaotic data sequence have been analyzed and have been found to be the most desirable for data encryption applications. The introductions of a transcendental sine function and a modulo operation render the system behavior extremely complex, leading to the removal of the periodic windows for large values of system parameters, the elimination of the correlation among the data sequences generated, and an almost ideally uniformly distributed random data series. The modulo function performs repeated folding operations, thereby generating truly uncorrelated random data. Subsequently, a new image encryption algorithm has been developed, which overcomes these problems associated with existing encryption algorithms based on chaotic dynamics: (i) limited key space since chaotic behavior only occurs under limited system parameters; (ii) presence of periodic behavior within key space which undermines security; and (iii) weakness in any possible known-plain-text attack due to the known correlation among the chaotic sequences. Extensive numerical experiments have revealed that the proposed image encryption algorithm offers advantages of unlimited key space and high-level security and is extremely robust against known-plain-text attacks. Coupled with system simplicity and computational efficiency, it is believed that the proposed algorithm will make truly efficient yet highly secure image encryptions a reality.

Conflicts of Interest

The authors declare that they have no conflict of interests.

References

J. Gleick, Chaos-Making a New Science, Penguin Books USA, 1988.

C. E. Shannon, “A mathematical theory of communication,” Bell Labs Technical Journal, vol. 27, pp. 379–423, 1948.

K. M. Cuomo and A. V. Oppenheim, “Circuit implementation of synchronized chaos with applications to communications,” Physical Review Letters, vol. 71, no. 1, pp. 65–68, 1993.

C. W. Wu and L. O. Chua, “A simple way to synchronize chaotic systems with applications to secure communication systems,” International Journal of Bifurcation and Chaos, vol. 3, no. 6, pp. 1619–1627, 1993.

L. Kocarev and U. Parlitz, “General approach for chaotic synchronization with applications to communication,” Physical Review Letters, vol. 74, no. 25, pp. 5028–5031, 1995.

A. A. Minai and T. D. Pandian, “Communicating with noise: how chaos and noise combine to generate secure encryption keys,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 8, no. 3, pp. 621–628, 1998.

X. Wang, M. Zhan, C.-H. Lai, and H. Gang, “Error function attack of chaos synchronization based encryption schemes,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 14, no. 1, pp. 128–137, 2004.

G. Chen, Y. Mao, and C. K. Chui, “A symmetric image encryption scheme based on 3D chaotic cat maps,” Chaos, Solitons & Fractals, vol. 21, no. 3, pp. 749–761, 2004.

R. Matthews, “On the derivation of a 'chaotic' encryption algorithm,” Cryptologia, vol. 13, no. 1, pp. 29–42, 1989.

A. N. Pisarchik, N. J. Flores-Carmona, and M. Carpio-Valadez, “Encryption and decryption of images with chaotic map lattices,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 16, no. 3, Article ID 033118, 2006.

C. Li, T. Xie, Q. Liu, and G. Cheng, “Cryptanalyzing image encryption using chaotic logistic map,” Nonlinear Dynamics, vol. 78, no. 2, pp. 1545–1551, 2014.

D. Arroyo, R. Rhouma, G. Alvarez, S. Li, and V. Fernandez, “On the security of a new image encryption scheme based on chaotic map lattices,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 18, no. 3, Article ID 033112, 2008.

G. Jakimoski and L. Kocarev, “Chaos and cryptography: block encryption ciphers based on chaotic maps,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 48, no. 2, pp. 163–169, 2001.

X. Wang and D. Xu, “Image encryption using genetic operators and intertwining logistic map,” Nonlinear Dynamics, vol. 78, no. 4, pp. 2975–2984, 2014.

X.-Y. Wang, Y.-Q. Zhang, and Y.-Y. Zhao, “A novel image encryption scheme based on 2-D logistic map and {DNA} sequence operations,” Nonlinear Dynamics, vol. 82, no. 3, pp. 1269–1280, 2015.

J. Fridrich, “Symmetric ciphers based on two-dimensional chaotic maps,” International Journal of Bifurcation and Chaos, vol. 8, no. 6, pp. 1259–1284, 1998.

E. Solak, C. Çokal, O. T. Yildiz, and T. Biyikoǧlu, “Cryptanalysis of fridrich's chaotic image encryption,” International Journal of Bifurcation and Chaos, vol. 20, no. 5, pp. 1405–1413, 2010.

D. Xiao, X. Liao, and P. Wei, “Analysis and improvement of a chaos-based image encryption algorithm,” Chaos, Solitons & Fractals, vol. 40, no. 5, pp. 2191–2199, 2009.

K. Wang, W. Pei, L. Zou, A. Song, and Z. He, “On the security of 3D Cat map based symmetric image encryption scheme,” Physics Letters A, vol. 343, no. 6, pp. 432–439, 2005.

S. Mazloom and A. M. Eftekhari-Moghadam, “Color image encryption based on coupled nonlinear chaotic map,” Chaos, Solitons & Fractals, vol. 42, no. 3, pp. 1745–1754, 2009.

S. Behnia, A. Akhshani, H. Mahmodi, and A. Akhavan, “A novel algorithm for image encryption based on mixture of chaotic maps,” Chaos, Solitons & Fractals, vol. 35, no. 2, pp. 408–419, 2008.

C. Li, D. Arroyo, and K.-T. Lo, “Breaking a chaotic cryptographic scheme based on composition maps,” International Journal of Bifurcation and Chaos, vol. 20, no. 8, pp. 2561–2568, 2010.

T. Xiang, K.-W. Wong, and X. Liao, “Selective image encryption using a spatiotemporal chaotic system,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 17, no. 2, Article ID 023115, 2007.

R. Hasimoto-Beltran, “High-performance multimedia encryption system based on chaos,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 18, no. 2, Article ID 023110, 2008.

X.-J. Tong, Z. Wang, M. Zhang, and Y. Liu, “A new algorithm of the combination of image compression and encryption technology based on cross chaotic map,” Nonlinear Dynamics, vol. 72, no. 1-2, pp. 229–241, 2013.

F. Chiaraluce, L. Ciccarelli, E. Gambi, P. Pierleoni, and M. Reginelli, “A new chaotic algorithm for video encryption,” IEEE Transactions on Consumer Electronics, vol. 48, no. 4, pp. 838–844, 2002.

S. Yang and S. Sun, “Video encryption method based on chaotic maps in DCT domain,” Progress in Natural Science, vol. 18, no. 10, pp. 1299–1304, 2008.

R. M. May, “Simple mathematical models with very complicated dynamics,” Nature, vol. 261, no. 5560, pp. 459–467, 1976.

C. Paar and J. Pelzl, Understanding Cryptography, Springer, 2010.