Keywords
|
Discrete Wavelet Transform, Singular Value Decomposition, Watermarking, Peak signal to noise ratio,Mean Square Error, Normalized Cross Correlation |
INTRODUCTION
|
The evolution of medical information systems has been supported by advances in information technology. Easier access to internet has evolved into new challenges regarding security of sensitive medical information [3,6]. DWT has been widely used in many fields, including image processing.The high spreading of broadband networks and new developments in digital technology has made ownership protection and authentication of digital data since it makes possible to identify the author of an image by embedding secret information to the host image. Robust watermarks are designed to resist any modifications and are designed for the copyright protection. Frequency (transform) domain methods are based on using some invertible transformations like discrete cosine transform (DCT), discrete Fourier transform (DFT), discrete wavelet transform (DWT). Watermarking is a robust method for hiding information on multimedia files. Digital watermarking is a process to embed some information called watermark into different kinds of media called Cover Work. It is used to protect the copyrights and restores the ownership. In spatial domain techniques, the watermark is applied in the pixel domain.[4]. In Frequency domain the watermark is embedded into transformed coefficients of image giving more robust images against watermarking attacks [2]Digital watermark is a perceptually invisible pattern embedded in a digital image. The watermark can carry information about the owner of the image or the recipient the image itself or some additional information with the image (image caption embedding). Typically, watermarking the original image is needed for watermark extraction techniques. In oblivious techniques, the watermark can be extracted from the watermarked or attacked image without access to the original image. |
LITERATURE SURVEY
|
Giakoumi et al have proposed a health data management system using wavelet- based multiple watermarking scheme.The robustness of the method is enhanced through a form of hybrid coding, which includes repetitive embedding of watermarks. Multiple other studies have demonstrated the effectiveness of wavelets in watermarking[4, 7 ]. Use of SVD is shown to give resilience to attacks such as signal processing operations (such as lossy compression, filtering etc) and common geometric transformations (such as cropping, scaling, translation and rotation etc) [5, 9]. A reliable watermark extraction scheme is developed for the extraction of watermark from distorted image. Modification in all frequencies allows the development of a watermarking scheme that is robust to a wide range of attacks[8]. |
In this study we have considered medical images of three popular categories namely Computed Tomography (CT), Magnetic Resonance Image (MRI) [1] and Positron Emission Tomography(PET). An optimal watermarking technique based on SVD and DWT domain for grey-scale images is evolved based on the performance parameters namely peak signal to noise ratio(PSNR), Mean Square Error(MSE) and Normalized Cross-correlation(NC) . The algorithm embeds the watermark by modifying the singular values of the host image with multiple scaling factors. |
Paper is organized as follows: Section 3 discusses about the proposed methodology. Section 4 deals with the algorithms for watermarking and extraction as well as the performance parameters considered. A discussion of the experimental results is done in Section 5. Section 6 discusses the conclusion and the future work. |
PROPOSED METHODOLOGY
|
The medical images considered are CT, MRI and PET images. Initially the watermark is embedded in the image by setting an initial value of scaling factor α. The original image is reconstructed by extraction of the watermark. Optimum value of scaling factor is found by iteration of the above and tabulating the results to obtain desirable values of PSNR, MSE and NC. |
Embedding of a watermark is made by modifications of the transform coefficients, accordingly to the watermark or its spectrum. The inverse transform is applied to obtain the marked image. The image is decomposed into four sub-bands namely LL,LH,HL,HH. |
Watermarking is being done by embedding the watermark in the first and second level DWT sub-bands of the host image, followed by the application on the selected DWT sub bands. A hybrid DWT-SVD-based watermarking scheme is developed that requires less computation effort yielding better performance. A sequence Haar wavelet is of rescaled "square-shaped" functions which together form a wavelet family or basis. |
The Haar wavelet's mother wavelet function áôê(t) can be described as |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e001.gif) |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e002.gif) |
Singular Value Decomposition (SVD) |
Singular Value Decomposition is a matrix factorization technique having many applications in image processing. SVD is based on a theorem from linear algebra which says that a rectangular matrix A can be broken down into the product of three matrices - an orthogonal matrix U, a diagonal matrix S, and the transpose of an orthogonal matrix V . The theorem is presented as: |
A = USVT where UTU = I, VTV = I; the columns of U are orthonormal eigenvectors of AAT , the columns of V are orthonormal eigenvectors of ATA, and S is a diagonal matrix containing the square roots of eigen values from U or V in descending order. |
An illustration: |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e003.gif) |
The diagonal elements in S are the singular values of A. |
The SVD of an image is computed to obtain two orthogonal matrices U, V and a diagonal matrix the watermark W is added to the matrix S. Then, a new SVD process is performed on the new matrix S+kW to get Uw, Sw and Vw, where k is the scale factor that controls the strength of the watermark embedded to the original image. Then, the watermarked image Fw is obtained by multiplying the matrices U, Sw, and V T . The steps of watermark embedding are summarized as follows: |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e004.gif) |
PROPOSED ALGORITHM
|
The medical images considered are CT, MRI and PET images. Initially the watermark is embedded in the image by setting an initial value of scaling factor α. The original image is reconstructed by extraction of the watermark. Optimum value of scaling factor is found by iteration of the above and tabulating the results to obtain desirable values of PSNR, MSE and NC. |
Algorithm to embed watermark into cover image:
|
1)Read the cover image & watermark image. |
2)Apply DWT to cover image to obtain approximation, horizontal, vertical, diagonal DWT coefficients i.e. ca1, ch1, cv1, cd1 (LL, HL, LH, HH) |
3)Calculate the approximate DWT coefficient by adding the watermark image as in equation Ca1 (i,j)=ca1(i,j)+(α*watermark) (14) |
Where Ca1 & ca1 are the modified & original approximation coefficients and α is a scaling value as set to10. |
4)Apply SVD to the decomposed sub-bands(LL) |
5)Find Inverse DWT |
6)Increment α, apply SVD and inverse DWT. |
Algorithm for Watermark Extraction:
|
1) Apply twoâÃâ¬ÃÂdimensional DWT, to obtain the first level decomposition of the watermarked cover image i.e. c1a1, c1h1, c1v1, c1d1 |
2) Extract the watermark from vertical & approximation DWT coefficient as per the equation WN= (SN-S)/ α; where α=10 |
3)Apply twoâÃâ¬ÃÂdimensional DWT, to obtain the first level decomposition of the watermarked image. LL1, HL1, LH1, HH1. |
Performance parameters:
|
The peak signal-to-noise ratio (PSNR) is used to measure the quality of reconstructed images. PSNR is the ratio between the maximum possible power of a signal and the power of corrupting noise that affects the fidelity of its representation. Because many signals have a very wide dynamic range, PSNR is usually expressed in terms of the logarithmic decibel scale. PSNR is most easily defined via the mean square error (MSE). Given a noisefree m×n monochrome image I and its noisy approximation K |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e005.gif) |
![equation](https://www.rroij.com/articles-images/IJAREEIE-120-e006.gif) |
Where f* denotes the complex conjugates of ‘f’ and’áôæ’ is the time lag. |
RESULTS
|
Fig 2. shows the sample original, watermarked and reconstructed images of the 3 categories of medical images considered. Then peak signal to noise ratio(PSNR), mean square error(MSE), and normalized cross-correlation (NC) of the extracted image are calculated. |
Then peak signal to noise ratio(PSNR), mean square error(MSE), and normalized cross-correlation (NC) of the extracted image for MRI image is shown in Table 1 peak signal to noise ratio(PSNR), mean square error(MSE), and normalized cross-correlation (NC) of the extracted image for CT images is shown in Table 2 and Table 3 shows the PSNR,MSE and NC for PET images. |
The proposed work has been implemented using Matlab R2013a.The medical images considered are MRI, CT and PET obtained from the database “MRI Images”. The images have been watermarked using a text message. The extraction is done by decomposing the original image using Haar wavelets and embedding the SVD coefficients in the Approximate and Vertical components of the decomposed image. Fig 2 shows the sample original, watermarked and reconstructed images of the 3 categories of medical images considered. Then peak signal to noise ratio(PSNR), mean square error(MSE), and normalized cross-correlation (NC) of the extracted image are calculated and tabulated as shown in Tables 1, 2 and 3. On experimentation the scaling factors below 8 and above 14 were found to give unrealistic PSNR values of 0 or infinity. Hence the scaling factor range is fixed in the range of 8 to 14. From Tables 1, 2 and 3, it can be observed that PSNR values are in the range of 62 to 80 dB while the MSE is in the range of 0.0015 to 0.048. With increasing values of scaling factor the PSNR is seen to reduce while MSE if found to increase as shown in the tables1,2,and 3.From the tabulation, optimum value for scaling factor is noted to be 10 which gives an average PSNR of 70 dB and MSE of 0.006. The results shows a good combination of PSNR(above 70dB) and MSE(below 0.0055). This justifies our selection of optimum scaling factor. |
CONCLUSION AND FUTURE WORK
|
The watermarking can be used to provide proof of the authenticity of medical images that is to say that the medical information of one patient has been issued from the right source.the watermark has been inserted with out interfering with the documents usefulness.The proposed methodology uses DWT +SVD techniques to calculate PSNR and MSE for different medical images .SVD is very aptly used with DWT. It has been tested for different images and PSNR,MSE,NC are calculated for each image. |
This work contributes in utilizing SVD generous properties that are unexploited in image processing. This paper also introduces new trends and challenges in using SVD in image processing applications.. This paper opens many tracks for future work in using SVD as an imperative tool in signal processing. |
Tables at a glance
|
|
Figures at a glance
|
![Click Here Figure](https://www.rroij.com/articles-images/IJAREEIE-120-g001.gif) |
![Click Here Figure](https://www.rroij.com/articles-images/IJAREEIE-120-g002.gif) |
Figure 1 |
Figure 2 |
|
References
|
- Nuha Omran Abokhdair, AzizahBt Abdul Manaf. “A Prediction-Based Reversible Watermarking for MRI Images.” International Journal of Computer, Information Science and Engineering Vol:7 pp 144-147,2013
- Verma Deepti, Mishra Keerti, Verma R.L. “ A Hybrid Watermarking Scheme For Medical Images” International Journal Of Engineering And Computer ScienceVolume 2 Issue 4 pp. 1250-1255, 2013
- P.Viswanathan, Dr.P.Venkata Krishna. “Text fusion watermarking in Medical image with Semi-reversible for Secure transfer and Authentication.” P r o c o f International Conference on Advances in Recent Technologies in Communication and Computing IEEE pp.585-589,2009
- Mohammad Abdullatif, Akram M. Zeki, JalelChebil. “ Wavelet Watermarking on Medical. Images” Proc of the International Conference on Artificial Intelligence in Computer Science ,Malaysia pp. 147-156, 2013
- Venkatakishore, Dr.T.Ramashri. “SVD based Robust Image Watermarking using Particle Swarm Optimization.” International Journal on Electronic & Electrical Engineering, Volume 05, pp. 1-6,2009
- Giakoumaki, S. Pavlopoulos D, Koutsouris. “Secure and efficient health data management through multiple watermarking on medical images.” Med Bio EngComput 44”,pp. 619-631, 2006
- S.Ramakrishnan, T.Gopalakrishnan, K.Balasamy. “SVD Based Robust Digital Watermarking For Still Images Using Wavelet Transform.” Computer Science Conference Proceedings in Computer Science & Information Technology (CS & IT) pp.155-157, 2011
- Manjit Thapa, Dr. Sandeep Kumar Sood, and A.P Meenakshi Sharma. “ Digital Image Watermarking Technique Based on Different Attacks” International Journal of Advanced Computer Science and Applications, Vol. 2, pp. 14-19, 2011
- Akshay Kumar Gupta and Mehul S Raval. “A robust and secure watermarking scheme based on singular values replacement.” Sadhana Vol. 37, Part 4 pp.425-440,2012
- www.mathworks.com
- www.ece.ncsu.edu/imaging/Archives/ImageDataBase
|