Skip to main content
Have a personal or library account? Click to login
AegisSentinel: Secure Medical Image Transmission with HFQE-S Encryption and COA Optimized HAPCNN Attack Detection Cover

AegisSentinel: Secure Medical Image Transmission with HFQE-S Encryption and COA Optimized HAPCNN Attack Detection

Open Access
|Aug 2026

Full Article

1. Introduction

Medical images such as CT, MRI, and X-ray scans are often shared over networks for remote analysis and consultation. When in transit, the images are susceptible to errors due to network noise, compression, disruptions, and hacking. Hackers can tamper images in an unnoticeable way, which may result in false diagnoses and put patients at risk.

There are many systems that are currently under development to make sure that images can be sent in a secure way (e.g., via encryption). However, at the moment, it is not possible to check whether the transmitted images have been tampered with after being sent. Currently, the techniques used do not have the capability of distinguishing between all kinds of attacks, and there is an excessive reliance on hyperparameters for almost all CNN models without any criteria for their selection. The aforementioned three constraints motivate the development of AegisSentinel.

AegisSentinel addresses each gap specifically: (1) HFQE-S provides image-level encryption so that the content is protected during transmission, (2) HAPCNN classifies the received encrypted image into one of six categories (clean or one of five attack types) to verify integrity after transmission, and (3) COA automatically tunes the HAPCNN hyperparameters to maximize detection accuracy without manual effort. The main contributions are:

  • HFQE-S technique is adopted to achieve CT image security through the confusion and diffusion stages, having high entropy value almost equal to ideal one at 7.9973 bits in lossless decoding.

  • HAPCNN architecture implementation to detect and classify different attacks including Gaussian noise attack, salt and pepper attack, intensity shift attack, occlusion attack, and FGSM attack on encrypted medical images.

  • COA technique employment to tune hyperparameters of HAPCNN model with accuracy of 96.52% and F1 score of 96.49%, outperforming other techniques such as PSO and WOA.

2. Related Work

A novel construction technique of four-dimensional augmented Lü hyperchaos system was used by Lin et al. [1] to design a new algorithm called CIEA-4DALHS for the encryption of color images using bidirectional spiral cross-scrambling, bit plane substitution, and hierarchical regional diffusion along with cascading through cross-channels. The method has an extensive key space, high entropy of the ciphertext, and strong ability to resist all kinds of attacks like statistical and differential attacks. Bidirectional spiral cross-scrambling destroys the spatial relationship of neighboring pixels and makes the local change of the ciphertext distribute evenly across the whole encrypted image by the help of hierarchical regional diffusion method. But the method has been designed for the encryption of color images only and does not have any mechanism for detection of any sort of tampering or transmission attacks.

The attack resistant Unet based watermarking framework was designed by Zou et al. [2], where the authors have incorporated an adaptive weighting mechanism to automatically determine the fusion weights for skip connection, along with residual interaction spatial feature transformation mechanism to maintain the fine detail information in the image during watermark embedding. The above framework has demonstrated very good imperceptibility when subjected to JPEG compression along with various standard attacks while maintaining bit accuracy. The application domain of this research is copyright protection via digital watermarking, but not encryption of medical images.

In [3], Zhou et al. have suggested quantum image encryption and watermarking techniques utilizing Quantum baker map and Quantum discrete cosine transformation. In their algorithm, 3D Henon Hyperchaotic system has been used for diffusion of the pixel value to authenticate the identity and encrypt the image. The addition of the quantum technique is one more step towards enhancing the security level compared to the traditional technique. This algorithm provides very high invisibility and robustness for the watermark verification but is dependent on the quantum computer primitives that are presently not available in real time for medical data transmission.

In Zou et al. [4], the HRFMS watermarking model was presented with attention-mask layer that extracts features through multiple receptive fields to effectively assist in guiding the process of watermark embedding. This method exhibited high PSNR and low bit-error rates. Multi-scale feature extraction enables this model to embed watermarking information in various image granularities in an adaptive manner. This contributes to the robustness of the model against geo-metric and photometric transformations. Unfortunately, this study is constrained only to watermarking robustness and invisibility as it does not perform any form of encryption or classification of transmission attacks.

A novel chaotic 4D memristor system based on the Sprott-C chaotic system has been designed by Zhang et al. [5] and analyzed in respect of dynamical complexity using Lyapunov exponent, bifurcation analysis, and hardware realization using PCB board. The 4D memristor system has been used for designing an image encryption scheme for color images using CNN as key generation scheme. The role of the memristor is to introduce flux-controlled non-linearity into the system making the proposed system more complex than other integer-order chaotic systems. The 4D memristor chaotic system is having higher degree of randomness than other Sprott chaotic systems; however, work on attack detection/classification has not yet been done.

The color image encryption and authentication scheme by Hu et al. [6] was developed through a combination of the ghost imaging technique, multi-parameter quaternion discrete fractional angular transform using two chaotic systems and a dual round fusion-permutation process with a high entropy of (7.9984 bits) and a huge key space. The use of dual chaotic architecture provides for a larger parameter space. It is resilient against brute-force attack as well as against chosen plaintext attacks. Though authentication has been integrated in the process of ghost imaging, yet the process requires higher computations owing to the ghost imaging process and the quaternion transform step. The ghost imaging system does not discriminate between various kinds of transmission attacks.

Gong et al. [7] have proposed an image encryption and authentication scheme where pixel adaptive diffusion algorithm is used along with cipherbook to provide security against differential attacks irrespective of use of key generated from plaintext using singular value decomposition ghost imaging to prove the origin of image. Pixel adaptive diffusion is an adaptive diffusion scheme which changes the extent of diffusion depending on the image statistics and, therefore, provides better uniformity in ciphertext distribution than static diffusion schemes. Nevertheless, as in case of [6], this approach fails to classify different types of attacks. In addition, authentication using ghost imaging is a costly process and may not be feasible in IoT-based healthcare applications.

Veerasekharreddy et al. [8] proposed hyperchaotic Fibonacci Polynomial CNN for efficient encryption and detection of attacks during the transmission of medical images in IoT-based healthcare network infrastructure. The combination of encryption using the chaos-based technique and the classification using a CNN was proposed, which forms the direct methodology used in this study. The introduction of the Fibonacci polynomials enhances the non-linearity of the chaotic map used and the classification of the anomalies through CNN takes place based on the feature extraction from the encrypted ciphertext domain. However, despite its effectiveness, the existing approach does not include the bio-inspired optimization phase for its hyperparameters and does not give any comparison of multiple optimization methods in this case, hence motivating the proposed COA optimization in this paper.

Reversible RGB medical image encryption scheme for IoMT using deep reinforcement learning was developed by Mahalakshmi and Nagarajan [9], in which Deep Q-Network makes encryption decisions based on the statistical properties derived from the intermediate encrypted images. The reinforcement learning approach provides the capability for the system to constantly improve the encryption policy based on the statistics of the image, thus providing the system with the advantage of adaptability, which is otherwise unavailable with rule-based static encryption techniques. This system increases the adaptability of the technique according to the changes in the properties of the image, but the system lacks the attack classification phase and relies upon the stability of the RL policy.

Subathra and Thanikaiselvan [10] introduced a framework for the encryption of medical images based on an innovative hyperchaotic system with five dimensions that was used in conjunction with a specialized U-Net network for segmenting important parts of the image along with dynamic DNA encryption and zigzag scrambling. The U-Net part of the algorithm detects diagnostically important features of the image, such as tumors and organ borders, ensuring that those are encrypted with highest priority. This improves encryption through the allocation of computing power to the most relevant areas of the image, but there is no post-transmission attack detection system in the proposed solution.

Based on the aforementioned analysis, the following are the current methods which include: purely chaotic/hyperchaotic encryption methods which are not aware of attacks [1,5]; watermarking-based techniques for copyright or integrity purposes other than encryption [24]; hybrid authentication-encryption techniques to verify the origin but not to classify the attack types [6,7]; and encryption together with the detection method lacking in hyperparameter optimization [810]. The gaps shown above lead us to conclude that there is no unified methodology which takes care of both strong encryption, intelligent multi-class attack detection, and hyperparameter optimization. The gap identified above leads to the development of the AegisSentinel framework which includes image-bound Fibonacci Q-Matrix encryption, multi-class attack detection, and hyperparameter optimization within a single computationally efficient pipeline. Table 1 shows the comparative summary of related work.

Table 1.

Comparative Summary of Related Work.

RefPaper/YearMethodology UsedEncryption/Security TechniqueAttack Detection ModelLimitation
[1]Lin et al., Chin. Phys. B, 2026CIEA-4DALHS color encryption4D augmented Lü hyperchaotic system, spiral scrambling, bit-plane substitutionNoneNo attack detection; not validated on medical images
[2]Zou et al., Signal Processing, 2026Attack-resilient Unet watermarkingAdaptive weighting + RISFT moduleNone (watermark robustness only)Watermarking only, no encryption/confidentiality, no attack classification
[3]Zhou et al., Int. J. Theor. Phys., 2024Quantum encryption + watermarkingQBM + QDCT + 3D Henon hyperchaosNoneQuantum primitives impractical for real-time clinical use; no attack classifier
[4]Zou et al., Neurocomputing, 2025HRFMS watermarking modelMulti-scale convolution, attention maskNoneWatermarking-only scope; no encryption or attack detection
[5]Zhang et al., Digital Signal Processing, 20264D memristor chaotic encryptionSprott-C-based memristor system, CNN key generationNoneEncryption-only; no attack-detection/classification stage
[6]Hu et al., Expert Systems with Applications, 2026CGI + QMPDFrAT encryption-authenticationDual chaotic systems, quaternion transformAuthentication only (no classification)Computationally heavy; no attack-type classification
[7]Gong et al., J. Modern Optics, 2026Pixel adaptive diffusion encryption-authenticationCipherbook-based diffusion, SVDGIAuthentication onlyNo attack-type classification; ghost-imaging overhead
[8]Veerasekharreddy et al., Neurocomputing, 2025 (Base Paper)Hyperchaotic Fibonacci polynomial CNNHyperchaotic encryptionCNN-based attack classifierNo bio-inspired hyperparameter optimization
[9]Mahalakshmi & Nagarajan, Frontiers in AI, 2026DRL-based reversible encryptionDQN-driven adaptive chaotic encryptionNoneNo attack classification; RL policy stability concerns
[10]Subathra & Thanikaiselvan, Scientific Reports, 20255D hyperchaotic + U-Net segmentation5D hyperchaos, dynamic DNA encoding, zig-zag scramblingNoneNo attack detection; segmentation adds processing overhead

3. Proposed System: AegisSentinel

3.1. System Overview

These three elements are integrated into AegisSentinel since they can each solve a specific problem which none of the other two can resolve. HFQE-S is employed for encrypting images, as it creates image-specific keys using SHA-256. Any alteration in the original image would cause a new key to be generated, which makes it impossible to decrypt the image. HAPCNN is utilized for detecting attacks since the polynomial convolution layers of the network are more sensitive to subtle pixel changes of encrypted (noisy) images compared to conventional CNN layers and therefore can detect attacks that cannot be detected visually. COA is applied to optimize the HAPCNN, since determining hyperparameters manually such as learning rate, dropout, and polynomial degree is unreliable, COA uses a simulation approach inspired by crayfish foraging.

The designed two-phase approach consists of two phases; namely, Training and Deployment. In the Training phase, the number of images from the TCGA-LUAD dataset used will be 475. These images will undergo preprocessing and encryption using the HFQE method. Afterward, we will conduct transmission attacks on these images with the help of five techniques, which give us a total number of 8,550 images. These images will then be optimized by the Crayfish Optimization Algorithm.

Figure 1 illustrates the sender-side HFQE-S encryption pipeline and the receiver-side HAPCNN-based attack detection and decryption workflow.

Figure 1.

Overall architecture of the proposed AegisSentinel system.

In deployment phase, the medical image acts as input at the sender’s end while HFQE-S module is used for encrypting this image and generation of encrypted image along with SHA-256 key. This encrypted image is sent over the network. At the receiver’s end, the prediction regarding attack on image is done with help of the pre-trained HAPCNN model. If there is no attack detected then the image is decrypted by using SHA-256 key and the real image is returned. Otherwise, the image is discarded.

3.2. HFQE-S Encryption with SHA Key

HFQE-S is the encryption module of AegisSentinel and turns a plain CT image into ciphertext through two stages: a deterministic, image-specific key derivation stage based on SHA-256, followed by a multi-round confusion–diffusion stage based on Fibonacci Q-Matrices. Expressing both stages mathematically, rather than only procedurally, makes explicit why the resulting cipher is image-bound and resistant to statistical attack.

The image is first flattened, row by row, into a one-dimensional pixel vector, and this vector is hashed to obtain a 256-bit digest (Using Equation (1)):

1
H = SHA-256(P)
Because SHA-256 is a one-way, avalanche-sensitive hash function, the digest depends on every pixel of the image [11]; changing even a single pixel value changes it unpredictably, which is what makes the key derived from it unique to that exact image.

The first sixteen bytes of the digest are then reinterpreted as an integer seed, which becomes the single source of randomness from which every structural parameter of the cipher is derived (Using Equation (2)):

2
S=Int16(H[0:16])
Deriving the seed directly from the hash, rather than sampling it independently, is what ties the parameters below to the image content instead of to an externally stored key. The seed is reduced, through modular arithmetic, into the block size used for diffusion (Using Equation (3)):
3
b=(Smod3)+2
and into the number of confusion–diffusion rounds applied to the image (Using Equation (4)):
4
R=S/28mod3+2

The block size takes a value of 2, 3, or 4, and the round count takes the same range independently. Both therefore vary from one image to the next, so no two images are ever encrypted with the same round structure.

Two independent generators are then seeded from the same value - one to produce the permutation tables used for confusion, and one to produce the keystream used later for diffusion - and the resulting key is assembled from every parameter derived so far (Using Equation (5)):

5
K={b,R, T1:R,RNG2}

Because the key is fully determined by the seed, and the seed is fully determined by the hash, it never has to be transmitted or stored separately - the receiver reconstructs it from the decrypted image's own hash once decryption succeeds. Encryption proceeds over the derived number of rounds, each combining a confusion step with a diffusion step.

Figure 2 depicts HFQE-S encryption engine showing the confusion and diffusion operations applied over R rounds with XOR keystream between rounds.

Figure 2.

Overall architecture of the proposed HFQE-S image encryption framework.

Algorithm 1. Encryption Key Generation from Image

Require: Plain CT image I ∈ R256×256 (grayscale)

Ensure: Encryption key K

1: Flatten image I into pixel vector

 P ← [p1, p2, …, p65536]

2: Compute hash: H ← SHA-256(P)

3: Extract seed: S ← Integer(H[0 : 16])

4: Compute block size: b ← (S mod 3) + 2

5: Compute number of rounds: R ← (S >> 8) mod 3 + 2

6: Initialize random number generators:

   RNG1 ← PCG64(S + 1)

   RNG2 ← PCG64(S + 2)

7: Generate permutation vectors:

 T1, T2, …, TR ← RNG1

8: Construct key:

 K ← {b, R, T1, T2, …, TR, RNG2}

9: return K

In the confusion step, that round's permutation table reorders the pixel vector so that spatially adjacent pixels are scattered across the image (Using Equation (6)):

6
δj[i]=P[Tj[i]]

Permutation alone provides no statistical protection - it only relocates existing pixel values - which is why it is always paired with the diffusion step that follows.

Diffusion is carried out using a Fibonacci Q-Matrix whose order matches the block size chosen for that image [12]. For a general order N, the entries follow a simple placement rule: every entry in the first row is one, every other row has a single one placed directly below-left of the diagonal, and every remaining entry is zero (Using Equation (7)):

7
QN(i,j)=1 if i=1 or j=i1, else 0

This is the direct generalization of the classical two-term Fibonacci recurrence to N terms: repeatedly multiplying by a matrix built this way advances an N-step Fibonacci-type sequence exactly the way repeated addition advances the ordinary sequence, which is what gives the diffusion step its mixing behaviour. For the block sizes used in HFQE-S, this rule produces the following three matrices (Using Equations (810)):

8
Q2=|1110|
9
Q3=|111100010|
10
Q4=|1111100001000010|

These three matrices are the only diffusion matrices HFQE-S ever uses, since the derived block size is always 2, 3, or 4; which one applies to a given image is decided entirely by that image's own hash value, not by any externally chosen setting.

Algorithm 2: Image Encryption using HFQE-S Scheme

Require: Plain image I;

Key K = {b, R, T1, …, TR, RNG2}

Ensure: Encrypted image E; auxiliary key H

1: Flatten image I into pixel vector

  P ← [p1, p2, …, p65536]

2:  for j = 1 to R do

3:   Confusion (Permutation):

4:    δj[i] ← P Tj[i], i ∈ {1, …, 65536}

5:    Reshape δj into matrix Mj ∈ R256×256

6:   Diffusion (Fibonacci Q-Matrix):

7:   for each b × b block B in Mj do

8:   B’ ← Qb·B mod 256

9:  end for

10:  XOR Diffusion:

11:  for each pixel Mj(i) do

12:   Mj(i) ← Mj(i) ⊕ RNG2.integers(0, 256)

13:  end for

14:  P ← Flatten(Mj)

15: end for

16: Reshape P into encrypted image E ∈ R256×256

17: return E, H

Diffusion itself is applied by multiplying each block of the permuted image by the matrix matching the chosen block size, under modulo-256 arithmetic so that every result remains a valid pixel intensity (Using Equation(11)):

11
B=Qb·B(mod256)

Because each of these matrices is invertible under modulo-256 arithmetic, this step can be undone exactly during decryption using the corresponding inverse matrix.

Finally, a keystream drawn from the second generator is combined with the diffused block using bitwise XOR, adding a layer of randomization that does not depend on the Fibonacci structure at all (Using Equation (12)):

12
Mj(i)=Mj(i)κi

Because this mask changes for every pixel and every round, two images that differ by only a single pixel produce ciphertexts that look completely unrelated once all rounds are complete.

Together, Equations (1)(12) show how a single hash of the input image determines every downstream choice made by HFQE-S - the key, the block size, the round count, the permutation, and the diffusion matrix - which is what gives the scheme its image-binding property without requiring any key to be stored or transmitted separately.

The immutability of the SHA-256 algorithm ensures that a particular image would always generate the same key. The slightest alteration of the picture results in an entirely new hash value.

Diffusion step is carried out using the use of the Fibonacci Q-Matrices in the sizes 2×2, 3×3 and 4×4. Encryption round involves two processes that occur alternately. Confusion process is where permutations take place. Permutations are done in a process involving multiplication of pixels by the Q-Matrix in modulo 256. Furthering the security of the process, application of XOR-based key-streams is done between rounds thus adding a new level of randomization. The combined effect of diffusion and confusion with XOR across multiple rounds ensures that even a single-pixel change in the image propagates unpredictably thus ensuring avalanche characteristics essential for robust security.

One of the main advantages of this method is the fact that the hash created with the help of the SHA-256 algorithm is created based on the initial picture. Therefore, any changes made to the picture will create a different hash and make decryption of the changed file impossible. Decryption is the reverse process of encryption.

3.3. Dataset Preparation and Attack Simulation

Pre-processing of data is performed before it is used for training and testing of the model. The process begins with extraction of images from the TCGA-LUAD dataset, which contains 475 sets of images, each set comprising DICOM images extracted from CT scans. Pre-processing of the images includes extraction of the middle slice, conversion of the image to PNG format, scaling of the image to 256×256, and graying of the image.

These images are encrypted using the HFQE-S encryption algorithm. After encryption, the images are classified as “clean” and assigned a label value of 0. Other forms of attacks considered are the Gaussian Noise Attack, Salt & Pepper Noise Attack, Intensity Shift Attack, Occlusion Attack, and FGSM Attack.

Because an adversary intercepting the channel only ever has access to ciphertext, all five attacks used to stress-test HAPCNN are simulated directly on the encrypted image Ic(x, y) rather than on the original plaintext. Pixel intensities are normalized to the range [0, 1] before a perturbation is applied and clipped back to a valid range afterward. The five attack models are defined below, with each attack first defining its underlying random or structural variable and then using that variable in the equation that produces the attacked image.

The Gaussian Noise Attack models additive channel or sensor noise picked up during transmission. A noise field η(x, y) is first drawn independently at every pixel from a zero-mean Gaussian distribution with standard deviation σ, resampled for each image as σ ∈ [0.05, 0.15] (Using Equation(13)):

13
η(x,y)~N(0,σ2)

This noise field is then added to the encrypted image and clipped to a valid range to obtain the attacked image (Using Equation(14)):

14
Ig(x,y)=clip[Ic(x,y)+η(x,y),0,1]

Table 2 depicts the clean and attacked image variations.The Salt & Pepper Noise Attack models impulsive bit-level corruption, in which a random subset of pixels is forced to the extreme intensity values. A binary corruption mask M(x, y) is first drawn from a Bernoulli distribution with density p ∈ [0.03, 0.08] (Using Equation(15)):

15
M(x,y)~Bernoulli(p)

Table 2.

Clean and Attacked Image Variations.

Medical Image0: Clean1: Gaussian2: Salt&Pepper3: Intensity4: Occlusion5: FGSM
graphic/j_ias-2026-0017_ingr_001.pnggraphic/j_ias-2026-0017_ingr_002.pnggraphic/j_ias-2026-0017_ingr_003.pnggraphic/j_ias-2026-0017_ingr_004.pnggraphic/j_ias-2026-0017_ingr_005.pnggraphic/j_ias-2026-0017_ingr_006.pnggraphic/j_ias-2026-0017_ingr_007.png

Each corrupted pixel is then assigned to 1 (salt) or 0 (pepper) with equal probability through V(x, y). Using this mask, the attacked image is composed as (Using Equation(16)):

16
Isp(x,y)=(1−M(x,y))·Ic(x,y)+M(x,y)·V(x,y)

The Intensity Shift Attack models a uniform gain or exposure drift applied to the whole image. A single scalar offset Δ is first drawn once per image from a uniform distribution over [-0.30, 0.30] (Using Equation(17)):

17
Δ~U(0.30,0.30)

This offset is then applied identically to every pixel of the encrypted image and clipped to a valid range (Using Equation(18)):

18
Iint (x,y)=clip[Ic(x,y)+Δ,0,1]

The Occlusion Attack models a physical obstruction or packet-loss blackout covering part of the image. An indicator mask R(x, y) is first defined, equal to 1 inside a randomly positioned rectangular region covering 10–25% of the image area and 0 elsewhere (Using Equation(19)):

19
R(x,y){0,1}

A randomly chosen fill value c ∈ {0, 1} is then substituted into the masked region to obtain the attacked image (Using Equation(20)):

20
Iocc(x,y)=(1-R(x,y))·Ic(x,y)+R(x,y)·c

The FGSM-style Perturbation Attack applies a bounded, sign-based perturbation modeled on the Fast Gradient Sign Method, to stress-test HAPCNN against adversarial-style inputs. A per-pixel sign field g(x, y) is first drawn from {−1, +1} (Using Equation(21)):

21
g(x,y){1,+1}

Because this perturbation is generated at the dataset-construction stage, before any classifier exists to differentiate against, g(x, y) is instantiated as a random sign field rather than the gradient of a trained model’s loss; this preserves FGSM’s characteristic bounded, high-frequency perturbation pattern while keeping the attack classifier-agnostic. Using a perturbation budget ε ∈ [0.02, 0.10], the attacked image is then obtained by (Using Equation(22)):

22
Ifgsm(x,y)=clip[Ic(x,y)+ε·sign(g(x,y)),0,1]

Equations (13)(22) generate the five attacked image classes (labels 1–5), which are used alongside the unperturbed clean class (label 0) for training and evaluation of the proposed model.

Algorithm 3: Dataset Generation Pipeline

Require: 475 DICOM CT scan series from TCGA-LUAD dataset

Ensure: Balanced dataset of 8,550 labeled images stored as dataset.npz

1: Initialize P ← ∅

2: for each DICOM series Si in dataset do

3:   Extract center slice Ii from Si

4:   Convert Ii to PNG format

5:   Resize Ii to 256 × 256

6:   Convert Ii to grayscale

7:   Ic ← HFQE(Ii)

8:   Append Ic to P with label 0

9: end for

10: Split P into Ptrain (333, 70%), Pval (71, 15%), Ptest (71, 15%)

11: // Process Training set

12: Xtrain ← ∅, ytrain ← ∅

13: for each clean image Ic in Ptrain do

14:   Append Ic to Xtrain; append label 0 to ytrain

15: I1 ← GaussianNoise(Ic); append I1 to Xtrain; append label 1 to ytrain

16: I2 ← SaltPepper(Ic); append I2 to Xtrain; append label 2 to ytrain

17: I3 ← IntensityShift(Ic); append I3 to Xtrain; append label 3 to ytrain

18: I4 ← Occlusion(Ic); append I4 to Xtrain; append label 4 to ytrain

19: I5 ← FGSM(Ic); append I5 to Xtrain; append label 5 to ytrain

20: end for

21: Xtrain_h ← HorizontalFlip(Xtrain); ytrain_h ← ytrain

22: Xtrain_v ← VerticalFlip(Xtrain); ytrain_v ← ytrain

23: Xtrain_total ← Xtrain ∪ Xtrain_h ∪ Xtrain_v

24: ytrain_total ← ytrain ∪ ytrain_h ∪ ytrain_v

25: Similarly process the Validation set (steps 12–24) to obtain (Xval_total, yval_total)

26: Similarly process the Test set (steps 12–24) to obtain (Xtest_total, ytest_total)

27: Save (Xtrain_total, ytrain_total), (Xval_total, yval_total), (Xtest_total, ytest_total) as dataset.npz

Figure 3 illustrates the DICOM CT scan series from the TCGA-LUAD dataset are processed by extracting the center slice, converting to PNG, resizing to 256×256, and converting to grayscale.

Figure 3.

Data acquisition and preprocessing pipeline.

3.4. HAPCNN: Hierarchical Auto-Associative Polynomial CNN

The attack detection stage of AegisSentinel is carried out by the Hierarchical Auto-Associative Polynomial CNN (HAPCNN). Its input is a 256×256 encrypted grayscale image, and its output is a six way classification of that same image into one of six categories: clean ciphertext or one of five attack types (Gaussian noise, salt-and-pepper noise, intensity shift, occlusion, or FGSM perturbation). The network architecture is organized hierarchically, consisting of three stacked polynomial convolution layers, an auto-associative bottleneck layer, and a fully connected output layer.

Instead of the conventional ReLU activation, each convolutional stage applies a polynomial expansion of the filtered response. This design choice is central to the model, because random-looking ciphertext pixels carry information not in their raw magnitude but in the higher-order statistical relationships between neighbouring pixels. A first-order (linear) filter is blind to these relationships, whereas a polynomial filter can retain quadratic and higher cross-terms that change measurably when noise, occlusion, or an adversarial perturbation is injected into the ciphertext. The response produced by the 1-th polynomial convolution layer is given by (Using Equation(23)):

23
i=1nwi·xin+b

where x is the feature map received from the previous layer, wi is the i-th learnable polynomial coefficient (kernel) of the layer, n is the polynomial order, and b is the layer bias. Raising the convolutional response to successive powers before summation allows the layer to learn feature interactions of order greater than one, which is what lets HAPCNN pick up the faint, structured distortions that an attack leaves behind in an otherwise noise-like encrypted image.

The bottleneck stage that follows the three polynomial convolution layers performs an auto-associative mapping: it compresses the polynomial feature maps to a lower-dimensional code and then reconstructs the original feature map from that code, in the same manner as an autoencoder. Because this reconstruction can only succeed if the compressed code preserves the fine-grained structure of the ciphertext, the bottleneck acts as an internal consistency check - a clean ciphertext reconstructs with low error, while an attacked ciphertext, whose local statistics have been disturbed, reconstructs with a measurably larger error. This reconstruction signal is combined with the ordinary classification signal, and the two are optimized jointly through a composite loss (Using Equation (24)):

24
L=Lcls+λ Lrec

where Lcls is the classification loss, Lrec is the reconstruction loss produced by the auto-association layer, and λ is a weighting coefficient that balances the two terms during training. The classification term is a standard categorical cross-entropy over the six output classes (Using Equation (25)):

25
c=16yclogy^c

where yc is the one-hot ground-truth label for class c and ŷc is the probability that HAPCNN assigns to class c. The reconstruction term is the mean squared error between the original encrypted image and the image reconstructed from the bottleneck code (Using Equation (26)):

26
1/Ni=1N(xix^i)2

where xi is the i-th pixel of the original encrypted image, x^i is the corresponding reconstructed pixel produced by the bottleneck decoder, and N is the total number of pixels. Minimizing Lrec alongside Lcls forces the network to keep enough ciphertext detail through the bottleneck that any attack-induced deviation remains detectable, rather than being smoothed away as the network learns to classify.

Training proceeds for a maximum of 80 epochs with a batch size of 16, using the hyperparameter set Θ* returned by the Crayfish Optimization Algorithm described in the next subsection. At inference time, an output of class 0 signals a clean ciphertext and triggers decryption, while any output in the range 1–5 signals a detected attack and raises an alarm instead of proceeding to decryption.

Algorithm 4: Training Procedure for HAPCNN Model

Require: Training set Dtrain (5985 images),

   Validation set Dval (1282 images),

   Optimal hyperparameters Θ*

Ensure: Trained model weights hapcnn_weights.pth

1: Initialize HAPCNN model:

  Polynomial Convolution layers

  Bottleneck layer

  Fully connected output layer

2: Set batch size = 16, max epochs = 80

3: for e = 1 to 80 do

4: for each mini-batch B ⊂ Dtrain do

5:   Forward pass: ŷ ← HAPCNN(B)

6:   Compute loss: L ← Lcls + Lrec

7:  Backpropagate and update model parameters

8:  end for

9:  Evaluate model on validation set Dval

10:  if validation accuracy improves then

11:   Save model weights

12:  end if

13:  if no improvement for 20 epochs then

14:   break

15:  end if

16: end for

17: Load best saved model

18: Evaluate on test set (1283 images)

3.5. Crayfish Optimization Algorithm (COA)

The Crayfish Optimization Algorithm (COA) is a nature-inspired meta-heuristic that models the feeding and competing behaviour of crayfish as they respond to changes in water temperature [13]. In AegisSentinel, COA is used as a pre-training search procedure to determine the optimal values of the hyperparameters that most strongly affect HAPCNN accuracy - namely the learning rate η, the regularization coefficient λ, the dropout rate d, and the polynomial/filter parameter p. Rather than hand-tuning these values, COA treats each candidate hyperparameter set as an individual “crayfish” and lets a population of such individuals compete and forage until the best-performing configuration emerges.

Each candidate Ci is scored by briefly training HAPCNN with that configuration and measuring how well it performs on the validation set. Because COA is framed as a minimization procedure, the validation accuracy is converted into a fitness value where a smaller number is better (Using Equation (27)):

27
f(Ci)=1Accval(Ci)

where Accval(Ci) is the validation accuracy obtained by HAPCNN after a short training run using the hyperparameter values encoded in candidate Ci. A candidate that yields higher validation accuracy therefore receives a lower (better) fitness value, and the candidate with the smallest f(Ci) across the whole population is retained as the current best solution C.

COA's distinguishing feature is that its search behaviour is not fixed but is instead driven by a temperature value T, resampled uniformly from the interval [20, 35] at every iteration. This temperature acts as a switch between two complementary search behaviours. When the temperature is high (T > 30), the crayfish behave competitively: each individual is pulled toward the current best solution, narrowing the search around the most promising region found so far, which is an exploitation move (Using Equation (28)):

28
Cit+1=Cit+r1·(C*Cit),T>30

where Cit is candidate i's position (hyperparameter values) at iteration t, C is the current best candidate, and r1 is a random number drawn from [0, 1] that controls the step size of the move toward C. When the temperature drops (T ≤ 30), the crayfish instead disperse to forage, which is an exploration move that widens the search away from the current best and helps the algorithm avoid settling into a local optimum (Using Equation(29)):

29
Cit+1=C*+r2·ε, T30

where r2 is a second random scalar in [0, 1] and ε is a small random perturbation vector applied around the current best solution C. Alternating between Equation (28) and Equation (29) according to the sampled temperature is what keeps COA from converging prematurely - a known weakness of gradient-free optimizers such as PSO and WOA - while still concentrating most of the search effort near the strongest candidates found so far.

In practice, the search space is defined first, and an initial population of five candidate crayfish is drawn at random. Each candidate is trained for five epochs to obtain an initial fitness score, and the best scoring candidate becomes C*. The population then undergoes ten rounds of temperature-driven updates alternating between Equations (28) and (29) as dictated by the freshly sampled temperature each round with every candidate retrained for five epochs and rescored after each update.

Algorithm 5: Hyperparameter Optimization using COA

Require: Training dataset D

Ensure: Optimal hyperparameters Θ* = {η, λ, d, p}

1: Define hyperparameter space: Θ = {η, λ, d, p}

2: for i = 1 to 5 do

3:  Ci ← Random(Θ)

4: end for

5: for each crayfish Ci do

6:  Train model for 5 epochs using Ci

7:  f(Ci) ← 1 - ValidationAccuracy

8: end for

9: C* ← arg mini f(Ci)

10: for t = 1 to 10 do

11:  Generate temperature T ∈ [20, 35]

12:  for each crayfish Ci do

13:   if T > 30 then

14:    Ci ← Ci + r1 · (C* - Ci)

15:   else

16:    Ci ← C* + r2 · ε

17:   end if

18:   Train model for 5 epochs using Ci

19:   f(Ci) ← 1 - ValidationAccuracy

20:  end for

21:  C* ← arg mini f(Ci)

22: end for

23: return Θ* ← C*

4. Experimental Results

4.1. Encryption Module Performance

The performance of the proposed HFQE-S algorithm is measured using the standard cryptographically defined performance criteria for the CT scan images of the TCGA-LUAD dataset, which includes entropy, PSNR, and encryption/decryption times [14].

Entropy measures the randomness of the pixel values of the encrypted image. An ideal encrypted image has pixel intensities uniformly distributed over all 256 gray levels, giving a theoretical maximum entropy of 8 bits per pixel. It is calculated as (Using Equation (30)):

30
H=i=0255p(i)log2p(i)

where p(i) is the probability of pixel intensity value i occurring in the image. The entropy of the image encrypted using the HFQE-S scheme is 7.9973 bits per pixel, which differs from the ideal value by only 0.0027 bits per pixel. This implies that the encrypted output image is statistically indistinguishable from random noise. This high entropy is obtained through the combination of multiple rounds of Fibonacci Q-Matrix diffusion and key-stream application via the XOR operator, which removes correlation among neighbouring pixels in all directions. A high entropy corresponds to strong resistance against statistical attacks.

HFQE-S scrambling and diffusion analysis was performed on images with different derived block sizes and round counts. For all cases tested, the entropy values remain closely clustered, ranging from 7.9969 to 7.9974 bits per pixel, which shows that the level of security offered by HFQE-S is consistent irrespective of the block size and round count derived from the SHA-256 key-generation parameters.

To quantify the visual dissimilarity between two images, the Peak Signal-to-Noise Ratio (PSNR) is used. It is defined as (Using Equation(31)):

31
PSNR=10×log10(MAX2MSE)

where MAX is the maximum pixel intensity value (255 for 8-bit images) and MSE is the Mean Squared Error between the two images being compared. MSE is defined as (Using Equation(32)):

32
MSE=1kl=1k(dl-(dl))2

where dl is the actual pixel value, d^l is the corresponding reconstructed/compared pixel value, and k is the total number of pixels. Between the original and the encrypted images, the PSNR values obtained range between 5.38 dB and 9.15 dB across the tested block-size and round configurations, which is sufficiently low to confirm that no visual information about the underlying CT image leaks into the ciphertext. The variation in PSNR across configurations is expected, since different block sizes cause different degrees of pixel alteration per round, even though the resulting entropy remains uniformly high in every case.

In the HFQE-S algorithm, the only role of SHA-256 is deterministic key generation from the image pixels [15]. SHA-256 does not take part in the encryption process itself; its task is to generate a 256-bit message digest of the original image pixel vector, which is then used to determine the block size b, the number of rounds R, and the seeds of the two random number generators. This design allows different encryption keys to be obtained for different images without relying on any additional key -management system. One SHA-256 hash computation takes very little time relative to the overall encryption process — execution of a single SHA-256 hash over a 256×256 pixel image takes less than 0.1 ms out of the total encryption time of 2.87 ms, remaining well below the real-time threshold of 5 ms. Using SHA-256 for key binding rather than for encryption itself eliminates computational bottlenecks while still ensuring plaintext sensitivity.

From the above results, the PSNR between the source image and the encrypted image is very low, confirming that the encrypted image bears no resemblance to the source image and that no diagnostic information can be extracted from it without the key. However, using the key generated through SHA-256, decryption recovers the original image exactly, so that the PSNR between the decrypted image and the source image becomes infinite (i.e., MSE = 0 in Eq. 31). This confirms perfect, lossless decryption, a property that is critical for medical-imaging applications, where any pixel-level distortion could obscure a diagnostic finding. Encryption performance across images is shown in Table 3. HFQE-S Encryption Module Performance Metrics is shown in Table 4.

Table 3.

Encryption Performance Across Images.

ImageBlock SizeRoundsPSNR (dB)Entropy (bits)
graphic/j_ias-2026-0017_ingr_008.png4×446.527.9973
graphic/j_ias-2026-0017_ingr_009.png2×249.157.9972
graphic/j_ias-2026-0017_ingr_010.png3×338.007.9971
graphic/j_ias-2026-0017_ingr_011.png3×336.687.9974
graphic/j_ias-2026-0017_ingr_012.png4×425.387.9969
Table 4.

HFQE-S Encryption Module Performance Metrics.

MetricAchieved ValueIdeal
Entropy (bits/pixel)7.99738.0000
PSNR (Original vs. Encrypted)Very Low (noisy)Low
PSNR (Decrypted vs. Original)Infinite (lossless)Infinite
Encryption Time2.87 ms< 5 ms
Decryption Time2.35 ms< 5 ms

4.2. Attack Detection Performance

The attack-detection capability of HAPCNN is evaluated on 1,293 test images spanning six classes (clean and five attack types), using accuracy, precision, recall and F1 score [16]. Accuracy is calculated as (Using Equation(33)):

33
 Accuracy = Number of correctly classified images  Total number of test images ×100%

Precision, recall and F1 score together characterize the trade-off between false alarms and missed detections and are calculated as (Using Equations (34,35, and 36)):

34
 Precision =TPTP+FP
35
 Recall =TPTP+FN
36
F1=2× Precision × Recall  Precision + Recall 

where TP, FP and FN denote the number of true positives, false positives and false negatives, respectively.

Four configurations were compared: HAPCNN without an optimizer, HAPCNN tuned with PSO, HAPCNN tuned with WOA, and the proposed HAPCNN tuned with COA.

HAPCNN without any optimizer achieves an accuracy of 83.94%, with a precision of 84.13%, a recall of 83.92%, and an F1 score of 83.80%. This establishes a strong baseline and demonstrates the effectiveness of the polynomial convolution layers and the auto-association mechanism in detecting attacks on encrypted images. The closeness of precision and recall at this stage indicates that the base architecture carries little inherent class bias, even before hyperparameter tuning.

When PSO is used to tune HAPCNN, accuracy improves to 89.94% and F1 score to 89.80%, confirming that systematic hyperparameter search yields substantial gains over manual or default configurations. A further improvement is observed with WOA, which raises accuracy to 92.51% and F1 score to 92.71%, reflecting the improved exploration–exploitation balance of WOA relative to PSO. The best performance is obtained with the proposed COA-tuned HAPCNN, which achieves 96.52% accuracy, 96.60% precision, 96.48% recall, and 96.49% F1 score. The high and closely matched precision and recall achieved by COA-tuned HAPCNN indicate that the model maintains both high sensitivity and high specificity simultaneously — a property that is essential in medical-image security, where either a missed attack or a false alarm can have serious downstream consequences. Comparative analysis of detection performance is shown in table 5.

Table 5.

Attack Detection Performance Comparison.

ConfigurationAcc. (%)Prec. (%)Rec. (%)F1 (%)
HAPCNN only (No optimizer)83.9484.1383.9283.80
HAPCNN + PSO89.9489.1389.9289.80
HAPCNN + WOA92.5192.5592.3892.71
HAPCNN + COA (Proposed)96.5296.6096.4896.49

In addition to detection accuracy, the practical feasibility of HAPCNN for IoT-healthcare deployment was evaluated in terms of model size and inference speed. HAPCNN contains 243,654 trainable parameters and requires approximately 1.276 GFLOPs per forward pass on a single 256×256 input image, resulting in an average inference time of 68.71 ms per image on CPU. This lightweight computational footprint confirms that HAPCNN is well suited for deployment on resource-constrained IoT-healthcare edge devices without requiring specialized hardware acceleration.

To visualize the class-wise behaviour of the proposed HAPCNN + COA detector after retraining on the corrected dataset, the confusion matrix over the 1,293 test images is shown in Figure 4. The corresponding per-class precision, recall and F1 score are macro-averaged across all six classes (Using Equation(37)),

37
 Macro F1=1Cc=1CF1c
where C is the number of classes (C = 6 in this study) and F1c is the F1 score obtained for class c. The macro-averaged precision, recall and F1 score obtained across the six classes are 0.9660, 0.9648 and 0.9652, respectively, consistent with the overall detection performance reported above. Table 6 presents the per-class performance metrics based on macro-averaged evaluation.

Figure 4.

Confusion matrix of the HAPCNN + COA model on the corrected dataset (1,293 test images, 6 classes).

Table 6.

Per-Class Metrics (Macro-Averaged).

ClassTPRowColPrecisionRecallF1
Clean2242242350.95321.00000.9760
Gaussian2082232250.92440.93270.9281
Salt & Pepper1771821771.00000.97250.9861
Intensity2252302251.00000.97830.9890
Occlusion2232232231.00001.00001.0000
FGSM1912112080.91830.90520.9117
Macro Average---0.96600.96480.9652

4.3. Comparative Analysis

The following analysis of the encryption quality is presented to justify the reasonability of the proposed framework against recent benchmarks. The HFQE-S encryption performance is compared with four recently reported encryption schemes in terms of entropy, encryption time, attack-detection capability and classification accuracy. The selected methods for comparison span different encryption paradigms, namely hyperchaotic, memristor chaos, ghost imaging via quaternion transforms, and pixel-adaptive diffusion.

The entropy provided by the HFQE-S scheme is 7.9973 bits/pixel, which is marginally lower than the entropy offered by Hu et al. [6] (7.9984 bits/pixel) and Zhang et al. [5] (7.9986 bits/pixel), within a margin of 0.0011–0.0013 bits/pixel, while HFQE-S performs light-weight block-wise processing on grayscale CT images without requiring quaternion transformation or hardware-based chaotic-circuit realization. Importantly, none of the compared techniques provide a multi-class attack-detection mechanism: Lin et al. [1] and Zhang et al. [5] offer encryption only, with no detection capability, and although Hu et al. [6] and Gong et al. [7] respectively provide authentication and source -verification schemes, neither classifies the type of attack. The comparison also includes the base Hyperchaotic Fibonacci Q-Matrix + HAPCNN-COA framework of Veerasekharreddy et al. [8], which reports an entropy of 7.9971 bits/pixel and an encryption time of 2.58 ms, together with a binary (tamper-vs-clean) detection accuracy of 91.7%, but does not report a per-class classification accuracy for the six-class attack-type setting used in this work. In contrast, the proposed HFQE-S with COA-tuned HAPCNN achieves 7.9973 bits/pixel entropy, 2.87 ms encryption time, and 96.52% six-class attack-classification accuracy, combining competitive encryption randomness with a substantially richer, multi-class detection capability. Table 7 indicates the comparative analysis of proposed work.

Table 7.

Comparative analysis of proposed work.

Encryption TechniqueEntropy (bits/pixel)Encrypt TimeAttack Detection
4D Memristor Chaotic + CNN Key Gen [17]7.99860.23 sNone
CGI + QMPDFrAT + Dual Chaotic [18]7.9984Not reportedAuthentication only
Pixel Adaptive Diffusion + SVDGI [19]~7.998Not reportedSource verification only
Hyperchaotic Fibonacci Q-Matrix + HAPCNN-COA [20]7.99712.58 msCOA-tuned HAPCNN (binary tamper check)
HFQE-S (SHA-256 + Fibonacci Q-Matrix)7.99732.87 msCOA-HAPCNN (6-class) [96.52%]

5. Discussion

These results may be accounted for in two ways. To begin with, the entropy value for the encryption method HFQE-S equals 7.9973, which is close to the theoretically expected value of 8.0. This means that the distribution of the pixels in the encrypted image is nearly uniform. Therefore, this image possesses a good level of resistance to any statistical attack. Moreover, the fact that the PSNR value is infinity means that there was no loss of any data while decrypting the image. This feature is highly required for any medical images.

For the purposes of detection, HAPCNN with no optimizers yields an accuracy of 83.94%, which demonstrates that the architecture itself is able to detect attack instances in noise-like images. With COA optimization, however, this accuracy increases to 96.52%, meaning that this method helps to improve accuracy by over 12 percentage points relative to the base case, over 6 percentage points relative to PSO (89.94%) and over 4 percentage points relative to WOA (92.51%). This is achieved because COA manages to select optimal values of such hyperparameters as learning rate, reconstruction loss weight, dropout and polynomial degree and not due to any changes in the model's architecture itself. Therefore, the value of HAPCNN is its capacity to detect attacks in noise-like images, while the value of COA is its ability to drive the model’s accuracy to higher level.

6. Conclusion

In this study, the authors proposed the AegisSentinel, which is the integrated technique used for the safe transmission of medical imaging information within the IoT-enabled healthcare environment. In this regard, the system is made up of the HFQE-S, HAPCNN, and COA algorithms, in which the HFQE-S algorithm manages to achieve the entropy value close to perfection (7.9973 bits), while encryption and decryption are performed in 2.87 ms and 2.35 ms, respectively, without any loss of data. Moreover, the COA-optimized HAPCNN algorithm achieved an accuracy of 96.52% and an F1 score of 96.49%.

This research paper will illustrate that the transportation of imaging data through the proposed architecture will be secure communication. Future work will involve extending this system for MRI, X-ray, and ultrasound imaging data with the new forms of attacks to the model (PGD, Carlini-Wagner). Different kinds of digital signatures will be analyzed for verifying the identities of both parties of communication.

Acknowledgment

The authors would like to thank Mepco Schlenk Engineering College, Sivakasi, for providing the computational resources and academic support for this research.

Notes

[1] Contributed by Author Contributions

Conceptualization, K.M.S.; methodology, K.M.S., S.S.S., N.G.K., and D.D.N.; software, S.S.S., N.G.K., and D.D.N.; validation, K.M.S. and S.S.S.; formal analysis, K.M.S. and N.G.K.; investigation, S.S.S., N.G.K., and D.D.N.; resources, K.M.S.; data curation, S.S.S., N.G.K., and D.D.N.; writing—original draft preparation, S.S.S., N.G.K., and D.D.N.; writing—review and editing, K.M.S.; visualization, N.G.K. and D.D.N.; supervision, K.M.S.; project administration, K.M.S. All authors have read and agreed to the published version of the manuscript.

[2] Conflicts of interest Conflict of Interest Statement

The authors declare no conflicts of interest.

[3] Data Availability Statement

The data supporting the findings of this study are available from the corresponding author upon reasonable request. The datasets are not publicly available as they are part of the ongoing research work.

8. Author Biographies (Optional)

graphic/j_ias-2026-0017_ingr_013.png K. Muthamil Sudar received the B.Tech. degree in computer science and engineering from Kalasalingam University, Srivilliputhur, Tamil Nadu, India, in 2014, the M.E. degree in computer science and engineering from RVS College of Engineering and Technology, Coimbatore, Tamil Nadu, India, in 2016, and the Ph.D. degree in computer science and engineering from Kalasalingam Academy of Research and Education, Srivilliputhur, Tamil Nadu, India, in 2021. He is currently an Assistant Professor with the Department of Computer Science and Engineering, Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India. His research interests include network security, software-defined networking, and machine learning. He has published more than 55 research articles in reputed international journals and leading conferences.

graphic/j_ias-2026-0017_ingr_014.png Sai Shobana Sri is currently pursuing the B.E. degree in computer science and engineering with the Department of Computer Science and Engineering, Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India. Her research interests include image processing, steganography, and cybersecurity.

graphic/j_ias-2026-0017_ingr_015.png Nigila G. K is currently pursuing the B.E. degree in computer science and engineering with the Department of Computer Science and Engineering, Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India.Her research interests include machine learning, network security, and web development.

graphic/j_ias-2026-0017_ingr_016.png Durga Devi N is currently pursuing the B.E. degree in computer science and engineering with the Department of Computer Science and Engineering, Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India. Her research interests include deep learning and secure data transmission.

DOI: https://doi.org/10.2478/ias-2026-0017 | Journal eISSN: 1554-1029 | Journal ISSN: 1554-1010
Language: English
Page range: 334 - 352
Published on: Aug 7, 2026
Published by: Cerebration Science Publishing Co., Limited
In partnership with: Paradigm Publishing Services
Publication frequency: Volume open

© 2026 K. Muthamil Sudar, Sai Shobana Sri, Nigila G K, Durga Devi N, published by Cerebration Science Publishing Co., Limited
This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License.