Distortion-free compression of medical images based on SPIHT algorithm

A coding / decoding scheme combining the second-generation wavelet transform and SPIHT algorithm in image lossless compression is proposed. The experimental results of this scheme and the compression performance comparison with other methods are given. Research shows that this method is suitable for medical images. Has a higher lossless compression rate.
Keywords: second-generation wavelet transform, lossless compression, SPIHT, medical images


1 Introduction With the development of society and the advancement of medical technology, people are more and more concerned about their health. Medical imaging is no longer information for doctors' reference only and has become an important basis for diagnosing diseases. Image compression coding under the condition of network transmission has become the key technology for establishing digital hospitals. At present, the compression standards for two-dimensional images include JPEG, GIF and JPEG2000 using wavelet transform. Medical images are special, and they generally do not allow the loss of useful details. The traditional DCT (Discrete Cosine Transform, discrete cosine transform) and the first generation wavelet will produce floating-point numbers after the image transformation, so the transformed data must be quantized, which will produce different degrees of distortion. It can be seen that the design of the quantizer is the key link that determines the image fidelity. Since the second-generation wavelet adopts the lifting method to realize the integer transform, it can realize the lossless compression of the image. Obviously, it is a compression method very suitable for medical images.
2 SPIHT algorithm Set ParTITIoning in Hierarchical Trees (SPIHT) based on hierarchical tree improves the embedded zero tree coding algorithm (EZW). After the wavelet transform is performed on the image, it more effectively utilizes the similarity between the important coefficients of sub-bands of different scales. It exhibits good characteristics: it constructs wavelets in the spatial domain without relying on the Fourier transform; a high PSNR (Peak Signal Noise RaTIo) guarantees a good quality of the reproduced image; integer arithmetic facilitates real-time rapid editing Decoding and network transmission; the gradual presentation of the image code stream is convenient for users to retrieve images of interest on the Internet.
The SPIHT algorithm uses the following coding steps for image information.
First, define three cohorts: insignificance coefficient queue LIP, saliency coefficient queue LSP, and insignificance set queue LIS.
Let O (i, j) denote the set of direct nodes of node (i, j); D (i, j) denote the set of child nodes of node (i, j); L (i, j) denote the exclusion of child nodes The collection after the direct node.
In the queue, each element is uniquely identified by a coordinate. It represents the isolation coefficient (root node without child nodes) in LIP and LSP, and D (i, j) or the second type of the first type element in LIS. L (i, j) of the element.
Significance test for a certain threshold T. Move elements greater than T into the LSP and remove the element in the LIP queue. The same test was performed on LIS, moving significant elements into the LSP, and then splitting the tree.

The SPIHT algorithm described in C ++-like language is as follows:
In the first step, the threshold T and the three queues (LSP, LIS, and LIP) are initialized.

coordinate of;


(2) if (x, y) is the second type of element, the significance test for L (i, j) if (L (i, j)) = = 1 all (k, l) ∈ O (i, j ) Move into LIS as the first type of element, and leave the team from LIS.
The third step is bit transmission / storage. Convert each coefficient in LSP into binary transmission / storage.
In the fourth step, the threshold value is updated and goes to the second step: T / = 2; gotostep2.
3 Lifting scheme and second-generation wavelet lifting method The construction of wavelet is divided into three steps: splitting, predicting and updating.
3.1 Split
Divide an original signal sequence Sj into two smaller, disjoint wavelet subsets Sj-1 and dj-1 according to even and odd numbers:

3.2 Prediction (predict)
Due to the correlation between the data, a prediction operator P can be defined so that dj-1 = P (Sj-1), so that adjacent even sequences can be used to predict odd sequences. If the difference between dj-1 and P (Sj-1) is used to replace dj-1, the data volume is much smaller than the original dj-1.

In the simplest case, the average value of the data of two adjacent even serial numbers is taken as the predicted value of the data of the odd serial numbers between them. which is,

3.3 Update
Since the above two processes generally cannot maintain some of the overall properties (such as brightness) in the original image, for this reason, we have to construct a U operator to update Sj-1 so that it maintains some characteristics of the original data set.

4 Encoding / decoding scheme In this paper, the front-end uses the second generation wavelet (lifTIng wavelet), then uses the SPIHT algorithm for the wavelet coefficients, and then uses Amir Said's adaptive arithmetic coding. Decoding is the inverse process of encoding, which includes three steps corresponding to forward SPIHT: recovery update, recovery prediction, and merge. The encoding / decoding scheme is shown in Figure 1.



If the front end uses the first generation wavelet for lossy compression, a higher compression ratio can be achieved. Obviously, the requirements of high fidelity and high compression ratio for data compression in the second generation wavelet transform are contradictory.
5. Experimental results and conclusions For the previous coding scheme, we tested the medical image and the Lena image separately, and the bit rate bbp used bit / pixel. Because of the lossless compression scheme, the three different encoding methods in Table 1 all have PSNR = ∞.

It can be seen from Table 1 that the performance difference when encoding the standard test image Lena is not very large, but because there are a large number of "zero pixels" at the edges of general medical images, a large number of "zero pixels" can be generated when encoding with SPIHT "Zero Tree" greatly reduces the amount of data. Therefore, when compressing medical images, the method of this article is more suitable.
Further analysis shows that, compared with the currently widely used JPEG, this compression scheme occupies less memory, has a high coding efficiency, and has no mosaic phenomenon. At low bit rates, the gap between the two is even more pronounced. If the scheme is implemented with parallel fast algorithms and hardware, its real-time performance will be further improved, so this medical image compression scheme has better application prospects.



references
1 A Said, WA Pearlman. A New Fast and Efficient ImageCode Based on Set Partitioning in Hierarchical Trees. IEEE Transactions on Circuits and Systems for VideoTechnology, 1996, 6 (3)
2 Erik Sjoblom. Compression of Medical Image Stacks Using Wavelets and Zero-Tree Coding. [Master Thesis]. Linkoping: Linkoping University, 2002
3 Xu Peixia. Wavelet analysis and application examples. Hefei: China University of Science and Technology Press, 2001
4 Chen Wufan. Wavelet analysis and its application in image processing. Beijing: Science Press, 2002
5 Hu Xuelong, Jiang Nan, Guo Zhenmin. Digital hospital image archiving and communication system (PACS). Journal of Engineering Graphics, 2001 (increasing)
6 Luo Wenbin, Zhang Liming. The application of the second generation wavelet in lossless compression of medical images. Computer Engineering, 1999 (10)
7 Tang Yan, Mo Yulong. The second generation wavelet transform is applied to lossless compression coding of images. Journal of Graphics and Images of China, 2000 (8)

7 Inch Tablet

When people talk about 7inch tablet, kids appears on their minds, especially for elementary students or kindergarten kids for playing intelligent exploitation games or online learning. Clients usually choose 7 inch tablet wifi only as 7 inch educational tablet for project, wifi one is much cheaper. Of course, 7 inch Android Tablet with 3G lite or 4G lite also optional. You can always see a right tablet at this store, no matter amazon 7 inch tablet, 8 inch android tablet, or 10.1 android tablet.

Except android tablet and window tablet, there are Education Laptop, Mini PC , All In One PC, which are is the main series at this store. Any other special configuration interest, just email or call us, you will receiving value information in 1-2 working days.

To meet clients` changing requirements, we contributes 10-20% profit to develop new designs according market research and client`s feedback.

So you are always welcome if can share your special demand or your clients opinion for the products.

7 Inch Tablet,7 Inch Android Tablet,Amazon 7 Inch Tablet,7 Inch Tablet Wifi Only,7 Inch Educational Tablet

Henan Shuyi Electronics Co., Ltd. , https://www.hnsylabtop.com