Open access peer-reviewed chapter

The DICOM Image Compression and Patient Data Integration using Run Length and Huffman Encoder

Written By

Trupti N. Baraskar and Vijay R. Mankar

Submitted: 18 June 2019 Reviewed: 12 August 2019 Published: 31 October 2019

DOI: 10.5772/intechopen.89143

From the Edited Volume

Coding Theory

Edited by Sudhakar Radhakrishnan and Muhammad Sarfraz

Chapter metrics overview

1,110 Chapter Downloads

View Full Metrics

Abstract

Maintaining human healthcare is one of the biggest challenges that most of the increasing population in Asian countries are facing today. There is an unrelenting need in our medical community to develop applications that are low on cost, with high compression, as huge number of patient’s data and images need to be transmitted over the network to be reviewed by the physicians for diagnostic purpose. This implemented work represents discrete wavelet-based threshold approach. Using this approach by applying N-level decomposition on 2D wavelet types like Biorthogonal, Haar, Daubechies, Coiflets, Symlets, Reverse Biorthogonal, and Discrete Meyer, various levels of wavelet coefficients are obtained. The lossless hybrid encoding algorithm, which combines run-length encoder and Huffman encoder, has been used for compression and decompression purpose. This work is proposed to examine the efficiency of different wavelet types and to determine the best. The objective of this research work is to improve compression ratio and compression gain.

Keywords

  • DICOM
  • discrete wavelet
  • N-level decomposition
  • threshold approach
  • data hiding algorithms

1. Introduction

Digital technology has, in the last few decades, entered in almost every aspect of medicine. There has been a huge development in noninvasive medical imaging equipment. Since there are multiple medical equipment manufacturers, there is a strong need to develop a standard for storage and exchange of medical images. DICOM (Digital Imaging and Communications in Medicine) makes medical image exchange easier and independent of the imaging equipment manufacturer. The DICOM standard has been developed by ACR-NEMA to meet the needs of manufacturers and users of medical imaging equipment for interconnection of devices on standard networks. The DICOM technology is suitable when sending images between different departments within hospitals and/or other hospitals and the consultant. DICOM file contains both a header, which include text information such as patient’s name, modality, image size, etc., and image data in the same file. Hence DICOM standards are widely used in the integration of digital imaging systems in medicine. Figure 1 shows the structure of DICOM image file. It has two main components. The first is header; it consists of 128 bytes of file preamble which is followed by string by 4-byte prefix, and it contains four-character string. The second is data set; it consists of multiple set of data elements. Each data element has four fields; these are tag, value representation, value length, and value field. The third is image pixel intensity data; it contains necessary medical image data display like number of frames, lines, columns, etc.

Figure 1.

The structure of a DICOM image file.

1.1 File format used by DICOM images

There are four major file formats in medical imaging, and they are Neuroimaging Informatics Technology Initiative (NIfTI), Analyze, DICOM, and MINC. The task of the image file format is to provide a standardized way to store the unique data in a much organized and systematic manner and showcase how the pixel data understood the correct loading, visualization, and analysis was derived by the software. The major file format currently useful in medical imaging is DICOM format. The DICOM format includes some information that can be useful for image registration, such as position and orientation of the image with respect to the data acquisition device and patient information with respect to voxel size. DICOM file format design consideration is based on the following concept such as pixel depth, photometric interpretation, metadata, and pixel data. The DICOM file format is created by addition of header size and pixel data. Mathematical equations are as follows:

DICOM File Format = Header Size + Pixel Data Size E1
Pixel Data Size = Rows × Columns × Pixel Depth × Number of Frames E2

The more popular formats used in daily practice are the JPEG, JPEG 2000, TIFF, GIF, PNG, and BMP formats. The images saved in these formats can be accessed on any personal computer without the need of specific viewers. File format are designed with the help of image conversion technique and coding schemes.

The rest of the chapter is organized as follows. Section 2 provides an overview of image file format and image standard used for medical image compression, and data hiding methods are described in the Section 3. Section 4 provides proposed work brief explanation. Section 5 discusses regarding results that are obtained after implementation of application. Finally, Section 6 concludes the chapter.

Advertisement

2. Related work

The related work is a comprehensive summery of previous research on image file format, standard image compression using transform coding, and patient information integration into image for DICOM images. The more popular formats used in daily practice are the JPEG, JPEG 2000, TIFF, GIF, PNG, and BMP formats. The images saved in these formats can be accessed on any personal computer without the need of specific viewers. File formats are designed with the help of image conversion technique and coding schemes [1, 2, 3]. Figure 2 shows the basic digital image file formats and its classification. The vector images are not commonly used in medical data processing.

Figure 2.

Basic digital classification image file format.

Table 1 gives the summary of various parameters of raster image file format [4, 5, 7, 8, 9, 10, 11]. Table 2 gives a characteristic overview of the major file formats currently used in medical imaging, i.e., NIfTI, Analyze, DICOM, and MINC [6, 12, 13].

Sr. no. File format Extension Bit depth Compression type Compression performance Name of supported free image viewer
Gray Color
1 BMP bmp 1,4,8 1,4,8,24 Lossless Very low IrfanView, XnView, Osiris, ImageJ
2 DICOM dcm 8, 16 8, 24, 48 Lossless Low IrfanView, XnView, Osiris, ImageJ
3 GIF gif 1, 4, 8 1, 4, 8 Lossless Medium IrfanView, XnView, Osiris, ImageJ
4 JPEG jpg 8 24 Lossy Average IrfanView, XnView, Osiris, ImageJ
5 JPEG 2000 Jp2 8, 16 24, 48 Lossless High IrfanView, XnView
6 PNG png 1, 4, 8, 16 1, 4, 8, 24, 48 Lossless High IrfanView, XnView
7 TIFF tiff 8, 16 8, 24, 48 Lossless or lossy Very high IrfanView, XnView

Table 1.

Summary of raster (bitmap) image file format.

Sr. no. File format Extension Data type Header Compression scheme supported
Integer Float Complex
1 DICOM .dcm Signed and unsigned (8 bit, 16 bit, 32 bit) Not supported Not supported Variable length binary format JPEG, RLE, JPEG-LS, MPEG2/MPEG4, JPEG XR
2 NIfTI .nii Signed and unsigned (8 bit to 64 bit) Signed and unsigned (32 bit to 128 bit) Signed and unsigned (64 bit to 256 bit) Fixed length (532 byte binary format) gzip (it is a software application used to store compressed and decompressed file)
3 MINC .mnc Signed and unsigned (8 bit to 32 bit) Signed and unsigned (32 bit to 64 bit) Signed and unsigned (32 bit to 64 bit) Extended binary format gzip (it is a software application used to store compressed and decompressed file)
4 Analyze .hdr and .img Signed (8 bit to 32 bit) and unsigned (8 bit) Signed (32 bit to 64 bit) Signed (64 bit) Fixed length (348 byte binary format) High dynamic range (HDR) imaging uses sub-band coding technique which is an example of lossy technique

Table 2.

Characteristics of medical image file format.

Two types of compression methods are classified. The lossless image has huge application in archival of medical and digital radiography document, where loss of information in original image could consider improper diagnosis. The medical imaging application required lossless image compression. Thus, medical image compression application development is a challenging problem. The survey paper [14] conveys that compression ratio 4:1 is possible using lossless compression. An increasing volume of data generated by new imaging modality, CT scan and MRI lossy compression technique are used to decrease the cost of storage and increase the efficiency of transmission over networks for teleradiology application [12]. There are two main categories of compression lossless (reversible) and lossy (irreversible). DICOM support lossless compression schemes like run-length encoding, Huffman coding, LZW coding, area coding, and arithmetic coding. The RLE is used for medical image compression in hybrid approach, where gray scale value gives certain interesting fact about the distribution in image. The background pixels of all the medical image are low values, and they differ by +3 or −3. The RLE is based on dynamic array implementation, no need to process whole image [8]. The paper served that the number of repeated zero count which is represented as “RUN” and appends the nonzero coefficients represented as “LEVEL” [15, 16]. The Huffman coding is a lossless compression technique, which is used for medical image compression. This works on variable length encoding principal, which includes calculation of length of unique codes. It generates a binary tree, which is also known as Huffman tree. Huffman algorithm gives higher compression ratio in the case of medical image compression. During the whole process of compression, there should not be any loss of information that will affect proper diagnosis [17, 18]. The Huffman code is designed to integrate the lowest probable symbols, and this integration is repeated until only two probabilities of two symbols are left. In this survey paper, certain improvements are discussed on the existing Huffman technique which will help to preserve any loss of information during compression that will affect proper diagnosis [19]. In lossy compression method, data are rejected during compression and cannot be recovered completely. This method reaches much greater compression performance than lossless compression. Wavelet and higher-level JPEG are the example of lossy compression technique where JPEG 2000 is a progressive lossless-to-lossy compression algorithm [20, 21, 22]. This article uses the concept of data hiding into image for data encryption. In order to enable large capacity of data hiding and maintaining good image quality, the data integration is applied on detail coefficients of high-frequency sub-bands. It works on transform domain of multilevel two-dimensional discrete wavelet transform. The objective of this implementation is to perform image compression as much as possible. It will help to reduce the redundancy of the image and to store or transmit data in an efficient form. As in telemedicine, the medical images are transmitted through advanced hyperlinks; medical image compression without any loss of useful information is of immense importance for the fast transfer of the medical data [23].

Advertisement

3. Proposed work

This proposed compression approach deals with .dcm file of DICOM format. It splits .dcm file into patient data with bmp.txt extension and gray scale image with .bmp extension. Then N-level DWT using various wavelet types is applied to a gray image. Firstly, this splits the image into n number of high-frequency sub-bands (HLn, LHn, HHn) where n = 1, 2, 3…..N and one low-frequency sub-band (LLn) where n = maximum level (N). The high-frequency sub-bands at levels 1, 2, 3, and 4 are threshold and quantized and find detail coefficients are encoded directly through run-length encoding. Secondly, the one low-frequency sub-band is also threshold and quantized and find high-level approximate coefficient. Lastly, both the coefficients (detail coefficients, high-level approximate coefficient) are encoded by Huffman coding.

In a gray scale image, each pixel is represented by 8-bit unsigned integer value. The minimum and maximum value of unsigned integer is 0 to 255. The 0 represented black and 255 represent white. In text file every text file is represented by ASCII value. The ASCII value is run between 0 and 128. The extended ASCII value is 8 bit, and it matched with 8-bit pixel intensity. So, both the entities are treated as normal integer. In this proposed system, we have practiced bitwise XOR approach for text integration into image. The proposed work deals ASCII conversion of patient data and multilevel two-dimensional discrete wavelet transform. Advantages of this work are high data integrity even with large patient data. Accepted levels of imperceptibility, excellent PSNR values, and high CR and good payload capacity are obtained. Figures 3 and 4 represented the block diagram of compression and decompression with data integration scheme. The decompression process is the inverse process of DICOM image compression as shown in block diagram in Figure 4 . The DICOM compressed image is recreated with acceptable quality through the abovementioned process.

Figure 3.

Block diagram of DICOM image compression and data integration method.

Figure 4.

Block diagram of DICOM image decompression and data and image extraction.

3.1 Algorithm for integration of patient information (text) into image file

  1. Select a proper Greyscale BMP image.

  2. Select patient text file with .txt extension

  3. Execute while loop (Number of Character Count in text file<= Number of pixels in the image)

  4. Convert Character vector and pixel value into an unsigned 16-bit integer using function str2num and unit 16.

  5. Integration of text in to Image Pixel = (Converted 16-bit unsigned pixel value) XORing_Bitwise (Converted 16-bit unsigned character value)

  6. Increment Pixel value until last value.

  7. Increment Character value until last value

  8. Check the Character Count in text file = Character Count +1

  9. End Loop

  10. Rest of Integration of text in Image Pixel = original pixels of image

3.2 Extraction of patient information (text) from the image file

  1. Open the original image and Integration of text in to Image Pixel

  2. Execute while loop (Number of Pixel Count in image file ≤ Number of pixels in the image)

  3. Convert pixel value into 16-bit integer value

  4. Integration of text in to Image Pixel = convert to 16-bit integer (Integration of text in to Image Pixel)

  5. A = (Original value of Pixel) XORing_Bitwise (Integration of text in to Image Pixel)

  6. if A = 0 then break, else extracted text file which is equal to A

  7. Extract original Pixel = Next original Pixel in image

  8. Integration of text in to Image Pixel = Next Integration of text in to Image Pixel

  9. End Loop

3.3 Multilevel 2D DWT decomposition on wavelet types

We implement an N-level 2D DWT decomposition. At each level of decomposition, the LL sub-band from the previous level is obtained, and each previous level is replaced with four new sub-bands. Each new sub-band is half the width and half the height of the LL sub-band from its parent sub-bands. The formula to calculate the total number of sub-bands depends on the number of level n. The number of sub-bands is therefore 3n + 1, where HHn represent high-frequency band, LLn is low-frequency band, and LHn and HLn are middle-frequency bands. The coefficients in LL are dominant. If any of the coefficients in LLn frequency band are changed, observer can observe that the corresponding spatial domain image has been modified.

Figure 5 shows the process of character integration in LHn sub-band, and it generates wavelet coefficients. It shows the scale and orientation selectivity of the DWT. Greatest energy is contained in the LLn sub-band, and the least energy is in the HHn sub-band. The HLn sub-band contains the vertical edges, and the LHn sub-band contains the horizontal edges. In this proposed work, we focus on data integration in LHn sub-bands because this band has high energy distribution as compared to other bands like HLn and HHn. The finest wavelet type and the appropriate coefficient selection method using threshold and quantization.

Figure 5.

The process of character integration in LHn sub-bands and generation of wavelet coefficients.

3.4 Apply thresholding and quantization technique on generated coefficients

3.4.1 Thresholding

Let n × n be matrix of an original image; noise observation can be written as s = x + n , where s = noise observation, x = original image, and n = noise. Let s (i), x(i), and n(i) denote ith sample of pixels. By applying discrete wavelet transform, the observed noised image obtained wavelet coefficients y = θ + z, where y = Ws, θ = Wx, and z = Wn, respectively. To recover θ and y, y is transformed into wavelet domain that decomposes y into many sub-bands [19, 20]. Then the coefficients with small value in the sub-bands are dominated by noise, thus replacing noise coefficients by zero. It is denoted by

y i = θ i + z i E3

If

y i ̂ = abs y i < λ E4
y i = 0 E5

where y (i) is the input and noise wavelet coefficients, λ is the Threshold Valu, y i ̂ is the Threshold output.

We define the PCDZ parameter; this parameter is required to calculate the percentage of nonzero DWT coefficients.

PCDZ = 100 NBz Ly E6

where NBz = number of zeros in DWT coefficients.

Ly = Number of Coefficients in DWT.

The proposed method used in the global threshold value that is derived by Donoho [19, 20] is given by the equation below. It is known to have a universal threshold.

λ = σ 2 log Ly E7

where Ly is the number of pixels in the medical image and σ is the noise variance.

3.4.2 DWT coefficient quantization

The quantization of each level permits to collect the set of nearest values. The uniform quantization on thresholded DWT coefficients in sub-bands will be transformed and to be contained in the interval width for quantization is between 0 to 2 Q . The quantized matrix will be computed as follows: to choose the quantization value Q (Q represents the interval width for quantization of the DWT coefficients in sub-bands), further determine the max(y(i)) and min(y(i)) values of the DWT coefficients which will represent as DWTmax and DWTmin. The uniform quantization on the resulting DWT coefficient sub-bands is formulated by the following equation:

DWTmin = min y i E8
DWTmax = max y i E9
DWTCQ = round ( 1 + 2 Q y i ̂ DWTmin DWTmax DWTmin E10
% ofIDWTCQ = round DWTmax DWTmin / round ( 1 + 2 Q DWTCQ + DWTmin ) E11

3.5 Encoding of wavelet coefficients using run-length encoder and Huffman encoder

In this implementation, hierarchical relationship of wavelet structure is explored to arrange wavelet coefficients into odd rows and even rows. The wavelet coefficients odd rows contain an ordering of wavelet coefficient that acts as approximate (smooth) value, and even rows contain different sign data that act as detail values. After evaluated many zeros in different orders of wavelets, by applying hard thresholding condition Eqs. (4), (5) on wavelet coefficients. In this whole process, separated approximate coefficients contain best information, while detail coefficients contain information like shapes and edges of image. The threshold condition chooses fixed threshold value to obtain desired quality of reconstructed image. After classifying threshold coefficients, need to transmit those coefficients using lossless method which will further be used for decompression purpose. Now encoded detail coefficients with run-length encoding excluding the highest approximate coefficient LL3 sub-bands, because LL3 sub-band does not have much long run of zeros. To convert repetitive data into bit stream, Huffman encoder has been used in this implementation. Huffman code is an example of optimum prefix code; these codes are generated using variable code length where a number of bits are essential. This will be helpful in average code length calculations, and thus the data compression is taking place where sometimes compressed image is smaller than original image.

Advertisement

4. Results and discussion

The performances of implemented method are based on few essential criteria: the obtained compression ratio (CR), compression gain, and the quality of the reconstructed image using PSNR, MSE (mean squared error), and SNR. Data compression equations are given below.

Data compression ratio = Uncompressed size/Compressed size

CR = X Y E12

Space saving (%) determines performance of transformation efficiency over storage of data bits for original bit size to unprocessed bit size. It is like compression ratio; however it reflects percentage of how much data space is saved following compression [21]. It is given by the following equation:

Percentage of Compression Gain = 1 – Compressed size/Uncompressed Size

Compression Gain = 1 Y X 100 E13

The calculated peak signal to noise ratio between maximum values is power of signal and power of distorting noise which affects the quality of its representation. The PSNR is generally expressed in terms of logarithmic decibel scale [22].

PSNR = 20 log 10 MAX f MSE E14

MAX f is the maximum signal value which exists in the original image, which is known to be good value of image?

where the MSE

MSE = 1 mn 0 m 1 0 n 1 f i j g i j 2 E15

where f is the matrix data of our original image. g is the matrix data of our degraded image. m is the numbers of rows of pixels of the images. i is the index of that row. n is the number of columns of pixels of the image. j is the index of that column.

where SNR is given as the ratio of the mean value of the signal and the standard deviation of the noise.

SNR = 20 log Intensity Signal Intensity Noise E16

To analyze the performance of our proposed method, we take Table 3 as an input MR image for evaluation of various parameters, and their information are as follows:

Image name Input size (in KB) Level of decomposition (N) Noise variance ( σ 2 ) The size of the DWT coefficient arrays (Ly) Threshold value ( λ )
1.2.840.113619.2.5.1762583153.215519.978957063.122.dcm 522 3,4 σ = 1 266,539 3.29416477

Table 3.

Input MR image for evaluation of various parameters.

Figure 6 shows that pop-up message is generated after compression of DICOM input image, and it displayed warning dialog that image is compressed successfully, and then after pressing the ok button, compressed image is stored in image folder with. Hdwt extension ( Figure 7 ).

Figure 6.

Graphical user interface of DICOM image 1.2.840.113619.2.5.1762583153.215519.978957063.122.dcm of 522 KB for N = 3 and after pressing pop-up button of image compression successfully with size.

Figure 7.

Graphical user interface of DICOM image 1.2.840.113619.2.5.1762583153.215519.978957063.122.dcm of 522 KB for N = 4 has displayed PSNR, MSE, and SNR parameters after pressing pop-up button of decompression.

Figure 8 shows that if input image size is 522 KB, then Biorthogonal DWT give highest compressed size, and Reverse Biorthogonal DWT gives lowest compressed size. The decomposition level = 4 gives better result than N = 3.

Figure 8.

The plot for size of compressed image in KB vs wavelet types for N = 3 and 4.

Figure 9 plot shows comparison of compression ratio when decomposition level is 3 or 4 for various wavelet types where size of input MRI is 522 KB. Biorthogonal DWT gives higher compression ratio, and Discrete Meyer DWT and Reverse Biorthogonal DWT give lower compression size. N = 4 gives better compression ratio (15:1 for Biorthogonal DWT and 10:1 Daubechies DWT) than N = 3.

Figure 9.

The plot for compression ratio vs wavelet types For image size of input(MRI) image size 522 KB, when N = 3 and 4 compare.

In Figure 10 , figure graph shows comparison of compression gain when decomposition level is 3 or 4 for various wavelet types where size of input MRI is 522 KB. Biorthogonal DWT gives higher compression gain, and Discrete Meyer DWT and Reverse Biorthogonal DWT give lower compression gain. N = 4 gives better compression gain (93.6781% for Biorthogonal DWT and 90.6130% for Daubechies DWT) than N = 3 ( Figure 11 ).

Figure 10.

The plot for compression gain vs wavelet types where image size of input (MRI) image size 522 KB for N = 3 and 4.

Figure 11.

The 36 KB of input .dcm file, compressed HDWT file, text file, and patient data integrated into the image file for Biorthogonal DWT.

Table 4 shows comparison of image quality parameters when application run at decomposition levels 3 and 4 for various wavelet types where size of input MRI is 522 KB. Daubechies DWT gives higher PSNR (42.0998db) where MSE is (35.2626db) and SNR is (28.3993 db). So, picture quality is good for N = 3. When we consider N = 4 decomposition level, we can choose Daubechies DWT for compression to achieve higher quality of medical image. This tool does not work on Haar DWT; there is no output for Haar wavelet type.

Initial size of image: 522 KB (image name: 1.2.840.113619.2.5.1762583153.215519.978957063.122.dcm) Image quality parameter for N = 3 (in decibels) Image quality parameter for N = 4 (in decibels)
Wavelet types PSNR MSE SNR PSNR MSE SNR
Biorthogonal DWT 35.3358 35.2626 21.6352 30.7777 35.3626 17.0771
Haar DWT
Daubechies DWT 42.0998 35.2626 28.3993 39.6729 35.3626 25.9723
Symlets DWT 39.0696 35.3626 25.3691 36.6956 35.3626 22.995
Coiflets DWT 40.4769 35.3626 26.7763 36.1026 35.3626 22.402
Reverse Biorthogonal DWT 32.0256 35.3626 18.325 26.3603 35.3626 12.6597
Discrete Meyer DWT 40.6572 35.3626 26.9566 32.5428 35.3626 18.8423

Table 4.

Input MRI .dcm file input having 522 KB size for evaluation of various image quality parameters.

Table 5 indicates that Daubechies, Symlets, and Coiflets DWT give higher compressed image size (9 KB) and highest compression ratio (4) and gain (75%). Decompression is working in this input size (36 KB). So, when the small-size input .dcm file, image reconstruction is required, this implementation work will be applicable. This application gives good PSNR (45.1486) and better MSE (29.88) and SNR (30.0124) for Daubechies DWT.

Sr. no. Wavelet type Size of compressed image (in KB) CR % of CG PSNR (in db) MSE (in db) SNR (in db)
1 Biorthogonal DWT 15 2.4 58.33333 35.2407 19.9722 39.0124
2 Haar DWT 11 3.272727 69.44444 44.8457 29.5771 30.0124
3 Daubechies DWT 9 4 75 45.1486 29.88 30.0124
4 Symlets DWT 9 4 75 44.9406 29.6721 39.0124
5 Coiflets DWT 9 4 75 38.9493 23.6808 39.0124
6 Reverse Biorthogonal DWT 18 2 50 35.7811 20.5126 39.0124
7 Discrete Meyer DWT 15 2.4 58.33333 44.6503 29.3818 39.0124

Table 5.

Compression and image quality performance of input (CT scan) image size 36 KB for different wavelet types.

Table 6 indicates that Biorthogonal DWT gives higher compressed image size (2 KB) and highest compression ratio (18:1) for threshold value = 50. But image quality is degraded due to PSNR value = 32.7496.

Parameters λ =3 λ =10 λ = 20 λ = 30 λ = 50
Compressed image in KB 9 5 4 3 2
Compression ratio 4 7.2 9 12 18
PSNR 35.1707 34.5081 33.673 33.1783 32.7496
MSE 19.9022 19.2396 18.4 17.9098 17.9098
SNR 39.0124 39.0124 39.0124 39.0124 39.0124

Table 6.

Compression and image quality performance of input (CT scan) image size 36 KB (gray scale) for different threshold values on Biorthogonal DWT using proposed method.

Table 7 shows that CT0081 gives good compression result for implemented method, but other input CT images give less compression ratio than .jpg file format. The obtained result shows that.

Sr. no. Image ID Dimension Depth in bit Input file size in KB DICOM size in KB .jpg file size in KB Implemented method .hdwt size in KB
1 CT0014 512 × 512 24 1030 769 25 29
2 CT0051 512 × 200 24 204 301 17 22
3 CT0052 250 × 512 24 254 376 17 26
4 CT0059 350 × 512 24 353 526 23 46
5 CT0074 512 × 512 24 5130 769 61 67
6 CT0081 888 × 733 24 2547 1908 167 165
7 CT0090 512 × 512 24 3591 769 93 152
8 CT0101 512 × 512 24 516 769 92 71
9 CT102 512 × 605 24 609 909 153 162
10 CT110 512 × 512 24 4616 769 61 68

Table 7.

Comparison of implemented method and .JPG format for 10 CT scan DICOM images.

Advertisement

5. Conclusions

Compression and decompression are necessary tasks in medical imaging applications. This implementation provides patient data integration within the medical image. It is very important to maintain the patient data security. This implementation work is very helpful to hide and recover patient information within the medical image and follow compression/decompression without any data loss. In this implanted work, 2D DWT and N-level decomposition are applied on medical image, and then the extracted detail coefficients are firstly encoded by RLE. Secondly, the extracted approximate coefficient and encoded detailed coefficients are encoded by Huffman encoder. The generated result shown that Biorthogonal DWT gives better compression size, compression ratio, and compression gain for higher decomposition level (N = 4), but image quality parameters like PSNR and MSE are degraded. After comparison with JPEG file format and implemented work, this work gives less compression size. We conclude that for medical image compression, we can select N = 3, decomposition level with λ = 3, threshold value, and Biorthogonal DWT for good image quality.

Advertisement

Acknowledgments

Authors thank Dr. S.V. Dudal, HOD, Department of Applied Electronics, SGBA University, Amravati, India, for providing all kind of facilities and support.

Advertisement

Conflict of interest

The authors declare that there is no conflict of interests, financial, potential, or otherwise associated with this manuscript.

References

  1. 1. NEMA Publications. DICOM Standard. Digital Imaging and Communications in Medicine (DICOM). 2008. Available from: ftp://medical.nema.org/medical/dicom/2008/
  2. 2. Cho K, Kim J, Jung SY, Kim K. Development of Medical Imaging Viewer Role in DICOM Standard. Hyun-Kook Kuhng Daegu University, School of Computer and Communication Engineering **ETRI, Broadcasting Media Research Group 0-7803-8940-91051 IEEE; 2005
  3. 3. Mousa WA, Shwehdi MH, Abdul-Malek MA. Conversion of DICOM System Images to Common Standard Image Format Using Matlab. Asia SENSE. SENSOR; 2003. pp. 251–255
  4. 4. Suapang P, Dejhan K, Yimmun S. Medical Image Compression and DICOM-Format Image Archive. Bangkok, Thailand: Department of Telecommunication Engineering, King Mongkut's Institute of Technology Ladkrabang, ICCAS-SICE; 2009
  5. 5. Larobina M, Murino L. Medical image file format. Journal of Digital Imaging. 2014;27(2):200-206. DOI: 10.1007/s10278-013-9657-9
  6. 6. Verma DR. Managing DICOM image: Tips and tricks for radiologist. Indian Journal of Radiology and Imaging. 2012;22(1):004-013. DOI: 10.4103/0971-3026.95396
  7. 7. Graham RNJ, Perriss RW, Scarsbrook AF. DICOM demystified: A review of digital file formats and their use in radiological practice. Clinical Radiology. 2005;60(11):1133-1140. DOI: 10.1016/j.crad.2005.07.003
  8. 8. Baraskar T, Pawar A. Conversion of DICOM image to common standard image formats. International Journal of Emerging Technology and Advanced Engineering. 2013;3:1-5
  9. 9. Ujgare NS, Baviskar SP. Conversion of DICOM image in to JPEG, BMP and PNG image format. International Journal of Computer Applications. 2013;62(11):22-26
  10. 10. Kaur S, Jindal G. Survey of databases used in image processing and their applications. International Journal of Scientific and Engineering Research. 2011;2(10):1-9
  11. 11. Chen P. Study on medical image processing technologies based on DICOM. Journal of Computers. 2012;7(10):1-8. DOI: 10.4304/jcp.7.10.2354-2361
  12. 12. Vinayak B, Gaikwad AN, Kanaskar M. DICOM medical data compression for telemedicine IN rural areas. Advances in Engineering Science. 2008;2:001-006
  13. 13. Zhifeng L, Changhong F, Xu F, Zhicong Q, Shunxiang W. An easy image compression method and its realization base on MATLAB, Information Engineering and Computer science International Conference; 2009
  14. 14. Dimitrovski I, Guguljanov P, Loskovska S. Implementation of web-based medical image retrieval system in Oracle. In: IEEE 2nd Intl. Conference on Adaptive Science & Technology. 2009
  15. 15. Xiaolei SHI, Wang M. Transformation of DICOM digital medical image format into BMP general image format. Microcomputer Informatics. 2010;26:195-197
  16. 16. Aliming H. The conversion between DICOM medical image format and common graphic format [master degree dissertation]. Chengdu, Sichuan: Chinqi Sichuan University; 2006
  17. 17. Cyriac M, Chellamuthu C. Medical image compression using visual quantization and modified run length encoding. Biomedical Imaging and Intervention Journal. 2013;7(2):1-8. DOI: 10.2349/biij9.2.e5
  18. 18. Hussain AJ, Al-Fayadh A, Radi N. Image compression technique: A survey in lossless and lossy algorithms. Neurocomputing. 2018;300:44-69
  19. 19. Akhtar MB, Qureshi AM, Qamar ul Islam. Optimum Run Length Coding for JPEG Image Compression used in Space Research Program of IST. IEEE; 2011
  20. 20. Janet J, Mohandass D, Meenalosini S. Lossless Compression Techniques for Medical Images in Telemedicine-Advances in Telemedicine Technologies, Enabling Factors, Scenarios. Austria: INTECH; 2011. pp. 111-130. ISBN 978- 953307-159-6
  21. 21. Salomon D. Data Compression-the Complete Reference. 2nd ed. Heidelberg: Springer; 2001
  22. 22. Shajun Nisha S, Kothar Mohideen S. Wavelet coefficients Thresholding techniques for Denoising MRI images. Indian Journal of Science and Technology. 2016;9(28):1-8. DOI: 10.17485/ijst/2016/v9i28/93872
  23. 23. Ray A. Performance evaluation of various image compression techniques using SVD, DCT and DWT. International Journal of Research in Engineering and Technology. 2017;6(6):1-6

Written By

Trupti N. Baraskar and Vijay R. Mankar

Submitted: 18 June 2019 Reviewed: 12 August 2019 Published: 31 October 2019