Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
OTE-SLAM: An Object Tracking Enhanced Visual SLAM System for Dynamic Environments
Previous Article in Journal
A Novel FDLSR-Based Technique for View-Independent Vehicle Make and Model Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel 3D Reversible Data Hiding Scheme Based on Integer–Reversible Krawtchouk Transform for IoMT

1
Dhar El Mahrez Faculty of Science, Sidi Mohamed Ben Abdellah-Fez University, Fez 30000, Morocco
2
National School of Applied Sciences, Sidi Mohamed Ben Abdellah-Fez University, Fez 30000, Morocco
3
Department of Computer Science, Faculty of Computer Science and Telecommunications, Cracow University of Technology, Warszawska 24, 31-155 Krakow, Poland
4
Institute of Theoretical and Applied Informatics, Polish Academy of Sciences, Bałtycka 5, 44-100 Gliwice, Poland
5
Information Countermeauser Technique Institute, Harbin Institute of Technology, School of Cyberspace Science, Faculty of Computing, Harbin 150001, China
6
Computer Science Department, Community College, King Saud University, Riyadh 11437, Saudi Arabia
7
Department of Mathematics and Computer Science, Faculty of Science, Menoufia University, Shebeen El-Kom 32511, Egypt
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(18), 7914; https://doi.org/10.3390/s23187914
Submission received: 7 August 2023 / Revised: 10 September 2023 / Accepted: 11 September 2023 / Published: 15 September 2023
(This article belongs to the Special Issue Trusted ICT in the Secure Human-Centered Computing (HCC))

Abstract

:
To avoid rounding errors associated with the limited representation of significant digits when applying the floating-point Krawtchouk transform in image processing, we present an integer and reversible version of the Krawtchouk transform (IRKT). This proposed IRKT generates integer-valued coefficients within the Krawtchouk domain, seamlessly aligning with the integer representation commonly utilized in lossless image applications. Building upon the IRKT, we introduce a novel 3D reversible data hiding (RDH) algorithm designed for the secure storage and transmission of extensive medical data within the IoMT (Internet of Medical Things) sector. Through the utilization of the IRKT-based 3D RDH method, a substantial amount of additional data can be embedded into 3D carrier medical images without augmenting their original size or compromising information integrity upon data extraction. Extensive experimental evaluations substantiate the effectiveness of the proposed algorithm, particularly regarding its high embedding capacity, imperceptibility, and resilience against statistical attacks. The integration of this proposed algorithm into the IoMT sector furnishes enhanced security measures for the safeguarded storage and transmission of massive medical data, thereby addressing the limitations of conventional 2D RDH algorithms for medical images.

1. Introduction

Integer-to-integer transforms (such as integer wavelet transforms [1], integer DCT [2], and integer color transform [3]) are an extension of traditional floating-point transforms. Unlike floating-point transforms that produce real-valued transformed coefficients, integer-to-integer transforms are specifically designed to directly generate integer-valued transformed coefficients. By exclusively operating on integer data throughout the entire transform process, including decomposition and reconstruction stages, these transforms eliminate round-off errors, computational overhead, and the need for quantization. Furthermore, they provide precise control over numerical accuracy. Integer-to-integer transforms have been studied and applied in lossless image compression [2,3,4,5,6] and image reversible data hiding (RDH) [7,8,9,10]. Such applications require accurate reconstruction of data without any loss of information. Various alternative variants of integer transforms have been investigated and developed, such as the integer discrete Fourier transform [11], integer fast Fourier transform [12], TS transform [5], the S + P transform [4], integer Tchebichef transform [6], and integer color transform [3].
In [13], the authors introduce for the first time another orthogonal transform called the Krawtchouk transform for digital image processing. This transform maps the input image data from the spatial domain to a different domain, called the Krawtchouk domain, defined by the Krawtchouk polynomials (kernel function of the transform). This transformed domain can provide valuable information about underlying patterns, structures, or features present in the data that may not be readily apparent in the original spatial representation. Furthermore, through the manipulation of control parameters within the framework of the Krawtchouk transform, it becomes possible to effectively capture and characterize local features present within an image. The Krawtchouk transform has demonstrated its utility in a range of image processing applications, including image reconstruction, image watermarking, and image zero-watermarking [13,14,15,16,17,18,19]. The authors in [13,14] have evaluated the Krawtchouk transform in image reconstruction and demonstrated that it outperforms several transforms (Tchebichef, Legendre, Zernike, and pseudo-Zernike transforms) in terms of reconstruction error. The authors in [15] proposed a geometrically invariant image watermarking scheme based on Krawtchouk transform. The authors in [17] proposed a Krawtchouk-transform-based watermarking scheme that automatically selects the optimum local region of the carrier image and embeds the watermark in the Krawtchouk domain. The authors in [19] proposed a robust blind zero-watermarking algorithm based on Krawtchouk radial transform, where the kernel function is the product of the Krawtchouk polynomials and an exponential function. The authors in [20] derived a new octonion transform called octonion Krawtchouk moment (OKM) transform that allows, in a compact and holistic way, global features to be simultaneously extracted from a set of seven images. And based on this transform, the authors introduced a new local zero-watermarking method for copyright protection of multiple medical images.
When applying the Krawtchouk transform on digital images that are represented as discrete pixel intensities (integer values), the output Krawtchouk coefficients are represented as floating-point numbers (i.e., real numbers) rather than integers (Figure 1). Consequently, the Krawtchouk transform is associated with some major problems, which are listed below:
  • The original image cannot be perfectly reconstructed from the Krawtchouk coefficients, even if the Krawtchouk transform itself is orthogonal, due to rounding errors associated with the limited representation of significant digits. This can have an impact on the performance of lossless applications requiring exact calculations;
  • Floating-point Krawtchouk coefficients require a more complex representation than integer image pixels (spatial domain), which means a higher cost in terms of memory and computing resources, especially when input images are large or voluminous (3D images). This can impact the performance of applications that require fast and efficient calculations;
  • The Krawtchouk transform is unsuitable for applications where exact integer arithmetic is desired, such as lossless compression, RDH, digital communication systems, and embedded systems with limited precision.
Our primary goal in this paper is to propose a framework for deriving the integer–reversible Krawtchouk transform (IRKT) by using a SERM factorization [21] of the polynomial Krawtchouk matrix. The Krawtchouk polynomial matrix undergoes a decomposition process, during which it is decomposed into a series of Single-row Elementary Reversible Matrices (SERMs) and a permutation matrix. This decomposition is the key step in calculating the IRKT and its inverse (iIRKT). Instead of directly using the Krawtchouk polynomial matrix itself, SERMs and the permutation matrix are used to implement the Krawtchouk transform. Thanks to this decomposition, IRKT and iIRKT can be applied without generating computational errors or losing information (Figure 2).
As an illustration of the IRKT applications, an RDH algorithm for secure storage and transmission is investigated in the context of the Internet of Medical Things (IoMT). In the IoMT field (see the related works in Section 4), there is a massive amount of information that requires secure transmission and storage. The substantial volume of data renders conventional 2D RDH algorithms unsuitable due to their limited embedding capacity. This is attributed to the fact that 2D RDH solely employs a 2D image as a carrier, within which additional/secret data are concealed inside pixel intensities. To achieve heightened capacity, the most straightforward approach involves segmenting the additional/secret data and embedding them into several 2D medical images via the application of a 2D RDH algorithm, such as [22,23]. However, this strategy brings about disadvantages, including the iterative application of computational processes, resulting in notable time overhead, as well as the absence of compact storage and/or unified transmission of the additional/secret data as a single entity. An effective solution emerges in the development of a 3D RDH algorithm that utilizes 3D medical images as carriers, particularly when significant information needs to be sent or stored [24]. Through embracing the third dimension, 3D RDH techniques harness the potential of 3D image stacks, thereby facilitating the seamless embedding of concealed data across the volume. This eliminates the necessity for repetitive computations and introduces a streamlined mechanism for accommodating and transmitting the additional/secret data as a cohesive unit.
In recent years, with the development of 3D technology, more and more 3D medical images have been created and applied to each IoMT application, and also spread in a large amount by the network. Numerous 3D irreversible methods [25,26,27,28] and 3D RDH methods [29,30,31] have been proposed based on 3D image models, specifically using polygon mesh representation. However, these methods are not well suited for 3D medical images due to the limitations of polygon mesh representation, which often lack the ability to explicitly capture the intricate interior information present in such images.
Our second goal in this paper is to propose a novel 3D RDH algorithm for 3D grayscale medical images based on data integrity preservation and the reversibility of the proposed IRKT. The histogram-modification-based spatial approach described in [32] has been applied in the IRKT domain for the 3D RDH of medical images. The IKRT-based 3D RDH algorithm exhibits significant potential in terms of high embedding capacity, high additional data imperceptibility, and high robustness against statistical attacks.
As a summary, our main contributions are shown as follows:
  • A precise, exact integer representation in the Krawtchouk domain is provided without rounding errors or accuracy limitations;
  • IRKT guarantees that transformed coefficients remain integers (without the need for quantization), enabling lossless data reconstruction;
  • The proposed IRKT can be easily generalized to accommodate 2D and 3D data representations;
  • A novel 3D RDH algorithm suitable for 3D medical images is proposed, to the best of our knowledge, for the first time;
  • Embedding data into 3D medical images does not increase their original size, thus optimizing infrastructure and maximizing resource utilization in the IoMT;
  • Medical images are recovered without any loss or damage after extracting additional data;
  • The embedding capacity and quality of the 3D stego image (image after data embedding) can be adjusted using a threshold-based embedding technique.
The rest of the paper is organized as follows. Some preliminaries on the Krawtchouk transform are provided and its limit is discussed in Section 2. In Section 3, the detailed procedure for deriving IRKT and some properties are provided. The 3D RDH algorithm for 3D grayscale medical images is presented in Section 4. Experiments designed to demonstrate its effectiveness are carried out in Section 5, and Section 6 concludes the paper.

2. Preliminaries

The Krawtchouk transform of a 1D signal f(x) of length N is defined as [13]
Q n = x = 0 N 1 K n ( x ; p , N 1 ) f ( x ) ,   n = 0 , 1 , , N 1
where Qn is the Krawtchouk transform of f(x), and Kn(x; p, N) is the Krawtchouk polynomials of degree n:
K n ( x ; p , N 1 ) = 2 F 1 ( n , x ; N + 1 ; 1 p )   ( 1 ) n ( N + 1 ) n ( N 1 x ) p n + x ( 1 p ) N x n 1 n !
where x, n = 0, 1, …, N − 1, p ( 0 , 1 ) is the Krawtchouk parameter controlling the shape of the polynomial, 2F1 represents the hypergeometric function, which can be defined as
2 F 1 ( a , b ; c ; z ) = k = 0 ( a ) k ( b ) k ( c ) k z k k !
and (a)k is the Pochhammer symbol, defined as
( a ) 0 = 1 ( a ) k = a ( a + 1 ) ( a + 2 ) ( a + k 1 ) ,   k 1
The Krawtchouk polynomials exhibit a distinct property of orthogonality, as stated in [13]
x = 0 N 1 K n ( x ; p , N 1 ) K m ( x ; p , N 1 ) = δ n , m
The property of orthogonality gives rise to the subsequent inverse Krawtchouk transform as follows [13]:
f ( x ) = n = 0 N 1 Q n K n ( x ; p , N 1 )
In the case of an image f(x, y) of size N × N, the two-dimensional Krawtchouk transform of order n + m and its inverse can be expressed mathematically as follows:
Q n , m = x = 0 N 1 y = 0 N 1 K n ( x ; p , N 1 ) K m ( y ; p , N 1 ) f ( x , y )
f ( x , y ) = n = 0 N 1 m = 0 N 1 Q n , m K n ( x ; p , N 1 ) K m ( y ; p , N 1 )
In the scenario of an N × N × N volumetric image f(x, y, z), the three-dimensional Krawtchouk transform and its inverse can be mathematically represented as follows [33]:
Q n , m , = x = 0 N 1 y = 0 N 1 z = 0 N 1 K n ( x ; p , N 1 ) K m ( y ; p , M 1 ) K ( z ; p , N 1 ) f ( x , y , z )
f ( x , y , z ) = n = 0 N 1 m = 0 N 1 = 0 N 1 Q n , m , K n ( x ; p , N 1 ) K m ( y ; p , N 1 ) K ( z ; p , N 1 )
When original integer data are transformed into the Krawtchouk domain, they are represented by a set of floating-point coefficients (real numbers) resulting from the projection of data onto the Krawtchouk polynomials. However, due to factors such as numerical precision or approximation errors during the Krawtchouk transform, the original data cannot be perfectly reconstructed from these coefficients. For example, Figure 1 shows an input data block, its Krawtchouk transform, the reconstructed block via the inverse Krawtchouk transform, and the absolute difference between the original and reconstructed blocks. The calculation is based on the IEEE Standard 754 for double precision. It is clear that the reconstructed block is not exactly the original input block, which means that information is lost in this process, even though the Krawtchouk transform itself is orthogonal (Equation (5)). In lossless image applications like reversible data hiding (RDH), lossless compression, or lossless watermarking, preserving the exact values and maintaining data integrity is crucial. Therefore, the Krawtchouk transform is not suitable for these applications where any information loss is not acceptable. In order to overcome this limitation, we propose in the next section an integer and reversible version of the Krawtchouk transform that is more suitable for lossless image applications.

3. Construction of the Integer–Reversible Krawtchouk Transform

In this section, we propose an integer–reversible Krawtchouk transform (IRKT), which can be applied on integer data and generate integer Krawtchouk coefficients (integer-based representation), and its inverse IRKT can exactly restore the original data from the Krawtchouk coefficients. The Krawtchouk polynomial matrix undergoes a decomposition process, during which it is decomposed into a series of Single-row Elementary Reversible Matrices (SERMs) and a permutation matrix. This decomposition is the key step in calculating the IRKT and its inverse (iIRKT). Instead of directly using the Krawtchouk polynomial matrix itself, SERMs and the permutation matrix are used to implement the Krawtchouk transform. Thanks to this decomposition, IRKT and iIRKT can be applied without generating computational errors or loss of information.

3.1. SERM Factorization of the Krawtchouk Polynomial Matrix

Based on [21], it has been established that for any nonsingular N × N matrix A, a unit SREM factorization of A = SNS2S1 exists if and only if all the minors of the leading principal submatrices of A are ones, where Sj(j = 0, 1, 2,…, N) represents the SERMs involved in the factorization. In particular, if the determinant of A is an integer number, then an N + 1 SERM factorization A = PSNS2S1S0 can be obtained. This factorization involves the use of permutation matrix P and unit SERMs Sj(j = 0, 1, 2, …, N).
The SERM factorization of the Krawtchouk polynomial matrix can be obtained as follows.
The matrix representation of the 1D Krawtchouk transform, which applies to a signal f(x) of length N as defined in Equation (1), can be expressed as
Q = K f
where K is the Krawtchouk polynomial matrix defined as follows:
K n , x = K n ( x ; p , N 1 ) , n , x = 0 , 1 , , N 1
Typical Krawtchouk polynomial matrices of size 4 × 4 and 8 × 8 for Krawtchouk transform and its inverse are as follows:
K 4 × 4 = ( 0.3536 0.6124 0.6124 0.3536 0.6124 0.3536 0.3536 0.6124 0.6124 0.3536 0.3536 0.6124 0.3536 0.6124 0.6124 0.3536 ) K 8 × 8 = ( 0.0884 0.2339 0.4050 0.5229 0.5229 0.4050 0.2339 0.0884 0.2339 0.4419 0.4593 0.1976 0.1976 0.4593 0.4419 0.2339 0.4050 0.4593 0.0884 0.3423 0.3423 0.0884 0.4593 0.4050 0.5229 0.1976 0.3423 0.2652 0.2652 0.3423 0.1976 0.5229 0.5229 0.1976 0.3423 0.2652 0.2652 0.3423 0.1976 0.5229 0.4050 0.4593 0.0884 0.3423 0.3423 0.0884 0.4593 0.4050 0.2339 0.4419 0.4593 0.1976 0.1976 0.4593 0.4419 0.2339 0.0884 0.2339 0.4050 0.5229 0.5229 0.4050 0.2339 0.0884 )
The matrix K satisfies the following properties:
  • Property 1: K is orthogonal, that is, KKT = KTK = I, KT is the transpose of K, and I is the identity matrix [13] of the appropriate size;
  • Property 2: K is nonsingular, that is, KK−1 = K−1K = I;
  • Property 3: The determinant of K is det(K) = 1;
  • Property 4: The minors of the leading principal submatrices of K are all ones.
According to the Corollary 6 in [21], the previous properties of K lead to the following SERM factorization:
K = P S N S 2 S 1 S 0
where P is the permutation matrix of K, and Sj(j = 0, 1,…, N) are SERMs. An SERM represents a square matrix of dimension N × N, wherein the diagonal is composed of integer factors, and solely one row possesses off-diagonal elements that are not entirely composed of zeros. The SERM matrix Sj can be expressed as
S 0 = I + e N s 0 T S j = I + e j s j T   ,   j = 1 , 2 , , N
The inverse of Sj can be expressed as follows:
S 0 1 = S 0 = I + e N s 0 T S j 1 = I e j s j T ,   j = 1 , 2 , , N
where I is the identity matrix of size N × N, ej refers to the column vector corresponding to the j-th column of I, and sj is a vector whose j-th element is zero.
The permutation matrices of K for sizes 4 × 4 and 8 × 8 are as follows:
P 4 × 4 = (   0   1   0   0   1   0   0   0   0   0   0   1   0   0   1   0 ) P 8 × 8 = (   0   0   0   1   0   0   0   0   0   0   1   0   0   0   0   0   0   0   0   0   0   0   1   0   0   0   0   0   1   0   0   0   1   0   0   0   0   0   0   0   0   0   0   0   0   1   0   0   0   1   0   0   0   0   0   0   0   0   0   0   0   0   0   1 )
The row vectors for computing the SERMs of K for sizes 4 × 4 and 8 × 8 are listed in Table 1 and Table 2, respectively.

3.2. Integer–Reversible Krawtchouk Transform

By utilizing the SERM factorization of the Krawtchouk polynomial matrix K, as discussed in the previous subsection, we can establish an integer–reversible version of the Krawtchouk transform (IRKT).
For an integer vector f = (f0, f1,…, fN)T of length N, the 1D IRKT can be defined as
Q = I R K T ( f ) = P S N S 1 S 0 f
where . is the rounding arithmetic, Sj(j = 0, 1,…, N) are SERMs of the N × N matrix K, and P is the permutation matrix of K.
The corresponding inverse IRKT (iIRKT) can be expressed as
f = i I R K T ( Q ) = S 0 1 S 1 1 S N 1 P T Q
To compute the 2D IRKT for an N × N block f = { f ( x , y ) } x , y = 0 x , y = N 1 , the 1D IRKT is firstly applied to each column of the block. This operation transforms each column individually. After that, the 1D IRKT is then applied to each row of the resulting transform obtained from the previous step. This process ensures that the transform is performed on both the columns and rows, completing the 2D IRKT for the N × N block. The mathematical formula of 2D-IRKT can be represented as
Q = I R K T ( ( I R K T ( f ) ) T )
and the 2-D inverse IRKT as
f = i I R K T ( ( i I R K T ( Q ) ) T )
Figure 2 shows a block of input data (the same block used in Figure 1), its IRKT, the block reconstructed via iIRKT, and the absolute difference between the original and reconstructed blocks. It is clear that the reconstructed block is exactly the original input block, meaning that the IRKT can provide an exact reconstruction of the original data.
The IRKT can be easily generalized to 3D situations as follows. To compute the 3D IRKT for an N × N × N volumetric cuboid, a stepwise approach is employed as follows: (1) First, the 1D IRKT is applied to each column of every 2D plane in the cuboid. This operation transforms each column within each 2D plane individually. (2) Next, the 1D IRKT is applied to each row of every 2D plane in the transformed cuboid obtained from the previous step. This operation transforms each row within each 2D plane. (3) Finally, the 1D IRKT is applied to each depth plane of the cuboid. This operation transforms each individual depth plane.
Detailed descriptions for implementing 3D IRKT and inverse 3D IRKT are given in Algorithm 1 and Algorithm 2, respectively.
Algorithm 1 3D Integer–Reversible Krawtchouk Transform.
Input: C a volumetric cuboid of size N × N × N, P the permutation matrix of K, Sj(j = 0, 1,…, N) the SERMs of K.
Output: Transformed cuboid Q of size N × N × N
Initialize empty matrices Q, Q1, Q2 of size N × N × N.
for i = 1 to N do
   A = PSN⌊…⌊ S1S0C(:,:,i) ⌋⌋…⌋;//C(:,:,i) is the i-th plane along the z-axis.
   Q1(:,:,i) = PSN⌊…⌊ S1S0AT⌋⌋…⌋;//AT is the transpose of A.
end
//Transpose the 3D matrix Q1 into another 3D matrix Q2.
for i = 1 to N do
   for j = 1 to N do
     Q2(j,:,i) = Q1(i,:,j);
   end
end
//Generate the IRKT matrix Q
for i = 1 to N do
   A = PSN⌊…⌊S1S0Q2(:,:,i) ⌋⌋…⌋;
   Q(:,:,i) = PSN⌊…⌊S1S0AT⌋…⌋;
end
Return the IRKT matrix Q.
Algorithm 2 Inverse 3D Integer–Reversible Krawtchouk Transform.
Input: IRKT matrixQ of size N × N × N,P the permutation matrix of K, Sj(j = 0, 1,…, N) the SERMs of K.
Output: Reconstructed cuboid R of size N × N×N.
Initialize empty matrices R, R1, R2 of size N × N × N.
Compute S j 1 the inverse of Sj.
for i = 1 to N do
   A= ⌊ S 0 1 S 1 1 …⌊ S N 1 PTQ(:,:,i)⌋…⌋⌋;//PT is the transpose of P.
   R1(:,:,i) = ⌊ S 0 1 S 1 1 …⌊ S N 1 PTAT⌋…⌋⌋;
end
for i = 1 to N do
   for j = 1 to N do
     R2(j,:,i) = R1(i,:,j);
   end
end
for i = 1 to N do
   A= ⌊ S 0 1 S 1 1 …⌊ S N 1 PTR2 (:,:,i)⌋…⌋⌋;
   R(:,:,i) = S 0 1 S 1 1 …⌊ S N 1 PTAT⌋…⌋⌋;
end
Return the Reconstructed cuboid R.

4. Application in RDH for 3D Medical Images

The proposed IRKT can be applied in many lossless applications such as lossless compression, lossless encryption, and lossless watermarking. Due to the limited paper length, the reversible data hiding (RDH) application is only discussed in the context of the Internet of Medical Things (IoMT).
Schematically, RDH is a practical lossless technique that can be mainly used to ensure secure transmission and storage of data [32]. This can be accomplished by embedding additional data (secret or sensitive data), in an imperceptible manner, into an image (carrier image) while ensuring the complete recovery of the original image without any loss after data extraction [34,35,36,37,38,39,40]. RDH techniques can be classified into three categories: space-based techniques, compression-based techniques, and transform-domain-based techniques [9]. Since space- and compression-based techniques embed the additional data into the carrier image directly, they are relatively weak under statistical attacks [9]. In contrast, transform-domain-based methods generally provide an additional layer of security and are relatively more resistant to statistical attacks [41].
Recently, a number of data hiding methods based on moment transforms [42,43,44,45,46] have been proposed. However, it is important to note that these methods are irreversible since they embed additional data into the real-valued moment transform domain. Consequently, these methods result in a loss of information from the original image, making it nonrecoverable. This is why the integer–reversible transforms have been more widely studied in RDH applications. Noteworthy examples of such transforms include the integer wavelet transforms [8,9] and integer DCT [10].
Among the RDH strategies that are simple to implement and require a short execution time, histogram modification is proposed by Ni et al. [32]. The authors embedded the additional data by modifying the histogram of the carrier image. In addition to its simplicity, this strategy produces a stego image of high quality. Ni’s strategy has been successfully adopted by several methods in the transform domain [9,10,47,48,49,50].
In the IoMT field (See Section 4.1), the amount of information requiring secure transmission and storage is gigantic/huge, making the conventional 2D RDH algorithm for medical images unsuitable due to its limited embedding capacity. Indeed, the data that can be embedded into the carrier medical image are the Electronic Patient Record (EPR), diagnostic annotations, measurements, or any other additional details that can contribute to more accurate and informed diagnoses and decisions.
To achieve high capacity, the most straightforward practice is to divide the additional data, and then embed them in several 2D medical images by applying the 2D RDH algorithm, such as [22,23]. However, this approach has the disadvantages of applying the computation scheme several times, adding significant time overhead, and the additional data are not stored and/or sent compactly in a single entity. An effective solution is to develop a 3D RDH algorithm that uses 3D medical images as carriers when a large amount of information needs to be sent or stored [24].
In recent years, with the development of 3D technology, more and more 3D medical images are created and applied to each IoMT application, and are also spread in a large amount by the network.
Numerous 3D irreversible methods [25,26,27,28] and 3D RDH methods [29,30,31] have been proposed based on 3D image models, specifically using polygon mesh representation. However, these methods are not well suited for 3D medical images. While polygonal models can effectively represent the surfaces of objects and scenes in 3D space, they often lack explicit capture of interior information. In contrast, voxel-based representation is commonly employed in 3D medical imaging, where each voxel represents a small volume element and stores information about the object at that specific location. This representation allows for the inclusion of both exterior and interior details, making it more suitable for accurately representing the complex structures and characteristics found in 3D medical images.
In this section, we propose a new 3D RDH scheme for 3D medical images based on data integrity preservation and reversibility of the proposed IRKT.

4.1. The Usefulness of 3D RDH in IoMT

The IoMT is a burgeoning field that encompasses the integration of medical devices, sensors, and healthcare systems with the internet. It enables the seamless exchange, analysis, and management of healthcare data, revolutionizing the delivery of medical services [51]. The IoMT has the potential to enhance patient care, facilitate remote monitoring, enable real-time interventions, and improve healthcare outcomes. However, the extensive connectivity and data sharing in IoMT also raise concerns regarding privacy, data security, data transmission and storage, and efficient data utilization. The RDH technique allows additional information to be hidden within medical images without compromising their visual quality, ensuring secure storage and transmission while preserving patient privacy.
The utilization of 3D RDH techniques, employing 3D medical images as carrier images, provides distinct advantages in IoMT applications:
  • Due to the substantial volume of data generated by IoMT, storage and transmission may encounter limitations. By storing and transmitting a significant amount of medical data within the same carrier image, efficient resource utilization and management can be achieved within IoMT applications;
  • Protecting patient privacy is crucial. The 3D RDH can selectively embed confidential patient identifiers or sensitive information within 3D medical images, ensuring that only authorized personnel can access these details;
  • By embedding additional diagnostic information into 3D medical images, the accuracy of diagnoses can be enhanced, thereby granting healthcare professionals extended access to relevant data during remote consultations;
  • In telemedicine scenarios, real-time interaction is limited. The 3D RDH can hide additional explanations, annotations, or visual cues within 3D medical images, offering remote healthcare professionals a more detailed understanding of the patient’s condition;
  • In the context of IoMT, medical data are frequently employed for research and analysis purposes. The adoption of 3D RDH alongside 3D medical images enables the integration of research data, metadata, or annotations into the images. This facilitates comprehensive data analysis, data exploration, and collaboration among researchers, thereby enabling them to gain deeper insights into medical conditions and treatment outcomes.
In the following, the histogram-modification-based spatial approach described in [32] has been applied in the IRKT domain for 3D RDH of medical images. The IRKT-based 3D RDH algorithm has the following features: (1) it is an integer-value-to-integer-value mapping in both data embedding and data extraction processes; (2) it exhibits significant potential in terms of high embedding capacity, high additional data imperceptibility, and high robustness against statistical attacks; (3) embedding data into 3D medical images does not increase their original size, thus optimizing infrastructure and maximizing resource utilization in the IoMT; (4) medical images are recovered without any loss or damage after extracting additional information; (5) it is well suited for high-quality 3D medical images; (6) the embedding capacity and quality of the 3D stego image (image after data embedding) can be adjusted using a threshold-based embedding technique; (7) partitioning the 3D carrier image into smaller cuboids offers practical advantages in terms of implementation and algorithm complexity reduction.
The embedding, extraction, and restoration procedures are described in the following subsections.

4.2. Embedding Procedure

Consider a 3D medical image f = { f ( x , y , z ) } x , y , z = 0 x , y , z = N 1 with N × N × N voxels. The additional data are embedded into the carrier image f, as depicted in Figure 3, through a series of steps, which are outlined below:
Step 1: Encode the additional data into a binary representation. We refer to the resulting binary vector of length L as V = { V ,   = 0 , 1 , , L 1 } . Here, the additional data are the Electronic Patient Record (EPR), diagnostic annotations, measurements, or any other additional details that can contribute to more accurate and informed diagnoses and decisions.
Step 2: The 3D carrier image f is divided into n × n × n nonoverlapping volumetric cuboids, and the 3D IRKT (Algorithm 1) is computed for each cuboid. Then, the IRKT cuboids are collected into a single 3D matrix of size N × N × N, denoted by Q.
Step 3: Create an empty space in the matrix Q to insert the additional data. Here, the term “empty space” refers to the values of the IRKT coefficients being modified by replacing their original values. The size of the empty space determines the amount of data that can be inserted. To control and adjust the amount of data to be inserted, we employ a threshold-based embedding technique. A positive integer value, T > 0 , is used to adjust the embedding capacity and the quality of the stego image.
Let us consider Q(i, j, k) as the value of the IRKT coefficient at coordinate (i, j, k), and its modified value denoted by Q′(i, j, k). The modification follows the rules
Q ( i , j , k ) = { Q ( i , j , k ) + T   if   Q ( i , j , k ) T Q ( i , j , k ) T   if   Q ( i , j , k ) < T Q ( i , j , k )   otherwise
where Q = { Q ( i , j , k ) ,   0 i , j , k < N } represents the modified IRKT matrix. From Equation (20), all the IRKT coefficients with values satisfying Q ( i , j , k ) T are increased by T, while those with values satisfying Q ( i , j , k ) < T are decreased by T. The values in the interval [−T, T) are not modified.
By applying this step, an empty space, in [−2T, −T) and [T, 2T), is created within the matrix Q′, allowing for the subsequent embedding of the additional data.
Step 4: In the proposed method, the number of additional data bits to be embedded is equal to the number of IRKT coefficients in the range [−T, T). The matrix Q′ is scanned, and when IRKT coefficients, whose value Q′(i, j, k) is in the interval [−T, T), are encountered, they are modified according to the embedded bit values. The rule for embedding { V } is as follows:
Q s ( i , j , k ) = { Q ( i , j , k ) + T   if   0 Q ( i , j , k ) < T   &   V = 1 Q ( i , j , k ) T   if   T Q ( i , j , k ) < 0   &   V = 1 Q ( i , j , k )   if   T Q ( i , j , k ) < T   &   V = 0 Q ( i , j , k )   if   Q ( i , j , k ) < T o r   Q ( i , j , k ) T
where Q s = { Q s ( i , j , k ) ,   0 i , j , k < N } represents the stego IRKT matrix.
Through this step, the empty space in the intervals [−2T, −T) and [T, 2T) is effectively filled within the matrix Qs. Negative IRKT coefficients in the interval [−T, T) are moved into the empty interval [−2T, −T) and positive ones into the empty interval [T, 2T) if V = 1 . If not, the IRKT coefficients remain unchanged.
Step 5: The 3D inverse IRKT (Algorithm 2) is applied to each stego cuboid of size n × n × n within the matrix Qs to obtain the stego image fs.
The output stego image, which contains the additional data, can be effectively transmitted or securely stored within the IoMT ecosystem. In both cases, the additional data remain in the stego image. Authorized parties can extract these data and recover the original image from the stego image by running the decoding algorithm presented in the next section.

4.3. Additional Data Extraction and Carrier Image Restoration

Step 1: The 3D stego image fs is divided into nonoverlapping cuboids of size n × n × n. The resulting IRKT cuboids are combined into a matrix denoted by Qs.
Step 2: During the extraction of additional data, Qs is scanned as in the embedding procedure, and then the sequence { V } is extracted from Qs using the following rule:
V = { 0   if   T Q s ( i , j , k ) < T   1   if   2 T Q s ( i , j , k ) < T   o r   T Q s ( i , j , k ) < 2 T
Step 3: During the restoration of the carrier image, the IRKT coefficients that were modified during embedding steps 3 and 4 are reverted to their original values. This is accomplished as follows:
Q ( i , j , k ) = { Q s ( i , j , k ) + T   if   Q s ( i , j , k ) < T Q s ( i , j , k ) T   if   Q s ( i , j , k ) 0 Q s ( i , j , k )   otherwise
where Q = { Q ( i , j , k ) ,   0 i , j , k < N } represents the original carrier IRKT matrix.
Step 4: Finally, the inverse IRKT is applied on Q to restore the original carrier image.
The efficiency of the IRKT-based 3D RDH algorithm in terms of embedding capacity and 3D stego image quality can be improved by adjusting the threshold T in the IRKT domain, the Krawtchouk parameter p ( 0 , 1 ) , and the cuboid size used for 3D IRKT computation. The following section discusses the influence of these parameters on the algorithm.

5. Experimental Results

A series of experiments were conducted to assess the efficiency of the 3D RDH algorithm based on IRKT for medical images. The algorithm was implemented using Matlab R2022b on a computing system equipped with an AMD Ryzen 5 5600U CPU, which has six cores running at 2.30 GHz and 8 GB of RAM. In the conducted experiments, a set of 208 grayscale 3D medical images obtained from the LIDC-IDRI [52], OpenNeuro [53], IRCAD [54], and BraTS 2021 [55] datasets served as the carrier images. Simultaneously, a randomly generated sequence of binary digits {0, 1} was employed as the additional data intended for embedding. More details on the test datasets are provided in Table 3. This dataset comprises 3D medical images encompassing both CT and MRI modalities, including both signed and unsigned images, with a depth of 16 bits (indicating high quality). The images are stored in DICOM and NIfTI formats. On average, it took 4.5 s, 8.7 s, 13.65 s, and 21.41 s to embed additional data bits into 3D medical images of sizes 96 × 96 × 96, 128 × 128 × 128, 160 × 160 × 160, and 192 × 192 × 192, respectively.
This IRKT-based RDH algorithm is compared with the spatial approach reported in [29,30,31,32,56], as well as the RDH using integer DCT [10], integer wavelet transforms [8,9], standard DCT [57], and standard DWT [58].
The algorithm’s performance in terms of additional data imperceptibility is quantified by the peak signal-to-noise ratio (PSNR) defined as follows.
PSNR ( dB ) = 10 log 10 ( M a x 2 MSE )
where Max is the maximum possible value for a voxel in the image, and MSE is the mean squared error between the carrier and stego images:
MSE = 1 M × N × K x = 1 M 1 y = 1 N 1 z = 1 K 1 [ f ( x , y , z ) f s ( x , y , z ) ] 2
where f and fs are the carrier image and stego one of size M × N × K, respectively.
The amount of information carried by the 3D stego image, known as embedding capacity (EC), is measured in bits per voxel (bpv), and defined as follows:
EC ( bpv ) = Additional   data   size   ( bits ) Total   number   of   voxels   in   3 D   carrier   image
In the following experiments, the cuboid size used for 3D IRKT computation is 4 × 4 × 4, unless otherwise specified in the paper.
The first test aimed to demonstrate the significant advantage of the proposed IRKT over the classical Krawtchouk transform. As explained in Section 3, the IRKT is a lossless and reversible transform, ensuring that images reconstructed using the IRKT are indistinguishable from the original images. Table 4 presents the average MSE and PSNR values for medical images reconstructed using the IRKT. These results hold true regardless of the value assigned to the Krawtchouk polynomial parameter p within the range of p ( 0 , 1 ) . Additionally, we compared the IRKT with the classical Krawtchouk transform when p = 0.5, as this value captures the maximum image information [13]. Notably, the calculations for these transforms were performed using the long format specified by the IEEE® standard for floating-point arithmetic. Table 4 clearly shows that medical images reconstructed using the proposed IRKT are exactly identical to the original images. In contrast, the classical Krawtchouk transform lacks reversibility, preventing an exact reconstruction of the original images.
In the second test, we investigated the impact of varying the Krawtchouk polynomial parameter p on the embedding capacity and quality of stego images. For each test image, the IRKT was applied, and then the additional data were inserted into the IRKT domain using the threshold T = 1. Subsequently, the stego image was generated by applying the inverse IRKT. In this test, the IRKT was implemented by increasing the Krawtchouk polynomial parameter p from 0.1 to 0.9 with an increment of 0.1. Figure 4 shows the average PSNRs and ECs of stego images for these different Krawtchouk polynomial parameters. It can be seen from Figure 4 that the parameters p = 0.1, 0.5, 0.8 and 0.9 are comparable and give high additional data imperceptibility (high PSNR) compared to the other p values. On the other hand, the Krawtchouk parameter p = 0.5 gives the highest EC. This is justified by the fact that p = 0.5 provided an IRKT histogram more concentrated on the interval [−T, T) with significantly high peaks in this interval. And as the capacity of our method is dependent on the number of IRKT coefficients in the interval [−T, T), we can embed more additional data with the parameter p = 0.5 and therefore a high embedding capacity. In the following tests, the Krawtchouk polynomial parameter p = 0.5 is used for calculating IRKT and inverse IRKT.
In the third test, we analyzed the influence of threshold T variation on IRKT-based RDH performance. The IRKT was applied to each test image, then the additional data were embedded in the IRKT domain with threshold T values ranging from 1 to 30. Figure 5 presents the average PSNR and EC of the stego images under these different threshold values. Figure 5 illustrates that as the threshold value T increases, stego images can carry more additional data, because a higher T value reserves more space in the IRKT domain for embedding the additional data bits. However, raising the threshold T results in a reduction in the PSNR value for the stego images. This is because the IRKT coefficients undergo a T-position shift, leading to more significant modifications when using higher T values. The figure also demonstrates that there is a notable variation in PSNR and EC for lower threshold values, but this variation gradually diminishes as the thresholds increase. Eventually, at a given threshold, the variation in PSNR and EC becomes nearly negligible. For example, when increasing the threshold from T = 1 to T = 5 for LIDC-IDRI images, the average PSNR and EC values changed from (93.06 dB, 0.2061 bpv) to (83.01 dB, 0.4008 bpv), respectively, resulting in a variation of 10.05 dB and 0.1947 bpv, respectively. Conversely, when increasing the threshold from T = 25 to T = 30, the average PSNR and EC values changed from (69.90 dB, 0.6037 bpv) to (68.39 dB, 0.6250 bpv), respectively, resulting in a variation of 1.51 dB and 0.0213 bpv, respectively.
Table 5, Table 6, Table 7 and Table 8 show the average results (PSNR and EC) over thresholds ranging from T = 1 to T = 30, for various stego medical images. These tables clearly show that the IRKT-based 3D HDR algorithm offers a high data embedding capacity while preserving a high PSNR.
In the fourth test, the robustness of the proposed IRKT-based RDH algorithm against statistical attacks was tested. Here, robustness against statistical attacks refers to the ability of the method to resist attempts at statistical analysis or detection of additional data that are hidden inside an image. Figure 6 shows a histogram-based comparison of many different types of medical images and their stego versions. We can clearly observe from this figure the great similarity between the two histograms, which is an indication of the robustness of our method against the statistical attacks.
In order to quantitatively measure the difference between the histogram distributions of the two images, we use the Bhattacharyya distance [59]. It quantifies the statistical distance or divergence between the distributions, indicating how much they overlap or differ from each other. The Bhattacharyya distance ranges from 0 to infinity. A smaller distance indicates a higher similarity between the histograms, while a larger distance indicates greater dissimilarity. Note that the Bhattacharyya distance is defined between two probability distributions. Therefore, it is important to normalize the histograms to ensure that they represent valid probability distributions before computing the distance. Table 9 presents the average Bhattacharyya distance between the histogram distributions of the original and stego images from four medical databases, using thresholds ranging from T = 1 to T = 30. The small Bhattacharyya distances shown in Table 9 clearly indicate the high similarity between the histograms of the stego images and the original images. This similarity implies that the statistical properties of the image remain largely unchanged after the additional data have been embedded, making it difficult for an observer or a system to detect the presence of hidden information through statistical analysis alone.
The fifth conducted test shows that the performance of the IRKT-based 3D RDH is improved by the cuboid size used when calculating IRKT. Figure 7 compares the PSNR and EC values of four medical databases, using cuboids of size 2 × 2 × 2, 4 × 4 × 4, and 8 × 8 × 8, with thresholds ranging from T = 1 to T = 30. From this figure, it is observed that the cuboid of size 4 × 4 × 4 offers better embedding capacity for all thresholds, followed by the cuboid of size 8 × 8 × 8 due to their more concentrated IRKT histograms in the interval [−T, T) and the presence of significantly higher peaks in this interval. On the other hand, the cuboid of size 2 × 2 × 2 shows a low embedding capacity for all thresholds due to its less-concentrated IRKT histogram in the interval [−T, T).
In terms of 3D stego image quality, the cuboid of size 4 × 4 × 4 offers a relatively high PSNR compared with the cuboid of size 8 × 8 × 8, especially for lower thresholds, while maintaining a comparable PSNR for higher thresholds. In contrast, the cuboid of size 2 × 2 × 2 has the highest PSNR due to its more limited embedding capacity. Overall, the cuboid of size 4 × 4 × 4 offers the maximum embedding capacity and a satisfactory PSNR for all thresholds.
The sixth test was performed to showcase the losslessness and reversibility of the IRKT-based 3D RDH algorithm. The results are presented in Table 10, which displays the average MSE and PSNR between the carrier 3D images and the restored images after the extraction of the additional data. Notably, the table demonstrates that the restored 3D images are exactly the same as the original carrier images. This significant advantage stems from the reversible nature of the proposed IRKT, ensuring the exact recovery of the original 3D carrier image upon extraction of the additional data. It is important to note that this achievement cannot be achieved using the classical Krawtchouk transform due to its nonreversible nature.
To show the superiority of the proposed IRKT-based RDH algorithm, we compare it with other algorithms from the perspective of embedding capacity and stego image quality. Table 11 presents a comparison between the proposed algorithm and other existing 3D algorithms based on histogram modification reported in [29,30,31]. It can be seen that the proposed algorithm is not only suitable for medical images, but also has a much better average embedding capacity than the existing 3D methods selected for this comparison.
Table 12 presents a comparison between the proposed algorithm and a similar lossless algorithm proposed by Huang et al. [22] for high-quality medical images. From this table, our method seems to have slightly better image quality based on the provided average PSNR values compared to Huang’s method. Moreover, our method achieves a higher average embedding capacity compared to Huang’s method, indicating that our method can potentially store more data within the same unit of space.
Table 13 provides a comprehensive comparison between the IRKT histogram-modification-based method and several other methods across different domains. The domains considered in the analysis include the spatial domain [32,56], the standard DCT domain [57], the standard DWT domain [58], the integer DCT domain [10], and the integer wavelet transform domains, such as integer CDF transform [8], integer 2/6 transform [9], integer Haar transform [9], and integer 9/7-F transform [9]. The comparison was conducted using four well-known 2D images: Lena, Mandrill, Airplane, and Boats. The average EC and PSNR values for each method are presented in the last columns. The results in Table 13 demonstrate that our proposed algorithm (2D IRKT, histogram modification, threshold-based embedding technique) achieves a very high embedding capacity and good stego image quality. The EC and PSNR values are higher than almost all other published algorithms selected for comparison, indicating a preference for our algorithm. Additionally, our algorithm’s average stego image quality is relatively lower than that of [32,56], integer Haar [9], and integer 2/6 [9], but it is relatively high compared to the other selected algorithms.

6. Conclusions

In this paper, we have successfully addressed the limitations of the conventional floating-point Krawtchouk transform by introducing the IRKT. By generating integer-valued coefficients in the Krawtchouk domain, the IRKT aligns effectively with the exact representation commonly utilized in lossless image applications. Expanding upon the IRKT framework, we have introduced a novel 3D RDH algorithm designed for the secure storage and transmission of extensive medical data in the IoMT sector. Our proposed algorithm employs histogram modification within the IRKT domain and a threshold-based embedding technique to embed additional/secret data. The histogram modification technique within the IRKT domain optimizes the utilization of the available embedding space, while the threshold-based embedding technique adjusts the embedding capacity and quality of the resulting 3D stego image based on a selected integer threshold.
Through comprehensive experimental assessments, we have demonstrated the effectiveness, imperceptibility, and robustness of our algorithm against statistical attacks. The algorithm showcases a substantial embedding capacity and upholds the quality of the host images. By integrating our proposed algorithm into the IoMT sector, we contribute enhanced security measures that guarantee the safe storage and transmission of extensive medical data.
Looking ahead, there are several exciting research avenues that warrant exploration. One particularly intriguing direction involves applying the IRKT in lossless compression, thereby further optimizing data storage and transmission efficiency. Furthermore, the incorporation of the IRKT into medical image encryption holds significant potential for reinforcing data privacy and security. Delving deeper into areas such as RDH within encrypted domains, stereo image RDH, and other relevant research domains could potentially yield valuable insights and drive notable advancements. Concurrently, we are dedicated to reducing the execution times of IRKT-based algorithms, particularly in scenarios where real-time or near-real-time processing is crucial. To achieve this, our strategy involves harnessing the capabilities of parallel computing techniques and leveraging the progress achieved in hardware acceleration methods [51,60], capitalizing on the advancements in these areas.

Author Contributions

Conceptualization, M.Y. and A.D.; methodology, M.Y. and A.D.; software, M.Y. and A.A.A.E.-L.; formal analysis, H.M., P.P. and A.A.A.E.-L.; investigation, O.A. and A.A.A.E.-L.; resources, H.M. and O.A.; supervision, P.P. and A.A.A.E.-L. writing—original draft, M.Y. and A.D.; writing—review and editing, P.P. and A.A.A.E.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Researchers Supporting Project Number (RSP2023R102) King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon reasonable request.

Acknowledgments

This work was funded by the Researchers Supporting Project Number (RSP2023R102) King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Calderbank, A.R.; Daubechies, I.; Sweldens, W.; Yeo, B.-L. Wavelet transforms that map integers to integers. Appl. Comput. Harmon. Anal. 1998, 5, 332–369. [Google Scholar] [CrossRef]
  2. Chen, Y.; Hao, P. Integer reversible transformation to make JPEG lossless. In Proceedings of the 7th International Conference on Signal Processing, 2004, Beijing, China, 31 August–4 September 2004; Volume 1, pp. 835–838. [Google Scholar]
  3. Gormish, M.J.; Schwartz, E.L.; Keith, A.F.; Boliek, M.P.; Zandi, A. Lossless and nearly lossless compression for high-quality images. In Proceedings of the Very High Resolution and Quality Imaging II, San Jose, CA, USA, 8–14 February 1997; Volume 3025, pp. 62–70. [Google Scholar]
  4. Said, A.; Pearlman, W.A. An image multiresolution representation for lossless and lossy compression. IEEE Trans. Image Process. 1996, 5, 1303–1310. [Google Scholar] [CrossRef] [PubMed]
  5. Zandi, A.; Allen, J.D.; Schwartz, E.L.; Boliek, M. CREW: Compression with reversible embedded wavelets. In Proceedings of the DCC’95 Data Compression Conference, Snowbird, UT, USA, 28–30 March 1995; pp. 212–221. [Google Scholar]
  6. Xiao, B.; Lu, G.; Zhang, Y.; Li, W.; Wang, G. Lossless image compression based on integer Discrete Tchebichef Transform. Neurocomputing 2016, 214, 587–593. [Google Scholar] [CrossRef]
  7. Xuan, G.; Zhu, J.; Chen, J.; Shi, Y.Q.; Ni, Z.; Su, W. Distortionless data hiding based on integer wavelet transform. Electron. Lett. 2002, 38, 1646–1648. [Google Scholar] [CrossRef]
  8. Xuan, G.; Yang, C.; Zhen, Y.; Shi, Y.Q.; Ni, Z. Reversible data hiding using integer wavelet transform and companding technique. In IWDW 2004: Digital Watermarking: Proceedings of the 3rd International Workshop, Seoul, Republic of Korea, 30 October–1 November 2004; Revised Selected Papers 3; Springer: Berlin/Heidelberg, Germany, 2005; pp. 115–124. [Google Scholar]
  9. Shaik, A.; Thanikaiselvan, V. Comparative analysis of integer wavelet transforms in reversible data hiding using threshold based histogram modification. J. King Saud Univ.-Comput. Inf. Sci. 2021, 33, 878–889. [Google Scholar] [CrossRef]
  10. Lin, Y.-K. High capacity reversible data hiding scheme based upon discrete cosine transformation. J. Syst. Softw. 2012, 85, 2395–2404. [Google Scholar] [CrossRef]
  11. Pei, S.-C.; Ding, J.-J. Integer discrete Fourier transform and its extension to integer trigonometric transforms. In Proceedings of the 2000 IEEE International Symposium on Circuits and Systems (ISCAS), Geneva, Switzerland, 28–31 May 2000; Volume 5, pp. 513–516. [Google Scholar]
  12. Oraintara, S.; Chen, Y.-J.; Nguyen, T.Q. Integer fast Fourier transform. IEEE Trans. Signal Process. 2002, 50, 607–618. [Google Scholar] [CrossRef]
  13. Yap, P.-T.; Paramesran, R.; Ong, S.-H. Image analysis by Krawtchouk moments. IEEE Trans. Image Process. 2003, 12, 1367–1377. [Google Scholar]
  14. Yap, P.T.; Raveendran, P.; Ong, S.H. Krawtchouk moments as a new set of discrete orthogonal moments for image reconstruction. In Proceedings of the 2002 International Joint Conference on Neural Networks, IJCNN’02 (Cat. No.02CH37290), Honolulu, HI, USA, 12–17 May 2002; Volume 1, pp. 908–912. [Google Scholar]
  15. Venkataramana, A.; Raj, P.A. Image Watermarking Using Krawtchouk Moments. In Proceedings of the 2007 International Conference on Computing: Theory and Applications (ICCTA’07), Kolkata, India, 5–7 March 2007; pp. 676–680. [Google Scholar]
  16. Papakostas, G.A.; Tsougenis, E.D.; Koulouriotis, D.E. Near optimum local image watermarking using Krawtchouk moments. In Proceedings of the 2010 IEEE International Conference on Imaging Systems and Techniques, Thessaloniki, Greece, 1–2 July 2010; pp. 464–467. [Google Scholar]
  17. Papakostas, G.A.; Tsougenis, E.D.; Koulouriotis, D.E. Moment-based local image watermarking via genetic optimization. Appl. Math. Comput. 2014, 227, 222–236. [Google Scholar] [CrossRef]
  18. Yamni, M.; Daoui, A.; Karmouni, H.; Elmalih, S.; Ben-fares, A.; Sayyouri, M.; Qjidaa, H.; Maaroufi, M.; Alami, B.; Jamil, M.O. Copyright protection of multiple CT images using Octonion Krawtchouk moments and grey Wolf optimizer. J. Frankl. Inst. 2023, 360, 4719–4752. [Google Scholar] [CrossRef]
  19. Yamni, M.; Karmouni, H.; Daoui, A.; Sayyouri, M.; Qjidaa, H. Blind image zero-watermarking algorithm based on radial krawtchouk moments and chaotic system. In Proceedings of the 2020 International Conference on Intelligent Systems and Computer Vision (ISCV), Fez, Morocco, 9–11 June 2020; pp. 1–7. [Google Scholar]
  20. Atakishiyev, N.M.; Wolf, K.B. Fractional Fourier–Kravchuk transform. JOSA A 1997, 14, 1467–1477. [Google Scholar] [CrossRef]
  21. Hao, P.; Shi, Q. Matrix factorizations for reversible integer mapping. IEEE Trans. Signal Process. 2001, 49, 2314–2324. [Google Scholar]
  22. Huang, L.-C.; Tseng, L.-Y.; Hwang, M.-S. A reversible data hiding method by histogram shifting in high quality medical images. J. Syst. Softw. 2013, 86, 716–727. [Google Scholar] [CrossRef]
  23. Kumar, C.V.; Natarajan, V.; Bhogadi, D. High capacity reversible data hiding based on histogram shifting for medical images. In Proceedings of the 2013 International Conference on Communication and Signal Processing, Melmaruvathur, India, 3–5 April 2013; pp. 730–733. [Google Scholar]
  24. Girdhar, A.; Kumar, V. Comprehensive survey of 3D image steganography techniques. IET Image Process. 2018, 12, 1–10. [Google Scholar] [CrossRef]
  25. Li, M.-T.; Huang, N.-C.; Wang, C.-M. A novel high capacity 3D steganographic algorithm. Int. J. Innov. Comput. Inf. Control 2011, 7, 1055–1074. [Google Scholar]
  26. Chao, M.-W.; Lin, C.; Yu, C.-W.; Lee, T.-Y. A high capacity 3D steganography algorithm. IEEE Trans. Vis. Comput. Graph. 2008, 15, 274–284. [Google Scholar] [CrossRef]
  27. Tsai, Y.-Y. An adaptive steganographic algorithm for 3D polygonal models using vertex decimation. Multimed. Tools Appl. 2014, 69, 859–876. [Google Scholar] [CrossRef]
  28. Anish, K.; Arpita, N.; Nikhil, H.; Sumant, K.; Bhagya, S.; Desai, S.D. Intelligence system security based on 3D image. In Proceedings of the 5th International Conference on Frontiers in Intelligent Computing: Theory and Applications, Bhubaneswar, India, 16–17 September 2016; Advances in Intelligent Systems and Computing; 2016. Springer: Singapore, 2017; Volume 1, pp. 159–167. [Google Scholar]
  29. Chuang, C.-H.; Cheng, C.-W.; Yen, Z.-Y. Reversible Data Hiding with Affine Invariance for 3D Models. In Proceedings of the IET International Conference on Frontier Computing. Theory, Technologies and Applications, Taichung, Taiwan, 4–6 August 2010. [Google Scholar]
  30. Jhou, C.-Y.; Pan, J.-S.; Chou, D. Reversible data hiding base on histogram shift for 3D vertex. In Proceedings of the Third International Conference on Intelligent Information Hiding and Multimedia Signal Processing (IIH-MSP 2007), Kaohsiung, Taiwan, 26–28 November 2007; Volume 1, pp. 365–370. [Google Scholar]
  31. Huang, Y.-H.; Tsai, Y.-Y. A reversible data hiding scheme for 3D polygonal models based on histogram shifting with high embedding capacity. 3D Res. 2015, 6, 20. [Google Scholar] [CrossRef]
  32. Ni, Z.; Shi, Y.-Q.; Ansari, N.; Su, W. Reversible data hiding. IEEE Trans. Circuits Syst. Video Technol. 2006, 16, 354–362. [Google Scholar]
  33. Mademlis, A.; Axenopoulos, A.; Daras, P.; Tzovaras, D.; Strintzis, M.G. 3D content-based search based on 3D Krawtchouk moments. In Proceedings of the Third International Symposium on 3D Data Processing, Visualization, and Transmission (3DPVT’06), Chapel Hill, NC, USA, 14–16 June 2006; pp. 743–749. [Google Scholar]
  34. Kumar, N.; Kumar, R.; Caldelli, R. Local moment driven PVO based reversible data hiding. IEEE Signal Process. Lett. 2021, 28, 1335–1339. [Google Scholar] [CrossRef]
  35. Wang, X.; Wang, X.; Ma, B.; Li, Q.; Wang, C.; Shi, Y. High-performance reversible data hiding based on ridge regression prediction algorithm. Signal Process. 2023, 204, 108818. [Google Scholar] [CrossRef]
  36. Singh, Y.V.; Khan, S.; Shukla, S.K.; Jung, K.-H. Local-Moment-Driven Robust Reversible Data Hiding. Appl. Sci. 2022, 12, 11826. [Google Scholar] [CrossRef]
  37. Shi, Y.-Q.; Li, X.; Zhang, X.; Wu, H.-T.; Ma, B. Reversible data hiding: Advances in the past two decades. IEEE Access 2016, 4, 3210–3237. [Google Scholar] [CrossRef]
  38. Mandal, P.C.; Mukherjee, I.; Chatterji, B.N. High capacity reversible and secured data hiding in images using interpolation and difference expansion technique. Multimed. Tools Appl. 2021, 80, 3623–3644. [Google Scholar] [CrossRef]
  39. Lin, J.-Y.; Horng, J.-H.; Chang, C.-C.; Li, Y.-H. Asymmetric Orientation Combination for Reversible and Authenticable Data Hiding of Dual Stego-images. Symmetry 2022, 14, 819. [Google Scholar] [CrossRef]
  40. Mehbodniya, A.; Douraki, B.K.; Webber, J.L.; Alkhazaleh, H.A.; Elbasi, E.; Dameshghi, M.; Abu Zitar, R.; Abualigah, L. Multilayer reversible data hiding based on the difference expansion method using multilevel thresholding of host images based on the slime mould algorithm. Processes 2022, 10, 858. [Google Scholar] [CrossRef]
  41. Ghebleh, M.; Kanso, A. A robust chaotic algorithm for digital image steganography. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 1898–1907. [Google Scholar] [CrossRef]
  42. Zhang, Y.; Luo, X.; Guo, Y.; Qin, C.; Liu, F. Zernike moment-based spatial image steganography resisting scaling attack and statistic detection. IEEE Access 2019, 7, 24282–24289. [Google Scholar] [CrossRef]
  43. Li, Q.; Wang, X.; Ma, B.; Wang, X.; Wang, C.; Xia, Z.; Shi, Y. Image steganography based on style transfer and quaternion exponent moments. Appl. Soft Comput. 2021, 110, 107618. [Google Scholar] [CrossRef]
  44. Soria-Lorente, A.; Berres, S.; Díaz-Nuñez, Y.; Avila-Domenech, E. Hiding data inside images using orthogonal moments. J. Inf. Secur. Appl. 2022, 67, 103192. [Google Scholar] [CrossRef]
  45. Radeaf, H.S.; Mahmmod, B.M.; Abdulhussain, S.H.; Al-Jumaeily, D. A steganography based on orthogonal moments. In Proceedings of the International Conference on Information and Communication Technology, Baghdad, Iraq, 15–16 April 2019; pp. 147–153. [Google Scholar]
  46. Tahiri, M.A.; Bencherqui, A.; Karmouni, H.; Amakdouf, H.; Mirjalili, S.; Motahhir, S.; Abouhawwash, M.; Askar, S.S.; Sayyouri, M.; Qjidaa, H. Implementation of a Steganography System Based on Hybrid Square Quaternion Moment Compression in IoMT. J. King Saud Univ.-Comput. Inf. Sci. 2023, 35, 101604. [Google Scholar] [CrossRef]
  47. Arsalan, M.; Malik, S.A.; Khan, A. Intelligent reversible watermarking in integer wavelet domain for medical images. J. Syst. Softw. 2012, 85, 883–894. [Google Scholar] [CrossRef]
  48. Jinna, S.K.; Ganesan, L. Reversible image data hiding using lifting wavelet transform and histogram shifting. arXiv 2010, arXiv:1004.1791. [Google Scholar]
  49. Muhammad, N.; Bibi, N.; Mahmood, Z.; Akram, T.; Naqvi, S.R. Reversible integer wavelet transform for blind image hiding method. PLoS ONE 2017, 12, e0176979. [Google Scholar] [CrossRef] [PubMed]
  50. Huang, H.-C.; Fang, W.-C. Integrity preservation and privacy protection for medical images with histogram-based reversible data hiding. In Proceedings of the 2011 IEEE/NIH Life Science Systems and Applications Workshop (LiSSA), Bethesda, MD, USA, 7–8 April 2011; pp. 108–111. [Google Scholar]
  51. Daoui, A.; Yamni, M.; Karmouni, H.; Sayyouri, M.; Qjidaa, H.; Motahhir, S.; Jamil, O.; El-Shafai, W.; Algarni, A.D.; Soliman, N.F. Efficient Biomedical Signal Security Algorithm for Smart Internet of Medical Things (IoMTs) Applications. Electronics 2022, 11, 3867. [Google Scholar] [CrossRef]
  52. Mirsky, Y.; Mahler, T.; Shelef, I.; Elovici, Y. CT-GAN: Malicious Tampering of 3D Medical Imagery using Deep Learning. In Proceedings of the USENIX Security Symposium, Santa Clara, CA, USA, 14–16 August 2019; Volume 2019. [Google Scholar]
  53. Grössinger, D.; Fischmeister, F.P.S.; Witte, M.; Koschutnig, K.; Ninaus, M.; Neuper, C.; Kober, S.E.; Wood, G. The role of superstition of cognitive control during neurofeedback training. bioRxiv 2021. [Google Scholar] [CrossRef]
  54. Soler, L.; Hostettler, A.; Agnus, V.; Charnoz, A.; Fasquel, J.; Moreau, J.; Osswald, A.; Bouhadjar, M.; Marescaux, J. 3D image reconstruction for comparison of algorithm database: A patient specific anatomical and medical image database. IRCAD Strasbg. Fr. Tech Rep 2010, 1. [Google Scholar]
  55. Baid, U.; Ghodasara, S.; Mohan, S.; Bilello, M.; Calabrese, E.; Colak, E.; Farahani, K.; Kalpathy-Cramer, J.; Kitamura, F.C.; Pati, S. The RSNA-ASNR-MICCAI BraTS 2021 benchmark on brain tumor segmentation and radiogenomic classification. arXiv 2021, arXiv:2107.02314. [Google Scholar]
  56. Kuo, W.-C.; Jiang, D.-J.; Huang, Y.-C. A reversible data hiding scheme based on block division. In Proceedings of the 2008 Congress on Image and Signal Processing, Sanya, China, 27–30 May 2008; Volume 1, pp. 365–369. [Google Scholar]
  57. Chang, C.-C.; Lin, C.-C.; Tseng, C.-S.; Tai, W.-L. Reversible hiding in DCT-based compressed images. Inf. Sci. 2007, 177, 2768–2786. [Google Scholar] [CrossRef]
  58. Huang, H.-Y.; Chang, S.-H. A 9/7 wavelet-based lossless data hiding. In Proceedings of the 2011 IEEE Symposium On Computational Intelligence For Multimedia, Signal And Vision Processing, Paris, France, 11–15 April 2011; pp. 1–6. [Google Scholar]
  59. Kailath, T. The divergence and Bhattacharyya distance measures in signal selection. IEEE Trans. Commun. Technol. 1967, 15, 52–60. [Google Scholar] [CrossRef]
  60. Daoui, A.; Yamni, M.; Chelloug, S.A.; Wani, M.A.; El-Latif, A.A.A. Efficient Image Encryption Scheme Using Novel 1D Multiparametric Dynamical Tent Map and Parallel Computing. Mathematics 2023, 11, 1589. [Google Scholar] [CrossRef]
Figure 1. (a) Input block, (b) Krawtchouk transform, (c) reconstructed block, (d) absolute difference between the input block and the reconstructed block, factored by 10−14.
Figure 1. (a) Input block, (b) Krawtchouk transform, (c) reconstructed block, (d) absolute difference between the input block and the reconstructed block, factored by 10−14.
Sensors 23 07914 g001
Figure 2. (a) Input block, (b) IRKT, (c) reconstructed block, (d) absolute difference input block and the reconstructed block.
Figure 2. (a) Input block, (b) IRKT, (c) reconstructed block, (d) absolute difference input block and the reconstructed block.
Sensors 23 07914 g002
Figure 3. General flowchart of the proposed 3D RDH in IRKT domain.
Figure 3. General flowchart of the proposed 3D RDH in IRKT domain.
Sensors 23 07914 g003
Figure 4. Average PSNR and EC obtained depending on the Krawtchouk parameter p.
Figure 4. Average PSNR and EC obtained depending on the Krawtchouk parameter p.
Sensors 23 07914 g004
Figure 5. Average PSNR and EC obtained depending on the embedding threshold T.
Figure 5. Average PSNR and EC obtained depending on the embedding threshold T.
Sensors 23 07914 g005
Figure 6. Histogram-based comparison of many different types of medical images and their stego versions.
Figure 6. Histogram-based comparison of many different types of medical images and their stego versions.
Sensors 23 07914 g006aSensors 23 07914 g006b
Figure 7. Average PSNR and EC using different cuboid sizes when calculating IRKT: (a) Results based on 80 3D images from LIDC-IDRI database; (b) results based on 70 3D images from OpenNeuro database; (c) results based on 20 3D images from IRCAD database; (d) results based on 38 3D images from BraTS 2021 database.
Figure 7. Average PSNR and EC using different cuboid sizes when calculating IRKT: (a) Results based on 80 3D images from LIDC-IDRI database; (b) results based on 70 3D images from OpenNeuro database; (c) results based on 20 3D images from IRCAD database; (d) results based on 38 3D images from BraTS 2021 database.
Sensors 23 07914 g007aSensors 23 07914 g007b
Table 1. Row vectors for computing the SERMs of 4 × 4 K matrix.
Table 1. Row vectors for computing the SERMs of 4 × 4 K matrix.
s 0 T = [0.63299−1.001131.001130]
s 1 T = [0−0.259510.25951−0.61237]
s 2 T = [0.1297600.224740.43301]
s 3 T = [0.68330−0.8165000.35355]
s 4 T = [0.93265−0.57735−1.096380]
Table 2. Row vectors for computing the SERMs of 8 × 8 K matrix.
Table 2. Row vectors for computing the SERMs of 8 × 8 K matrix.
s 0 T = [0.912366−1.293027−1.738221−1.8570551.9232384.1647484.7875270]
s 1 T = [0−0.478498−1.251264−1.2362421.2708502.5201262.305815−0.522913]
s 2 T = [0.03549600.8368630.453748−1.166437−1.687982−1.5617380.423608]
s 3 T = [0.025099−0.12975700.320848−0.824795−0.785335−1.8114220.299536]
s 4 T = [0.0229130.5270460.2051630−0.045825−1.629780−1.6535950.273437]
s 5 T = [−0.0219090.1507010.1817910.44909001.5583681.581139−0.261456]
s 6 T = [0.750597−0.7099150.8486610.782428−0.52262302.956796−0.286411]
s 7 T = [0.6128600.1420440.4429290.0798320.234718−0.2391460−0.233854]
s 8 T = [−1.775545−2.571369−0.0042220.503776−2.003696−0.9413286.1660020]
Table 3. Summary of 208 grayscale 3D medical images used as carrier images.
Table 3. Summary of 208 grayscale 3D medical images used as carrier images.
DatabaseNumber of 3D ImagesBody Part ExaminedModalityFormat and SizeClass and Bit Depth
LIDC-IDRI [52]80Human lungsCTDICOMSigned
160 × 160 × 16016 bit
OpenNeuro [53]70HeadMRINIfTISigned
192 × 192 × 19216 bit
IRCAD [54]20LiverCTNIfTISigned
128 × 128 × 12816 bit
BraTS 2021 [55]38BrainMRIDICOMUnsigned
96 × 96 × 9616 bit
Table 4. Average error between reconstructed 3D images and original ones using IRKT and classical Krawtchouk.
Table 4. Average error between reconstructed 3D images and original ones using IRKT and classical Krawtchouk.
3D Medical Images DatabaseMetric Proposed   IRKT   for   p ( 0 , 1 ) Classical Krawtchouk Transform for p = 0.5
80 images from LIDC-IDRI [52]PSNR (dB)341.5722
MSE02.9904 × 10−25
70 images from OpenNeuro [53]PSNR (dB)312.5929
MSE02.3640 × 10−22
20 images from IRCAD [54]PSNR (dB)343.7595
MSE01.8072 × 10−25
38 images from BraTS 2021 [55]PSNR360.6442
MSE03.7028 × 10−27
Table 5. Average results for seven signed 16-bit 3D CT images from LIDC-IDRI database for all thresholds.
Table 5. Average results for seven signed 16-bit 3D CT images from LIDC-IDRI database for all thresholds.
Image Name from the DatabaseAverage EC (bits) Average EC (bpv)Average PSNR (dB)
10032,156,8090.526675.7744
10052,125,9350.519075.7660
10092,207,0800.538875.8119
10452,191,7400.535175.7590
10672,190,8360.534975.7844
12012,037,0640.497375.7113
12192,115,3370.516475.7513
Table 6. Average results for seven signed 16-bit 3D IRMs from OpenNeuro database for all thresholds.
Table 6. Average results for seven signed 16-bit 3D IRMs from OpenNeuro database for all thresholds.
Image Name from the Database Average EC (bits)Average EC (bpv)Average PSNR (dB)
sub-012,434,1810.343975.3041
sub-022,271,3340.320975.2472
sub-082,944,9460.416175.4962
sub-092,649,7080.374475.3851
sub-102,538,4290.358675.3416
sub-112,641,6530.373275.3831
sub-122,664,4510.376475.3923
Table 7. Average results for seven signed 16-bit 3D CT images from IRCAD database for all thresholds.
Table 7. Average results for seven signed 16-bit 3D CT images from IRCAD database for all thresholds.
Image Name from the Database Average EC (bits)Average EC (bpv)Average PSNR (dB)
ircad_e011,127,5320.537675.8147
ircad_e021,216,5110.580175.9162
ircad_e031,084,8960.517375.7673
ircad_e041,222,9730.583275.8694
ircad_e051,149,0710.547975.8471
ircad_e061,138,6790.543075.8361
ircad_e071,217,3450.580575.9177
Table 8. Average results for seven signed 16-bit 3D IRMs from BraTS 2021 database for all thresholds.
Table 8. Average results for seven signed 16-bit 3D IRMs from BraTS 2021 database for all thresholds.
Image Name from the Database Average EC (bits)Average EC (bpv)Average PSNR (dB)
00114759,5840.858576.9145
00119707,0220.799176.6859
00125786,9240.889477.0169
00153784,9080.887277.0111
00161716,1770.809576.7315
00174704,2890.796076.6777
00181736,1300.832076.8149
Table 9. The difference between the histogram distributions of carrier images and the stego images using Bhattacharyya distance.
Table 9. The difference between the histogram distributions of carrier images and the stego images using Bhattacharyya distance.
IRKT-Based 3D RDH Method with80 Images from LIDC-IDRI70 Images from OpenNeuro20 Images from IRCAD38 Images from BraTS 2021
For T = 11.4416 × 10−77.8335 × 10−72.8873 × 10−52.2724 × 10−6
For all thresholds3.49 × 10−21.2647 × 10−48.87 × 10−24.9249 × 10−4
Table 10. Average MSE and PSNR of the images restored after additional data extraction.
Table 10. Average MSE and PSNR of the images restored after additional data extraction.
3D Medical Image DatabaseMSEPSNR (dB)
80 3D images from LIDC-IDRI database0
70 3D images from OpenNeuro database0
20 3D images from IRCAD database0
38 3D images from BraTS 2021 database0
Table 11. Comparison between the proposed algorithm and histogram-modification-based 3D RDH algorithms [29,30,31].
Table 11. Comparison between the proposed algorithm and histogram-modification-based 3D RDH algorithms [29,30,31].
AlgorithmEmbedding DomainType of Carrier DataSuitable for Medical ImageAverage EC
[31]Spatial3D image modelsNo0.2460 bits/vertex
[30]Spatial3D image modelsNo0.1725 bits/vertex
[29]Spatial3D image modelsNo0.2118 bits/vertex
OurIRKT3D grayscale medical imageYes0.4228 bits/voxel
Table 12. Comparison between the proposed algorithm and a similar lossless algorithm [22] for high-quality medical images.
Table 12. Comparison between the proposed algorithm and a similar lossless algorithm [22] for high-quality medical images.
MethodEmbedding DomainUnsign 16-bit Grayscale Medical ImageSign 16-bit Grayscale Medical Image
Average ECAverage PSNRAverage ECAverage PSNR
[22]Spatial0.248 bits/pixel74.660.124 bits/pixel74.77
Proposed IRKM0.803 bits/voxel76.700.422 bits/voxel75.51
Table 13. Comparison between our proposed method and other RDH methods [8,9,10,32,56,57,58] for some commonly used images.
Table 13. Comparison between our proposed method and other RDH methods [8,9,10,32,56,57,58] for some commonly used images.
MethodEmbedding DomainLenaMandrillAirplaneBoatAverage
EC (bits)PSNREC (bits)PSNREC (bits)PSNREC (bits)PSNREC (bits)PSNR
[32]Spatial546048.20542148.2016,17148.3730148.208588.348.2250
[56]Spatial883548.2542348.2023,19948.310,21748.3011,91948.25
[8]CDF transform52,42946.2326,21443.45--------------------39,32244.8400
[9]Integer Haar 40,10848.3913,10748.2157,14748.5327,78748.3334,53748.3650
[9]Integer 2/6 transform44,82748.2514,68048.0361,60448.3631,19548.2238,07748.2150
[9]Integer 9/7-F50,07045.2517,56445.4667,63345.3535,91445.2642,79545.3300
[10]Integer DCT64,34541.78723,90441.50593,95142.03950,84841.69758,26241.7570
[57]Standard DCT36,85730.3435,38926.4636,80529.9836,70029.7536,43829.1325
[58]Standard DWT39,87230.08----------39,45327.8439,13831.4139,48829.7767
ProposedIRKT80,56945.97819,98246.03063,37445.98670,26346.01658,54746.0031
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yamni, M.; Daoui, A.; Pławiak, P.; Mao, H.; Alfarraj, O.; El-Latif, A.A.A. A Novel 3D Reversible Data Hiding Scheme Based on Integer–Reversible Krawtchouk Transform for IoMT. Sensors 2023, 23, 7914. https://doi.org/10.3390/s23187914

AMA Style

Yamni M, Daoui A, Pławiak P, Mao H, Alfarraj O, El-Latif AAA. A Novel 3D Reversible Data Hiding Scheme Based on Integer–Reversible Krawtchouk Transform for IoMT. Sensors. 2023; 23(18):7914. https://doi.org/10.3390/s23187914

Chicago/Turabian Style

Yamni, Mohamed, Achraf Daoui, Paweł Pławiak, Haokun Mao, Osama Alfarraj, and Ahmed A. Abd El-Latif. 2023. "A Novel 3D Reversible Data Hiding Scheme Based on Integer–Reversible Krawtchouk Transform for IoMT" Sensors 23, no. 18: 7914. https://doi.org/10.3390/s23187914

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