Next Article in Journal
A Review of Remote Sensing Image Dehazing
Next Article in Special Issue
Artificial Neural Networks to Solve the Singular Model with Neumann–Robin, Dirichlet and Neumann Boundary Conditions
Previous Article in Journal
Weak Calibration of a Visible Light Positioning System Based on a Position-Sensitive Detector: Positioning Error Assessment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Medical Images Encryption Based on Adaptive-Robust Multi-Mode Synchronization of Chen Hyper-Chaotic Systems

1
Faculty of Electrical Engineering, Zabol Branch, Islamic Azad University, Zabol 1939598616, Iran
2
Electrical and Computer Engineering Faculty, Semnan University, Semnan 3513119111, Iran
3
Faculty of Electrical Engineering, Biomedical Data Acquisition Lab (BDAL), K. N. Toosi University of Technology, Tehran 1631714191, Iran
4
Faculty of Electrical Engineering, Gonabad Branch, Islamic Azad University, Gonabad 6518115743, Iran
5
Faculty of Engineering, Mashhad Branch, Islamic Azad University, Mashhad 91735413, Iran
6
Computer Engineering Department, Ferdowsi University of Mashhad, Mashhad 9177948974, Iran
7
Institute for Intelligent Systems Research and Innovation (IISRI), Deakin University, Waurn Ponds, VIC 3217, Australia
8
Department of Signal Theory, Networking and Communications, Universidad de Granada, 52005 Granada, Spain
*
Author to whom correspondence should be addressed.
Submission received: 10 May 2021 / Revised: 1 June 2021 / Accepted: 2 June 2021 / Published: 7 June 2021
(This article belongs to the Special Issue Explainable AI in Medical Sensors)

Abstract

:
In this paper, a novel medical image encryption method based on multi-mode synchronization of hyper-chaotic systems is presented. The synchronization of hyper-chaotic systems is of great significance in secure communication tasks such as encryption of images. Multi-mode synchronization is a novel and highly complex issue, especially if there is uncertainty and disturbance. In this work, an adaptive-robust controller is designed for multimode synchronized chaotic systems with variable and unknown parameters, despite the bounded disturbance and uncertainty with a known function in two modes. In the first case, it is a main system with some response systems, and in the second case, it is a circular synchronization. Using theorems it is proved that the two synchronization methods are equivalent. Our results show that, we are able to obtain the convergence of synchronization error and parameter estimation error to zero using Lyapunov’s method. The new laws to update time-varying parameters, estimating disturbance and uncertainty bounds are proposed such that stability of system is guaranteed. To assess the performance of the proposed synchronization method, various statistical analyzes were carried out on the encrypted medical images and standard benchmark images. The results show effective performance of the proposed synchronization technique in the medical images encryption for telemedicine application.

1. Introduction

Recently, telemedicine systems have been introduced to assist the diagnosis and treatment of various diseases [1,2,3]. These systems, through telecommunication systems and clinical data registration devices, have been able to yield good medical services between specialists and patients easily; the significance of this issue arises when some patients for any reason cannot attend medical specialty centers [4,5]. Specialist physicians can apply telemedicine to diagnose the disease using clinical data such as medical images and signals.
When sending and receiving clinical data in telecommunication channels, information security is one of the main dilemmas of telemedicine systems [6]; because telecommunication channels are not secure enough to send medical information. Various medical information of patients is confidential, and at the time of sending and receiving from the patient to the specialist must have a sufficient security communication channel and cannot be accessed by other individuals and organizations [7,8]. In order to enhance the security of patients’ information at the time of sending, various approaches are exploited that the most important of which are cryptographic techniques.
One of the most important fields of cryptography is chaotic systems [9,10]. The most significant feature of chaotic systems is the high reaction to small changes in the initial conditions [11,12]. Much research has been conducted in this field for the last two decades since Correll et al. [13,14,15] introduced the synchronization of chaotic systems. In recent years, the employment of chaotic systems in clinical data encryption has grown remarkably [16]. Because chaos systems mainly guarantee medical information confidentiality when sending and receiving from telecommunication channels [16].
Data encryption is usually done by chaotic or hyper-chaotic methods [17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37]. In some researches, chaotic or hyper-chaotic systems are used to encrypt peripheral data [38,39,40,41,42,43]. In recent years, researchers have focused on the encryption of medical data using chaotic methods [17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]. Medical data contains important information about patients [41,44,45,46]. Therefore, the confidentiality of medical information is essential when sending it through telecommunication channels. Chaotic theory, as a nonlinear and complex phenomenon, plays an important role in increasing the security of medical data. In the following, a number of conducted investigations using cryptographic procedures in clinical data have been examined.
Lin et al. [17] used chaos theory to the encryption of electroencephalogram (EEG) signals. In this paper, the chaos-based EEG encryption system involves three levels of encryption and is implemented in the C# programming language.
A chaotic theory-based cryptographic method for clinical signals is presented in [18]. The proposed encryption algorithm is based on a logistic map with double chaotic layer encryption (DCLE). In this study, EEG, electrocardiograms (ECG), and blood pressure (BP) data from the PhysioBank database were acquired.
The idea of the optical chaos method for secure EEG signals transmission in telemedicine applications has been discussed in Shahzadi et al.’s research [19]. In the proposed scheme, a semiconductor laser source is utilized to generate optical chaos and conceals the EEG signal in the chaotic signal before transmission over the optical fiber. The generated data is then sent over optical fiber.
Applying chaotic modulation in intrinsic mode functions to encrypt medical ECG and EEG data is done in [20]. The results prove that the proposed method for encrypting and decrypting EEG and ECG signals containing chaotic behavior has satisfactory results.
In the research of Ibrahim et al. [21], medical image encryption was conducted using dynamic S-boxes and chaotic maps. Experimental results explicate that the proposed scheme has successfully passed all security tests in medical image encryption.
Gafsi et al. [22] proposed an improved chaotic system for medical image encryption. In this paper, a complex chaotic-based pseudorandom number generator (PRNG) has been designed to generate a high-quality encryption key. The generated key indicates the high complexity behavior of the approach.
In [23], a medical image encryption method has been proposed based on a chaotic map and fractional discrete cosine transform (FrDCT) coefficients. The proposed method comprises two steps of applying FrDCT on an image and, after that, a chaotic map on the coefficients.
The fourth-order chaotic system for clinical data encryption introduced in reference [24]. They calculated a variety of significant statistical parameters, such as correlation for the encrypted images in the proposed scheme.
An encryption algorithm based on the properties of block encryption, 4-dimensional logistics map and, DNA systems has been introduced by Stalin et al. [25]. In this method, multiple key sequences are generated first. In the second step, the pixel is substituted by a 4D nonlinear logistic map. Finally, encryption is performed using DNA rules to ensure that the various blocks are securely encrypted.
A chaotic system based on the Arnold cat method for medical images encryption in reference [26] is presented. The results demonstrate the high security and robustness of the proposed approach for secure medical image transmission.
In another research, Choi et al. [27] presented a secure and robust medical images encryption method. The proposed method include two permutations based on plain image and Latin square (PPILS) and bi-directional adaptive diffusion.
Choi et al. [28] proposed the NCA-based color medical image encryption algorithm and the generalized 3D chaotic cat map in another study. In this work, the NCA is an effective PRNG generating nonlinear pseudorandom sequences. Generalized 3D chaotic cat mapping is also adopted to create effective shuffle in encryption of color medical images.
A medical image encryption system based on multiple chaotic and message-digest algorithm 5 (MD5) has been discussed in [29]. The system employs two rounds of chaotic diffusion using XOR bit operation to encrypt an image. The results show that the proposed method is highly efficient in medical image encryption.
Designing a cryptographic algorithm based on pseudorandom sequence generation using a generalized double humped logistic map was done in [30], which ensures the transmission of magnetic resonance imaging (MRI) and X-ray images with high security. In this work, various statistical analyzes have been exploited.
Chaotic map-based authentication for telecare medicine information systems (TMISs) was examined by Li et al. in [31]. The results show that the proposed method presented in this paper can be applied in e-health services.
A method based on the integration of fused coupled chaotic map (FCCM) to encrypt ECG signals has been proposed in Pandey et al.’s study [32]. The FCCM generates four different chaotic maps based on control parameters. The effective-ness of the proposed approach has been confirmed using various analyzes, including key sensitivity, key space, correlation evaluation, and histogram analysis.
In the research of Sangavi et al. [33], the method (MMIE) using Rossler dynamical system and sine map (Sine) for medical image encryption has been proposed. For the method, statistical metrics such as histogram analysis, entropy, differential, and correlation have been computed.
Presented by researchers in [34] is an improved technique for medical image encryption based on discrete wavelet transform (DWT), discrete cosine transform (DCT) and, singular value decomposition (SVD) transformations. Additionally, the security of the method has been improved by exploiting 2D logistic maps based on chaotic cryptography.
In Zheng et al.’s [35] study, the data encryption standard (DES) and elliptic curves cryptography (ECC) algorithms are used for EEG signals. The results reveal that the ECC method has superior performance in encryption and decryption of EEG signals.
In this paper, a novel medical images encryption based on synchronization of hyper chaotic systems is presented. The main purpose of providing synchronization methods is to ad-just the controller parameters in chaotic systems properly. In order to synchronize, various control methods such as adaptive control [47,48], sliding mode [49,50], fuzzy [51,52], optimal [53,54], predictive [55,56], etc. have been proposed so far. Chaotic synchronization is such that the state paths of the following system follow the state paths of a base system [57,58].
Due to the numerous applications of synchronization, various other methods have been proposed that have different applications based on the existence of chaos and uncertainty in them. The following are some of the newest synchronization methods. In [59], the idea of synchronization based on the sliding mode control method and developing the advanced encryption standard algorithm is presented. Synchronization of chaotic systems based on additional Fourier modes is one of the newest methods proposed by [60]. Ding et al. [38] performed the encryption based on fractional-order Henon chaotic map and the discrete-time wavelet transform for a hyper-chaotic system. A new encryption idea based on the discrete-time sliding mode method and the improved El-Gamal encryption system is presented in [61]. Synchronization with the optimized proportional integral derivative (PID) controller with particle swarm optimization (PSO) algorithm is another new method in this field that has been introduced by [62].
In this paper, the main idea of medical image encryption is based on robust adaptive control. Also, the synchronization method is multi-mode. The multi-state synchronization with the presence of uncertainty, disturbance, unknowns, and variable parameters has been carried out. The proposed synchronization method is used alongside the Chen hyper chaotic system to encrypt medical images. In the provided example, one Chen hyper-chaotic system is considered as the master system, and two Chen hyper chaotic systems are considered as the slave systems. According to the simulations, the robustness of the proposed method and the convergence of all error types to zero verify its capabilities. In the following, taking the advantage of masking technique, it is possible to use the proposed synchronization method for medical images encryption. For encryption, first, a variety of standard benchmark images have been employed. Then, the chest CT and X-ray images of COVID-19 patients with various noises have been examined. Also, to accurately verify the performance of the proposed method, a variety of important cryptographic parameters including histogram analysis, correlation, differential attack, PSNR and information entropy of CT images with different noises have been calculated to represent the performance of the proposed method correctly. The results reveal that the synchronization technique of the proposed hyper-chaotic system for encrypting standard images and medical images has shown successful results and can therefore be adopted as an encryption section in a COVID-19 detection telemedicine system. The proposed approach is capable of encrypting other medical images.
In the following, other sections of the paper are presented. The second part is devoted to the theories of the proposed synchronization method. In the Section 3, the idea of encrypting medical images and standard benchmarks based on the chaotic masking method is expressed. The implementation of the proposed synchronization method on the Chen hyper-chaotic system is introduced in Section 4. The types of statistical parameters to show the efficiency of the proposed synchronization method in the encryption of different images are provided in the Section 5. The Section 6 is devoted to the results of the proposed method in encrypting medical images and standard benchmarks. The advantages and disadvantages of the proposed method are described in Section 7. Finally, the discussion and conclusion are presented in Section 8.

2. Formulation of the Theorem

This section first describes the problem of synchronization of multi-mode transitive chaotic systems. The multiple loop synchronization is then checked. In both cases, the adaptive laws and controllers are designed using the adaptive control method.

2.1. Comparative Synchronization between More Response Systems and a Master System

Figure 1 shows the synchronization model between a master system and several slave systems.
The master chaotic system is defined as follows [48]:
x ˙ 1 = f 1 ( x 1 ) + H 1 ( x 1 ) θ 1 ( t )
where x 1 ( t ) = [ x 11 , x 12 , , x 1 n ] T is the state vector of the system. f 1 ( x 1 ( t ) ) = [ f 11 , f 12 , , f 1 n ] T is the continuous function, H 1 ( x 1 ( t ) ) = [ H 11 , H 12 , , H 1 n ] T is matrix function and θ 1 ( t ) = [ θ 11 , θ 12 , , θ 1 n ] T is the main parameters of the master system, which are piecewise constant.
There are (N − 1) chaotic slave systems with control function as given by (2) [63]:
x ˙ i = f i ( x i ) + H i ( x i ) θ i ( t )   i = 2 ,   3 ,   ,   N
where x i ( t ) = [ x i 1 , x i 2 , , x i n ] T is the state vector of the ith system, f i ( x 1 ( t ) ) = [ f i 1 , f i 2 , , f i n ] T   is the continuous function, H i ( x i ( t ) ) = [ H i 1 , H i 2 , , H i n ] T is matrix function, θ i ( t ) = [ θ i 1 , θ i 2 , , θ i n ] T are the main parameters which are piecewise constant of the ith slave system, and u i 1 ( t ) = [ u i 1.1 ( t ) , u i 1.2 ( t ) , , u i 1. n ( t ) ] T , is the control function of ith system. Therefore, according to Equations (1) and (2), synchronization of the chaotic system with the control function is as follows:
x 1 ˙ = f 1 ( x 1 ) + H 1 ( x 1 ) θ 1 ( t ) x 2 ˙ = f 2 ( x 2 ) + H 2 ( x 2 ) θ 2 ( t ) + u 1 ( t ) x N ˙ = f N ( x N ) + H N ( x N ) θ N ( t ) + u N 1 ( t )
In the form of multiple synchronization mode, the synchronization error is defined as follows:
e i 1 ( t ) = x i ( t ) x 1 ( t )   i = 2 , 3 , , N  
Definition 1.
Chaotic system (3) using the controllers u i ( t ) and unknown parameters have multi-state transformative synchronization if the following condition is held:
lim t || e i 1 ( t ) || = l i m t   || x i ( t ) x 1 ( t ) ||   0   i = 2 , 3 , , N  
where the error dynamic defined as follows:
e ˙ i 1 ( t ) = f i ( x i ) f 1 ( x 1 ) + H i ( x i ) θ i ( t ) H 1 ( x 1 ) θ 1 ( t ) + u i 1   i = 2 , 3 ,   , N 1
is met, then the adaptive transmission synchronization between N chaotic systems with unknown parameters is realized.
The design of controllers and adaptive laws to achieve the above goal is based on Lyapunov’s function, and the synchronization under the mode of transmission synchronization is performed. The control law for u 1 ( t ) , u 2 ( t ) , , u N 1 ( t ) can be designed as follows:
u i 1 ( t ) = f i ( x i ) + f 1 ( x 1 ) H i ( x i ) θ ^ i ( t ) + H 1 ( x 1 ) θ ^ 1 ( t ) + K i 1 e i 1   i = 2 , 3 ,   , N 1  
in which:
K i 1 = K i 1 1 + K i 1 2 , ( K i 1 1 ) T = K i 1 1 ,   K i 1 2 = d i a g ( k i 1.1 . ,   k i 1.2 . , ,   k i 1. n . ) , k i 1 , j . > 0 ,   j = 1 , 2 , , n
Therefore, error dynamics are given as follows:
e ˙ i 1 ( t ) = H i ( x i ) θ ˜ i ( t ) H 1 ( x 1 ) θ ˜ 1 ( t ) + K i 1 e i 1   i = 2 , 3 ,   , N 1
where θ ^ i ( t ) is an approximation of θ i ( t ) and θ ˜ i ( t ) = θ i ( t ) θ ^ i ( t ) is an approximation error.
Theorem 1.
Transmission synchronization of chaotic systems (3) with controls (5), error dynamics (6) and update laws (9a), (9b), (10a), and (10b) is guaranteed.
Proof. 
Consider the following Lyapunov function:
V ( e . θ ˜ ) = 1 2 ( i = 2 N [ e i 1 T e i 1 + θ ˜ i T θ ˜ i ] + θ ˜ 1 T θ ˜ 1 )
The derivative of Lyapunov’s function (7) is obtained according to Equation (8)
d d t V = i = 2 N [ e i 1 T ( H i ( x i ) θ ˜ i H 1 ( x 1 ) θ ˜ 1 + 1 2 ( K i 1 + K i 1 T ) e i 1 ) + θ ˜ i T θ ˜ ˙ i ] + θ ˜ 1 T θ ˜ ˙ 1
Therefore, the parameters laws are chosen as follows:
θ ˜ ˙ i = ( H i ( x i ) e i 1 + σ i θ ˜ i )   σ i > 0   i = 2 , 3 , , N
θ ˜ ˙ 1 = i = 2 N 1 H 1 ( x 1 ) e i 1 σ 1 θ ˜ 1   σ 1 > 0
If θ i are constant then θ ˙ i = 0 and the parameters estimation laws are determined as follows:
θ ^ ˙ i = θ ˜ ˙ i = H i ( x i ) e i 1 + σ i θ ˜ i   σ i > 0   i = 2 , 3 , , N  
θ ^ ˙ 1 = θ ˜ ˙ 1 = i = 2 N 1 H 1 ( x 1 ) e i 1 + σ 1 θ ˜ 1   σ 1 > 0
By setting the estimation laws in (8), we will have:
d d t V = i = 2 N [ e i 1 T K i 1 2 e i 1 σ i θ ˜ i T θ ˜ i ] σ 1 θ ˜ 1 T θ ˜ 1   φ   V
where φ = min i , j ( σ i , k i 1. j . ) > 0 and K i 1 2   are Hurwitz. So system (6) is stable and || e i ( t ) || 0. Therefore, obtain an adaptive synchronization between (N − 1) slave systems and a master system. □
Note 1.
If the systems are time varying: θ i = θ i ( t ) Equations (9a) and (9b) are still stablished, then in these cases the Theorem 1 is also established. Since the synchronization error converges to zero, the Equations (10a) and (10b) are converted as follows for relatively long periods of time:
θ ˜ ˙ i = σ i θ ˜ i   θ ˜ i ( t ) = θ ˜ i ( 0 ) exp ( σ i t ) 0   σ i > 0   i = 1 , 2 , 3 , , N
Therefore, if θ i ( t ) is a vector with piecewise constant element and the appropriate time interval between changes in any element, the updating rule (10) is established and the parameters are estimated accurately. If | d d t θ i ( t ) | < s i < 1 then the update Equation (10) are also well established.

2.2. Circular Synchronization of Multiple Chaotic Systems with Unknown Parameters

Figure 2 displays circular multi-mode synchronization. In this type of synchronization, systems 1 and N have the role of master and slave respectively and other systems both have the role of slave and master at the same time (i-th system is master role of i + 1-th system), which makes it more complex.
According to Figure 2, the systems are synchronized with a circular structure. In this case, the first chaotic system is described by the following expression:
x ˙ 1 ( t ) = f 1 ( x 1 ) + H 1 ( x 1 ) θ 1 ( t )  
The (N − 1) chaotic systems are as follows:
{ x ˙ 1 ( t ) = f 1 ( x 1 ) + H 1 ( x 1 ) θ 1 ( t ) x ˙ 2 ( t ) = f 2 ( x 2 ) + H 2 ( x 2 ) θ 2 ( t ) + m 1 ( t ) x ˙ N ( t ) = f N ( x N ) + H N ( x N ) θ N ( t ) + m N 1 ( t )
The control input is
m i 1 ( t ) = [ m i 1.1 , m i 1.2 , , m i 1. n ] T
For N chaotic system described by Equation (13), if there are adaptive controllers, so that the error dynamics systems are defined as follows:
r i ( t ) = x i + 1 ( t ) x i ( t )   i = 1 , 2 , 3 , , N 1 r ˙ i 1 ( t ) = f i ( x i ) f i 1 ( x i 1 ) + H i ( x i ) θ i H i 1 ( x i 1 ) θ i 1 + m i 1 m i 2 i = 2 ,   3 ,   , N 1
and condition:
lim t || r i 1 ( t ) || = lim t || x i ( t ) x i 1 ( t ) ||   0   i = 2 , 3 , , N
is met then a circular synchronization between the N chaotic system and the unknown parameters is realized.
Theorem 2.
The following statements are used to synchronize in both circular and transmission modes.
(A) 
If there is a transmission and circular synchronization with the u(t) and m(t) controllers. Then:
i , j   :   l i m t || x i ( t ) x j ( t ) ||   0
(B) 
If transmission synchronization is established, circular synchronization is also realized and vice versa.
Proof. 
Assume that a transmission synchronization is established, so:   i   :   e i ( t ) 0 so:
  i , j : || x i ( t ) x j ( t ) || = || ( x i ( t ) x 1 ( t ) ) ( x j ( t ) x 1 ( t ) ) || || e i ( t ) || + || e j ( t ) || 0
If a circular synchronization is established (assuming i > j ):
  i , j : || x i ( t ) x j ( t ) || = || ( x i ( t ) x i 1 ( t ) ) + ( x i 1 ( t ) x i 2 ( t ) ) + + ( x j + 1 ( t ) x j ( t ) ) || || r i 1 ( t ) || + + || r j ( t ) || k = j i 1 || r k ( t ) || 0
 □
Proof. 
The following relationships between errors are established in two synchronization modes:
r i ( t ) = { e i ( t ) e i 1 ( t )   i = 2 , 3 , , N 1   e i ( t )   i = 1 ,       e i ( t ) = k = 1 i r k ( t )
If a transmission synchronization is established, e i ( t ) 0   therefore:
|| r i ( t ) || = { || e i ( t ) e i 1 ( t ) || || e i ( t ) || + || e i 1 ( t ) || 0 i 2 || e 1 ( t ) || 0 i = 1
So circular synchronization is established. Conversely, suppose a circular synchronization is established, r i ( t ) 0 therefore:
|| e i ( t ) || = || k = 1 i r k ( t ) || k = 1 i || r k ( t ) || 0
Therefore, transmission synchronization is also established. So the two types of synchronization are equivalent. □
Theorem 3.
The control law is same for transmission synchronization and circular synchronization. In other words: m i ( t ) = u i ( t ) .
Proof. 
Using the relationship between two types of errors:
e ˙ i ( t ) = e ˙ i 1 ( t ) + r ˙ i ( t )
Which is obtained by placing in Equations (2) and (14):
u i 1 ( t ) = u i ( t ) ( m i ( t ) m i 1 ( t ) )   i 2
For i = 1, we have:
e 1 ( t ) = r 1 ( t ) e ˙ 1 ( t ) = r ˙ 1 ( t ) m 1 ( t ) = u 1 ( t )
By placing in Equation (16), we have: m i ( t ) = u i ( t ) so the proof is completed. □

2.3. Synchronization with the Presence of Disturbance and Uncertainty in the System

In this case, the master and slave systems are accompanied by disturbance and uncertainty as follows:
{ x ˙ 1 = f 1 ( x 1 ) + H 1 ( x 1 ) θ 1 ( t ) + Δ f 1 ( x 1 ) + D 1 ( t ) x ˙ 2 = f 2 ( x 2 ) + H 2 ( x 2 ) θ 2 ( t ) + Δ f 2 ( x 2 ) + D 2 ( t ) + u 1 ( t ) x ˙ N = f N ( x N ) + H N ( x N ) θ N ( t ) + Δ f N ( x N ) + D N ( t ) + u N 1 ( t )
where it is assumed that uncertainties and disturbances are bounded but with unknown bound.
| Δ f i ( x i ) | γ i g i ( x i ) ,   | D i ( t ) |   d i   i = 1 , 2 , , N
where γ i and di are constant but indefinite and g i ( x i ) is a definite function and in the special case   g i ( x i ) = | x i | . The dynamics of the errors are as follows:
e ˙ i 1 ( t ) = f i ( x i ) f 1 ( x 1 ) + H i ( x i ) θ i ( t ) H 1 ( x 1 ) θ 1 ( t ) + Δ f i ( x i ) Δ f 1 ( x 1 ) + D i ( t ) D 1 ( t ) + u i 1 ( t )   i = 2 , 3 ,   , N 1
By defining the control function as follows:
u i 1 ( t ) = f i ( x i ) + f 1 ( x 1 ) H i ( x i ) θ ^ i ( t ) + H 1 ( x 1 ) θ ^ 1 ( t ) + K i 1 e i 1 + u ¯ i 1 ( t ) i = 2 , 3 ,   , N 1
where θ ^ i ( t ) is an estimate of θ i ( t ) and u ¯ i 1 ( t ) is part of the control function which is introduced below. By placing the control function in (18), the dynamics of errors are as follows:
e ˙ i 1 ( t ) = H i ( x i ) θ ˜ i ( t ) H 1 ( x 1 ) θ ˜ 1 ( t ) + Δ f i ( x i ) Δ f 1 ( x 1 ) + D i ( t ) D 1 ( t ) + K i 1 e i 1 + u ¯ i 1 ( t )   i = 2 , 3 ,   , N 1
Theorem 4.
The dynamic system of error Equation (20) under control Equation (19), Equation (25) and update laws Equations (29) and (31) is stable and synchronization errors are convergent to zero despite uncertainty and disturbance.
Proof. 
by defining the Lyapunov function as follows:
V = 1 2 ( V e + V θ + V γ + V d )
where:
V e = i = 2 N e i 1 T e i 1   V θ = i = 2 N θ ˜ i T θ ˜ i + θ ˜ 1 T θ ˜ 1 .
V γ = i = 2 N γ ˜ i 2 + γ ˜ 1 2   V d = i = 2 N d ˜ i 2 + d ˜ 1 2 .
By calculating the derivative of Lyapunov’s function:
d d t V = i = 2 N [ e i 1 T ( H i ( x i ) θ ˜ i H 1 ( x 1 ) θ ˜ 1 + Δ f i ( x i ) Δ f 1 ( x 1 ) + D i ( t ) D 1 ( t ) + 1 2 ( K i 1 + K i 1 T ) e i 1 + u ¯ i 1 ( t ) ) + θ ˜ i T θ ˜ ˙ i + γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ] + θ ˜ 1 T θ ˜ ˙ 1 + γ ˜ 1 γ ˜ ˙ 1 + d ˜ 1 d ˜ ˙ 1 .
Note that for θ ˜ i ( t ) the same laws (9) and (10) are established, so Equation (22) becomes the following:
d d t V = i = 2 N [ e i 1 T ( Δ f i ( x i ) Δ f 1 ( x 1 ) + D i ( t ) D 1 ( t ) + 1 2 ( K i 1 + K i 1 T ) e i 1 + u ¯ i 1 ( t ) ) + γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ] + d ˜ 1 d ˜ ˙ 1 i = 1 N σ i θ ˜ i T θ ˜ i
If Δ f i j   , D i j   ,   e i 1 j   and u ¯ i 1 j are the component j of the vectors   Δ f i ,   D i , e i 1 and u ¯ i 1 ( t ) then:
d d t V = i = 2 N j = 1 n e i 1 j ( Δ f i j Δ f 1 j + D i j D 1 j + u ¯ i 1 j ) + i = 2 N ( γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ) + i = 2 N e i 1 T K i 1 2 e i 1 + γ ˜ 1 γ ˜ ˙ 1 + d ˜ 1 d ˜ ˙ 1 i = 1 N σ i θ ˜ i T θ ˜ i
Therefore:
d d t V i = 2 N j = 1 n [ | e i 1 j | ( | Δ f i j | + | Δ f 1 j | + | D i j | + | D 1 j | ) + e i 1 j u ¯ i 1 j ] + i = 2 N ( γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ) + i = 2 N e i 1 T K i 1 2 e i 1 + γ ˜ 1 γ ˜ ˙ 1 + d ˜ 1 d ˜ ˙ 1 i = 1 N σ i θ ˜ i T θ ˜ i
The bounded condition of disturbance and uncertainty can be developed on the components Δ f i and D i as follows:
| Δ f i j | max j | Δ f i j | | Δ f i ( x i ) | γ i g i ( x i )
| D i j ( t ) | max j | D i j ( t ) | | D i ( t ) |   d i
Which we have placed in the following equation.
d d t V i = 2 N j = 1 n [ | e i 1 j | ( γ i   g i ( x i ) + γ 1   g 1 ( x 1 ) + d i + d 1 ) + e i 1 j u ¯ i 1 j ] + i = 2 N (   γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ) + i = 2 N e i 1 T K i 1 2 e i 1 + γ ˜ 1 γ ˜ ˙ 1 + d ˜ 1 d ˜ ˙ 1 i = 1 N σ i θ ˜ i T θ ˜ i
If u ¯ i 1 j ( t ) is defined as follows:
u ¯ i 1 j ( t ) = ( γ ^ i   g i ( x i ) + γ ^ 1   g 1 ( x 1 ) + d ^ i + d ^ 1 ) · sgn ( e i 1 j )  
Then:
d d t V i = 2 N j = 1 n [ | e i 1 j | ( γ ˜ i   g i ( x i ) + γ ˜ 1 g 1 ( x 1 ) + d ˜ i + d ˜ 1 ) ] + i = 2 N e i 1 T K i 1 2 e i 1 + i = 2 N (   γ ˜ i γ ˜ ˙ i + d ˜ i d ˜ ˙ i ) + γ ˜ 1 γ ˜ ˙ 1 + d ˜ 1 d ˜ ˙ 1 i = 1 N σ i θ ˜ i T θ ˜ i
By defining the estimation error: γ ˜ i = γ i γ ^ i ,   d ˜ i = d i d ^ i . The update laws are as follows:
γ ˜ i γ ˜ i + j = 1 n | e i 1 j | γ ˜ i   g i ( x i ) = α i γ ˜ i 2 γ ˜ ˙ i = g i ( x i ) j = 1 n | e i 1 j | α i γ ˜ i   i = 2 , 3 , , N .  
γ ˜ 1 γ ˜ ˙ 1 + i = 2 N j = 1 n | e i 1 j | γ ˜ 1 g 1 ( x 1 ) = α 1 γ ˜ 1 2 γ ˜ ˙ 1 = g 1 ( x 1 ) i = 1 N j = 1 n | e i 1 j | α 1 γ ˜ 1 .
d ˜ i d ˜ ˙ i + i = 2 N j = 1 n | e i 1 j | d ˜ i = β i d ˜ i 2 d ˜ ˙ i = j = 1 n | e i 1 j | β i d ˜ i   i = 2 , 3 , , N .  
d ˜ 1 d ˜ ˙ 1 + i = 2 N j = 1 n | e i 1 j | d ˜ 1 = β 1 d ˜ 1 2 d ˜ ˙ 1 = i = 1 N j = 1 n | e i 1 j | β 1 d ˜ 1 .
where α i   and β i   are positive values. By placing the above update laws in Equation (28) we have:
d d t V i = 2 N 1 2 e i 1 T K i 1 2 e i 1 i = 1 N ( α i γ ˜ i 2 + β i d ˜ i 2 ) i = 1 N σ i θ ˜ i T θ ˜ i < μ V  
where: μ = min i . j ( α i , β i ,   σ i , k i 1. j . ) > 0 .
Therefore, since K i 1 2   are Hurwitz, the stability of the system is proved. Also, the convergence of synchronization errors to zero is guaranteed despite the uncertainty and disturbance. The laws for updating estimates will also be as follows:
d d t γ ^ i = g i ( x i ) j = 1 n | e i 1 j | + α i γ ˜ i   i = 2 , 3 , , N  
d d t γ ^ 1 = g 1 ( x 1 ) i = 1 N j = 1 n | e i 1 j | + α 1 γ ˜ 1
d d t d ^ i = j = 1 n | e i 1 j | + β i d ˜ i   i = 2 , 3 , , N
d d t d ^ 1 = i = 1 N j = 1 n | e i 1 j | + β 1 d ˜ 1  
Note 2.
Theorems 2 and 3 are established in spite of uncertainty and disturbance, because the nature of proving them does not depend on the existence or non-existence of uncertainty and disturbance. Therefore, the problem of circular synchronization is also solved in the presence of uncertainty and disturbance.
Note 3.
If systems are time-varying: θ i = θ i ( t ) Equations (9a) and (9b) are still established, so in these circumstances the Theorem 4 is also established. Therefore, as stated in note 1: circular synchronization of time-varying chaotic systems is also proven in the presence of uncertainty and disturbance.
Note 4.
The following relation can be used to make the control function continuous:
u ¯ i 1 j ( t ) = ( γ ^ i g i ( x i ) + γ ^ 1 g 1 ( x 1 ) + d ^ i + d ^ 1 ) · tan h ( λ e i 1 j ) )   λ 10  
Note 5.
If uncertainty forms are: | Δ f i ( x i ) | γ i | x i |   , i = 1 , 2 , , N then in update and control laws, it should be   g i ( x i ) = | x i | .
Note 6.
The final control function is as follows:
u i 1 ( t ) = f i ( x i ) + f 1 ( x 1 ) H i ( x i ) θ ^ i ( t ) + H 1 ( x 1 ) θ ^ 1 ( t ) + K i 1 e i 1 ( γ ^ i g i ( x i ) + γ ^ 1 g 1 ( x 1 ) + d ^ i + d ^ 1 ) · tan h ( λ e i 1 ) )

3. Application in Secure Communication Based on Chaotic Masking

In this scheme, the information signal is added to the linear composition of the master chaotic system state signals. Assuming that S(t) is the information signal carried by the master system, then the transmission message q ( t ) is defined as follows [64,65]:
q ( t ) = S ( t ) + i = 1 n η i z i ( t )  
In the given relation, z i ( t ) denotes i-th component of the master system and the signal q ( t ) is masked using the chaotic signal. This signal is transmitted from the sender to the receiver through the communication channel. Adopting the proposed controller, multi-state chaotic synchronization is conducted on one of its states. Using Equation (34), the received signal can be recovered [66]:
S ^ ( t ) = q ( t ) i = 1 n η i y i ( t ) .
S ^ ( t ) = S ( t ) + i = 1 n η i z i ( t ) i = 1 n η i y i ( t )
= S ( t ) + i = 1 n η i ( z i ( t ) y i ( t ) ) = S ( t ) + i = 1 n η i e i ( t ) S ( t )  
The multistate system for image masking is shown in Figure 3. According to Figure 3, the model entails a master system and several slave systems. Consequently, here the message is first transmitted to the master system. In the next step, the image is encrypted and synchronized using the slave system. Finally, after synchronization in the receiver, the image is decrypted and the original messages are recovered.
Utilizing the masking method, encryption and decryption of two signals are fulfilled using the proposed multi-state synchronization method. The simulation is executed with MATLAB software. Decryption and encryption are implemented after the synchronization of the master and slave chaotic systems in the presence of disturbance and uncertainty. To that end, the given image is first converted into a signal. In the next step, this signal is added to one of the components of the master system state vector and is transmitted to the receiver through a communication channel. In the receiver part, first using the synchronization, the message signal is extracted from the chaotic signal and then converted into an image.

4. Implementation of the Proposed Synchronization Method on Chen Hyper-Chaotic System

A Chen hyper- chaotic system is considered as a master system and two Chen hyper chaotic systems are considered as slaves. They are defined as follows:
{ x ˙ 11 = θ 11 ( x 12 x 11 ) + x 14 x ˙ 12 = θ 12 x 11 x 11 x 13 + θ 13 x 12 x ˙ 13 = x 11 x 12 θ 14 x 13 x ˙ 14 = x 12 x 13 + θ 15 x 14
f 1 ( x 1 ) = [ x 14 x 11 x 13 x 11 x 12 x 12 x 13 ]   H 1 ( x 1 ) = [ x 12 x 11 0 0 0 0 x 11 0 0 0 x 12 0 0 0 0 x 13 0 0 0 0 x 14 ]  
Slave systems:
{ x ˙ 21 = θ 21 ( x 22 x 21 ) + x 24 + u 11 x ˙ 22 = θ 22 x 21 x 21 x 23 + θ 23 x 22 + u 12 x ˙ 23 = x 21 x 22 θ 24 x 23 + u 13 x ˙ 24 = x 22 x 23 + θ 25 x 24 + u 14
f 2 ( x 2 ) = [ x 24 x 21 x 23 x 21 x 22 x 22 x 23 ]   H 2 ( x 2 ) = [ x 22 x 21 0 0 0 0 x 21 0 0 0 x 22 0 0 0 0 x 23 0 0 0 0 x 24 ]
{ x ˙ 31 = θ 31 ( x 32 x 31 ) + x 34 + u 21 x ˙ 32 = θ 32 x 31 x 31 x 33 + θ 33 x 32 + u 22 x ˙ 33 = x 31 x 32 θ 34 x 33 + u 23 x ˙ 34 = x 32 x 33 + θ 35 x 34 + u 24
f 3 ( x 3 ) = [ x 34 x 31 x 33 x 31 x 32 x 32 x 33 ]
H 3 ( x 3 ) = [ x 32 x 31 0 0 0 0 x 31 0 0 0 x 32 0 0 0 0 x 33 0 0 0 0 x 34 ]  
In all examples, the parameters are set as follows:
K i 1 2 = D i a g ( 10 , 10 , 10 , 10 ) ,   σ i = 10 ,   α i = 20 , β i = 20 ,   i = 1 , 2 , 3 , 4
Case A. Multi-mode synchronization with time-varying parameters (without disturbance and uncertainty).
For parameter values, all systems are chaotic as follows:
θ i 1 = 35   θ i 2 = 3   θ i 3 = 12   θ i 4 = 0.3   θ i 5 = 7 ,   i = 1 , 2 , 3 , 4 , 5
In the simulation of the basic conditions, the master system and two slave systems are selected as follows:
(x11(0), x12(0), x13(0), x14(0)) = (10, 10, 10, 10)
(x21(0), x22(0), x23(0), x24(0)) = (2, 2, 2, 2)
(x31(0), x32(0), x33(0), x34(0)) = (3, 3, 3, 3)
The initial values of the parameters are assumed to be as follows:
θ ^ 1 ( 0 ) = [ 3 3 3 3 ]   θ ^ 2 ( 0 ) = [ 3 3 3 3 ]   θ ^ 3 ( 0 ) = [ 4 4 4 4 ]
The main parameters over time are in the form of variable steps and as follows:
For a time from zero to 2 s:
θ11 = 35, 12θ = 7, 21θ = 35, 22θ = 7, 31θ = 35, 32θ = 7
For a time from 2 to 6 s:
13θ = 12, 14θ = 3, 15θ = 0.3, 23θ = 12, 24θ = 3, θ25 = 0.3, 33θ = 12, 34θ = 3, 35θ = 0.3,
For a time from 2 to 5 s:
11θ = 33, 12θ = 6, 21θ = 32, 22θ = 5, 31θ = 34, 32θ = 5,
For a time from 6 to 10 s:
13θ = 13 14θ = 1, 15θ = 0.1, 23θ = 12, 24θ = 5, θ25 = 0.2, 33θ = 10, 34θ = 2, 35θ = 0.2,
Figure 4 shows that synchronization errors in presence of uncertainty and disturbances are quickly converged to zero. On the other hand, parameter changes have little effect on synchronization errors and are quickly attenuated.
Figure 5 shows that control laws are continuous functions prevent the phenomenon of chattering.
As displayed in Figure 4, synchronization of chaotic systems as the main objective is well carried out and synchronization errors converge to zero despite changes in parameters.
As shown in Figure 5, the obtained control signals are continuous. In addition, changes in the parameters are met with a rapid reaction in the control signals to converge the synchronization errors to zero.
Case B. Multi-mode synchronization with time-varying parameters despite disturbance and uncertainty. In this case, disturbance and uncertainties are applied to master and slave systems as follows:
Δ f 1 = [ 0.2 x 11 sin ( x 11 + x 14 ) 0.1 x 13 sin ( x 11 + 2 x 12 ) 0.2 ( x 14 + x 12 ) cos ( x 11 x 13 ) x 12 sin ( x 14 + x 12 ) ]
Δ f 2 = [ 0.5 x 22 sin ( x 21 x 22 ) 0.2 x 24 sin ( x 21 3 x 22 ) 0.5 ( x 24 x 22 ) sin ( x 21 + 2 x 23 ) 0.4 ( x 21 x 22 ) cos ( x 24 x 22 ) ]
Δ f 3 = [ 0.3 x 32 c o s ( x 31 + x 34 ) 0.9 x 31 s i n ( x 31 x 32 ) 0.4 x 32 s i n ( x 31 + x 33 ) 0.3 ( x 31 x 32 ) s i n ( x 34 2 x 32 ) ]
| Δ f i j |   max j | Δ f i j | | Δ f i ( x i ) | γ i | x i |   γ 1 = 1   γ 2 = 0.5   γ 3 = 0.9
D 1 = [ 0.2 sin ( π 3 t ) 0.1 sin ( π 10 t ) 0.25 sin ( π 4 t ) 0.3 sin ( π 5 t ) ]   D 2 = [ 0.25 sin ( π 2 t ) 0.15 sin ( π 20 t ) 0.20 sin ( π 10 t ) 0.35 sin ( π 3 t ) ]   D 3 = [ 0.3 sin ( π 30 t ) 0.2 sin ( π 20 t ) 0.15 sin ( π 10 t ) 0.1 sin ( π 25 t ) ]
| D i j ( t ) | max j | D i j ( t ) | | D i ( t ) |   d i   d 1 = 0.3   d 2 = 0.35   d 3 = 0.3
In this case, the change in parameters is similar to case (A). Figure 6 shows the estimation of variable parameter of the system. It can be noted that, estimating time-varying parameters in the presence of uncertainty and disturbance is correctly done.
Figure 6 shows the multi-mode synchronization errors despite the structural uncertainty and disturbance. We have obtained good performance as the synchronization errors are close to zero. Multi-mode synchronization errors for both master and slave systems, despite disturbance and uncertainties, have quickly reached to zero. Subsequently, change in the parameter causes small errors, which quickly converges to zero (Figure 7).
As can be seen in Figure 6, the synchronization of chaotic systems as the main aim is well done, and the synchronization errors converge to zero despite changes in parameters, presence of disturbance, and uncertainty.
Case C. Multi-mode synchronization, taking into account the disturbance and uncertainties bounded with the function:
For the case C   | Δ f i j |   max j | Δ f i j | | Δ f i ( x i ) | γ i g i ( x i )
Δ f 1 = [ 2 x 11 2 cos ( x 11 + x 14 ) x 14 2 sin ( x 11 + 2 x 12 ) 0.2 sin ( x 11 x 13 ) x 13 2 cos ( x 14 + x 12 ) ] .   g 1 ( x 1 ) = | x 1 | 2 + 0.1 | x 1 |
Δ f 2 = [ 0.5 x 24 2 sin ( x 21 x 22 ) 0.2 x 23 2 cos ( x 21 3 x 22 ) 0.5 ( x 24 2 + 0.6 x 22 2 ) cos ( x 21 + 2 x 23 ) 0.4 ( x 24 2 + 0.6 x 22 2 ) sin ( x 24 x 22 ) ] .   g 2 ( x 2 ) = | x 2 | 2  
Δ f 3 = [ 0.3 x 31 2 c o s ( x 31 + x 34 ) 0.5 x 33 2 s i n ( x 31 x 32 ) 0.4 ( x 31 x 32 ) 2 s i n ( x 31 + x 33 ) 0.3 s i n ( x 34 2 x 32 ) ] .   g 3 ( x 3 ) = | x 3 | 2 + 1.2 | x 3 |
In this case, the change in parameters is similar to case (A), the disturbances are similar to case (B) and uncertainties are not bounded with norm of states but bounded with a function of state norms. This increases the control problem complexity. In Figure 8, the multi-mode synchronization error in case (C) shows that despite the time-varying parameters, the functional bounded uncertainty and disturbance converge to zero.
As displayed in Figure 8, the synchronization of chaotic systems as the main objective is well fulfilled, and in spite of changes in parameters, the presence of disturbance, and uncertainty, synchronization errors are converged to zero.
Figure 9 shows that the estimation errors of uncertainties are more fluctuating than in case (B) but are close to zero. Estimation errors of disturbance and uncertainty bounds have reached zero with low fluctuations. The uncertainties are complex functions. The peaks appear in the figures only when the parameters have changed. However, over a short period of time, the magnitude of curves are decreased and eventually approached to zero. The proposed approach is able to reject uncertainties bounded with polynomial of state norms and unknown non-negative coefficients.
The chaotic behavior of systems is shown in Figure 10. The above figures exhibit the chaotic behavior of the master (green) and slave systems (blue and red). The presence of disturbance and uncertainty as well as changes in system parameters make the system behavior more complex.
Chaotic behavior leads to the complicatedness of secure communications method and the possibility of detecting it is so low because in chaotic systems, in addition to uncertainty and disturbance, its parameters are assumed to be variable, which promotes the security of the approach.

5. Statistical Metrics

5.1. Histogram Analysis

In image processing, the histogram represents the distribution of the pixel values of an image. Histogram variance is employed in image encryption. In this manner, low variance values indicate higher uniformity in encrypted images. The variance of histograms is expressed by Equation (36) [24]:
V a r ( I ) = 1 n 2 i = 1 n j = 1 n 1 2 ( I i I j ) 2  
According to Equation (36), I denotes the vector of histogram values, Ii and Ij indicate the number of pixels whose gray values are equal to i and j, respectively [24].

5.2. Correlation Analysis

Correlation distributions and correlation coefficients play a crucially significant role in the analysis of encrypted images, which are discussed in this section. The following statements are taken into account to calculate these values [24]:
{ E ( x ) = 1 N i = 1 N x i D ( x ) = 1 N i = 1 N ( x i E ( x ) ) 2 C o v ( x , y ) = 1 N i = 1 N ( x i E ( x ) ) ( y i E ( y ) ) γ x , y = C o v ( x , y ) D ( x ) D ( y )

5.3. Differential Attack Analysis

Pixel change rate (NPCR) and unified average changing Intensity (UACI) are two other important measurable parameters in image encryption, whose relationships of each are shown below [24]:
D ( i , j ) = { 1 I e 1 I e 2 0 I e 1 = I e 2
N P C R = i = 1 M j = 1 N D ( i , j ) N × M × 100 %
U A C I = i = 1 M j = 1 N | I e 1 ( i , j ) I e 2 ( i , j ) | 255 N × M × 100 %
Evaluation of encrypted images by NPCR and UACI parameters is depicted in Table 1, Table 2, Table 3, Table 4 and Table 5. The results reveal that a swift change in the original image leads to a change in the ciphered image. This signifies that the proposed scheme has a high ability to resist differential attack [24].

5.4. PSNR Analysis

PSNR (peak signal to noise ratio) is the ratio of peak signal power to noise power. An encryption method achieves successful performance when the encrypted image has a low PSNR. PSNR relationships are described below [24]:
{ M S E = 1 M * N i = 0 m 1 j = 0 n 1 [ I ( i , j ) K ( i , j ) ] 2 P S N R = 20 * log 10 ( M A X I ) 10 * log 10 ( M S E )
where MSE is the mean squared error and MAXI is the maximum possible pixel value of the image [24].

5.5. Information Entropy Analysis

Information entropy is another image analysis procedure in cryptographic applications and is defined as follows [24]:
H ( s ) = i = 0 2 N 1 P ( s i ) log P ( s i )

6. Experiment Results

The results of images encryption based on the synchronization method of the proposed chaotic systems are yielded in this section. In the following, first, the images used to perform the experiments are introduced. Then, the results of the proposed approach are expressed by application.

6.1. Image Benchmarks

In this paper, in order to evaluate the proposed cryptographic method, multiple medical images and benchmarks have been adopted. The images tested in this section contain a variety of standard benchmark images and medical images including CT and X-ray. In Figure 11, 10 images of standard benchmarks for conducting experiments are displayed. Also, shown in Figure 12 are X-ray images of COVID-19 patients. Finally, 10 CT images of patients with COVID-19 are used for the final tests, as shown in Figure 13.

6.2. Simulation

In this section, the results of the proposed synchronization method of chaotic systems for various images encryption are discussed. All images utilized in this research have a size of 300 × 300 in png format. Figure 14, Figure 15 and Figure 16 show encrypted images of standard benchmark, X-ray, and CT patients with COVID-19, respectively.
The results of the proposed chaotic synchronization method for famous benchmark images encryption are shown in Figure 14. As it is manifest, the ten selected images have different contrasts, and the results displayed in Figure 14 verify the high effectiveness of the proposed techniques in the encryption of these images.
In Figure 15, X-ray images of patients with COVID-19 for cryptography are exploited. According to Figure 13, it can be seen that the synchronization method of the proposed chaotic system also achieves highly successful outcomes in X-ray images encryption.
Figure 16 demonstrates the various CT encrypted images. According to the figure, the CT images are first fed to the input of the proposed chaotic synchronization method. CT images encrypted are then generated. Finally, reconstructed CT images are represented. According to Figure 16, it can be seen that the cryptographic procedure based on synchronization of the proposed chaotic system achieves successful results.
As can be observed in Figure 16, the recommended method has also been very successful in CT images encryption of people with COVID-19. Additionally, with a conclusion, it is perceived that the cryptographic method based on synchronization of chaotic systems presented in this paper for various medical images and standard benchmarks achieves satisfactory results.
In Figure 17, a number of standard X-ray and CT images are randomly selected and their histogram diagrams are drawn. According to Figure 17, histograms related to input and decoded images are displayed.
In Figure 17, histogram diagrams of the number of the original input images and encrypted images are illustrated. The results show that the histograms of the input and decoded images are very close to each other. This means that the proposed synchronization method of the chaotic system has performed highly successfully in encrypting and recovering various images. For further experiments, Gaussian noises with different variances have been applied to CT images. Then, the types of statistical parameters of Section 5 are calculated for them. CT images have been corrupted by Gaussian noise with 0.001, 0.003, 0.006, and 0.009 variances, and the results of their statistical analysis are shown in Table 1, Table 2, Table 3, Table 4 and Table 5.
In Table 1, the results of the proposed methodology have been yielded on CT images without noise. As can be seen, this method has statistically satisfactory results in CT cryptography of patients with COVID-19.
In Table 2, the results of the synchronization method of chaotic systems for CT images with Gaussian noise with 0.001 variance are shown. The results in this table prove that the proposed method for encrypting CT images with Gaussian noise with 0.001 variance has been successful.
Table 3 indicates the statistical results of the proposed synchronization approach for CT images cryptography of COVID-19 patients with Gaussian noise (variance 0.003). It can be perceived that the proposed technique has been functioned successfully in encrypting these noise-corrupted CT images.
The findings of CT images encryption corrupted by Gaussian noise with 0.006 variance are given in Table 4. It is obvious that the proposed scheme for CT images for Gaussian noise with 0.006 variance compared to the previous modes has been able to report good results.
Finally, in Table 5 the results of the proposed method cryptography with CT images for Gaussian noise with 0.009 variance have been illustrated. From the results in this table, it can be seen that the proposed synchronization method is also robust to this type of Gaussian noise in CT images encryption.

7. Advantages and Disadvantages

The proposed method has the ability to deal with uncertainty, disturbance, and changes in system parameters. On the other hand, the resulting control functions are continuous, which prevents the chattering problems. The only problem with our proposed algorithm is the large magnitude of few control functions. In the future, we intend to develop the proposed method for synchronization of chaotic and hyper-chaotic fractional-order systems. Using an objective function as the total of control norms can reduce the magnitude of control functions.
The proposed method has strengths from various perspectives, the most important of which are discussed in the following. Ensuring the synchronization error converges to zero, which is the most mattering objective in the synchronization issue. These errors converge to zero in a short period of time, despite uncertainty and disturbance, as well as changes in the master and slave systems. Providing a continuous control function that prevents the occurrence of the chattering phenomenon. Proving the equivalence of control functions in two synchronization methods is another significant result. Finally, according to the conducted analysis, the resulting procedure has good robustness against noise, which is one of the most important outcomes of this paper.
Also, the disadvantages of the proposed method are expressed in this paragraph. In the proposed approach, there are some constraints, which incorporate: large control signal that increases the cost of synchronization. On the other hand, there is no physical factor and a significant delay in the master and slave systems. Ideal communication channel, as well as parameter changes that are considered stepwise, are other limitations of this method.

8. Discussion and Conclusions

Today, telemedicine is employed in a variety of medical applications. In all telemedicine systems, the issue of information security is vitally significant. To this end, various approaches for the information security of these systems have been presented so far [12]. Cryptography based on chaotic theory is one of the best approaches to protect the information in secure communications.
Encryption techniques based on chaos theory in medicine have also attracted the attention of many. In the control field, it has been proven that synchronization techniques of hyper chaotic systems to enhance information security have reported remarkably successful outcomes [54,55,56,66]. In this paper, for the first time, a novel multi-mode synchronization of a chaotic system has been utilized to encrypt medical images of CT and X-ray for COVID-19 patients. This method can be employed as an encryption section in the various telemedicine systems.
In this paper, the multi-mode synchronization of two hyper chaotic systems in two forms of multi slave, and one master system synchronization and circular synchronization are examined. Firstly, it is proved that two forms of synchronizations are equivalent and control laws are exactly equal. The adaptive control method and the definition of appropriate Lyapunov function are used for synchronization. The convergence of synchronization errors to zero is ensured despite variable parameters with time and disturbance, and uncertainty is bounded with a known function. Adaptive laws to estimate the system-time-varying parameters, disturbance, and uncertainty bounds, are also determined to ensure that the system stability is guaranteed. For preventing the phenomenon of chattering, the law of controls is established as a continuous function. To test the effectiveness of the proposed method, simulations are performed in the presence of structural uncertainty and disturbance. The results show that the proposed controller is able to perform efficiently by yielding zero errors for synchronization, estimations of disturbance, and uncertainty bounds.
After implementing the proposed multi-mode synchronization method, subsequently, medical image encryption of COVID-19 patients has been fulfilled. To verify the efficiency of the recommended method, CT images were corrupted by Gaussian noise with different variances and the statistical analysis results revealed the effectiveness of the proposed approach.
The advantages of the method involve synchronization in the presence of disturbance and uncertainty, which indicates its robustness. Using the Lyapunov function, which guarantees the convergence of all types of errors, including synchronization error and the estimation error to zero. On the other hand, the variables of the system are another great issue that makes the synchronization problem more complicated, and finally, the images are encrypted with high security, which is very difficult to detect. In Table 6, the advantages of the proposed method are compared with other researches in this field.
According to Table 6, the proposed method has more advantages than other proposals in the field of medical data encryption using chaotic methods.
For future work, the proposed method can be utilized to encrypt EEG signals of epileptic seizures [67,68], attention deficit hyperactivity disorder (ADHD) [69,70], schizophrenia [71,72], and multiple sclerosis [73]. Additionally, the technique presented in this study can be adopted to encrypt MRI [74,75], X-ray [76], mammography [77,78], and other medical data [79,80,81,82,83,84].
In the continuation of this paper in order to develop and eliminate deficiencies such as applying optimal control to synchronize with less control effort and considering the important factor of delay in master and slave systems as well as method extension for fractional-order systems in which case the complexity of the control method increases and the security of the method is enhanced. For another future works, due to the fact that medical data has uncertainty, neuro fuzzy, fuzzy type1 or fuzzy type 2 systems can be used [85,86,87,88,89,90].

Author Contributions

Conceptualization, A.A.K.J. and A.S.; methodology, A.Z., A.A.K.J., A.S. and J.M.G.; software, A.A.K.J. and A.S.; validation, A.Z. and A.A.K.J.; formal analysis, A.S., R.A., M.K. and N.G.; resources, A.Z. and A.A.K.J.; writing—original draft preparation, A.A.K.J., A.S., R.A. and M.J.; writing—review and editing, A.S., A.Z., M.K., J.M.G., N.G. and M.J.; visualization, A.Z., A.A.K.J., and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the MINECO/ FEDER under the RTI2018-098913-B100 CV20-45250 and A-TIC- 080-UGR18 projects.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fadhil, A. Beyond patient monitoring: Conversational agents role in telemedicine & healthcare support for home-living elderly individuals. arXiv 2018, arXiv:1803.06000. [Google Scholar]
  2. Wijesinghe, I.; Gamage, C.; Perera, I.; Chitraranjan, C. A smart telemedicine system with deep learning to manage diabetic retinopathy and foot ulcers. In Proceedings of the 2019 Moratuwa Engineering Research Conference (MERCon), Moratuwa, Sri Lanka, 3–5 July 2019; pp. 686–691. [Google Scholar]
  3. Leite, H.; Hodgkinson, I.R.; Gruber, T. New development: ‘Healing at a distance’—Telemedicine and COVID-19. Public Money Manag. 2020, 40, 483–485. [Google Scholar] [CrossRef] [Green Version]
  4. Portnoy, J.; Waller, M.; Elliott, T. Telemedicine in the era of COVID-19. J. Allergy Clin. Immunol. Pract. 2020, 8, 1489–1491. [Google Scholar] [CrossRef]
  5. Angelucci, A.; Kuller, D.; Aliverti, A. A home telemedicine system for continuous respiratory monitoring. IEEE J. Biomed. Health Inform. 2020, 25, 1247–1256. [Google Scholar] [CrossRef] [PubMed]
  6. Garg, V.; Brewer, J. Telemedicine security: A systematic review. J. Diabetes Sci. Technol. 2011, 5, 768–777. [Google Scholar] [CrossRef] [PubMed]
  7. Abou El Kalam, A.; Ferreira, A.; Kratz, F. Bilateral teleoperation system using QoS and secure communication networks for telemedicine applications. IEEE Syst. J. 2015, 10, 709–720. [Google Scholar] [CrossRef]
  8. Rezaeibagha, F.; Mu, Y. Practical and secure telemedicine systems for user mobility. J. Biomed. Inform. 2018, 78, 24–32. [Google Scholar] [CrossRef]
  9. Murillo-Escobar, M.A.; Meranza-Castillón, M.O.; López-Gutiérrez, R.M.; Cruz-Hernández, C. Suggested integral analysis for chaos-based image cryptosystems. Entropy 2019, 21, 815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. 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]
  11. Mirrezapour, S.Z.; Zare, A. A new fractional sliding mode controller based on nonlinear fractional-order proportional integral derivative controller structure to synchronize fractional-order chaotic systems with uncertainty and disturbances. J. Vib. Control. 2021. [Google Scholar] [CrossRef]
  12. Kekha Javan, A.A.; Shoeibi, A.; Zare, A.; Hosseini Izadi, N.; Jafari, M.; Alizadehsani, R.; Moridian, P.; Mosavi, A.; Acharya, U.R.; Nahavandi, S. Design of Adaptive-Robust Controller for Multi-State Synchronization of Chaotic Systems with Unknown and Time-Varying Delays and Its Application in Secure Communication. Sensors 2021, 21, 254. [Google Scholar] [CrossRef] [PubMed]
  13. Zare, A.; Mirrezapour, S.Z.; Hallaji, M.; Shoeibi, A.; Jafari, M.; Ghassemi, N.; Alizadehsani, R.; Mosavi, A. Robust Adaptive Synchronization of a Class of Uncertain Chaotic Systems with Unknown Time-Delay. Appl. Sci. 2020, 10, 8875. [Google Scholar] [CrossRef]
  14. Pecora, L.M.; Carroll, T.L. Synchronization in chaotic systems. Phys. Rev. Lett. 1990, 64, 821. [Google Scholar] [CrossRef] [PubMed]
  15. Sweetha, S.; Sakthivel, R.; Harshavarthini, S. Finite-time synchronization of nonlinear fractional chaotic systems with stochastic actuator faults. Chaos Solitons Fractals 2021, 142, 110312. [Google Scholar] [CrossRef]
  16. Yang, F.; Mou, J.; Ma, C.; Cao, Y. Dynamic analysis of an improper fractional-order laser chaotic system and its image encryption application. Opt. Lasers Eng. 2020, 129, 106031. [Google Scholar] [CrossRef]
  17. Lin, C.F.; Shih, S.H.; Zhu, J.D. Chaos based encryption system for encrypting electroencephalogram signals. J. Med. Syst. 2014, 38, 1–10. [Google Scholar] [CrossRef] [PubMed]
  18. Murillo-Escobar, M.A.; Cardoza-Avendaño, L.; López-Gutiérrez, R.M.; Cruz-Hernández, C. A double chaotic layer encryption algorithm for clinical signals in telemedicine. J. Med. Syst. 2017, 41, 59. [Google Scholar] [CrossRef] [PubMed]
  19. Shahzadi, R.; Anwar, S.M.; Qamar, F.; Ali, M.; Rodrigues, J.J.; Alnowami, M. Secure EEG signal transmission for remote health monitoring using optical chaos. IEEE Access 2019, 7, 57769–57778. [Google Scholar] [CrossRef]
  20. Chen, J.; Ling, B.W.K.; Feng, P.; Lei, R. Computer cryptography through performing chaotic modulation on intrinsic mode functions with non-dyadic number of encrypted signals. IET Signal. Process. 2018, 13, 7–13. [Google Scholar] [CrossRef]
  21. Ibrahim, S.; Alhumyani, H.; Masud, M.; Alshamrani, S.S.; Cheikhrouhou, O.; Muhammad, G.; Hossain, M.S.; Abbas, A.M. Framework for Efficient Medical Image Encryption Using Dynamic S-Boxes and Chaotic Maps. IEEE Access 2020, 8, 160433–160449. [Google Scholar] [CrossRef]
  22. Gafsi, M.; Abbassi, N.; Hajjaji, M.A.; Malek, J.; Mtibaa, A. Improved Chaos-Based Cryptosystem for Medical Image Encryption and Decryption. Sci. Program. 2020, 2020, 6612390. [Google Scholar]
  23. Kumar, S.; Panna, B.; Jha, R.K. Medical image encryption using fractional discrete cosine transform with chaotic function. Med. Biol. Eng. Comput. 2019, 57, 2517–2533. [Google Scholar] [CrossRef]
  24. Liu, J.; Tang, S.; Lian, J.; Ma, Y.; Zhang, X. A novel fourth order chaotic system and its algorithm for medical image encryption. Multidimens. Syst. Signal. Process. 2019, 30, 1637–1657. [Google Scholar] [CrossRef]
  25. Stalin, S.; Maheshwary, P.; Shukla, P.K.; Maheshwari, M.; Gour, B.; Khare, A. Fast and secure medical image encryption based on non linear 4D logistic map and DNA sequences (NL4DLM_DNA). J. Med. Syst. 2019, 43, 1–17. [Google Scholar] [CrossRef]
  26. Banday, S.A.; Mir, A.H.; Malik, S. Multilevel medical image encryption for secure communication. In Advances in Computational Techniques for Biomedical Image Analysis; Academic Press: Cambridge, MA, USA, 2020; pp. 233–252. [Google Scholar]
  27. Chai, X.; Zhang, J.; Gan, Z.; Zhang, Y. Medical image encryption algorithm based on Latin square and memristive chaotic system. Multimed. Tools Appl. 2019, 78, 35419–35453. [Google Scholar] [CrossRef]
  28. Choi, U.S.; Cho, S.J.; Kang, S.W. Color medical image encryption using 3D Chaotic Cat Map and NCA. In Proceedings of the 2019 10th IFIP International Conference on New Technologies, Mobility and Security (NTMS), Canary Islands, Spain, 24–26 June 2019; pp. 1–5. [Google Scholar]
  29. Dagadu, J.C.; Li, J.; Shah, F. An efficient di-chaotic diffusion based medical image cryptosystem. In Proceedings of the 2017 14th International Computer Conference on Wavelet Active Media Technology and Information Processing (ICCWAMTIP), Chengdu, China, 15–17 December 2017; pp. 206–210. [Google Scholar]
  30. Ismail, S.M.; Said, L.A.; Radwan, A.G.; Madian, A.H.; Abu-Elyazeed, M.F. Generalized double-humped logistic map-based medical image encryption. J. Adv. Res. 2018, 10, 85–98. [Google Scholar] [CrossRef]
  31. Li, X.; Wu, F.; Khan, M.K.; Xu, L.; Shen, J.; Jo, M. A secure chaotic map-based remote authentication scheme for telecare medicine information systems. Future Gener. Comput. Syst. 2018, 84, 149–159. [Google Scholar] [CrossRef]
  32. Pandey, A.; Singh, B.; Saini, B.S.; Sood, N. A novel fused coupled chaotic map based confidential data embedding-then-encryption of electrocardiogram signal. Biocybern. Biomed. Eng. 2019, 39, 282–300. [Google Scholar] [CrossRef]
  33. Sangavi, V.; Thangavel, P. An exotic multi-dimensional conceptualization for medical image encryption exerting Rossler system and Sine map. J. Inf. Secur. Appl. 2020, 55, 102626. [Google Scholar]
  34. Thakur, S.; Singh, A.K.; Ghrera, S.P.; Elhoseny, M. Multi-layer security of medical data through watermarking and chaotic encryption for tele-health applications. Multimed. Tools Appl. 2019, 78, 3457–3470. [Google Scholar] [CrossRef]
  35. Zheng, L.; Wang, Z.; Tian, S. Comparative study on electrocardiogram encryption using elliptic curves cryptography and data encryption standard for applications in Internet of medical things. Concurr. Comput. Pract. Exp. 2020, e5776. [Google Scholar] [CrossRef]
  36. Murillo-Escobar, M.A.; Cruz-Hernández, C.; Abundiz-Pérez, F.; López-Gutiérrez, R.M.; Del Campo, O.A. A RGB image encryption algorithm based on total plain image characteristics and chaos. Signal. Process. 2015, 109, 119–131. [Google Scholar] [CrossRef]
  37. Ye, G.; Jiao, K.; Wu, H.; Pan, C.; Huang, X. An asymmetric image encryption algorithm based on a fractional-order chaotic system and the RSA public-key cryptosystem. Int. J. Bifurc. Chaos 2020, 30, 2050233. [Google Scholar] [CrossRef]
  38. Ding, L.; Ding, Q. A Novel Image Encryption Scheme Based on 2D Fractional Chaotic Map, DWT and 4D Hyper-chaos. Electronics 2020, 9, 1280. [Google Scholar] [CrossRef]
  39. Zhang, S.; Gao, T.; Gao, L. A novel encryption frame for medical image with watermark based on hyperchaotic system. Math. Probl. Eng. 2014, 2014, 240749. [Google Scholar] [CrossRef]
  40. Ye, G.; Wong, K.W. An image encryption scheme based on time-delay and hyperchaotic system. Nonlinear Dyn. 2013, 71, 259–267. [Google Scholar] [CrossRef]
  41. Alizadehsani, R.; Roshanzamir, M.; Hussain, S.; Khosravi, A.; Koohestani, A.; Zangooei, M.H.; Abdar, M.; Beykikhoshk, A.; Shoeibi, A.; Zare, A.; et al. Handling of uncertainty in medical data using machine learning and probability theory techniques: A review of 30 years (1991–2020). Ann. Oper. Res. 2021, 1–42. [Google Scholar] [CrossRef]
  42. Ye, G.; Pan, C.; Dong, Y.; Shi, Y.; Huang, X. Image encryption and hiding algorithm based on compressive sensing and random numbers insertion. Signal. Process. 2020, 172, 107563. [Google Scholar] [CrossRef]
  43. Zhou, J.; Zhou, N.R.; Gong, L.H. Fast color image encryption scheme based on 3D orthogonal Latin squares and matching matrix. Opt. Laser Technol. 2020, 131, 106437. [Google Scholar] [CrossRef]
  44. Sharifrazi, D.; Alizadehsani, R.; Joloudari, J.H.; Shamshirband, S.; Hussain, S.; Sani, Z.A.; Hasanzadeh, F.; Shoaibi, A.; Dehzangi, A.; Alinejad-Rokny, H. CNN-KCL: Automatic Myocarditis Diagnosis using Convolutional Neural Network Combined with K-means Clustering. Preprints 2020. [Google Scholar] [CrossRef]
  45. Bashashati, A.; Fatourechi, M.; Ward, R.K.; Birch, G.E. A survey of signal processing algorithms in brain–computer interfaces based on electrical brain signals. J. Neural Eng. 2007, 4, R32. [Google Scholar] [CrossRef] [Green Version]
  46. Chow, L.S.; Paramesran, R. Review of medical image quality assessment. Biomed. Signal. Process. Control. 2016, 27, 145–154. [Google Scholar] [CrossRef]
  47. Jin, X.; Wang, Z.; Yang, H.; Song, Q.; Xiao, M. Synchronization of Multiplex Networks with Stochastic Perturbations via Pinning Adaptive Control. J. Frankl. Inst. 2021, 358, 3994–4012. [Google Scholar] [CrossRef]
  48. Ji, Y.; Gong, Y.; Su, S.; Bai, X. Fixed-Time Synchronization for Different Dimensional Complex Network Systems with Unknown Parameters via Adaptive Control. Complexity 2021, 2021, 6680287. [Google Scholar] [CrossRef]
  49. Mofid, O.; Momeni, M.; Mobayen, S.; Fekih, A. A Disturbance-Observer-Based Sliding Mode Control for the Robust Synchronization of Uncertain Delayed Chaotic Systems: Application to Data Security. IEEE Access 2021, 9, 16546–16555. [Google Scholar] [CrossRef]
  50. Yao, Q. Synchronization of second-order chaotic systems with uncertainties and disturbances using fixed-time adaptive sliding mode control. Chaos Solitons Fractals 2021, 142, 110372. [Google Scholar] [CrossRef]
  51. Balootaki, M.A.; Rahmani, H.; Moeinkhah, H.; Mohammadzadeh, A. Non-singleton fuzzy control for multi-synchronization of chaotic systems. Appl. Soft Comput. 2021, 99, 106924. [Google Scholar] [CrossRef]
  52. Soleimanizadeh, A.; Nekoui, M.A. Optimal type-2 fuzzy synchronization of two different fractional-order chaotic systems with variable orders with an application to secure communication. Soft Comput. 2021, 25, 6415–6426. [Google Scholar] [CrossRef]
  53. Rigatos, G.; Abbaszadeh, M. Nonlinear optimal control and synchronization for chaotic electronic circuits. J. Comput. Electron. 2021, 20, 1050–1063. [Google Scholar] [CrossRef]
  54. Motallebzadeh, F.; Motlagh MR, J.; Cherati, Z.R. Synchronization of different-order chaotic systems: Adaptive active vs. optimal control. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 3643–3657. [Google Scholar] [CrossRef]
  55. Khan, A.; Kumar, S. TS fuzzy modeling and predictive control and synchronization of chaotic satellite systems. Int. J. Model. Simul. 2019, 39, 203–213. [Google Scholar] [CrossRef]
  56. Kaheni, H.R.; Yaghoobi, M. A new approach in anti-synchronization of a fractional-order hyper-chaotic Duffing system based on new nonlinear predictive control. Int. J. Dyn. Control. 2020, 8, 917–931. [Google Scholar] [CrossRef]
  57. Gupta, S.; Varshney, P.; Srivastava, S. Whale optimization based synchronization and control of two identical fractional order financial chaotic systems. J. Intell. Fuzzy Syst. 2021, 1–14. [Google Scholar] [CrossRef]
  58. Xiong, P.Y.; Jahanshahi, H.; Alcaraz, R.; Chu, Y.M.; Gómez-Aguilar, J.F.; Alsaadi, F.E. Spectral Entropy Analysis and Synchronization of a Multi-Stable Fractional-Order Chaotic System using a Novel Neural Network-Based Chattering-Free Sliding Mode Technique. Chaos Solitons Fractals 2021, 144, 110576. [Google Scholar] [CrossRef]
  59. Lin, C.H.; Hu, G.H.; Chan, C.Y.; Yan, J.J. Chaos-Based Synchronized Dynamic Keys and Their Application to Image Encryption with an Improved AES Algorithm. Appl. Sci. 2021, 11, 1329. [Google Scholar] [CrossRef]
  60. Moon, S.; Baik, J.J.; Seo, J.M. Chaos synchronization in generalized Lorenz systems and an application to image encryption. Commun. Nonlinear Sci. Numer. Simul. 2021, 96, 105708. [Google Scholar] [CrossRef]
  61. Wan, P.Y.; Liao, T.L.; Yan, J.J.; Tsai, H.H. Discrete sliding mode control for chaos synchronization and its application to an improved El-Gamal cryptosystem. Symmetry 2019, 11, 843. [Google Scholar] [CrossRef] [Green Version]
  62. Yau, H.T.; Hung, T.H.; Hsieh, C.C. Bluetooth based chaos synchronization using particle swarm optimization and its applications to image encryption. Sensors 2012, 12, 7468–7484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Chen, X.; Park, J.H.; Cao, J.; Qiu, J. Sliding mode synchronization of multiple chaotic systems with uncertainties and disturbances. Appl. Math. Comput. 2017, 308, 161–173. [Google Scholar] [CrossRef]
  64. Ren, H.P.; Bai, C.; Huang, Z.Z.; Grebogi, C. Secure communication based on hyperchaotic Chen system with time-delay. Int. J. Bifurc. Chaos 2017, 27, 1750076. [Google Scholar] [CrossRef] [Green Version]
  65. Lin, T.C.; Huang, F.Y.; Du, Z.; Lin, Y.C. Synchronization of fuzzy modeling chaotic time delay memristor-based Chua’s circuits with application to secure communication. Int. J. Fuzzy Syst. 2015, 17, 206–214. [Google Scholar] [CrossRef]
  66. Kwon, O.M.; Park, J.H.; Lee, S.M. Secure communication based on chaotic synchronization via interval time-varying delay feedback control. Nonlinear Dyn. 2011, 63, 239–252. [Google Scholar] [CrossRef]
  67. Shoeibi, A.; Ghassemi, N.; Alizadehsani, R.; Rouhani, M.; Hosseini-Nejad, H.; Khosravi, A.; Panahiazar, M.; Nahavandi, S. A comprehensive comparison of handcrafted features and convolutional autoencoders for epileptic seizures detection in EEG signals. Expert Syst. Appl. 2021, 163, 113788. [Google Scholar] [CrossRef]
  68. Ghassemi, N.; Shoeibi, A.; Rouhani, M.; Hosseini-Nejad, H. Epileptic seizures detection in EEG signals using TQWT and ensemble learning. In Proceedings of the 2019 9th International Conference on Computer and Knowledge Engineering (ICCKE), Mashhad, Iran, 24–25 October 2019; pp. 403–408. [Google Scholar]
  69. Ahmadi, A.; Kashefi, M.; Shahrokhi, H.; Nazari, M.A. Computer aided diagnosis system using deep convolutional neural networks for ADHD subtypes. Biomed. Signal. Process. Control. 2021, 63, 102227. [Google Scholar] [CrossRef]
  70. Chen, H.; Song, Y.; Li, X. A deep learning framework for identifying children with ADHD using an EEG-based brain network. Neurocomputing 2019, 356, 83–96. [Google Scholar] [CrossRef]
  71. Sadeghi, D.; Shoeibi, A.; Ghassemi, N.; Moridian, P.; Khadem, A.; Alizadehsani, R.; Teshnehlab, M.; Gorriz, J.M.; Nahavandi, S. An Overview on Artificial Intelligence Techniques for Diagnosis of Schizophrenia Based on Magnetic Resonance Imaging Modalities: Methods, Challenges, and Future Works. arXiv 2021, arXiv:2103.03081. [Google Scholar]
  72. Górriz, J.M.; Ramírez, J.; Ortíz, A.; Martínez-Murcia, F.J.; Segovia, F.; Suckling, J.; Leming, M.; Zhang, Y.-D.; Álvarez-Sánchez, J.R.; Bologna, G.; et al. Artificial intelligence within the interplay between natural and artificial computation: Advances in data science, trends and applications. Neurocomputing 2020, 410, 237–270. [Google Scholar] [CrossRef]
  73. Shoeibi, A.; Khodatars, M.; Jafari, M.; Moridian, P.; Rezaei, M.; Alizadehsani, R.; Khozeimeh, F.; Gorriz, J.M.; Heras, J.; Panahiazar, M.; et al. Applications of Deep Learning Techniques for Automated Multiple Sclerosis Detection Using Magnetic Resonance Imaging: A Review. arXiv 2021, arXiv:2105.04881. [Google Scholar]
  74. Ghassemi, N.; Shoeibi, A.; Rouhani, M. Deep neural network with generative adversarial networks pre-training for brain tumor classification based on MR images. Biomed. Signal. Process. Control. 2020, 57, 101678. [Google Scholar] [CrossRef]
  75. Manni, F.; van der Sommen, F.; Fabelo, H.; Zinger, S.; Shan, C.; Edström, E.; Elmi-Terander, A.; Ortega, S.; Callicó, G.M.; de With, P.H.N. Hyperspectral imaging for glioblastoma surgery: Improving tumor identification using a deep spectral-spatial approach. Sensors 2020, 20, 6955. [Google Scholar] [CrossRef]
  76. Ahsan, M.; Based, M.; Haider, J.; Kowalski, M. COVID-19 Detection from Chest X-ray Images Using Feature Fusion and Deep Learning. Sensors 2021, 21, 1480. [Google Scholar]
  77. Mohammadpoor, M.; Shoeibi, A.; Shojaee, H. A hierarchical classification method for breast tumor detection. Iran. J. Med. Phys. 2016, 13, 261–268. [Google Scholar]
  78. Mambou, S.J.; Maresova, P.; Krejcar, O.; Selamat, A.; Kuca, K. Breast cancer detection using infrared thermal imaging and a deep learning model. Sensors 2018, 18, 2799. [Google Scholar] [CrossRef] [Green Version]
  79. Shoeibi, A.; Khodatars, M.; Alizadehsani, R.; Ghassemi, N.; Jafari, M.; Moridian, P.; Khadem, A.; Sadeghi, D.; Hussain, S.; Zare, A.; et al. Automated detection and forecasting of covid-19 using deep learning techniques: A review. arXiv 2020, arXiv:2007.10785. [Google Scholar]
  80. Alizadehsani, R.; Sharifrazi, D.; Izadi, N.H.; Joloudari, J.H.; Shoeibi, A.; Gorriz, J.M.; Hussain, S.; Arco, J.E.; Sani, Z.A.; Khozeimeh, F.; et al. Uncertainty-aware semi-supervised method using large unlabelled and limited labeled COVID-19 data. arXiv 2021, arXiv:2102.06388. [Google Scholar]
  81. Ayoobi, N.; Sharifrazi, D.; Alizadehsani, R.; Shoeibi, A.; Gorriz, J.M.; Moosaei, H.; Khosravi, A.; Nahavandi, S.; Chofreh, A.G.; Goni, F.A.; et al. Time Series Forecasting of New Cases and New Deaths Rate for COVID-19 using Deep Learning Methods. arXiv 2021, arXiv:2104.15007. [Google Scholar]
  82. Ghassemi, N.; Shoeibi, A.; Khodatars, M.; Heras, J.; Rahimi, A.; Zare, A.; Pachori, R.B.; Gorriz, J.M. Automatic Diagnosis of COVID-19 from CT Images using CycleGAN and Transfer Learning. arXiv 2021, arXiv:2104.11949. [Google Scholar]
  83. Khozeimeh, F.; Sharifrazi, D.; Izadi, N.H.; Joloudari, J.H.; Shoeibi, A.; Alizadehsani, R.; Gorriz, J.M.; Hussain, S.; Sani, Z.A.; Moosaei, H.; et al. CNN AE: Convolution Neural Network combined with Autoencoder approach to detect survival chance of COVID 19 patients. arXiv 2021, arXiv:2104.08954. [Google Scholar]
  84. Alizadehsani, R.; Khosravi, A.; Roshanzamir, M.; Abdar, M.; Sarrafzadegan, N.; Shafie, D.; Khozeimeh, F.; Shoeibi, A.; Nahavandi, S.; Panahiazar, M.; et al. Coronary Artery Disease Detection Using Artificial Intelligence Techniques: A Survey of Trends, Geographical Differences and Diagnostic Features 1991–2020. Comput. Biol. Med. 2020, 128, 104095. [Google Scholar] [CrossRef] [PubMed]
  85. Fereydouneyan, F.; Zare, A.; Mehrshad, N. Using a fuzzy controller optimized by a genetic algorithm to regulate blood glucose level in type 1 diabetes. J. Med. Eng. Technol. 2011, 35, 224–230. [Google Scholar] [CrossRef] [PubMed]
  86. Namadchian, Z.; Zare, A. Stability analysis of dynamic nonlinear interval type-2 TSK fuzzy control systems based on describing function. Soft Comput. 2020, 24, 14623–14636. [Google Scholar] [CrossRef]
  87. Baydokhty, M.E.; Zare, A.; Balochian, S. Performance of optimal hierarchical type 2 fuzzy controller for load–frequency system with production rate limitation and governor dead band. Alex. Eng. J. 2016, 55, 379–397. [Google Scholar] [CrossRef] [Green Version]
  88. Bajestani, N.S.; Kamyad, A.V.; Zare, A. A piecewise type-2 fuzzy regression model. Int. J. Comput. Intell. Syst. 2017, 10, 734–744. [Google Scholar] [CrossRef] [Green Version]
  89. Zare, A.; Okauti, M. Automatic road extraction based on neuro-fuzzy algorithm. In Proceedings of the 10th WSEAS International Conference on Robotics, Control and Manufacturing Technology, Hangzhou, China, 11–13 April 2010. [Google Scholar]
  90. Bajestani, N.S.; Kamyad, A.V.; Esfahani, E.N.; Zare, A. Nephropathy forecasting in diabetic patients using a GA-based type-2 fuzzy regression model. Biocybern. Biomed. Eng. 2017, 37, 281–289. [Google Scholar] [CrossRef]
Figure 1. Structure of transmission synchronization.
Figure 1. Structure of transmission synchronization.
Sensors 21 03925 g001
Figure 2. Structure of circular synchronization.
Figure 2. Structure of circular synchronization.
Sensors 21 03925 g002
Figure 3. Block diagram of chaotic masking with the aid of multi-state synchronization.
Figure 3. Block diagram of chaotic masking with the aid of multi-state synchronization.
Sensors 21 03925 g003
Figure 4. Multi-state synchronization error curves (case A).
Figure 4. Multi-state synchronization error curves (case A).
Sensors 21 03925 g004
Figure 5. Multi-mode synchronization control effort curves (case A).
Figure 5. Multi-mode synchronization control effort curves (case A).
Sensors 21 03925 g005
Figure 6. Synchronization error curves for disturbance and structural uncertainty (case B).
Figure 6. Synchronization error curves for disturbance and structural uncertainty (case B).
Sensors 21 03925 g006
Figure 7. Estimation error curves of bounds related to uncertainties (right) and disturbances (left)—(case B).
Figure 7. Estimation error curves of bounds related to uncertainties (right) and disturbances (left)—(case B).
Sensors 21 03925 g007
Figure 8. Curves of synchronization errors in case of disturbance and structural uncertainty of case (C).
Figure 8. Curves of synchronization errors in case of disturbance and structural uncertainty of case (C).
Sensors 21 03925 g008
Figure 9. Estimation error curves for uncertainty (right) and disturbance (left) case (C).
Figure 9. Estimation error curves for uncertainty (right) and disturbance (left) case (C).
Sensors 21 03925 g009
Figure 10. Phase curves for master chaotic system (green), Slave 1 (blue) and slave 2 (red) chaotic systems.
Figure 10. Phase curves for master chaotic system (green), Slave 1 (blue) and slave 2 (red) chaotic systems.
Sensors 21 03925 g010
Figure 11. Standard benchmark images used in the proposed method.
Figure 11. Standard benchmark images used in the proposed method.
Sensors 21 03925 g011
Figure 12. X-ray images used in the proposed method.
Figure 12. X-ray images used in the proposed method.
Sensors 21 03925 g012
Figure 13. CT images used in the proposed method.
Figure 13. CT images used in the proposed method.
Sensors 21 03925 g013
Figure 14. Standard benchmarks images encryption using the proposed method.
Figure 14. Standard benchmarks images encryption using the proposed method.
Sensors 21 03925 g014
Figure 15. X-ray images encryption using the proposed method.
Figure 15. X-ray images encryption using the proposed method.
Sensors 21 03925 g015
Figure 16. CT images encrypted using the proposed method.
Figure 16. CT images encrypted using the proposed method.
Sensors 21 03925 g016
Figure 17. Displayed histograms for various images encrypted using the proposed method.
Figure 17. Displayed histograms for various images encrypted using the proposed method.
Sensors 21 03925 g017
Table 1. Statistical metrics for CT modality without noise.
Table 1. Statistical metrics for CT modality without noise.
ImagesHistogramCorrelationDifferential AttackPSNRInformation Entropy
StandardEncryptedNPCR (%)UACI (%)
Image 11,555,1643,375,5080.997099.61133.461 34.6235.5081
Image 21,942,7775,880,5770.996499.61033.460 33.1865.2665
Image 32,926,3666,511,8500.997699.60933.463 34.9824.8647
Image 41,969,8295,299,0180.995499.61133.459 33.1825.077
Image 52,776,0875,805,2490.996799.61033.460 34.9454.6828
Image 61,983,2215,394,1680.995699.61033.462 33.1985.0558
Image 73,108,7046,811,4270.997499.61133.460 35.0534.5683
Image 82,783,6718,248,9270.995899.61033.460 33.4834.5248
Image 92,887,1186,433,4370.997699.61233.461 34.9794.7115
Image 101,168,3853,379,1350.995699.61033.463 32.8505.6934
Table 2. Statistical metrics for CT Modality (Gaussian noise with 0.001 variance).
Table 2. Statistical metrics for CT Modality (Gaussian noise with 0.001 variance).
ImagesHistogramCorrelationDifferential AttackPSNRInformation Entropy
StandardEncryptedNPCR (%)UACI (%)
Image 1371,790.95414,854.680.996599.61133.45934.1466.4966
Image 2637,672.20775,806.220.995899.60933.46132.5995.9531
Image 3704,174.55770,799.320.996999.61033.46234.2845.8406
Image 4655,672.19823,316.520.994199.61233.45932.5835.7021
Image 5711,500.39800,070.880.995999.60833.460 34.2865.5608
Image 6658,916.62837,271.820.994499.61133.461 32.6035.6814
Image 7764,189.21849,713.850.996699.61233.463 34.3045.5387
Image 8910,672.691,120,143.860.994299.61133.462 32.7285.3159
Image 9712,326.25774,208.570.996999.61333.458 34.2765.6761
Image 10377,977.28471,233.220.994899.61133.461 32.4266.4675
Table 3. Statistical metrics for CT Modality (Gaussian noise with 0.003 variance).
Table 3. Statistical metrics for CT Modality (Gaussian noise with 0.003 variance).
ImagesHistogramCorrelationDifferential AttackPSNRInformation Entropy
StandardEncryptedNPCR (%)UACI (%)
Image 1239,642.29273,559.380.996399.61133.461 34.0526.6729
Image 2399,917.39493,824.290.995499.60933.462 32.3766.2428
Image 3439,060.82492,967.110.996799.61333.458 34.1416.1468
Image 4453,864.31566,619.140.993799.61033.460 32.4165.9343
Image 5490,271.52542,562.670.995699.61133.461 34.1155.8208
Image 6450,636.56557,774.790.993999.61233.459 32.3995.9312
Image 7505,486.28558,333.490.996399.61333.463 34.1565.8280
Image 8585,431.25723,173.690.993699.60933.462 32.4815.6384
Image 9456,812.35501,600.980.996699.61133.462 34.1205.9823
Image 10244,061.98307,606.740.994499.61233.460 32.2636.6424
Table 4. Statistical metrics for CT Modality (Gaussian noise with 0.006 variance).
Table 4. Statistical metrics for CT Modality (Gaussian noise with 0.006 variance).
ImagesHistogramCorrelationDifferential AttackPSNRInformation Entropy
StandardEncryptedNPCR (%)UACI (%)
Image 1183,606.33207,963.360.996199.61133.461 33.9756.7881
Image 2287,175.67347,957.530.995199.61033.459 32.2716.4632
Image 3318,969.51351,484.970.996599.60833.458 34.0276.3680
Image 4353,759.37420,480.950.993499.61133.462 32.3026.1234
Image 5379,975.17413,399.530.995499.60933.457 34.0546.0266
Image 6347,921.57414,116.530.993699.61133.456 32.3046.1326
Image 7381,476.11422,535.340.996299.61233.461 34.0976.0471
Image 8435,904.80517,108.820.993299.61033.460 32.3545.8833
Image 9333,557.81368,586.440.996599.61133.459 34.0736.2279
Image 10185,554.28227,767.280.994299.60733.454 32.1856.7663
Table 5. Statistical metrics for CT Modality (Gaussian noise with 0.009 variance).
Table 5. Statistical metrics for CT Modality (Gaussian noise with 0.009 variance).
ImagesHistogramCorrelationDifferential AttackPSNRInformation Entropy
StandardEncryptedNPCR (%)UACI (%)
Image 1160,406.07181,599.020.996199.60933.459 33.9656.8530
Image 2234,952.50276,693.140.995099.61133.461 32.2376.6058
Image 3262,347.29294,388.460.996499.60833.458 34.0076.5026
Image 4300,806.96345,570.370.993199.61033.457 32.2046.2599
Image 5319,370.10347,968.940.995399.61233.462 34.0466.1673
Image 6295,654.72341,419.340.993499.60733.457 32.2136.2688
Image 7320,298.48349,954.050.996099.61133.459 33.9956.1951
Image 8364,283.38423,633.800.993099.60933.462 32.2906.0396
Image 9276400.42307161.050.996499.61033.460 34.0006.3806
Image 10162565.94191476.220.994099.61233.461 32.1396.8390
Table 6. Comparison of the proposed method with other related works.
Table 6. Comparison of the proposed method with other related works.
Properties Encryption MethodEncryption MethodData TypesWorks
DisturbanceUnknown ParameterUncertainty
Chaos Logic Map EEG Signals [17]
Double Chaotic Layer Encryption (DCLE)EEG, ECG Signals[18]
Optical Chaos (Additive Chaos Masking)EEG Signals[19]
Chaotic Modulation on the Intrinsic Mode FunctionsEEG, ECG Signals[20]
Dynamic S-Boxes and Chaotic Maps Medical Images[21]
Improvement Chaotic SystemMedical Images[22]
chaotic Map + Fractional Discrete Cosine Transform (FrDCT) CoefficientsMedical Images[23]
Fourth Order Chaotic SystemMedical Images[24]
Non Linear 4D Logistic Map and DNA Sequences (NL4DLM_DNA)Medical Images [25]
Chaotic Method Based on Arnold’s Cat MapMRI Images[26]
Latin Square + Memristive Chaotic SystemMedical Images[27]
3D Chaotic Cat Map + NCAMedical Images[28]
Multiple Chaotic Systems + MD5Medical Images[29]
Double-Humped Logistic MapMRI, X-ray Images[30]
Chaotic Map-based Remote Authentication SchemeMedical Informatics [31]
Fused Coupled Chaotic Map (FCCM)ECG Signals[32]
Rossler Dynamical Chaotic system + Sine MapMedical Images[33]
DWT, DCT, SVD + Chaotic SystemMRI Images[34]
Data Encryption Standard (DES), and Elliptic Curves Cryptography (ECC)EEG Signals[35]
Chen Hyper-Chaotic System + Adaptive-Robust Multi-Mode SynchronizationMedical Images (CT, X-ray), Standard Benchmarks Proposed Method
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Javan, A.A.K.; Jafari, M.; Shoeibi, A.; Zare, A.; Khodatars, M.; Ghassemi, N.; Alizadehsani, R.; Gorriz, J.M. Medical Images Encryption Based on Adaptive-Robust Multi-Mode Synchronization of Chen Hyper-Chaotic Systems. Sensors 2021, 21, 3925. https://0-doi-org.brum.beds.ac.uk/10.3390/s21113925

AMA Style

Javan AAK, Jafari M, Shoeibi A, Zare A, Khodatars M, Ghassemi N, Alizadehsani R, Gorriz JM. Medical Images Encryption Based on Adaptive-Robust Multi-Mode Synchronization of Chen Hyper-Chaotic Systems. Sensors. 2021; 21(11):3925. https://0-doi-org.brum.beds.ac.uk/10.3390/s21113925

Chicago/Turabian Style

Javan, Ali Akbar Kekha, Mahboobeh Jafari, Afshin Shoeibi, Assef Zare, Marjane Khodatars, Navid Ghassemi, Roohallah Alizadehsani, and Juan Manuel Gorriz. 2021. "Medical Images Encryption Based on Adaptive-Robust Multi-Mode Synchronization of Chen Hyper-Chaotic Systems" Sensors 21, no. 11: 3925. https://0-doi-org.brum.beds.ac.uk/10.3390/s21113925

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