# Novel structures of chaos-based parallel multiple image encryption and FPGA implementation

Thang Manh Hoang, Pham Quang Anh, Manh-Hai Hoang, Do Quang Tran, Safwan El Assad

► To cite this version:

Thang Manh Hoang, Pham Quang Anh, Manh-Hai Hoang, Do Quang Tran, Safwan El Assad. Novel structures of chaos-based parallel multiple image encryption and FPGA implementation. Scientific Reports, 2025, 16 (1), pp.954. 10.1038/s41598-025-30471-x. hal-05449749

HAL Id: hal-05449749

https://hal.science/hal-05449749v1

Submitted on 15 Jan 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# Novel structures of chaos-based parallel multiple image encryption and FPGA implementation

Thang Manh Hoang<fnref n="1" />✉, Pham Quang Anh<fnref n="1" />, Manh-Hai Hoang², Do Quang Tran<fnref n="1" /> & Safwan El Assad³

Image data has been generated massively by devices in medical imaging modalities, cameras, and even by artificial intelligence. Encryption is the powerful method to keep the image content confidential, in which an encryption algorithm must include the confusion and diffusion properties. For massive images, an efficient method of encryption must be chosen to meet the demands of encryption speed and confidentiality. So far, chaos-based image encryption has been an active topic of research because it is considered an effective method to remove the correlation in image data as well as to keep confidential by the involvement of chaotic system in the encryption process. Besides, multiple image encryption algorithms encrypt multiple images in parallel, and it provides highly efficient performance in terms of speed if it is implemented on a parallel computing platform such as multiple core processing as well as digital hardware design. Chaos-based multiple image encryption is constructed by integrating a chaotic system into multiple image encryption. Recently, many algorithms of chaos-based multiple image encryption have been proposed, and they are proved to have high efficiency in terms of both speed and confidentiality. However, all the existing algorithms of chaos-based multiple image encryption require images of the same size and of the same number of bits representing pixels. Further, they encrypt a cohort of plain images at the same time, and all ciphertext images of a cohort must also be decrypted at the same time. It means that if it does not allow to decrypt one or some selected ciphertext images from a cohort separately; and as a result, it wastes time and energy to decrypt unwanted images. In this paper, three novel structures of chaos-based multiple image encryption are proposed which overcome the drawbacks of existing algorithms. That is, the proposed cryptosystems accept cohort images of different sizes; pixels of images can be represented by different numbers of bits; and any selected ciphertext images from a cohort can be decrypted separately. The security is improved by using the session keys of image-content dependency. The proposed structures of multiple image encryption consist of permutation, substitution, and diffusion processes. The difference between three structures is the order of such processes. A perturbed chaotic map and a linear-feedback shift register are employed to generate pseudo-random bit sequences for session keys. The simulation results for the exemplar designs using the proposed structures show the effectiveness by means of the statistical analysis for the session keys using the NIST randomness test, information entropy, histogram, and correlation coefficients of adjacent pixels in ciphertext images, and security analysis by means of space and sensitivity of the secret key. The hardware implementation on the FPGA platform demonstrates the feasibility of the proposed structures by means of throughput and hardware efficiency.

Today, modern devices such as cameras, medical imaging modalities, and even artificial intelligence generate massive images. In order to keep image data confidential, encryption algorithms must be efficient in terms of speed and confidentiality in the context of massive images. So far, chaos-based image encryption is considered an effective method to remove the correlation in image data as well as to keep confidential by the involvement of chaotic system in the encryption process.

Since R. Matthews introduced the first chaos-based encryption<fnref n="1" />, many structures and algorithms have been proposed for chaos-based image encryption. In the early stage of development, a chaotic system was used as a

<footnote n="1">¹School of Electrical and Electronic Engineering, Hanoi University of Science and Technology, 1 Dai Co Viet, 11600 Hanoi, Vietnam. ²Department of Radiology and Diagnostic Imaging, Faculty of Medicine and Dentistry, University of Alberta, Edmonton, AB T5K 0Z4, Canada. ³Nantes University, IETR Laboratory, Nantes 44000, France. ✉email: thang.hoangmanh@hust.edu.vn</footnote>
pseudo-random number generator (PRNG) and chaotic values were mixed with image data to obtain ciphered images, for example<sup>2,3</sup>. In the next stage of development, chaotic systems were utilized in structured cryptography models, i.e., substitution-permutation networks (SPN)<sup>4</sup> and Feistel networks<sup>5</sup>. Recently, chaos has also been used with other techniques to construct chaotic image encryption, e.g., DNA encoding<sup>6,7</sup> and quantum<sup>8,9</sup>. However, most of the existing chaotic encryption have dealt with a single image, namely single image encryption (SIE), e.g.<sup>10,11</sup>.

Rather than SIE, more efficient performance is achieved by multiple image encryption (MIE), in which multiple images, called a cohort, are encrypted simultaneously. Since 2005, MIE has been implemented in the optical domain<sup>12–16</sup>, and in the digital data domain, for example<sup>17,18</sup>. After the first chaos-based MIE was introduced<sup>19</sup>, many works on MIE have been reported, for example<sup>6,20–31</sup>. The input of chaos-based MIE algorithms is multiple images and multiple ciphertext images are achieved at the output. So far, there are several approaches to chaos-based MIE algorithms, and they are categorized based on different purposes, i.e., parallelism, techniques, image-content dependency, and structural configuration. In the following literature review, these approaches are considered only for chaos-based MIE.

*   **Parallelism:** The purpose of designing with parallelism operations is to strive for high speed when it is implemented on parallel computing platforms such as multiple core computing or digital hardware. With respect to the level of parallelism, chaos-based MIE algorithms can be categorized into three groups depending on the processing progress of pixels from images of a cohort being encrypted. The first group of chaos-based MIE algorithms has the lowest level of parallelism, in which all original images are merged to become a large image, then the data of the big image are encrypted pixel by pixel sequentially throughout the big image. For example, the algorithm in<sup>22,24,32–34</sup> is completely sequential, in which all original images are merged, and then permutation and diffusion are performed pixel by pixel through the big image. The encryption scheme in<sup>34</sup> is the kind of lowest level of parallelism, in which four RGB images are merged before permutation and substitution for pixels of the big image. The second group of chaos-based MIE algorithms has the highest level of parallelism, e.g.<sup>7</sup>, in which pixels at the same locations from all original images are encrypted with the same progress. For example, the algorithm in<sup>28,29</sup> is completely parallel if the intra-image permutation scheme is used. The third group of chaos-based MIE algorithms has a moderate level of parallelism, in which some encryption processes of algorithms can be implemented in parallel, and some others are performed sequentially, e.g.<sup>23,31,35,36</sup>. For example, in Karawia<sup>23</sup>, multiple original images with the same size and the same pixel depth are merged to become a big image; then small pieces of images that are divided from the big image are shuffled with each other; after that, small pieces of images are converted into binary vectors, which are scrambled by XORing with binary vectors obtained by chaotic values; scrambled pieces of image are merged to become the big image, before all pixels are shuffled in the range of the big image; the big image is divided into images of the same size of original images. In practice, only the step where small pieces of images are scrambled by XORing with binary vectors generated by the chaotic map can be performed in parallel, while the shuffling operations of small pieces of images and of pixels must be implemented sequentially.

*   **Techniques:** So far, many techniques have been used in the design of chaos-based MIE. Techniques is of variety, ranged from cellular automa (CA), transforms to the deoxyribonucleic acid (DNA), and others. The technique is used with the purpose to ensure the security, that is, to improve with the confusion and diffusion properties for chaos-based MIE. First of all, the cellular automa in<sup>22</sup> is used for generating rules for the permutation process which changes the permutation key from a pixel to pixel by means of their indexes. Secondly, in the DNA technique, bits are mapped by the rules of based pairing<sup>37</sup> by four DNA bases, i.e. Adenine (A), Thymine (T), Guanine (G) and Cytosine (C), and there are eight sets of encoding and decoding rules<sup>38</sup>. Accordingly, the computing operations based on DNA such as XOR, addition, and subtraction, are employed for the permutation and diffusion processes<sup>6,7,22,25,39,40</sup>. After that, the image data in the DNA bases are decoded to represent image in binary. In fact, the DNA technique is a form of binary image encryption, but it deals with groups of two bits. Thirdly, in the technique with transforms, image data is converted into transform domains and the encryption process is performed on the data in the transform domains. Then, the encrypted data is inverse transformed to obtained ciphertext images. Several transforms are used for this purpose, i.e., discrete fractional random transform<sup>20,40,41</sup>, Rubik cube transform<sup>40,42</sup>, discrete fractional cosine transform<sup>43,44</sup>, discrete cosine transform<sup>32,45</sup>, fractional discrete Fourier transform<sup>46,47</sup>, Krawtchouk fractional transform<sup>48</sup>, fractional Hartley transformation<sup>35</sup>, etc. However, some techniques are complicated and slow down the encryption process. For example, the transform requires a lot of computational effort, so it is not suitable for encryption with high-speed priority. The most speed is achieved when the permutation is performed in spatial domain and the diffusion process is carried out with pixel values.

*   **Image-content dependency:** In most existing chaos-based MIE algorithms, the private key and session key are independent from the image content. Recently, there are some schemes of chaos-based MIE that are designed so the the image content is involve in the generation of the keys. There are several two main categories, i.e., secret key of image-content dependency and session keys of image-content dependency. As the first approach, the secret key of image-content dependency is the simplest way, in which some values are computed from the image, e.g. hash values, and used as the secret key. It is is kept constant during encryption. For example, the initial values and/or parameters of chaotic map are hashed values of plaintext images<sup>26,27,33,36,40,42,49,50</sup>. With this way, the secret key at the encryption side must be sent to the decryption side for each plaintext image. This creates the problem in the key management. In the second approach, the session keys are generated and dependent on the image content of intermediate ciphertexts during the encryption. In fact, the orbit of chaotic map is regularly perturbed by pixel values of intermediate ciphertext image while it produces the session keys as presented in<sup>28,40,51–57</sup>. However, most existing algorithms with the session keys of image-content dependency were designed for the chaos-based SIE. To the best our knowledge, there is only one algorithm
presented in<sup>29</sup> being for the chaos-based MIE. In summary, the image encryption with image-content dependency provides superior security, especially for the session keys of image-content dependency, because the encryption is sensitive to changes in the plaintext and ciphertext, in other words, it can resist from the plaintext- and ciphertext-related attacks. In addition, the session keys of image-content dependency requires more computation than the secret key of image-content dependency, so it is only suitable for the encryption algorithms that are implemented in the parallel computing platform.

*   **Structural configuration:** So far, there are two main architectures of chaos-based MIE. The first one is based on the substitution-permutation networks (SPN) and the second one is the Feistel structure<sup>5</sup>. Most of algorithms of chaos-based image employs the SPN, and very few schemes use the Feistel structure<sup>58-60</sup>. For the SPN, the substitution provides the confusion property by a substitution box (S-box) while the diffusion property is achieved by the permutation process. It is well known that the S-box is very simple but it makes the encryption resisted from most typical attacks. Most of existing algorithms of chaos-based encryption, that are successfully analyzed, are of permutation-diffusion structure, and are without any S-box<sup>61-65</sup>. This means that the permutation-diffusion structure does not provide strong security enough, and an S-box should be included in the image encryption algorithm. In addition, one of the weakness in most existing chaos-based MIE algorithms is that they are designed for performing only a single round of encryption, e.g.<sup>7,23,25,49</sup>. It is known that the combination of diffusion and permutation has been proved efficient in saving time for chaotic iterations as presented in<sup>33,66,67</sup> but, there are few chaos-based MIE algorithms designed using the combination of diffusion and permutation<sup>29,33</sup>. Besides, most of existing chaos-based MIE algorithms are with the structures that limit themselves to the number of grayscale images with the same size to be encrypted, i.e., two<sup>19,20,32</sup>, three<sup>7,35,39,42,46,48,68</sup>, six<sup>40,41,44</sup>, twelve<sup>34</sup>, and multiple<sup>24,26,27,33,36,47,49,50</sup>. In case of different sizes, the plaintext images are required to reshape or pad pixels to become the same size before being encrypted. It needs overhead time to reshape and pad additional pixels. It also wastes time and energy to encrypt, and decrypt padded pixels. In practice, there are some areas of applications such as medical imaging modalities<sup>69</sup>, aerial imagery<sup>70</sup>, etc, in which various types of images are of different sizes and different pixel depths. Even the encryption scheme presented in<sup>33</sup> allows multiple images of a cohort with any sizes, it does not accept variety of pixel depths and it encrypts images sequentially, but not in parallel.

In summary, existing chaos-based MIE algorithms still suffer from limited security, performance, and generality. For generality, most existing MIE algorithms lack generality: they rely on specific chaotic maps and particular image parameterizations, employ non-standard architectures (neither SPN nor Feistel) and perform only a single encryption round. In other words, these algorithms are designed for specific scenarios rather than as a general model.

For performance, although parallel MIE algorithms have been proposed, all encrypted images within a cohort must still be decrypted to retrieve only some of them. This implies that undesired images also need to be decrypted to access the desired ones, thereby wasting computational time and energy. Furthermore, all existing parallel MIE algorithms require plaintext images of the same size. Padding pixels with specific values is therefore necessary to make all images uniform in size, which again increases time and energy consumption during encryption and decryption. Moreover, all existing MIE algorithms assume that pixels in every image are represented with the same bit depth. However, in practice certain application domains, such as medical imaging modalities<sup>69</sup> and aerial imagery<sup>70</sup>, involve heterogeneous images of different sizes and bit depths.

Regarding security, most existing MIE algorithms use static session keys that are independent of image content. Ideally, MIE algorithms should employ session keys that depend on the content of the image during the encryption process. This dependency makes the schemes more resistant to attacks that exploit relationships between plaintext and ciphertext, such as chosen-plaintext, known- plaintext, and chosen-ciphertext attacks. In addition, many chaotic image encryption schemes have been successfully cryptanalyzed, as reported in<sup>61-65</sup>. The main reason is that their diffusion processes are based on simple linear transformations between plaintext and previous ciphertext pixels. Although such linear equations provide high processing speed, they do not guarantee adequate security. Therefore, nonlinear substitution mechanisms, such as S-boxes, should be incorporated to enhance robustness.

These gaps in the existing MIE algorithms are addressed in our work. In this work, three MIE structures are proposed to overcome the disadvantages, and the main contributions and novelty are as follows.

*   The proposed structures of chaos-based MIE are generic models in the form of substitution-permutation network, so that allow for employing any types of chaotic map, applying any number of encryption rounds to meet the security requirement. Three processes, i.e., permutation, diffusion, and substitution (S-box), are included in the proposed models, and those are combined in the same step and are performed pixel by pixel within individual images. The proposed models can encrypt any number of images, and work on spatial and pixel value domains. The images of a cohort are encrypted in parallel, while any ciphertext images from a cohort can be decrypted to retrieve individual images. Images can be of different sizes and different pixel depths.

*   The session keys in the proposed structures are dependent on the pixel values of the intermediate ciphertext images. The chaotic map is frequently perturbed by a random number generator, and chaotic values, together with random numbers, are used to generate session keys in encryption. In addition, the substitution box (S-box) is included to enhance the confusion property of the encryption. As a result, ciphers that employ the proposed structures can withstand any type of attack that exploits the difference in the plaintexts and ciphertexts.

*   The digital hardware is designed and implemented on the FPGA platform to demonstrate the effectiveness and feasibility of the proposed structures in terms of throughput and efficiency. The resultant performance is compared to other existing cryptographic hardware designs.
The remainder of the paper is organized as follows. Section Proposed MIE structures presents the details of the encryption, decryption, and operation of the proposed structures. Next, section Exemplar simulation describes the details about the cohort of images, the chaotic map, and the values of the simulation parameters. The simulation results and statistical and security analysis are discussed in this section. Then, section Digital hardware design presents details of the hardware architecture for Structure 1, including operation, and architecture of the MIE Intellectual Property and FPGA Implementation results. The throughput and efficiency of the hardware design are calculated and compared with those in other reported designs. Finally, sections Discussion and Conclusion highlight the highlights of the proposed structures in terms of operation, security, and hardware performance.

## Proposed MIE structures

Let us define $\mathbb{I} = \{I_k | k = 1, ..., K\}$ to be a set of K images, called a cohort, in which image $I_k = \{p(i_k, j_k) | i_k = 1, ..., M_k, j_k = 1, ..., N_k\}$ has a size of $M_k \times N_k$; $p(i_k, j_k)$ denotes the pixel value of image k at coordinates $(i_k, j_k)$, and is represented by $Q_k$ bits, called the depth of pixels. It is assumed that the images of a cohort can be of different sizes ($M_i \neq M_j$ and/or $N_i \neq N_j$ for $i \neq j$) and that the pixels of different images in a cohort have different pixel depths ($Q_i \neq Q_j$ for $i \neq j$). $I_k$ is simply considered as a matrix of values.

## Encryption

As shown in Fig. 1, each of the proposed structures consists of three processes, that is, permutation, diffusion, and substitution. Here, the diffusion property is achieved by the permutation and diffusion processes, while the confusion property is obtained by the substitution process. The diffusion process depends on the session key; therefore, it also enhances the confusion property. In fact, there can be several encryption structures depending on the relative position of these processes. In this research, the permutation is fixed and performed first, and then the others to ensure that the values of all pixels are changed by the encryption. These processes are performed separately and in parallel on all images. The encryption procedure, details on processes, and the function of blocks of the proposed structures are presented in the following.

Consider encryption for K images at the same time. The plaintext image k with the size $M_k \times N_k$, $k \in \{1, 2, ..., K\}$, is denoted by $I_k = \{p(i_k, j_k) | 1 \le i_k \le M_k, 1 \le j_k \le M_k\}$. The source and destination locations are respectively denoted by $XY_k = (i_k, j_k)$ and $XY_k = (i'_k, j'_k)$, where i and j are row numbers, and k and j' are column numbers of the image k, with $k \in \{1, 2, ..., K\}$. Intermediate-ciphertext pixels at the source and destination locations are denoted with a hat as $\hat{c}(i_k, j_k)$ and $\hat{c}(i'_k, j'_k)$, respectively. The ciphertext pixels in an encryption round at the source and destination locations are denoted by $c(i_k, j_k)$ and $c(i'_k, j'_k)$, respectively. The final ciphertext image k after $N_e$ rounds of encryption is denoted by $C_k$, which is constructed by ciphertext pixels $c(i_k, j_k)$ with $1 \le i_k \le M_k$ and $1 \le j_k \le M_k$.

In the encryption procedure, pixels are scanned from left to right, row by row until all pixels in an image are encrypted. All images of a cohort are encrypted separately and in parallel for $N_e$ rounds of encryption.

### Structure 1

The first proposed structure, named Structure 1, is illustrated in Fig. 1a, in which the order of the processes is permutation, substitution, and then diffusion.

*   **Permutation:** The permutation is to move the plain source pixel $p(i_k, j_k)$ and the destination pixel $p(i'_k, j'_k)$ in the same image k. Here, $i_k$ and $j_k$ of the source pixel are incrementally scanned from left to right and row by row of the image $I_k$, and $i'_k$ and $j'_k$ of the destination pixel are induced from chaotic values. Details about the way in which the values of $i_k$ and $j'_k$ are calculated will be shown in the next subsection.

*   **Substitution:** The values of the permuted pixels at the source and destination locations, $p(i_k, j_k)$ and $p(i'_k, j'_k)$, are then changed by substitution. The substitution is carried out by the S-box to get substituted values as

$$ \begin{aligned} \hat{c}(i_k, j_k) &= \text{S-box}(p(i_k, j_k)), \\ \hat{c}(i'_k, j'_k) &= \text{S-box}(p(i'_k, j'_k)), \end{aligned} \quad (1) $$

where S-box(.) is the function that represents the substitution process and returns the substitution values.

*   **Diffusion:** After being replaced by Eq. (1), the pixels at the source and destination locations, $\hat{c}(i_k, j_k)$ and $\hat{c}(i'_k, j'_k)$, are diffused to obtain ciphertext pixels as

$$ \begin{aligned} c(i_k, j_k) &= F_k^{(Diff)} \left( \hat{c}(i_k, j_k), c(i_k, j_k - 1), \phi_k^{(S)} \right), \\ c(i'_k, j'_k) &= F_k^{(Diff)} \left( \hat{c}(i'_k, j'_k), c(i'_k, j'_k - 1), \phi_k^{(D)} \right), \end{aligned} \quad (2) $$

where $F_k^{(Diff)}(.)$ is the diffusion function that operates bit-wise; $c(i_k, j_k - 1)$ and $c(i'_k, j'_k - 1)$ are the values of the ciphertext pixels that stand immediately before the source and destination locations, respectively; and $\phi_k^{(S)}$ and $\phi_k^{(D)}$ are pseudo-random values for the source and destination pixels. Note that the superscripts (S) and (D) are denoted for the source and destination in the following text.

In Structure 1, every source pixel is subsequently performed by three processes as described above. It is clear that the processes of permutation, substitution, and diffusion are combined and performed pixel-by-pixel. All images of a cohort are encrypted separately and in parallel for $N_e$ encryption rounds.
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** The detail of proposed structures.</figcaption></figure>

### Structure 2

The second proposed structure, named Structure 2, is presented in Fig. 1b. After permuting pixels, the source and destination pixels are diffused and then substituted as specified below.

* Permutation: The permutation of Structure 2 is the same as that in Structure 1.

* Diffusion: Pixels at the source and destination locations $p(i_k, j_k)$ and $p(i'_k, j'_k)$ are diffused by
$$ 
\begin{aligned}
\hat{c}(i_k, j_k) &= F_k^{(\text{Diff})} \left( p(i_k, j_k), c(i_k, j_k - 1), \phi_k^{(S)} \right), \\
\hat{c}(i'_k, j'_k) &= F_k^{(\text{Diff})} \left( p(i'_k, j'_k), c(i'_k, j'_k - 1), \phi_k^{(D)} \right),
\end{aligned} 
\quad (3) $$

where function $F_k^{(\text{Diff})}(.)$, pixels $c(i_k, j_k - 1)$, $c(i'_k, j'_k - 1)$, and chaotic values $\phi_k^{(S)}$ and $\phi_k^{(D)}$ are defined in Eq. (2).

* Substitution: The diffused pixels at the source and destination locations, $\hat{c}(i_k, j_k)$ and $\hat{c}(i'_k, j'_k)$, are then changed by substitution. The substitution is carried out by the S-box to get ciphertext pixels as

$$ 
\begin{aligned}
c(i_k, j_k) &= \text{S-box}(\hat{c}(i_k, j_k)), \\
c(i'_k, j'_k) &= \text{S-box}(\hat{c}(i'_k, j'_k)).
\end{aligned} 
\quad (4) $$

Compared with Structure 1, the diffusion process of Structure 2 is carried out before substitution. The encryption procedure is also performed pixel by pixel.

### Structure 3

The third proposed structure, named Structure 3, is illustrated in Fig. 1c, in which permutation and diffusion are combined and performed pixel by pixel, the same as those in Structure 2. After all pixels of an image are permuted and diffused, a substitution is then performed on all pixels separated from permutation and diffusion. Specifically, the encryption processes are as follows.

* Permutation: The permutation of Structure 3 is the same as that in Structure 1.

* Diffusion: The permutation of Structure 3 is the same as that in Structure 2 as given in Eq. (3). After all pixels of the image are permuted and diffused, the intermediate image $\hat{C}_k = \{\hat{c}(i_k, j_k) | 1 \le i_k \le N_k, 1 \le j_k \le M_k\}$ is obtained.

* Substitution: The pixels in the intermediate image $\hat{C}_k$ are substituted by the S-box as

$$ c(i_k, j_k) = \text{S-box}(\hat{c}(i_k, j_k)). \quad (5) $$

where $\hat{c}(i_k, j_k)$ is the pixel in location $XY_k = (i_k, j_k)$ of $\hat{C}_k$.

### Remarks

The highlight of the proposed structures can be derived from the configuration as follows.

* More than one process is combined and performed at the same time. Those processes are permutation, substitution, and diffusion in Structure 1; permutation, diffusion, and substitution in Structure 2; and permutation and diffusion in Structure 3. The combination of more than one process is the difference in the proposed structures from most existing encryption algorithms. As a result, chaotic values are efficiently used for the processes of permutation and diffusion, and the speed is improved compared to that of most existing encryption algorithms in which all pixels of an image are completed for a process before the next process is performed, where it requires chaotic iterations for separate processes.

* As illustrated in Fig. 1, the proposed structures encrypt images separately and in parallel. Therefore, the various sizes of the images can be accepted. In addition, because the session keys are generated separately for multiple images, the decryption can be done separately for any selected images. In addition, encryption separation allows separation of the decryption, which is the unique feature of the proposed structures.

### Decryption

The decryption structures have the same block functions except that the order of processes is reversed compared to that in the encryption. Specifically, the permutation is the same as in the encryption while inverse substitution and inverse diffusion are used. In addition, in every individual image, the direction of the pixel scanning must be reversed in compared with that in the encryption. The decryption is also performed in $N_e$ rounds.

All blocks of the proposed structures that provide parameters for the encryption processes are described in detail below.

### Functions of blocks in the proposed structures

A chaotic map is used to generate pseudo-random values, and its dynamics of a chaotic map is perturbed by bits from the block Session Key Scheduling (SKS), the so-called perturbed chaotic map (PCM). The PCM and SKS are initiated and then iterated to generate bit sequences for permutation and diffusion during encryption. As shown in Fig. 1, Bit Pre-processing of image k, $BPP_k$, uses bits of next plaintext, $p(i_k, j_k + 1)$, previous ciphered pixels, $c(i_k, j_k - 1)$, and the SKS bit sequence, $sk_k$, to generate the coordinates of the destination pixel for permutation ($XY_k$) and pseudo-random values ($\Phi_k$) for diffusion to encrypt the present pixel of image K. The functions of blocks in the proposed structures are explained as follows.

#### Perturbed chaotic map (PCM)

The general model of PCM is employed in the proposed structures expressed by
$$ 
\left\{ 
\begin{aligned}
X_{n+1} &= F(X_n, \hat{\Gamma}_n), \\
\hat{\Gamma}_n &= \Gamma_0 \oplus \Delta\Gamma_n, \\
\Delta\Gamma_n &= \left\{ \begin{array}{ll} Y_1(E) & \text{for } n=1, \\ Y_2(X_n) & \text{for } n>1, \end{array} \right.
\end{aligned} 
\right. 
\qquad (6) 
$$

where $F(.)$ is a chaotic function; $Y_i$ with $i = 1, ..., 4$ are bit arrangements in the construction of bit sequences of perturbation; $X_n = [x_n^{(D)} \ x_n^{(D-1)} \ ... \ x_n^{(2)} \ x_n^{(1)}]^T$ is the vector of state variables; $\Gamma_0 = [\gamma_0^{(G)} \ \gamma_0^{(G-1)} \ ... \ \gamma_0^{(2)} \ \gamma_0^{(1)}]^T$ is the vector of initial control parameters and is regularly perturbed by $\Delta\Gamma_n = [\delta_{\gamma_n^{(G)}} \ \delta_{\gamma_n^{(G-1)}} \ ... \ \delta_{\gamma_n^{(2)}} \ \delta_{\gamma_n^{(1)}}]^T$. $D$ and $G$ are the numbers of dimensions and those of the control parameters of the chaotic map, respectively. $\hat{\Gamma}_n = [\hat{\gamma}_n^{(G)} \ \hat{\gamma}_n^{(G-1)} \ ... \ \hat{\gamma}_n^{(2)} \ \hat{\gamma}_n^{(1)}]^T$ is the vector of perturbed control parameters$^{53}$.

Figure 2 illustrates the general PCM configuration, in which $Y_i$, $i = 1, 2$, performs the bit arrangement so that the bit of $X_n$ is extracted and used to perturb the chaotic map. The bit sequence $E$ with $L_E = L_X + L_\Gamma$ bits in length is used as an external perturbation. The values of $x_n^{(d)}$ and $\gamma_n^{(g)}$, with $d = 1, ..., D$ and $g = 1, ..., G$, are represented by the $m_1$ and $m_2$ bits, respectively. The values of $x_n^{(d)}$ and $\gamma_n^{(g)}$ are real, so the numbers of bits for the integer and fraction portions of $x_n^{(d)}$ and $\gamma_n^{(g)}$ are denoted by $\langle m_1^{(\text{int})}, m_1^{(\text{frac})} \rangle$ and $\langle m_2^{(\text{int})}, m_2^{(\text{frac})} \rangle$, respectively, with $m_1 = m_1^{(\text{int})} + m_1^{(\text{frac})}$ and $m_2 = m_2^{(\text{int})} + m_2^{(\text{frac})}$ as

$$ 
\begin{aligned}
x_n^{(d)} &= \sum_{i=-m_1^{(\text{frac})}-1}^{m_1^{(\text{int})}-1} a_i \cdot 2^i, \\
\gamma_n^{(g)} &= \sum_{i=-m_2^{(\text{frac})}-1}^{m_2^{(\text{int})}-1} b_i \cdot 2^i.
\end{aligned} 
\qquad (7) 
$$

where $a_i$ and $b_i$ are binary. The amounts of perturbation, $\delta_{x_n^{(d)}}$ and $\delta_{\gamma_n^{(g)}}$, are

$$ 
\begin{aligned}
\delta_{x_n^{(d)}} &= \sum_{i=-m_1^{(\text{frac})}-1}^{m_1^{(\text{int})}-1} c_i \cdot 2^i, \\
\delta_{\gamma_n^{(g)}} &= \sum_{i=-m_2^{(\text{frac})}-1}^{m_2^{(\text{int})}-1} d_i \cdot 2^i.
\end{aligned} 
\qquad (8) 
$$

where $c_i$ and $d_i$ are binary. Similarly, the values of $\hat{X}_n$ and $\hat{\Gamma}_n$ are represented by

$$ 
\begin{aligned}
\hat{x}_n^{(d)} &= \sum_{i=-m_1^{(\text{frac})}-1}^{m_1^{(\text{int})}-1} e_i \cdot 2^i \quad \text{for } j = 1, \dots, D \\
\hat{\gamma}_n^{(g)} &= \sum_{i=-m_2^{(\text{frac})}-1}^{m_2^{(\text{int})}-1} f_i \cdot 2^i \quad \text{for } k = 1, \dots, G.
\end{aligned} 
\qquad (9) 
$$

where $e_i$ and $f_i$ are also binary. The perturbation is the XORing operation as follows.

$$ 
\begin{aligned}
e_i &= a_i \oplus c_i \quad \text{for } i = -m_1^{(\text{frac})}, \dots, m_1^{(\text{int})} - 1, \\
f_i &= b_i \oplus d_i \quad \text{for } i = -m_2^{(\text{frac})}, \dots, m_2^{(\text{int})} - 1.
\end{aligned} 
\qquad (10) 
$$

The value ranges of state variables and control parameters must be constrained in those as in the original model. So, a certain number of bits in the values of $x_n^{(d)}$ and $\gamma_n^{(g)}$ are fixed with specified values to ensure that the PCM exhibits chaotic behavior, while the other bits can be flippable. The positions of the perturbed bits are specified by $Y_i$.

In this work, PCM is initiated by $X_0$ and $\Gamma_0$, and iterates for $R$ times to get $X_R$ to encrypt a set of $K$ pixels; each pixel comes from an image. This process is repeated until all pixels in every image are encrypted for $N_e$ rounds of encryption. More details on PCM can be found in$^{53}$.

<figure><img src="image_4.png" /><figcaption>Fig. 2. The configuration of the PCM$^{53}$.</figcaption></figure>
### Session key scheduling (SKS)

The SKS produces session keys to encrypt individual pixels and perturb the PCM. The core of SKS is the keystream generator based on the linear feedback shift register (LFSR)⁷¹.

The input of the keystream generator is the initial vector $S_0$ and the predefined private key ID. The SKS outputs are the session keys $sk_k$ with $k = 1, ..., K$ and the bit sequence $E$ for PCM perturbation. The bit sequence $sk_k$ will be used to identify the location of destination pixels ($XY'_k$) and to generate diffusion values ($\Phi_k$), so the length of $sk_k$ must be

$$L_{sk} = \lceil \log_2(N_{pix\_max}) \rceil + 2Q_{max} \quad (\text{bits}) \qquad (11)$$

where $N_{pix\_max} = M_{max} * N_{max}$ is with $M_{max} = \max_{k \in \{1, 2, ..., K\}} M_k$ and $N_{max} = \max_{k \in \{1, 2, ..., K\}} N_k$; $Q_{max}$ is the maximum number of bits that represent a pixel of images, $Q_{max} = \max_{k \in \{1, 2, ..., K\}} Q_k$; and $\lceil \cdot \rceil$ is the round-up function.

The length $L_E$ requires ($K * L_{sk} + L_E$) bits to encrypt a group of K pixels from K images at the same time. Therefore, the sequence requires at least ($K * L_{sk} + L_E$) $N_{pix\_max}$ bits to encrypt all pixels of K images. Consequently, the number of states of LFSR⁷² must be at least

$$N_{LFSR} = \lceil \log_2 ((K * L_{sk} + L_E) * N_{pix\_max}) \rceil \quad (\text{states}). \qquad (12)$$

As shown in Fig. 3a, each frame is of $K * L_{sk} + L_E$ bits to encrypt K pixels from K images. So, it requires $N_{pix\_max}$ frames to encrypt all pixels of K images. In detail, each sequence $sk_k$ is used to produce pseudo-random values, as shown in Fig. 3b.

The bit sequence $E$ is used to construct $\delta_{x_k}^{(d)}$ and $\delta_{y_k}^{(g)}$ to perturb the PCM as shown in Eq. (8). The length of $E$, $L_E$, is considered using the numbers of bits represented for the values of the state variables, and the control parameters are $D * m_1$ and $G * m_2$, respectively. However, some constraints are made for the range of state variables and control parameters of the PCM to ensure the existence of chaotic behavior. In terms of security, the number of bits that are extracted to construct $XY'_k$ and $\Phi_k$ should be lower than $D * m_1$, so the attacker is challenged to recover chaotic dynamics and infer values of $X_0$ and $\Gamma_0$.

### Bit Pre-processing (BPP)

The BPP takes bits from $X_R$, $sk_k$, and $pc_k$ to produce the coordinate of the destination pixel, $XY'_k = (i'_k, j'_k)$, for permutation with the pixel at $(i_k, j_k)$ of the image k as

$$ \begin{aligned} i'_{k} &= \text{mod}(F_{k}^{(IJ)}(bs_{k}^{(I)}, sk_{k}^{(I)}, pc_{k}), M_{k}) + 1, \\ j_{k} &= \text{mod}(F_{k}^{(IJ)}(bs_{k}^{(J)}, sk_{k}^{(J)}, pc_{k}), N_{k}) + 1, \end{aligned} \qquad (13) $$

and pseudo-random values, $\Phi_k = \{\phi_k^{(S)}, \phi_k^{(D)}\}$, for diffusion of source and destination pixels as

$$ \begin{aligned} \phi_k^{(S)} &= F_k^{(\Phi)}(bs_k^{(S)}, sk_k^{(S)}, pc_k), \\ \phi_k^{(D)} &= F_k^{(\Phi)}(bs_k^{(D)}, sk_k^{(D)}, pc_k), \end{aligned} \qquad (14) $$

where $F_k^{(IJ)}(.)$ and $F_k^{(\Phi)}(.)$ are bitwise functions of image k. bit sequences $bs_k^{(I)}, bs_k^{(J)}, bs_k^{(S)}$, and $bs_k^{(D)}$ are extracted from the values of $X_R$; bit sequence $pc_k$ is formed by interleaving bits of pixels $p(i_k, j_k + 1)$ and $c(i_k, j_k - 1)$ as given in Eq. (16); $M_k$ and $N_k$ are the actual size of the image $I_k$, and each pixel is represented by $Q_k$ bits; and $\text{mod}(.)$ is the modulus function.

<figure><img src="image_5.png" /><figcaption>Fig. 3. Arrangement of the keystream in the SKS.</figcaption></figure>
Bit interleaving

The block Bit interleaving is used to interleave bits of the next plain pixel $p(i_k, j_k + 1)$ and the past ciphertext pixel $c(i_k, j_k - 1)$ of the image $I_k$, and a bit sequence $pc_k$ is obtained. Specifically, the representation of $p(i_k, j_k + 1)$ and $c(i_k, j_k - 1)$ is

$$ 
\begin{aligned}
p(i_k, j_k + 1) &= \sum_{i=0}^{q_k-1} b_i^{(p)} \cdot 2^i, \\
c(i_k, j_k - 1) &= \sum_{i=0}^{q_k-1} b_i^{(c)} \cdot 2^i,
\end{aligned} 
\qquad (15) 
$$

where $b_i^{(p)}$ and $b_i^{(c)}$ are binary. The bitwise operation of interleaving is as follows.

$$ pc_k = b_{q_k-1}^{(p)} b_{q_k-1}^{(c)} b_{q_k-2}^{(p)} b_{q_k-2}^{(c)} \cdots b_0^{(p)} b_0^{(c)} = \left\|_{i=0}^{q_k-1} b_i^{(p)} b_i^{(c)} \right\|, \qquad (16) $$

where $\|$ is the bit concatenation operation. The initial ciphered pixels $c_k^{(0)}$ and the plain pixels $p(M_k, N_k + 1)$ denoted by $p_k^{(0)}$ for $\forall k \in \{1, 2, \dots, K\}$ are used to encrypt the pixels at (1, 1) and the last pixel at $(M_k, N_k)$, respectively. Those are part of the secret key.

Operation and security of the proposed structures

In this subsection, the operation of the structures proposed in Fig. 1 is described step by step in the form of a flow chart as in Fig. 4. The security in terms of the secret key space and the image-content dependency is discussed.

The operation

The flow charts in Fig. 4 show that the encryption is separated for each image and in parallel for multiple images. The inputs are the parameters' values, including the number of encryption rounds ($N_e$), the initial vectors and control parameters of the chaotic map ($X_0$ and $\Gamma_0$), the number of cohort images ($K$) and the size of images ($M_k, N_k$ and $Q_k$), the number of chaotic iterations ($R$), the substitution box (S-box), initial secret values ($p_k^{(0)}$, $c_k^{(0)}$). From that, running values are set for the encryption of $K$ images.

Each encryption session is separate for an individual image. For each image encryption session, $Pause_k$ is the flag that indicates that the encryption process for image $I_k$ is in progress ($Pause_k = 0$) or complete ($Pause_k = 1$). In an encryption progress for the image $I_k$, the processes of permutation, diffusion, and substitution are performed in the orders described in Subsection Description for the encryption.

<figure><img src="image_6.png" /><figcaption>Fig. 4. The flowchart of the encryption.</figcaption></figure>
As shown in the flowchart for Structures 1 and 2 in Fig. 4a and 4b, pixels are chosen to encrypt sequentially, that is, scanning from left to right and row by row. Each pair of pixels, source and destination pixels, in Structures 1 and 2 is processed by combining the permutation, diffusion, and substitution in one step. Then, the next pair of pixels is processed as indicated by the increment of $i_k$ and $j_k$. During each process, appropriate equations as described for Structures 1 and 2 in Subsection Description for the encryption are used. If all pixels of $I_k$ are scanned and encrypted for rounds $N_e$, that is, $n_k > N_e$, the encryption of the image $I_k$ is complete by setting the flag $Pause_k = 1$.

The flowchart in Fig. 4c for Structure 3 shows that the permutation and diffusion processes are carried out the same as those in Structures 1 and 2, except that the substitution process is separated from the permutation and diffusion. The substitution process is performed after all pixels of $I_k$ are permuted and diffused.

The encryption for all images is finished if all the flags are set to '1', that is, $Pause_k = 1$ for $\forall k \in \{1, 2, \dots, K\}$. Specifically, K encryption sessions are performed in parallel for K images. The chaotic map is iterated while the LFSR is shifted at the same time. The session keys are generated in parallel for K encryption sessions.

It is clear from the structural configuration and flow charts of the encryption algorithms that all images are encrypted separately and in parallel. Therefore, three features are achieved; that is, it allows images with different sizes and different pixel depths, and it allows decrypting some ciphertext images of a cohort, instead of decryption for all ciphertext images as in existing algorithms. In addition, these structures can be easily implemented on the parallel computing platform, such as multiple cores and hardware design, to achieve high-speed encryption.

### The security

Consider the security of the proposed structures in terms of secret key space and image-content dependency. Firstly, the secret key space of the proposed structures is considered by several factors, i.e. the initial vectors of state variables and control parameters of the chaotic map, the LFSR's seed, initial values of ciphertexts $c_k^{(0)}$ and plaintexts $p_k^{(0)}$ with $k = 1, 2, \dots, K$ as in Eq. (15), and the user key of the block SKS. Specifically, the number of bits for the LFSR's seed and the user key can be chosen as desired, but it must satisfy the number of states given in Eq. (12) to ensure randomness of session keys. The initial values of $c_k^{(0)}$ and $p_k^{(0)}$ with $k = 1, 2, \dots, K$ have a fixed number of bits. For example, 64 bits for $c_k^{(0)}$ and 64 bits for $p_k^{(0)}$ with $K = 8$ and 8-bit grayscale pixels. The initial vectors of the chaotic map consist of the initial value of the state variables, $X_0$, and the initial value of the control parameters, $\Gamma_0$. Here, the chaotic map is implemented on a digital platform, so $X_0$ and $\Gamma_0$ are represented by a number of bits. The chaotic map is perturbed as shown in Fig. 2, therefore, the bits that are changed during the iteration of the chaotic map are considered the secret key. The number of bits of the secret key can be as long as desired and depending on the representation for the values of state variables and control parameters. In addition, a chosen chaotic map whose state variables and control parameters have a large value range of chaotic behavior also provides a large number of bits for the secret key. Moreover, depending on the chaotic map adopted, higher dimensions of state space can provide a large number of bits for the secret key. For example, the Logistic map has only one state variable, and the value range of chaotic variable is in the range 0 to 1. Therefore, the number of bits for the secret key is dependent on the number of bits representing the values of the state variable and control parameters. In contrast, the Cat map and the Standard map are two-dimensional state variables, and the value ranges of state variables and control parameters that exhibit chaos are larger than the Logistic map. Therefore, the Cat map and Standard maps can provide a number of bits for the secret key larger than that from the Logistic map. In summary, the number of bits for the secret key or the key space can be as large as desired in several ways, as given above. However, adoptions of a large number of bits representing values of state variables and control parameters and/or higher-dimensional chaos pay the speed of encryption. It takes longer to compute for an iteration of the chaotic map.

Secondly, the proposed MIE algorithms have the image-content dependency. Specifically, the current session keys are the coordinates of destination pixels $XY_k$ and pseudo-random values $\Phi_k$ with $k = 1, 2, \dots, K$. For image $I_k$, $XY_k$ and $\Phi_k$ in Eqs. (13) and (14) are directly related to the values of the previous ciphertext pixel $c(i_k, j_k - 1)$ and the next plaintext pixel $p(i_k, j_k + 1)$ by means of $p_{C_k}$ as given in Eq. (16). This means that the permutation and diffusion processes depend on the image content. The dependency on image content provides better security and resistance to the types of attacks such as chosen-plaintext and chosen-ciphertext attacks. However, the algorithms of proposed structures continually encrypt pixel-by-pixel throughout images, so they cannot resist the ciphertext-crop attack.

### Exemplar simulation

In this section, examples of all proposed structures are simulated using Python v3.12.6, and then statistical and security analyzes are performed using Matlab R2016a for the simulation results to confirm effectiveness.

### Chosen chaotic maps, images and values of parameters for simulation

In this example, a set of eight images as shown in Fig. 5 is encrypted using each of the three proposed designs. Table 1 shows the size of the image and the number of bits that represent pixels in different images. As given in Table 2, the critical values for encryption are to specify its capability to adopt the values of other parameters. $K_{max}$, $M_{max}$, $N_{max}$, and $Q_{max}$ are, respectively, the maximum numbers of images, rows, and columns of pixels, and of bits representing a pixel. Here, the chosen values are greater than or at least equal to the actual numbers characterized by the images.

* Chosen chaotic maps: The Cat map is chosen and its perturbation is described in the form of Eq. (6) as
<figure><img src="image_7.png" /><figcaption>(a)
Clock</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(b) Flowers</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(c) Baboon</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(d) Peppers</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(e) Boats</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(f) Goldhill</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>(g)
Black</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>(h)
White</figcaption></figure>

<table><caption>**Fig. 5.** Exemplar set of images with actual ratios.</caption>
  <thead>
    <tr>
      <th>Image</th>
      <th>Color</th>
      <th>Size (M<sub>k</sub> × N<sub>k</sub>)</th>
      <th>#bits per pixel (Q<sub>k</sub>)</th>
    </tr>
  </thead>
  <tbody>
    <tr>
      <td>Clock</td>
      <td>GrayScale</td>
      <td>256 × 256</td>
      <td>8</td>
    </tr>
    <tr>
      <td>Flowers</td>
      <td>RGB</td>
      <td>362 × 500</td>
      <td>24</td>
    </tr>
    <tr>
      <td>Baboon</td>
      <td>RGB</td>
      <td>480 × 500</td>
      <td>24</td>
    </tr>
    <tr>
      <td>Peppers</td>
      <td>RGB</td>
      <td>512 × 512</td>
      <td>24</td>
    </tr>
    <tr>
      <td>Boats</td>
      <td>GrayScale</td>
      <td>576 × 720</td>
      <td>8</td>
    </tr>
    <tr>
      <td>Goldhill</td>
      <td>RGB</td>
      <td>532 × 720</td>
      <td>24</td>
    </tr>
    <tr>
      <td>Black</td>
      <td>GrayScale</td>
      <td>790 × 616</td>
      <td>8</td>
    </tr>
    <tr>
      <td>White</td>
      <td>GrayScale</td>
      <td>492 × 874</td>
      <td>8</td>
    </tr>
  </tbody>
</table>

<table><caption>**Table 1.** Parameters of chosen images.</caption>
  <thead>
    <tr>
      <th>Highest parameter</th>
      <th>Value</th>
      <th>Unit</th>
    </tr>
  </thead>
  <tbody>
    <tr>
      <td>K<sub>max</sub></td>
      <td>16</td>
      <td>images</td>
    </tr>
    <tr>
      <td>M<sub>max</sub></td>
      <td>1024</td>
      <td>rows</td>
    </tr>
    <tr>
      <td>N<sub>max</sub></td>
      <td>1024</td>
      <td>columns</td>
    </tr>
    <tr>
      <td>N<sub>pix_max</sub></td>
      <td>2<sup>10</sup> × 2<sup>10</sup></td>
      <td>pixels</td>
    </tr>
    <tr>
      <td>Q<sub>max</sub></td>
      <td>24</td>
      <td>bits</td>
    </tr>
  </tbody>
</table>

<caption>**Table 2.** Critical values for the exemplar encryption.</caption>

$$
\left\{
\begin{array}{l}
\begin{bmatrix} x_{n+1}^{(2)} \\ x_{n+1}^{(1)} \end{bmatrix} = \text{mod} \left( \left[ \hat{\gamma}_n^{(1)} x_n^{(2)} + (1 + \hat{\gamma}_n^{(2)} \hat{\gamma}_n^{(1)}) x_n^{(1)} \right], 1 \right), \\
\\
\begin{bmatrix} \hat{\gamma}_n^{(2)} \\ \hat{\gamma}_n^{(1)} \end{bmatrix} = \begin{bmatrix} \gamma_0^{(2)} \\ \gamma_0^{(1)} \end{bmatrix} \oplus \begin{bmatrix} \delta \gamma_n^{(2)} \\ \delta \gamma_n^{(1)} \end{bmatrix}, \\
\\
\begin{bmatrix} \delta \gamma_n^{(2)} \\ \delta \gamma_n^{(1)} \end{bmatrix} = \left\{
\begin{array}{ll}
\begin{bmatrix} Y_1^{(a)}(E) \\ Y_1^{(b)}(E) \end{bmatrix} & \text{for } n=1, \\
\\
\begin{bmatrix} Y_2^{(a)}(x_n^{(2)}) \\ Y_2^{(b)}(x_n^{(1)}) \end{bmatrix} & \text{for } 1 < n \le R,
\end{array}
\right.
\end{array}
\right.
\tag{17}
$$

where $X_n = [x_n^{(2)} \ x_n^{(1)}]^T$ are state variables and $\Delta_{\Gamma_n} = [\gamma_n^{(2)} \ \gamma_n^{(1)}]^T$ are control parameters; these are in the ranges of $x_n^{(2)}, x_n^{(1)} \in (0, 1)$ and $\gamma_n^{(2)}, \gamma_n^{(1)} > 0$; $Y_2^{(a)}, Y_2^{(b)}, Y_1^{(a)}$ and $Y_1^{(b)}$ are to arrange bits of either E or $X_n$ to construct permutation amounts. The values of the state variables and control parameters are represented as fixed-point numbers. The formats, bit patterns, value ranges, and bit positions to be perturbed are given in Table 3. $Y_1^{(a)}$ and $Y_1^{(b)}$ are chosen the same format as $Y_2^{(a)}$ and $Y_2^{(b)}$, respectively, but bits are from E. Bit arrangements $Y_2^{(a)}$ and $Y_2^{(b)}$ are chosen to construct perturbation amounts $\Delta_{\Gamma_n}$ as shown in Table 4; where $B_0$ at specific positions are fixed at '0', and the group $(i, b_j)$ is denoted for bit $b_j$ of the state variable $x_n^{(i)}$. These are bits used as $d_i$ in Eq.(8), and the perturbation is based on the XOR operation as given in Eq. (10). The initial values in decimal for the PCM are chosen as in Table 5. It is noted that even though the value range of control parameters are separated and the perturbation is applied, the chaotic behavior certainly exists with nonstationary dynamics.
<table><caption>**Table 3.** Representation of state variables and control parameters, their value ranges, and initially perturbed bits.</caption><thead><tr><th>Parameter</th><th>Format</th><th>Pattern in binary</th><th>Value range</th><th>Position of 16 perturbed bits</th></tr></thead><tbody><tr><td>γ<sub>n</sub><sup>(2)</sup></td><td>⟨52, 49⟩</td><td>xx1(.)xxx...xx</td><td>[1.0,2.0), [3.0,4.0), [5.0,6.0), and [7.0,8.0)</td><td>b<sub>2</sub>b<sub>1</sub>b<sub>-1</sub>...b<sub>-13</sub>b<sub>-14</sub></td></tr><tr><td>γ<sub>n</sub><sup>(1)</sup></td><td>⟨52, 49⟩</td><td>xx1(.)xxx...xx</td><td>[1.0,2.0), [3.0,4.0), [5.0,6.0), and [7.0,8.0)</td><td>b<sub>2</sub>b<sub>1</sub>b<sub>-1</sub>...b<sub>-13</sub>b<sub>-14</sub></td></tr><tr><td>x<sub>n</sub><sup>(2)</sup></td><td>⟨50, 49⟩</td><td>0(.)xxx...x1</td><td>[2<sup>-49</sup>,1.0)</td><td>b<sub>-1</sub>b<sub>-2</sub>...b<sub>-15</sub>b<sub>-16</sub></td></tr><tr><td>x<sub>n</sub><sup>(1)</sup></td><td>⟨50, 49⟩</td><td>0(.)xxx...x1</td><td>[2<sup>-49</sup>,1.0)</td><td>b<sub>-1</sub>b<sub>-2</sub>...b<sub>-15</sub>b<sub>-16</sub></td></tr></tbody></table>

<table><caption>**Table 4.** Chosen bit arrangement of $Y_2$ to construct amounts of perturbation for control parameters of the PCM.</caption><thead><tr><th>Parameter</th><th>Bit arrangement Y<sub>i</sub></th></tr></thead><tbody><tr><td>δ<sub>γ<sub>n</sub></sub><sup>(2)</sup></td><td>Y<sub>2</sub><sup>(a)</sup> = [(2, 49)(1, 48)B<sub>0</sub>(2, 47)(1, 46)...(1, 34)B<sub>0</sub>B<sub>0</sub>...B<sub>0</sub>]</td></tr><tr><td>δ<sub>γ<sub>n</sub></sub><sup>(1)</sup></td><td>Y<sub>2</sub><sup>(b)</sup> = [(1, 49)(2, 48)B<sub>0</sub>(1, 47)(2, 46)...(2, 34)B<sub>0</sub>B<sub>0</sub>...B<sub>0</sub>]</td></tr></tbody></table>

<table><caption>**Table 5.** Chosen initial values.</caption><thead><tr><th>Entity</th><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>Cat map</td><td>γ<sub>n</sub><sup>(2)</sup></td><td>5.1122334455</td></tr><tr><td>Cat map</td><td>γ<sub>n</sub><sup>(1)</sup></td><td>3.1122334455</td></tr><tr><td>Cat map</td><td>x<sub>0</sub><sup>(2)</sup></td><td>0.5544332211</td></tr><tr><td>Cat map</td><td>x<sub>0</sub><sup>(1)</sup></td><td>0.1122334456</td></tr><tr><td>SKS</td><td>S<sub>0</sub></td><td>0x12345678</td></tr><tr><td>BPP</td><td>p<sub>k</sub><sup>(0)</sup>, k = 1, ..., 8</td><td>[0xC, 0x53BD01, 0x7F05F2, 0x9EC7AD, 0x1A, 0x05D4CB, 0xF0, 0xE8]</td></tr><tr><td>BPP</td><td>c<sub>k</sub><sup>(0)</sup>, k = 1, ..., 8</td><td>[0xE4, 0xDE2F91, 0x57CD1A, 0xA0D14F, 0xC8, 0x1B6AA3, 0委员会, 0x9B]</td></tr></tbody></table>

* The SKS: From Table 2, the maximum number of pixels in each image is $N_{pix\_max} = 2^{20}$ or 2048 pixels. According to Eq. (11), the length of $sk_k$ is $L_{sk} = 68$ bits. As assumed that each of the state variables and control parameters of the PCM is perturbed by 16 bits as given in the leftmost column of Table 3, therefore, the length of bit sequences to perturb the PCM is $L_E = 64$ bits for the perturbed Cat map. In this example, the LFSR is chosen and, according to Eq. (12), the number of states of LFSR is $N = 32$ with the primitive polynomial $X^{32} + X^{22} + X^2 + 1$ as suggested in $^{72}$. The seed of LFSR is $S_0$ in hexadecimal form, as shown in Table 5.

* The BPP: Table 6 shows bits of chaotic values, $b_{-17}...b_{-48}$, being extracted and used to construct values of coordinates of destination pixels and pseudo-random values in Eqs. (13) and (14), respectively. More specifically, bits sequences $bs_k^{(t)}$ with $t \in \{I, J, S, D\}$ in the four leftmost columns of Table 6 are taken bits of chaotic values in the second column. Then, bit sequences $bs_k^{(t)}$ with $t \in \{I, J, S, D\}$ are XORed with portions of $sk_k$ (see Fig. 3b) to produce pseudo-random values for permutation and diffusion. The modulus is applied as given in Eqs. (13) and (14) for smaller images and smaller number of bits representing the pixels. The XOR operation is chosen for functions $F_k^{(IJ)}$ and $F_k^{(Φ)}$ in Eqs. (13) and (14) as $F_k^{(IJ)} = bs_k^{(t_1)} \oplus sk_k^{(t_1)} \oplus \bigoplus_{i=1}^{T} pc_k^{(i)}$ and $F_k^{(Φ)} = bs_k^{(t_2)} \oplus sk_k^{(t_2)} \oplus \bigoplus_{i=1}^{T} pc_k^{(i)}$, with $t_1 \in \{I, J\}$ and $t_2 \in \{S, D\}$. Here, the bit sequences $sk_k^{(t_1)}$ and $sk_k^{(t_2)}$ can be cut to shorten and fit the lengths of $bs_k^{(t_1)}$ and $bs_k^{(t_2)}$, respectively. Also, bit sequence $pc_k$ is padded with bits 1's to have the length equal to an integer number of times of the length of $bs_k^{(t)}$ with $t \in \{I, J, S, D\}$, then it is divided into T sub-sequences before XORing all together.
<table><caption>**Table 6.** Bit positions of $X_R$ to produce pseudo-random values.</caption><thead><tr><th rowspan="2">Value of</th><th rowspan="2">Bits of x<sub>n</sub><sup>(i)</sup></th><th rowspan="2">k</th><th colspan="4">Chosen</th></tr><tr><th>24 bits for bs<sub>k</sub><sup>(S)</sup></th><th>24 bits for bs<sub>k</sub><sup>(D)</sup></th><th>10 bits for bs<sub>k</sub><sup>(I)</sup></th><th>10 bits for bs<sub>k</sub><sup>(J)</sup></th></tr></thead><tbody><tr><td rowspan="4">x<sub>n</sub><sup>(2)</sup></td><td>b<sub>-17...b-47</sub></td><td>1</td><td>b<sub>-17...b-40</sub></td><td>b<sub>-18...b-41</sub></td><td>b<sub>-31...b-40</sub></td><td>b<sub>-32...b-41</sub></td></tr><tr><td></td><td>2</td><td>b<sub>-19...b-42</sub></td><td>b<sub>-20...b-43</sub></td><td>b<sub>-33...b-42</sub></td><td>b<sub>-34...b-43</sub></td></tr><tr><td></td><td>3</td><td>b<sub>-21...b-44</sub></td><td>b<sub>-22...b-45</sub></td><td>b<sub>-35...b-44</sub></td><td>b<sub>-36...b-45</sub></td></tr><tr><td></td><td>4</td><td>b<sub>-23...b-46</sub></td><td>b<sub>-24...b-47</sub></td><td>b<sub>-37...b-46</sub></td><td>b<sub>-38...b-47</sub></td></tr><tr><td rowspan="4">x<sub>n</sub><sup>(1)</sup></td><td>b<sub>-18...b-48</sub></td><td>5</td><td>b<sub>-18...b-41</sub></td><td>b<sub>-19...b-42</sub></td><td>b<sub>-32...b-41</sub></td><td>b<sub>-33...b-42</sub></td></tr><tr><td></td><td>6</td><td>b<sub>-20...b-43</sub></td><td>b<sub>-21...b-44</sub></td><td>b<sub>-34...b-43</sub></td><td>b<sub>-35...b-44</sub></td></tr><tr><td></td><td>7</td><td>b<sub>-22...b-45</sub></td><td>b<sub>-23...b-46</sub></td><td>b<sub>-36...b-45</sub></td><td>b<sub>-37...b-46</sub></td></tr><tr><td></td><td>8</td><td>b<sub>-24...b-47</sub></td><td>b<sub>-25...b-48</sub></td><td>b<sub>-38...b-47</sub></td><td>b<sub>-39...b-48</sub></td></tr></tbody></table>

<figure><img src="image_15.png" /><figcaption>**Fig. 6.** Ciphertext images in the example using Structure 1 in Fig. 1a encryption rounds $N_e = 1, ..., 5$.</figcaption></figure>

* The diffusion: Similarly to the functions in the BPP, the function $F_k^{(Diff)}(.)$ is the XOR operation among its entries in Eqs. (2) and (3) as $F_k^{(Diff)} = \hat{c}(i_k, j_k) \oplus c(i_k, j_k - 1) \oplus \phi_k^{(t)}$ and $F_k^{(Diff)} = p(i_k, j_k) \oplus c(i_k, j_k - 1) \oplus \phi_k^{(t)}$, respectively, with $t \in \{S, D\}$.

* The S-box: In this exemplar simulation, the 8-bit AES S-box$^{73}$ and its corresponding inverse S-box are chosen. In fact, a new S-box can be generated by using chaotic values. For pixels represented by 24-bit data, i.e. red, green, and blue colors, three consecutive substitutions using the same S-box are performed, such that each is for 8 bit data of pixels. Then, the outputs of the S-box are merged to return 24-bit data of pixels.
<figure><img src="image_16.png" /><figcaption>Fig. 7. Ciphertext images in the example using Structure 2 in Fig. 1b encryption rounds $N_e = 1, ..., 5$.</figcaption></figure>

## Simulation results

Here, the set of eight images is encrypted for several rounds of encryption, $n_k = 1..N_e$, with $k = 1...K$ and $N_e = 5$. Figures 6, 7 and 8 display the ciphertext images produced by the exemplar designs using Structures 1, 2, and 3, respectively. The ciphertext images look like random images in every first round of encryption. In addition, these ciphertext images were successfully decrypted for any selected individual ciphered images.

It should be noted that the three-color image is considered as a grayscale image with a pixel depth of 24 bits. However, the color image is separated into three 8-bit grayscale images for analysis.

## Statistical analysis

Firstly, the randomness test is performed using the NIST SP 800-22 test standard to confirm the quality of session keys. Then, a statistical analysis of ciphertext images is performed with two main metrics, namely information entropy and correlation coefficients. Statistical analysis is measured for individual channels of color images, whereas bits of color channels are concatenated for encryption. Table 7 shows the information entropy and correlation coefficients of the original images.

### NIST randomness test for session keys

Here, the randomness test is performed for the bit sequences generated by the SKS blocks to confirm the quality of the session keys. As shown in Fig. 1, the session keys are $sk_k$ with $k = 1, 2, ..., 8$, and more than 2 million bits of each session key $sk_k$ are recorded for the test. In this study, the NIST SP 800-22 test standard is used for the session key randomness test. Each session key $sk_k$ with a length of 165.457.600 bits is collected during encryption for the test. The test result is judged by $P-value$, and it is passed if $P-value > 0.01$. Table 8 shows the test result of $sk_k$ with $k = 1, 2, ..., 8$, and it is clear that all $P-value$ are greater than 0.01 or pass. In other words, all the session keys in the exemplar simulation have good randomness and are qualified for encryption.

### Information entropy

The information entropy (IE) of an image measures the randomness of the values of the pixels, and the ideal value of IE is equal to the number of bits that represent the pixels. An image is mostly random if its IE is close to the ideal value.
<figure><img src="image_17.png" /><figcaption>Fig. 8. Ciphertext images in the example using Structure 3 in Fig. 1c encryption rounds $N_e = 1, ..., 5$.</figcaption></figure>

<table><caption>Table 7. Statistical measures of original images.</caption><thead><tr><th rowspan="2">Channel</th><th colspan="3">Entropy</th><th colspan="3">Chi-square test (χ²)</th><th colspan="9">Correlation coefficients in</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>6.501</td><td>6.446</td><td>6.273</td><td>579,798.39</td><td>595,896.90</td><td>742,573.87</td><td>0.924</td><td>0.869</td><td>0.907</td><td>0.874</td><td>0.776</td><td>0.884</td><td>0.867</td><td>0.750</td><td>0.854</td></tr><tr><td>Flowers</td><td>7.383</td><td>7.235</td><td>7.364</td><td>206,427.64</td><td>272,706.08</td><td>192,953.29</td><td>0.972</td><td>0.951</td><td>0.953</td><td>0.972</td><td>0.950</td><td>0.953</td><td>0.955</td><td>0.922</td><td>0.926</td></tr><tr><td>Peppers</td><td>7.339</td><td>7.496</td><td>7.059</td><td>213,187.22</td><td>318,382.93</td><td>491,428.18</td><td>0.964</td><td>0.981</td><td>0.967</td><td>0.966</td><td>0.982</td><td>0.966</td><td>0.956</td><td>0.969</td><td>0.948</td></tr><tr><td>Goldhill</td><td>7.471</td><td>7.429</td><td>7.476</td><td>231,332.11</td><td>268,762.60</td><td>278,677.24</td><td>0.968</td><td>0.973</td><td>0.977</td><td>0.961</td><td>0.977</td><td>0.979</td><td>0.940</td><td>0.958</td><td>0.962</td></tr><tr><td>Black</td><td>0.000</td><td></td><td></td><td>124,093,200.00</td><td></td><td></td><td>NaN</td><td></td><td></td><td>NaN</td><td></td><td></td><td>NaN</td><td></td><td></td></tr><tr><td>Boats</td><td>7.088</td><td></td><td></td><td>586,046.22</td><td></td><td></td><td>0.969</td><td></td><td></td><td>0.978</td><td></td><td></td><td>0.950</td><td></td><td></td></tr><tr><td>Clock</td><td>6.706</td><td></td><td></td><td>282,061.56</td><td></td><td></td><td>0.956</td><td></td><td></td><td>0.974</td><td></td><td></td><td>0.939</td><td></td><td></td></tr><tr><td>White</td><td>0.000</td><td></td><td></td><td>109,630,537.22</td><td></td><td></td><td>NaN</td><td></td><td></td><td>NaN</td><td></td><td></td><td>NaN</td><td></td><td></td></tr></tbody></table>

In this example, *IE* of the color images is computed for each channel, so the ideal value of *IE* is 8 bits. Table 9 displays *IE* of encrypted images with different numbers of encryption rounds, $N_e = 1, ..., 5$. The values of *IE* of all encrypted images are almost equal to the ideal value (8 bits) in every encryption round. The information entropy of all plaintext images in the simulation changes quickly from small values in the original images as given in Table 7 to almost 8 after the first round of encryption, at $N_e = 1$. The average information entropy is measured in all ciphertext images and is unchanged with increasing number of encryption rounds. In all simulations, the average value of 7.9992 is very close to the ideal value of 8, which means that the values of the pixels in the ciphertext images are highly random and uniformly distributed.
<table><caption>**Table 8.** The result of NIST test.</caption><thead><tr><th>Sub-tests</th><th>sk<sub>1</sub></th><th>sk<sub>2</sub></th><th>sk<sub>3</sub></th><th>sk<sub>4</sub></th><th>sk<sub>5</sub></th><th>sk<sub>6</sub></th><th>sk<sub>7</sub></th><th>sk<sub>8</sub></th><th>Result</th></tr></thead><tbody><tr><td>Frequency</td><td>0.3505</td><td>0.9114</td><td>0.3505</td><td>0.0179</td><td>0.3505</td><td>0.7399</td><td>0.5341</td><td>0.1223</td><td>All Pass</td></tr><tr><td>Block Frequency (m=128)</td><td>0.7399</td><td>0.5341</td><td>0.3505</td><td>0.7399</td><td>0.3505</td><td>0.5341</td><td>0.1223</td><td>0.3505</td><td>All Pass</td></tr><tr><td>Cusum-Forward</td><td>0.5341</td><td>0.3505</td><td>0.3505</td><td>0.1223</td><td>0.2133</td><td>0.1223</td><td>0.5341</td><td>0.7399</td><td>All Pass</td></tr><tr><td>Cusum-Reverse</td><td>0.9114</td><td>0.3505</td><td>0.5341</td><td>0.1223</td><td>0.3505</td><td>0.3505</td><td>0.5341</td><td>0.1223</td><td>All Pass</td></tr><tr><td>Runs</td><td>0.9915</td><td>0.1223</td><td>0.7399</td><td>0.7399</td><td>0.3505</td><td>0.3505</td><td>0.2133</td><td>0.3505</td><td>All Pass</td></tr><tr><td>Longest Runs of Ones</td><td>0.7399</td><td>0.1223</td><td>0.5341</td><td>0.5341</td><td>0.2133</td><td>0.5341</td><td>0.9114</td><td>0.7399</td><td>All Pass</td></tr><tr><td>Rank</td><td>0.2133</td><td>0.1223</td><td>0.2133</td><td>0.7399</td><td>0.3505</td><td>0.0179</td><td>0.5341</td><td>0.0179</td><td>All Pass</td></tr><tr><td>FFT</td><td>0.5341</td><td>0.7399</td><td>0.3505</td><td>0.1223</td><td>0.5341</td><td>0.9915</td><td>0.9114</td><td>0.3505</td><td>All Pass</td></tr><tr><td>NonOverlapping Templates (m = 9)</td><td>0.4923</td><td>0.5030</td><td>0.5424</td><td>0.4888</td><td>0.4684</td><td>0.5091</td><td>0.4672</td><td>0.4534</td><td>All Pass</td></tr><tr><td>Overlapping Templates (m = 9)</td><td>0.3505</td><td>0.7399</td><td>0.9114</td><td>0.0669</td><td>0.7399</td><td>0.0352</td><td>0.5341</td><td>0.5341</td><td>All Pass</td></tr><tr><td>Universal</td><td>0.3353</td><td>0.4584</td><td>0.3655</td><td>0.5486</td><td>0.2255</td><td>0.5423</td><td>0.0313</td><td>0.4896</td><td>All Pass</td></tr><tr><td>Approximate Entropy (m = 10)</td><td>0.7399</td><td>0.9114</td><td>0.5341</td><td>0.5341</td><td>0.5341</td><td>0.3505</td><td>0.9915</td><td>0.3505</td><td>All Pass</td></tr><tr><td>Serial (m = 16) P - value<sub>1</sub></td><td>0.0669</td><td>0.9114</td><td>0.1223</td><td>0.2133</td><td>0.3505</td><td>0.5341</td><td>0.2133</td><td>0.1223</td><td>All Pass</td></tr><tr><td>Serial (m = 16) P - value<sub>2</sub></td><td>0.5341</td><td>0.5341</td><td>0.2133</td><td>0.2133</td><td>0.7399</td><td>0.7399</td><td>0.9114</td><td>0.3505</td><td>All Pass</td></tr><tr><td>Linear Complexity (M = 500)</td><td>0.2133</td><td>0.2133</td><td>0.9114</td><td>0.7399</td><td>0.1223</td><td>0.0089</td><td>0.1223</td><td>0.2133</td><td>All Pass</td></tr></tbody></table>

### Histogram analysis

The histogram of an image reflects the distribution of pixel values, and the statistical histogram can be measured using the chi-square test with a significance level of 0.05. The equation of $\chi^2$ is

$$ \chi^2 = \sum_{i=0}^{V-1} \frac{(O_i - E_i)^2}{E_i} \qquad (18) $$

where V is the number of gray level ($V = 256$ for 8-bit grayscale images), and $O_i$ and $E_i$ are respectively observed and expected occurrence frequencies of gray level $i$, with $0 \le i \le V - 1$. The expected occurrence frequency of 8-bit grayscale images is $E_i = \frac{M \times N}{V}$; $M \times N$ is the size of the images. The hypothesis test considers the significance of the histogram if it complies with a uniform distribution, and the test is accepted (or the histogram is uniformly distributed) if $\chi^2 \le \chi^2_\alpha(V - 1)$. In this example, the significance level $\alpha = 0.05$ is considered, and $\chi^2_{0.05}(255) = [293.247]$.

Table 10 gives the test results of the ciphertext images for the exemplar simulation. The value of $\chi^2$ is less than [!tb]. In other words, the histogram statistics of all ciphertext images pass the test, or the histogram of ciphertext images is uniformly distributed. It is clear that the chi-square values of the plain text images in Table 7 are large, but are quickly reduced to and smaller than $\chi^2_{0.05}(255)$. The average values are calculated on the overall ciphertext images, and the chi-square test with significant level $\alpha = 0.05$ and passed, even after the first encryption round, means that the histogram of ciphertext images are uniformly distributed.

### Correlation coefficients between adjacent pixels

It is well known that a high similarity in the values of adjacent pixels in an image makes the visual structure. In this example, the correlation coefficients (CCs) between pairs of adjacent pixels in three directions, that is, horizontal (H), vertical (V) and diagonal (D), is measured, and the Pearson correlation coefficient is used for measurements for grayscale images and every channel of color images. The closer the CCs are to zero, the more random the image is.

Tables 11, 12 and 13 show the CCs of the ciphertext images in the exemplar simulation of each design. The results show that the CC values are close to zero in every encryption round, which means that the ciphertext images. However, the values of Clock's CCs are the biggest in all examples. The average values are calculated on all the ciphertext images of the cohort. Even the average values are small, but they fluctuate around zero and regardless to number of encryption rounds. This means that the visual structure in the plaintext images, with large value of CCs as seen in Table 7, is destroyed in the ciphertext images, with small values in all directions of horizontal, vertical and diagonal, as seen in Table 11, 12 and 13.

Compared with the values of the original images given in Table 7, the IE, the Chi-square test and the CC of the plain images are significantly changed compared to those of the ciphertext images given in Table 9, 10, 11, 12 and 13.

### Security analysis

Let us consider the security aspects of the examples using the proposed structures using two crucial metrics, namely the secret key space and the sensitivity of the secret key.

### Space of the secret key

The space of the secret key is considered by the number of flippable bits which are initial values of parameters. In this simulation, the flippable bits of initial values of the Cat map ($X_0$ and $\Gamma_0$) and of the SKS ($S_0$) as well
<table><caption>**Table 9.** Information entropy of encrypted images.</caption><caption>Exemplar design using Structure 1 in Fig. 1a</caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th colspan="16">Channel</th></tr><tr><td>Baboon</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9991</td><td>7.9992</td><td>7.9993</td><td>7.9992</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9993</td></tr><tr><td>Flowers</td><td>7.9989</td><td>7.9990</td><td>7.9991</td><td>7.9990</td><td>7.9990</td><td>7.9989</td><td>7.9989</td><td>7.9990</td><td>7.9990</td><td>7.9989</td><td>7.9989</td><td>7.9990</td><td>7.9989</td><td>7.9989</td><td>7.9990</td></tr><tr><td>Peppers</td><td>7.9994</td><td>7.9993</td><td>7.9992</td><td>7.9994</td><td>7.9992</td><td>7.9993</td><td>7.9993</td><td>7.9994</td><td>7.9992</td><td>7.9994</td><td>7.9994</td><td>7.9993</td><td>7.9992</td><td>7.9994</td><td>7.9992</td></tr><tr><td>Goldhill</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9995</td><td>7.9996</td></tr><tr><td>Black</td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Boats</td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9995</td><td></td><td></td></tr><tr><td>Clock</td><td>7.9974</td><td></td><td></td><td>7.9972</td><td></td><td></td><td>7.9971</td><td></td><td></td><td>7.9970</td><td></td><td></td><td>7.9970</td><td></td><td></td></tr><tr><td>White</td><td>7.9995</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Average</td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td></tr><tr><th colspan="16">Exemplar design using Structure 2 in Fig. 1b</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr><tr><th colspan="16">Channel</th></tr><tr><td>Baboon</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9993</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9991</td><td>7.9991</td><td>7.9993</td><td>7.9993</td></tr><tr><td>Flowers</td><td>7.9989</td><td>7.9992</td><td>7.9988</td><td>7.9990</td><td>7.9989</td><td>7.9990</td><td>7.9988</td><td>7.9991</td><td>7.9990</td><td>7.9989</td><td>7.9991</td><td>7.9990</td><td>7.9990</td><td>7.9989</td><td>7.9991</td></tr><tr><td>Peppers</td><td>7.9993</td><td>7.9992</td><td>7.9994</td><td>7.9993</td><td>7.9994</td><td>7.9994</td><td>7.9993</td><td>7.9992</td><td>7.9993</td><td>7.9994</td><td>7.9993</td><td>7.9993</td><td>7.9992</td><td>7.9994</td><td>7.9994</td></tr><tr><td>Goldhill</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9995</td><td>7.9995</td><td>7.9996</td><td>7.9996</td></tr><tr><td>Black</td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9997</td><td></td><td></td><td>7.9997</td><td></td><td></td></tr><tr><td>Boats</td><td>7.9996</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Clock</td><td>7.9974</td><td></td><td></td><td>7.9976</td><td></td><td></td><td>7.9973</td><td></td><td></td><td>7.9972</td><td></td><td></td><td>7.9970</td><td></td><td></td></tr><tr><td>White</td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Average</td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td></tr><tr><th colspan="16">Exemplar design using Structure 3 in Fig. 1c</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr><tr><th colspan="16">Channel</th></tr><tr><td>Baboon</td><td>7.9992</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9993</td><td>7.9993</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9992</td><td>7.9992</td></tr><tr><td>Flowers</td><td>7.9991</td><td>7.9988</td><td>7.9991</td><td>7.9991</td><td>7.9990</td><td>7.9989</td><td>7.9988</td><td>7.9991</td><td>7.9990</td><td>7.9992</td><td>7.9991</td><td>7.9991</td><td>7.9989</td><td>7.9990</td><td>7.9991</td></tr><tr><td>Peppers</td><td>7.9993</td><td>7.9994</td><td>7.9993</td><td>7.9992</td><td>7.9993</td><td>7.9994</td><td>7.9993</td><td>7.9994</td><td>7.9994</td><td>7.9993</td><td>7.9994</td><td>7.9993</td><td>7.9992</td><td>7.9993</td><td>7.9993</td></tr><tr><td>Goldhill</td><td>7.9996</td><td>7.9996</td><td>7.9997</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9996</td><td>7.9997</td><td>7.9996</td><td>7.9996</td><td>7.9997</td><td>7.9996</td><td>7.9996</td><td>7.9996</td></tr><tr><td>Black</td><td>7.9997</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Boats</td><td>7.9995</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9995</td><td></td><td></td><td>7.9995</td><td></td><td></td></tr><tr><td>Clock</td><td>7.9974</td><td></td><td></td><td>7.9972</td><td></td><td></td><td>7.9976</td><td></td><td></td><td>7.9979</td><td></td><td></td><td>7.9974</td><td></td><td></td></tr><tr><td>White</td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td><td>7.9996</td><td></td><td></td></tr><tr><td>Average</td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9992</td><td></td><td></td><td>7.9993</td><td></td><td></td><td>7.9992</td><td></td><td></td></tr></tbody></table>

as of plaintext and ciphertext pixels ($c_k^{(0)}$ and $p_k^{(0)}$ with $k = 1..K$) in the context of encrypting eight images.

Table 14 shows 486 flippable bits that contribute by parameters. In other words, the secret key space is $2^{486}$, and the ciphers in this exemplar simulation based on the proposed structures can resist brute-force attacks using modern computers.

## Sensitivity of the secret key

The sensitivity of the secret key of a MIE algorithm can be assessed by measuring the difference between two sets of ciphertext images that are obtained from the same set of plain images using two almost identical secret keys. Two secret keys are said to be almost identical because one of them is produced by using the other with slight modification. In this simulation, the least significant bits of the parameters are changed in turn to make the smallest difference between pairs of secret keys, as shown in Table 15.

The number of pixels change rate (NPCR) and the unified averaged changed intensity (UACI)$^{74,75}$ are used to measure the sensitivity of the secret key in two sets of ciphertext images. Specifically, a pair of ciphertext images $C_1$ and $C_2$ are used to calculate NPCR and UACI as

$$NPCR = \frac{\sum_{x,y} D(x, y)}{M \times N} \times 100\%, \quad (19)$$

and
<table><caption>χ<sup>2</sup> of encrypted images in the exemplar simulation using Structure 1 in Fig. 1a</caption><thead><tr><th>N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>225.024</td><td>274.968</td><td>274.831</td><td>283.294</td><td>278.586</td><td>238.726</td><td>270.159</td><td>236.399</td><td>276.024</td><td>266.733</td><td>272.489</td><td>244.002</td><td>278.942</td><td>256.361</td><td>231.366</td></tr><tr><td>Flowers</td><td>287.042</td><td>241.166</td><td>229.438</td><td>248.549</td><td>261.338</td><td>287.673</td><td>266.248</td><td>251.737</td><td>252.365</td><td>264.305</td><td>277.905</td><td>251.847</td><td>266.132</td><td>276.033</td><td>250.014</td></tr><tr><td>Peppers</td><td>271.232</td><td>246.325</td><td>226.382</td><td>272.219</td><td>241.743</td><td>227.020</td><td>212.012</td><td>246.864</td><td>245.899</td><td>256.662</td><td>267.964</td><td>239.662</td><td>269.368</td><td>245.315</td><td>271.356</td></tr><tr><td>Goldhill</td><td>256.172</td><td>223.146</td><td>269.022</td><td>255.493</td><td>268.755</td><td>258.332</td><td>225.218</td><td>269.318</td><td>278.298</td><td>227.722</td><td>253.505</td><td>252.721</td><td>241.462</td><td>284.386</td><td>258.382</td></tr><tr><td>Black</td><td>258.190</td><td></td><td></td><td>275.968</td><td></td><td></td><td>266.930</td><td></td><td></td><td>242.293</td><td></td><td></td><td>266.213</td><td></td><td></td></tr><tr><td>Boats</td><td>232.903</td><td></td><td></td><td>244.562</td><td></td><td></td><td>246.557</td><td></td><td></td><td>264.909</td><td></td><td></td><td>299.651</td><td></td><td></td></tr><tr><td>Clock</td><td>237.016</td><td></td><td></td><td>251.578</td><td></td><td></td><td>260.742</td><td></td><td></td><td>272.211</td><td></td><td></td><td>270.281</td><td></td><td></td></tr><tr><td>White</td><td>267.814</td><td></td><td></td><td>262.724</td><td></td><td></td><td>283.123</td><td></td><td></td><td>258.653</td><td></td><td></td><td>215.717</td><td></td><td></td></tr><tr><td>Average</td><td>251.292 (Pass)</td><td></td><td></td><td>259.785 (Pass)</td><td></td><td></td><td>255.493 (Pass)</td><td></td><td></td><td>257.099 (Pass)</td><td></td><td></td><td>261.311 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>χ<sup>2</sup> of encrypted images in the exemplar simulation using Structure 2 in Fig. 1b</caption><thead><tr><th>N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>226.044</td><td>260.122</td><td>260.563</td><td>236.239</td><td>277.043</td><td>263.742</td><td>248.465</td><td>235.881</td><td>274.929</td><td>253.086</td><td>251.480</td><td>301.811</td><td>286.812</td><td>227.428</td><td>234.404</td></tr><tr><td>Flowers</td><td>267.646</td><td>205.408</td><td>292.405</td><td>248.334</td><td>287.783</td><td>259.227</td><td>295.285</td><td>231.308</td><td>252.419</td><td>273.642</td><td>233.257</td><td>262.540</td><td>243.087</td><td>282.856</td><td>227.670</td></tr><tr><td>Peppers</td><td>263.356</td><td>269.438</td><td>255.138</td><td>234.802</td><td>251.378</td><td>280.429</td><td>220.387</td><td>260.395</td><td>249.479</td><td>255.202</td><td>249.676</td><td>290.457</td><td>270.369</td><td>251.735</td><td>220.364</td></tr><tr><td>Goldhill</td><td>266.360</td><td>269.895</td><td>268.893</td><td>280.378</td><td>249.177</td><td>232.245</td><td>280.120</td><td>252.197</td><td>253.164</td><td>273.825</td><td>235.658</td><td>294.107</td><td>304.127</td><td>255.396</td><td>238.926</td></tr><tr><td>Black</td><td>241.076</td><td></td><td></td><td>251.418</td><td></td><td></td><td>276.914</td><td></td><td></td><td>230.526</td><td></td><td></td><td>216.151</td><td></td><td></td></tr><tr><td>Boats</td><td>233.224</td><td></td><td></td><td>262.757</td><td></td><td></td><td>250.315</td><td></td><td></td><td>261.549</td><td></td><td></td><td>230.815</td><td></td><td></td></tr><tr><td>Clock</td><td>240.406</td><td></td><td></td><td>222.055</td><td></td><td></td><td>244.469</td><td></td><td></td><td>255.430</td><td></td><td></td><td>273.727</td><td></td><td></td></tr><tr><td>White</td><td>211.514</td><td></td><td></td><td>240.997</td><td></td><td></td><td>238.493</td><td></td><td></td><td>261.427</td><td></td><td></td><td>245.593</td><td></td><td></td></tr><tr><td>Average</td><td>251.968 (Pass)</td><td></td><td></td><td>254.875 (Pass)</td><td></td><td></td><td>254.014 (Pass)</td><td></td><td></td><td>261.480 (Pass)</td><td></td><td></td><td>250.591 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Table 10. Chi-square test results.</caption><caption>χ<sup>2</sup> of encrypted images in the exemplar simulation using Structure 3 in Fig. 1c</caption><thead><tr><th>N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>258.603</td><td>235.392</td><td>280.506</td><td>253.197</td><td>276.450</td><td>265.444</td><td>259.392</td><td>251.960</td><td>233.137</td><td>243.712</td><td>258.878</td><td>274.037</td><td>269.745</td><td>261.978</td><td>270.443</td></tr><tr><td>Flowers</td><td>235.692</td><td>296.968</td><td>220.527</td><td>219.684</td><td>245.689</td><td>273.235</td><td>295.500</td><td>219.682</td><td>262.359</td><td>205.923</td><td>232.756</td><td>231.313</td><td>267.801</td><td>254.817</td><td>233.418</td></tr><tr><td>Peppers</td><td>265.548</td><td>249.135</td><td>266.452</td><td>220.354</td><td>227.686</td><td>267.629</td><td>255.489</td><td>275.327</td><td>256.135</td><td>247.028</td><td>255.613</td><td>239.657</td><td>228.358</td><td>273.516</td><td>265.902</td></tr><tr><td>Goldhill</td><td>231.418</td><td>249.874</td><td>218.995</td><td>222.326</td><td>247.528</td><td>237.083</td><td>267.491</td><td>278.048</td><td>203.367</td><td>241.728</td><td>271.877</td><td>219.586</td><td>245.952</td><td>265.682</td><td>253.891</td></tr><tr><td>Black</td><td>203.720</td><td></td><td></td><td>245.266</td><td></td><td></td><td>267.113</td><td></td><td></td><td>266.415</td><td></td><td></td><td>256.545</td><td></td><td></td></tr><tr><td>Boats</td><td>281.506</td><td></td><td></td><td>257.516</td><td></td><td></td><td>273.678</td><td></td><td></td><td>286.019</td><td></td><td></td><td>263.912</td><td></td><td></td></tr><tr><td>Clock</td><td>237.133</td><td></td><td></td><td>255.656</td><td></td><td></td><td>218.313</td><td></td><td></td><td>186.922</td><td></td><td></td><td>234.008</td><td></td><td></td></tr><tr><td>White</td><td>235.127</td><td></td><td></td><td>247.135</td><td></td><td></td><td>258.162</td><td></td><td></td><td>240.050</td><td></td><td></td><td>259.806</td><td></td><td></td></tr><tr><td>Average</td><td>247.912 (Pass)</td><td></td><td></td><td>247.617 (Pass)</td><td></td><td></td><td>254.697 (Pass)</td><td></td><td></td><td>243.845 (Pass)</td><td></td><td></td><td>256.611 (Pass)</td><td></td><td></td></tr></tbody></table>

$$ UACI = \frac{1}{N^2} \left[ \sum_{x,y} \frac{C_1(x,y) - C_2(x,y)}{255} \right] \times 100\%, \quad (20) $$

where $D(x, y) = 1$ if $C_1(x, y) \neq C_2(x, y)$, and $D(x, y) = 0$ if $C_1(x, y) = C_2(x, y)$. The values of *NPCR* and *UACI* are used to determine the random tests. A random test with the critical value $\alpha = 0.05$ is passed if $NPCR \ge 99.569\%$ and $UACI \in [33.282\%, 33.644\%]$ as mentioned in$^{75}$.

As given in Eq. (16), the initial plain pixels, $p_k^{(0)}$ with $k = 1, \dots, K$, are used to perturb the chaotic map in encryption for the last pixels in the final encryption round. So, the sensitivity of the initial plain pixels is not used to consider the sensitivity for encryption. Similarity, sensitivity of the initial ciphertext pixels is only considered for the decryption.

Tables 16, 17, 18, 19, 20 and 21 present *NPCR* and *UACI* ciphertexts to analyze the sensitivity of each secret key parameter. The tests are measured for individual images and for all images using average values. It is clear that most of the tests for individual images have been passed, and average values indicate that ciphertext images of the cohort have been passed for sensitivity of every parameter. Passing the test of average values is very important for the case of MIE. Here, the average values of *NPCR* and *UACI* fluctuate slightly around 99.6 and 33.4 for all changes in the secret key and regardless of the number of encryption rounds. The sensitivity test is passed for every parameter of the secret key; therefore, any parameters can be used equivalently to lengthen the
<table><caption>Correlation coefficients in horizontal</caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>0.0010</td><td>-0.0008</td><td>-0.0027</td><td>0.0025</td><td>0.0013</td><td>-0.0002</td></tr><tr><td>Flowers</td><td>0.0011</td><td>0.0005</td><td>-0.0023</td><td>0.0018</td><td>0.0005</td><td>-0.0011</td></tr><tr><td>Peppers</td><td>0.0012</td><td>-0.0002</td><td>0.0003</td><td>-0.0011</td><td>-0.0009</td><td>-0.0035</td></tr><tr><td>Goldhill</td><td>0.0017</td><td>0.0004</td><td>-0.0028</td><td>-0.0022</td><td>-0.0002</td><td>-0.0010</td></tr><tr><td>Black</td><td>-0.0007</td><td>-0.0007</td><td>-0.0008</td><td>-0.0011</td><td>-0.0008</td><td>-0.0011</td></tr><tr><td>Boats</td><td>-0.0004</td><td>0.0000</td><td>0.0000</td><td>0.0000</td><td>-0.0003</td><td>0.0001</td></tr><tr><td>Clock</td><td>0.0033</td><td>0.0013</td><td>0.0043</td><td>0.0043</td><td>-0.0054</td><td>-0.0050</td></tr><tr><td>White</td><td>0.0018</td><td>0.0001</td><td>-0.0003</td><td>-0.0003</td><td>-0.0003</td><td>0.0021</td></tr><tr><td>Average</td><td>0.0001</td><td>-0.0002</td><td>0.0006</td><td>0.0001</td><td>-0.0001</td><td>-0.0001</td></tr></tbody></table>

<table><caption>Correlation coefficients in vertical</caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>-0.0009</td><td>0.0017</td><td>-0.0052</td><td>-0.0022</td><td>0.0018</td><td>-0.0020</td></tr><tr><td>Flowers</td><td>0.0016</td><td>0.0004</td><td>-0.0022</td><td>0.0009</td><td>0.0011</td><td>0.0017</td></tr><tr><td>Peppers</td><td>-0.0013</td><td>-0.0006</td><td>-0.0002</td><td>-0.0003</td><td>-0.0008</td><td>-0.0005</td></tr><tr><td>Goldhill</td><td>0.0012</td><td>-0.0003</td><td>0.0002</td><td>0.0025</td><td>-0.0006</td><td>-0.0003</td></tr><tr><td>Black</td><td>0.0016</td><td>-0.0011</td><td>-0.0053</td><td>0.0041</td><td>-0.0020</td><td>0.0001</td></tr><tr><td>Boats</td><td>-0.0011</td><td>-0.0013</td><td>-0.0005</td><td>-0.0009</td><td>-0.0029</td><td>0.0005</td></tr><tr><td>Clock</td><td>0.0035</td><td>0.0000</td><td>-0.0053</td><td>-0.0022</td><td>-0.0029</td><td>0.0001</td></tr><tr><td>White</td><td>0.0000</td><td>0.0005</td><td>0.0001</td><td>0.0000</td><td>0.0000</td><td>0.0002</td></tr><tr><td>Average</td><td>-0.0005</td><td>-0.0001</td><td>0.0001</td><td>0.0000</td><td>-0.0007</td><td>-0.0005</td></tr></tbody></table>

<table><caption>Table 11. Correlation coefficients of encrypted images in the exemplar design using Structure 1 in Fig. 1a.</caption><caption>Correlation coefficients in diagonal</caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>-0.0018</td><td>0.0019</td><td>-0.0035</td><td>-0.0026</td><td>0.0012</td><td>-0.0012</td></tr><tr><td>Flowers</td><td>0.0043</td><td>-0.0011</td><td>-0.0018</td><td>0.0025</td><td>-0.0005</td><td>0.0006</td></tr><tr><td>Peppers</td><td>0.0003</td><td>-0.0002</td><td>-0.0012</td><td>0.0002</td><td>-0.0004</td><td>-0.0005</td></tr><tr><td>Goldhill</td><td>0.0002</td><td>-0.0004</td><td>0.0005</td><td>-0.0016</td><td>0.0020</td><td>0.0016</td></tr><tr><td>Black</td><td>-0.0021</td><td>-0.0004</td><td>-0.0007</td><td>-0.0004</td><td>0.0000</td><td>0.0000</td></tr><tr><td>Boats</td><td>-0.0002</td><td>-0.0014</td><td>-0.0007</td><td>-0.0014</td><td>-0.0004</td><td>0.0000</td></tr><tr><td>Clock</td><td>0.0013</td><td>0.0027</td><td>0.0027</td><td>-0.0022</td><td>-0.0015</td><td>-0.0003</td></tr><tr><td>White</td><td>-0.0044</td><td>-0.0007</td><td>-0.0007</td><td>0.0009</td><td>0.0012</td><td>0.0012</td></tr><tr><td>Average</td><td>-0.0005</td><td>0.0001</td><td>0.0001</td><td>-0.0009</td><td>-0.0005</td><td>-0.0004</td></tr></tbody></table>

number of bits of the secret key. In other words, the exemplar ciphers in this simulation that use the proposed structures are highly sensitive to the smallest changes in the secret key.

## Comparison to the simulation results of some algorithms of chaos-based MIE

Table 22 shows the comparison between the exemplar simulation using the proposed structures and some existing MIE algorithms based on chaos. Here, comparisons are made only on the statistical properties of ciphertexts and the security properties of algorithms, regardless of structural configuration, parallelism, techniques, and image-content dependency. It is clear that the statistical properties of the ciphertexts in this simulation, i.e. information entropy and correlation coefficients, are comparable to those of existing algorithms. Although the key space of the exemplar designs (2$^{486}$) is moderate compared to that of other algorithms, it is large enough to resist the brute-force attack using a modern computer. The randomness of ciphertexts and sensitivity of secret key in this simulation are measured with pass by means of the chi-square test and of NCPR and UAIC, while some other algorithms are not quantitatively reported.

The simulation shows by means of the exemplar ciphers that the proposed structures are effective for practical applications. Next, the hardware for Structure 1 is designed, simulated, verified, and implemented on the FPGA platform to confirm the feasibility of the proposed designs in hardware applications.
<table><caption>**Table 12.** Correlation coefficients of encrypted images in the exemplar design using Structure 2 in Fig. 1b.</caption><thead><tr><th colspan="15">Correlation coefficients in horizontal</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>-0.0007</td><td>-0.0006</td><td>-0.0022</td><td>-0.0036</td><td>0.0001</td><td>0.0024</td><td>-0.0008</td><td>0.0029</td><td>0.0025</td><td>0.0048</td><td>-0.0020</td><td>0.0004</td><td>-0.0010</td><td>-0.0006</td><td>-0.0020</td></tr><tr><td>Flowers</td><td>0.0046</td><td>0.0017</td><td>0.0016</td><td>-0.0026</td><td>0.0012</td><td>0.0009</td><td>0.0040</td><td>-0.0003</td><td>0.0013</td><td>0.0002</td><td>0.0003</td><td>0.0016</td><td>-0.0027</td><td>0.0020</td><td>-0.0022</td></tr><tr><td>Peppers</td><td>0.0007</td><td>0.0014</td><td>-0.0006</td><td>-0.0009</td><td>-0.0015</td><td>0.0016</td><td>0.0005</td><td>-0.0020</td><td>0.0016</td><td>-0.0029</td><td>-0.0001</td><td>-0.0006</td><td>-0.0001</td><td>0.0002</td><td>0.0002</td></tr><tr><td>Goldhill</td><td>0.0004</td><td>-0.0007</td><td>0.0007</td><td>-0.0014</td><td>0.0008</td><td>-0.0004</td><td>0.0001</td><td>-0.0037</td><td>-0.0002</td><td>0.0028</td><td>-0.0027</td><td>-0.0007</td><td>0.0010</td><td>0.0029</td><td>0.0017</td></tr><tr><td>Black</td><td>0.0007</td><td></td><td></td><td>-0.0010</td><td></td><td></td><td>-0.0009</td><td></td><td></td><td>-0.0012</td><td></td><td></td><td>0.0002</td><td></td><td></td></tr><tr><td>Boats</td><td>-0.0024</td><td></td><td></td><td>-0.0018</td><td></td><td></td><td>-0.0006</td><td></td><td></td><td>-0.0015</td><td></td><td></td><td>-0.0005</td><td></td><td></td></tr><tr><td>Clock</td><td>0.0013</td><td></td><td></td><td>0.0084</td><td></td><td></td><td>0.0018</td><td></td><td></td><td>0.0071</td><td></td><td></td><td>-0.0012</td><td></td><td></td></tr><tr><td>White</td><td>0.0021</td><td></td><td></td><td>0.0019</td><td></td><td></td><td>-0.0033</td><td></td><td></td><td>0.0027</td><td></td><td></td><td>-0.0005</td><td></td><td></td></tr><tr><td>Average</td><td>0.0005</td><td></td><td></td><td>0.0003</td><td></td><td></td><td>0.0002</td><td></td><td></td><td>0.0005</td><td></td><td></td><td>-0.0002</td><td></td><td></td></tr></tbody><thead><tr><th colspan="15">Correlation coefficients in vertical</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>-0.0002</td><td>0.0000</td><td>0.0028</td><td>0.0009</td><td>-0.0040</td><td>-0.0005</td><td>-0.0029</td><td>0.0008</td><td>0.0022</td><td>-0.0053</td><td>-0.0003</td><td>0.0017</td><td>-0.0011</td><td>-0.0027</td><td>-0.0022</td></tr><tr><td>Flowers</td><td>0.0020</td><td>-0.0021</td><td>-0.0024</td><td>0.0009</td><td>0.0037</td><td>-0.0047</td><td>0.0011</td><td>-0.0013</td><td>0.0046</td><td>0.0000</td><td>-0.0004</td><td>0.0023</td><td>0.0016</td><td>0.0051</td><td>-0.0015</td></tr><tr><td>Peppers</td><td>-0.0003</td><td>-0.0004</td><td>0.0005</td><td>-0.0002</td><td>-0.0005</td><td>0.0004</td><td>0.0009</td><td>0.0014</td><td>0.0003</td><td>-0.0010</td><td>0.0026</td><td>0.0009</td><td>-0.0001</td><td>0.0023</td><td>-0.0029</td></tr><tr><td>Goldhill</td><td>-0.0005</td><td>0.0008</td><td>-0.0009</td><td>0.0003</td><td>0.0010</td><td>0.0002</td><td>-0.0019</td><td>0.0002</td><td>-0.0010</td><td>-0.0020</td><td>-0.0013</td><td>0.0014</td><td>0.0025</td><td>-0.0023</td><td>0.0001</td></tr><tr><td>Black</td><td>0.0008</td><td></td><td></td><td>0.0009</td><td></td><td></td><td>-0.0002</td><td></td><td></td><td>0.0000</td><td></td><td></td><td>0.0001</td><td></td><td></td></tr><tr><td>Boats</td><td>0.0020</td><td></td><td></td><td>0.0018</td><td></td><td></td><td>0.0008</td><td></td><td></td><td>-0.0002</td><td></td><td></td><td>-0.0017</td><td></td><td></td></tr><tr><td>Clock</td><td>-0.0017</td><td></td><td></td><td>0.0011</td><td></td><td></td><td>0.0041</td><td></td><td></td><td>-0.0025</td><td></td><td></td><td>0.0086</td><td></td><td></td></tr><tr><td>White</td><td>0.0007</td><td></td><td></td><td>0.0022</td><td></td><td></td><td>-0.0021</td><td></td><td></td><td>0.0004</td><td></td><td></td><td>-0.0007</td><td></td><td></td></tr><tr><td>Average</td><td>0.0001</td><td></td><td></td><td>0.0002</td><td></td><td></td><td>0.0004</td><td></td><td></td><td>-0.0002</td><td></td><td></td><td>0.0003</td><td></td><td></td></tr></tbody><thead><tr><th colspan="15">Correlation coefficients in diagonal</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>0.0012</td><td>0.0009</td><td>-0.0013</td><td>-0.0007</td><td>-0.0005</td><td>0.0013</td><td>-0.0002</td><td>-0.0018</td><td>0.0007</td><td>0.0003</td><td>-0.0028</td><td>-0.0005</td><td>-0.0005</td><td>0.0004</td><td>-0.0013</td></tr><tr><td>Flowers</td><td>0.0012</td><td>0.0030</td><td>-0.0021</td><td>-0.0022</td><td>0.0006</td><td>0.0012</td><td>0.0018</td><td>-0.0031</td><td>0.0050</td><td>0.0049</td><td>-0.0057</td><td>0.0016</td><td>0.0036</td><td>-0.0012</td><td>-0.0021</td></tr><tr><td>Peppers</td><td>0.0019</td><td>0.0015</td><td>0.0022</td><td>-0.0004</td><td>0.0001</td><td>-0.0009</td><td>-0.0002</td><td>-0.0017</td><td>0.0002</td><td>0.0004</td><td>-0.0015</td><td>0.0009</td><td>0.0007</td><td>-0.0002</td><td>0.0004</td></tr><tr><td>Goldhill</td><td>-0.0006</td><td>-0.0031</td><td>0.0011</td><td>0.0008</td><td>-0.0003</td><td>0.0007</td><td>0.0014</td><td>0.0007</td><td>-0.0007</td><td>-0.0010</td><td>0.0003</td><td>-0.0038</td><td>0.0000</td><td>-0.0039</td><td>0.0011</td></tr><tr><td>Black</td><td>-0.0026</td><td></td><td></td><td>-0.0019</td><td></td><td></td><td>0.0016</td><td></td><td></td><td>-0.0023</td><td></td><td></td><td>0.0006</td><td></td><td></td></tr><tr><td>Boats</td><td>-0.0006</td><td></td><td></td><td>-0.0007</td><td></td><td></td><td>-0.0015</td><td></td><td></td><td>-0.0005</td><td></td><td></td><td>-0.0011</td><td></td><td></td></tr><tr><td>Clock</td><td>-0.0122</td><td></td><td></td><td>0.0027</td><td></td><td></td><td>-0.0021</td><td></td><td></td><td>-0.0013</td><td></td><td></td><td>-0.0032</td><td></td><td></td></tr><tr><td>White</td><td>-0.0002</td><td></td><td></td><td>0.0005</td><td></td><td></td><td>0.0012</td><td></td><td></td><td>0.0024</td><td></td><td></td><td>0.0023</td><td></td><td></td></tr><tr><td>Average</td><td>-0.0006</td><td></td><td></td><td>0.0000</td><td></td><td></td><td>0.0001</td><td></td><td></td><td>-0.0005</td><td></td><td></td><td>-0.0003</td><td></td><td></td></tr></tbody></table>

## Digital hardware design

In this section, the hardware for MIE Structure 1 as illustrated in Fig. 1a is designed and implemented on the Intel FPGA device. The reason for choosing the FPGA platform is that the pixels of an individual image are permuted in the range of the image. The GPU only outperforms the FPGA if an image is processed in separate portions of pixels. In addition, the pipeline is implemented in FPGA to parallelize the encryption process for every individual image.

The hardware design uses the value of parameters for in the simulation mentioned above. The results on resource utilization, throughput, and hardware performance will also be evaluated and discussed. The hardware architecture of decryption is analogous to that of encryption, but the difference is only the reverse order of pixel processing compared to that of encryption.

## Hardware architecture of MIE

Figure 9 shows the MIE system architecture that incorporates MIE hardware as a core component, which is named MIE Intellectual Property (MIE-IP). The MIE-IP receives instructions from the processor through the AXI4-Lite communication interface to manage control information such as parameter setting for the encryption system, encryption requests, and read/write requests of data. The DDR4 SDRAM module is used to accommodate the image data for encryption or decryption. Both the processor and MIE-IP interface with the Memory Controller to access data in the DDR4 SDRAM. Furthermore, the application connects with the MIE system through high-speed data communication protocols such as PCIe or Ethernet. The operation of MIE System can be tracked through the JTAG communication interface.
<table><caption>Table 13. Correlation coefficients of encrypted images in the exemplar design using Structure 3 in Fig. 1c.</caption><thead><tr><th colspan="12">Correlation coefficients in horizontal</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>0.0019</td><td>0.0003</td><td>0.0015</td><td>0.0017</td><td>-0.0014</td><td>0.0005</td><td>0.0001</td><td>-0.0031</td><td>-0.0010</td><td>-0.0017</td><td>0.0001</td><td>-0.0007</td><td>-0.0033</td><td>0.0016</td><td>0.0004</td></tr><tr><td>Flowers</td><td>0.0040</td><td>0.0000</td><td>0.0008</td><td>-0.0007</td><td>0.0005</td><td>0.0048</td><td>0.0000</td><td>-0.0019</td><td>-0.0039</td><td>-0.0002</td><td>0.0025</td><td>-0.0040</td><td>-0.0015</td><td>-0.0024</td><td>-0.0045</td></tr><tr><td>Peppers</td><td>0.0002</td><td>-0.0017</td><td>0.0009</td><td>-0.0019</td><td>0.0003</td><td>0.0011</td><td>0.0003</td><td>0.0009</td><td>-0.0028</td><td>-0.0008</td><td>0.0016</td><td>-0.0042</td><td>0.0006</td><td>-0.0007</td><td>0.0003</td></tr><tr><td>Goldhill</td><td>0.0032</td><td>0.0017</td><td>0.0003</td><td>-0.0016</td><td>-0.0002</td><td>-0.0021</td><td>-0.0018</td><td>-0.0003</td><td>0.0005</td><td>0.0004</td><td>0.0001</td><td>-0.0007</td><td>0.0028</td><td>0.0010</td><td>0.0002</td></tr><tr><td>Black</td><td>-0.0007</td><td></td><td></td><td>0.0024</td><td></td><td></td><td>0.0023</td><td></td><td></td><td>0.0006</td><td></td><td></td><td>-0.0017</td><td></td><td></td></tr><tr><td>Boats</td><td>0.0003</td><td></td><td></td><td>0.0007</td><td></td><td></td><td>-0.0003</td><td></td><td></td><td>0.0015</td><td></td><td></td><td>-0.0041</td><td></td><td></td></tr><tr><td>Clock</td><td>-0.0027</td><td></td><td></td><td>-0.0025</td><td></td><td></td><td>0.0072</td><td></td><td></td><td>0.0067</td><td></td><td></td><td>-0.0027</td><td></td><td></td></tr><tr><td>White</td><td>-0.0012</td><td></td><td></td><td>0.0000</td><td></td><td></td><td>-0.0024</td><td></td><td></td><td>-0.0001</td><td></td><td></td><td>0.0004</td><td></td><td></td></tr><tr><td>Average</td><td>0.0006</td><td></td><td></td><td>0.0001</td><td></td><td></td><td>-0.0004</td><td></td><td></td><td>0.0001</td><td></td><td></td><td>-0.0009</td><td></td><td></td></tr><tr><th colspan="15">Correlation coefficients in vertical</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>0.0003</td><td>-0.0028</td><td>0.0003</td><td>0.0031</td><td>-0.0011</td><td>0.0018</td><td>0.0007</td><td>0.0014</td><td>0.0014</td><td>0.0020</td><td>0.0013</td><td>-0.0012</td><td>-0.0025</td><td>-0.0006</td><td>0.0002</td></tr><tr><td>Flowers</td><td>0.0013</td><td>0.0047</td><td>0.0046</td><td>0.0013</td><td>-0.0032</td><td>0.0003</td><td>-0.0043</td><td>-0.0008</td><td>0.0049</td><td>-0.0009</td><td>0.0022</td><td>0.0011</td><td>-0.0002</td><td>-0.0008</td><td>-0.0035</td></tr><tr><td>Peppers</td><td>0.0015</td><td>-0.0011</td><td>-0.0001</td><td>0.0008</td><td>-0.0019</td><td>0.0016</td><td>0.0025</td><td>0.0014</td><td>0.0018</td><td>-0.0004</td><td>0.0011</td><td>-0.0020</td><td>0.0017</td><td>0.0001</td><td>0.0002</td></tr><tr><td>Goldhill</td><td>0.0001</td><td>0.0014</td><td>-0.0024</td><td>0.0010</td><td>-0.0021</td><td>0.0026</td><td>0.0029</td><td>-0.0022</td><td>-0.0026</td><td>-0.0003</td><td>0.0005</td><td>-0.0023</td><td>-0.0008</td><td>0.0001</td><td>0.0007</td></tr><tr><td>Black</td><td>0.0006</td><td></td><td></td><td>0.0000</td><td></td><td></td><td>-0.0009</td><td></td><td></td><td>-0.0010</td><td></td><td></td><td>-0.0007</td><td></td><td></td></tr><tr><td>Boats</td><td>0.0002</td><td></td><td></td><td>-0.0006</td><td></td><td></td><td>0.0014</td><td></td><td></td><td>-0.0003</td><td></td><td></td><td>-0.0021</td><td></td><td></td></tr><tr><td>Clock</td><td>0.0069</td><td></td><td></td><td>-0.0030</td><td></td><td></td><td>0.0003</td><td></td><td></td><td>0.0031</td><td></td><td></td><td>0.0056</td><td></td><td></td></tr><tr><td>White</td><td>-0.0010</td><td></td><td></td><td>0.0020</td><td></td><td></td><td>-0.0005</td><td></td><td></td><td>-0.0014</td><td></td><td></td><td>-0.0005</td><td></td><td></td></tr><tr><td>Average</td><td>0.0009</td><td></td><td></td><td>0.0002</td><td></td><td></td><td>0.0005</td><td></td><td></td><td>0.0001</td><td></td><td></td><td>-0.0002</td><td></td><td></td></tr><tr><th colspan="15">Correlation coefficients in diagonal</th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr><tr><th colspan="15">Channel</th></tr><tr><td>Baboon</td><td>0.0008</td><td>0.0034</td><td>0.0015</td><td>-0.0015</td><td>0.0000</td><td>-0.0002</td><td>0.0003</td><td>0.0019</td><td>-0.0009</td><td>0.0012</td><td>-0.0018</td><td>0.0012</td><td>-0.0031</td><td>0.0006</td><td>0.0003</td></tr><tr><td>Flowers</td><td>-0.0024</td><td>0.0033</td><td>0.0055</td><td>0.0006</td><td>-0.0011</td><td>0.0000</td><td>-0.0003</td><td>-0.0018</td><td>0.0000</td><td>0.0007</td><td>0.0042</td><td>0.0056</td><td>0.0025</td><td>0.0015</td><td>-0.0057</td></tr><tr><td>Peppers</td><td>-0.0011</td><td>-0.0026</td><td>-0.0005</td><td>-0.0019</td><td>0.0004</td><td>0.0018</td><td>-0.0002</td><td>-0.0014</td><td>-0.0004</td><td>0.0017</td><td>-0.0020</td><td>-0.0016</td><td>0.0017</td><td>-0.0005</td><td>0.0003</td></tr><tr><td>Goldhill</td><td>-0.0026</td><td>0.0004</td><td>-0.0016</td><td>0.0009</td><td>-0.0006</td><td>0.0019</td><td>-0.0006</td><td>-0.0008</td><td>0.0006</td><td>-0.0008</td><td>-0.0013</td><td>-0.0007</td><td>-0.0006</td><td>-0.0004</td><td>0.0007</td></tr><tr><td>Black</td><td>0.0015</td><td></td><td></td><td>0.0008</td><td></td><td></td><td>-0.0027</td><td></td><td></td><td>0.0027</td><td></td><td></td><td>0.0005</td><td></td><td></td></tr><tr><td>Boats</td><td>0.0013</td><td></td><td></td><td>-0.0007</td><td></td><td></td><td>0.0010</td><td></td><td></td><td>0.0005</td><td></td><td></td><td>0.0029</td><td></td><td></td></tr><tr><td>Clock</td><td>-0.0067</td><td></td><td></td><td>0.0027</td><td></td><td></td><td>0.0003</td><td></td><td></td><td>-0.0002</td><td></td><td></td><td>-0.0017</td><td></td><td></td></tr><tr><td>White</td><td>-0.0011</td><td></td><td></td><td>-0.0020</td><td></td><td></td><td>0.0000</td><td></td><td></td><td>0.0004</td><td></td><td></td><td>0.0014</td><td></td><td></td></tr><tr><td>Average</td><td>-0.0001</td><td></td><td></td><td>0.0001</td><td></td><td></td><td>-0.0003</td><td></td><td></td><td>0.0006</td><td></td><td></td><td>0.0000</td><td></td><td></td></tr></tbody></table>

<table><caption>Table 14. Number of bits for the secret key.</caption><thead><tr><th>Entity</th><th>Parameters</th><th>Number of flipplable bits (bits)</th><th>Reference</th></tr></thead><tbody><tr><td>Cat map</td><td>γ<sub>0</sub><sup>(2)</sup></td><td>51</td><td>Eq. (17) and Table 3 and 5</td></tr><tr><td>Cat map</td><td>γ<sub>0</sub><sup>(1)</sup></td><td>51</td><td>Eq. (17) and Table 3 and 5</td></tr><tr><td>Cat map</td><td>x<sub>0</sub><sup>(2)</sup></td><td>48</td><td>Eq. (17) and Table 3 and 5</td></tr><tr><td>Cat map</td><td>x<sub>0</sub><sup>(1)</sup></td><td>48</td><td>Eq. (17) and Table 3 and 5</td></tr><tr><td>Initial value of SKS</td><td>S<sub>0</sub></td><td>32</td><td>Eq. (12) and Table 5</td></tr><tr><td>Initial values of ciphertext</td><td>{c<sub>1</sub><sup>(0)</sup>, c<sub>2</sub><sup>(0)</sup>, ..., c<sub>K</sub><sup>(0)</sup>}</td><td>∑<sub>k=1</sub><sup>K=8</sup> q<sub>k</sub>=128</td><td>Table 1</td></tr><tr><td>Initial values of plaintext</td><td>{p<sub>1</sub><sup>(0)</sup>, p<sub>2</sub><sup>(0)</sup>, ..., p<sub>K</sub><sup>(0)</sup>}</td><td>∑<sub>k=1</sub><sup>K=8</sup> q<sub>k</sub>=128</td><td>Table 1</td></tr><tr><td>Sum</td><td></td><td>486</td><td></td></tr></tbody></table>
<table><caption>**Table 15.** Modified secret key for the sensitivity analysis of secret key.</caption><thead><tr><td>Parameters</td><td>Bit being changed</td><td>amount of difference</td><td>Reference/Note</td></tr></thead><tbody><tr><td>γ<sub>0</sub><sup>(2)</sup></td><td>b<sub>-49</sub></td><td>2<sup>-49</sup></td><td>Table 3</td></tr><tr><td>γ<sub>0</sub><sup>(1)</sup></td><td>b<sub>-49</sub></td><td>2<sup>-49</sup></td><td>Table 3</td></tr><tr><td>x<sub>0</sub><sup>(2)</sup></td><td>b<sub>-48</sub></td><td>2<sup>-48</sup></td><td>Table 3</td></tr><tr><td>x<sub>0</sub><sup>(1)</sup></td><td>b<sub>-48</sub></td><td>2<sup>-48</sup></td><td>Table 3</td></tr><tr><td>S<sub>0</sub></td><td>b<sub>0</sub></td><td>1</td><td>Table 5</td></tr><tr><td>{c<sub>1</sub><sup>(0)</sup>, c<sub>2</sub><sup>(0)</sup>, ..., c<sub>K</sub><sup>(0)</sup>}</td><td>b<sub>0</sub></td><td>1</td><td>- Only one of initial ciphertext pixels changed</td></tr><tr><td>{p<sub>1</sub><sup>(0)</sup>, p<sub>2</sub><sup>(0)</sup>, ..., p<sub>K</sub><sup>(0)</sup>}</td><td>b<sub>0</sub></td><td>1</td><td>- Only one of initial plain pixels changed</td></tr></tbody></table>

## Operation of digital hardware

The MIE-IP hardware performs the encryption function described above. The execution time is balanced to have optimal speed for the entire system. So, the MIE system performs on blocks of $K \times N_{co}$ pixels in the internal memory, SRAM, of the FPGA chip. A block of pixels is illustrated in Fig. 10. A cohort of K pixels, each from an individual image, is encrypted at the same time. The encryption is repeated block by block until all images are encrypted. Encryption of a block is called an iteration. In each iteration, the hardware performs reading a block of data from the DDR4 SDRAM into the SRAM, encrypting blocks' pixels, and writing encrypted pixels from the SRAM into the DDR4 SDRAM.

### Memory reads/write

The scheme of consecutive read/write commands is used for reading/writing blocks of pixels between the SRAM and the DDR4 SDRAM. As an example for two consecutive read commands, the detail of the time is shown in Fig. 11. Specifically, the duration is $T_2 = t_{acc} + t_{read1} + t_{read2} + t_{diff}$ with $t_{diff} \ll t_{acc}$ and $t_{diff} \ll t_{read}$; where $t_{acc}$ is the time access of the read commands; $t_{read}$ is the time for all data to be read, and $t_{diff}$ is the delay between two read commands or between two groups of data from consecutive read commands. Multiple consecutive read commands can be performed to save time. This is also applicable to write commands.

### Operation phases of digital hardware

In this work, the allocation of time for the operating phases of the hardware is illustrated in Fig. 12. The encryption process is divided into multiple iterations and each iteration encrypts a block of pixels. An iteration consists of a series of phases, that is, Initial Key Phase and Normal Computation Phases, as shown in Fig. 12.

At the beginning of MIE, the SKS performs Initial Key Phase, and then the PCM iterates for the first value in Phase 1. Blocks of pixels are read from the DDR4 SDRAM into the SRAM, encrypted in SDRAM, and written back to the DDR4 SDRAM in Phases 2, 3, and 4. In order to optimize computation time, Session Key Scheduling, called Key Phase, for the next block is performed concurrently with Phases 2, 3, and 4 of the current block.

## MIE-IP architecture

The abstract architecture of MIE-IP is shown in Fig. 13. The design consists of five key components, as detailed in the following.

### Register block

Register Block (RB) is a collection of registers that are used to store parameters, operational information, and system status information. Processor (CPU) writes control information to and reads status information from the RB via the AXI4-Lite communication interface.

### Global computing

Global Computing (GC) consists of three main components: Session Key Scheduling, PCM Computing, and Position Generator.

Session Key Scheduling (SKS) works as an LFSR-based keystream generator. The number of states for the LFSR in the SKS is calculated based on Eq. (12), i.e., $N = 32$ states as in the simulation. However, the number of bits required to encrypt a cohort of K pixels is equal to the length of a frame given in Fig. 3, i.e., $L_F + K * L_{sk} = 608$ bits in this exemplar design. With an LFSR of 32 states, it requires 608 clock cycles to obtain 608 bits. To reduce the time delay in generating a sequence of 608 bits, this design implements the pipeline of 19 parallel-running LFSRs, and each LSFR is 32 bits in length. The seed $S_0$ of 32 bits is modified to become 19 seeds for 19 LFSRs. The bit sequences in the outputs of 19 LFSRs are merged to become 608 bits. In other words, the sequence of 608 bits is generated in only 32 clock cycles.

PCM Computing (PC) performs R chaotic iterations to get a value vector of $X_R$. So, a block of $N_{co}$ cohorts requires $N_{co}$ values of $X_R$. In order to optimize the speed of the hardware, the pipeline technique is applied to the hardware design for the PC, in which R PCMs are implemented and run in parallel using the same clock signal. However, such the R PCMs are put together in chain to reduce the computation time that the output of a PCM is the input of another. Therefore, instead of $R * N_{co}$ clock cycles to obtain $N_{co}$ values of $X_R$, it requires
<table><caption>Sensitivity on γ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.606</td><td>99.620</td><td>99.617</td><td>99.607</td><td>99.587</td><td>99.615</td><td>99.617</td><td>99.595</td><td>99.613</td><td>99.602</td><td>99.601</td><td>99.602</td><td>99.615</td><td>99.611</td><td>99.597</td></tr><tr><td>Flowers</td><td>99.606</td><td>99.603</td><td>99.610</td><td>99.610</td><td>99.604</td><td>99.587</td><td>99.610</td><td>99.617</td><td>99.586</td><td>99.622</td><td>99.599</td><td>99.601</td><td>99.597</td><td>99.603</td><td>99.610</td></tr><tr><td>Peppers</td><td>99.626</td><td>99.617</td><td>99.616</td><td>99.605</td><td>99.621</td><td>99.614</td><td>99.616</td><td>99.605</td><td>99.613</td><td>99.614</td><td>99.599</td><td>99.600</td><td>99.597</td><td>99.614</td><td>99.621</td></tr><tr><td>Goldhill</td><td>99.601</td><td>99.618</td><td>99.609</td><td>99.620</td><td>99.599</td><td>99.602</td><td>99.612</td><td>99.607</td><td>99.623</td><td>99.605</td><td>99.600</td><td>99.608</td><td>99.606</td><td>99.602</td><td>99.607</td></tr><tr><td>Black</td><td>99.613</td><td></td><td></td><td>99.619</td><td></td><td></td><td>99.612</td><td></td><td></td><td>99.592</td><td></td><td></td><td>99.614</td><td></td><td></td></tr><tr><td>Boats</td><td>99.602</td><td></td><td></td><td>99.631</td><td></td><td></td><td>99.613</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.618</td><td></td><td></td></tr><tr><td>Clock</td><td>99.635</td><td></td><td></td><td>99.661</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.599</td><td></td><td></td><td>99.622</td><td></td><td></td></tr><tr><td>White</td><td>99.606</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.607</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.600</td><td></td><td></td></tr><tr><td>Average</td><td>99.613 (Pass)</td><td></td><td></td><td>99.612 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.603 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on γ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.625</td><td>99.592</td><td>99.625</td><td>99.608</td><td>99.602</td><td>99.607</td><td>99.590</td><td>99.639</td><td>99.610</td><td>99.604</td><td>99.598</td><td>99.627</td><td>99.614</td><td>99.593</td><td>99.615</td></tr><tr><td>Flowers</td><td>99.608</td><td>99.595</td><td>99.585</td><td>99.636</td><td>99.626</td><td>99.635</td><td>99.636</td><td>99.577</td><td>99.610</td><td>99.614</td><td>99.591</td><td>99.583</td><td>99.624</td><td>99.634</td><td>99.606</td></tr><tr><td>Peppers</td><td>99.615</td><td>99.618</td><td>99.620</td><td>99.600</td><td>99.605</td><td>99.625</td><td>99.617</td><td>99.620</td><td>99.618</td><td>99.609</td><td>99.616</td><td>99.634</td><td>99.626</td><td>99.617</td><td>99.595</td></tr><tr><td>Goldhill</td><td>99.602</td><td>99.616</td><td>99.615</td><td>99.608</td><td>99.604</td><td>99.631</td><td>99.605</td><td>99.619</td><td>99.614</td><td>99.606</td><td>99.601</td><td>99.596</td><td>99.617</td><td>99.595</td><td>99.612</td></tr><tr><td>Black</td><td>99.600</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.596</td><td></td><td></td></tr><tr><td>Boats</td><td>99.603</td><td></td><td></td><td>99.616</td><td></td><td></td><td>99.588</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.615</td><td></td><td></td></tr><tr><td>Clock</td><td>99.647</td><td></td><td></td><td>99.641</td><td></td><td></td><td>99.576</td><td></td><td></td><td>99.586</td><td></td><td></td><td>99.641</td><td></td><td></td></tr><tr><td>White</td><td>99.607</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.598</td><td></td><td></td><td>99.618</td><td></td><td></td></tr><tr><td>Average</td><td>99.611 (Pass)</td><td></td><td></td><td>99.617 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td><td>99.605 (Pass)</td><td></td><td></td><td>99.614 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on α<sub>0</sub><sup>(2)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.605</td><td>99.591</td><td>99.607</td><td>99.615</td><td>99.585</td><td>99.590</td><td>99.624</td><td>99.615</td><td>99.615</td><td>99.598</td><td>99.584</td><td>99.605</td><td>99.612</td><td>99.614</td><td>99.633</td></tr><tr><td>Flowers</td><td>99.584</td><td>99.632</td><td>99.593</td><td>99.616</td><td>99.608</td><td>99.627</td><td>99.613</td><td>99.605</td><td>99.606</td><td>99.601</td><td>99.615</td><td>99.604</td><td>99.598</td><td>99.626</td><td>99.619</td></tr><tr><td>Peppers</td><td>99.613</td><td>99.620</td><td>99.617</td><td>99.618</td><td>99.599</td><td>99.603</td><td>99.615</td><td>99.590</td><td>99.616</td><td>99.615</td><td>99.607</td><td>99.601</td><td>99.616</td><td>99.597</td><td>99.602</td></tr><tr><td>Goldhill</td><td>99.603</td><td>99.610</td><td>99.617</td><td>99.606</td><td>99.611</td><td>99.597</td><td>99.619</td><td>99.611</td><td>99.610</td><td>99.615</td><td>99.608</td><td>99.623</td><td>99.588</td><td>99.592</td><td>99.616</td></tr><tr><td>Black</td><td>99.625</td><td></td><td></td><td>99.617</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.607</td><td></td><td></td><td>99.607</td><td></td><td></td></tr><tr><td>Boats</td><td>99.612</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.595</td><td></td><td></td><td>99.622</td><td></td><td></td><td>99.629</td><td></td><td></td></tr><tr><td>Clock</td><td>99.576</td><td></td><td></td><td>99.635</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.612</td><td></td><td></td><td>99.617</td><td></td><td></td></tr><tr><td>White</td><td>99.612</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.600</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.603</td><td></td><td></td></tr><tr><td>Average</td><td>99.607 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on α<sub>0</sub><sup>(1)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.606</td><td>99.591</td><td>99.628</td><td>99.605</td><td>99.605</td><td>99.602</td><td>99.594</td><td>99.611</td><td>99.611</td><td>99.624</td><td>99.602</td><td>99.598</td><td>99.610</td><td>99.605</td><td>99.598</td></tr><tr><td>Flowers</td><td>99.604</td><td>99.608</td><td>99.636</td><td>99.633</td><td>99.596</td><td>99.603</td><td>99.590</td><td>99.615</td><td>99.609</td><td>99.586</td><td>99.601</td><td>99.610</td><td>99.601</td><td>99.614</td><td>99.619</td></tr><tr><td>Peppers</td><td>99.621</td><td>99.609</td><td>99.610</td><td>99.622</td><td>99.618</td><td>99.624</td><td>99.609</td><td>99.602</td><td>99.608</td><td>99.612</td><td>99.623</td><td>99.597</td><td>99.593</td><td>99.620</td><td>99.611</td></tr><tr><td>Goldhill</td><td>99.617</td><td>99.609</td><td>99.614</td><td>99.619</td><td>99.612</td><td>99.609</td><td>99.619</td><td>99.626</td><td>99.605</td><td>99.612</td><td>99.628</td><td>99.592</td><td>99.610</td><td>99.596</td><td>99.612</td></tr><tr><td>Black</td><td>99.601</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.622</td><td></td><td></td><td>99.605</td><td></td><td></td></tr><tr><td>Boats</td><td>99.612</td><td></td><td></td><td>99.584</td><td></td><td></td><td>99.621</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.613</td><td></td><td></td></tr><tr><td>Clock</td><td>99.590</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.629</td><td></td><td></td><td>99.649</td><td></td><td></td><td>99.600</td><td></td><td></td></tr><tr><td>White</td><td>99.610</td><td></td><td></td><td>99.609</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.596</td><td></td><td></td><td>99.615</td><td></td><td></td></tr><tr><td>Average</td><td>99.610 (Pass)</td><td></td><td></td><td>99.612 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.610</td><td>99.606</td><td>99.608</td><td>99.602</td><td>99.610</td><td>99.615</td><td>99.620</td><td>99.623</td><td>99.603</td><td>99.623</td><td>99.604</td><td>99.622</td><td>99.617</td><td>99.613</td><td>99.602</td></tr><tr><td>Flowers</td><td>99.599</td><td>99.609</td><td>99.588</td><td>99.602</td><td>99.599</td><td>99.602</td><td>99.613</td><td>99.622</td><td>99.628</td><td>99.620</td><td>99.617</td><td>99.607</td><td>99.608</td><td>99.570</td><td>99.622</td><td>99.592</td></tr><tr><td>Peppers</td><td>99.618</td><td>99.610</td><td>99.590</td><td>99.607</td><td>99.620</td><td>99.599</td><td>99.612</td><td>99.625</td><td>99.629</td><td>99.620</td><td>99.617</td><td>99.608</td><td>99.622</td><td>99.621</td><td>99.610</td></tr><tr><td>Continued</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Goldhill</td><td>99.597</td><td>99.617</td><td>99.607</td><td>99.628</td><td>99.613</td><td>99.612</td><td>99.603</td><td>99.616</td><td>99.613</td><td>99.610</td><td>99.607</td><td>99.615</td><td>99.609</td><td>99.602</td><td>99.635</td></tr></tbody></table>
<table><caption>**Table 16.** Exemplar simulation for Structure 1 in Fig. 1a: Sensitivity of secret key by NPCR.</caption><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Black</td><td>99.618</td><td></td><td></td><td>99.618</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.590</td><td></td><td></td><td>99.604</td><td></td><td></td></tr><tr><td>Boats</td><td>99.616</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.609</td><td></td><td></td><td>99.596</td><td></td><td></td><td>99.603</td><td></td><td></td></tr><tr><td>Clock</td><td>99.602</td><td></td><td></td><td>99.657</td><td></td><td></td><td>99.594</td><td></td><td></td><td>99.574</td><td></td><td></td><td>99.580</td><td></td><td></td></tr><tr><td>White</td><td>99.617</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.608</td><td></td><td></td><td>99.621</td><td></td><td></td><td>99.620</td><td></td><td></td></tr><tr><td>Average</td><td>99.607 (Pass)</td><td></td><td></td><td>99.614 (Pass)</td><td></td><td></td><td>99.614 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td></tr></tbody></table>

only $R + N_{co} - 1$ clock cycles as Fig. 14 illustrates the timing diagram of PCM Computing, where values of $X_R$ are obtained at $t_R$ in different PCMs. In this exemplar design, $N_{co} = 8$ and $R = 5$ are chosen, so it only requires 13 clock cycles, instead of 39 clock cycles.

Position Generator performs the computation using Eqs. (13) and (14), and is implemented on hardware using only combinational logic gates, so it does not require any clock cycle in the timing diagram of the system.

### Image computing core

Each Image Computing Core (ICC) is responsible for encrypting pixels of an individual image. The more ICCs integrated into the system, the greater the number of images that can be encrypted simultaneously.

Each ICC composed of four primary components, i.e. Bit Pre-processing (BPP), Bit Interleaving (BI), Permutation-S-box-Diffusion (PSD), and Read-Write Image Controller (RWIC). The specific roles of BPP, BI, and PSD are detailed above. The RWIC component oversees the reading and writing processes of blocks of pixels from the Memory Read/Write Coordinator. The ICCs perform its computations in parallel, each ICC operates independently without interference from others.

### Memory read/write coordinator

Memory Read/Write Coordinator (MRWC) is responsible for coordinating data transfer between the Memory Controller via the AXI4 interface and the ICCs. The MRWC uses the SRAM on the FPGA chip to temporarily store image data retrieved from the DDR4 SDRAM in Phase 2 and to temporarily hold encrypted image data before being written to the DDR4 SDRAM by MIE-IP. In addition, this block also performs the crucial function of converting pixel coordinates into corresponding storage addresses in the DDR4 SDRAM, and vice versa. There are two specific data organization methods as shown in Fig. 15, where pixels with the same highlighted color correspond to the same image.

The first method in Fig. 15a organizes the images in a continuous way. The second method in Fig. 15b allocates pixels of the same image to a separate memory region where pixels of the same image are non-adjacent addresses in the DDR4 SDRAM. In fact, the first method provides better optimization for memory space compared to the second. However, the first method creates complexity in the coordinate-address transformation process, i.e., the conversion of pixel coordinates into storage addresses in the DDR4 SDRAM. Consequently, the computational complexity of the first method of data organization reduces the overall encryption speed. The second method requires less hardware resources, compared to that in the first one, for the coordinate-address transformation. Since the amount of image data is relatively modest in comparison to the space of the DDR4 SDRAM, the second method is chosen to optimize the speed of the hardware.

### Encryption controller

Encryption Controller (EC) oversees the entire operation of Global Computing, Image Computing Cores, and Memory Read/Write Coordinator. Specifically, the EC consists of subsystems with different responsibilities, that is, the Register Controller manages the read and write operations of Register Block; Global Controller coordinates the functions of Global Computing; Image Controller directs the activities of all Image Computing Cores; and the AXI4 Controller generates the control signals for the Memory Read/Write Coordinator.

Global Controller operates on the basis of a finite state machine (FSM) as depicted in Fig. 16. The AXI4 controller consists of two FSMs, as illustrated in Fig. 17 to manage the reading and writing processes of MIE-IP using the AXI4 protocol.

### FPGA implementation

The MIE system in Fig. 9 is implemented on the Han Pilot Platform and tested as shown in Fig. 18a to confirm the operation of the proposed design. The Han Pilot Platform is a development kit based on the Intel Arria 10 SoC FPGA from Intel-Altera. In this work, the code is written in SystemVerilog, then synthesized and implemented using Quartus Prime 21.1 Standard Edition. Functionality testing and performance evaluation are conducted for the implementation.

### FPGA-based system architecture

Figure 18b shows the details of the MIE system on the Han Pilot Platform for testing; where the Nios II CPU is used for general control and the JTAG-UART is utilized for console communication.
<table><caption>Sensitivity on γ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.603</td><td>99.598</td><td>99.586</td><td>99.610</td><td>99.600</td><td>99.609</td></tr><tr><td>Flowers</td><td>99.599</td><td>99.627</td><td>99.609</td><td>99.596</td><td>99.609</td><td>99.641</td></tr><tr><td>Peppers</td><td>99.625</td><td>99.604</td><td>99.619</td><td>99.622</td><td>99.608</td><td>99.627</td></tr><tr><td>Goldhill</td><td>99.610</td><td>99.605</td><td>99.589</td><td>99.599</td><td>99.608</td><td>99.606</td></tr><tr><td>Black</td><td>99.615</td><td></td><td></td><td>99.607</td><td></td><td>99.614</td></tr><tr><td>Boats</td><td>99.607</td><td></td><td></td><td>99.604</td><td></td><td>99.630</td></tr><tr><td>Clock</td><td>99.576</td><td></td><td></td><td>99.620</td><td></td><td>99.615</td></tr><tr><td>White</td><td>99.625</td><td></td><td></td><td>99.610</td><td></td><td>99.594</td></tr><tr><td>Average</td><td>99.606 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td>99.615 (Pass)</td></tr></tbody></table>

<table><caption>Sensitivity on γ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.601</td><td>99.618</td><td>99.622</td><td>99.615</td><td>99.617</td><td>99.602</td></tr><tr><td>Flowers</td><td>99.606</td><td>99.603</td><td>99.613</td><td>99.590</td><td>99.605</td><td>99.592</td></tr><tr><td>Peppers</td><td>99.621</td><td>99.609</td><td>99.615</td><td>99.587</td><td>99.621</td><td>99.615</td></tr><tr><td>Goldhill</td><td>99.607</td><td>99.599</td><td>99.612</td><td>99.595</td><td>99.619</td><td>99.614</td></tr><tr><td>Black</td><td>99.596</td><td></td><td></td><td>99.609</td><td></td><td>99.618</td></tr><tr><td>Boats</td><td>99.614</td><td></td><td></td><td>99.590</td><td></td><td>99.605</td></tr><tr><td>Clock</td><td>99.657</td><td></td><td></td><td>99.580</td><td></td><td>99.596</td></tr><tr><td>White</td><td>99.602</td><td></td><td></td><td>99.612</td><td></td><td>99.626</td></tr><tr><td>Average</td><td>99.612 (Pass)</td><td></td><td></td><td>99.604 (Pass)</td><td></td><td>99.611 (Pass)</td></tr></tbody></table>

<table><caption>Sensitivity on χ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.590</td><td>99.612</td><td>99.606</td><td>99.596</td><td>99.612</td><td>99.621</td></tr><tr><td>Flowers</td><td>99.605</td><td>99.591</td><td>99.609</td><td>99.622</td><td>99.610</td><td>99.619</td></tr><tr><td>Peppers</td><td>99.621</td><td>99.623</td><td>99.599</td><td>99.617</td><td>99.609</td><td>99.625</td></tr><tr><td>Goldhill</td><td>99.609</td><td>99.605</td><td>99.621</td><td>99.623</td><td>99.608</td><td>99.607</td></tr><tr><td>Black</td><td>99.588</td><td></td><td></td><td>99.613</td><td></td><td>99.612</td></tr><tr><td>Boats</td><td>99.602</td><td></td><td></td><td>99.606</td><td></td><td>99.595</td></tr><tr><td>Clock</td><td>99.559</td><td></td><td></td><td>99.599</td><td></td><td>99.657</td></tr><tr><td>White</td><td>99.598</td><td></td><td></td><td>99.623</td><td></td><td>99.624</td></tr><tr><td>Average</td><td>99.602 (Pass)</td><td></td><td></td><td>99.613 (Pass)</td><td></td><td>99.615 (Pass)</td></tr></tbody></table>

<table><caption>Sensitivity on χ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.608</td><td>99.624</td><td>99.592</td><td>99.600</td><td>99.611</td><td>99.612</td></tr><tr><td>Flowers</td><td>99.596</td><td>99.617</td><td>99.599</td><td>99.619</td><td>99.643</td><td>99.592</td></tr><tr><td>Peppers</td><td>99.623</td><td>99.607</td><td>99.616</td><td>99.603</td><td>99.614</td><td>99.625</td></tr><tr><td>Goldhill</td><td>99.592</td><td>99.619</td><td>99.625</td><td>99.615</td><td>99.617</td><td>99.593</td></tr><tr><td>Black</td><td>99.604</td><td></td><td></td><td>99.615</td><td></td><td>99.608</td></tr><tr><td>Boats</td><td>99.611</td><td></td><td></td><td>99.598</td><td></td><td>99.608</td></tr><tr><td>Clock</td><td>99.594</td><td></td><td></td><td>99.619</td><td></td><td>99.625</td></tr><tr><td>White</td><td>99.628</td><td></td><td></td><td>99.600</td><td></td><td>99.610</td></tr><tr><td>Average</td><td>99.610 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td>99.612 (Pass)</td></tr></tbody></table>

<table><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th>N<sub>e</sub></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Channel</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>99.624</td><td>99.605</td><td>99.615</td><td>99.614</td><td>99.619</td><td>99.597</td></tr><tr><td>Flowers</td><td>99.588</td><td>99.609</td><td>99.599</td><td>99.613</td><td>99.606</td><td>99.592</td></tr><tr><td>Peppers</td><td>99.589</td><td>99.600</td><td>99.597</td><td>99.605</td><td>99.623</td><td>99.616</td></tr><tr><td>Goldhill</td><td>99.592</td><td>99.618</td><td>99.603</td><td>99.623</td><td>99.608</td><td>99.601</td></tr></tbody></table>
<table><caption>**Table 17.** Exemplar simulation for Structure 2 in Fig. 1b: Sensitivity of secret key by NPCR.</caption><thead><tr><th colspan="11">Sensitivity on S<sub>0</sub></th></tr><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Black</td><td>99.620</td><td></td><td></td><td>99.615</td><td></td><td></td><td>99.621</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.620</td><td></td><td></td></tr><tr><td>Boats</td><td>99.611</td><td></td><td></td><td>99.607</td><td></td><td></td><td>99.595</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.606</td><td></td><td></td></tr><tr><td>Clock</td><td>99.632</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.657</td><td></td><td></td><td>99.600</td><td></td><td></td><td>99.611</td><td></td><td></td></tr><tr><td>White</td><td>99.620</td><td></td><td></td><td>99.599</td><td></td><td></td><td>99.593</td><td></td><td></td><td>99.594</td><td></td><td></td><td>99.603</td><td></td><td></td></tr><tr><td>Average</td><td>99.608 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.621 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.615 (Pass)</td><td></td><td></td></tr></tbody></table>

### Hardware metrics

The encryption and decryption processes of the FPGA-based system are simulated using Questa *-Intel® FPGA Edition software. After verifying the MIE-IP through functional simulation, the hardware metrics are assessed to demonstrate the utilization of the design resources, including the number of adaptive logic modules (ALM), digital signal processing blocks (DSP), registers, and block memory bits. The performance of the design is evaluated, including speed (maximum frequency, throughput), efficiency (throughput-to-ALM ratio), and power consumption.

* The throughput: Here, the throughput is computed for the maximum image data that goes through the hardware during the encryption time. It depends on, at maximum, the number of images that are encrypted at the same times, the size of the images, and the depth of the pixels as

$$Throughput = \frac{Q_{max} * K_{max} * M_{max} * N_{max}}{T_{total}} \quad (Mbps), \quad (21)$$

where $T_{total}$ is the encryption time for all data by only the ICCs. $T_{total}$ is counted by $T_{total} = \frac{N_{iterations} * T_{Nco} * N_{co}}{f_{max}}$, in which $f_{max}$ is the maximum frequency of the MIE-IP core; $N_{iterations}$ is the number of iterations needed to encrypt all pixels of the largest image, $N_{iterations} = \lceil \frac{M_{max} * N_{max}}{N_{co}} \rceil$; and $T_{Nco}$ is the number of clock cycles to encrypt plain pixels $p(i_k, j_k)$ and $p(i'_k, j'_k)$ for $k = 1..K$ as

$$T_{Nco} = T_{SDRAM\_to\_SRAM} + T_{SRAM\_to\_ICCs} + T_{ICCs\_to\_SRAM} + T_{SRAM\_to\_SDRAM}, \quad (22)$$

where $T_{SDRAM\_to\_SRAM}$ takes six clock cycles to write $p(i_k, j_k)$, $p(i'_k, j'_k)$, $c(i_k, j_k - 1)$, $c(i'_k, j'_k - 1)$, $p(i_k, j_k - 1)$, $p(i'_k, j'_k - 1)$ for $k = 1..K$ from the DDR4 SDRAM to the SRAM; $T_{SRAM\_to\_ICCs}$ is one clock cycle to read data from the SRAM to the ICCs; $T_{ICCs\_to\_SRAM}$ is one clock cycle to write the output of the ICCs to the SRAM; $T_{SRAM\_to\_SDRAM}$ are two clock cycles to write $c(i_k, j_k)$ and $c(i'_k, j'_k)$ back to the DDR4 SDRAM. Consequently, $T_{Nco}$ is ten clock cycles, and substituting $T_{Nco}$ into Eq.(21), the throughput is

$$Throughput = \frac{Q_{max} * K_{max} * M_{max} * N_{max} * f_{max}}{10N_{co} \lceil \frac{M_{max} * N_{max}}{N_{co}} \rceil} \quad (Mbps), \quad (23)$$

where $f_{max}$ is the maximum frequency of MIE hardware, and it has been found that the critical path lies in the fixed-point multiplication of the chaotic map implementation (PCM Computing in Fig. 13) in this exemplar design; $Q_{max}$ is the maximum number of bits representing a pixel; and $K$ is the number of images encrypted simultaneously. In this work, all hardware components in the design were optimized for the highest speed.

* The hardware efficiency: The hardware efficiency is computed as

$$Efficiency = \frac{Throughput}{ALMs} \quad (Mbps/ALM). \quad (24)$$

Tables 23 and 24 summarize hardware metrics and power consumption. The throughput is 998.2 Mbps measured with the value of the simulation parameters as $f_{max} = 51.99$ MHz, $K_{max} = 8$, $M_{max} = 1024$, $N_{max} = 1024$, $Q_{max} = 24$, and $N_{co} = 8$. These metrics indicate that the exemplar design requires modest resources on the FPGA chip, while it achieves high throughput, as given in Table 24. The low efficiency is obtained because there is hardware sacrifice to optimize the speed with the use of pipeline technique. Note from Table 23 that the amount of 96.1% RAM blocks is required for the entire MIE system, which is composed of several hardware modules such as NiOS, DRAM Controller, etc., for performing an end-to-end application. In fact, BRAM utilization is only counted for MIE-IP solely at 0.516%. The amount of hardware utilization for MIE-IP is less than 10% of the resource of the Intel Arria 10 SoC FPGA chip. The power consumption has been measured using the software tools by the FPGA chip manufacturer after complete compilation, placement, and routing.
<table><caption>Sensitivity on γ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Channel</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Baboon</td><td>99.606</td><td>99.620</td><td>99.617</td><td>99.607</td><td>99.587</td><td>99.615</td><td>99.617</td><td>99.595</td><td>99.613</td><td>99.602</td><td>99.601</td><td>99.602</td><td>99.615</td><td>99.611</td><td>99.597</td></tr><tr><td>Flowers</td><td>99.606</td><td>99.603</td><td>99.610</td><td>99.610</td><td>99.604</td><td>99.587</td><td>99.610</td><td>99.617</td><td>99.586</td><td>99.622</td><td>99.599</td><td>99.601</td><td>99.597</td><td>99.603</td><td>99.610</td></tr><tr><td>Peppers</td><td>99.619</td><td>99.620</td><td>99.617</td><td>99.601</td><td>99.621</td><td>99.624</td><td>99.623</td><td>99.609</td><td>99.607</td><td>99.621</td><td>99.596</td><td>99.607</td><td>99.612</td><td>99.621</td><td>99.626</td></tr><tr><td>Goldhill</td><td>99.601</td><td>99.618</td><td>99.609</td><td>99.620</td><td>99.599</td><td>99.602</td><td>99.612</td><td>99.607</td><td>99.623</td><td>99.605</td><td>99.600</td><td>99.608</td><td>99.606</td><td>99.602</td><td>99.607</td></tr><tr><td>Black</td><td>99.613</td><td></td><td></td><td>99.619</td><td></td><td></td><td>99.612</td><td></td><td></td><td>99.592</td><td></td><td></td><td>99.614</td><td></td><td></td></tr><tr><td>Boats</td><td>99.602</td><td></td><td></td><td>99.631</td><td></td><td></td><td>99.613</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.618</td><td></td><td></td></tr><tr><td>Clock</td><td>99.635</td><td></td><td></td><td>99.661</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.599</td><td></td><td></td><td>99.622</td><td></td><td></td></tr><tr><td>White</td><td>99.606</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.607</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.600</td><td></td><td></td></tr><tr><td>Average</td><td>99.613 (Pass)</td><td></td><td></td><td>99.612 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.604 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on γ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Channel</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Baboon</td><td>99.625</td><td>99.592</td><td>99.625</td><td>99.608</td><td>99.602</td><td>99.607</td><td>99.590</td><td>99.639</td><td>99.610</td><td>99.604</td><td>99.598</td><td>99.627</td><td>99.614</td><td>99.593</td><td>99.615</td></tr><tr><td>Flowers</td><td>99.608</td><td>99.595</td><td>99.585</td><td>99.636</td><td>99.626</td><td>99.635</td><td>99.636</td><td>99.577</td><td>99.610</td><td>99.614</td><td>99.591</td><td>99.583</td><td>99.624</td><td>99.634</td><td>99.606</td></tr><tr><td>Peppers</td><td>99.618</td><td>99.617</td><td>99.622</td><td>99.599</td><td>99.611</td><td>99.624</td><td>99.631</td><td>99.607</td><td>99.615</td><td>99.609</td><td>99.617</td><td>99.621</td><td>99.622</td><td>99.617</td><td>99.598</td></tr><tr><td>Goldhill</td><td>99.602</td><td>99.616</td><td>99.615</td><td>99.608</td><td>99.604</td><td>99.631</td><td>99.605</td><td>99.619</td><td>99.614</td><td>99.606</td><td>99.601</td><td>99.596</td><td>99.617</td><td>99.595</td><td>99.612</td></tr><tr><td>Black</td><td>99.600</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.596</td><td></td><td></td></tr><tr><td>Boats</td><td>99.603</td><td></td><td></td><td>99.616</td><td></td><td></td><td>99.588</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.615</td><td></td><td></td></tr><tr><td>Clock</td><td>99.647</td><td></td><td></td><td>99.641</td><td></td><td></td><td>99.576</td><td></td><td></td><td>99.586</td><td></td><td></td><td>99.641</td><td></td><td></td></tr><tr><td>White</td><td>99.607</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.598</td><td></td><td></td><td>99.618</td><td></td><td></td></tr><tr><td>Average</td><td>99.611 (Pass)</td><td></td><td></td><td>99.617 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td><td>99.604 (Pass)</td><td></td><td></td><td>99.614 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on α<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Channel</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Baboon</td><td>99.605</td><td>99.591</td><td>99.607</td><td>99.615</td><td>99.585</td><td>99.590</td><td>99.624</td><td>99.615</td><td>99.615</td><td>99.598</td><td>99.584</td><td>99.605</td><td>99.612</td><td>99.614</td><td>99.633</td></tr><tr><td>Flowers</td><td>99.584</td><td>99.632</td><td>99.593</td><td>99.616</td><td>99.608</td><td>99.627</td><td>99.613</td><td>99.605</td><td>99.606</td><td>99.601</td><td>99.615</td><td>99.604</td><td>99.598</td><td>99.626</td><td>99.619</td></tr><tr><td>Peppers</td><td>99.618</td><td>99.620</td><td>99.599</td><td>99.619</td><td>99.600</td><td>99.611</td><td>99.614</td><td>99.592</td><td>99.618</td><td>99.631</td><td>99.614</td><td>99.604</td><td>99.622</td><td>99.604</td><td>99.593</td></tr><tr><td>Goldhill</td><td>99.603</td><td>99.610</td><td>99.617</td><td>99.606</td><td>99.611</td><td>99.597</td><td>99.619</td><td>99.611</td><td>99.610</td><td>99.615</td><td>99.608</td><td>99.623</td><td>99.588</td><td>99.592</td><td>99.616</td></tr><tr><td>Black</td><td>99.625</td><td></td><td></td><td>99.617</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.607</td><td></td><td></td><td>99.607</td><td></td><td></td></tr><tr><td>Boats</td><td>99.612</td><td></td><td></td><td>99.603</td><td></td><td></td><td>99.595</td><td></td><td></td><td>99.622</td><td></td><td></td><td>99.629</td><td></td><td></td></tr><tr><td>Clock</td><td>99.576</td><td></td><td></td><td>99.635</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.612</td><td></td><td></td><td>99.617</td><td></td><td></td></tr><tr><td>White</td><td>99.612</td><td></td><td></td><td>99.604</td><td></td><td></td><td>99.600</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.603</td><td></td><td></td></tr><tr><td>Average</td><td>99.607 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.609 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on α<sub>0</sub><sup>(1)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Channel</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Baboon</td><td>99.606</td><td>99.591</td><td>99.628</td><td>99.605</td><td>99.605</td><td>99.602</td><td>99.594</td><td>99.611</td><td>99.611</td><td>99.624</td><td>99.602</td><td>99.598</td><td>99.610</td><td>99.605</td><td>99.598</td></tr><tr><td>Flowers</td><td>99.604</td><td>99.608</td><td>99.636</td><td>99.633</td><td>99.596</td><td>99.603</td><td>99.590</td><td>99.615</td><td>99.609</td><td>99.586</td><td>99.601</td><td>99.610</td><td>99.601</td><td>99.614</td><td>99.619</td></tr><tr><td>Peppers</td><td>99.615</td><td>99.599</td><td>99.609</td><td>99.620</td><td>99.617</td><td>99.619</td><td>99.621</td><td>99.600</td><td>99.591</td><td>99.601</td><td>99.621</td><td>99.623</td><td>99.607</td><td>99.608</td><td>99.619</td></tr><tr><td>Goldhill</td><td>99.617</td><td>99.609</td><td>99.614</td><td>99.619</td><td>99.612</td><td>99.609</td><td>99.619</td><td>99.626</td><td>99.605</td><td>99.612</td><td>99.628</td><td>99.592</td><td>99.610</td><td>99.596</td><td>99.612</td></tr><tr><td>Black</td><td>99.601</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.610</td><td></td><td></td><td>99.622</td><td></td><td></td><td>99.605</td><td></td><td></td></tr><tr><td>Boats</td><td>99.612</td><td></td><td></td><td>99.584</td><td></td><td></td><td>99.621</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.613</td><td></td><td></td></tr><tr><td>Clock</td><td>99.590</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.629</td><td></td><td></td><td>99.649</td><td></td><td></td><td>99.600</td><td></td><td></td></tr><tr><td>White</td><td>99.610</td><td></td><td></td><td>99.609</td><td></td><td></td><td>99.606</td><td></td><td></td><td>99.596</td><td></td><td></td><td>99.615</td><td></td><td></td></tr><tr><td>Average</td><td>99.609 (Pass)</td><td></td><td></td><td>99.612 (Pass)</td><td></td><td></td><td>99.610 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Channel</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Baboon</td><td>99.610</td><td>99.606</td><td>99.608</td><td>99.602</td><td>99.610</td><td>99.615</td><td>99.620</td><td>99.623</td><td>99.603</td><td>99.623</td><td>99.604</td><td>99.622</td><td>99.617</td><td>99.613</td><td>99.602</td></tr><tr><td>Flowers</td><td>99.599</td><td>99.609</td><td>99.588</td><td>99.602</td><td>99.599</td><td>99.602</td><td>99.613</td><td>99.622</td><td>99.628</td><td>99.617</td><td>99.607</td><td>99.608</td><td>99.570</td><td>99.622</td><td>99.592</td></tr><tr><td>Peppers</td><td>99.618</td><td>99.621</td><td>99.591</td><td>99.621</td><td>99.617</td><td>99.602</td><td>99.599</td><td>99.608</td><td>99.614</td><td>99.597</td><td>99.603</td><td>99.621</td><td>99.623</td><td>99.619</td><td>99.611</td></tr><tr><td>Continued</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Goldhill</td><td>99.597</td><td>99.617</td><td>99.607</td><td>99.628</td><td>99.613</td><td>99.612</td><td>99.603</td><td>99.616</td><td>99.613</td><td>99.610</td><td>99.607</td><td>99.615</td><td>99.609</td><td>99.602</td><td>99.635</td></tr></tbody></table>
<table><caption>**Table 18.** Exemplar simulation for Structure 3 in Fig. 1c: Sensitivity of secret key by NPCR.</caption><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub><br>Channel</th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Black</td><td>99.618</td><td></td><td></td><td>99.618</td><td></td><td></td><td>99.605</td><td></td><td></td><td>99.590</td><td></td><td></td><td>99.604</td><td></td><td></td></tr><tr><td>Boats</td><td>99.616</td><td></td><td></td><td>99.626</td><td></td><td></td><td>99.609</td><td></td><td></td><td>99.596</td><td></td><td></td><td>99.603</td><td></td><td></td></tr><tr><td>Clock</td><td>99.602</td><td></td><td></td><td>99.657</td><td></td><td></td><td>99.594</td><td></td><td></td><td>99.574</td><td></td><td></td><td>99.580</td><td></td><td></td></tr><tr><td>White</td><td>99.617</td><td></td><td></td><td>99.611</td><td></td><td></td><td>99.608</td><td></td><td></td><td>99.621</td><td></td><td></td><td>99.620</td><td></td><td></td></tr><tr><td>Average</td><td>99.608 (Pass)</td><td></td><td></td><td>99.615 (Pass)</td><td></td><td></td><td>99.611 (Pass)</td><td></td><td></td><td>99.607 (Pass)</td><td></td><td></td><td>99.608 (Pass)</td><td></td><td></td></tr></tbody></table>

* The latency: Considering the hardware latency, the analysis is based on the configurations of the proposed structures shown in Fig. 1 and their flowcharts in Fig. 4. The latency is defined as the time duration between the input of the first pixel of an image and the output of the first ciphertext pixel. It should be noted that the Bit Interleaving, Bit Pre-Processing (BPP), and Diffusion blocks are implemented with combinational logic (without flip-flops); therefore, these blocks operate combinationally and do not consume any clock cycles. Only the Permutation block requires memory read and write operations for pixel data. The S-box is also implemented with combinational logic and functions as a codebook; however, its delay depends on how it is utilized in each structure. Specifically, the S-boxes in Structures 1 and 2 do not incur any delay, whereas the S-boxes in Structure 3 require time for reading and writing pixel values from and to memory during the substitution process. In Structure 3, after the permutation and diffusion processes are completed for all pixels of an image, the substitution process is performed sequentially by reading the pixel values from memory and writing them back after substitution. Therefore, the latency of the proposed structures can be analyzed as follows. For Structures 1 and 2: The pixels in image $I_k$ are processed sequentially. In the exemplar design, both pixels $p(i_k, j_k + 1)$ and $c(i_k, j_k - 1)$ can be read simultaneously because they are stored in the same word of the SRAM. Therefore, only one SRAM read cycle is required. As described above, each permutation requires one SRAM read and one SRAM write cycle, denoted as $T_{reading\_SRAM}$ and $T_{writing\_SRAM}$, respectively. The latency is calculated as the time interval from when the first pixel $p(1, 1)$ enters the encryption core to when the first ciphertext pixel $c(1, 1)$ is produced. However, since the permutation operates over the entire image-space range, the overall latency corresponds to the time needed to complete encryption of all pixels in the smallest image across $N_e$ encryption rounds. The size of the smallest image is used in this calculation because its ciphertext pixels are returned first. Therefore, the latency of Structures 1 and 2 for $N_e$ rounds of encryption is given by

$$Latency_{1,2} = N_e * M_{min} * N_{min} * (2 * T_{reading\_SRAM} + T_{writing\_SRAM}) \quad (\text{clock cycles}) \quad (25)$$

where $M_{min}$ and $N_{min}$ are the size of smallest image in terms of the numbers of pixels. For Structure 3: Encryption is divided into two stages. In the first stage, all pixels are completely permuted and diffused. In the second stage, the pixels are substituted using the S-box. The total latency is determined by the time required to complete the permutation and diffusion processes over $N_e$ rounds, plus the time required for ($N_e - 1$) rounds of substitution. Therefore, the latency of Structure 3 is given by

$$Latency_3 = (2 * N_e - 1) * M_{min} * N_{min} * (T_{reading\_SRAM} + T_{writing\_SRAM}) \quad (\text{clock cycles}) \quad (26)$$

In FPGA, it only needs one clock cycle to read/write to SRAM, which means $T_{reading\_SRAM} = T_{writing\_SRAM} = 1/f_{max}$. Equations (25) and (26) become

$$Latency_{1,2} = \frac{3 * N_e * M_{min} * M_{min}}{f_{max}} \quad (\text{seconds}) \quad (27)$$

and

$$Latency_3 = \frac{2 * (2 * N_e - 1) * M_{min} * N_{min}}{f_{max}} \quad (\text{seconds}) \quad (28)$$

In the exemplar hardware design with $N_e = 5$, $M_{min} = N_{min} = 256$ (see Table 1), the estimated latencies using Eq. (27) for Structures 1 and 2 and Eq. (28) for Structure 3 are $Latency_{1,2} \approx 18.908$ milliseconds and $Latency_3 \approx 22.689$ milliseconds, respectively.

* Comparison to other works: Regardless of single- and multiple-image encryption, structure, or block and stream ciphers, Table 25 shows the results of throughput and efficiency between the exemplar design and other non-chaotic and chaos-based ciphers for reference. The referenced works are ordered by the implementation target and by years of publication. The results indicate that the proposed hardware achieves throughput comparable to that of recent ciphers, i.e., LSPT-SC<sup>77</sup> and LED (128 bit key) in<sup>78</sup>. The disadvantage is that the
<table><caption>Sensitivity on γ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>33.378</td><td>33.502</td><td>33.427</td><td>33.487</td><td>33.428</td><td>33.368</td><td>33.427</td><td>33.487</td><td>33.476</td><td>33.395</td><td>33.451</td><td>33.473</td><td>33.454</td><td>33.498</td><td>33.452</td></tr><tr><td>Flowers</td><td>33.514</td><td>33.456</td><td>33.465</td><td>33.461</td><td>33.462</td><td>33.519</td><td>33.469</td><td>33.518</td><td>33.406</td><td>33.461</td><td>33.552</td><td>33.519</td><td>33.383</td><td>33.424</td><td>33.529</td></tr><tr><td>Peppers</td><td>33.440</td><td>33.455</td><td>33.478</td><td>33.506</td><td>33.540</td><td>33.499</td><td>33.410</td><td>33.454</td><td>33.351</td><td>33.450</td><td>33.421</td><td>33.432</td><td>33.444</td><td>33.421</td><td>33.500</td></tr><tr><td>Goldhill</td><td>33.449</td><td>33.475</td><td>33.422</td><td>33.460</td><td>33.439</td><td>33.473</td><td>33.437</td><td>33.512</td><td>33.468</td><td>33.420</td><td>33.435</td><td>33.452</td><td>33.499</td><td>33.440</td><td>33.470</td></tr><tr><td>Black</td><td>33.492</td><td></td><td></td><td>33.457</td><td></td><td></td><td>33.446</td><td></td><td></td><td>33.532</td><td></td><td></td><td>33.511</td><td></td><td></td></tr><tr><td>Boats</td><td>33.443</td><td></td><td></td><td>33.485</td><td></td><td></td><td>33.510</td><td></td><td></td><td>33.480</td><td></td><td></td><td>33.460</td><td></td><td></td></tr><tr><td>Clock</td><td>33.580</td><td></td><td></td><td>33.568</td><td></td><td></td><td>33.447</td><td></td><td></td><td>33.474</td><td></td><td></td><td>33.360</td><td></td><td></td></tr><tr><td>White</td><td>33.490</td><td></td><td></td><td>33.401</td><td></td><td></td><td>33.525</td><td></td><td></td><td>33.505</td><td></td><td></td><td>33.468</td><td></td><td></td></tr><tr><td>Average</td><td>33.467 (Pass)</td><td></td><td></td><td>33.472 (Pass)</td><td></td><td></td><td>33.459 (Pass)</td><td></td><td></td><td>33.466 (Pass)</td><td></td><td></td><td>33.457 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on γ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>33.485</td><td>33.417</td><td>33.565</td><td>33.533</td><td>33.441</td><td>33.486</td><td>33.526</td><td>33.505</td><td>33.466</td><td>33.395</td><td>33.477</td><td>33.537</td><td>33.466</td><td>33.353</td><td>33.495</td></tr><tr><td>Flowers</td><td>33.419</td><td>33.503</td><td>33.422</td><td>33.531</td><td>33.471</td><td>33.470</td><td>33.458</td><td>33.539</td><td>33.435</td><td>33.474</td><td>33.491</td><td>33.509</td><td>33.428</td><td>33.606</td><td>33.337</td></tr><tr><td>Peppers</td><td>33.467</td><td>33.354</td><td>33.532</td><td>33.431</td><td>33.513</td><td>33.420</td><td>33.415</td><td>33.499</td><td>33.459</td><td>33.464</td><td>33.529</td><td>33.489</td><td>33.541</td><td>33.430</td><td>33.424</td></tr><tr><td>Goldhill</td><td>33.492</td><td>33.448</td><td>33.422</td><td>33.467</td><td>33.507</td><td>33.532</td><td>33.495</td><td>33.448</td><td>33.453</td><td>33.459</td><td>33.436</td><td>33.466</td><td>33.450</td><td>33.487</td><td>33.424</td></tr><tr><td>Black</td><td>33.493</td><td></td><td></td><td>33.477</td><td></td><td></td><td>33.463</td><td></td><td></td><td>33.451</td><td></td><td></td><td>33.441</td><td></td><td></td></tr><tr><td>Boats</td><td>33.490</td><td></td><td></td><td>33.494</td><td></td><td></td><td>33.483</td><td></td><td></td><td>33.489</td><td></td><td></td><td>33.465</td><td></td><td></td></tr><tr><td>Clock</td><td>33.643</td><td></td><td></td><td>33.481</td><td></td><td></td><td>33.398</td><td></td><td></td><td>33.564</td><td></td><td></td><td>33.425</td><td></td><td></td></tr><tr><td>White</td><td>33.476</td><td></td><td></td><td>33.454</td><td></td><td></td><td>33.438</td><td></td><td></td><td>33.446</td><td></td><td></td><td>33.469</td><td></td><td></td></tr><tr><td>Average</td><td>33.477 (Pass)</td><td></td><td></td><td>33.482 (Pass)</td><td></td><td></td><td>33.468 (Pass)</td><td></td><td></td><td>33.480 (Pass)</td><td></td><td></td><td>33.453 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on ω<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>33.443</td><td>33.405</td><td>33.436</td><td>33.426</td><td>33.391</td><td>33.423</td><td>33.515</td><td>33.496</td><td>33.520</td><td>33.519</td><td>33.483</td><td>33.371</td><td>33.411</td><td>33.475</td><td>33.512</td></tr><tr><td>Flowers</td><td>33.462</td><td>33.555</td><td>33.549</td><td>33.460</td><td>33.421</td><td>33.405</td><td>33.503</td><td>33.488</td><td>33.467</td><td>33.445</td><td>33.475</td><td>33.531</td><td>33.476</td><td>33.531</td><td>33.429</td></tr><tr><td>Peppers</td><td>33.495</td><td>33.531</td><td>33.488</td><td>33.387</td><td>33.549</td><td>33.455</td><td>33.448</td><td>33.511</td><td>33.520</td><td>33.459</td><td>33.470</td><td>33.467</td><td>33.433</td><td>33.424</td><td>33.445</td></tr><tr><td>Goldhill</td><td>33.469</td><td>33.463</td><td>33.445</td><td>33.453</td><td>33.448</td><td>33.392</td><td>33.497</td><td>33.512</td><td>33.444</td><td>33.489</td><td>33.504</td><td>33.450</td><td>33.420</td><td>33.437</td><td>33.449</td></tr><tr><td>Black</td><td>33.461</td><td></td><td></td><td>33.492</td><td></td><td></td><td>33.464</td><td></td><td></td><td>33.489</td><td></td><td></td><td>33.465</td><td></td><td></td></tr><tr><td>Boats</td><td>33.458</td><td></td><td></td><td>33.396</td><td></td><td></td><td>33.449</td><td></td><td></td><td>33.472</td><td></td><td></td><td>33.465</td><td></td><td></td></tr><tr><td>Clock</td><td>33.515</td><td></td><td></td><td>33.541</td><td></td><td></td><td>33.386</td><td></td><td></td><td>33.385</td><td></td><td></td><td>33.519</td><td></td><td></td></tr><tr><td>White</td><td>33.484</td><td></td><td></td><td>33.485</td><td></td><td></td><td>33.468</td><td></td><td></td><td>33.449</td><td></td><td></td><td>33.449</td><td></td><td></td></tr><tr><td>Average</td><td>33.479 (Pass)</td><td></td><td></td><td>33.445 (Pass)</td><td></td><td></td><td>33.481 (Pass)</td><td></td><td></td><td>33.466 (Pass)</td><td></td><td></td><td>33.459 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on ω<sub>0</sub><sup>(1)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>33.520</td><td>33.447</td><td>33.503</td><td>33.472</td><td>33.469</td><td>33.395</td><td>33.428</td><td>33.500</td><td>33.526</td><td>33.484</td><td>33.541</td><td>33.491</td><td>33.430</td><td>33.395</td><td>33.450</td></tr><tr><td>Flowers</td><td>33.429</td><td>33.483</td><td>33.398</td><td>33.396</td><td>33.450</td><td>33.453</td><td>33.456</td><td>33.415</td><td>33.544</td><td>33.583</td><td>33.418</td><td>33.517</td><td>33.411</td><td>33.405</td><td>33.288</td></tr><tr><td>Peppers</td><td>33.469</td><td>33.499</td><td>33.440</td><td>33.459</td><td>33.438</td><td>33.487</td><td>33.518</td><td>33.379</td><td>33.500</td><td>33.400</td><td>33.529</td><td>33.449</td><td>33.488</td><td>33.411</td><td>33.445</td></tr><tr><td>Goldhill</td><td>33.517</td><td>33.447</td><td>33.444</td><td>33.455</td><td>33.476</td><td>33.489</td><td>33.423</td><td>33.452</td><td>33.406</td><td>33.420</td><td>33.476</td><td>33.451</td><td>33.414</td><td>33.461</td><td>33.440</td></tr><tr><td>Black</td><td>33.516</td><td></td><td></td><td>33.490</td><td></td><td></td><td>33.441</td><td></td><td></td><td>33.485</td><td></td><td></td><td>33.532</td><td></td><td></td></tr><tr><td>Boats</td><td>33.454</td><td></td><td></td><td>33.390</td><td></td><td></td><td>33.505</td><td></td><td></td><td>33.530</td><td></td><td></td><td>33.457</td><td></td><td></td></tr><tr><td>Clock</td><td>33.532</td><td></td><td></td><td>33.471</td><td></td><td></td><td>33.494</td><td></td><td></td><td>33.548</td><td></td><td></td><td>33.258</td><td></td><td></td></tr><tr><td>White</td><td>33.496</td><td></td><td></td><td>33.492</td><td></td><td></td><td>33.416</td><td></td><td></td><td>33.451</td><td></td><td></td><td>33.532</td><td></td><td></td></tr><tr><td>Average</td><td>33.475 (Pass)</td><td></td><td></td><td>33.455 (Pass)</td><td></td><td></td><td>33.463 (Pass)</td><td></td><td></td><td>33.486 (Pass)</td><td></td><td></td><td>33.426 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Baboon</td><td>33.466</td><td>33.523</td><td>33.434</td><td>33.520</td><td>33.432</td><td>33.502</td><td>33.528</td><td>33.487</td><td>33.525</td><td>33.458</td><td>33.445</td><td>33.446</td><td>33.498</td><td>33.473</td><td>33.435</td></tr><tr><td>Flowers</td><td>33.450</td><td>33.456</td><td>33.519</td><td>33.409</td><td>33.457</td><td>33.484</td><td>33.462</td><td>33.487</td><td>33.525</td><td>33.478</td><td>33.388</td><td>33.551</td><td>33.441</td><td>33.579</td><td>33.440</td></tr><tr><td>Peppers</td><td>33.489</td><td>33.470</td><td>33.366</td><td>33.439</td><td>33.433</td><td>33.412</td><td>33.435</td><td>33.420</td><td>33.449</td><td>33.513</td><td>33.499</td><td>33.417</td><td>33.458</td><td>33.537</td><td>33.540</td></tr><tr><td>Continued</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Goldhill</td><td>33.527</td><td>33.423</td><td>33.525</td><td>33.469</td><td>33.419</td><td>33.474</td><td>33.501</td><td>33.483</td><td>33.474</td><td>33.422</td><td>33.478</td><td>33.442</td><td>33.478</td><td>33.453</td><td>33.512</td></tr></tbody></table>
<table><caption>**Table 19.** Exemplar simulation for Structure 1 in Fig. 1a: Sensitivity of secret key by UACI.</caption><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><td>Black</td><td>33.426</td><td></td><td></td><td>33.551</td><td></td><td></td><td>33.431</td><td></td><td></td><td>33.446</td><td></td><td></td><td>33.526</td><td></td><td></td></tr><tr><td>Boats</td><td>33.497</td><td></td><td></td><td>33.449</td><td></td><td></td><td>33.436</td><td></td><td></td><td>33.465</td><td></td><td></td><td>33.446</td><td></td><td></td></tr><tr><td>Clock</td><td>33.501</td><td></td><td></td><td>33.633</td><td></td><td></td><td>33.383</td><td></td><td></td><td>33.573</td><td></td><td></td><td>33.471</td><td></td><td></td></tr><tr><td>White</td><td>33.494</td><td></td><td></td><td>33.492</td><td></td><td></td><td>33.518</td><td></td><td></td><td>33.496</td><td></td><td></td><td>33.521</td><td></td><td></td></tr><tr><td>Average</td><td>33.447 (Pass)</td><td></td><td></td><td>33.448 (Pass)</td><td></td><td></td><td>33.483 (Pass)</td><td></td><td></td><td>33.466 (Pass)</td><td></td><td></td><td>33.460 (Pass)</td><td></td><td></td></tr></tbody></table>

efficiency of the exemplar design is significantly lower than the others, except for RC4. The exemplar design has low efficiency because the entire system in this work is implemented on hardware, rather than a combination of hardware and software as in AES-128⁷⁹. In addition, the throughput is also dependent on the level of security expected by means of number of encryption rounds. Therefore, these numbers of metrics in Table 25 are meant only for reference purposes.

## Discussion

The MIE structures have been proposed and verified through exemplar simulation and hardware design. Some discussion of structure, performance, security, and hardware improvements is given below.

In terms of structural configuration, the proposed structures are generic models, so that they can be scalable for as many images with different sizes and different pixel depths as possible. Here, the processes of reshaping and padding pixels are not required. In addition, any chaotic map can be employed in the proposed structures. Images of a cohort are encrypted separately and in parallel. Due to the separate encryption, it allows any ciphertext images from the cohort to be decrypted separately.

In terms of performance, the simulation shows that all images have the same number of encryption rounds, $N_e$. So, encryption for small images will be completed earlier than for larger ones. In fact, the proposed structures can encrypt different images with different numbers of encryption rounds by setting the condition for flags $Pause_k$. Unlike most existing algorithms of chaos-based MIE, the permutation and diffusion are combined in the same step, so it saves the time of chaotic iterations in generation of session keys.

In terms of security, with the inclusion of an S-box and session keys of image content dependency, the ciphers employing the proposed structures are sensitive to the image content and can resist any type of attack that exploits the difference in the plaintexts and ciphertexts. In addition, the key space can be expanded by lengthening the secret key in several ways, i.e., increasing the number of bits representing chaotic values, choosing a higher-dimensional chaotic map, and lengthening the size of the LFSR. However, these solutions require increasing computational complexity.

In terms of hardware performance, the hardware design for the example using Structure 1 has been synthesized and implemented on the FPGA platform, in which the value of parameters is chosen as given in the simulation. The timing of the hardware design is optimized for the speed of the MIE core. That is, the pipeline technique is used for the design of PCM, SKS, and MRWCs. The speed of the hardware is improved because the SRAM in the FPGA chip is used to temporarily store pixels for encryption. According to Eq. (23), for a certain maximum frequency of the MIE-IP core, $f_{max}$, the throughput can be scalable by the number of images to be encrypted at the same time, $K_{max}$, as well as by the number of bits representing pixels, $Q_{max}$. In fact, the trade-off between the throughput and the utilization of hardware resource must be considered in this case. It is noted that the other hardware designs listed in Table 25 cannot be scaled to have the desired throughput value. In addition, the hardware design based on Structures 2 and 3 is almost the same as that given in the exemplar implementation of Structure 1, except for the position of the S-box relative to the diffusion. The S-box is implemented as a codebook by logic gates, so it does not take any clock cycle.

## Conclusion

In summary, three chaos-based MIE models have been proposed that allow encrypting multiple images in parallel, and the images of the cohort can be various in sizes and pixel depths. Any selected individual ciphertext images from the cohort can be decrypted separately. The simulation results and hardware implementation for the example show the effectiveness and feasibility of the proposed structures. The proposed models are suitable for the application of diverse image data.

In these models, the encryption is performed pixel by pixel sequentially for an entire image, rather than blocks of pixels, therefore, it is sensitive to any change in the value of ciphertext pixels. This means that there are challenges for the types of attack based on the analysis of pairs of plaintexts and ciphertexts. However, the proposed structures are not robust to the ciphertext-crop attack. In future work, the proposed structures will be further developed to improve both the speed and efficiency of the hardware, and the security to resist the ciphertext-crop attack.
<table><caption>Sensitivity on γ<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th>Channel</th><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>Baboon</th><td>33.375</td><td>33.499</td><td>33.403</td><td>33.513</td><td>33.421</td><td>33.462</td><td>33.549</td><td>33.462</td><td>33.426</td><td>33.391</td><td>33.476</td><td>33.471</td><td>33.422</td><td>33.364</td><td>33.465</td></tr><tr><th>Flowers</th><td>33.432</td><td>33.414</td><td>33.429</td><td>33.483</td><td>33.447</td><td>33.495</td><td>33.475</td><td>33.578</td><td>33.567</td><td>33.459</td><td>33.535</td><td>33.415</td><td>33.369</td><td>33.431</td><td>33.440</td></tr><tr><th>Peppers</th><td>33.419</td><td>33.420</td><td>33.511</td><td>33.547</td><td>33.509</td><td>33.455</td><td>33.389</td><td>33.531</td><td>33.460</td><td>33.477</td><td>33.466</td><td>33.507</td><td>33.478</td><td>33.445</td><td>33.446</td></tr><tr><th>Goldhill</th><td>33.432</td><td>33.470</td><td>33.445</td><td>33.436</td><td>33.455</td><td>33.453</td><td>33.432</td><td>33.444</td><td>33.471</td><td>33.463</td><td>33.424</td><td>33.498</td><td>33.454</td><td>33.495</td><td>33.495</td></tr><tr><th>Black</th><td>33.510</td><td></td><td></td><td>33.416</td><td></td><td></td><td>33.531</td><td></td><td></td><td>33.535</td><td></td><td></td><td>33.492</td><td></td><td></td></tr><tr><th>Boats</th><td>33.442</td><td></td><td></td><td>33.504</td><td></td><td></td><td>33.451</td><td></td><td></td><td>33.480</td><td></td><td></td><td>33.409</td><td></td><td></td></tr><tr><th>Clock</th><td>33.572</td><td></td><td></td><td>33.403</td><td></td><td></td><td>33.456</td><td></td><td></td><td>33.513</td><td></td><td></td><td>33.497</td><td></td><td></td></tr><tr><th>White</th><td>33.490</td><td></td><td></td><td>33.396</td><td></td><td></td><td>33.441</td><td></td><td></td><td>33.511</td><td></td><td></td><td>33.513</td><td></td><td></td></tr><tr><th>Average</th><td>33.454 (Pass)</td><td></td><td></td><td>33.462 (Pass)</td><td></td><td></td><td>33.479 (Pass)</td><td></td><td></td><td>33.476 (Pass)</td><td></td><td></td><td>33.451 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on γ<sub>0</sub><sup>(1)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th>Channel</th><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>Baboon</th><td>33.420</td><td>33.398</td><td>33.506</td><td>33.502</td><td>33.445</td><td>33.379</td><td>33.566</td><td>33.486</td><td>33.497</td><td>33.416</td><td>33.438</td><td>33.554</td><td>33.489</td><td>33.536</td><td>33.472</td></tr><tr><th>Flowers</th><td>33.500</td><td>33.466</td><td>33.451</td><td>33.472</td><td>33.515</td><td>33.421</td><td>33.418</td><td>33.445</td><td>33.435</td><td>33.485</td><td>33.520</td><td>33.496</td><td>33.436</td><td>33.455</td><td>33.479</td></tr><tr><th>Peppers</th><td>33.413</td><td>33.357</td><td>33.466</td><td>33.460</td><td>33.423</td><td>33.519</td><td>33.406</td><td>33.398</td><td>33.429</td><td>33.457</td><td>33.429</td><td>33.492</td><td>33.489</td><td>33.527</td><td>33.440</td></tr><tr><th>Goldhill</th><td>33.461</td><td>33.442</td><td>33.435</td><td>33.468</td><td>33.363</td><td>33.499</td><td>33.464</td><td>33.406</td><td>33.492</td><td>33.495</td><td>33.466</td><td>33.429</td><td>33.380</td><td>33.460</td><td>33.500</td></tr><tr><th>Black</th><td>33.545</td><td></td><td></td><td>33.503</td><td></td><td></td><td>33.541</td><td></td><td></td><td>33.474</td><td></td><td></td><td>33.521</td><td></td><td></td></tr><tr><th>Boats</th><td>33.417</td><td></td><td></td><td>33.486</td><td></td><td></td><td>33.456</td><td></td><td></td><td>33.436</td><td></td><td></td><td>33.419</td><td></td><td></td></tr><tr><th>Clock</th><td>33.418</td><td></td><td></td><td>33.328</td><td></td><td></td><td>33.418</td><td></td><td></td><td>33.281</td><td></td><td></td><td>33.537</td><td></td><td></td></tr><tr><th>White</th><td>33.496</td><td></td><td></td><td>33.444</td><td></td><td></td><td>33.481</td><td></td><td></td><td>33.458</td><td></td><td></td><td>33.437</td><td></td><td></td></tr><tr><th>Average</th><td>33.449 (Pass)</td><td></td><td></td><td>33.452 (Pass)</td><td></td><td></td><td>33.459 (Pass)</td><td></td><td></td><td>33.458 (Pass)</td><td></td><td></td><td>33.474 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Sensitivity on ω<sub>0</sub><sup>(2)</sup></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th>Channel</th><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>Baboon</th><td>33.455</td><td>33.417</td><td>33.542</td><td>33.451</td><td>33.407</td><td>33.397</td><td>33.431</td><td>33.423</td><td>33.399</td><td>33.501</td><td>33.520</td><td>33.486</td><td>33.412</td><td>33.492</td><td>33.447</td></tr><tr><th>Flowers</th><td>33.398</td><td>33.522</td><td>33.477</td><td>33.501</td><td>33.491</td><td>33.475</td><td>33.513</td><td>33.402</td><td>33.525</td><td>33.451</td><td>33.387</td><td>33.398</td><td>33.446</td><td>33.481</td><td>33.469</td></tr><tr><th>Peppers</th><td>33.424</td><td>33.460</td><td>33.489</td><td>33.406</td><td>33.418</td><td>33.492</td><td>33.466</td><td>33.543</td><td>33.432</td><td>33.551</td><td>33.428</td><td>33.564</td><td>33.460</td><td>33.470</td><td>33.463</td></tr><tr><th>Goldhill</th><td>33.471</td><td>33.450</td><td>33.490</td><td>33.452</td><td>33.443</td><td>33.406</td><td>33.467</td><td>33.436</td><td>33.448</td><td>33.479</td><td>33.435</td><td>33.413</td><td>33.483</td><td>33.485</td><td>33.517</td></tr><tr><th>Black</th><td>33.460</td><td></td><td></td><td>33.469</td><td></td><td></td><td>33.484</td><td></td><td></td><td>33.510</td><td></td><td></td><td>33.409</td><td></td><td></td></tr><tr><th>Boats</th><td>33.446</td><td></td><td></td><td>33.461</td><td></td><td></td><td>33.451</td><td></td><td></td><td>33.472</td><td></td><td></td><td>33.505</td><td></td><td></td></tr><tr><th>Clock</th><td>33.570</td><td></td><td></td><td>33.325</td><td></td><td></td><td>33.557</td><td></td><td></td><td>33.473</td><td></td><td></td><td>33.364</td><td></td><td></td></tr><tr><th>White</th><td>33.425</td><td></td><td></td><td>33.438</td><td></td><td></td><td>33.530</td><td></td><td></td><td>33.503</td><td></td><td></td><td>33.450</td><td></td><td></td></tr><tr><th>Average</th><td>33.469 (Pass)</td><td></td><td></td><td>33.440 (Pass)</td><td></td><td></td><td>33.469 (Pass)</td><td></td><td></td><td>33.473 (Pass)</td><td></td><td></td><td>33.460 (Pass)</td><td></td><td></td></tr></tbody></table>

<table><caption>Table 20. Exemplar simulation for Structure 2 in Fig. 1b: Sensitivity of secret key by UACI.</caption><caption>Sensitivity on S<sub>0</sub></caption><thead><tr><th rowspan="2">N<sub>e</sub></th><th colspan="3">1</th><th colspan="3">2</th><th colspan="3">3</th><th colspan="3">4</th><th colspan="3">5</th></tr><tr><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th><th>R</th><th>G</th><th>B</th></tr></thead><tbody><tr><th>Channel</th><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>Baboon</th><td>33.521</td><td>33.452</td><td>33.544</td><td>33.473</td><td>33.478</td><td>33.452</td><td>33.518</td><td>33.519</td><td>33.416</td><td>33.527</td><td>33.499</td><td>33.463</td><td>33.491</td><td>33.498</td><td>33.395</td></tr><tr><th>Flowers</th><td>33.468</td><td>33.430</td><td>33.471</td><td>33.483</td><td>33.528</td><td>33.482</td><td>33.421</td><td>33.434</td><td>33.533</td><td>33.491</td><td>33.437</td><td>33.475</td><td>33.439</td><td>33.388</td><td>33.569</td></tr><tr><th>Peppers</th><td>33.499</td><td>33.439</td><td>33.459</td><td>33.432</td><td>33.547</td><td>33.555</td><td>33.429</td><td>33.443</td><td>33.394</td><td>33.401</td><td>33.425</td><td>33.433</td><td>33.321</td><td>33.468</td><td>33.412</td></tr><tr><th>Goldhill</th><td>33.448</td><td>33.523</td><td>33.439</td><td>33.474</td><td>33.494</td><td>33.429</td><td>33.499</td><td>33.413</td><td>33.469</td><td>33.440</td><td>33.478</td><td>33.406</td><td>33.487</td><td>33.473</td><td>33.497</td></tr><tr><th>Black</th><td>33.452</td><td></td><td></td><td>33.479</td><td></td><td></td><td>33.520</td><td></td><td></td><td>33.485</td><td></td><td></td><td>33.519</td><td></td><td></td></tr><tr><th>Boats</th><td>33.481</td><td></td><td></td><td>33.449</td><td></td><td></td><td>33.500</td><td></td><td></td><td>33.406</td><td></td><td></td><td>33.399</td><td></td><td></td></tr><tr><th>Clock</th><td>33.502</td><td></td><td></td><td>33.628</td><td></td><td></td><td>33.463</td><td></td><td></td><td>33.514</td><td></td><td></td><td>33.451</td><td></td><td></td></tr><tr><th>White</th><td>33.458</td><td></td><td></td><td>33.500</td><td></td><td></td><td>33.517</td><td></td><td></td><td>33.516</td><td></td><td></td><td>33.520</td><td></td><td></td></tr><tr><th>Average</th><td>33.470 (Pass)</td><td></td><td></td><td>33.474 (Pass)</td><td></td><td></td><td>33.453 (Pass)</td><td></td><td></td><td>33.477 (Pass)</td><td></td><td></td><td>33.490 (Pass)</td><td></td><td></td></tr></tbody></table>