arXiv:1708.01862v1 [cs.IT] 6 Aug 2017

A combination chaotic system and
application in color image encryption
R. Parvaza∗, M. Zarebniaa†

a

Department of Mathematics, University of Mohaghegh Ardabili, 56199-11367
Ardabil, Iran.
Abstract
In this paper, by using Logistic, Sine and Tent systems we define a
combination chaotic system. Some properties of the chaotic system are
studied by using figures and numerical results. A color image encryption
algorithm is introduced based on new chaotic system. Also this encryption
algorithm can be used for gray scale or binary images. The experimental
results of the encryption algorithm show that the encryption algorithm is
secure and practical.

Keywords: Encryption; Color image; Chaotic system; Cyclic shift.

1

Introduction

Image plays an important role in the data transfer. With rapid development network communication, image security has become increasingly important. The first step in chaotic encryption was introduced by Matters [1]. In
recent years, much attention has been given in the literature to the development, analysis and implementation of chaotic system for the image and the
data encryption; see, for example, [2, 3, 4, 5, 6]. Chaos maps as Logistic map,
Sine map and Tent map are used in image encryption algorithm because chaotic
maps have high sensitivity to their initial values and control parameters [7, 8].
Logistic, Sine and Tent maps have some disadvantages. These maps for some
values of r have chaotic behavior. Also these maps have non-uniform distribution over output. Many methods have been proposed to solve these problems,
for example see [9]. To overcome these problems, in this paper by using different
functions as sin(x), cos(x), . . ., we combine Logistic (or Sine) map and Tent map
then by using this combination the different chaotic systems can be found. In
the next step by using combination map, XOR operation and circ shift the color
∗ Corresponding

author: rparvaz@uma.ac.ir

† zarebnia@uma.ac.ir

1

image encryption algorithm is introduced. In this encryption algorithm in the
first step, color image is divided into twelve parts and then by using combination
map the encryption process for each of the parts is done, in the last step, we
combine parts and then the encryption process combination image is repeated.
The organization of this paper is as follows: In Section 2, Logistic-Tent combination map is explained. In Section 3, we present the color image encryption
algorithm. Simulation results and security analysis are given in Section 4. A
summary is given at the end of the paper in Section 5.

2

A combination chaotic system

In this section, we describe our combination of chaotic systems. Logistic,
Sine and Tent maps are defined as follows
xn+1 = L(r, xn ) := rxn (1 − xn ),
xn+1 = S(r, xn ) := r sin(πxn )/4,

xn+1 = T (r, xn ) :=


 rxn /2,


(2.1)
(2.2)

when xn < 0.5,
(2.3)

r(1 − xn )/2, when xn ≥ 0.5.

where parameter r ∈ (0, 4]. It is known that the Logistic system (Sine or
Tent system) for some values of r ∈ (0, 4] has chaotic behavior. Fig 2(a) and
Fig 4(a) show that the Logistic system for r = 2 has not chaotic behavior.
To overcome the above problem, the combination of chaotic system as Logistic
Tent system (LTS) introduced in [9]. Histogram of the Logistic Tent system is
showed in Fig 3(a). Chaotic range is not limited for the Logistic Tent system
but from Fig 3(a) we can see that the histogram of the Logistic Tent system is
not flat enough. Non-uniform distribution over output series leads to weakness
in the statistical attack. For solving this problem we add weights and functions
in the Logistic Tent system or the Sine Tent system as follows
xn+1 = Gr (xn ) :=

(β −r)x
 ω1 f1 ◦ F (r, xn ) + α1 g1 (rxn ) + ξ1 1 2 n mod 1,


when xn < 0.5,

n)
mod 1, when xn ≥ 0.5,
ω2 f2 ◦ F (r, xn ) + α2 g2 (rxn ) + ξ2 (β2 −r)(1−x
2
(2.4)

where F (r, xn ) is Logistic or Sine map. In (2.4), fi (x) and gi (x) (i = 1, 2) can be
considered as ax, sin(ax), cos(ax), tan(ax), cot(ax), exp(ax), log(ax) (where a is
real constant) and any other appropriate function. Also in the above formula,
ωi , αi , ξi and βi (i = 1, 2) are real numbers and parameter r ∈ (0, 4]. In next
step to investigate the properties of the new system, we consider the following
cases
2

Figure 1:

The combination chaotic system.

• i: 2ω1 = ω2 = 20, −α1 = α2 = −2, 2ξ1 = ξ2 = 4, β1 = 4, β2 =
−20, f1 (x) = sin(x), f2 (x) = exp(x), g1 (t) = cot(x), g2 = cos(πx), F (r, xn ) =
L(r, xn ).
• ii: ω1 = ω2 = 20, −α1 = α2 = 0, 2ξ1 = ξ2 = 1, β1 = β2 = 80, f1 (x) =
sin(x), f2 (x) = exp(x), F (r, xn ) = L(r, xn ).
• iii: ω1 = ω2 = 1, α1 = α2 = 1, ξ1 = 7, ξ2 = 15, β1 = 2β2 = 40, f1 (x) =
cos(x), f2 (x) = tan(x), g1 (x) = tan(x), g2 (x) = x, F (r, xn ) = S(r, xn ).
For discrete time system as (2.4), the Lyapunov exponent for an orbit starting
with x0 is defined as follows
n−1
1X
ln |G′r (xi )|.
n→∞ n
i=0

LE(x0 , r) := lim

(2.5)

The degree of ”sensitivity to initial conditions” can be measured by using the
Lyapunov exponent. In [10] for the Lyapunov exponent the following theorem
is mentioned.
Theorem 2.1. If at least one of the average Lyapunov exponents is positive,
then the system is chaotic; if the average Lyapunov exponent is negative, then the
orbit is periodic and when the average Lyapunov exponent is zero, a bifurcation
occurs.

3

The Lyapunov exponent plot for Case (i) is given in Fig 2(b). From this
figure we can see that for some values of r the Lyapunov exponent are negative.
The Bifurcation diagram for the Case (i) has been shown in Fig. 5(b). From
this figure we can see that white lines appear in places where the Lyapunov
exponent are negative. Also for Case (ii) and Case (iii) from Fig.s 2(c)-2(d)
we can say that for all values of r ∈ (0, 4] the Lyapunov exponent are positive.
Also the Cobweb plots 4(b), 4(c) and 4(d) show chaotic behavior for the Case
(i), Case (ii) and the Case (iii). The Fig.s 3(c)-3(d) show that the Case (ii) and
Case (iii) have uniform distribution over output range. Also from Fig. 4(b) we
can see that the Case (i) has not uniform distribution over its output range.
Two orbits of the Case (i) and Case(ii) are shown in Figs. 5(a)-5(b). By using
this figure we can see that the Case (i) and Case(ii) are much more sensitive
to the starting points. By using the Lyapunov Exponents figures we can say
that the Lyapunov Exponents of two examples are all larger than the Logistic
map. Also the distribution of new combination chaotic system is more uniform
than the distribution Logistic Tent system. Then by using suitable functions
and parameters, appropriate chaotic system can be found.
In the next section, as application of the proposed chaotic system, we introduce
an image encryption algorithm.
Lyapunov exponent

Lyapunov exponent
6

0.5

0
5
−0.5

−1

4

LE

LE

−1.5

−2

3

−2.5
2
−3

−3.5

1

−4
0

−4.5
−1

0

1

2

3

4

5

−2

−1

0

1

2

3

4

5

6

r

r

(a)

(b)

Lyapunov exponent

Lyapunov exponent

5.5
6.5

5
6

4.5

LE

LE

5.5

4

5

3.5
4.5

3
4

2.5
0.5

1

1.5

2

2.5

3

3.5

4

0.5

r

1.5

2

2.5

3

3.5

4

r

(c)

Figure 2:

1

(d)

Lyapunov exponent plot for the (a) Logistic map, (b) Case (i), (c) Case (ii), (d)

Case (iii).

3

Proposed encryption and decryption process

In this section, we introduce some details about the proposed encryption and
decryption algorithm.

4

2000
1800

1800

1600

1600

1400

Frequency

Frequency

1400

1200

1000

800

1200

1000

800

600

600

400

400

200

200

0

0
0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0

0.1

0.2

0.3

0.4

0.5

0.7

0.8

0.9

1

(b)

1800

1800

1600

1600

1400

1400

Frequency

Frequency

(a)

1200

1000

800

1200

1000

600

800

600

400

400

200

200

0

0
0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

0

0.1

0.2

0.3

0.4

0.5

xn

0.6

0.7

0.8

0.9

1

xn

(c)

Figure 3:

0.6

xn

xn

(d)

Histogram of the (a) Logistic Tent map, (b) Case (i), (c) Case (ii), (d) Case(iii).
Cobweb plot

0.8

0.8

0.7

0.7

0.6

0.6

xn+1

1

0.9

0.5

x

n+1

Cobweb plot
1

0.9

0.5

0.4

0.4

0.3

0.3

0.2

0.2

0.1

0

0.1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

0

1

0

0.1

0.2

0.3

0.4

0.5

xn

(a)

0.7

0.8

0.9

1

(b)

Cobweb plot

Cobweb plot
1

0.9

0.9

0.8

0.8

0.7

0.7

0.6

0.6

n+1

1

0.5

0.5

x

xn+1

0.6

xn

0.4

0.4

0.3

0.3

0.2

0.2

0.1

0

0.1

0

0.1

0.2

0.3

0.4

0.5

x

0.6

0.7

0.8

0.9

0

1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

x

n

n

(c)

(d)

Figure 4:

Cobweb plot for the (a) Logistic map, (b) Case (i), (c) Case (ii), (d) Case(iii)
(blue: r=2, red: r=4).

3.1

Encryption process

We assume that the size of the input color image is m × n. The encryption
process is written in the following steps.

Step 1.

For the color plain image A, we determine yk , yk0 (k = 1, 2, 3) as
5

(a)

(b)

(c)

Figure 5:

(d)

Bifurcation diagram for the (a) Logistic Tent map, (b) Case (i), (c) Case (ii), (d)

1

1

0.9

0.9

0.8

0.8

0.7

0.7

0.6

0.6

xn

x

n

Case(iii).

0.5

0.5

0.4

0.4

0.3

0.3

0.2

0.2

0.1

0

0.1

0

50

100

150

200

250

300

350

400

450

0

500

0

50

100

150

200

Iteration

250

300

350

400

450

500

Iteration

(a)

(b)

Figure 6:

(a) Case (i) with x0 = 0.5 (blue) and x0 = 0.5 + 10−15 (red); (b) Case (ii) with
x0 = 0.5 (blue) and x0 = 0.5 + 10−15 (red).

follows
yk :=

X

A(i, j, k), k = 1, 2, 3,

i,j

yk0 :=

yk
, k = 1, 2, 3.
m × n × 255

Therefore in this step we can find the different secret keys for different plain
images. Also to generate different keys in each iteration, we consider y00 and y40
as random
gray scale image or binary image, we consider
P numbers in [0, 1]. For
y1
,y00 := rand() and y20 := rand(). Also for
y1 := i,j A(i, j), y10 := m×n×255
gray scale image or binary image yk , yk0 (k = 3, 4) are not used in the encryption
algorithm.

6

Figure 7:

The encryption process of the proposed algorithm.

Step 2. By using cyclic shift operation we find
y1
y2
y3
⌋ − ⌊ ⌋ ⌊ ⌋]),
4
4
4
y2 + y3
y3 + y1
y1 + y2
⌋ ⌊
⌋ −⌊
⌋]),
A = circshif t(A, [⌊
10
10
10
A = circshif t(A, [⌊

where circshif t(A, r) circularly shifts the values in array A by r positions. For
gray scale image or binary image, this step has been changed as follows
y1
y1
⌋ − ⌊ ⌋]),
4
3
2y1
2y1
A = circshif t(A, [⌊
⌋ ⌊
⌋]).
10
5

A = circshif t(A, [⌊

Step 3. In this step, by using Div(A) image A is divided into twelve parts.
In the Div(A) function, image is divided into three color (red, green, blue) then
each part is divided into four equal parts. However in the encryption algorithm,
other division functions can be used. We can find
{Ari , Agi , Abi }4i=1 = Div(A).
Step 4. We define Λ(y 0 , r, n) as follows
Λ(y 0 , r, n) := (x0 , x1 , . . . , xn−1 ),

(3.1)

where xi (i = 1, . . . , n−1) are defined by using (2.4) with x0 = y 0 and parameter
r. Also we consider Ω(Ki , r) for all i ∈ {1, 2, 3, 4} as a matrix with the following

7

elements
Ω(Ki , r)(1, :) := Λ(yi0 , ri , ⌊n/2⌋), i = 1, 3,
Ω(Ki , r)(1, :) :=

Λ(yi0 , ri , n

(3.2)

− ⌊n/2⌋), i = 2, 4,

(3.3)
T

Ω(K1 , r)(2 : ⌊m/2⌋, j) := Λ(Ω(K1 , r1 )(1, j), r1 , ⌊n/2⌋) , j = 1, . . . ⌊n/2⌋,
(3.4)
T

Ω(K2 , r)(2 : ⌊m/2⌋, j) := Λ(Ω(K2 , r2 )(1, j), r2 , n − ⌊n/2⌋) , j = 1, . . . , n − ⌊n/2⌋,
(3.5)
Ω(K3 , r)(2 : m − ⌊m/2⌋, j) := Λ(Ω(K3 , r3 )(1, j), r3 , ⌊n/2⌋)T , j = 1, . . . ⌊n/2⌋,
(3.6)
Ω(K4 , r)(2 : m − ⌊m/2⌋, j) := Λ(Ω(K4 , r4 )(1, j), r4 , n − ⌊n/2⌋)T , j = 1, . . . , n − ⌊n/2⌋.
(3.7)
(a) For i = 1, . . . , 4 we find
n
Ki := Λ(yi0 , ri , ⌊ ⌋),
2
n
ti := mod(Ki (1, ⌊ ⌋) × 105 , m × n).
2
Also for i = 0 we define
P4

i=1 ti

⌋,
4
K0 := Λ(y00 , r0 , m × n).

t0 := ⌊

(b) We define Zi1 , Zi2 (i = 0, . . . , 4) as
Zi1 := Ω(Ki , ri ) × 255,
Zi2 := circshif t(Zi1 , ti ).
Also Z0j (j = 1, 2) are defined as follows
f or i = 1 : m do
f or j = 1 : n do
Z01 (i, j) := K0 (1, j + (i − 1)n) × 255,
end
end
Z02 := circshif t(Z01 , −t0 ),
(c) By using bitxor operation matrixs Aji (i = 1, 2, 3, 4, j = r, g, b) are encrypted as follows

8

f or i = 1 : 4do
f or j = r, g, b do
f or w = 1 : 2do
DAji := bitxor(Aji , Ziw ),
end
end
end
Step 6. In the encryption algorithm, the inverse function of the Div(·) is
shown by EDiv(·).
(a) By using the inverse function of the Div(·) we find
EA := EDiv({DAri , DAgi , DAbi }4i=1 ),
where EA denotes the joined encrypted image.
(b) We transform matrix EA as follows
EA(:, :, 1) = circshif t(EA(:, :, 1), [−y1 + t0 t0 ]),
EA(:, :, 2) = circshif t(EA(:, :, 2), [t0 − y2 + t0 ]),
EA(:, :, 3) = circshif t(EA(:, :, 3), [t0 − y3 + t0 ]).
(c) The final encrypted image is found as follows
f or i = 1 : 3
EA(:, :, i) = bitxor(EA(:, :, i), Z01 ),
EA(:, :, i) = bitxor(EA(:, :, i), Z02 ),
end

EA = circshif t(EA, [y1 + y2 y2 + y3 y1 + y3 ]),
y1
y2
y3
EA = circshif t(EA, [⌊ ⌋ ⌊ ⌋ − ⌊ ⌋]).
50
50
50

(3.8)
(3.9)

For gray scale image or binary image, Eq.s (3.8)-(3.9) has been changed as
follows
EA = circshif t(EA, [y1 − 2y1 ),
y1
y1
EA = circshif t(EA, [⌊ ⌋ − ⌊ ⌋]).
50
35

3.2

Decryption process

The decryption process is the inverse process of the encryption process, so we
remove the decryption details.
9

4

Simulation results and security analysis

4.1

Simulation results

In order to illustrate the encryption algorithm results and security analysis, we
consider some test problems. In this section in all figures we use Case (ii) in
the encryption process. We have computed the numerical results by MATLAB
7.11.0 programming. As plain images, we use color images “Sailboat on lake” ,
“Airplane (F-16)” (512×512×3 pixels), “Lena” (256×256×3, 512×512×3 pixels)
and gray scale image “Lena” (256 × 256 , 512 × 512 pixels), “girl” (256 × 256
pixels). Also “Horse” (444×455 pixels) is used as binary image. The simulations
results are given in Figs. 8-14. In Figs. 8-9 we show the encryption algorithm
results. The binary image has only two pixel values 0 and 1, therefore it is a
difficult case for encryption. As can be seen in the Fig.9, by using the proposed
algorithm for the binary image, the image is changed to a noise-like encrypted
image with flat histogram. The running time of encryption and decryption are
shown in Table 1.

Table 1:

Encryption and decryption times of the color images by using Case (ii).

Image
Lena (512 × 512)
Sailboat on lake
Airplane (F-16)

4.2

Encryption time
1.130 s
1.105 s
1.139 s

Decryption time
1.044 s
1.044 s
1.024 s

Security analysis

In this subsection, the different security parameters are discussed, for details
see [11, 12].
4.2.1

Security key analysis

From the proposed algorithm, we can say that the security keys of the algorithm
are composed of ten parameters {ri }4n=0 and {yi0 }4n=0 , The range for ri (i =
0, . . . 4) are (0, 4] Also yi0 (i = 0, . . . 4) are in range of [0, 1]. If in the image
encryption algorithm we consider the precision as 10−15 the key space is almost
10150 , and this space is sufficiently large to resist the brute force attack.
To show key sensitivity, in the encryption process we consider key as follows
(this key is used to obtain all results)
k := (r0 , r1 , r2 , r3 , r4 ) = (2, 1, 2, 3.5, 1.75).
10

(4.1)

Figure 8:

(a) Original image, (b-c-d) Histograms of the R, G, B components of the original
image, (e) Encrypted image, (f-g-h) Histograms of the R, G, B components of the encrypted
image (Sailboat on lake).

Figure 9: (a-b) Original and encrypted image, (e-f) Histograms of the original image and
encrypted image (Lena -256 × 256), respectively, (c-d) Original and encrypted image, (g-h)
Histograms of the original image and encrypted image (Horse), respectively.

11

In the decryption process, we use a small change for r0 as r0 = 2 + 10−15 ,
we consider new key as K1 . Simulation results in Fig. 10 show that by using
new key we can not reconstruct the original image. Therefore, the proposed
algorithm has high key sensitivity.

Figure 10:

(a) Encrypted image (Sailboat on lake), (b) Decrypted image by using K, (c)

Decrypted image by using K1 .

4.2.2

Statistical analysis

In this section, we study statistical analysis as part of the security analysis. It
is clear that a good encrypted image should be unrecognized hence the correlation values of a good encrypted image are close to zero. Table 3 shows the
correlation values for the original images and encrypted images. We have found
the correlation values using formula [13]
Cxy =

E[(x − µx )(y − µy )]
,
σx σy

(4.2)

where E[·] denotes the expectation value, µ is the mean value and σ represents
standard deviation. From Table 3, we can see that the original images have
high correlation values while encrypted images have very low correlation values. Also correlation distributions for the original image and encrypted image
of Lena image are shown in Fig.11.

12

Figure 11:

Correlation of neighborhood pixels at different directions before and after

encryption of Lena ( 512 × 512 pixels).

13

Table 2: Correlation coefficients of the encrypted color image of ”Sailboat on lake”.
Chaos map
Component Horizontal Vertical
Diagonal

Case (i)

Case (ii)

Case (iii)

Logistic Tent

Original image

Table 3:

R
G
B
R
G
B
R
G
B
R
G
B
R
G
B

-0.0027
0.0009
-0.0017
0.0007
0.0029
-0.0007
0.0003
0.0029
-0.0007
-0.0029
-0.0009
-0.0009
0.9558
0.9715
0.9710

-0.0003
-0.0019
-0.0005
-0.0003
-0.0005
0.0009
0.0018
-0.0005
-0.0021
0.0011
0.0033
0.0025
0.9541
0.9663
0.9694

Diagonal

(lower left to top right)

(lower right to top left)

-0.0016
-0.0008
-0.0019
-0.0006
-0.0003
0.0010
0.0020
-0.0028
0.0007
0.0028
0.0035
0.0012
0.9373
0.9510
0.9512

0.0019
0.0027
0.0024
-0.0002
-0.0003
-0.0023
0.0027
0.0030
-0.0021
-0.0005
-0.0002
0.0005
0.9420
0.9530
0.9530

Comparison of Correlation coefficients of image of “Lena” (512 × 512).

Method
The proposed
Method in [14]

Horizontal
-0.0036
0.0139

Vertical
-0.0020
0.0073

Diagonal
-0.0026
0.0104

The other test is the information entropy. The values of entropy are in range
of [0, 8]. This test is used for evaluating the randomness of an image. If the
value of entropy of encrypted image is close to the maximum value means the
excellent random property. The information entropy is defined as follows [12]
H(k) = −

w−1
X

P (ki ) log2 P (ki ),

(4.3)

i=0

where w is the gray level and P (·) denotes the probability of symbol. Results
for the information entropy are tabulated in Table 4. From Table 4, we can say
that test results for the encryption algorithm are close to the maximum value.
Also Figs. 8-9 show that the histograms of plain images are not flat while the
histograms of encrypted images are in flat distributions.
From the above discussion, it can be concluded that the proposed algorithm
has stronger ability to resist statistical attacks.

14

Table 4: Informaion entropies of the encrypted
Chaos map
Image
R
Lena (512 × 512 × 3)
7.9970
Case (i)
Sailboat on lake
7.9993
Airplane (F-16)
7.9994
Lena (512 × 512 × 3)
7.9967
Case (ii)
Sailboat on lake
7.9992
Airplane (F-16)
7.9994
Lena (512 × 512 × 3)
7.9970
Case (iii)
Sailboat on lake
7.9994
Airplane (F-16)
7.9991
Lena (512 × 512 × 3)
7.9975
Logistic Tent
Sailboat on lake
7.9993
Airplane (F-16)
7.9993
4.2.3

color image.

G
7.9972
7.9993
7.9992
7.9967
7.9993
7.9993
7.9973
7.9993
7.9992
7.9970
7.9992
7.9994

B
7.9970
7.9992
7.9993
7.9972
7.9993
7.9994
7.9975
7.9993
7.9994
7.9970
7.9993
7.9994

Sensitivity analysis

NPCR (Number of Pixels change Rate) denotes the number of pixels change rate
while one pixel of plain image changed. Also UACI (Unified Average Changing
Intensity) measures the average intensity of differences between the plain image
and encrypted image. The ideal value for NPCR is 100% while the ideal value
for UACI is 33.33%. When the value of NPCR gets closer to 100%, the encryption algorithm is more sensitive to the changing of plain image, therefore the
algorithm can effectively resist plaintext attack. Also when the value of NPCR
gets closer to 33.33%, the proposed algorithm can effectively resist differential
attack. In computing, we consider NPCR and UACI as follows
P
i,j D(i, j)
N P CR =
× 100%,
(4.4)
m ×P
n
1  i,J |C1 (i, j) − C2 (i, j)| 
× 100%,
(4.5)
U ACI =
m×n
255
where
D(i, j) :=


 1, when C1 (i, j) 6= C2 (i, j),


(4.6)

0, when C1 (i, j) = C2 (i, j).

In above formulae, C1 (i, j) and C2 (i, j) are denoted the cipher image before and
after one pixel of the plain image is changed. In the Lena, Sailboat on lake and
Airplane images, A(200, 150, 1), A(100, 250, 2) and A(180, 334, 3) are changed
to 0, respectively. As can be seen in Table 5 , results are close to ideal values.

15

Table 5:
Image
Lena (512 × 512 × 3)
Sailboat on lake
Airplane (F-16)

Figure 12:

UACI and NPCR of the encrypted color image.

R
33.4457
33.4617
33.4744

UACI
G
33.5589
33.3928
33.4482

B
33.5243
33.5061
33.4813

R
99.6078
99.6082
99.6353

NPCR
G
99.6140
99.6231
99.6059

B
99.6033
99.6048
99.6021

Cropped attack experiment: (a) Encrypted image, (b) cropped attack image, (c)

Decryption result, (d-e-f) Decryption result for the R, G, B components, respectively.

4.2.4

Noise and Data loss attacks

In continuation of discussion, by using the simulation results we study noise
and data loss attacks. In real applications, a part of the encrypted image may
be lost during transmission. A proper encryption algorithm should resist the
data loss and noise attacks. In Fig. 12, we remove 200 × 200 of encrypted
image. In fact we consider A(80 : 280, 80 : 280) = 0200×200 . Also in Fig. 13, for
encrypted image we add Gaussian noise with zero-mean and different variance.
After the decryption process, we can see that the reconstructed images contain
most of original visual in formation and we can recognize the original image
from decrypted image.

16

Figure 13:

Noise attack experiment with Gaussian noise (a-b-c) noise attack image with

var=0.1,0.4,0.8, respectively, (d-e-f) decryption with var=0.1,0.4,0.8, respectively.

4.2.5

Chosen-plain text attack

In the encryption algorithm, in Step 1, random numbers are used. So this algorithm create different encrypted image each time when the encryption algorithm
is applied to the same image. In Fig. 14 we use K as key and we run the encryption algorithm twice. The first and second encrypted image are consider as
C1 (Fig. 14(b)) and C2 (Fig. 14(c)) , respectively. To illustrate the difference
between the two images, we use pixel-to-pixel difference as |C1 − C2 |. As can
be seen in Fig. 10, two encrypted images are different. Then the our algorithm
is able to withstand the chosen-plain text attack.

5

Conclusion

In this paper, we have constructed a new combination chaotic system based
on Logistic and Tent systems. By using this new combination chaotic system
a large number of chaotic map can be produced. Also we have proposed a
new image encryption algorithm based on combination chaotic system. It is
shown that the proposed encryption algorithm can effectively resist differential,
statistical, noise, data loss, chosen-plain text attacks.

17

Figure 14:

(a-e) The original image and its histogram, (b-f) the first encrypted image and

its histogram, (c-g) the second encrypted image and its histogram, (d-h) the pixel-to-pixel
difference and its histogram.

References
[1] R. Matthews, On the derivation of a chaotic encryption algorithm, Cryptologia 4 (1989) 29-42.
[2] Ü. Çavuşoğlu, S. Kaçarb, I. Pehlivanb, A. Zengina, Secure image encryption algorithm design using a novel chaos based S-Box, Chaos, Solitons &
Fractals 95 (2017) 92-101.
[3] A. Kanso, M. Ghebleh, An algorithm for encryption of secret images into
meaningful images, Optics and Lasers in Engineering 90 (2017) 196-208.
[4] X.Y. Wang, L. Yang, R.Liu, A. Kadir, A chaotic image encryption algorithm based on perceptron model, Nonlinear Dynamics 62 (2010) 615-621
[5] G. Gu, J. Ling, A fast image encryption method by using chaotic 3D cat
maps, Optik - International Journal for Light and Electron Optics 125
(2014) 4700-4705.
[6] M. Kumar, A. Vaish, Encryption of color images using MSVD in DCST
domain, Optics and Lasers in Engineering 3 (1990) 278-285.
[7] M. Ausloos, M. Dirickx, The Logistic Map and the Route to Chaos: From
the Beginnings to Modern Applications, Springer 2006.

18

[8] R. Hilborn, Chaos and Nonlinear Dynamics: An Introduction for Scientists
and Engineers, Oxford University Press 2000.
[9] S.C. Satapathy, A. Govardhan, K.S. Raju, J.K. Mandal, Emerging ICT
for Bridging the Future - Proceedings of the 49th Annual Convention of
the Computer Society of India (CSI), Springer International Publishing
Switzerland 2015.
[10] S. Lynch, Dynamical Systems with Applications using MATLAB r, Second
Edition, Birkhäuser Boston 2014.
[11] A. Uhl, A. Pommer, Image and Video Encryption From Digital Rights
Management to Secured Personal Communication, Springer 2004.
[12] F.A El-Samie, H.H. Ahmed, I.F. Elashry, M.H. Shahieen, O.S. Faragallah,
E.M. El-Rabaie, S.A. Alshebeili, Image Encryption: A Communication
Perspective, CRC Press 2014.
[13] N. Zheng, J. Xue, Statistical learning and pattern analysis for image and
video processing, Springer Science & Business Media, 2009.
[14] G. Gu, J. Ling, A fast image encryption method by using chaotic 3D cat
maps, Optik-International Journal for Light and Electron Optics, 17 (2014)
4700-4705.

19

