Next Article in Journal
Buck-Boost/Flyback Hybrid Converter for Solar Power System Applications
Previous Article in Journal
Deadlock-Free Planner for Occluded Intersections Using Estimated Visibility of Hidden Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

JPEG Images Encryption Scheme Using Elliptic Curves and A New S-Box Generated by Chaos

by
Erick Hernández-Díaz
1,
Héctor Pérez-Meana
1,*,
Víctor Silva-García
2 and
Rolando Flores-Carapia
2
1
SEPI ESIME Culhuacan, Instituto Politécnico Nacional, Mexico City 04400, Mexico
2
CIDETEC, Instituto Politécnico Nacional, Mexico City 07700, Mexico
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(4), 413; https://doi.org/10.3390/electronics10040413
Submission received: 26 December 2020 / Revised: 26 January 2021 / Accepted: 2 February 2021 / Published: 8 February 2021
(This article belongs to the Section Computer Science & Engineering)

Abstract

:
This paper proposes a new symmetric encryption system based on an elliptical curve and chaos, where the encryption is done in a single block and runs for 14 rounds. Here, the 15 encryption keys have the same size as the image and are generated using a solution point of a strong elliptic curve. Using a string of random numbers obtained with a logistic map, a permutation and its inverse are generated, which improve the encryption level and add diffusion to the cryptosystem. Another important contribution to this research is the generation of a substitution box with a non-linearity of 100, which strengthens the cryptosystem against differential and linear attacks that add confusion to the scheme. Moreover, the cryptographic properties of the proposed S-Box are compared with those of the S-Box of the Advanced Encryption Standard (AES) to ascertain that it is a suitable alternative that it is resistant to differential power analysis (DPA) attacks. To verify the robustness of proposed cryptosystem against cryptanalysis and the quality of the cipher text obtained, this system is subjected to different tests, such as entropy, correlation coefficient, χ 2 , Changing Pixel Rate (NPCR), and Unified Averaged Changing Intensity (UACI). The results are shown and compared with recently proposed systems.

1. Introduction

A large amount of sensitive information is transmitted through the internet, and such information can be easily accessed through public as well as private networks. The images, which represent an important amount of transmitted information, can be photographs, agreements, contracts, identification documents, account statements, or other kind of scanned documents, with a high intrinsic value. Any person with a mobile device and an internet connection will be able to send images through the internet and store them in physical devices or in the cloud. The security and confidentiality of such information have led to the development of several efficient cryptographic systems.
In some practical applications, involving scientific and engineering fields such as astronomy or medicine, the use of images without compression and loss information is commonly required. Among the image formats that meet these requirements, we have the Microsoft-designed Bitmap (BMP) format, which contains all the image information in a simple format operating in the spatial domain. It can be modified and easily edited because it can be debugged and viewed without special tools; besides that, it is used to display image files, including the color of each pixel, and although it was designed to be a Windows standard, it is currently supported by a variety of operating systems [1]. In several countries, such as Mexico, where this research was carried out, it is prohibited to lose information in documents with sensitive information that need to be encrypted [2]; the use of BMP or another lossless scheme is imperative when it is required to store digital documents. To encrypt the BMP images, reference [3] proposed an efficient cryptographic algorithm based on elliptic curves and chaos. However, the size of BMP format images is too large for many applications; then, it is necessary to use some compressing format. Thus, this paper proposes an encryption algorithm for use in JPEG (Joint Photographic Experts Group) format [4], which is widely used in a variety of electronic devices such as smartphones, computers, or digital cameras. This task is carried out in two stages: the baseline, which has no loss of information, and DCT (Discrete Cosine Transform), which considerably reduces the file size due to the loss of information, which is achieved by discarding redundant pixels using a quantization matrix, reducing in such way the image size and doing it in an easier manner than file transmission or storing [4]. It is important to mention that the picture quality will depend on the compression level; usually, a level from 95% to 0.70% will give excellent quality but low compression, while if the maximum level is applied, 0% to 20%, the resulting image quality will degrade and sometimes become indistinguishable. Moreover, most of the other image formats are based on the RGB color models (Red, Green, Blue), but the JPEG format employs the color model YUV, where “Y” means luminance while “U” and “V” means chrominance, which works under the principle that the human eye is unable to identify certain bands of colors [5]. Thus, although there is a loss of information, it is imperceptible to the human eye.
Protecting images with sensitive information transmitted through public communication channels is an important task that has been a topic of active research during the last several decades, leading to the development of several efficient cryptographic algorithms. To this end, recently, researchers have developed some efficient cryptosystems based on chaotic systems and/or in elliptic curves that avoid discrete logarithm attack, since it has been proven that they are robust against differential [6,7], linear [8], or even statistical cryptanalysis.
Several efficient cryptographic schemes have been proposed recently in the literature based on elliptic curves and chaos, such as those described below. In [9], the authors propose two image encryption schemes in which their encryption keys are generated using the logistic equation that is presented in Equation (13) and Henon maps. Then, the image encryption is performed using the solution points of an elliptic curve and applying bitwise XOR operations in three different stages. In [10], the authors propose a three-step scheme where first, it uses two chaotic maps known as LTM and TSM to permute the pixels of the input image; the input parameters of these are obtained using the SHA-512 of the input image, from which the encryption keys are also generated. Then, using an asymmetric cryptosystem known as Elliptic Curves-ElGamal, the image is encrypted. Finally, using a genetic algorithm based on DNA sequences, diffusion is applied to the previously obtained output. In [11], the authors propose an encryption scheme in two scenarios; both use an Elliptic Curves-Diffie-Hellman protocol and another key generated from the SHA-256 of the input image. Furthermore, in the encoding process, a triple chaotic STH map is used, and a digital signature is added to verify its origin before being decrypted. In [12], the authors propose a scheme that uses ECC to share a point of an elliptic curve that will serve to obtain the initial parameters ( x 0 ,   μ ) of a logistic map and the number of rounds of Arnold’s map that will serve to permute the pixels from an image. The parameter k of the elliptic curve used is obtained with the SHA-512 and using the multiplication of points between k and a generator, the encryption keys are created. In [13], the authors propose a four-step scheme. Here, firstly, they reduce the input image, which in this case is an RGB image, and extend it in the grayscale for the entire field E (F_p). Secondly, they apply an enhanced 4D Arnold cat map t times to permute the image pixels. Next, they encrypt the image with an ECC. Finally, they apply a 3D Lorenz chaotic map to add diffusion to the final encrypted image. Finally, in [14], the authors propose a new method called EC-GRP permutation operation, which combines the cryptographic properties of the ECC and the operation group. This procedure requires a pseudo-random bit generator with good attributes according to the cryptography field, which is used to encrypt the input image. According to the researchers, the results obtained are highly resistant to cryptanalysis tests.
Despite the efficient algorithms described above, several issues remain that must be improved. To this end, this paper proposes a cryptographic algorithm for improving the above-mentioned schemes, whose main contributions are summarized as follows: (a) It uses elliptic curves with a constant l equal to zero for the generation of the set of encryption keys. (b) It proposes an algorithm to generate the elliptic curves in addition to being mandatory that they comply with certain characteristics. (c) It uses a chaotic logistic equation to generate permutations and a substitution box plus its inverse with a non-linearity level of 100. (d) It implements a fingerprint for the receiver to identify if the received file corresponds to the one sent by the issuer. (e) It encrypts the images in a single block, which allows obtaining an adequate encryption speed. On the other hand, the cryptosystem was evaluated using several tests. After analyzing the evaluation of the results reported in Section 5, it follows that the proposed structure is robust and capable of withstanding linear, differential, statistical, brute force, or modification attacks, as well as some of the better known as the discrete logarithm and the MOV.
The rest of this paper is organized as follows. In Section 2, the preliminaries of this research are stated. Section 3 describes how the encryption keys, the permutations, the substitution box, and the complete encryption algorithm are developed. Section 4 provides several essential information such as the images that are encrypted, the explanation of every test necessary to demonstrate how robust is our proposal against cryptanalysis, and the results obtained. In Section 5, the analysis of the data obtained after applying all the tests is provided. Finally, Section 6 presents the conclusions of this research.

2. Preliminaries

This section provides a background and a detailed explanation about the two pillars on which this research is based: the elliptic curves and chaos.

2.1. Elliptic Curves in the Field of Cryptography

The elliptic curves have been used in many fields of science and engineering since their discovery. However, they were firstly introduced in the cryptography field in 1985 and 1987 by the mathematicians, Victor S. Miller in [15] and N. Koblitz [16]. An elliptic curve E is a projective geometric shape that is defined in a field F p in which a set of solutions # E ( F p ) is constructed that has two variables that satisfy the Weierstrass’ mathematical expression given by
y 2 x 3 + k x + l   m o d   p .
This set of solutions is Abelian, in which it is possible to define the addition operation ( E ,   + ) . In # E ( F p ) , a cyclic subset is built, where q is a prime factor that represents the number of elements (solutions) of this and the order of the elliptic curve. The prime factor q is calculated as follows:
q = p + 1 + 2 a 4   m o d   p .
The element of the subset that generates all other solutions is known as the generator and is denoted with the Greek letter α . The element ( q 1 ) α is the additive inverse of α ; that satisfies ( q 1 ) α = ( x 0 , y 0 ) . So, q α is the null element, which is written as infinity [17]. As it has already been established, with the addition of the other elements, it can be calculated—that is 2 α , 3 α . To calculate their coordinates ( x ,   y ) , it is necessary to obtain the slope of the line and there can be three cases, which are described below.
  • When x 0 x 1 , Equations (3)–(5) are used.
    λ = ( y 1 y 0 ) ( x 1 x 0 ) 1   m o d   p ,
    x 2 = ( y 2 x 0 x 1 )   m o d   p ,
    y 2 = λ ( x 0 x 2 ) y 0   m o d   p .
  • When x 0 = x 1 and, y 0 = y 1 , Equations (3)–(5) are used.
    λ = ( 3 x 0 2 + a ) ( 2 y 0 ) 1   m o d   p ,
    x 2 = ( λ 2 2 x 0 )   m o d   p .
  • When x 0 = x 1 and y 0 = y 1 , the null element is obtained, such that ( x 0 ,   y 0 ) + ( x 0 ,   y 0 ) = , it means that both are inverses concerning the elliptic curve addition operation. Furthermore, α + = α .
In this research, we use the elliptic curves that satisfy Equation (8), whose palpable difference with (1) is that the constant l is equal to zero; thus, for practical reasons, it is omitted in the mathematical expression,
y 2 x 3 k x   m o d   p .
The elements k ,   p ,   q are prime numbers; to ensure that they have this characteristic, it is necessary to apply a primality test such as Miller Rabin’s [18]. To calculate k, Equation (9) is used in addition to complying with the characteristics described in Theorem 1 [19], which is also useful to explain how # E ( F p ) is calculated, that is
k = ( x 0 3 y 0 2 ) ( x 0 ) 1   m o d   p .
Theorem 1.
Let p = a 2 + b 2 and p 1 mod 4 ; where p is a prime number, a a positive odd integer, b a positive even integer, and # E ( F p ) the number of solutions of the elliptic curve shown in (8). Then, the number of solutions is # E ( F p ) =   p +   1 +   2 a , provided that k is not the fourth power modulo p of any element of the field F p ; however, k must be the power squared of some element of the field F p [19].
In addition to having this expression, they must meet four requirements:
  • Being Non-Singular [17]. This happens when an elliptic curve meets the condition: 4 ( k ) 3 mod p   0 .
  • To be Non-Supersingular [17]. Non-Supersingular elliptic curves resist the attack designed by Menezes, Okamoto, and Vanstone [20]. They fulfill the following condition: q mod p   1 .
  • Not be an Elliptic Curve of Trace One [17]. The elliptic curves of Trace One are weak. If they do not belong to this type, they must satisfy that p q .
  • The size of q is at least 2 160 bits [17]. With this characteristic in the elliptic curves, it is impossible to solve the problem of the discrete logarithm that is capable of finding m with Q = m P when P and Q are known [6]. Even though the existence of the Pohlig–Hellman algorithm can solve this problem [21], it requires a high computational cost and exponential processing time.
There are research works that show that the elliptic curves such as those described in this paper are robust and adequate for the development of safe systems [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. Figure 1 presents the complete algorithm for the search of the elliptic curves to be used in this research.

2.2. Logistic Map

The Australian mathematician Robert May [18] proposed a logistic map to explain the abnormal behavior in the growth of populations; however, in the field of Cryptography, it is very useful to generate strings of pseudo-random numbers. This equation is obtained from the discrete expression shown in Equation (10), where e , f > 0 .
d P d t = e P f P 2 ,
After this, the variable P takes discrete values in time such that P ( t 0 ) ,   P ( t 1 ) , P ( t 2 ) ,   P ( t n ) . Afterward, Euler’s Algorithm is applied so that the result obtained is shown in Equation (11) [23].
P n + 1 = P n + ( e P n f P n 2 ) h .
The next step is to rewrite the equation as shown in Equation (12), where s = 1 + e h and, t = f h .
P n + 1 = s P n t P n 2 .
Lastly, the iterative equation shown in Equation (13) is obtained. Where 0 <   x n < 1 and 0 <   s < 4 [24],
x n + 1 = s x n ( 1 x n ) ,
where 0 <   x n < 1 and 0 <   s < 4 [24]. Two behaviors can be observed when iterating the logistic equation after 1000 cycles: a limit with a tendency to infinity ( x = lim n X n ), in which case it is stated that chaos occurs, or stabilization of the digits at the right of the decimal point, which implies that chaos does not take place, since there is a pattern to follow. Both are exemplified in Table 1.
This logistic equation has three main characteristics:
  • It is deterministic; if the initial parameters are the same, the same string will be obtained.
  • It is highly sensitive to changes, as seen in Table 1; by varying one of the parameters a little, the chain obtained is different.
  • The output string is impossible to predict, since it is not the result of any algebraic equation with rational coefficients, unless the input parameters are known beforehand, in which case it might be found.

3. Development

In this section, we provide a complete description of how the proposed cryptosystem is developed. Firstly, the algorithm used to generate the rounds keys is explained, which is followed by the way to get permutations and the proposed S-Box, and finally, all the previous pieces are joined to explain the complete encryption scheme.

3.1. Encryption Keys

The security of any cryptosystem relies on the encryption key; anyone with access to it will be able to decrypt the information sent. In this study, a private key cryptosystem (also called symmetric) is proposed in which the main key is private; that is, it will be exclusively known to the sender and receiver. From this one and using a keys schedule, all the round keys are generated. In this type of scheme, the set of encryption keys are the same for the encryption and the decryption process, but to decrypt, they are employed strictly in reverse order.
The proposed symmetric cryptosystem generates the main key K from a point of an elliptic curve that meets the characteristics described in section II and that must be known to all the people involved in the secure communication scheme. In total, 15 round keys are required.
To begin with, the sender and the receiver will choose an integer r , which is the private key, and which will refer to a solution point of an elliptic curve; also, it will have to fulfill with 1 < r < p 1 , and its size in the binary expression must be equal to 256 bits. From r, a new subset solution will be constructed, that is, r α = ( x 0 ,   y 0 )     α = ( x 0 ,   y 0 ) . Starting from , all the coordinates x ,   y of every point will be concatenated to get a string of the same size as the input image ( m   ×   n ) ; if its length is overcome, it will be adjusted suppressing the surplus bits. For example, if the input image has the dimensions 512   ×   512 , that is, 262,144 pixels, the length of the array K must be 262,144 bits after any adaptation.
K =   α   | |   2 α   | |   3 α   | |     | |   n α     K =   x 0   | |   y 0   | |   x 1   | |   y 1   | |     x 2   | |   y 2   | |     | |   x n   y n
Once K has been obtained, the procedure to obtain the first round key k 1 will be to permute the bits of each key with the S-Box that is presented in Section 3.3 to increase its level of randomness; then, at the output obtained, a 5-bit shift is applied to the right, and this process will be repeated to obtain k 2 but this time starting from k 1 , and so on until leading to k 15 whose origin will be k 14 .
k 1 = ( S B O X ( K ) ) 5 ;
k 2 = ( S B O X ( K 1 ) ) 5 ;
.
k 15 = ( S B O X ( K 14 ) ) 5
Since each coordinate of a point is considered as pseudo-random, K is considered as a pseudo-random string, too. The only way an attacker could infer K is to know the value r and the elliptic curve chosen to start the concatenation of values, which can be solved by protecting the process of sending this information through insecure means by encrypting it with some asymmetric cryptosystem. In addition, it would face the problem that an elliptic curve similar to the one proposed in this research and presented in Section 4.2 can have a solution set with an extension of at least 2 160 [12].

3.2. Permutations and the Proposed Substitution Box

One of the objectives of any cryptosystem is to obtain an encrypted text with a high level of randomness; this characteristic ensures that a scheme is resistant to differential and linear cryptanalysis. Having said the above, two techniques help to strengthen any encryption system: these are diffusion and confusion [25]. Diffusion disperses the elements of the plain text so that one may hide the relationship that exists between it and the ciphertext obtained at the output; to achieve it, permutations are used. The confusion confuses an attacker and makes it difficult for him to establish a relationship between the ciphertext and the key. Substitution boxes (S-Box) are used to achieve this. Both are used in symmetric cryptosystems that encrypt by blocks such as DES, Triple-DES, or AES, the proposal presented in this research work falls into that category, which is why they are also included.
In the proposed scheme, different permutations are performed in each execution. To build them requires a pseudo-random number obtained from Equation (13) in which it is proven that chaos has occurred. It is important to mention that for this task, only the digits to the right of the decimal point are taken into account. Then, a bijective function described in [22] is used; this algorithm adjusts to construct a permutation of N elements depending on the dimensions m   ×   n of the input image.
The same algorithm is used to build the S-Box; however, their design is a more complex task, since it must satisfy various criteria to be considered safe and resistant to DPA (differential power analysis) attacks that use statistical techniques to obtain information that helps to infer the encryption key of a cryptosystem and are usually more efficient than linear and differential attacks [26].
Perhaps the most important parameter to consider is non-linearity, which is the number of bits that must be modified in the truth table of the Boolean function to get closer to the closest affine function [27]. The non-linearity is represented in G F ( 2 N ) , that is, N =   2 m 1 2 m 2 1 , for m = 8 ; thus, in theory, the upper bound of N is 120 [28]. The S-Box proposed in this research work has a non-linearity of 100 equivalent to 78%, whose value is within the range of expected parameters [29], and it is measured using the Walsh function [22], which is presented in Equation (15).
L N = min N L i ,
where
N L i = 2 7 1 2 m a x a 2 8 | W f i ( a ) | .
Table 2 shows a complete list of the cryptographic properties of the proposed S-Box, and they are compared with those of AES. Although our proposal only exceeds the second S-Box in two parameters—Robustness to Differential Cryptanalysis and Transparency order—the others are found within the ranges of expected values for which it is considered adequate for the research being presented and robust enough against DPA attacks.
The input parameters used to build the substitution box are listed below. In addition, the proposed S-Box and its inverse are shown in Figure 2.

3.3. Complete Cryptosystem

This cryptosystem is symmetric; it is of the Substitution Permutation Network type, since the image encryption process is carried out through several encryption rounds, using XOR bitwise operations, a permutation, and the substitution box proposed in Section 3.2 [25]. Symmetric cryptosystems usually perform block encryption as with Triple-DES or AES [30,31], whose characteristic consists of dividing the information into blocks of fixed size (64 or 128). However, in this proposal, the encryption is carried out in one whose size is determined by the dimensions of the image, as already explained in Section 3.1. It was also explained that using a solution point of an Elliptic Curve, 15 encryption keys are generated, which will be applied one at a time for each of the 14 rounds of this proposal, except for the last one, which needs two. Starting from the fact that rounds 1 and 14 are the only different ones and rounds 2 to 13 are identical, we proceed to explain the complete operation of the algorithm. To describe the encryption and decryption process, the following nomenclature is used:
  • I for the input image.
  • I C for the encrypted image.
  • I D for the decrypted image.
  • R A ,   B , for each round, where A is equal to the round number (1–14), and, B for each intermediate operation having a maximum of 3.
  • O A , which denotes the output obtained after round A .
  • k 1 to k 15 for every encryption key from one to 15.
  • P   ( R A ,   B ) and P 1   ( R A ,   B ) for the permutations that modify the previous output.
  • S B O X   ( R A ,   B ) and S B O X I N V   ( R A ,   B ) for the replacement box that modifies the previous output.

3.3.1. Encryption Process

Step 1. Reading and Modifying the input image. An input image I of dimensions m   ×   n is read. If m = n = 512 , the dimension m is made to grow by one whereby m = 513 and n = 512 . Then, the pixel values of the new row from position 0 to 255 will match with every bit of the integer r used to generate the main key; that is, if r = 1110 0101 , pp (pixel in position) 0 = 1, pp 1 = 1, pp 2 = 1, pp 3 = 0, …, pp 252 = 0, pp 253 = 1, pp 254 = 0, pp 255 = 1. It is important to mention that these additional 256 pixels can serve as a fingerprint that verifies the origin of the received image in the decryption process. The pixels from position 256 to the last one can be filled with pseudo-random values.
Step 2. Round one of encryption. The operations carried out are summarized as follows:
  • R 1 ,   1 = I     k 1
  • R 1 ,   2 = S B O X   ( R 1 ,   1 )
  • O 1 = P   ( R 1 ,   1 ) .
Step 3. Rounds two through thirteen of encryption. The operations carried out are summarized as follows, and since all these rounds are the same, only two are described:
  • R 2 ,   1 = O 1     k 2
  • O 2 = S B O X   ( R 2 ,   1 ) .
Step 4. Round fourteen of encryption. The operations carried out are summarized as follows:
  • R 14 ,   1 = O 13     k 14
  • R 14 ,   2 = S B O X   ( R 14 ,   1 )
  • R 14 ,   3 = P 1   ( R 14 ,   2 )
  • I C = R 14 ,   3     k 15 .

3.3.2. Decryption Process

Step 1. Round one of decryption. IC is read, and the operations carried out are summarized as follows:
  • R 1 ,   1 = I C     k 15
  • R 1 ,   2 = P   ( R 1 ,   1 )
  • R 1 ,   3 = S B O X I N V   ( R 1 ,   2 )
  • O 1 = R 1 ,   3     k 14 .
Step 2. Rounds two through thirteen of decryption. The operations carried out are summarized as follows, and since all these rounds are the same, only two are described:
  • R 2 ,   1 = S B O X I N V   ( O 1 )
  • O 2 = R 2 ,   1     R 2 ,   1 .
Step 3. Round fourteen of decryption. The operations carried out are summarized as follows:
  • R 14 ,   1 = P 1   ( O 14 )
  • R 14 ,   2 = S B O X I N V   ( R 14 ,   1 )
  • I D = R 14 ,   2     k 15 .
Step 4. Reading and Modifying the output image. I D is read; thus, the last row of pixels is deleted; in this case, the final dimensions of I D will be m = 512 and n = 512 . Finally, it is easy to conclude that I D = I .

4. Experiments and Security Analysis

In this section, we carry out all the required experiments to prove if this cryptosystem is resistant against the main types of cryptanalysis and modifications. Another piece of essential important information is provided, too. To facilitate understanding, graphs and tables are added.

4.1. Images for Experiments

The proposed cryptosystem has undergone several tests to demonstrate its robustness to the differential, lineal, and differential attacks. To carry out each experiment, JPEG images of different sizes have been chosen; these are:
  • Lena.jpeg, Baboon.jpeg, and Boat.jpeg, with dimensions of 512 × 512 pixels.
  • Barbara.jpeg, with dimensions of 720 × 576.
  • City.jpeg, with dimensions of 1280 × 720.
  • House.jpeg, with dimensions of 1920 × 1080.
  • Security.jpeg, with dimensions of 4900 × 3464 pixels.
The first four images are commonly used in peer-reviewed papers, and the others are proposed in this research. The desktop application in which the proposed algorithm was implemented in Java programming language and the BufferedImage library was used to support certain tasks such as the reading of the pixels of each image [32]. All pictures are shown in Figure 3.

4.2. Elliptic Curve to Generate the Encryption Keys

To perform all the experiments, the elliptic curve has been chosen with the following data:
a = 1bdb5c33c01c799169c62dae71c4176c5d1b0d
b = 1b1dbcab63715d81db85c8db589685e3b5f1ec
k = 59196553ed14991dbb65913217a9e4cffd676c103fa9df60df6b4b1681806d3571a8bfba36a
p = 5e74aaa3d208d1c5322ef4c2327fa7e23143aa8eab32b38e0dd1705f83c7cd9909626a3f039
q = 179d2aa8f48234714c8bbd308c9fe9f88c50eb8285ae4ae46740e766125f66f46313fd78995
r α = 8e28e3a5cef3c15189a8b1d9a8d87a2ffc9c79146fde29b7984fd997ea4b2fbc
x 0 =   742fae6e43b6b33e4d4b49672c2097c87bbcbdba4f7b28b8969e844f19f3941acf646d88ae
y 0 =   55546d0927fc1d205516fbbe4ceb24ed2495e01adc421a814a5b404167a3674e8f2c15d9bf
Next, fast analysis of the chosen elliptic curve is conducted:
  • If the calculation 4 ( k ) 3 mod p is performed, the result obtained will be: 5a67ecc091bc618e74d07f4905361dd6c7e67fc44d6d26c76fbfbd9625a8b43a4eddb4fee4a; therefore, it is shown that this is a Non-Singular Curve.
  • Later, if the calculation q mod p is performed, the result obtained will be: 179d2aa8f48234714c8bbd308c9fe9f88c50eb8285ae4ae46740e766125f66f46313fd78995; therefore, it is verified that this elliptic curve is not Supersingular.
On the other hand, it is easy to observe that p and q are different; in this way, it is stated that this elliptic curve is not of Trace One.
  • Finally, this elliptic curve has a solution set q of a size of 2 256 .
Thus, it is concluded that the chosen elliptic curve fulfills the four requirements described in Section 2.1 and is safe and suitable to generate the set of round keys.

4.3. Encrypted Images

Figure 4 shows the encrypted results of images shown in Figure 3. They prove that it is visibly impossible to find a pattern that allows inferring its origin. Nevertheless, it is important to understand that a visual inspection is not enough to demonstrate the cipher text is impossible to reverse; for that reason, in the next sections, several tests will be applied to ascertain it mathematically.

4.4. Statistical Cryptanalysis

An important aspect in any cryptosystem is to quantify the quality of the encryption, that is, the level of randomness, which will determine its resistance to statistical attacks that can determine the encryption key or the plaintext through existing biases or patterns in the ciphertext. Statistical tests such as those presented below are useful for this purpose.

4.4.1. Entropy

Entropy is one of the most important works of the French mathematician Claude E. Shannon [33]. In an image, it must be applied in each band of this depending on its type of color space. This test analyzes the histogram of the figure after it has been coded and determines if its frequency distribution is much more uniform than that of the original image. It is known that 255 is the maximum value of a pixel, and for its binary expression, it requires 8 bits; this implies that the perfect distribution of an image after being encrypted must be 8, which is unlikely to happen in practice. Therefore, any value greater than 7.9 indicates a high entropy level [15]. Equation (17) is used to calculate entropy.
H ( x ) = i = 0 255 P ( x i ) l o g 2 P ( x i )
The results obtained in this test are shown in Table 3. The evaluation results show that the entropy provided by the proposed scheme is quite close to 8, which is the maximum theoretical value. Moreover, Figure 5 shows the histograms per every RGB channel of the Lena image used in the experiments to visualize its distribution before and after encryption. Here, it can be shown that after encryption, the estimated histogram is almost flat, independently of the histogram shape of the image before encryption. These types of graphs are widely used to carry out a visual inspection in this test.

4.4.2. Correlation Coefficient

The correlation coefficient is based on the problem of analyzing the relationship between two variables x ,   y [34]. In the specific case of image encryption, it analyzes whether the position of two contiguous pixels x ,   y is determined by some given pattern or if there is a dependency between two pixels in each image. Otherwise, it is said that there is no correlation, and both are randomly positioned. This measure is obtained with Equation (18).
r = [ ( x ( x ) ) ( y ( y ) ] i n i = 1 n ( x i ( x ) ) 2 i n i = 1 n ( y i ( y ) ) 2
If the measurement yields a 0, it is said that both images are completely different, and if a 1 or −1 is obtained, the conclusion would be that both images are equal. In practical terms, it is almost impossible to get zero, but anything value close to this indicates a high level of randomness between the pixels of an encrypted image. The most appropriate way to carry out this measurement is in three directions: Horizontal, that is, a pixel x and its neighbour to the right; Vertical, that is, a pixel x and its neighbor below; and Diagonal, that is, a pixel x and its neighbor to below shifted by one space to the right.
Table 4, Table 5 and Table 6 show the correlation coefficients obtained; only absolute values are registered. The evaluation results show that the values of horizontal, vertical, and diagonal cross-correlation of the images under analysis closely approach zero, which means that it is not possible to infer one pixel of the image under analysis. Thus, if the correlation between the pixels of the encrypted image approach to zero, it is not possible to infer the original image using only information of the encrypted one. In addition to the experimental data given in Table 2, Table 3, Table 4, Table 5 and Table 6, Figure 6 shows the scatter plot per every RGB channel of one of the images used in the experiments to visualize the position of all the pixels before and after encryption; these types of graphs are widely used to carry out a visual inspection in this test. From this figure, it follows that in the encrypted image, it is not possible to estimate the value of pixel (x, y) using the value of pixel on the other locations of the same image.

4.4.3. χ 2 Test

This statistical test is based on the statement of two hypotheses: the first of them is known as null and is identified as H 0 [35]; in this case, its role is to affirm that the pixels of an image after being encrypted have a random distribution. The second one is known as the alternative and is called H 1 ; its function is to reject the assertion proposed by the null. Using Equation (19), χ 2 is calculated, and it determines which of the two hypotheses is accepted and which is rejected. f o refers to the observed frequencies, that is, those of the encrypted image, and f e talks about the number of expected frequencies, which in this case is 256.
χ 2 = i = 1 n ( f o i f e ) 2 f e .
It is known that in tests based on hypotheses, there are two errors: a type I error that is the most important, that is, unequivocally rejecting, H 0 ; and, type II, that is, accepting H 0 wrongly. For this research work, the type I error is used, whose value is α = 0.01 . In practice, any threshold ( χ 2 ) less than 308 is enough for H 0 to be accepted [36]. Table 7 shows the results obtained from the test after the experiments.

4.5. Differential and Linear Cryptanalysis

Cryptanalysis is the antagonistic science of cryptography whose objective is to develop attacks capable of compromising or breaking encryption algorithms. Symmetric and block encryption cryptosystems must be robust to at least two types of cryptanalysis, differential and linear [37].
The differential attacks were first proposed by Eli Biham and Adi Shamir with the aim of breaking the DES cryptosystem [7]. Attacks of this type exploit the high probability of the existence of differences in the plain text ∆X and differences in the last round of encryption ∆Y, which is known as the differential. This attack selects inputs and analyzes outputs to find the encryption key.
The linear attack was proposed for the first time by Mitsuru Matsui, who sought to exploit the DES cryptosystem through the known plaintext at entry [6]. This attack works considering the linear correlations between some of the plaintext bits (input block) and the output bits (cipher block) to infer the cipher key. There are various tests that a cryptosystem must undergo to verify that it is not vulnerable to the attacks described above; this process is carried out below.

4.5.1. NPCR and UACI

Both standards serve to test the resistance of any cryptosystem against differential cryptanalysis; they function as follows. Starting from the fact that there are two encrypted images, C 1 and C 2 , which come from two images whose only difference is a pixel, and the encryption process has been carried out with the same keys. If the proposed cryptosystem is robust, images C 1 and C 2 must be practically different, which can be measured with the NPCR and UACI standards [38]. The first is defined by Equation (20) and the second is defined by (21) where, C 1 and C 2 refer to the images, T refers to the number of pixels in each image ( n   ×   m   × number of planes), and D is defined in Equation (22).
N P C R : N ( C 1 ,   C 2 ) = i , j D ( i , j ) T   ×   100 % ,
U A C I : U ( C 1 ,   C 2 ) = i , j | C 1 ( i , j ) C 2 ( i , j ) | 255     T   ×   100 % ,  
D ( i , j ) = { 0 ,   i f   C 1 ( i , j ) = C 2 ( i , j ) ,   1 ,   i f   C 1 ( i , j ) C 2 ( i , j ) .
In the case of this experiment, both images are encrypted with the set of encryption keys, which are generated from the K corresponding to the SHA-1 of both. In all the images used, pixel 3750 of the blue channel has been modified. In practical terms, in these tests, it is expected to obtain percentages between the range of 99.5% and 99.6% for the NPCR and between 33.4% and 33.5% for the UACI. The results obtained in this test are shown in Table 8.

4.5.2. Avalanche Effect

The avalanche effect, also known as avalanche attack, is another important manner to test the resistance of any cryptosystem against differential attacks [39]. This standard works with the same principle that was observed in NPCR and UACI tests, where a tiny change made in an image will produce a practically different ciphered image, although in this case, the change is made on a bit level. Having said that, let I 1 and I 2 be two plain images with only one different bit; then, they are ciphered using a group of round keys generated from a main key K with just a distinct bit among each other. Thus, the resulting images C 1 and C 2 must be shown a bit rate of changing approached to 50% [40]. Using Equations (22) and (23), it is possible to obtain such a measure.
A v a l a n c h e = i , j D ( i , j ) T o t a l   b i t s   ×   100 %
The results obtained in the avalanche test are given in Table 9.
The evaluation results show that the avalanche effect resulting from the proposed scheme is quite close to the ideal value, which is equal to 50%.

4.5.3. Chosen/Known Plain-Text Attacks

These types of attacks fit into the category of linear cryptanalysis, and there four of them in total, which are the known-plaintext attack, the chosen-plaintext attack, the ciphertext-only attack, and the chosen-ciphertext attack; the first two are the most important, and it is stated that any cryptosystem capable of supporting them will do the same with the last two [39]. The procedure to test if a cryptosystem can resist them consists of encrypting two images: one white and one black, and then measuring their entropies and correlation coefficients and verifying if they fall within the parameters described in Section 4.4.1 and Section 4.4.2. In this experiment, five images of each are used, corresponding to the dimensions of those chosen for the experimentation stage. Table 10 shows the average results obtained in this test.

4.6. Key Sensitivity Test

The objective of this test is to analyze the percentage of different pixels between two equal images, C 1 and C 2 . Both will be encrypted with a set of encryption keys generated from a specific K , and in the case of the second image, the same K will be used as for C 1 but with a small modification that could be the change of one bit. If we want to decrypt C 1 with the keys of C 2 , it must be impossible, and vice versa. The percentage of change is calculated with Equation (24), and in practical terms, it is expected that minimum values of DiffImg equal to 99% will be obtained. Table 11 shows the results obtained from the test after the experiments.
D i f f I m g =   D i f f P i x e l s T o t a l P i x e l s   ×   100

4.7. Keyspace Analysis

The brute force attack on an encryption key is inescapable, where its objective is to try all possible combinations until finding the one that matches. In this case, the keys are generated from a solution point of an elliptic curve, which implies that there are q different possibilities, because q is a prime number. Furthermore, for a key to be secure, it must have a minimum size of 2 100 bits. Since the keys depend on the integer r which is 256 bits, the size of the keys used in this proposal is about 2 256 , which is far beyond the minimum required.

4.8. Occlusion and Noise Attacks

When an encrypted image is transmitted over an insecure medium, it takes the risk of being intentionally or incidentally modified or distorted by an attacker. Any clipping, obstruction, or change in pixels that can be interpreted as added noise implies information that has been lost and is impossible to recover. Nevertheless, if the image encryption quality of the image is high and all the pixels have been perfectly distributed in a way considered as random, when the image is decrypted, much of the original information will still be possible to display. Therefore, it is very important to measure the resistance of any cryptosystem, mainly to two attacks: that is the occlusion and the added noise.
Furthermore, in this test, an image that displays a written message is used, which is shown in Figure 7. Sometimes, these types of pictures are sent to share notes or news and, commonly, attackers can try to interfere or damage them to avoid that the communication can be completed.
Regarding the first attack, Lena.jpeg and Text.jpeg are encrypted, and then, various sections of an amplitude of 25%, 50%, and 75% of the total surface are intentionally cut out. Thus, when the image is deciphered, it is visually analyzed if the pixels remaining rearranged are enough to be able to infer which was the original image. The second attack consists of adding salt and pepper noise to the encrypted image, which in practice consists of scattering black or white pixels pseudo-randomly on the surface; for this case, the noise is added in densities of 25%, 50%, and 75%, and then, the same visual analysis is done as in the first one. The obtained results are shown in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13.

4.9. Time Encryption

The computer used for all the tests had the next hardware resources:
  • Processor: Intel Core 13 7350k, 4.00 GHz, dual-core.
  • RAM Memory: 8 GB, 1600 Mhz, DDR3L.
  • Video Card: Intel 630.
  • Hard disk: SATA 7200 RPM.
  • Operative System: Windows 10 Pro.
Having explained this, the speed encryption was measured five times for every image; thus, in Table 12 shows the average speed calculated.

4.10. Comparison with Other Articles of the State of the Art

Table 13 shows the average results obtained in this research in four tests for entropy and correlation coefficient, while Table 14 shows a comparison of the NPCR, UACI, and the avalanche obtained by proposed scheme compared with those obtained using other similar papers recently published in the literature.
Table 13 and Table 14 shows that the evaluation results provided by the proposed are quite competitive with other previously proposed schemes.

5. Analysis and Results

This section provides an analysis of the evaluation results obtained by every test performed in a strictly sequential order as they were performed. The first evaluated parameter is the Entropy; according to the results of Section 4.4.1, this test seeks to find a value that means that the encrypted image has almost a uniform frequency distribution. If a value greater than 7.9 is obtained, the entropy is high and complies with the previous approach. It can be seen from Table 3 that the lowest calculated value is 7.9992, which is very close to the ideal value for which an image encoded with 8 bit/pixel is equal to 8. The next evaluated parameter is the Correlation coefficient. In Section 4.4.2, it was explained that this test seeks to measure the level of dependence between the contiguous pixels of an image encrypted in three directions: horizontal, vertical, and diagonal. There is a high non-linearity when the values approach zero, which happens in all cases. It should be remembered that the results of Table 4, Table 5 and Table 6 are expressed in absolute values. These tables show that the cross-correlation values of encrypted signals approach zero, which means that the knowledge of some pixel values does not allow the estimate the encrypted images. The χ 2 Test proposed in Section 4.4.3 is based on the proposition of two hypotheses; the first indicates that the encryption carried out is random, and the second contradicts it; the acceptance of each one depends on a threshold ( χ 2 ) that must be less than 308. It is observed that all the values recorded in Table 7 are less than 300 in all cases.
The resistance to differential and linear cryptoanalysis was also evaluated. To this end, the NPCR and UACI were estimated. As mentioned in Section 4.5.1, this test seeks to measure using both standards of the quantity and percentage of different pixels between two images whose only difference is one pixel. It was also explained that for the NPCR, minimum values of 99.5% are expected, and for the UACI, we expect values not less than 33.4%; as seen in Table 8, the lowest value of the first case is 99.60%, and for the second, it is 33.44%. Other test is the avalanche test, in which is evaluated that a tiny change made in an image will produce a practically different ciphered image. The evaluation results shown in Table 9 show that the obtained results are very close to 50%, which is the ideal value. The Chosen/Known plain-text attacks are the other recommended evaluation. In this test, black and white images of different dimensions are coded, and subsequently, their entropies and correlation coefficients are measured, which are recorded in Table 10. All measurements obtained are within the ranges specified in Section 4.5.1 and Section 4.5.2.
The Keyspace and key sensitivity was also evaluated. In the case of the first, it was already explained in Section 4.7 that the minimum size for a key to be considered secure is 2 100 , which is far exceeded, since in our proposal, it is 2 256 . On the other hand, the second test measures the percentage of distinct pixels when an image is encrypted with two keys with minimal variation. Table 11 shows percentages greater than 99% that correspond to what is expected in practice. Another evaluated attack is the Occlusion and noise attacks. The objective of both is to determine with a visual inspection how much information is recovered from an encrypted image that is covered or cropped to a certain extent as in the occlusion attack, or when it is randomly covered with white or black pixels after adding salt and pepper noise. For both cases, it is carried out in extensions or intensities of 25%, 50%, and 75%. Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 show that although the level of visibility after decryption is lower with higher intensities, it is possible to perceive the message. Following the foregoing, it is determined that the results obtained are good and all the tests have been satisfactorily passed.
Finally, Table 13 and Table 14 compare the results of the entropy, correlation coefficient, NPCI, UACI, and avalanche tests of the proposal presented with other articles of the state of the art. It is easy to see that the results obtained are close to those reported by other recently proposed papers. In some of the cases, they are equal or slightly lower, and in others, better measurements are obtained; therefore, it can be said that the system proposed in this research is competitive concerning the existing state of the art.
Finally, the future works that can be carried out for enhancing the performance of this research work are mentioned as follows. Researchers could develop substitution boxes with a higher level of non-linearity and better cryptography properties or use other logistic equations such as Lawrence’s equation. The authors consider that it would also be of great value to verify the effectiveness of this proposal or others in the encryption of audio or video. Finally, the introduction of a steganography scheme could increase the security of the system and propose a scheme that includes a key distribution system.

6. Conclusions

In this article, an original symmetric cryptosystem has been presented whose purpose is the encryption of images in JPEG format, which stands out for its high quality, its lightness, and its ease of transmission, which makes it one of the most used at present. The contributions of this work that make it different from others that were mentioned in the state of the art are summarized as follows. Firstly, the use of elliptic curves with a constant l equal to zero for the generation of the set of encryption keys. Next, an algorithm is proposed to generate the elliptic curves in addition to being mandatory that they comply with certain characteristics. The proposed scheme uses a chaotic logistic equation to generate permutations and a substitution box plus its inverse with a non-linearity level of 100. This algorithm also proposes the implementation of a fingerprint for the receiver to identify if the file received corresponds to the one sent by the issuer. The proposed algorithm carries out the encryption of the images in a single block, obtaining an adequate encryption speed. On the other hand, the cryptosystem has been subjected to several tests; after analyzing the results obtained, it is stated that this proposal is robust and capable of withstanding linear, differential, statistical, brute force, or modification attacks, as well as some of the better known as the discrete logarithm and the MOV. The proposed scheme provides quite competitive results compared with other previously proposed schemes.

Author Contributions

Conceptualization, V.S.-G., R.F.-C., H.P.-M.; methodology and software, E.H.-D.; formal analysis, E.H.-D., V.S.-G.; investigation, H.P.-M.; data analysis R.F.-C., V.S.-G.; writing and original draft preparation, E.H.-D., H.P.-M.; review, V.S.-G., R.F.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the Instituto Politécnico Nacional (National Polytechnic Institute), ESIME Culhuacan, CIDETEC, and CONACyT for their support to develop this researching work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Azmi, T.N.T.; Azma Abdullah, N.; Ab Rahman, N.H.; A Hamid, I.R.; Chuah, C.W. Image size variation influence on corrupted and non-viewable BMP image. Mater. Sci. Eng. 2017, 226, 1–12. [Google Scholar] [CrossRef] [Green Version]
  2. Norma Oficial Mexicana NOM-151-SCFI. Prácticas Comerciales, Requisitos que Deben Observarse para la Conservación de Mensajes de Datos; Diario Oficial de la Federación: Mexico City, Mexico, 2002; (In SpanicJ. Nashh). [Google Scholar]
  3. Hernández-Díaz, E.; Pérez-Meana, H.; Silva-García, V.M. Encryption of RGB images by means of a novel cryptosystem using elliptic curves and chaos. IEEE Lat. Am. Trans. 2020, 18, 1407–1415. [Google Scholar] [CrossRef]
  4. Abed Elgabar, E.E. Comparison of LSB steganography in BMP and JPEG images. Int. J. Soft Comput. Eng. 2013, 3, 91–95. [Google Scholar]
  5. Stanescu, D.; Stratulat, M.; Groza, V.; Ghergulescu, I.; Borca, D. Steganography in YUV color space. In Proceedings of the IEEE International Workshop on Robotic and Sensors Environments, Ottawa, ON, Canada, 12–13 October 2007; pp. 1–4. [Google Scholar]
  6. Schoof, R. The discrete logarithm problem. In Open Problems in Mathematics; Nash, J., Rassias, M., Eds.; Springer: Zurich, Switzerland, 2016; pp. 403–416. [Google Scholar]
  7. Biham, E.; Shamir, A. Differential cryptanalysis of DES-like cryptosystems. J. Cryptol. 1991, 14, 3–72. [Google Scholar] [CrossRef]
  8. Matsui, M. Linear cryptanalysis method for DES cipher. In Advances in Cryptology—EUROCRYPT ‘93; Goos, G., Hartmanis, J., Eds.; Springer: Lofthus, Norway, 1994; pp. 386–397. [Google Scholar]
  9. Abd-Elhafiez, W.M.; Reyad, O.; Mofaddel, M.A.; Fathy, M. Image encryption algorithm methodology based on multi-mapping image pixel. In Proceedings of the International Conference on Advanced Machine Learning Technologies and Applications (AMLTA2019), Cairo, Egypt, 28–30 March 2019; Hassanien, A.E., Azar, A.T., Gaber, T., Bhatnagar, R., Tolba, M.F., Eds.; Springer: Cairo, Egypt, 2019; Volume 921, pp. 645–655. [Google Scholar]
  10. Luo, Y.; Ouyang, X.; Liu, J.; Cao, L. An image encryption method based on elliptic curve elgamal encryption and chaotic systems. IEEE Access 2019, 7, 38507–38522. [Google Scholar] [CrossRef]
  11. Abdelfatah, E. Secure image transmission using chaotic-enhanced elliptic curve cryptography. IEEE Access 2019, 8, 3875–3890. [Google Scholar] [CrossRef]
  12. Laiphrakpam, D.S.; Khumanthem, M.S. A robust image encryption scheme based on chaotic system and elliptic curve over a finite field. Multimed. Tools Appl. 2018, 77, 8629–8652. [Google Scholar] [CrossRef]
  13. Wu, J.; Liao, X.; Yang, B. Color image encryption based on chaotic systems and elliptic curve ElGamal scheme. Signal Process. 2017, 141, 109–124. [Google Scholar] [CrossRef]
  14. Reyad, O.; Khalifa, H.S.; Kharabsheh, R. Image pixel permutation operation based on elliptic curve cryptography. Appl. Math. Inf. Sci. 2019, 13, 183–189. [Google Scholar]
  15. Miller, V.S. Use of elliptic curves in cryptography. In Advances in Cryptology—CRYPTO ‘85; Goos, G., Hartmanis, J., Eds.; Springer: Berlin, Germany, 1985; pp. 417–426. [Google Scholar]
  16. Koblitz, N. Elliptic curve cryptosystems. Math. Comput. 1987, 48, 203–209. [Google Scholar] [CrossRef]
  17. Stinson, D.R. Public-key cryptography and discrete logarithms. In Cryptography Theory and Practice, 3rd ed.; Rosen, K.H., Ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2006; pp. 54–56, 254–267. [Google Scholar]
  18. Liskov, M. Miller–rabin probabilistic primality test. In Encyclopedia of Cryptography and Security, 2nd ed.; Van Tilborg, H.C.A., Jajodia, S., Eds.; Springer: Boston, MA, USA, 2011; Volume 1, p. 784. [Google Scholar]
  19. Washington, L.C. The basic theory. In Elliptic Curves Number Theory and Cryptography, 2nd ed.; Rosen, K.H., Ed.; Chapman & Hall/CRC: College Park, MD, USA, 2008; pp. 47–58. [Google Scholar]
  20. Menezes, A.J.; Okamoto, T.; Vanstone, S.A. Reducing elliptic curve logarithms to logarithms in a finite field. IEEE Trans. Inf. Theory 1993, 39, 1639–1646. [Google Scholar] [CrossRef]
  21. Pohlig, S.C.; Hellman, M.E. An improved algorithm for computing logarithms over GF(p) and its cryptographic significance. IEEE Trans. Inf. Theory 1978, 24, 106–110. [Google Scholar] [CrossRef]
  22. Silva-García, V.M.; Flores-Carapia, R.; Rentería-Márquez, C.; Luna-Benoso, B.; Chimal-Eguía, J.C. Image cipher applications using the elliptical curve and chaos. Int. J. Appl. Math. Comput. Sci. 2020, 30, 377–391. [Google Scholar]
  23. Tiessen, T.; Knudsen, L.R.; Kölbl, S.; Lauridsen, M.M. Security of the AES with a secret S-box. In Proceedings of the 22nd International Workshop on Fast Software Encryption, Istanbul, Turkey, 8–11 March 2015; pp. 175–189. [Google Scholar]
  24. Edwards, C.H.; Penney, D.E.; Calvis, D.T. Nonlinear systems and phenomena. In Differential Equations and Boundary Value Problems: Computing and Modeling, 5th ed.; Hoffman, W., Ed.; Pearson: London, UK, 2015; pp. 426–436. [Google Scholar]
  25. Feistel, H.; Notz, W.A.; Smith, J.L. Some cryptographic techniques for machine-to-machine data communications. Proc. IEEE 1975, 63, 1545–1554. [Google Scholar] [CrossRef]
  26. Prouff, E. DPA Attacks and S-Boxes; Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2005; Volume 3557, pp. 424–441. [Google Scholar]
  27. Farwa, S.; Shah, T.; Idrees, L. A highly nonlinear S-box based on a fractional linear transformation. Springerplus 2016, 5, 1568. [Google Scholar] [CrossRef] [Green Version]
  28. Anees, A.; Phoebe-Chen, Y.-P. Designing secure substitution boxes based on permutation of symmetric group. Neural Comput. Appl. 2020, 32, 7045–7056. [Google Scholar] [CrossRef]
  29. Mohamed, K.; Hj Mohd-Ali, F.H.; Ariffin, S.; Hafiza-Zakaria, N.; Mohammed-Pauzi, M.N. An Improved AES S-box based on Fibonacci numbers and prime factor. Int. J. Netw. Secur. 2018, 20, 1206–1214. [Google Scholar]
  30. FIPS PUB 46-3. Data Encryption Standard (DES); FIPS 46-3, Data Encryption Standard (DES)|CSRC (nist.gov); National Institute of Standards and Technology: Gaithersburg, MD, USA, 1999.
  31. FIPS PUB 197. Advanced Encryption Standard (AES); FIPS 197, Advanced Encryption Standard (AES) (nist.gov); National Institute of Standards and Technology: Gaithersburg, MD, USA, 2001.
  32. BufferedImage Library. Available online: https://docs.oracle.com/javase/8/docs/api/java/awt/image/BufferedImage.html (accessed on 1 December 2020).
  33. Shannon, E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  34. Walpole, R.E.; Myers, R.H.; Myers, S.L.; Ye, K. Simple linear regression and correlation. In Probability and Statistics for Engineering and Sciences, 9th ed.; Pearson Education: Boston, MA, USA, 2012; pp. 430–435. [Google Scholar]
  35. Steel, R.G.D.; Torrie, J.H. Biostatistical, Principles and Procedures, 1st ed.; Mc. Graw Hill: New York, NY, USA, 1985; pp. 56–57. [Google Scholar]
  36. Silva-García, V.M.; González-Ramírez, M.D.; Flores-Carapia, R.; Vega-Alvarado, E.; Rodríguez-Escobar, E. A novel method for image encryption based on chaos and transcendental numbers. IEEE Access 2019, 7, 163729–163739. [Google Scholar] [CrossRef]
  37. Heys, H.M. A tutorial on linear and differential cryptanalysis. Cryptologia 2006, 22, 189–221. [Google Scholar] [CrossRef]
  38. Wu, Y.; Noonan, J.P.; Agaian, S. NPCR and UACI randomness tests for image encryption. J. Sel. Areas Telecommun. 2011, 1, 31–38. [Google Scholar]
  39. Diffie, W.; Hellman, M.E. New directions in cryptography. IEEE Trans. Inf. Theory 1976, 22, 644–654. [Google Scholar] [CrossRef] [Green Version]
  40. Alvarez, G.; Li, S. Some basic cryptographic requirements for chaos-based cryptosystems. Int. J. Bifurc. Chaos 2006, 16, 2129–2151. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow diagram that describes the algorithm for the search for the elliptic curves used in this research.
Figure 1. Flow diagram that describes the algorithm for the search for the elliptic curves used in this research.
Electronics 10 00413 g001
Figure 2. (a) Proposed S-Box and (b) Inverse of proposed S-Box.
Figure 2. (a) Proposed S-Box and (b) Inverse of proposed S-Box.
Electronics 10 00413 g002aElectronics 10 00413 g002b
Figure 3. Images for tests: (a) Lena.jpeg, (b) Babbon.jpeg, (c) Boat.jpeg, (d) Barbara.jpeg, (e) City.jpeg, (f) House.jpeg, (g) Security.jpeg.
Figure 3. Images for tests: (a) Lena.jpeg, (b) Babbon.jpeg, (c) Boat.jpeg, (d) Barbara.jpeg, (e) City.jpeg, (f) House.jpeg, (g) Security.jpeg.
Electronics 10 00413 g003
Figure 4. Images after encryption: (a) Lena_encrypted.jpeg, (b) Babbon_encrypted.jpeg, (c) Boat_ encrypted.jpeg, (d) Barbara_ encrypted.jpeg, (e) City_ encrypted.jpeg, (f) House_ encrypted.jpeg, (g) Security_ encrypted.jpeg.
Figure 4. Images after encryption: (a) Lena_encrypted.jpeg, (b) Babbon_encrypted.jpeg, (c) Boat_ encrypted.jpeg, (d) Barbara_ encrypted.jpeg, (e) City_ encrypted.jpeg, (f) House_ encrypted.jpeg, (g) Security_ encrypted.jpeg.
Electronics 10 00413 g004
Figure 5. (a) Histograms of Lena.jpeg, (b) Histograms of Lena_encrypted.jpeg.
Figure 5. (a) Histograms of Lena.jpeg, (b) Histograms of Lena_encrypted.jpeg.
Electronics 10 00413 g005
Figure 6. (a) Scatter plots of Lena.jpeg, (b) Scatter plots of Lena_encrypted.jpeg.
Figure 6. (a) Scatter plots of Lena.jpeg, (b) Scatter plots of Lena_encrypted.jpeg.
Electronics 10 00413 g006
Figure 7. Message.jpeg.
Figure 7. Message.jpeg.
Electronics 10 00413 g007
Figure 8. (a) Lena_encrypted.jpeg with 25% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg. after 25% occlusion attack, (c) Message_encrypted.jpeg with 25% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 25% occlusion attack.
Figure 8. (a) Lena_encrypted.jpeg with 25% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg. after 25% occlusion attack, (c) Message_encrypted.jpeg with 25% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 25% occlusion attack.
Electronics 10 00413 g008
Figure 9. (a) Lena_encrypted.jpeg with 50% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg after 50% occlusion attack, (c) Message_encrypted.jpeg with 50% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 50% occlusion attack.
Figure 9. (a) Lena_encrypted.jpeg with 50% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg after 50% occlusion attack, (c) Message_encrypted.jpeg with 50% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 50% occlusion attack.
Electronics 10 00413 g009
Figure 10. (a) Lena_encrypted.jpeg with 75% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg after 75% occlusion attack, (c) Message_encrypted.jpeg with 75% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 75% occlusion attack.
Figure 10. (a) Lena_encrypted.jpeg with 75% occlusion in different zones, (b) Decryption of Lena_encrypted.jpeg after 75% occlusion attack, (c) Message_encrypted.jpeg with 75% occlusion in different zones, (d) Decryption of Message_encrypted.jpeg after 75% occlusion attack.
Electronics 10 00413 g010
Figure 11. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 25%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 25%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 25%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 25%.
Figure 11. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 25%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 25%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 25%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 25%.
Electronics 10 00413 g011
Figure 12. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 50%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 50%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 50%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 50%.
Figure 12. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 50%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 50%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 50%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 50%.
Electronics 10 00413 g012
Figure 13. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 75%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 75%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 75%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 75%.
Figure 13. (a) Lena_encrypted.jpeg under a salt and pepper noise attack with a density of 75%, (b) Decryption of Lena_encrypted.jpeg after noise attack with a density of 75%, (c) Message_encrypted.jpeg under a salt and pepper noise attack with a density of 75%, (d) Decryption of Message_encrypted.jpeg after noise attack with a density of 75%.
Electronics 10 00413 g013
Table 1. Examples when chaos is done or not after 1000 iterations.
Table 1. Examples when chaos is done or not after 1000 iterations.
xsResultObservation
0.013.10970.5528279577750637Chaos occurs.
0.013.10980.5527757512022222Chaos does not take place.
0.0012.4750.595959595959595Chaos does not take place.
0.0012.4760.5961227786752827Chaos occurs.
Table 2. Comparison of cryptographic properties between the proposed S-Box and AES S-Box.
Table 2. Comparison of cryptographic properties between the proposed S-Box and AES S-Box.
Cryptographic PropertiesProposed S-BoxAES S-BoxExpected Value (EV) [28]
Algebraic Degree77 E V 10
Algebraic Immunity44 E V 4
Balance00 E V = 0
Confusion Coefficient Variance0.1428880.111 E V 0
Correlation Immunity00 E V 0
Nonlinearity100112 100 E V 120
Fixed and Opposite Fixed Points 00 E V 4
Propagation Characteristic00 E V 0
Robustness to Differential Cryptanalysis0.9610.984 E V < 0.98
Signal to Noise Ratio (SNR) Diff. Power Analysis8.8339.600 E V > 0.98
Transparency Order7.8017.860 E V < 7.8
Table 3. Entropy results.
Table 3. Entropy results.
ImageRed BandGreen BandBlue BandAverage Entropy
Lena7.99928255117.99939982157.99944521457.9993758623
Baboon7.99949438267.99935384777.99931386337.9993873646
Boat7.99927124027.99928187987.99925045617.9992678587
Barbara7.99950838977.99954477447.99944077857.9994979809
City7.99977596617.99979669707.99981503987.9997959010
House7.99992547277.99992259257.99993188897.9999266513
Security7.99998892777.99998880957.99998774607.9999884944
Table 4. Horizontal correlation coefficient results.
Table 4. Horizontal correlation coefficient results.
ImageRed BandGreen BandBlue BandAverage H.C.C.
Lena0.00279109110.00130756940.00126782270.0017888277
Baboon0.00251464950.00142692710.01716318870.0070349218
Boat0.00281135950.01088956960.00013505900.0046119960
Barbara0.00110784570.02452728850.00100703370.0088807226
City0.00143789430.01052968730.00593162780.0059664031
House0.00137076500.00269463360.00828138350.0041155940
Security0.00085862660.00501128730.00531378040.0037278981
Table 5. Vertical correlation coefficient results.
Table 5. Vertical correlation coefficient results.
ImageRed BandGreen BandBlue BandAverage V.C.C.
Lena0.00432276170.01958628900.00143252720.0084471926
Baboon0.03033894830.00291450880.00267126570.0119749076
Boat0.00927023970.00610607800.01689400170.0107567731
Barbara0.00274810280.00056385560.00464978500.0026539145
City0.02125114180.00375709190.00137356650.0087939334
House0.00646186210.01904887960.00055572700.0086888229
Security0.00357725130.00526164760.00342034550.0040864148
Table 6. Diagonal correlation coefficient results.
Table 6. Diagonal correlation coefficient results.
ImageRed BandGreen BandBlue BandAverage D.C.C.
Lena0.00401610970.00816569370.00286427750.0050153603
Baboon0.00618493340.01760824770.00212364950.0086389435
Boat0.00769699240.00197183010.00626195880.0053102604
Barbara0.00670012260.00175907210.00594051030.0047999017
City0.00621773860.00490870840.00530747830.0054779751
House0.00314398670.00438495220.00005074670.0025265619
Security0.00984614660.00160379190.00543118310.0056270405
Table 7. χ 2 test results.
Table 7. χ 2 test results.
ImageRed BandGreen BandBlue Band Average   χ 2
Lena248.3114536285.4836539252.6187869262.1379648
Baboon184.2003975234.2884898249.1635853222.5508242
Boat264.0449749261.5555760272.5446485266.0483998
Barbara273.4035741262.9516374221.0672448252.4741521
City286.2245508260.9730802236.9778295261.3918201
House254.2933231264.0568020232.1029117250.1510123
Security260.6685332263.5265517188.4832084237.5594311
Table 8. NPCR and UACI results.
Table 8. NPCR and UACI results.
ImageNPCR (%)UACI (%)
Lena99.603410290833.4703504432
Baboon99.617877883433.4530838915
Boat99.609248091333.4670956081
Barbara99.614144040133.4441091361
City99.623678051333.4647289134
House99.611346929533.4587882284
Security99.612574726933.4640385760
Table 9. Avalanche effect result.
Table 9. Avalanche effect result.
ImageAvalanche Effect (%)
Lena49.9823570251
Baboon50.0461260477
Boat49.9830722808
Barbara50.0082585841
City49.9927435980
House50.0056389702
Security50.0402729690
Table 10. Average entropy and correlation coefficient results of a black and white image.
Table 10. Average entropy and correlation coefficient results of a black and white image.
ImageDimensionsAverage
Entropy
Average
H.C.C.
Average
V.C.C.
Average
D.C.C.
White512 × 5127.99934582780.01020433710.02174384210.0117693713
Black512 × 5127.99928520490.01771579350.01605361080.0088100118
White720 × 5767.99959958150.00788892360.00151651200.0051731527
Black720 × 5767.99952395520.00223548850.00795108110.0029753841
White 1280 × 7207.99980726720.00218396640.01044614630.0246513849
Black 1280 × 7207.99979508480.00132470860.01513928630.0018480579
White 1920 × 10807.99990261770.01298273770.01442222510.0066555182
Black 1920 × 10807.99990083700.00721951630.01873552240.0035342527
White 4900 × 34647.99999052110.00152791810.00487574590.0081609654
Black 4900 × 34647.99998899520.00285200190.02316564320.0013768905
Table 11. Key sensitivity test results.
Table 11. Key sensitivity test results.
ImageDifferent PixelsPercentage Change
Lena260,83599.50065613%
Baboon260,90399.52659607%
Boat260,51799.37934875%
Barbara412,40599.44179205%
City917,22399.52506510%
House2,061,20999.40244020%
Security16,892,44199.52185158%
Table 12. Average encryption time measured.
Table 12. Average encryption time measured.
ImageDimensionsEncryption Time
Lena512 × 512 pixels0.123 s
Baboon512 × 512 pixels0.157 s
Boat512 × 512 pixels0.125 s
Barbara720 × 576 pixels0.189 s
City1280 × 720 pixels0.408 s
House1920 × 1080 pixels1.131 s
Security4900 × 3464 pixels7.946 s
Table 13. Performance comparison between proposed scheme and others recently proposed.
Table 13. Performance comparison between proposed scheme and others recently proposed.
Peer-Reviewed ArticleEntropyH.C.C.V.C.C.D.C.C.
Ours-Lena7.9993750.0017880.0084470.005015
Ours-Baboon7.9993870.0070340.0119740.008638
Ref. [10]-Lena7.9993000.0019000.0024000.001100
Ref. [10]-Baboon7.9993000.0024000.0011000.000800
Ref. [11]-Lena7.9998600.0015000.0006000.000900
Ref. [11]-Baboon7.9988400.0015000.0006400.000100
Ref. [12]-Lena7.9992000.0034750.0040130.007655
Ref. [12]-Baboon7.9991000.0003800.0052880.002308
Ref. [13]-Lena7.9912260.0001000.0089000.009100
Ref. [14]-Lena7.9962000.0031100.0542000.002310
Table 14. Performance comparison between proposed scheme and others recently proposed.
Table 14. Performance comparison between proposed scheme and others recently proposed.
Peer-Reviewed ArticleNPCRUACIAvalanche
Ours-Lena99.603433.470349.9824
Ours-Baboon99.617833.453050.0461
Ref. [10]-Lena99.611333.468250.0334
Ref. [10]-Baboon99.611233.491950.0411
Ref. [11]-Lena99.620433.489850.0419
Ref. [11]-Baboon99.610033.460050.0824
Ref. [12]-Lena--------49.9322
Ref. [12]-Baboon--------50.0324
Ref. [13]-Lena10033.4464----
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hernández-Díaz, E.; Pérez-Meana, H.; Silva-García, V.; Flores-Carapia, R. JPEG Images Encryption Scheme Using Elliptic Curves and A New S-Box Generated by Chaos. Electronics 2021, 10, 413. https://doi.org/10.3390/electronics10040413

AMA Style

Hernández-Díaz E, Pérez-Meana H, Silva-García V, Flores-Carapia R. JPEG Images Encryption Scheme Using Elliptic Curves and A New S-Box Generated by Chaos. Electronics. 2021; 10(4):413. https://doi.org/10.3390/electronics10040413

Chicago/Turabian Style

Hernández-Díaz, Erick, Héctor Pérez-Meana, Víctor Silva-García, and Rolando Flores-Carapia. 2021. "JPEG Images Encryption Scheme Using Elliptic Curves and A New S-Box Generated by Chaos" Electronics 10, no. 4: 413. https://doi.org/10.3390/electronics10040413

APA Style

Hernández-Díaz, E., Pérez-Meana, H., Silva-García, V., & Flores-Carapia, R. (2021). JPEG Images Encryption Scheme Using Elliptic Curves and A New S-Box Generated by Chaos. Electronics, 10(4), 413. https://doi.org/10.3390/electronics10040413

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

Article Metrics

Back to TopTop