Next Article in Journal
Multifractal and Chaotic Properties of Solar Wind at MHD and Kinetic Domains: An Empirical Mode Decomposition Approach
Next Article in Special Issue
A Chaotic Electromagnetic Field Optimization Algorithm Based on Fuzzy Entropy for Multilevel Thresholding Color Image Segmentation
Previous Article in Journal / Special Issue
Kapur’s Entropy for Color Image Segmentation Based on a Hybrid Whale Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Image Encryption Based on Pixel-Level Diffusion with Dynamic Filtering and DNA-Level Permutation with 3D Latin Cubes

1
School of Economic Information Engineering, Southwestern University of Finance and Economics, Chengdu 611130, China
2
College of Computer Science, Sichuan University, Chengdu 610064, China
3
College of Electronics and Information Engineering, Sichuan University, Chengdu 610064, China
*
Authors to whom correspondence should be addressed.
Entropy 2019, 21(3), 319; https://doi.org/10.3390/e21030319
Submission received: 12 March 2019 / Revised: 19 March 2019 / Accepted: 21 March 2019 / Published: 24 March 2019
(This article belongs to the Special Issue Entropy in Image Analysis)

Abstract

:
Image encryption is one of the essential tasks in image security. In this paper, we propose a novel approach that integrates a hyperchaotic system, pixel-level Dynamic Filtering, DNA computing, and operations on 3D Latin Cubes, namely DFDLC, for image encryption. Specifically, the approach consists of five stages: (1) a newly proposed 5D hyperchaotic system with two positive Lyapunov exponents is applied to generate a pseudorandom sequence; (2) for each pixel in an image, a filtering operation with different templates called dynamic filtering is conducted to diffuse the image; (3) DNA encoding is applied to the diffused image and then the DNA-level image is transformed into several 3D DNA-level cubes; (4) Latin cube is operated on each DNA-level cube; and (5) all the DNA cubes are integrated and decoded to a 2D cipher image. Extensive experiments are conducted on public testing images, and the results show that the proposed DFDLC can achieve state-of-the-art results in terms of several evaluation criteria.

1. Introduction

As one of the most important information carriers, hundreds of millions of images are generated, stored, and transmitted every day. How to ensure image security has become a very hot topic of research in recent years. Image encryption is one of the most important image security methods. Encryption algorithms can be roughly classified into two categories: symmetric key (private key) and asymmetric key (public key) algorithms. The former uses the same key for both encryption and decryption, while the latter uses a key for encryption and another key for decryption. Typical private key algorithms include data encryption standard (DES), international data encryption algorithm (IDEA), advanced encryption standard (AES) and so on. Rivest-Shamir-Adleman (RSA) and Elliptic-curve cryptography (ECC) are among the very popular public key algorithms. The symmetric key algorithms are fast, efficient, but difficult to manage keys, while the asymmetric encryption algorithms are slow but have higher security [1,2]. Due to the inherent characteristics of images such as strong correlation, high redundancy and bulky data capacity, the above mentioned encryption algorithms are usually not suitable for direct applications in images. To address this issue, a variety of image encryption algorithms have been proposed in recent years [3,4,5,6].
There are many kinds of operations for the purpose of image encryption, such as shuffling, permutation, rotation, substitution, confusion, diffusion, transposition, and so on [7]. Among the operations, diffusion and permutation are very popular ones because they can achieve good results and are easy to implement. The diffusion is to change the values of the pixels in images, while the permutation aims at changing the positions of the pixels. Some practical image encryption algorithms are capable of handling diffusion and permutation jointly. Due to the characteristics of ergodicity, pseudorandomness, unpredictability, and extreme sensitivity to initial values and parameters, chaos-based image encryption has become increasingly popular in recent years. The main idea of chaos-based image encryption is to conduct diffusion and/or permutation according to the pseudorandom sequences generated from chaotic systems [8,9,10,11,12,13]. Very recently, Flores-Vergara et al. have implemented a chaotic cryptosystem on embedded systems with multiprocessors. The NIST statistical test and the security analysis have confirmed the proposed cryptosystem is very secure and robust for image encryption [14]. Wang et al. used a spatiotemporal chaotic system to generate a pseudorandom sequence, and then used the sequence to conduct permutation and diffusion simultaneously [15]. Pareek et al. employed two chaotic logistic maps and eight different types of operations to encrypt the pixels of images, and the experimental results demonstrated the proposed scheme was real-time, efficient and secure [16]. Hua et al. put forward a new 2D Logistic-Sine-coupling map that has more complex behavior, better ergodicity, and larger chaotic range than some other 2D chaotic maps, for image encryption scheme. The experiments showed that the proposed scheme had better security performance than several state-of-the-art encryption approaches [17]. Sahari and Boukemara proposed a novel 3D chaotic map by integrating the piecewise and logistic maps for color image encryption, the experimental results showed the efficiency and safety of the proposed scheme [18]. Zhou et al. proposed a novel image encryption scheme by combining quantum 3D Arnold transform and quantum XOR operations with scaled 3D Zhongtang chaotic system [19]. Low-dimensional chaotic systems have the advantages of simple forms, only a few parameters, and easy implementation. However, such properties may make it easy to estimate the orbits and the initial parameters of the low-dimensional chaotic systems and hence the security of encryption is limited.
In a dynamical system, the Lyapunov exponent (LE) is used to measure the rate of separation of infinitesimally close trajectories [20]. If a chaotic system has at least two positive LEs, the system is said to be hyperchaotic. The image encryption algorithms with hyperchaotic systems have been demonstrated more secure [2,6,21,22,23,24,25,26]. Chai et al. used a 4D memristive hyperchaotic system to encrypt 4 compound bit planes recombined from the 24-bit planes of components R, G, and B [27]. Li et al. proposed a quantum image compression-encryption approach with quantum cosine transform and a 5D hyperchaotic system, and the experiments demonstrated that the proposed compression-encryption approach outperformed some classical image encryption approaches [28]. Zhou et al. used a 5D hyperchaotic system for quantum color image encryption. Some researchers also applied 6D or 7D hyperchaotic systems to generate hyperchaotic sequences for image encryption [6,29].
Like other tasks in signal processing, image encryption can also be conducted in both spatial or transform domain [30,31,32,33,34,35]. The encryption in spatial domain is very direct, which changes the values and/or the positions of pixels. To improve the efficiency or the effectiveness of image encryption, sometime the algorithms can be conducted on higher-level data (blocks of pixels) or lower-level data (DNA-level data and bit-level data) [36,37,38]. Generally speaking, for the same processing power, the lower the data level, the more pixels will be involved in encryption. Therefore, the encryption processing lower-level data usually achieves better performance of encryption [6]. In the field of image encryption, the introduction of transform domain is for the purpose of compressing images. Typical transform methods include discrete cosine transform (DCT) [39,40,41,42,43], Fourier transform [44,45,46], and wavelet transform [47,48,49,50]. With these transform methods, the image encryption can focus on the high-energy parts of the images only and discard some low-energy (zero coefficients) parts. Then the image can be recovered by decryption and corresponding reverse transform.
Some recent progress has improved the performance of image encryption. Regarding diffusion, Hua and Zhou introduced filtering, a very popular technique in image processing, into image encryption. The authors make the filtering reversible by setting the right-bottom point of the filter to “1”, and they proposed an image encryption algorithm using block-based scrambling and image filtering (BSIF) with a fixed filter for all pixels [51]. Very recently, Hua et al. have extended image encryption with Josephus scrambling and filtering diffusion, where the filter is a 2 × 2 square with fixed values [52]. Li et al. used a 1 × 3 or 3 × 1 filter with dynamically variable values decided by a 7D hyperchaotic system for filtering (so-called dynamic filtering), and bit cuboid operations, namely DFBC, for image encryption, and the experiments demonstrated the DFBC could achieve state-of-the-art results [6]. As far as permutation is concerned, in theory, any reversible position transform can be used for image encryption. Latin squares are such popular transforms which help to achieve good results of permutation [53,54,55]. Xu et al. extended the use of Latin squares in image encryption, and they treated the pixel-level image as a 3D bit matrix and then conducted operations of Latin cubes on the 3D matrix, and the experimental results showed that the proposed image encryption achieves both a desirable level of security and high efficiency [56].
Motivated by the diffusion with filtering and the permutation with Latin cubes, in this paper, we propose a novel approach that integratings a hyperchaotic system, Dynamic Filtering, DNA computing, and Latin Cubes, termed as DFDLC, for image encryption. Specifically, the DFDLC consists of five stages: (1) A 5D hyperchaotic system with 2 positive LEs is applied to generate the chaotic sequences for subsequent diffusion and permutation. (2) Filters with variable values are generated from the chaotic sequences, and filtering is conducted on each pixel of the image with a different filter. That is to say, the value of each pixel is changed by a different filter. This is called pixel-level diffusion with dynamic filtering. (3) The 2D pixel plane is transformed into several DNA cubes via DNA encoding rules determined by the chaotic sequence. (4) For each DNA cube, we generate a Latin cube with the same size and then change the position of each element in the DNA cube via the Latin cube. This operation is called DNA-level permutation with 3D Latin cubes. (5) All the DNA cubes are integrated and decoded to a 2D pixel image. The main contributions of this paper are three-aspect: (1) We propose a novel image encryption using a newly found 5D hyperchaotic system. (2) Pixel-level dynamic filtering and DNA-level permutation with Latin cubes are used to improve the performance of image encryption. (3) Extensive experiments on several public images show that the DFDLC is promising for image encryption.
The remainder of this paper is structured as the following. A brief description of a 5D hyperchaotic system with two positive LEs, filtering, DNA computing and Latin square is given in Section 2. In Section 3, a novel image encryption algorithm with dynamic filtering and Latin cube transformation, namely DFDLC, is proposed in detail. Experimental results are reported and analyzed in Section 4. Finally, the paper is concluded in Section 5.

2. Preliminaries

2.1. Hyperchaotic Systems

As one of the most popular chaotic systems, the Lorenz chaotic system and its extensions are very popular in image encryption. Most recently, Wang et al. have found a new 5D autonomous hyperchaotic system with 2 positive LEs by adding feedback controllers to the Lorenz system, formulated as Equation (1) [57]:
x 1 ˙ = x 2 x 2 ˙ = x 2 + a x 1 + b x 1 3 + c x 1 x 5 x 3 ˙ = x 4 x 4 ˙ = x 4 + d x 3 + e x 3 3 + f x 3 x 5 x 5 ˙ = g x 5 + h x 1 2 + i x 3 2
where x j ( j = 1 , 2 , , 5 ) are state variables, and ( a , b , c , d , e , f , g , h , i ) are constant parameters. There are several numerical methods to solve this system, such as Forward Euler (FE), 4th order Runge-Kutta (RK) and newly proposed trigonometric polynomials [58]. In this paper, we use the 4th order RK method with a step size of h = 0.001 to solve the hyperchaotic system. When the parameters ( a , b , c , d , e , f , g , h , i ) = ( 4 , 1 , 1 , 2 , 1 , 2 , 0.0 , 6 , 1 ) and initial values ( x 1 0 , x 2 0 , x 3 0 , x 4 0 , x 5 0 ) = ( 1.618 , 3.14 , 2.718 , 4.6692 , 0.618 ) × 10 2 , the attractors of the 5D hyperchaotic system are shown in Figure 1.

2.2. Filtering

Filtering, also termed as convolution, is a very popular operation in the field of image processing, which can be applied to denoising, smoothing, and sharpening images by changing the values of pixels. Typically, the operation of filtering is to do convolution between a mask, also known as a kernel/filter/window, and an image. The values of pixels in an image are changed and hence it seems that filtering can be used for diffusion directly. However, since traditional filtering cannot be reversible, the cipher image with such diffusion cannot be recovered. To cope with this issue, Hua and Zhou set the right-bottom point of the filter to “1” and then align this point to the processed pixel in the image for convolution, and they proposed a novel image encryption algorithm with block-based scrambling and such image filtering (BSIF) [51]. However, the BSIF used a fixed filter for all pixel when doing convolution, limiting the encryption performance. An ideal scheme should use a variable/dynamic filter for convolution with each pixel.

2.3. DNA Computing

DNA computing, invented by Leonard Adleman, is a type of parallel computing technique that the information is expressed by four nucleic acids, i.e., adenine (A), cytosine (C), guanine (G), and thymine (T) [59]. The key factors of DNA for encryption are encoding and decoding rules, and algebraic operations for DNA sequences. Like 0 and 1 are complementary pairs in binary, 00 (0) and 11 (3), and 01 (1) and 10 (2) are also complementary pairs in DNA computing. Although there are 4 ! = 24 combinations in total for DNA encoding, there are only 8 kinds of DNA bases are capable of meeting the DNA complementary rules, as listed in Table 1. With the encoding rule, an 8-bit pixel in grayscale image can be expressed by 4 letters. For example, following Rule 5 and Rule 8 in Table 1, the decimal gray-level 156 (’10011100’ in binary) can be transformed into a 4-letter DNA sequence ’TAGC’ orand ’ATCG’, respectively. It can be seen that for a fixed binary sequence, different rules lead to totally different DNA sequences.
In image encryption, several algebraic operations, such as addition (++), subtraction (--) and exclusive OR (XOR, ), as listed in Table 2, Table 3 and Table 4, can be used to change the values of nucleic acids [2].

2.4. Latin Square

A Latin square of order N is an N × N matrix which includes a set S with N different symbol elements, and each symbol shows only once in each row and each column [53]. For instance, L is a Latin square of order N, i and j represent the row and column index of an element in L respectively, and S k is the k - th element in set S. We can draw a formula as follows:
f ( i , j , k ) = 1 , L ( i , j ) = S k 0 , o t h e r w i s e
Given S = { 0 , 1 , , N 1 } , Figure 2 shows an example of Latin square of order 4.

3. The Proposed Image Encryption Approach

3.1. Hyperchaotic Sequence Generation

In this paper, we used the 5D hyperchaotic system described in Section 2.1 to generate the hyperchaotic sequence for encryption. Specifically, the generating process has three steps:
Step 1:
The sequences generated by the first N 0 iterations are discarded to eliminate the adverse effects.
Step 2:
The 5D hyperchaotic system continues to iterate to generate sequences long enough for image encryption. In the j - th iteration, we can obtain five state values denoted as s j = { x 1 j , x 2 j , , x 5 j } .
Step 3:
When the iteration completes, a hyperchaotic sequence S can be obtained by contacting all the s j ( j = 1 , 2 , , N ) as
S = { s 1 , s 2 , , s N } = { x 1 1 , x 2 1 , , x 5 1 , , x 1 N , x 2 N , , x 5 N } = { s 1 , s 2 , s 3 , , s 5 N 2 , s 5 N 1 , s 5 N } .
The real value sequence S is further mapped to an integral sequence as Equation (4):
k i = m o d ( m o d ( ( k i k i ) × 10 15 ) , 10 8 ) , 256 ) ,
where mod, · and · denote the operations of modulo, absolute value, and flooring, respectively [2,6].

3.2. Dynamic Filtering

The modified filtering can be applied to image encryption, according to the very recent work BSIF by Hua and Zhou [51]. However, the BSIF does convolution on each pixel in an image with a fixed kernel generated from a random sequence. Li et al. used a 1 × 3 or 3 × 1 variable kernel to convolute each pixel in an image, that is to say, the kernels associated with each pixel for convolution are different, so-called dynamic filtering [6]. The experimental results demonstrated the effectiveness of dynamic filtering. A reasonable assumption is that a dynamic kernel with larger size (e.g., 3 × 3 or 5 × 5 ) will lead to better encryption. An example of dynamic filtering with two 3 × 3 filters is shown in Figure 3, where the 3 × 3 red kernel and the 3 × 3 blue kernel are conducted on the pixels of 34 and 178 in the plain image, and the results of dynamic filtering will be 140 and 214 in the cipher image, respectively. We can see that with dynamic filtering, the values of pixels in the plain image are changed, and this procedure can be reversible [51]. Therefore, we can use this operation for diffusion.

3.3. Image to Cubes

Since 3D Latin cube transformation can be conducted on cubes only, the image for encryption must be reshaped to one or several cubes. The pseudocode of such transformation algorithm (I2C) can be described as below.
Step 1:
Given an image with size h × w × d , where h , w , and d represent the height, width, and depth, respectively, calculate the number of the pixels N = h × w × d .
Step 2:
Let L = N 3 , if L is an integer, jump to Step 3, else jump to Step 4.
Step 3:
Get a cube with size L × L × L , return.
Step 4:
Define K = 2 n , n N , find the biggest K that meets K L ; then we get a cube with size K × K × K .
Step 5:
Update N = N K 3 , if N = 0 , return; else jump to Step 2.
For instance, a DNA-level image with size 512 × 512 × 4 can be transformed into 4 cubes with size 64 × 64 × 64 , while a DNA-level image with size 256 × 256 × 4 can be transformed into 8 32 × 32 × 32 cubes. Unlike the previous work that can only encrypts images of specified sizes [56], the proposed DFDLC can handle images of any sizes with such transformation.
Accordingly, one or several cubes can be merged into a plain image with the reverse procedure of the I2C.

3.4. 3D Latin Cube

Latin cube is a generalized version of the Latin square. A Latin cube of order N is an N × N × N cube which includes a set S with N different symbol elements, and each symbol occurs only once in each row, each column, and each file [56]. Given a chaotic sequence x = { x 0 , x 1 , , x q n 1 } (q is a prime and q n is the order of the Latin cubes to generate), we can sort the sequence by ascending to get an index sequence y = { y 0 , y 1 , , y q n 1 } and then construct a finite field F q n on y via redefining “+” and “×”. With three distinct nonzero elements p 1 , p 2 and p 3 in F q n , the element of L t ( i , j , s ) can be obtained by Equation (5):
L t ( i , j , s ) = y s + p t × y j + p t 2 × y i ,
where t = { 1 , 2 , 3 } is the index of the Latin cube, and “+” and “×” are the addition and multiplication in F q n , respectively [56]. Figure 4 shows three Latin cubes of order 3 on the set S = { 0 , 1 , 2 } , named as L 1 , L 2 , and L 3 . When we superimpose the same position of three Latin cubes on the set S, if each combination occurs only once, we can say these three Latin cubes are orthogonal. For example, when we combine the three Latin cubes L 1 , L 2 and L 3 , each of the 27 combinations 000 , 001 , 002 , , 222 occurs only once, so they are orthogonal. By combining L 1 , L 2 and L 3 , we can get a new cube K shown in Figure 5. Then a spatial permutation is obtained: ( 0 , 0 , 0 ) ( 0 , 0 , 0 ) , ( 0 , 1 , 0 ) ( 1 , 1 , 2 ) , ( 0 , 2 , 0 ) ( 2 , 2 , 1 ) , , ( 2 , 2 , 2 ) ( 0 , 2 , 1 ) , i.e., the element in the left position is transferred to the right position. More generally, K s ( i , j ) = ( L 1 ( i , j , s ) , L 2 ( i , j , s ) , L 3 ( i , j , s ) ) , where s is the index of K (or L), and i and j are the indices of the row and the column, respectively.

3.5. DFDLC: The Proposed Image Encryption Approach with Dynamic Filtering and Latin Cubes

The DFDLC is conducted on pixel-level diffusion and DNA-level permutation. Specifically, regarding pixel-level diffusion, we mainly apply dynamic filtering on each pixel in a plain 2D image. We also used the ciphertext diffusion in crisscross pattern (CDCP) to improve the diffusion results [60]. For DNA-level permutation, we mainly use Latin cube to change the position of each nucleic acid. In addition, a kind of global DNA permutation similar to the global bit permutation is adopted for DNA permutation [6]. The proposed DFDLC is illustrated in Figure 6. With the hyperchaotic sequence generated by the 5D chaotic system, the main steps of the DFDLC are as the following: hyperchaotic sequence generation, pixel-level diffusion, pixel-to-DNA transformation, DNA permutation, and DNA-to-pixel transformation.
The details of the DFDLC are described as follows:
Step 1:
Given the keys, generate a hyperchaotic sequence with Equations (1), (3) and (4).
Step 2:
Conduct CDCP with pixels of the image. This operation expands a little change in one pixel of the plain image to very large changes in a variety of pixels of the cipher image.
Step 3:
Dynamic filtering on the image. For each pixel, firstly, generate a 3 × 3 kernel with the hyperchaotic sequence and set the right-bottom grid to 1. Secondly, do convolution with the kernel and corresponding sub-region of the image associated with the pixel. Thirdly, use the result of the convolution as the new value of the pixel in the cipher image.
Step 4:
Transform the pixel image to a DNA image. For each pixel, use an encoding rule decided by the hyperchaotic sequence to encode one pixel into a string with 4 nucleic acids. The DNA encoding rule (Rule N) can be formulated as: N = 1 + m o d ( x , 8 ) , where x is a corresponding value in the hyperchaotic sequence regarding the pixel.
Step 5:
Transform the DNA image into one or several cubes using I2C.
Step 6:
Conduct DNA-level Latin cube permutation. For each DNA-level cube, generate a Latin cube and then change the position of each nucleic acid according to the Latin cube. In addition, the DNA XOR operation is conducted on the DNA-level cube with a generated DNA cube from the hyperchaotic sequence.
Step 7:
Integrate all the DNA-level cubes into a DNA image.
Step 8:
Conduct global DNA permutation as described in [2].
Step 9:
Decode the DNA image into a pixel image. For each nucleic acid, the DNA encoding rule is decided as the encoding rule in Step 6. The pixel image is the cipher image.
The proposed DFDLC consists of five stages: hyperchaotic sequence generation (Step 1), pixel-level diffusion (Step 2-3), a transformation from a plane image to cubes (Step 4-5), DNA-level Latin cube permutation (Step 6-8) and a transformation from cubes to a plane image (Step 9). The keys of the DFDLC are Step 3 and Step 6, i.e., pixel-level diffusion with dynamic filtering and DNA-level permutation with 3D Latin cubes, respectively. Although the main objective of Latin cubes in the DFDLC is for permuting the DNA, it also results in diffusion because the change of the position of DNA can change the corresponding value of the pixel naturally [6].
The cipher image can be easily decrypted by the inverse steps as listed above.

4. Experimental Results

4.1. Experimental Settings

To validate the performance of the proposed DFDLC, we compare it with some state-of-the-art image encryption schemes, such as the image encryption with a fractional-order hyperchaotic system and DNA computing (FOHCDNA) [2], the hyperchaotic and DNA sequence-based method (HC-DNA) [61], CDCP [60], BSIF [51] and DFBC [6]. We set the parameters for the DFDLC as following. For the 5D hyperchaotic system, we set ( x 1 0 , x 2 0 , x 3 0 , x 4 0 , x 5 0 ) = ( 1.618 , 3.14 , 2.718 , 4.6692 , 0.618 ) × 10 2 and 1000 as the initial values and the preiterating times, respectively. For the compared methods, we use the parameters as set by the corresponding references.
We used ten publicly accessed images for validating the proposed DFDLC, and the details of the images are listed in Table 5.
The experiments were conducted using MATLAB 2016b (MathWorks, Natick, MA, USA) on a 64-bit Windows 7 Ultimate (Microsoft, Redmond, WA, USA) with 32 GB memory and a 3.6 GHz I7 CPU.

4.2. Security Key Analysis

A feasible image encryption algorithm should have a large enough key space and extreme sensitivity to the key to resist brute force attacks. In this subsection, we will analyze the key space and the sensitivity of the security key.

4.2.1. Key Space

The key space is the set of all possible security keys that can be used in a system of image encryption. It was reported that the size of a key space larger than 2 100 can provide enough security [62]. Basically, the 5 initial values of the 5D hyperchaotic systems, i.e., ( x 1 0 , x 2 0 , x 3 0 , x 4 0 , x 5 0 ) for Equation (1), can be constructed as the security keys. If each initial value has the same precision of 10 15 , the DFDLC has a key space with size of 10 15 5 = 10 75 2 249 , which is much larger than 2 100 . Therefore, the DFDLC can resist all types of brute-force attacks from current computers. Besides, the distinct nonzero elements in the finite filed for Latin cubes can be used as security keys to improve the key space.

4.2.2. Sensitivity to Security Key

An ideal image encryption approach should be sensitive enough to the security key, that is to say, a very little change in the security keys will lead to a completely different decrypted image.
We use two groups of slightly different keys to validate the sensitivity to the security keys of the proposed DFDLC. The first group keys are the initial values of the hyperchaotic system, i.e., g 1 = ( x 1 0 , x 2 0 , x 3 0 , x 4 0 , x 5 0 ) = ( 1.618 , 3.14 , 2.718 , 4.6692 , 0.618 ) × 10 2 , while the second groups are almost the same as the first group except x 1 0 = 0.0168 + 10 15 , i.e., g 2 = ( x 1 0 + 10 13 , x 2 0 , x 3 0 , x 4 0 , x 5 0 ) = ( 1.618 + 10 13 , 3.14 , 2.718 , 4.6692 , 0.618 ) × 10 2 . We apply g 1 and g 2 to decrypt the first five images in Table 5, and the results are shown in Figure 7. It is clear that even the security keys are changed very little such as 10 15 , the cipher images cannot be recovered correctly, demonstrating the high sensitivity to security keys of the proposed DFDLC [6].

4.3. Statistical Analysis

Statistical analysis, including histogram analysis, entropy analysis, and correlation analysis are essential for a cryptosystem. An ideal image encryption algorithm should have the ability to resist kinds of statistical attacks.

4.3.1. Histogram Analysis

Histogram describes the distribution of pixels for an image. The histogram of a natural image usually shows an irregular (unevenly distributed) shape. A good image encryption approach should change the irregular shape of a plain image as evenly distributed as possible, leading to a completely random-like cipher image. Regarding evaluating the image encryption approach with histogram, the more uniform the histogram is, the better the encryption approach is [2]. The histograms of the plain images and the corresponding cipher images are shown in Figure 8.
It can be seen that the histograms of the plain images except Bw look like mountains, including peaks and valleys. However, the histograms of their corresponding cipher images are so flat that they are very close to uniform distributions. It is worth pointing out that regarding the image Bw, it has only two values of grayscale level, i.e., 0 and 255, and its histogram looks like two needles. However, the histogram of its cipher image is still very uniform similar to histograms of other cipher images. Although the plain images are very different, the histograms of their corresponding cipher images are so uniform and so close that it looks like that each grayscale level appears about 1000 times in all cipher images. This characteristic of cipher images can be easily found in the last column in Figure 8. The experiments indicate that the proposed DFDLC can obtain very uniform histograms for different types of images and hence it can resist histogram attacks very well.

4.3.2. Information Entropy

Information entropy (IE), originally proposed by Shannon, is one of the key measures to quantify the degree of uncertainty (randomness) of a given system in information theory [63]. It can be applied to measure the randomness of an image encryption system. Given an 8-bit grayscale level that has 2 8 = 256 possible pixel values, i.e., 0 , 1 , , 255 , the IE can be formulated as Equation (6)
IE ( I ) = i = 0 255 p ( I i ) l o g 2 p ( I i ) ,
where p ( I i ) is the probability of the i - th gray value I i appears in an image I. For a cipher image, when each gray value I i appears with equal probability, i.e., 1 256 , the IE obtains the maximum 8. Therefore, an ideal image encryption approach should have an IE close to 8.
The IEs of the test images and corresponding cipher images with the DFDLC and the compared approaches are shown in Table 6. It can be seen that the testing natural images in this experiment have close IEs around 7, while the image of Bw has the lowest entropy 1, showing that the distribution of pixel values is irregular, as indicated by their histograms in Figure 8. It can be seen that the IEs of all cipher images are very close to the ideal value 8. Specifically, all encryption approaches except for HCDNA achieve very stable IEs, i.e., 7.9992 7.9994 , which are also very close to 8, indicating that these approaches are secure enough to resist entropy attacks. Although the IEs achieved by the HCDNA are slightly worse than those by the other approaches, they are still very close to the ideal value except that the IE of Bw by HCDNA is as low as 7.9158 . Among the approaches, the BSIF obtains the highest IEs with 6 out of 10 cases, followed by DFDLC, FHDNA and DFBC, which all achieve the highest IEs 4 out of 10 times. However, the HCDNA achieves the highest IE only once. The experimental results demonstrate that the DFDLC are advantageous over or comparable to other approaches in terms of IE.
As mentioned above, the IEs reflect the randomness of the grayscale values in an image. The IEs achieved by DFDLC are very close to 8, indicating that the pixel values are distributed very uniformly, as the histograms shown in the last column in Figure 8. Therefore, the results of histograms are consistent with the analysis of IEs, confirming that the proposed DFDLC has good statistical properties in terms of image encryption.

4.3.3. Correlation Analysis

Natural images usually show high correlation, that is, neighboring pixels have very close grayscale levels. When an image is permutated, the neighboring pixels will be randomly distributed in the whole image and hence the high correlation in plain image is broken. An ideal image encryption approach should decrease the correlation to zero in the cipher image. One of the popular ways to measure the correlation in images is the correlation coefficient γ defined as Equation (7) [6,64]
E ( x ) = 1 M i = 1 M x i , S ( x ) = 1 M i = 1 M ( x i E ( x ) ) 2 , c o v ( x , y ) = 1 M i = 1 M ( x i E ( x ) ) ( y i E ( y ) ) , γ = c o v ( x , y ) S ( x ) S ( y ) ,
where x and y are grayscale levels of two adjacent pixels in an image, and M denotes the number of pairs of involved pixels, and E ( x ) , S ( x ) and c o v ( x , y ) are the expectation of x, the standard deviation of x and the covariance of x and y, respectively.
To analyze the correlation, we firstly use all the pairs of adjacent pixels from each plain image and the corresponding cipher image in the horizontal direction, the vertical direction, and the diagonal direction to compute the correlation coefficients, denoted by γ h , γ v and γ d , respectively. The results are shown in Table 7. We can see that the correlation coefficients of all plain images in all directions are very high, especially the γ h of the image Bw equals to the maximum value of 1. However, all the correlation coefficients of the encrypted images decrease to close to zero, showing that the high correlation in plain images is broken. Regarding the encryption approaches, each outperforms others in several cases, indicating they are comparable in terms of reducing the correlation in images. If we consider the range of the γ achieved by the approaches, we can see that the ranges by DFDLC, FHDNA, HCDNA, CDCP, BSIF and DFBC are [−0.0023,0.0030], [−0.0049, 0.0057], [−0.0032, 0.0038], [−0.0032, 0.0028], [−0.0032, 0.0034] and [−0.0029, 0.0027], respectively. Accordingly, the interval widths of γ by the approaches are 0.0053 , 0.0106 , 0.0070 , 0.0060 , 0.0066 and 0.0056 . Among the interval widths, the DFDLC achieves the narrowest one, indicating that the DFDLC is the most stable approach in terms of γ .
Then, we randomly select 4000 pairs of horizontally adjacent pixels from each plain image and its corresponding cipher image to plot the distribution maps of the grayscale levels of the adjacent pixels, as shown in Figure 9. It can be seen that the correlation of natural images is so strong that the grayscale levels of the adjacent pixels are concentrated near the diagonal line. The figure of the plain Bw is a special case because its distribution has only two possible combinations, i.e., ( 0 , 0 ) and ( 255 , 255 ) . The strong correlation of all the plain images is thoroughly destroyed by the proposed DFDLC so that the grayscale levels of adjacent pixels are evenly distributed over the entire plane. It further demonstrates that the DFDLC has good performance regarding correlation.

4.4. Analysis of Resisting Differential Attacks

Differential attack is to study how a tiny change in a plain image can affect the corresponding cipher image. A good encryption approach should have the ability to resist differential attacks, that is to say, any small changes (even if changing a bit) in a plain image will result in a completely different cipher image. Two of the most popular indices to quantify the performance of resisting differential attacks in image encryption are the number of pixels change rate (NPCR) and the unified average changing intensity (UACI), as defined by Equations (8) and (9), respectively [65]
NPCR = 1 W × H i = 1 W j = 1 H d i j × 100 % ,
UACI = 1 255 × W × H i = 1 W j = 1 H C i j 1 C i j 2 × 100 % ,
where W and H denote the width and the height of the cipher images respectively, C 1 and C 2 are two cipher images, and d i j is defined as Equation (10)
d i j = 0 , C i j 1 = C i j 2 , 1 , C i j 1 C i j 2 .
As far as the two indices are concerned, the NPCR focuses on the variation ratio of two cipher images whose plain images are slightly changed while the UACI defines the mean intensity of the two cipher images. Wu et al. proposed a threshold and a range for NPCR and UACI respectively to evaluate if an encryption approach can pass the differential attack test for a given specified size image at a significance level α . Specifically, for a 512 × 512 8-bit grayscale image, if the NPCR score is bigger than the threshold N 0.05 * = 99.5893 % , it passes the NPCR test at α = 0.05 . In addition, if the UACI score falls into the interval U 0.05 * l , U 0.05 * u = ( 33.3730 % , 33.5541 % ) , it is said to pass the UACI test at α = 0.05 [65].
We add 1 to the value of a randomly selected pixel to compute one score of the NPCR and the UACI. The computation is repeated 10 times and then the mean, standard deviation, and times of passing the test of NPCR and UACI are reported in Table 8 and Table 9, respectively. The mean scores that pass the NPCR or the UACI tests at a significance level α = 0.05 are shown in bold. One can see that both DFDLC and BSIF can pass both tests on all images in terms of the mean scores of NPCR and UACI, while CDCP and DFDC can pass most tests. In contrast, the FHDNA and the HCDNA failed the tests with all images, although the mean scores by the FHDNA are very close to N 0.05 * and U 0.05 * l , U 0.05 * u . If we look at the times of passing the NPCR test, both the DFDLC and the BSIF can pass the test in 99 out of 10 × 10 = 100 times and they are far superior to other methods. However, regarding times of passing the UACI test, the DFDLC is slightly worse than the BSIF, but it outperforms other methods. The experimental results demonstrate that the proposed DFDLC is capable of resisting differential attacks.

4.5. Discussion

The proposed DFDLC conducts encryption on pixel-level and DNA-level, with dynamic filtering for diffusion and Latin cubes for permutation. From the above analysis, we can see that the DFDLC can resist brute force attacks, statistical attacks as well as differential attacks, and the experiments have also demonstrated that DFDLC is superior or comparable to the compared state-of-the-art image encryption methods. In addition, the proposed I2C allows the DFDLC to handle images with any sizes, making it more practical.
One limitation of the DFDLC is the running time. It takes about 0.84s and 3.15s to encrypt an image of size 256 × 256 and 512 × 512 respectively in our experimental environment. The DFDLC is time consuming because the DNA operations (DNA encoding, decoding and algebraic operation) are actually operations on strings. This can be resolved by introducing lookup tables of DNA operations. Another possible way is to use GPU to accelerate DNA operations.

5. Conclusions

Image encryption is one of the core tasks of image security. To improve image security, in this paper, a novel image encryption algorithm that uses a 5D hyperchaotic system with 2 positive LEs, pixel-level dynamic filtering, DNA computing, and 3D Latin cubes, namely DFDLC, is proposed. The novelty of the DFDLC is introducing a new type of dynamic filtering to conduct pixel-level diffusion and permutating images with DNA-level data via Latin cubes. Extensive experiments on ten public test images have indicated that the proposed DFDLC has a large key space, is very sensitive to security keys, has good statistical characteristics, and can resist types of attacks. In the future, we will extend the proposed DFDLC in several aspects. First, we will apply trigonometric polynomials to generate the hyperchaotic sequence for the DFDLC. Second, we will try a variety of shapes of the filters for dynamic filtering. Third, we may use GPU or lookup tables to speed up the encoding and decoding of DNA and corresponding arithmetic operations. Finally, we can apply the DFDLC to color image encryption.

Author Contributions

Investigation, T.L., J.W. and F.P.; Methodology, T.L. and X.L.; Software, T.L., J.S. and X.L.; Supervision, T.L.; Writing—original draft, T.L. and J.S.; Writing—review & editing, T.L. and J.W.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (Grant No. JBK1902029, No. JBK1802073 and No. JBK170505), the Ministry of Education of Humanities and Social Science Project (Grant No. 19YJAZH047), Sichuan Science and Technology Program (Grant No. 2019YFG0117) and the Scientific Research Fund of Sichuan Provincial Education Department (Grant No. 17ZB0433).

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (Grant No. JBK1902029, No. JBK1802073 and No. JBK170505), the Ministry of Education of Humanities and Social Science Project (Grant No. 19YJAZH047), Sichuan Science and Technology Program (Grant No. 2019YFG0117) and the Scientific Research Fund of Sichuan Provincial Education Department (Grant No. 17ZB0433).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmad, J.; Hwang, S.O.; Ali, A. An experimental comparison of chaotic and non-chaotic image encryption schemes. Wirel. Pers. Commun. 2015, 84, 901–918. [Google Scholar] [CrossRef]
  2. Li, T.; Yang, M.; Wu, J.; Jing, X. A novel image encryption algorithm based on a fractional-order hyperchaotic system and DNA computing. Complexity 2017, 2017, 9010251. [Google Scholar] [CrossRef]
  3. Abd El-Latif, A.A.; Abd-El-Atty, B.; Talha, M. Robust encryption of quantum medical images. IEEE Access 2018, 6, 1073–1081. [Google Scholar] [CrossRef]
  4. Guillén-Fernández, O.; Meléndez-Cano, A.; Tlelo-Cuautle, E.; Núñez-Pérez, J.C.; de Jesus Rangel-Magdaleno, J. On the synchronization techniques of chaotic oscillators and their FPGA-based implementation for secure image transmission. PLoS ONE 2019, 14, e0209618. [Google Scholar] [CrossRef] [PubMed]
  5. Flores-Vergara, A.; García-Guerrero, E.; Inzunza-González, E.; López-Bonilla, O.; Rodríguez-Orozco, E.; Cárdenas-Valdez, J.; Tlelo-Cuautle, E. Implementing a chaotic cryptosystem in a 64-bit embedded system by using multiple-precision arithmetic. Nonlinear Dyn. 2019, 1–20. [Google Scholar] [CrossRef]
  6. Li, X.; Xie, Z.; Wu, J.; Li, T. Image encryption based on dynamic filtering and bit cuboid operations. Complexity 2019, 2019, 7485621. [Google Scholar] [CrossRef]
  7. Praveenkumar, P.; Thenmozhi, K.; Rayappan, J.B.B.; Amirtharajan, R. Inbuilt Image Encryption and Steganography Security Solutions for Wireless Systems: A Survey. Res. J. Inf. Tech. 2017, 9, 46–63. [Google Scholar]
  8. Chen, G.; Mao, Y.; Chui, C.K. A symmetric image encryption scheme based on 3D chaotic cat maps. Chaos Solitons Fractals 2004, 21, 749–761. [Google Scholar] [CrossRef]
  9. Guan, Z.H.; Huang, F.; Guan, W. Chaos-based image encryption algorithm. Phys. Lett. A 2005, 346, 153–157. [Google Scholar] [CrossRef]
  10. Huang, X. A designed image encryption algorithm based on chaotic systems. J. Comput. Theor. Nanosci. 2012, 9, 2130–2135. [Google Scholar] [CrossRef]
  11. Wu, Y.; Yang, G.; Jin, H.; Noonan, J.P. Image encryption using the two-dimensional logistic chaotic map. J. Electron. Imaging 2012, 21, 013014. [Google Scholar] [CrossRef]
  12. Hua, Z.; Zhou, Y. Image encryption using 2D Logistic-adjusted-Sine map. Inf. Sci. 2016, 339, 237–253. [Google Scholar] [CrossRef]
  13. Liu, X.; Xiao, D.; Xiang, Y. Quantum image encryption using intra and inter bit permutation based on logistic map. IEEE Access 2019, 7, 6937–6946. [Google Scholar] [CrossRef]
  14. Flores-Vergara, A.; Inzunza-González, E.; García-Guerrero, E.E.; López-Bonilla, O.R.; Rodríguez-Orozco, E.; Hernández-Ontiveros, J.M.; Cárdenas-Valdez, J.R.; Tlelo-Cuautle, E. Implementing a Chaotic Cryptosystem by Performing Parallel Computing on Embedded Systems with Multiprocessors. Entropy 2019, 21, 268. [Google Scholar] [CrossRef]
  15. Wang, Y.; Wong, K.W.; Liao, X.; Chen, G. A new chaos-based fast image encryption algorithm. Appl. Soft Comput. 2011, 11, 514–522. [Google Scholar] [CrossRef]
  16. Pareek, N.K.; Patidar, V.; Sud, K.K. Image encryption using chaotic logistic map. Image Vis. Comput. 2006, 24, 926–934. [Google Scholar] [CrossRef]
  17. Hua, Z.; Jin, F.; Xu, B.; Huang, H. 2D Logistic-Sine-coupling map for image encryption. Signal Process. 2018, 149, 148–161. [Google Scholar] [CrossRef]
  18. Sahari, M.L.; Boukemara, I. A pseudo-random numbers generator based on a novel 3D chaotic map with an application to color image encryption. Nonlinear Dyn. 2018, 94, 723–744. [Google Scholar] [CrossRef]
  19. Zhou, N.; Yan, X.; Liang, H.; Tao, X.; Li, G. Multi-image encryption scheme based on quantum 3D Arnold transform and scaled Zhongtang chaotic system. Quantum Inf. Process. 2018, 17, 338. [Google Scholar] [CrossRef]
  20. Cvitanović, P.; Artuso, R.; Mainieri, R.; Tanner, G.; Vattay, G.; Whelan, N.; Wirzba, A. Chaos: Classical and Quantum; Niels Bohr Institute: Copenhagen, Denmark, 2005. [Google Scholar]
  21. Gangadhar, C.; Rao, K.D. Hyperchaos based image encryption. Int. J. Bifurcation Chaos. 2009, 19, 3833–3839. [Google Scholar] [CrossRef]
  22. Ye, G.; Wong, K.W. An image encryption scheme based on time-delay and hyperchaotic system. Nonlinear Dyn. 2013, 71, 259–267. [Google Scholar] [CrossRef]
  23. Chai, X.; Gan, Z.; Yang, K.; Chen, Y.; Liu, X. An image encryption algorithm based on the memristive hyperchaotic system, cellular automata and DNA sequence operations. Signal Process. Image Commun. 2017, 52, 6–19. [Google Scholar] [CrossRef]
  24. Zhang, L.M.; Sun, K.H.; Liu, W.H.; He, S.B. A novel color image encryption scheme using fractional-order hyperchaotic system and DNA sequence operations. Chin. Phys. B 2017, 26, 100504. [Google Scholar] [CrossRef]
  25. Bouslehi, H.; Seddik, H. Innovative image encryption scheme based on a new rapid hyperchaotic system and random iterative permutation. Multimed. Tools Appl. 2018, 77, 30841–30863. [Google Scholar] [CrossRef]
  26. Zhou, H.; Wilke, V.S. Research on image selective encryption and compression algorithm under hyperchaotic system. J. Intell. Fuzzy Syst. 2018, 35, 4329–4337. [Google Scholar] [CrossRef]
  27. Chai, X.; Gan, Z.; Lu, Y.; Zhang, M.; Chen, Y. A novel color image encryption algorithm based on genetic recombination and the four-dimensional memristive hyperchaotic system. Chin. Phys. B 2016, 25, 100503. [Google Scholar] [CrossRef]
  28. Li, X.; Chen, W.; Wang, Y. Quantum image compression-encryption scheme based on quantum discrete cosine transform. Int. J. Theor. Phys. 2018, 57, 2904–2919. [Google Scholar] [CrossRef]
  29. Wu, X.; Wang, D.; Kurths, J.; Kan, H. A novel lossless color image encryption scheme using 2D DWT and 6D hyperchaotic system. Inf. Sci. 2016, 349, 137–153. [Google Scholar] [CrossRef]
  30. Alfalou, A.; Brosseau, C.; Abdallah, N.; Jridi, M. Simultaneous fusion, compression, and encryption of multiple images. Opt. Express 2011, 19, 24023–24029. [Google Scholar] [CrossRef]
  31. Li, T.; Zhou, M. ECG classification using wavelet packet entropy and random forests. Entropy 2016, 18, 285. [Google Scholar] [CrossRef]
  32. Annaby, M.H.; Rushdi, M.A.; Nehary, E.A. Color image encryption using random transforms, phase retrieval, chaotic maps, and diffusion. Opt. Laser Eng. 2018, 103, 9–23. [Google Scholar] [CrossRef]
  33. Li, T.; Hu, Z.; Jia, Y.; Wu, J.; Zhou, Y. Forecasting Crude Oil Prices Using Ensemble Empirical Mode Decomposition and Sparse Bayesian Learning. Energies 2018, 11, 1882. [Google Scholar] [CrossRef]
  34. Zhou, Y.; Li, T.; Shi, J.; Qian, Z. A CEEMDAN and XGBOOST-based approach to forecast crude oil prices. Complexity 2019, 2019, 4392785. [Google Scholar] [CrossRef]
  35. Deng, W.; Zhang, S.; Zhao, H.; Yang, X. A novel fault diagnosis method based on integrating empirical wavelet transform and fuzzy entropy for motor bearing. IEEE Access 2018, 6, 35042–35056. [Google Scholar] [CrossRef]
  36. Zhang, Q.; Wang, Q.; Wei, X. A novel image encryption scheme based on DNA coding and multi-chaotic maps. Adv. Sci. Lett. 2010, 3, 447–451. [Google Scholar] [CrossRef]
  37. Wu, X.; Kan, H.; Kurths, J. A new color image encryption scheme based on DNA sequences and multiple improved 1D chaotic maps. Appl. Soft. Comput. 2015, 37, 24–39. [Google Scholar] [CrossRef]
  38. Faragallah, O.S.; Alzain, M.A.; El-Sayed, H.S.; Al-Amri, J.F.; El-Shafai, W.; Afifi, A.; Naeem, E.A.; Soh, B. Block-based optical color image encryption based on double random phase encoding. IEEE Access 2019, 7, 4184–4194. [Google Scholar] [CrossRef]
  39. Naeem, E.A.; Abd Elnaby, M.M.; Soliman, N.F.; Abbas, A.M.; Faragallah, O.S.; Semary, N.; Hadhoud, M.M.; Alshebeili, S.A.; Abd El-Samie, F.E. Efficient implementation of chaotic image encryption in transform domains. J. Syst. Softw. 2014, 97, 118–127. [Google Scholar] [CrossRef]
  40. Qian, Z.; Zhang, X.; Ren, Y. JPEG encryption for image rescaling in the encrypted domain. J. Vis. Commun. Image Represent. 2015, 26, 9–13. [Google Scholar] [CrossRef]
  41. Lima, J.B.; da Silva, E.S.; Campello de Souza, R.M. Cosine transforms over fields of characteristic 2: Fast computation and application to image encryption. Signal Process. Image Commun. 2017, 54, 130–139. [Google Scholar] [CrossRef]
  42. Wu, J.; Guo, F.; Liang, Y.; Zhou, N. Triple color images encryption algorithm based on scrambling and the reality-preserving fractional discrete cosine transform. Optik 2014, 125, 4474–4479. [Google Scholar] [CrossRef]
  43. Wu, J.; Zhang, M.; Zhou, N. Image encryption scheme based on random fractional discrete cosine transform and dependent scrambling and diffusion. J. Mod. Opt. 2017, 64, 334–346. [Google Scholar] [CrossRef]
  44. Chen, B.; Yu, M.; Tian, Y.; Li, L.; Wang, D.; Sun, X. Multiple-parameter fractional quaternion Fourier transform and its application in colour image encryption. IET Image Process. 2018, 12, 2238–2249. [Google Scholar] [CrossRef]
  45. Liu, X.; Xiao, H.; Li, P.; Zhao, Y. Design and implementation of color image encryption based on qubit rotation about axis. Chin. J. Electron. 2018, 27, 799–807. [Google Scholar] [CrossRef]
  46. Liansheng, S.; Xiao, Z.; Chongtian, H.; Ailing, T.; Asundi, A.K. Silhouette-free interference-based multiple-image encryption using cascaded fractional Fourier transforms. Opt. Laser Eng. 2019, 113, 29–37. [Google Scholar] [CrossRef]
  47. Fan, C.; Ding, Q. A novel image encryption scheme based on self-synchronous chaotic stream cipher and wavelet transform. Entropy 2018, 20, 445. [Google Scholar] [CrossRef]
  48. Lv, X.; Liao, X.; Yang, B. A novel scheme for simultaneous image compression and encryption based on wavelet packet transform and multi-chaotic systems. Multimed. Tools Appl. 2018, 77, 28633–28663. [Google Scholar] [CrossRef]
  49. Vaish, A.; Kumar, M. Color image encryption using MSVD, DWT and Arnold transform in fractional Fourier domain. Optik 2017, 145, 273–283. [Google Scholar] [CrossRef]
  50. Raja, S.P. Joint medical image compression-encryption in the cloud using multiscale transform-based image compression encoding techniques. Sadhana-Acad. Proc. Eng. Sci. 2019, 44, 28. [Google Scholar] [CrossRef]
  51. Hua, Z.; Zhou, Y. Design of image cipher using block-based scrambling and image filtering. Inf. Sci. 2017, 396, 97–113. [Google Scholar] [CrossRef]
  52. Hua, Z.; Xu, B.; Jin, F.; Huang, H. Image encryption using josephus problem and filtering diffusion. IEEE Access 2019, 7, 8660–8674. [Google Scholar] [CrossRef]
  53. Wu, Y.; Zhou, Y.; Noonan, J.P.; Agaian, S. Design of image cipher using latin squares. Inf. Sci. 2014, 264, 317–339. [Google Scholar] [CrossRef]
  54. Panduranga, H.T.; Kumar, S.K.N.; Kiran. Image encryption based on permutation-substitution using chaotic map and Latin Square Image Cipher. Eur. Phys. J. Spec. Top. 2014, 223, 1663–1677. [Google Scholar] [CrossRef]
  55. Xu, M.; Tian, Z. A novel image encryption algorithm based on self-orthogonal Latin squares. Optik 2018, 171, 891–903. [Google Scholar] [CrossRef]
  56. 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]
  57. Wang, H.; Li, X. A novel hyperchaotic system with infinitely many heteroclinic orbits coined. Chaos Solitons Fractals 2018, 106, 5–15. [Google Scholar] [CrossRef]
  58. Pano-Azucena, A.; Tlelo-Cuautle, E.; Rodriguez-Gomez, G.; de la Fraga, L. FPGA-based implementation of chaotic oscillators by applying the numerical method based on trigonometric polynomials. AIP Adv. 2018, 8, 075217. [Google Scholar] [CrossRef]
  59. Adleman, L.M. Molecular computation of solutions to combinatorial problems. Science 1994, 266, 1021–1024. [Google Scholar] [CrossRef]
  60. Zhu, C.; Hu, Y.; Sun, K. New image encryption algorithm based on hyperchaotic system and ciphertext diffusion in crisscross pattern. J. Electron. Inf. Tech. 2012, 34, 1735–1743. [Google Scholar] [CrossRef]
  61. Zhan, K.; Wei, D.; Shi, J.; Yu, J. Cross-utilizing hyperchaotic and DNA sequences for image encryption. J. Electron. Imaging 2017, 26, 013021. [Google Scholar] [CrossRef]
  62. Stinson, D.R. Cryptography: Theory and Practice; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  63. Shannon, C.E. Communication theory of secrecy systems. Bell Syst. Tech. J. 1949, 28, 656–715. [Google Scholar] [CrossRef]
  64. Wang, Z.; Huang, X.; Li, Y.; Song, X. A new image encryption algorithm based on the fractional-order hyperchaotic Lorenz system. Chin. Phys. B 2013, 22, 010504. [Google Scholar] [CrossRef]
  65. Wu, Y.; Noonan, J.P.; Agaian, S. NPCR and UACI randomness tests for image encryption. J. Sel. Areas Telecommun. 2011, 31–38. [Google Scholar]
Figure 1. The attractors of the 5D hyperchaotic system.
Figure 1. The attractors of the 5D hyperchaotic system.
Entropy 21 00319 g001
Figure 2. An example of Latin square of order 4.
Figure 2. An example of Latin square of order 4.
Entropy 21 00319 g002
Figure 3. An example of dynamic filtering.
Figure 3. An example of dynamic filtering.
Entropy 21 00319 g003
Figure 4. Three examples of Latin cube of order 3.
Figure 4. Three examples of Latin cube of order 3.
Entropy 21 00319 g004
Figure 5. A new cube K constructed by L 1 , L 2 , L 3 . K 1 , K 2 and K 3 are the 1st, 2nd and 3rd squares of K respectively.
Figure 5. A new cube K constructed by L 1 , L 2 , L 3 . K 1 , K 2 and K 3 are the 1st, 2nd and 3rd squares of K respectively.
Entropy 21 00319 g005
Figure 6. The framework of the proposed DFDLC.
Figure 6. The framework of the proposed DFDLC.
Entropy 21 00319 g006
Figure 7. Decrypted images of Lena, Cameraman, Barbara, Mandril and Bw with security keys g 1 and g 2 . The first and the second row is with g 1 and g 2 , respectively.
Figure 7. Decrypted images of Lena, Cameraman, Barbara, Mandril and Bw with security keys g 1 and g 2 . The first and the second row is with g 1 and g 2 , respectively.
Entropy 21 00319 g007
Figure 8. Histograms of the plain images and their corresponding cipher images. The first and the second columns are the plain images and their corresponding histograms, respectively. The third and the fourth columns are the cipher images and their corresponding histograms, respectively.
Figure 8. Histograms of the plain images and their corresponding cipher images. The first and the second columns are the plain images and their corresponding histograms, respectively. The third and the fourth columns are the cipher images and their corresponding histograms, respectively.
Entropy 21 00319 g008
Figure 9. The adjacent-pixel distribution maps of the plain images and the corresponding cipher images in horizontal direction.
Figure 9. The adjacent-pixel distribution maps of the plain images and the corresponding cipher images in horizontal direction.
Entropy 21 00319 g009
Table 1. Encoding and decoding rules of DNA.
Table 1. Encoding and decoding rules of DNA.
RULERule 1Rule 2Rule 3Rule 4Rule 5Rule 6Rule 7Rule 8
00ATTACGCG
01CGCGAATT
10GCGCTTAA
11TAATGCGC
Table 2. Addition (++) operation.
Table 2. Addition (++) operation.
++ACGT
ACATG
CACGT
GTGCA
TGTAC
Table 3. Subtraction (--) operation.
Table 3. Subtraction (--) operation.
--ACGT
ACGTA
CACGT
GTACG
TGTAC
Table 4. XOR ( ) operation.
Table 4. XOR ( ) operation.
ACGT
AACGT
CCATG
GGTAC
TTGCA
Table 5. Testing images.
Table 5. Testing images.
ImageSize ( h × w )ImageSize ( h × w )
Lena 512 × 512 Cameraman 512 × 512
Barbara 512 × 512 Mandril 512 × 512
Bw 512 × 512 Pirate 512 × 512
Couple 512 × 512 Finger 512 × 512
Peppers 512 × 512 Houses 512 × 512
Table 6. The IEs of the testing images.
Table 6. The IEs of the testing images.
ImageInputCipher Images
DFDLCFHDNA [2]HCDNA [61]CDCP [60]IC-BSIF [51]DFBC [6]
Lena7.44557.99937.99937.99947.99937.99947.9994
Cameraman7.04807.99927.99937.99817.99937.99937.9992
Barbara7.63217.99937.99947.99937.99927.99937.9993
Mandril7.29257.99947.99927.99927.99937.99937.9993
Bw1.00007.99937.99927.91587.99927.99937.9993
Pirate7.23677.99947.99937.99887.99937.99947.9993
Couple7.05727.99937.99927.99927.99937.99927.9993
Finger6.72797.99937.99947.99907.99927.99947.9993
Peppers7.59257.99937.99947.99917.99937.99937.9994
Houses7.65487.99927.99937.99937.99947.99947.9993
Table 7. The correlation coefficients γ of the testing images.
Table 7. The correlation coefficients γ of the testing images.
Image γ InputCipher Images
DFDLCFHDNA [2]HCDNA [61]CDCP [60]IC-BSIF [51]DFBC [6]
Lena γ h 0.96910.00230.0000−0.0015−0.0004−0.00320.0002
γ v 0.98410.0009−0.0022−0.00200.00280.00130.0010
γ d 0.96390.00080.00040.00240.0016−0.00090.0006
Cameraman γ h 0.98300.00110.00130.0004−0.0001−0.0015−0.0008
γ v 0.98870.00090.00330.00030.00190.0010−0.0013
γ d 0.9746−0.0002−0.0000−0.00130.0010−0.0012−0.0002
Barbara γ h 0.8940−0.0003−0.00220.0010−0.0026−0.00020.0027
γ v 0.95720.0030−0.00020.00040.0006−0.0004−0.0029
γ d 0.8942−0.0029−0.0000−0.00090.00050.0010−0.0005
Mandril γ h 0.93220.00220.0016−0.00070.00120.0026−0.0006
γ v 0.91000.00050.0035−0.00010.0009−0.0001−0.0018
γ d 0.8647−0.0023−0.0025−0.0017−0.00040.00010.0016
Bw γ h 1.00000.00190.00060.0004−0.00040.00030.0000
γ v 0.9922−0.00060.00090.00130.0001−0.0005−0.0002
γ d 0.9961−0.0012−0.0012−0.00020.00050.0002−0.0016
Pirate γ h 0.9593−0.00000.0015−0.0023−0.0012−0.0026−0.0012
γ v 0.96750.00090.0057−0.0000−0.0008−0.00060.0013
γ d 0.94320.00150.00010.00110.00060.00050.0005
Couple γ h 0.94510.00120.00130.0014−0.0001−0.0006−0.0009
γ v 0.95140.0025−0.00260.00080.00010.00230.0022
γ d 0.91160.0017−0.0011−0.00070.0005−0.0008−0.0024
Finger γ h 0.9343−0.00010.00020.0007−0.00230.0004−0.0025
γ v 0.91680.0002−0.00250.0029−0.0032−0.00090.0004
γ d 0.86640.00170.0005−0.0022−0.00100.0030−0.0006
Peppers γ h 0.97330.0003−0.00450.0000−0.0003−0.00310.0008
γ v 0.9763−0.0010−0.0049−0.00050.0003−0.0010−0.0003
γ d 0.96500.0011−0.0012−0.0005−0.00250.0017−0.0010
Houses γ h 0.90770.00200.00060.00040.00260.0001−0.0002
γ v 0.91730.00150.0004−0.00320.00020.00170.0006
γ d 0.84390.00200.00210.0038−0.00110.00340.0002
Range[0.8439,1.000][−0.0023,0.0030][−0.0049, 0.0057][−0.0032, 0.0038][−0.0032, 0.0028][−0.0032, 0.0034][−0.0029, 0.0027]
Interval Width0.15610.00530.01060.00700.00600.00660.0056
Table 8. The mean/standard deviation/times of passing the test of NPCR (%) of running the schemes 10 times ( α = 0.05 ).
Table 8. The mean/standard deviation/times of passing the test of NPCR (%) of running the schemes 10 times ( α = 0.05 ).
ImageDFDLCFHDNA [2]HCDNA [61]CDCP [60]BSIF [51]DFBC [6]
Lena99.6103/0.0129/1099.5814/0.0119/643.5948/16.8360/099.6201/0.2837/599.6166/0.0109/1099.5995/0.0002/10
Cameraman99.6055/0.0126/1099.5795/0.0137/464.6306/31.1442/099.6146/0.2372/699.6057/0.0121/1099.6143/0.0002/10
Barbara99.6171/0.0075/1099.5842/0.0099/837.8473/19.6663/099.6048/0.2136/699.6165/0.0144/1099.5833/0.0002/10
Mandril99.6047/0.0117/1099.5774/0.0125/351.2024/28.3679/099.5697/0.2107/499.6070/0.0117/1099.5998/0.0001/10
Bw99.6030/0.0094/1099.3196/0.2433/147.5142/15.7628/099.6362/0.2321/599.6180/0.0142/1099.6033/0.0000/10
Pirate99.6176/0.0133/1099.5812/0.0127/435.8150/27.9995/099.6403/0.3222/799.6116/0.0116/1099.5751/0.0002/0
Couple99.6089/0.0133/1099.5779/0.0076/558.1698/27.5116/099.5718/0.1939/399.6079/0.0112/1099.5586/0.0001/0
Finger99.6097/0.0158/1099.5792/0.0152/460.3329/29.8886/099.5984/0.1420/699.6171/0.0096/1099.6132/0.0002/10
Peppers99.6099/0.0154/999.5800/0.0119/545.6316/38.2206/099.5493/0.2509/499.6099/0.0150/999.6166/0.0001/10
Houses99.6130/0.0126/1099.5795/0.0083/663.0733/19.2267/099.6039/0.2053/799.6135/0.0119/1099.6151/0.0001/10
Table 9. The mean/standard deviation/times of passing the test of UACI (%) of running the schemes 10 times ( α = 0.05 ).
Table 9. The mean/standard deviation/times of passing the test of UACI (%) of running the schemes 10 times ( α = 0.05 ).
ImageDFDLCFHDNA [2]HCDNA [61]CDCP [60]BSIF [51]DFBC [6]
Lena33.4504/0.0466/933.2700/0.0490/018.5974/9.5490/033.5212/0.0775/633.4714/0.0339/1033.4818/0.0005/10
Cameraman33.4909/0.0457/933.3010/0.0320/027.0047/13.9227/033.4222/0.0658/733.4755/0.0485/1033.4406/0.0005/10
Barbara33.4451/0.0350/1033.2533/0.0431/013.6480/8.2289/033.4464/0.1075/733.4722/0.0476/933.4808/0.0007/10
Mandril33.4704/0.0334/1033.2988/0.0336/022.2006/14.1993/033.4467/0.0928/533.4449/0.0423/1033.5136/0.0003/10
Bw33.4334/0.0471/1032.0705/1.0272/018.5654/6.1072/033.4555/0.1144/433.4500/0.0468/1041.6585/0.0010/0
Pirate33.4736/0.0275/1033.3021/0.0431/114.8888/14.1945/033.4664/0.0766/833.4644/0.0328/1033.4668/0.0003/10
Couple33.4282/0.0385/933.2796/0.0381/017.8782/7.3362/033.4293/0.1011/933.4632/0.0439/1033.4717/0.0003/10
Finger33.4311/0.0504/833.2907/0.0413/026.0775/14.9380/033.4911/0.1004/633.4856/0.0399/933.5263/0.0006/10
Peppers33.4618/0.0432/1033.2735/0.0347/019.8106/18.6275/033.4626/0.0752/633.4301/0.0379/1033.4525/0.0009/10
Houses33.4634/0.0358/1033.3322/0.0273/122.1296/7.4892/033.4721/0.0467/1033.4448/0.0343/1033.4545/0.0004/10

Share and Cite

MDPI and ACS Style

Li, T.; Shi, J.; Li, X.; Wu, J.; Pan, F. Image Encryption Based on Pixel-Level Diffusion with Dynamic Filtering and DNA-Level Permutation with 3D Latin Cubes. Entropy 2019, 21, 319. https://doi.org/10.3390/e21030319

AMA Style

Li T, Shi J, Li X, Wu J, Pan F. Image Encryption Based on Pixel-Level Diffusion with Dynamic Filtering and DNA-Level Permutation with 3D Latin Cubes. Entropy. 2019; 21(3):319. https://doi.org/10.3390/e21030319

Chicago/Turabian Style

Li, Taiyong, Jiayi Shi, Xinsheng Li, Jiang Wu, and Fan Pan. 2019. "Image Encryption Based on Pixel-Level Diffusion with Dynamic Filtering and DNA-Level Permutation with 3D Latin Cubes" Entropy 21, no. 3: 319. https://doi.org/10.3390/e21030319

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