Next Article in Journal
Development of Novel Lightweight Dual-Phase Al-Ti-Cr-Mn-V Medium-Entropy Alloys with High Strength and Ductility
Next Article in Special Issue
An Insider Data Leakage Detection Using One-Hot Encoding, Synthetic Minority Oversampling and Machine Learning Techniques
Previous Article in Journal
Energy Disaggregation Using Elastic Matching Algorithms
Previous Article in Special Issue
An Algorithm of Image Encryption Using Logistic and Two-Dimensional Chaotic Economic Maps
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Image Encryption Algorithm Based on Random Hamiltonian Path

Software College, Northeastern University, No.11, Lane 3, Wenhua Road, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(1), 73; https://doi.org/10.3390/e22010073
Submission received: 7 December 2019 / Revised: 30 December 2019 / Accepted: 3 January 2020 / Published: 6 January 2020

Abstract

:
In graph theory, Hamiltonian path refers to the path that visits each vertex exactly once. In this paper, we designed a method to generate random Hamiltonian path within digital images, which is equivalent to permutation in image encryption. By these means, building a Hamiltonian path across bit planes can shuffle the distribution of the pixel’s bits. Furthermore, a similar thought can be applied for the substitution of pixel’s grey levels. To ensure the randomness of the generated Hamiltonian path, an adjusted Bernoulli map is proposed. By adopting these novel techniques, a bit-level image encryption scheme was devised. Evaluation of simulation results proves that the proposed scheme reached fair performance. In addition, a common flaw in calculating correlation coefficients of adjacent pixels was pinpointed by us. After enhancement, correlation coefficient becomes a stricter criterion for image encryption algorithms.

1. Introduction

When computers and the internet came on the scene, here came the era of information, accompanied by the formidable challenge of information security. Among complicated information, the vivid multimedia information is preferred by people, especially digital images. Consequently, such information involves both collective interests and personal interests. For instance, images of military affairs are related to the safety of whole country. Privacy and copyright of images influence everyone’s peace of mind. To protect the rights of image’s owners, methods like steganography, watermarking, and encryption are frequently utilized [1]. Among these techniques, encryption is a direct and thorough means. Nowadays, image encryption is an inviting and fruitful field, and many imaginative image encryption algorithms are proposed.
One picture is worth more than ten thousand words, and there indeed are tens of thousands of pixels in a digital image. To encrypt the bulk data of images, traditional cryptosystems are not efficient enough. Among specific image encryption schemes, the permutation–diffusion structure is widely used. Essentially, permutation is to rearrange image pixels on different dimensions. In [2], 2D CMT (chaotic magic transform) was proposed for permutation. In [3], image scrambling was performed by a parametric 2D Sudoku matrix. In [4], horizontal and vertical wave lines were utilized to realize row rotation and column rotation. This is also a 2D method. In [5], spatial permutation was performed on a 3D bit matrix by using orthogonal Latin cubes. Moreover, file-based algorithms like [6] deem images as 1D binary files when realizing permutation. Considering the features of bit distribution in digital images, encryption schemes [7,8] with bit-level permutation are proposed.
Sometimes, the permutation phase is accompanied by a sort operation, such as [2]. However, time complexity of sorting is usually nonlinear. To obtain high efficiency, the additional operation should be avoided. Regarding an image as a 1D pixel array, permutation can be depicted as an arrangement of pixels, which is represented in the bijective map from plain image to permuted image. If we connect the pixels by order of the arrangement, all the pixels are traversed exactly once. Deeming pixels as the vertices of a graph, such a path of traversing is known as a Hamiltonian path in graph theory. Conversely, a Hamiltonian path corresponds to an arrangement of permutations. Following this thought, the method of building Hamiltonian path is equal to a permutation scheme. As a Hamiltonian path can be generated without sort operation, the corresponding permutation algorithm has the advantage of efficiency. In cryptography, substitution is a classical method of cipher schemes. To substitute pixel’s grey value, arrangement of all the possible grey levels is requisite. Hence, the thought of random Hamiltonian path is also suitable for the substitution of grey levels.
It is common knowledge that chaotic systems have conspicuous advantages for cryptosystems. High dimension chaotic systems possess complex chaotic behavior, while 1D chaotic systems are convenient for implementation. Under synthetical consideration, some combined chaotic maps have been explored in recent encryption schemes [2,9,10,11,12]. Just resembling the series-parallel connection of resistors in circuits, these chaotic maps are a combination or adjustment of the original chaotic maps. Multiple chaotic maps can be coupled as CML (chaotic map lattice) [13,14]. By these means, chaotic behavior is magnified, leading to better chaotic performance. In this paper, a new chaotic map was also explored.
There are some innovative works in this paper:
  • A method of building a random Hamiltonian path within digital images was designed, which is equivalent to permutation. On this basis, bit-level permutation of high efficiency was achieved.
  • Following the thought of the random Hamiltonian path, arrays for grey levels’ substitution can be generated.
  • An adjusted Bernoulli map is proposed, which is suitable for image encryption schemes.
  • The ambiguous definition of diagonal direction is normalized to two orthogonal directions when calculating correlation coefficients.
The rest of this paper is organized as following: Section 2 explains the Hamiltonian path and the procedures to generate such paths within images. Section 3 expounds the adjusted Bernoulli map. In Section 4 the proposed scheme is thoroughly introduced. The results of simulation experiments are exhibited in Section 5. Section 6 is the summary of the entire paper.

2. Hamiltonian Path

As was mentioned earlier, the scheme of generating a random Hamiltonian path is tantamount to permutation. In this section, the relevant theories are presented, while the method of generating a Hamiltonian path is proposed.

2.1. Basic Theory of Hamiltonian Path

Graph theory is a classical branch of mathematics. The term graph refers to the figures composed by points and the connecting lines between the points. Commonly, the points are called vertices, and the lines are called edges. The definition of graph is G = (V, E). Here V is a nonempty set of finite vertices and the set of edges E = {(x, y)|x, yV}. If the vertex pair (x, y) in E is ordered, the graph is named a direct graph. Otherwise, it is named an undirected graph.
In an undirected graph, a path P is a sequence of vertices v1 v2vk, and there exists anedge between each of the vertex pairs vi vi+1. The k is the number of vertices that P contains, in other words, the length of P.
There are two special categories of path, Euler path and Hamiltonian path. Euler path refers to the paths that traverse each edge once and only once. A famous instance is the problem of Konigsberg bridges [15]. In 1736, Leonhard Euler had proved that there is no solution for the problem. This is known as the beginning of graph theory. Hamiltonian path refers to the paths that traverse each node once and only once, or an arrangement of all vertices in which every adjacent vertex pair is connected by at least one edge. The problem of Hamiltonian path can be traced back to 1859, when Willian Hamilton talked about a mathematical game: traverse all the vertices of a dodecahedron and pass by the vertices exactly once. Figure 1 is the graphic illustration of the two famous problems.
Graphs are generally intricate. The problem of finding a Hamiltonian path is a nondeterministic polynomial complete problem (NP-C problem), one of the most burdensome challenges in mathematics [16,17,18,19]. Off the beaten track, DNA computing [20] and light-based computers [21] have been developed to solve this problem efficiently. However, generating a Hamiltonian path within digital images can be much easier.

2.2. Hamiltonian Path Within Digital Images

For an undirected graph of N vertices, there are, at most, N × (N + 1)/2 edges. Under this condition, any two vertices are connected by an edge. Such a graph is called a complete graph. Some instances of complete graphs are shown in Figure 2. The complete graphs with three nodes, four nodes, and five nodes are shown in Figure 2a–c, respectively.
There are varieties of theorems to measure whether there are Hamiltonian paths in a graph. One of the theorems is as below:
Dirac theorem: In a graph G of N vertices, if for each vertex vi there always is d(vi) ≥ N/2, then at least one Hamiltonian path exists in G.
The d(vi), otherwise called the degree of vi, represents the quantity of edges connected with vi.
In our scheme, digital images were regarded as complete graphs. Hereof, the pixels are the vertices, and there is an edge between every two pixels. According to Dirac theorem, there always exist Hamiltonian paths in such graphs.
To build a Hamiltonian path within an image, pixels are divided into two parts. One is composed by the pixels that have been added into the path. The other one is composed by the rest of the pixels. Firstly, a pixel is chosen to be the path’s outset. Then, the other pixels are added to the path one by one. If the image’s size is M × N and its pixels are {P1, P2, …, PM×N}, this progress can be generalized as the following steps:
Step 1: Choose a pixel from {P1, P2, …, PM×N} and put it in the position of PM×N.
Step 2: Choose a pixel from {P1, P2, …, PM×N-1} and put it in the position of PM×N-1.
Step 3: Choose a pixel from {P1, P2, …, PM×N-2} and put it in the position of PM×N-2.
Step M×N–1: Choose a pixel from {P1, P2} and put it in the position of P2.
In the above process, pixels that have been added into the path are insulated at the back of image’s pixel array. Among the whole image, these pixels are deemed as permutated pixels. The complete process is shown in Figure 3. Figure 4 illustrates the generated Hamiltonian path from a graph perspective.

2.3. Hamiltonian Path across Bit Planes

In [7], the intrinsic features of bit distribution in digital images were revealed. Higher bits of pixels hold higher weight of an image’s information, and there are strong correlations among the higher bit planes. In the instance of Figure 5, the 8th bit plane and the 7th bit plane tend to have opposite values. These features shall not be neglected in a secure cryptosystem.
To build a Hamiltonian path across bit planes, the strategy of [7] was extended to greyscale images in this paper. By these means, a plain image of size M × N is expanded to 2M × 2N. All the bit planes of a plain image’s pixels were placed to the 1st bit plane and the 2nd bit plane of the expanded image’s pixels. After generating a Hamiltonian path, the bit planes were restored, and a permutated image of size M × N was formed. The whole procedure can be generalized into Figure 6.
Figure 7 is the illustration of the generated Hamiltonian path across bit planes.

3. Adjusted Bernoulli Map

To ensure the randomness of the Hamiltonian path, chaotic maps can serve as pseudo random number generators. Theoretically, any 1D chaotic map is compatible. In this section, an adjusted Bernoulli map is proposed.

3.1. Bernoulli Map

The original definition of Bernoulli map [22,23] is given by:
x n + 1 = 2 x n mod 1 = { 2 x n 0 < x n < 0 . 5 2 x n 1 0 . 5 x n < 1 .
The piecewise linear property of a Bernoulli map is demonstrated in Figure 8. When implemented into discrete computer systems, the map resembles bit shifting of floating numbers. Such degradation means that the original Bernoulli map is seldom applied to encryption algorithms directly.

3.2. Adjusted Bernoulli Map

To amplify the limited nonlinear property of original Bernoulli map, cascaded modulus operations are adopted in the adjusted Bernoulli map (ABM).
x n + 1 = β ( α x n mod 1 ) mod 1 .
The parameters α and β can be many of the floating-point numbers that are bigger than two. Though the multiplication operation is linear in mathematics, the multiplication operation in computer systems involves the conversion between decimal number and binary number. The ABM possesses fair chaotic behavior in practice, especially when the parameters α and β are random. Owing to the finite precision of computers, the ABM does not work well when its parameters are big numbers, and special values such as 2N and 10N should be avoided. Here the N is the set of natural numbers. Part of the parameters’ value range is shown in Figure 9.
To examine the randomness of the pseudo-random numbers generated by the ABM, the NIST SP800-22 test suite [24] was utilized. In our experiment, 300 bitstreams of length 106 were generated and tested. The α = 10.45678 and β = 10.123 in these bitstreams. The initial value of x was increased by 0.0033, ranging from 0.001 to 0.991. The test results are listed in Table 1.

4. Proposed Scheme

After the progress of Section 2, a bit-level permutation was completed. To obtain fair diffusion properties, XOR operations were performed on pixels, and their grey values were substituted dynamically. The whole encryption scheme is detailed in this section.

4.1. Encryption Algorithm

As is illustrated in Figure 10, the whole cryptosystem is handled by ABM. The inputs of the algorithm are the plain image P of size M × N and the parameters of ABM. The output is the cipher image C.
The whole encryption process is as below:
Step 1: Read in P. Iterate ABM to avoid transient effect.
Step 2: Decompose P’s bit planes. Make a montage of these bit planes to obtain an image B of size 2M × 2N.
Step 3: For i = 2M × 2N, 2M × 2N − 1, 2M × 2N − 2, …, 3, 2, iterate ABM to generate pseudo random number ri and use Equation (3) to quantize it. Swap M’s ith pixel Mi and jth pixel Mj.
j = r o u n d ( r i × 10 14 ) mod i + 1 .
Step 4: Merge the decomposed bit planes to obtain the permutated image H of size M × N.
Step 5: Initialize two 1D arrays or vectors S and T by Equation (4). For i = 0, 1, 2, …, 255,
S i = T i = i .
Step 6: For i = 255, 254, …, 2, 1, iterate ABM to generate pseudo random number ui and use Equation (5) to quantize it. Swap Si and Sj.
j = r o u n d ( u i × 10 14 ) mod i .
Step 7: For i = 255, 254, …, 2, 1, iterate ABM to generate pseudo random number pi and use Equation (6) to quantize it. Swap Ti and Tj.
j = r o u n d ( p i × 10 14 ) mod i .
Step 8: For i = 1, 2, …, M × N − 1, use Equation (7) to diffuse H’s pixel Hi+1. Here, ai is the pseudo random number generated by ABM.
H i + 1 = H i + 1 S T H i ( r o u n d ( a i × 10 14 ) mod 256 ) .
Step 9: For i = M × N, M × N − 1, …, 3, 2, use Equation (8) to diffuse H’s pixel Hi−1. Here, bi is the pseudo random number generated by ABM.
H i 1 = H i 1 T S H i ( r o u n d ( b i × 10 14 ) mod 256 ) .
Step 10: Save H as the C.

4.2. Discussion

In Step 2, Step 3, and Step 4, the permutation phase that works on a bit-level was performed. According to the method of building a Hamiltonian path, two arrays were generated for grey value’s substitution in Step 5, 6, and 7. The arrays were arrangements of integers from 0 to 255, in accordance with pixel’s grey levels. As the arrays were randomly generated, there were 256! ≈ 8.578 × 10506 possible arrangements. In this way, the modification of plain images could be amplified and transmitted in Step 8 and Step 9, causing an avalanche effect.

4.3. Decryption Algorithm

The decryption algorithm is the reverse progress of the encryption algorithm, as can be seen from Figure 11.
In the encryption process, the substitution is realized by arrays S and T. The reverse operation of substitution needs the inverse map of S and T, which can be generated by Equation (9).
S S i = T T i = i .
In the formula, i = 255, 254, …, 2, 1, 0. The S’ and T’ are the inverse maps of S and T, respectively.

5. Simulation Experiments

To check the performance of the proposed scheme, the results of simulation experiments were evaluated by several criteria in this section. Our experimental environment was a desktop PC with 64-bit Windows 10 OS, Intel i7-2600 CPU, and 8GB RAM. The programming language was C++, and the developing environments were Visual Studio 2019 and OpenCV 4.1.0. The test images were chosen from SIPI image database [25].

5.1. Secret Key Analysis

Secret key is an indispensable component of a cryptosystem. The key space is suggested to be no less than 2100 [26]. The secret keys of the proposed scheme are parameters of ABM. In our simulation experiment, the data type of the keys was double precision floating point numbers. According to IEEE 754 standard, each key occupies 8 Bytes and owns significant digit of 52 bits. The structure of secret key is as shown in Figure 12, and the key space is bigger than 2100.
To examine the key sensitivity in the encryption process and decryption process, a strict test for bit change rate—NBCR (the number of bit change rate) [27]—was utilized:
N B C R ( C 1 , C 2 ) = H a m ( C 1 , C 2 ) M × N × d
In the above formula, C1 and C2 are two images of size M × N and bit-depth d. Ham(C1,C2) represents the Hamming distance between C1 and C2; in other words, the number of different bits between the two images. The calculation results of NBCR should be close to 50%, which indicates that around 50% of the bits are different between C1 and C2.
In our work, three groups of modified keys were set as the illegal keys. These illegal keys are utilized to encrypt the plain images in the encryption process and decrypt the cipher images in the decryption process. The obtained encrypted images and decrypted images were made in comparison with the original plain images and cipher images. The NBCRs are listed in Table 2.

5.2. Histograms

The histogram is the foundation of various spatial image processing techniques, e.g., image enhancement. Moreover, the inherent information of histograms is useful in image compression and segmentation. For an image of size M × N and bit-depth d, the histogram is a discrete function:
h ( r i ) = q i .
Here i = 0, 1, 2, …, 2d − 1, qi is the pixels’ quantity of grey value ri. The variance of histogram can be calculated by Equation (12).
V a r ( h i ) = 1 2 d i = 0 2 d 1 ( q i μ h ) 2 .
The μ h is the arithmetic mean value of qi. The histogram of a cipher image should be relatively uniform. After encryption, the variance of an image’s histogram should be reduced. In Table 3, the variance of several images’ histograms are listed.
In practice, histograms are often normalized by Equation (13).
p ( r i ) = q i M × N .
After normalization, p(ri) represents the emergence probability of the ith grey value. The normalized histograms of plain images and cipher images are shown in Figure 13.

5.3. Information Entropy

Information entropy was proposed by C. E. Shannon, which is a measurement of information’s randomness. For a digital image, it is hard to predict the content if its information entropy is high. The calculation formula of information entropy is as shown in Equation (14). Here, the p(ri) is identical to Equation (13).
H = i = 0 2 d 1 p ( r i ) log 2 p ( r i ) .
The ideal value of a cipher image’s information entropy is its bit-depth d. In Table 4, the information entropy of plain images and cipher images are listed.

5.4. Differential Attack

To resist differential attacks, tiny modification in plain images should cause massive changes in the cipher image. This is known as diffusion property in cryptography. NPCR (number of pixel change rate) and UACI (unified averaged changed intensity) are two common indicators for an algorithm’s ability of resisting differential attacks [29]. If C1 and C2 are two images of size M × N and bit-depth d, then
N P C R = [ 1 M × N i = 1 M j = 1 N D ( i , j ) ] × 100 %
U A C I = [ 1 M × N × 255 i = 1 M j = 1 N | C 1 ( i , j ) C 2 ( i , j ) | ] × 100 % .
Here,
D ( i , j ) = { 0 , C 1 ( i , j ) = C 2 ( i , j ) 1 , C 1 ( i , j ) C 2 ( i , j ) .
In our experiment, the plain image boat of size 512 × 512 was utilized for evaluating the diffusion effect. Some pixels were chosen in the image, and the last bit of these pixels were reversed, respectively. Then, the modified images were encrypted. As can be seen from Table 5, the NPCRs and UACIs were close to theoretical values after two encryption rounds.

5.5. Correlation Coefficients

Plain images usually are redundant in the spatial domain, which means that adjacent pixels are highly correlated. Whereas, in cipher images, such a correlation should be broken. To measure the correlation between adjacent pixels, we calculated correlation coefficients as below:
r x y = E [ ( x μ x ) ( y μ y ) ] σ x σ y .
The x and y are pixel vectors of the same length. The μx and μy are their arithmetic mean values, and the σx and σy are their standard deviations. The range of correlation coefficients is [−1, 1]. If x and y are not correlated, rxy shall be close to 0.
Commonly, the adjacent pixels of three directions are calculated, respectively horizontal, vertical, and diagonal. Whereas, there exist two orthogonal diagonal directions in 2D matrices of pixels—the principal diagonal direction (from upper-left to lower-right) and the minor diagonal direction (from upper-right to lower-left). For instance, in a pixel block [ p 1 p 2 p 3 p 4 ] , p1 and p4 are adjacent in the principal diagonal direction, while p2 and p3 are adjacent in the minor diagonal direction. In the field of image encryption, the definition of diagonal direction is usually ambiguous. However, the two diagonal directions are nonequivalent for some image processing techniques and image encryption algorithms [30,31,32,33,34]. Under the extreme circumstances in Figure 14 and Figure 15, it is not enough to calculate only three of the four directions.
For all the plain images and cipher images in our experiments, correlation coefficients of 10,000 adjacent pixel pairs in each of the four directions were calculated. The results are listed in Table 6.

5.6. Efficiency

In the proposed scheme, bit-level permutation is performed in linear time complexity. Meanwhile, the diffusion phase is also linear. If the encrypted image is of size M × N, the algorithm’s time complexity is O(MN). The time complexity of the algorithm in [4] is also O(MN). However, our bit-level scheme is slower than the pixel-level scheme of [4]. In [2], permutation was companied by sorting operation. Thus, the scheme’s efficiency was related to the adopted sorting algorithm. In [28], the algorithm’s time complexity is O(MN(M + N)). The comparison between these algorithms’ efficiency is presented in Table 7.

6. Conclusions

In this paper, a 1D adjusted Bernoulli map is proposed, which is suitable for encryption systems. Based on the new chaotic map, an innovative image encryption algorithm was designed. The permutation phase was realized by generating a random Hamiltonian path, which was performed across different bit planes. Then, the idea of random Hamiltonian path was extended for substitution of grey levels in the diffusion phase. Various criterions indicate that our scheme had a pretty good performance. Besides, for measuring the correlation of adjacent pixels more reasonably, both the principal diagonal direction and the minor diagonal direction are involved when calculating correlation coefficients.

Author Contributions

Project administration, W.Z.; writing—original draft preparation, S.W.; writing—review and editing, W.H.; supervision, H.Y.; Project Administration, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The work described in this paper was supported by the National Natural Science Foundation of China (Grant No. 61402092, 61374178, 61603182), the Fundamental Research Funds for the Central Universities (N171704004).

Acknowledgments

We thank the anonymous reviewers for their helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, M.; Shah, T. A literature review on image encryption techniques. 3D Res. 2014, 5, 29. [Google Scholar] [CrossRef]
  2. Hua, Z.Y.; Zhou, Y.C.; Pun, C.M.; Chen, C.L.P. 2D sine logistic modulation map for image encryption. Inf. Sci. 2015, 297, 80–94. [Google Scholar] [CrossRef]
  3. Wu, Y.; Zhou, Y.; Agaian, S.; Noonan, J.P. 2D Sudoku associated bijections for image scrambling. Inf. Sci. 2016, 327, 91–109. [Google Scholar] [CrossRef]
  4. Ye, G.; Zhao, H.; Chai, H. Chaotic image encryption algorithm using wave-line permutation and block diffusion. Nonlinear Dyn. 2016, 83, 2067–2077. [Google Scholar] [CrossRef]
  5. Xu, M.; Tian, Z. A novel image cipher based on 3D bit matrix and latin cubes. Inf. Sci. 2019, 478, 1–14. [Google Scholar] [CrossRef]
  6. Diaconu, A.V. Circular inter-intra bit-level permutation and chaos-based image encryption. Inf. Sci. 2016, 355, 314–327. [Google Scholar] [CrossRef]
  7. Zhang, W.; Wong, K.; Yu, H.; Zhu, Z. A symmetric color image encryption algorithm using the intrinsic features of bit distributions. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 584–600. [Google Scholar] [CrossRef]
  8. Cao, C.; Sun, K.; Liu, W. A novel bit-level image encryption algorithm based on 2D-LICM hyperchaotic map. Signal Process. 2018, 143, 122–133. [Google Scholar] [CrossRef]
  9. Hua, Z.; Zhou, Y.; Chen, C.L.P. A new series-wound framework for generating 1D chaotic maps. In Proceedings of the 2013 IEEE Digital Signal Processing and Signal Processing Education Meeting, Napa, CA, USA, 11–14 August 2013; pp. 118–123. [Google Scholar]
  10. Lan, R.; He, J.; Wang, S.; Gu, T.; Luo, X. Integrated chaotic systems for image encryption. Signal Process. 2018, 147, 133–145. [Google Scholar] [CrossRef]
  11. Pak, C.; Huang, L. A new color image encryption using combination of the 1D chaotic map. Signal Process. 2017, 138, 129–137. [Google Scholar] [CrossRef]
  12. Zhou, Y.; Bao, L.; Chen, C.L.P. Image encryption using a new parametric switching chaotic system. Signal Process. 2013, 93, 3039–3052. [Google Scholar] [CrossRef]
  13. Chapaneri, S.; Chapaneri, R.; Sarode, T. Evaluation of Chaotic Map Lattice Systems for Image Encryption. In Proceedings of the 2014 International Conference on Circuits, Systems, Communication and Information Technology Applications (CSCITA), Mumbai, India, 4–5 April 2014; pp. 59–64. [Google Scholar]
  14. Tang, Y.; Wang, Z.; Fang, J. Image encryption using chaotic coupled map lattices with time-varying delays. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 2456–2468. [Google Scholar] [CrossRef]
  15. Bellman, R.; Cooke, K.L. The Konigsberg bridges problem generalized. J. Math. Anal. Appl. 1969, 25, 1–7. [Google Scholar] [CrossRef] [Green Version]
  16. Bax, E.T. Inclusion and exclusion algorithm for the Hamiltonian path problem. Inf. Process. Lett. 1993, 47, 203–207. [Google Scholar] [CrossRef]
  17. Bertossi, A.A. The edge Hamiltonian path problem is NP-complete. Inf. Process. Lett. 1981, 13, 157–159. [Google Scholar] [CrossRef]
  18. Conrad, A.; Hindrichs, T.; Morsy, H.; Wegener, I. Solution of the knight’s Hamiltonian path problem on chessboards. Discret. Appl. Math. 1994, 50, 125–134. [Google Scholar] [CrossRef] [Green Version]
  19. Schiermeyer, I. Problems remaining NP-complete for sparse or dense graphs. Discuss. Math. Gr. Theory 1995, 15, 33–41. [Google Scholar] [CrossRef]
  20. Baumgardner, J.; Acker, K.; Adefuye, O.; Crowley, T.S.; DeLoache, W.; Dickson, J.O.; Heard, L.; Martens, A.; Morton, N.; Ritter, M.; et al. Solving a hamiltonian path problem with a bacterial computer. J. Biol. Eng. 2009, 3, 11. [Google Scholar] [CrossRef] [Green Version]
  21. Oltean, M. Solving the Hamiltonian path problem with a light-based computer. Nat. Comput. 2007, 7, 57–70. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, W.; Zhu, Z.; Yu, H. A symmetric image encryption algorithm based on a coupled logistic–bernoulli map and cellular automata diffusion strategy. Entropy 2019, 21, 504. [Google Scholar] [CrossRef] [Green Version]
  23. Saito, A.; Yamaguchi, A. Pseudorandom number generator based on the Bernoulli map on cubic algebraic integers. Chaos Interdiscip. J. Nonlinear Sci. 2017, 28, 1054–1500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Dong, L.; Yong, Z.; Ji, L.; Han, X. Study on the Pass Rate of NIST SP800-22 Statistical Test Suite. In Proceedings of the 2014 Tenth International Conference on Computational Intelligence and Security (CIS), Kunming, China, 15–16 November 2014; pp. 402–404. [Google Scholar]
  25. The USC-SIPI Image Database. Available online: http://sipi.usc.edu/database/database.php (accessed on 2 December 2019).
  26. Alvarez, G.; Li, S. Some basic cryptographic requirements for chaos-based cryptosystems. Int. J. Bifurc. Chaos 2006, 16, 2129–2151. [Google Scholar] [CrossRef] [Green Version]
  27. Castro, J.C.H.; Sierra, J.M.; Seznec, A.; Izquierdo, A.; Ribagorda, A. The strict avalanche criterion randomness test. Math. Comput. Simul. 2005, 68, 1–7. [Google Scholar] [CrossRef]
  28. Zhang, Y. The unified image encryption algorithm based on chaos and cubic S-Box. Inf. Sci. 2018, 450, 361–377. [Google Scholar] [CrossRef]
  29. Wu, Y.; Noonan, J.P.; Agaian, S. NPCR and UACI randomness tests for image encryption. Cyber J. Multidiscip. J. Sci. Technol. J. Sel. Areas Telecommun. 2011, 7714, 31–38. [Google Scholar]
  30. Ji, X.; Bai, S.; Guo, Y.; Guo, H. A new security solution to JPEG using hyper-chaotic system and modified zigzag scan coding. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 321–333. [Google Scholar] [CrossRef]
  31. Maniccam, S.S.; Bourbakis, N.G. Image and video encryption using SCAN patterns. Pattern Recognit. 2004, 37, 725–737. [Google Scholar] [CrossRef]
  32. Ramasamy, P.; Ranganathan, V.; Kadry, S.; Damaševičius, R.; Blažauskas, T. An image encryption scheme based on block scrambling, modified zigzag transformation and key generation using enhanced logistic—Tent map. Entropy 2019, 21, 656. [Google Scholar] [CrossRef] [Green Version]
  33. Richter, T. Lossless coding extensions for JPEG. In Proceedings of the Data Compression Conference, Snowbird, UT, USA, 7–9 April 2015; pp. 143–152. [Google Scholar]
  34. Chai, X.; Zheng, X.; Gan, Z.; Han, D.; Chen, Y. An image encryption algorithm based on chaotic system and compressive sensing. Signal Process. 2018, 148, 124–144. [Google Scholar] [CrossRef]
Figure 1. (a) The problem of Konigsberg bridges and (b) Hamiltonian path.
Figure 1. (a) The problem of Konigsberg bridges and (b) Hamiltonian path.
Entropy 22 00073 g001
Figure 2. Complete graphs. (a) Three nodes, (b) Four nodes, (c) Five nodes.
Figure 2. Complete graphs. (a) Three nodes, (b) Four nodes, (c) Five nodes.
Entropy 22 00073 g002
Figure 3. Process of generating a Hamiltonian path in an image of size 3 × 3. The blue pixels are chosen, and the grey pixels have been added into the path.
Figure 3. Process of generating a Hamiltonian path in an image of size 3 × 3. The blue pixels are chosen, and the grey pixels have been added into the path.
Entropy 22 00073 g003
Figure 4. Generated Hamiltonian path. The red pixel is the beginning of the path and the blue pixel is the rear of the path.
Figure 4. Generated Hamiltonian path. The red pixel is the beginning of the path and the blue pixel is the rear of the path.
Entropy 22 00073 g004
Figure 5. Bit planes of Lena, (ah) are from the 8th bit plane to the 1st bit plane, respectively.
Figure 5. Bit planes of Lena, (ah) are from the 8th bit plane to the 1st bit plane, respectively.
Entropy 22 00073 g005
Figure 6. Modified expand–shrink strategy.
Figure 6. Modified expand–shrink strategy.
Entropy 22 00073 g006
Figure 7. Hamiltonian path across bit planes.
Figure 7. Hamiltonian path across bit planes.
Entropy 22 00073 g007
Figure 8. Bernoulli map.
Figure 8. Bernoulli map.
Entropy 22 00073 g008
Figure 9. Bifurcation diagrams of ABM. (a) β = 3. The value of α is increased by 0.1, ranging from 2.1 to 202.1. (b) α = 3. The value of β is increased by 0.1, ranging from 2.1 to 202.1.
Figure 9. Bifurcation diagrams of ABM. (a) β = 3. The value of α is increased by 0.1, ranging from 2.1 to 202.1. (b) α = 3. The value of β is increased by 0.1, ranging from 2.1 to 202.1.
Entropy 22 00073 g009
Figure 10. Encryption progress.
Figure 10. Encryption progress.
Entropy 22 00073 g010
Figure 11. Decryption progress.
Figure 11. Decryption progress.
Entropy 22 00073 g011
Figure 12. Encryption process.
Figure 12. Encryption process.
Entropy 22 00073 g012
Figure 13. Histograms. (a) Plain image boat; (b) plain image male; (c) plain image clock; (d) histogram of plaintext boat; (e) histogram of plaintext male; (f) histogram of plaintext clock; (g) cipher image boat; (h) cipher image male; (i) cipher image clock; (j) histogram of cyphertext boat; (k) histogram of cyphertext male; (l) histogram of cyphertext clock.
Figure 13. Histograms. (a) Plain image boat; (b) plain image male; (c) plain image clock; (d) histogram of plaintext boat; (e) histogram of plaintext male; (f) histogram of plaintext clock; (g) cipher image boat; (h) cipher image male; (i) cipher image clock; (j) histogram of cyphertext boat; (k) histogram of cyphertext male; (l) histogram of cyphertext clock.
Entropy 22 00073 g013aEntropy 22 00073 g013b
Figure 14. One example. (a) An image in which all adjacent pixels of minor diagonal direction are equal; (b) its scatter plots.
Figure 14. One example. (a) An image in which all adjacent pixels of minor diagonal direction are equal; (b) its scatter plots.
Entropy 22 00073 g014
Figure 15. Another example. (a) An image in which all adjacent pixels of principal diagonal direction are equal; (b) its scatter plots.
Figure 15. Another example. (a) An image in which all adjacent pixels of principal diagonal direction are equal; (b) its scatter plots.
Entropy 22 00073 g015
Table 1. Randomness test using NIST SP800-22 test suite.
Table 1. Randomness test using NIST SP800-22 test suite.
Statistical TestsP-valuePass Rate (%)
Frequency0.798139100.00
Block frequency0.10879199.33
Cumulative Sums *0.28280499.83
Runs0.58865299.67
Longest run0.24507299.33
Rank0.319084100.00
FFT0.28030699.00
Non overlapping template *0.46813998.95
Overlapping template0.42505998.00
Universal0.44967299.33
Approximate entropy0.56122799.67
Random excursions *0.53300598.95
Random excursions variant *0.41954299.27
Serial *0.46463298.83
Linear complexity0.91574599.33
∗ Average value of multiple tests.
Table 2. Key sensitivity (Δ = 0.00000000000001).
Table 2. Key sensitivity (Δ = 0.00000000000001).
Boat
(512 × 512)
Couple
(512 × 512)
Tank
(512 × 512)
Male
(1024 × 1024)
Clock
(256 × 256)
Key sensitivity in encryption process(x0 + Δ, α, β)0.4993210.4999970.5000570.4999040.501156
(x0, α + Δ, β)0.4999230.5001550.4997650.4999370.500164
(x0, α, β + Δ)0.499940.5000760.5004990.5000780.500856
Key sensitivity in decryption process(x0 + Δ, α, β)0.5002890.4997410.5000820.4998580.500328
(x0, α + Δ, β)0.5001540.4994150.50010.499920.501308
(x0, α, β + Δ)0.4997470.5003370.4997420.499860.499378
Table 3. Variance of histograms.
Table 3. Variance of histograms.
Plain ImageCipher Image
Chemical plant (256 × 256)50,326.4248.469
Clock (256 × 256)282,062248.328
Moon surface (256 × 256)135,688248.094
Boat (512 × 512)1,535,8801137.66
Couple (512 × 512)1,195,4601002.11
Lena (512 × 512)632,254986.281
Tank (512 × 512)8,103,6001043.73
Airplane (1024 × 1024)115,199,0003783.7
Airport (1024 × 1024)31,596,4003832.03
Male (1024 × 1024)11,349,4004412.8
Table 4. Information entropy.
Table 4. Information entropy.
Plain ImageProposed Scheme[2][4][28]
Chemical plant(256 × 256)7.342437.997257.997167.996927.99683
Clock(256 × 256)6.705677.997277.997267.996927.99705
Moon surface(256 × 256)6.709317.997257.997387.99747.9972
Boat(512 × 512)7.191377.999227.999347.99947.99921
Couple(512 × 512)7.201017.999317.999347.999317.99936
Lena(512 × 512)7.445517.999327.999297.999347.99932
Tank(512 × 512)5.495747.999287.999347.999237.99934
Airplane(1024 × 1024)5.641457.999847.999847.999837.99981
Airport(1024 × 1024)6.830337.999847.999837.999817.99983
Male(1024 × 1024)7.523747.999817.999787.999817.99981
Table 5. Results of NPCR and UACI.
Table 5. Results of NPCR and UACI.
Index of Modified PixelNPCR
(1 Round)
UACI
(1 Round)
NPCR
(2 Rounds)
UACI
(2 Rounds)
00.9969830.33490.9959410.335169
2550.997330.3352240.9960630.3338
5110.9966160.3347270.9961430.333902
65,1510.998970.33550.9960780.333911
65,4070.9983330.3356410.9962540.334896
130,5600.9963650.3337190.9958610.334763
130,8160.999950.3356140.9961470.334734
131,0710.9999850.3358760.9962160.335797
196,0960.9999430.3361460.9958040.333875
196,3520.9990310.3354290.9962690.334645
261,6320.99810.3352040.9961810.333829
261,8880.9992490.3355510.9960370.334519
262,1430.9976080.3353610.9959980.334605
Theoretical value0.9960940.3346350.9960940.334635
Table 6. Correlation coefficients.
Table 6. Correlation coefficients.
BoatMaleClockFigure 14aFigure 15a
Plain imageHorizontal0.9365020.9780160.956658−0.0104305−0.0347679
Vertical0.9701650.9817110.973594−0.0262352−0.0253942
Principal diagonal0.9221030.9656810.9409880.03090251
Minor diagonal0.9242850.9677240.93422510.00261237
Proposed
scheme
Horizontal−0.00790818−0.00530914−0.006273260.01062150.00724055
Vertical−0.0032019−0.00593338−0.007879230.00002616870.00380631
Principal diagonal−0.007532230.0180877−0.005196520.000707789−0.0161754
Minor diagonal0.0012620.00994515−0.00729377−0.00609870.00262413
[2]Horizontal0.02727320.004440880.01055370.00498476−0.00189389
Vertical−0.0321433−0.000856047−0.00733777−0.01261750.0129397
Principal diagonal−0.00603878−0.00964336−0.01381180.0118652−0.0067004
Minor diagonal−0.00132560.0046903−0.005019110.00299615−0.0185178
[4]Horizontal0.00361182−0.00595886−0.002368480.000340202−0.0171794
Vertical0.00145023−0.0103426−0.004370460.005203040.00879099
Principal diagonal0.003954350.00305054−0.000705693−0.0120762−0.00874923
Minor diagonal−0.0001653270.002324920.000369637−0.00743531−0.00299425
[28]Horizontal0.008994910.007547750.0004110310.0057195−.00967912
Vertical−0.00416340.000629605−0.00538419−0.002668450.0105734
Principal diagonal−0.004636510.000007108760.011115−0.0034216−0.00922371
Minor diagonal0.01277110.006773950.006712560.01211950.00781511
Table 7. Efficiency of algorithms.
Table 7. Efficiency of algorithms.
Proposed Tcheme[2][4][28]
Chemical plant(256 × 256)0.008 s0.02 s0.004 s0.03 s
Clock(256 × 256)0.008 s0.019 s0.003 s0.028 s
Moon surface(256 × 256)0.008 s0.019 s0.003 s0.029 s
Boat(512 × 512)0.029 s0.093 s0.011 s0.236 s
Couple(512 × 512)0.03 s0.094 s0.012 s0.234 s
Lena(512 ×512) 0.028 s0.09 1s0.009 s0.239 s
Tank(512 × 512)0.029 s0.091 s0.011 s0.243 s
Airplane(1024 × 1024)0.129 s0.441 s0.032 s2.392 s
Airport(1024 × 1024)0.121 s0.44 7s0.033 s2.388 s
Male(1024 × 1024)0.116 s0.434 s0.033 s2.398 s
Average throughput66.062 Mbps21.884 Mbps194.571 Mbps9.542 Mbps

Share and Cite

MDPI and ACS Style

Zhang, W.; Wang, S.; Han, W.; Yu, H.; Zhu, Z. An Image Encryption Algorithm Based on Random Hamiltonian Path. Entropy 2020, 22, 73. https://doi.org/10.3390/e22010073

AMA Style

Zhang W, Wang S, Han W, Yu H, Zhu Z. An Image Encryption Algorithm Based on Random Hamiltonian Path. Entropy. 2020; 22(1):73. https://doi.org/10.3390/e22010073

Chicago/Turabian Style

Zhang, Wei, Shuwen Wang, Weijie Han, Hai Yu, and Zhiliang Zhu. 2020. "An Image Encryption Algorithm Based on Random Hamiltonian Path" Entropy 22, no. 1: 73. https://doi.org/10.3390/e22010073

APA Style

Zhang, W., Wang, S., Han, W., Yu, H., & Zhu, Z. (2020). An Image Encryption Algorithm Based on Random Hamiltonian Path. Entropy, 22(1), 73. https://doi.org/10.3390/e22010073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop