Medical Image Blind Integrity Verification with Krawtchouk Moments
A new blind integrity verification method for medical image is proposed in this paper. It is based on a new kind of image features, known as Krawtchouk moments, which we use to distinguish the original images from the modified ones. Basically, with our scheme, image integrity verification is accomplished by classifying images into the original and modified categories. Experiments conducted on medical images issued from different modalities verified the validity of the proposed method and demonstrated that it can be used to detect and discriminate image modifications of different types with high accuracy. We also compared the performance of our scheme with a state-of-the-art solution suggested for medical images—solution that is based on histogram statistical properties of reorganized block-based Tchebichef moments. Conducted tests proved the better behavior of our image feature set.
With the development of eServices such as eHealth, eCommerce, and eLearning, a huge amount of digital information, such as images and videos, is transmitted over the internet. However, as data sharing is facilitated and accelerated, data security needs also increase. In particular, data reliability is of major concern. It is the key of trust one can have in the data he or she received. Taking eHealth as an example, a lot of medical information and images are exchanged through the internet between health professionals . Herein, medical images already play an important role in teleradiology and telesurgery applications, for the identification of potential diseases as well as for therapy. Indeed, medical images convey many details and specific pieces of image information that can be raised up to physicians via the tools of image processing techniques . However, such frameworks are often sensitive to various kinds of threats. Particularly, data can be intercepted and modified for illegal and malevolent purposes ranging from the patient life endangerment to the wrongful accusation of health professionals or medical institutions . Therefore, data should be proven trustworthy before being exploited.
Regarding the trustworthy of medical images, two aspects are usually considered. One is the image integrity and the other is image authenticity. The image integrity verification consists in the detection (and even prevention or correction) of image degradation or alteration, while the authenticity verification allows for determination of authorship of medical images . Various strategies have been applied to verify the integrity and authenticity of images. One class of techniques is based on image signature [5–7]. To verify the integrity, one just has to compare the signature shared along with the image with the one computed by the user from the image he received. Any difference would alert the user about an integrity loss. An aspect of consideration for such an approach is where to store the signature. One strategy stands in using image file headers, as proposed in the DICOM standard Part 15 . However, commercial implementation of the DICOM standard is not always in compliance with Part 15 specification. Moreover, trustworthy verification might fail if image signature is deleted from the image DICOM file . Another approach is based on watermarking, which allows embedding the signature into the image itself by imperceptibly modifying the gray values in images. The main drawback of watermarking methods is the possible induced image degradation that can endanger the image’s diagnostic value . The last class employs blind forensics techniques. These techniques do not require any additional image information besides the image itself to detect whether it has been modified or not . Many blind forensics methods have been proposed for the detection of various modifications of natural images including noise addition , median filtering , copy-move modification [12, 13], JPEG compression . Notice that blind forensics techniques can also be used for image steganalysis purpose, that is, to detect whether an image has been modified so as to dissimulate a secret message in between spies [15–18]. These methods can thus detect subtle image modifications if correctly designed. For verifying medical image integrity and further identifying the type of modification, very few methods have been proposed for medical images. Huang et al.  recently proposed a scheme in which a set of image features that are sensitive to modification is generated before being submitted to a classifier for image integrity verification. These features are derived from the histogram statistical properties of reorganized block-based Tchebichef moments (HRBT).
It is known that Krawtchouk moments have better image reconstruction performance than Tchebichef moments [20, 21] and have already find applications in image recognition  and fractional transform domain construction . Therefore, in this work, we propose to take advantage of them so as to build a new feature set in which the features are extracted from the histogram statistical properties of reorganized block-based Krawtchouk moments (HRBK). The proposed features can be applied for medical image integrity verification and further applied for image modification classification with better accuracy than with Tchebichef moments. The rest of this paper is organized as follows. Section 2 reviews the definition and some properties of the Krawtchouk moments. The HRBK feature generation procedure is described in Section 3. Some experimental results on image integrity verification with HRBK and its comparison with HRBT are given in Section 4. Section 5 concludes the paper.
2. Krawtchouk Moments
Let be one-dimensional signal in length , the 1D Krawtchouk transform in terms of weighted Krawtchouk polynomial is defined as in :
where is the nth order weighted Krawtchouk polynomial, defined as
Here, the weighting coefficient (x;p,N-1) and normalization coefficient ρ(n;p,N-1) are defined as
is the classical Krawtchouk polynomial:
is the hypergeometric function:
is the Pochhammer symbol:
The weighted Krawtchouk polynomial satisfies the following orthogonality property:
With the above orthogonality property, the signal can be reconstructed with the following inverse transform:
For an N×N image , the two-dimensional Krawtchouk moments are defined as
and its inverse transform is given by
3. Proposed Blind Image Forensics
In order to blindly detect whether a medical image has been modified by some global image processing techniques (e.g., filtering, lossy image compression) and further identify the type of modification, two issues should be addressed: a set of image features should be designed to distinguish original images from modified ones; a classifier model needs to be built to predict the integrity of an input image. In our work, the support vector machine (SVM) is used as classifier mainly because it has been shown that, for classification problems, SVM based solutions outperform traditional neural network approaches, such as multilayer perceptron (MLP) and radial basis functions [24, 25]. The remaining issue is thus to construct an image feature set that is sensitive to the image modifications one wants to detect and identify. In this work, histogram statistical properties of reorganized block-based Krawtchouk moments are exploited to design this feature set. This approach is developed based on the fact that low order Krawtchouk moments often have better image reconstruction performance than Tchebichef moments [21, 22] do, indicating as a consequence that most of the image energy is compacted into these low order moments. As a result, it is reasonable to extend the strategy proposed in  to Krawtchouk moments with some improvements taking advantage of the specificities of Krawtchouk moments. Figure 1 shows the schematic diagram of the extraction process of our features, which we elaborate in detail in the following part of this section.
Firstly, the image is divided into n×n nonoverlapping blocks and the Krawtchouk moments of each block are computed. In each block, Krawtchouk moments are partitioned into 3L+1 (n = 2L) subbands. Then the coefficients of the same subband in each block are grouped together to generate an L-scale coefficient tree for the whole image. In our experiments, we considered n = 8. As shown in Figure 2, one moment block is thus divided into 10 subbands. Figure 3 illustrates the corresponding L-scale coefficient tree for the whole image, where the coefficients of the same subband i (i = 0, 1, …, 9) from all blocks are reorganized into the group Gi.
Contrarily to , in order to emphasize the Krawtchouk moment statistical variations under modifications, we further divide coefficients groups , , and as illustrated in Figure 4. At last, all the block-based Krawtchouk moments are reorganized into 25 groups. The detailed way of reorganizing the Krawtchouk moment block can be found in .
Once Krawtchouk moments are reorganized, the image features can be generated. The first class of features we use corresponds to the statistical moments of the discrete Fourier transform (DFT) of histogram of one Krawtchouk moment subband :
where is the DFT coefficient at frequency in the histogram of one subband of reorganized Krawtchouk moment transform coefficients. As it can be seen, features defined in (12)-(14) act as high-pass filters for the histogram. In order to take advantage of the rest of frequency information of the histogram, the second class of features obtained is defined as 
In addition to these two classes of features generated from the original image, features are also generated from the prediction error image. The basic idea is to achieve a second set of features that is more image content independent. Let be an image pixel at position . In this work, its predicted value is defined as 
Afterwards, the prediction error image is constructed by subtracting the predicted image from the original one, that is, .
Finally, the feature set are extracted from the 25 groups of the reorganized Krawtchouk moments of the original image, giving access to a first set of 150 features. Regarding the prediction error image, the twelve subbands, ~, ~, ~, are not exploited because the coefficients in these subbands are of very small values and are insignificant for image description purpose, so only 78 image features are extracted from the prediction error image. To conclude, one image will be represented or summarized by a feature vector of 228 components which are used to train the SVM classifier for discriminating the original images from the modified ones, as well as to identify the kind of modification (e.g., filtering, lossy compression).
4. Experimental Results and Discussion
In order to verify the validity of the proposed image feature set for medical image blind verification, two medical image datasets in different imaging modalities were used. The performance of our features is compared with performance of those features proposed in  which are based on Tchebichef moments. We detail and discuss some experimental results in this section.
4.1. Test Data, Modifications, and Test Parameterization
Our image datasets consist of medical images issued from two modalities: 100 computed tomography (CT) images of size 512×512, 12 bits encoded; 100 magnetic resonance (MR) images of size 181×181, encoded onto 12 bits [28, 29]. Some samples of these two test datasets are illustrated in Figure 5. Notice that the features are generated from the 128×128 image block centered in each image. To make a fair comparison with the method proposed in , we used the same types of image modifications. That is, seven types of common image modifications were considered: JPEG2000 compression, Gaussian filtering, Laplacian filtering, brightening, scaling up, histogram equalization, and JPEG compression. For each type of modifications, five different modification intensities were considered. The detailed parameters that correspond to each modification type are listed in Table 1.
In the following experiments, the 100 original images of each modality along with their modified versions were divided into two groups. One group was used as the training set of SVM classifier and the other used for SVM testing. Training and testing sets were randomly generated. It is important to notice that the following experimental results are given in average after ten rounds of training and testing.
4.2. Detection of Image Modification
In this experiment, the objective was to discriminate the original images from the modified ones using the proposed HRBK features. In the first test, each modification type described in Section 4.1 was treated separately. For each modification type, fifty original test images and their modified ones, that is to say, three hundred images (50 originals and 250 modified images corresponding to 5 different intensities; see Table 1), were used to train an SVM classifier. The rest of the original images and their corresponding modified images were then exploited as the test set. The detection rates for each type of modifications are provided in Table 2. As it can be seen, the proposed scheme was able to detect all types of modifications with higher accuracy than the cases with HRBT features.
In the second test, all types of modifications were considered together, the objective was to discriminate the original images from the modified ones. Again, fifty original test images and all their modified versions were used so as to generate the SVM model for classification. This training set was thus composed of the features of 1800 images. Then the rest of the original images and their corresponding modified images were fed to the trained SVM model to determine whether an input image was modified or not. The modification detection rates for this setup are recorded in the last row of Table 2. As it can be seen from this table, HRBK detection rates are better than HRBT ones.
The third test was conducted to compare the performance of the two feature sets under different modification parametrization settings. In this test, for each modification type and parameter, the training dataset was composed of fifty original images and their modified ones, which are modified according to the corresponding modification type and parameter; the test set included the remaining fifty original images and their modified ones. The training and testing of SVM model were similarly done as in aforementioned tests. The results of this test are recorded in Tables 3 and 4 for the two image datasets. As it can be seen from these two tables, HRBK outperforms HRBT in most cases.
4.3. Identification of the Image Modification Type
In this second experiment, the objective was not only to detect whether an image had been modified, but also to determine the type of the image tampering. To do so, a multiclass SVM was constructed based on one-versus-one binary classifier, and pairwise coupling  was employed to combine results from all binary classifiers. As previously, fifty original images along with all their modified versions, that is to say, 1800 images, were used as training set. The rest of the images, original and modified, were used as the test set for the evaluation of the trained classifier. Detection rates are provided in Table 5 for different modification types considering CT images and MR images. As it can be seen, for each modification type, the detection rates of HRBK features are higher than those of HRBT features in most situations. This is due to the better image representation of the Krawtchouk moments than that of Tchebichef moments. It can be observed that HRBK achieved much better results for CT images with JPEG compression than HRBT did, but worse results for MR images with JPEG compression. This could be due to the fact that MR images have more image details than CT images do. After JPEG compression, more image discriminative information is lost for MR images, which leads to lower modification detection rates.
4.4. Evaluation of the Influence of Block Size on Integrity Verification Performance
In this last experiment, the objective was to determine the influence of the centered block size (size of the image block used to generate features) on detection rates for the different modification types described in Section 4.1. Let us recall that the number of Krawtchouk moments in a subband depends on the block size. In this test, three different block sizes were considered; i.e., 128×128, 64×64, and 32×32 were chosen to generate the image features while considering the same ten subbands (see Section 3). As before, for each modification type, fifty original images and their modified images were used as training set and the others as test set, and multiclass SVM was used as classifier. The plots in Figures 6 and 7 show the relations between detection rates and block sizes for CT and MR images, respectively. As it can be seen, detection rates decrease as the block sizes reduce. However, the HRBK method has better detection rates in most cases, and this superiority is much obvious for small block sizes. This phenomenon can be explained as the result of better description power of Krawtchouk moments compared to Tchebichef moments and their ability to keep essential discriminative image information as block size decreases.
In previous experiments, the ratio between the size of training set and that of the testing set was 1:1. The influence of the ratio value on the performance of the proposed method was also tested. In this test, we used seventy original images along with all their modified versions as the training set. The rest of the images, original and modified, were used as the test set to evaluate the performance of the trained classifier for modification classification. This leads the ratio between training set size and testing set size to be 7:3. Table 6 shows the detection rates of HRBK with 128×128 block size for different modification types. As it can be seen, the detection rates with the HRBK features increased in most cases. However, the choice of the optimal value of the ratio needs further study; this will be part of our further work.
For the purpose of verifying the integrity of medical images, a new set of image features was proposed, which was developed based on the histogram statistical properties of reorganized block-based Krawtchouk moments. Higher image modification detection rates are achieved when the proposed method is used in most modification situations. However, our method can only be used to detect modifications with “a priori knowledge.” More clearly, the possible types of modification an image may undergo are identified before the SVM training. In the case of an “unforeseen attack,” the proposed method will identify the modification as the closest or most similar a priori known modification type. Moreover, we considered only global image modification types; further study should be done to investigate the image integrity verification problem with local image modification types. To further improve the performance of our method, one way could be to use other prevalent machine learning methods, such as deep leaning, to construct the classification model.
In this paper, we have proposed a new set of image features based on Krawtchouk moments. These HRBK features are used to detect medical image modifications and also to distinguish the types of the image modifications through a classification-based strategy. Compared with the existing HRBT features, HRBK achieves better detection rates for almost all kinds of image modification types and is more robust with respect to feature extraction area size. The proposed integrity verification method relies on the ability of the image features to describe the differences between unmodified and modified images. Future works include exploiting image moment properties to construct more selective image features and speeding up the feature extraction process by taking advantage of parallelism in the computation.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
A. P. Dhawan, Medical Image Analysis, IEEE Press, Hoboken, NJ, USA, 2nd edition, 2011.View at: Publisher Site
D. Bouslimi, G. Coatrieux, M. Cozic, and C. Roux, “A telemedicine protocol based on watermarking evidence for identification of liabilities in case of litigation,” in Proceedings of the IEEE 14th International Conference on e-Health Networking, Applications and Services (Healthcom), pp. 506–509, October 2012.View at: Publisher Site | Google Scholar
H. Huang, G. Coatrieux, H. Z. Shu, L. M. Luo, and C. Roux, “Medical image integrity control and forensics based on watermarking – approximating local local modifications and identifying global image alterations,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, pp. 8062–8065, Boston, MA, USA, August 2011.View at: Publisher Site | Google Scholar
G. Coatrieux, H. Maitre, and B. Sankur, “Strict integrity control of biomedical images,” in Proceedings of SPIE, vol. 4314, pp. 229–240, 2011.View at: Google Scholar
“Part 15: Security and System Management Profiles,” in Digital Imaging and Communications in Medicine, DICOM PS 3.15, 2004.View at: Google Scholar
Q. Wang and R. Zhang, “Double JPEG compression forensics based on a convolutional neural network,” EURASIP Journal on Information Security, vol. 2016, no. 1, 2016.View at: Google Scholar
Y. Q. Shi, G. Xuan, D. Zou et al., “Image steganalysis based on moments of characteristic functions using wavelet decomposition, prediction-error image,” in Proceedings of International Conference on Multimedia Expo, p. 1, 2005.View at: Google Scholar
S. K. Jena and G. V. V. Krishmna, “Blind steganalysis: estimation of hidden message length,” Journal of Computers, Communications, and Control, vol. 2, pp. 149–158.View at: Google Scholar
P. T. Yap, R. Paramesran, and S. H. Ong, “Krawtchouk moments as a new set of discrete orthogonal moments for image reconstruction,” in International Joint Conference on Neural Networks, vol. 1, pp. 908–912, 2002.View at: Google Scholar
B. Kaur and G. Joshi, “Lower Order Krawtchouk Moment-Based Feature-Set for Hand Gesture Recognition,” Advances in Human Computer Interaction, vol. 2016, 2016.View at: Google Scholar
C. A. Cocosco, V. Kollokian, R. K.-S. Kwan, A. C. Evans, and BrainWeb, “Online interface to a 3D MRI simulated brain database,” in Proceedings of the Proceeding of 3rd International Conference on Functional Mapping of the Human Brain, vol. 5, p. 475, 1997.View at: Google Scholar