Image Encryption Based on Dynamic Filtering and Bit Cuboid Operations

As one of the most widely used media types, images play an important role in the era of the Internet. And hence how to enhance the security of images has become a hot topic in the field of information security. However, due to some intrinsic characteristics of images, image security is still a challenging task. For the purpose of coping with this issue, in this paper, we propose a novel algorithm that combines a hyperchaotic system, dynamic filtering, and bit cuboid operations, namely, DFBC, for image encryption. Specifically, the proposed DFBC consists of four steps: firstly, a 7D Lorenz hyperchaotic system is utilized to generate a pseudorandom sequence; secondly, variable 1D filters are derived from the pseudorandom sequence, and dynamic filtering is conducted on each pixel of an image; thirdly, a diffusion scheme is performed and then the image is transformed to a bit cuboid; and, finally, various types of permutation (rearranging, symmetry, rotation, zigzag, and global bit permutation) are performed on the bit cuboid. The experiments on several testing images demonstrate that the DFBC achieves state-of-the-art results in terms of several evaluation criteria, showing that the DFBC is promising for image encryption.


Introduction
Image security is one of the most important types of information security.In recent years, image encryption has become a hot topic in the field of image security.However, since images usually show some intrinsic properties such as bulky data capacity, high redundancy, and strong correlation, traditional encryption algorithms such as data encryption standard (DES), advanced encryption standard (AES), and international data encryption algorithm (IDEA) for common data are not able to be applied to images directly to achieve good results [1,2].To cope with this issue, various researchers have been devoted to proposing image encryption algorithms to enhance image security in recent years.
Among the image encryption algorithms, chaos-based ones have become more and more popular in recent years because chaotic systems have some special properties for image encryption, such as pseudorandomness, ergodicity, unpredictability, and extreme sensitivity to system parameters and initial values [3][4][5][6][7].Generally speaking, there are two ways to perform image encryption, i.e., permutation and diffusion.The former means changing the position of image data, while the latter means changing the values of image data.Many practical image encryption algorithms use these two approaches together.When chaotic systems are applied to image encryption, the common ways are to perform permutation and diffusion using the index and the values of chaotic sequences generated from chaotic systems, respectively.Chen et al. proposed a novel real-time secure symmetric image encryption scheme that generalized a 2D chaotic map to 3D and then used the 3D cat map and another chaotic map to conduct permutation and diffusion, respectively [3].Zhang proposed a novel fast image cryptosystem that used the piecewise linear chaotic map and cubic S-box to generate the sequence for image encryption [6].Typically, low-dimensional chaotic systems have simple forms and are 2 Complexity easy to implement, and hence some scholars applied them to image encryption [8,9].
However, low-dimensional chaotic maps usually have only a few variables and parameters, along with simple structures and chaotic orbits, making it easy to estimate the orbits and the initial parameters, and hence the security of image is reduced [2,10].Therefore, higher-dimensional chaotic systems have been widely applied to image encryption.Some researchers apply 3D chaotic maps to image encryption [11,12].In a dynamical system, the Lyapunov exponent (LE) is a quantity that characterizes the rate of separation of infinitesimally close trajectories and estimates the chaos of the system [13], and if a chaotic system has two or more positive LEs, the system is defined to be hyperchaotic [14].Recent research has demonstrated that image encryption algorithms associated with hyperchaotic systems show greater security [15][16][17][18][19][20][21][22][23].Among the approaches, 4-6D hyperchaotic systems are widely used to enhance image encryption.Zhu and Sun proposed an image encryption algorithm with two-round encryption operations based on a 4D hyperchaotic system [16].Xue et al. applied a 5D hyperchaotic system with 3 positive LEs for region of interest encryption for color images [21].Wu et al. presented a new lossless encryption algorithm that used the 2D discrete wavelet transform (DWT) and a 6D hyperchaotic system in both frequency and spatial domains for color images, and the experimental results indicated that the proposed algorithm was effective and efficient [23].Generally speaking, the more the number of positive LEs or the higher the dimension of a hyperchaotic system is, the more secured the hyperchaotic system-based image encryption algorithms are.
When an image encryption algorithm is performed on an image, it usually processes block-level data (a block of pixels) [24][25][26][27], pixel-level data [28,29], deoxyribonucleic acid-(DNA-) level data (2 bits) [30][31][32][33], bit-level data [34][35][36], or bit plane-level data [37,38].For fixed processing power, the lower the processing level is, the more pixels the algorithm can handle in one time.For example, if the algorithm can process 8 bits in one time for 256-level gray images, it means that it can process 1 pixel, 4 nucleic acid bases, or 8 bits, which involves 1/4/8 pixels at the most, respectively.Therefore, encryption algorithms associated with low-level data usually show better performance in image encryption.
Image filtering, also known as convolution, is one of the most important operations in image processing, and its typical applications include denoising, smoothing, and edge detection.Most recently, Hua and Zhou have applied this technique to image encryption for the first time, and they proposed an approach called block-based scrambling and image filtering (IC-BSIF) for image encryption.The experimental results showed that it could achieve better performance than some state-of-the-art encryption schemes [26].In this scheme, the authors use a fixed template of filter for all pixels, losing the diversity of the template and hence not maximizing the performance of the scheme.In the field of chaos, Yang et al. have found a new 7D autonomous hyperchaotic system with 5 positive LEs recently, which has very simple algebraic structure but can show complex dynamical behaviors [39].
Therefore, it has good potential for image encryption.Motivated by the above analysis, in this paper, we aim to propose a novel image encryption scheme that combines a 7D hyperchaotic system with 5 positive LEs, dynamic filtering, and bit cuboid operation, namely, DFBC, for image encryption.The main steps of the DFBC are as follows.(1) A 7D hyperchaotic system is used to generate a pseudorandom sequence for subsequent encryption operations.(2) 1D dynamic filtering is conducted on each pixel with random filters derived from the pseudorandom sequence.(3) The 2D image of pixel plane is transformed to a 3D bit cuboid.(4) Various types of permutation, such as rearranging, symmetry, rotation, zigzag, and global bit permutation, are performed on the bit cuboid.
The rest of this paper is structured as follows.A brief description of the 7D hyperchaotic system and filtering is given in Section 2. In Section 3, the novel image encryption algorithm, DFBC, is proposed in detail.Experimental results are reported in Section 4. Finally, we conclude the paper in Section 5.
The advantages of this hyperchaotic system are two aspects: (1) the algebraic structure is very simple, so it is easy to implement and (2) the system exhibits complex dynamical behaviors, and hence it is suitable for image encryption.Therefore, in this paper, we use this hyperchaotic system for image encryption.

Filtering Operation.
It is well-known that image filtering by convolution is widely used in many digital image processing fields, such as adding soft blur, sharpening details, accentuating edges, or denoising.Some of filtering is linear, which involves weighted combinations of pixels in small neighborhoods.The simplest filter to implement is the moving average or box filter, which simply averages the pixel values in an × window.Usually, the gradient, Laplacian, and difference of Gaussians (DoG) filters work very well for general purpose visual enhancement or edge detection.This is equivalent to convolution of the image with a mask of all ones.
Traditional image filtering cannot be reversed; that is to say, it usually cannot recover the original image with the filtered image and the filters.To overcome this shortcoming, very recently, Hua and Zhou have designed a special filter with which the original image can be recovered from the filtered image.And they have proposed a novel approach called IC-BSIF for image encryption, where image filtering Complexity is applied to image encryption for the first time [26].The experimental results demonstrated the effectiveness of the proposed IC-BSIF.

Bit-Level Operation.
Permutation conducted on pixellevel data or higher level data (e.g., blocks of pixels) cannot change the statistical characteristics of an image [45], while permutation on bit-level data can change such characteristics.Furthermore, for a fixed process width (bits processed in one run), the bit-level image encryption can involve more pixels than higher level data (e.g., DNA-level, pixel-level, or blocks of pixels), for both permutation and diffusion.Therefore, to enhance the security of image encryption and to speed up the algorithm, in this paper, we will perform bit-level encryption.Specifically, we first transform an image into a cuboid that is very close to a cube, and then various bit plane operations are conducted on the cuboid.

The Proposed Image Encryption Approach
3.1.Hyperchaotic Sequence Generation.Since the 7D hyperchaotic system described in Section 2.1 has good properties for image encryption, we use it to generate the hyperchaotic sequence applied in the proposed image encryption approach.The generating process consists of three steps.
Step 1.To eliminate the adverse effects, the 7D hyperchaotic system is firstly iterated  0 times and then the sequence generated is removed.
Step 2. The 7D hyperchaotic system continues to iterate until the generated sequences are long enough to image encryption.For the −th iteration, seven state values denoted by   = { Step 3.After the whole iteration, the hyperchaotic sequences  can be obtained by contacting all the   ( = 1, 2, ⋅ ⋅ ⋅ , ) as The purposes of the generated sequence  for encryption are two aspects: (1) sorting subsequence of  to get the index of original data for permutation; (2) using subsequence of  to change the values of images for diffusion.In our approach, we map the −th point in  to the integral range of [0,255] by (3) and to the real range of [0, 1) by (4) for the first purpose and the second purpose, respectively.
where mod is the modulo operation, | ⋅ | is the absolute value operation, and ⌊⋅⌋ denotes flooring operation [2].

Image Filtering.
Image filtering is the convolution operation between an image and a template , also known as a mask.This template is usually a 2D matrix whose size is 3 × 3, 5 × 5, 7 × 7, etc. Suppose the size of a 2D template is (2 + 1) × (2 + 1), the convolution operation of image filtering can be defined as If the template is a one-row or one-column matrix, it is called 1D convolution, which can be defined as the following: where the sizes of  in ( 6) and ( 7) are 1×(2+1) and (2+ 1) × 1, respectively.

1D Dynamic
Filtering.Hua and Zhou used 2D filtering to do diffusion in image encryption.They also introduced Proposition 1 and its proof to identify that image filtering operation satisfying some conditions can be reversible [26].
The key difference to the ordinary image filter is to set the mask center to 1; i.e.,  +1,+1 = 1.
Here we can prove that if the center of 1D image filter is set to 1, Proposition 1 still holds.And it is rewritten as Proposition 2.
Proposition 2. For any given 1D mask  of size (2 + 1) × 1, image  of size  × , and P's grayscale level , the operation is reversible and its inverse operation is if  ∈ N,  ∈ N and  +1 = 1.

Complexity 5
For the mask with direction along with coordinate , or, in other words, 1D mask  of size 1 × (2 + 1), 1D image filtering has the same property.
By now, the proof of Proposition 2 has been completed.It has the similar formation with (9).As we can see that, for Proposition 2, even if every pixel  , is filtered by different template  , ++1 , Proposition 2 still holds and then is detailed as Proposition 3 as follows.
Proposition 3.For any given mask  , of size (2 + 1) × 1 at pixel  , of P's grayscale level , the operation is reversible where image size is size  ×  and its inverse operation is The logic of setting our masks  , ++1 is described in Figure 3.
The colored  , +1 is the center of the mask of filtering which is 1 exactly.It pairs with the pixel  , which needs to be filtered.The other weights in the mask after the colored  , +1 are set to 0. Then this 2 + 1 mask is simplified to  + 1 and the weight with respect to the pixel  , is the last weight in the mask.In the proposed approach, the mask  , is said to be a dynamic one because, for each pixel  , , we use different values from the hyperchaotic sequence to fill the mask as the weights.
When filtering is conducted on the edge of an image, the original image has to be expanded [26].Since the filter is 1D, we only need to expand the image on horizontal or vertical direction.When horizontal filter is working, the most right  columns are copied to the left side of the image.Similarly, when vertical filter is working, the most bottom  rows are copied to the top side of the image.Of course, the number of rows or columns which need to be copied in filtering process has to change accordingly at the row or column between 1 and .This process is depicted by an example in Figure 4.The mask size in the example is set as 3 × 1 and 1 × 3.
Our dynamic diffusion mask of filtering comes from the hyperchaotic sequence.And the indicator of direction of filtering mask is determined by one number in the hyperchaotic sequence too.If the number is even, a horizontal mask is applied to the image.Otherwise, an odd number means a vertical mask.
After the dynamic masks are set, the convolution begins as the normal convolution does except that the mask has to be renewed for each pixel.

Global Pixel Diffusion.
Although the dynamic filtering is capable of diffusing images, to further enhance the performance of the diffusion of our approach, we conduct a simple two-step diffusion for image with pixels after dynamic filtering, following the diffusion scheme in [2].Given an image of size  =  × , where  and  are the height and width, respectively, the scheme firstly transforms the image into a 1D pixel sequence  = {  },  = 1, 2, ⋅ ⋅ ⋅ , .Given an initial key  0 and a chaotic key sequence  = {  ∈ [0, 255]},  = 1, 2, ⋅ ⋅ ⋅ , , the diffusion can be conducted in two steps, as formulated in the following equations, respectively: where ⊗ is the XOR operation and   is −th pixel in the diffused image.
Besides the dynamic filtering and global pixel diffusion, we also conduct the XOR operation on the pixels of image with a random mask to change the value of the pixels directly.

Bit Cuboid Operations
3.4.1.Pixel Plane to Bit Cuboid.Given a plain grayscale image with  rows (height) and  columns (width), it is a typical  pixel plane with  ×  pixels, each of which is one 256-level value.The total bits of the image are ××8, so it is naturally a cuboid.To perform bit-level operations better, the cuboid should be transformed into a new one whose width, height, and depth can be as uniform as possible.To achieve this, we transform the original cuboid with size of  ×  × 8 into a new one with size of   ×   × , as Equation ( 19) can ensure that both   and   are integers.For the simplicity, in this paper, we do not consider the cases when  is not integral.With this equation, the original pixel plane can be transformed into a cuboid which has as uniform a side length as possible.For example, grayscale images having sizes of 512 × 512 and 256 × 256 can be transformed into two cuboids with sizes of 128 × 128 × 128 and 64 × 64 × 128, respectively.

Operations on Bit-Planes
(1) Bit Plane Rearranging.Bit plane rearranging can be conducted at different directions for a bit cuboid.Since the operations at each direction are almost the same, here we use the operation at the direction of depth as an illustration for simplicity.The rearranging steps are as follows.
(4) Zigzag.Zigzag can disturb the high correlation for images to enhance the security of encrypted images [46].In this paper, we use zigzag to permutate the bit-planes in the bit cuboid one by one, as illustrated in Figure 5.In each bit plane, the zigzag confusion path starts from the upper left bit and ends at the bottom right bit.When scanning a plane is completed, the scan starts from the upper left bit of the next plane, and so on.When all the bit-planes are scanned, a bit string is obtained.Zigzag can not only permute the bits, but also diffuse corresponding pixel-level data.
(5) Global Bit Permutation.After Zigzag, the encryption scheme will obtain a bit sequence.The global bit permutation is to permute all bits in the bit sequence with a hyperchaotic sequence [2,30].The steps are as follows.
Step 2. According to the index sequence, rearrange the original -th bit in the bit sequence as the   -th one in the permuted bit sequence.
With the operation of global bit permutation, the bit sequence can be transformed to a pixel plane as the encrypted image.It also can be transformed into a bit cuboid for further bit-level permutation.

The DFBC: The Proposed Image Encryption Using Hyperchaotic System, Dynamic Filtering and Bit Cuboid Operation.
The framework of the DFBC is shown in Figure 6.After generating the hyperchaotic sequence, the encryption procedure consists of three steps: pixel-level diffusion (dynamic filtering and global pixel diffusion), transformation of pixel plane to bit cuboid, and bit-level permutation.The diffusion changes the gray value of every pixel while the permutation rearranges the positions of bits.
The detailed encryption procedure of the DFBC is described as follows.
Step 2. Use the dynamic masks to diffuse the image  by the method in Section 3.2.1.Firstly, the first number of the hyperchaotic sequence,  0 , is chosen to determine the direction of mask.Then,  ×  ×  integers are chosen from the hyperchaotic sequence.Secondly, reshape these integers into  × masks of 1× or ×1,  , , where the direction of mask is determined by the parity of  0 .Set the last element of the mask to one; i.e.,  ,  = 1.Thirdly, convolve each mask  , with its corresponding pixels to get a diffused image.
Step 3. Conduct global pixel diffusion by ( 17)-( 18) as described in Section 3.3.The operation will expand a little change in one pixel to other pixels.
Step 4. Conduct XOR operation on the image with a hyperchaotic subsequence to change the values of pixels directly.
Step 5. Transform the pixel plane to a bit cuboid as described in Section 3.4.1.The transformed bit cuboid has as uniform width, height, and depth as possible.
Step 6. Rearrange the bits in the bit cuboid based on bit plane at height, width, and depth direction as described in Section 3.4.2.(1).
Step 7. Randomly rotate each bit-plain based on chaos sequence as described in Section 3.4.2.(3).
Step 8. Randomly symmetrize each bit-plain based on chaos sequence as described in Section 3.4.2.(2).
Step 10.Conduct global bit permutation on the bit sequence as described in Section 3.4.2.(5).
Step 11.Transform the bit sequence to a pixel plane of the original size as the final encrypted image.
The procedure of the DFBC can be categorized into five stages: hyperchaotic sequence generation (Step 1), pixel-level diffusion (Step 2-4), pixel plane to bit cuboid transformation (Step 5), bit-level permutation (Steps 6-10), and bit sequence to pixel-plane transformation (Step 11).The DFBC diffuses images with dynamic filtering, global pixel diffusion, and XOR operation on pixel-level data.And it permutes the images on bit-level data.It is worth pointing out that since the permutation on bits can change the values of pixels, the bit cuboid operations in the DFBC are capable of diffusing the images; that is to say, the bit cuboid operations have the effects of both permutation and diffusion.
The decryption process of DFBC is the inverse procedure of encryption as listed above.

Experimental Settings.
In order to evaluate the performance of the proposed DFBC, some state-of-the-art schemes, such as the hyperchaotic and DNA sequence-based method (HC-DNA) [30], the image encryption using encrypted diffusion in crisscross pattern (CDCP) [16], a class hyperchaosbased scheme (CHC) [17], an image encrypted using blockbased scrambling, and image filtering (IC-BSIF) [26], are compared with some common statistical indexes.The parameters for the compared schemes are set as the corresponding papers.The parameters of the DFBC are set as follows.The initial values for the 7D hyperchaotic system are  0 1 = 0,  0 2 = 0.1,  0 3 = 0.2,  0 4 = 0.3,  0 5 = 0.5,  0 6 = 0.6, and  0 7 = 0.7.And the preiterating times  0 is set to 1000.The original hyperchaotic integer sequence generated from the 7D hyperchaotic equations is adopted from the beginning, which means the start position of hyperchaotic sequence is 1.Then our diffusion masks of filtering are created from this hyperchaotic integer sequence.The first byte of the hyperchaotic sequence determines the direction of 1D filter.Now that the image size is  ×  and mask size is 3 × 1 or 1 × 3, then  ×  × 3 integers of one byte are cut out from the hyperchaotic sequence to shape  ×  masks to filter every pixel  , iteratively.The round of 1D filtering and cuboid operation are both set to 1.
Six publicly accessed images and two images we collected, Monkey and Parking, with different sizes are used to test the proposed DFBC, as listed in Table 1.
The publicly accessed images Lena, Cameraman, Barbara, and Airfield are widely used in the field of image processing, including image encryption.The image Histogram is a graphical one which has an almost similar grayscale background and histogram which means most of its pixels are highly correlated.This image is quite different from other real photographs from cameras.The image Terrain is an image downloaded from the Internet which has a very single sharp mountain shape histogram.The image Parking is captured by the author.It has a big area of gray cloud background which has a very high correlation between pixels.The image Monkey is a real image captured from a video which is different from these well-known images.It is to validate that our DFBC can be applied to any real images robustly and easily.
All the experiments were conducted by Matlab 8.3 (Mathworks, Natick, MA, USA) on a 64-bit Windows 7 (Microsoft, Redmond, WA, USA) with 8 GB memory and 3.4 GHz I3 CPU.

Security Key Analysis.
A good encryption scheme should have an enough large key space and be extremely sensitive to any small changes in its security key.Both a large key space and extreme sensitivity can resist brute-force attacks.So key space and sensitivity to secret key are two essential points in encryption.

Key Space.
Basically, the security keys of the proposed DFBC are composed of 7 initial values, i.e., ( 0 1 ,  0 2 ,  0 3 ,  0 4 ,  0 5 ,  0 6 ,  0 7 ).If the precision of each initial value is 10 −15 , the key space size is 10 15×7 = 10 105 ≈ 2 348 .From the view of cryptology, the size of the key space larger than 2 100 can provide a high-level security [1,47].Therefore, the key space of the DFBC is large enough to resist all kinds of brute-force attacks from current computers.In addition, the start position of chaos sequence to form the random filter mask can also be used as a key to further enhance the key space.

Sensitivity to Security Key.
The extreme sensitivity of an image encryption algorithm implies that even one bit changed in the keys will lead to a completely different encrypted image.In other words, if the security key changes a little, the decrypted image will be completely different from the input image.
To demonstrate the sensitivity to secret key of the DFBC, we decrypt the encrypted images twice.In the first run, we use the exact encryption keys ( 0 1 = 0,  0 2 = 0.1,  0 3 = 0.2,  0 4 = 0.3,  0 5 = 0.5,  0 6 = 0.6, and  0 7 = 0.7) to decrypt the encrypted images, while, in the second run, we attempt to decrypt the encrypted images with slightly different keys ( 1  1 = 0 + 10 −15 ,  1 2 = 0.1,  1 3 = 0.2,  1 4 = 0.3,  1 5 = 0.5,  1 6 = 0.6, and  1 7 = 0.7).We conduct the experiments on the images of Lena, Cameraman, Monkey, and Parking, and the results are shown in Figure 7.As we can see, even if we only change the key extremely little such as 10 −15 , the decrypted images are completely different from the ones decrypted with the correct keys, validating that the proposed DFBC is highly sensitive to secret key.

Statistical Analysis. The statistical analysis is another widely used and effective way to analyze a cryptosystem.
We adopt some typical statistical analysis, such as histogram analysis, information entropy, and correlation analysis, to evaluate the performance of the proposed DFBC.A cryptosystem with good performance of statistical analysis can resist all kinds of statistical attacks.The experimental analysis demonstrates that our method DFBC also can handle images especially like the graphical Histogram very well.

Histogram Analysis.
Histogram is a popular and effective way to measure the distribution of pixel values in the plain image and the encrypted image for image encryption.The histogram of a plain image is usually unevenly distributed while that of an encrypted image by a good encryption scheme has a uniform distribution as much as possible.A uniform distribution of histogram indicates a totally randomlike image and the least probability of recovering its plain image.In other words, in terms of the performance of encryption schemes, the flatter of the histogram of the encrypted image is, the better the encryption scheme is.
The histograms of the plain images and their corresponding encrypted images are shown in Figure 8.The histograms of the plain images are shaped like some mountains or valleys while all the histograms of encrypted images are very close to a uniform distribution.These histograms of encrypted images indicate that the proposed DFBC have the ability to resist histogram attacks.The fifth image Histogram especially has a highly unusual histogram, where most of the pixels have the highest gray scale levels so that the rectangles of low gray scales are very short.And our method can encrypt the image Histogram to be an even and flat histogram like the other photographic images.It confirms that our method works very well for man-made graphical images like Histogram.
The diffused images by dynamic filtering are listed in the third column of Figure 8, and the fourth column is their corresponding histograms.It is worth noting that the diffused images by dynamic filtering have almost uniformly distributed histograms like the encrypted images.It can be seen that their histograms are just a little less even than the encrypted images.Compared with the diffusion in [26], our diffusion is as good as [26] at the same level with regard to the histogram analysis.Theoretically, our 1D image filtering is more efficient and can compute faster than 2D filtering with a fixed mask.

Information Entropy.
Information entropy (IE) is the average rate at which information is produced by a stochastic source of data.Here it is used to reflect the complexity or orderliness of the encryption system.The intensity of an 8bit grayscale image has 2 8 possible values ([0, 255]).The IE is defined as where (  ) is the probability that the pixel value   appears [30].When   of encrypted image has the same probability, i.e., 1/256, IE reaches maximum of 8.It matches the uniform distribution of the encrypted image perfectly.The IEs of input images and encrypted images by different encryption methods are listed in Table 2.The third column under DF (Dynamic Filtering) is the entropies of diffusion by dynamic filtering.The fourth column under DFBC is the entropies of final encrypted images.All the entropies of diffused images except Barbara are improved after our bit cuboid operations, which indicates that this random permutation is very effective.It can be seen that the IEs of input images are far below 8, while those of encrypted images are very close to the ideal value.The IEs of DFBC are within [7.9972, 7.9993], demonstrating that the DFBC is secure enough to resist entropy attack.Our DFBC achieves 5 out of 8 optimal values while HC-DNA and IC-BSIF achieve only 2 and 1 out of 8 optimal values, respectively.Although the IEs by DF are not the highest ones, they are very close to the ideal values, showing that the 1D dynamic filtering can achieve very good results of diffusion.

Correlation Analysis.
Two neighboring pixels in a natural image usually have similar grayscale levels and thus are highly correlated.A good image encryption algorithm should result in encrypted image with extremely low correlation.Correlation coefficient  is a widely used metric to measure the correlation and it can be formulated as [44] where  and  are grayscale levels of two neighboring pixels in an image and  denotes the total number of pixels in an image that satisfies  =  × .
As listed in Table 3, we calculate the correlation coefficients for all input images and encrypted images at three different directions, i.e., horizontal  ℎ , vertical  V , and diagonal   , respectively [30].It can be seen from this table that the correlation coefficients of all the input images are big number and close to 1 in all directions, while all the encrypted images are around 0, showing that the encryption schemes can effectively reduce the correlation of the adjacent pixels to a very low level after encryption.The DFBC outperforms the remaining schemes on 6 out of 24 correlation coefficients, whereas the CHC also achieves the optimal value 6 times and IC-BSIF achieves 5 times, showing that the DFBC has good performance for correlation of the encrypted images.When we further investigate the range of , however,  of DFBC is within [−0.0041, 0.0040], while those of HC-DNA, CDCP, CHC, and IC-BSIF are within [−0.0091, 0.0076], [−0.0054, 0.0048], [−0.0069, 0.0048], and [−0.0067, 0.0104], respectively.From the ranges of the algorithms, we can get the intervals of , i.e., 0.0081, 0.0167, 0.0102, 0.0117, and 0.0171, for  The number of pixels change rate (NPCR) and the unified average changing intensity (UACI) are two important metrics for differential attack analysis.The NPCR is defined as the variation ratio of two encrypted images whose plain images are slightly different, meaning the dissimilitude between two encrypted images.The UACI indicates the average intensity of the differences between the same encrypted images.In most cases, only one pixel or even one bit chosen is randomly changed to compute the NPCR and UACI.Mathematically, NPCR and UACI between two encrypted images  1 and  2 can be formulated as follows, respectively.
where  and  are the height and the width of the image, respectively, and   is defined as NPCR concentrates on the absolute number of pixels which change values in differential attack, while the UACI focuses on the average difference between the corresponding encrypted images.For an image with 256 grayscale levels, the expectations of NPCR and UACI of an ideally encrypted image are 99.6094% and 33.4635%, respectively [48].Generally speaking, the more NPCR gets close to 100% and the bigger UACI is, the more effective it is for the encryption scheme to resist differential attacks.
We randomly change one pixel and then add 1 on gray value in the plain images to compute one value of NPCR and UACI.We repeat the process 10 times and the average scores of NPCR and UACI of the corresponding image encryption schemes are recorded in Tables 4 and 5, respectively, as the final NPCR and UACI.
In Table 4, although the NPCR scores by the DFBC are not as good as those by CDCP, they are very close to the maximum theoretical scores.DFBC apparently outperforms HC-DNA in terms of NPCR, and it achieves comparable results with CHC and IC-BSIF.In Table 5, DFBC obviously outperform all the other schemes in 6 out of 8 cases, whereas HC-DNA still gets the poorest results in all cases.The values of NPCR and UACI indicate that the DFBC is able to resist differential attacks.

Robustness Analysis.
A digital image is inevitably contaminated by noise or has data loss during storage and transmission.Image encryption should have the ability to resist both noise and data loss.The experimental results of the DFBC for robustness analysis are shown in Figure 10.We firstly add 0.5%, 1%, and 2% salt and pepper noise to the encrypted Lena, and the corresponding decrypted images are shown in Figures 10(b)-10(d), respectively.It can be seen that, for 0.5% and 1% salt and pepper noise, although the decrypted image contains noise, it can clearly recover the original image very well.However, when the noise increases to 2%, we can only see the outline of Lena.When the encrypted image has 0.4% and 1.56% data loss, the proposed DFBC can recover Lena, as shown in Figures 10(e)-10(f), respectively.For 4% data loss, the decrypted image retains some information for us to recognize Lena, as shown in Figure 10(g).
From the analysis, we can see that the proposed DFBC is robust to a certain extent.However, as pointed out by Hua and Zhou, resisting noise, rotation, and cropping attacks are   current limitations for image encryption by filtering [26].The proposed DFBC will be significantly improved when such limitations are resolved.

Conclusions
Image security is one of the most important branches of information security.To enhance image security, in this paper, we have proposed a novel image encryption algorithm that combines a 7D hyperchaotic system with 5 positive LEs, dynamic filtering operation, and bit cuboid operations, namely, DFBC, for image encryption.Extensive experiments have shown that the proposed DFBC can outperform some state-of-the-art image encryption schemes in terms of several evaluation criteria, indicating that the DFBC is effective for image encryption.In the future, we will study construction of more complex key for the DFBC.Furthermore, since the dynamic filtering can be performed in parallel, we will also implement the algorithm on the framework of CUDA and use GPU to accelerate it.

Figure 3 :
Figure 3: The horizontal and vertical masks are deduced from the original ones with size 2 + 1 to the simplified ones with size  + 1.(a) Horizontal 1D mask.(b) Vertical 1D mask.

Figure 4 :
Figure 4: Demonstration of handling edge of an image.(a) Horizontally expanding when the first column is processing.(b) Vertically expanding when the first row is processing.The mask size is 3.

Figure 8 :
Figure 8: Histograms of the plain images and their corresponding encrypted images.The first column is plain images.The second column is the histograms of the plain images.The third column is diffusion result after 1D dynamic filtering.The fourth column is the histograms of the filtered image in the third column.The fifth column is encrypted images.And the sixth column is the histograms of the encrypted images.

Figure 9 :
Figure 9: The adjacent-pixel distribution maps of the input images and the corresponding encrypted images in horizontal direction.

Figure 10 :
Figure 10: Robustness analysis results of noise and data loss.(a) The encrypted Lena and its decrypted image.(b)-(d) The encrypted Lena with 0.5%, 1%, and 2% salt and pepper noise and its decrypted image, respectively.(e)-(g) The encrypted Lena with 0.4%, 1.56% and 4% data loss and its decrypted image, respectively.

Table 3 :
The correlation coefficients  of the testing images.

Table 4 :
The average NPCR (%) of running the schemes 10 times.
I m a g e D F B C H C -D N A [ 3 0 ] C D C P [ 1 6 ] C H C [ 1 7 ] I C -B S I F [ 2 6 ]

Table 5 :
The average UACI (%) of running the schemes 10 times.