Open access peer-reviewed chapter

VLSI Implementation of Medical Image Fusion Using DWT-PCA Algorithms

Written By

Surya Prasada Rao Borra, Rajesh K. Panakala and Pullakura Rajesh Kumar

Submitted: 15 August 2019 Reviewed: 21 January 2020 Published: 15 April 2020

DOI: 10.5772/intechopen.91298

From the Edited Volume

Field Programmable Gate Arrays (FPGAs) II

Edited by George Dekoulis

Chapter metrics overview

756 Chapter Downloads

View Full Metrics

Abstract

Nowadays, the usage of DIP is more important in the medical field to identify the activities of the patients related to various diseases. Magnetic Resonance Imaging (MRI) and Computer Tomography (CT) scan images are used to perform the fusion process. In brain medical image, MRI scan is used to show the brain structural information without functional data. But, CT scan image is included the functional data with brain activity. To improve the low dose CT scan, hybrid algorithm is introduced in this paper which is implemented in FPGA. The main objective of this work is to optimize performances of the hardware. This work is implemented in FPGA. The combination of Discrete Wavelet Transform (DWT) and Principle Component Analysis (PCA) is known as hybrid algorithm. The Maximum Selection Rule (MSR) is used to select the high frequency component from DWT. These three algorithms have RTL architecture which is implemented by Verilog code. Application Specified Integrated Chips (ASIC) and Field Programmable Gate Array (FPGA) performances analyzed for the different methods. In 180 nm technology, DWT-PCA-IF architecture achieved 5.145 mm2 area, 298.25 mW power, and 124 ms delay. From the fused medical image, mean, Standard Deviation (SD), entropy, and Mutual Information (MI) performances are evaluated for DWT-PCA method.

Keywords

  • application specified integrated chips
  • discrete wavelet transform
  • field programmable gate array
  • principle component analysis
  • maximum selection rule

1. Introduction

In recent years, Image Fusion (IF) importance has increased rapidly. The process of combining two or more images into one image is called as IF. Through this, all kinds of information are possible to take from the different images [1]. Based on the image stage, the fusion has been classified into two types, those are transform domain and spatial domain fusion [2]. IF is used in so many applications like medical, automated industry, engineering field, military, etc. [3]. Among all those fields, medical field application is more important in IF which helps to identify the human problems [4]. In medical, two major models like MRI and CT scan helps to analyze the normal and abnormal tissue and internal structure of human body because both MRI and CT contain some different information of the human brain [5]. MRI scan is used for soft tissue which detects the skull problems as well as CT scan is used for hard tissue to identify the bone structure [6]. Earlier many techniques used in IF like pixel level based, decision level, and feature level based [7]. Many of the existing algorithm has been used for IF process such as Electrical Capacitance Tomography (ECT) algorithm [8], Non-Subsample Contour let Transform (NSCT) [9], sparse representation and decision [10], Curvelet transform [11], hybrid Entropy concept [12], hybrid Dual tree complex wavelet transform [13], and hybrid IF and image registration [14]. The main problem with these methods is information loss. To check the hardware utilization and improve the efficiency, the IF has been implemented in FPGA. The way of implementation is also different in FPGA. In FPGA, DWT [15], multi model method [16], and configurable pixel level [17] methods have been implemented for IF process. The hardware utilization of these methods is high. To overcome these problems, hybrid algorithm with the maximum selection rule is implemented in this paper. From the DWT, high frequency component signal only processes the MSR and output of this is given to the Inverse DWT. The combination of DWT and PCA is named as hybrid algorithm. The PCA output gives the IF output. These methods implemented in FPGA architecture to improve the efficiency of the IF. At last, FPGA and ASIC performances improved in proposed method compared to conventional methods. Mean, Standard Deviation (SD), Entropy, and Mutual Information (MI) performances also calculated for all the algorithms. The rest of the paper is organized as follows: Section-2 elaborates the literature survey, Section-3 describes the proposed method, Section-4 discusses the experimental results, and Conclusion is summarized in section-5.

Advertisement

2. Literature review

Mishra et al. [18] presented Modified Frei-chen based image fusion method. This method was utilized in Structural Similarity (SS), and contrast in Night Vision (NV) based two-scale decomposition. This method achieved 48%, 15%, and 100% of improvements in total edge transfer, SS, and NV. This architecture was implemented in the Xilinx tool which consumes 4% of resources. This proposed method was analyzed in synopsis tool with 90 nm CMOS technology. This algorithm provides less accuracy and less fusion efficiency.

Bavirisetti and Dhulli [19] proposed two scale image fusion using saliency detection. This method was used for Saliency extraction process, which can highlight the significant information. This works gave better results compared to multi-scale fusion technique. This method failed to process the medical images perfectly.

Pemmaraju et al. [20] presented wavelet based image fusion using FPGA. This proposed method was implemented in Xilinx EDK 10.1 using Spartan 3E. This FPGA contains combinational blocks which are flexible for high speed application. This architecture contains memory, flip flops, and LUT. This proposed method was applied to multi focus image fusion. DWT does not provide stationary outputs and low frequency component has less efficiency.

Yang et al. [21] proposed multi model based image fusion based on fuzzy logic. With the help of type 2 fuzzy, NSCT was analyzed using pre-registered source image for getting low and high bands. Low frequency bands are used by local energy algorithm. The proposed fused image was taken with the help of inverse NSCT with all sub bands. The accuracy, contrast, and versatility was also evaluated. The main drawback of this method is low spatial resolution.

Bhaskar and Munde [22] proposed image fusion using Non-Subsampled Shear-let Transform (NST) in FPGA implementation. Input image was separated into individual image co-efficient using NST. Different rules were applied to fuse the high and low bands. With the help of inverse NST, the fused image was taken. This proposed method was implemented in Xilinx system generator and MATLAB. The power value was reduced in proposed method. But, the hardware utilization of this proposed method is high.

Agarwal and Bedi [23] presented hybrid image fusion for medical diagnosis. In this paper, wavelet and Curvelet transforms were used to perform the IF. The segmented blocks were fused into sub bands using Curvelet transform. The resolution of the fused image is too less which affects quality of the image.

Sanjay et al. [24] proposed IF based on DWT and type-2 fuzzy logic. In this paper, CT and MRI images were fused with the help of hybrid method. The fused low level bands and high level bands were reconstructed to perform the IDWT. This hybrid algorithm fails to use more logic function and analyses the hardware utilization.

Advertisement

3. DWT-PCA-IF architecture

Image Fusion is one of the important processes for obtaining more information from different images. The overall process of image fusion is shown in Figure 1.

  • The input CT image is read into MATLAB and the pixel is converted to binary value. These binary values are stored in a text file.

  • The same process is applied to MRI images also.

  • Both CT and MRI images binary values perform the DWT process which gives four frequency components such as Low High (LH), High Low (HL), High High (HH), and Low Low (LL).

  • These frequency components perform MSR. In this operation, high frequency component only required.

  • So, HH, HL, and LH frequency components performed MSR operation which gives three results.

    These three results are given to the IDWT process along with low frequency component (LL).

  • After performing IDWT, both results are given to the PCA component which gives the fused image.

  • DWT, MSR, IDWT and PCA are implemented in Verilog and the final output is written in text file.

  • With the help of MATLAB, that binary values are converted to pixel which shows the fused image.

Figure 1.

Block diagram of entire process.

3.1 DWT architecture

For analyzing the signal, wavelet converts the time domain to frequency domain. The DWT is implemented using two major blocks namely Filter Bank (FB) and Lifting Scheme (LS). The DWT is a decimated wavelet transform, where the size of the image reduces by half at each scale. It is easy to convert the spatial domain inputs into frequency domain in wavelet transform [25]. High pass and low pass coefficient series are obtained from the input series y0, y1, …. yn. The high pass and low pass coefficients are represented by using the following two Eqs. (1) and (2).

Hi=n=0l12jn.snzE1
Li=n=0l12jn.tnzE2

where, the wavelet filters are represented as sn(z) and tn(z), and length of the filter is denoted as l and j = 0, …, [n/2] − 1.

The spatial domain DWT is applied in two directions. First, 1D-DWT is applied on the horizontal axis and that results are applied to the vertical axis of 1D-DWT. There are four parts named as LL, LH, HL and HH obtained from the 2D-DWT.

The two-dimensional DWT applies to all the rows and columns of an image. If the input image is of size 2k × 2k pixels at level L + 1 its size will be 2k/2 × 2k/2 The various kinds of decomposition methods are used in wavelets over an image. The DWT is applied to the input image, which is decomposed into four sub image. These sub images are named as sub bands. The LL sub band is the coarse level sub image, HH, LH, and HL are the diagonal, vertical and horizontal components of the image respectively. Finally, the input image is decomposed into four major components that is shown in Figure 2. A high level 2D-DWT is developed by LL frequency and low pass components for multi resolution analysis.

Figure 2.

Discrete wavelet transform.

Let assume input image is Y.

Here, Y is splitting into two different bands such as Yo and Ye.

Yo=[Y(1), Y(3), Y(5)Y(2n1)]E3
Ye=[Y(2), Y(4), Y(6)Y(2n)]E4
Q1(n)=Yo(n)+a(Ye(n)+Ye(n+1))E5
V1(n)=Ye(n)+b(Q1(n)+Ye(n+1))E6

dc(n)=L.Q1(n) Here, Q1(n) is scaled by L.

The 2D-DWT architecture and 1D-DWT are shown in Figures 3 and 4. The control signals represent as clk and rst. The odd input and even input are mentioned as odd_in [7:0] and even_in [7:0]. These two inputs are given to the line buffer to perform even and odd extraction which outputs are given to the PIPO for capturing the data. From that block, four outputs are generated which is given to the lifting block. After processing the lifting block, the final output is generated as detailed co-efficient dc_out [27:0] and significance co-efficient sc_out [23:0].

Figure 3.

2D-DWT architecture.

Figure 4.

1D-DWT architecture.

3.2 Maximum selection rule

The MSR diagram is shown in Figure 5. This rule is applicable for the high frequency component. So that HH, HL, LH frequency values perform the MSR operation. Both DWT output values are connected to the MUX for choosing maximum value.

Figure 5.

Maximum selection rule diagram.

These outputs are given to the IDWT for changing the frequency domain into the time domain.

3.3 PCA architecture

The architecture of PCA is shown in Figure 6 which contains control engine, covariance matrix, MUX, multiplier, adder, and comparator. With the help of detected spike waveform, the covariance matrix is calculated. The covariance matrix is called as PC spike waveform. The MAC address is used for distilling and orthogonalization process to improve the PCA efficiency. Comparator and right shift are used to shift the procedure and level checking. The entire algorithm split into four processing units and the data is stored in register files. Finite State Machine (FSM) is used for scheduling and allocating the resources during the PCA processing. FSM is very effective for controlling the remaining signal [26]. These outputs are helpful to perform the image fusion. The fused architecture binary output is read in MATLAB for showing fused image.

Figure 6.

PCA architecture.

Advertisement

4. Experimental results and discussion

In this section, the experimental simulation results and discussion of the proposed methodology is detailed effectively in terms of performance measure. The performance of the proposed methodology was evaluated by ASIC and FPGA performances.

4.1 Discussion

The input images (CT and MRI) are shown in Figures 7 and 8. These images are converted to binary which are shown in Figures 9 and 10.

Figure 7.

Input CT image.

Figure 8.

Input MRI image.

Figure 9.

Binary value of CT image.

Figure 10.

Binary value of MRI image.

The ASIC performance of the different methods are tabulated in Table 1. In this table, values of ASIC performance of the Existing-I [18], existing-II [20], existing-III [22], and DWT-PCA-IF are compared.

TechnologyMethodArea (mm2)Power (mW)Delay (ms)
180 nmExisting-I [l8]8.471387.1180
Existing-II [20]7.321345.71158
Existing-III [22]6.214314.21143
DWT-PCA-IF5.145298.25124
45 nmExisting-I [l8]3.014198.25104
Existing-II [20]2.987168.12101
Existing-III [22]2.158148.68798
DWT-PCA-IF1.982111.2191

Table 1.

Comparison of area, power, and delay for different methods.

The comparison of ASIC performances is tabulated in Table 1. Here, all the methods are implemented and the results are tabulated. All the methods are implemented in the cadence RTL compiler with 180 and 45 nm technology. From this table, it’s clear that DWT-PCA-IF provides better performances when compared to previous existing architectures.

4.2 Comparative analysis

In this work, three papers have been compared with proposed method. A. Mishra, S. Mahapatra, and S. Banerjee [18], applied modified Frei-chen operator based IF for real time applications. Scalable decomposition was used to perform the fusion operation which was implemented in Virtex 4 FPGA. The overall architecture RTL was too complex to perform the IF algorithm which caused more area. Pemmaraju et al. [20], implemented IF based on DWT using FPGA. This algorithm was implemented in Xilinx EDK 10.1 FPGA Spartan 3E hardware. There is no explanation of RTL architecture, and .ucf file. Due to use of wavelet, the power consumption is too high. Bhaskar and Munde [22] performed IF based on non-subsampled shearlet transform. Xilinx system generator was used to implement this design with MATLAB tool. The fused image affected by more noise and it require more hardware utilization.

The comparison graph of area, power, and delay are shown in Figures 1113. The dark blue bar graph is represented as DWT-PCA-IF architecture. All the ASIC performance is reduced due to the hybrid algorithm.

Figure 11.

Comparison of area performance for different methods.

Figure 12.

Comparison of power performance for different methods.

Figure 13.

Comparison of delay performance for different methods.

The FPGA performances are tabulated in Table 2. In this table, Virtex 4 and Virtex 5 devices are used to evaluate LUT, flip flop, slices, and frequency. These values are shows that the DWT-PCA-IF architecture achieves better FPGA performance parameters.

DevicesMethodLUTFlip flopslicesFrequency
Virtex 4Existing-I [l8]403848522857250.3
Existing-II [20]400246572654289.64
Existing-III [22]354142142011314.21
DWT-PCA-IF301439871968355.14
Virtex 5Existing-I [l8]310441251964185.41
Existing-II [20]301440321847193.21
Existing-III [22]298739871752255.14
DWT-PCA-IF274137891648287.96

Table 2.

Comparison of FPGA performances for different methods.

The comparison graph of LUT, Flip flop, slices, and frequency are shown in Figures 1417. The hardware utilizations are evaluated from this FPGA performance. The top module and 2D DWT and 1D DWT RTL schematic diagram are shown in Figures 18 and 19.

Figure 14.

Comparison of LUT for different methods.

Figure 15.

Comparison of flip flop for different methods.

Figure 16.

Comparison of slices for different methods.

Figure 17.

Comparison of frequency for different methods.

Figure 18.

Top module of DWT architecture.

Figure 19.

Internal schematic of 2D-DWT.

The performance evaluation for different methods is given in Table 3. Here, some of the performances are evaluated such as Mean, Standard Deviation (SD), Entropy, and Mutual information (MI). This performance evaluated for fused medical image. From this table, it is clears that DWT-PCA gives better performances than existing methods. Finally, the fused image is shown in Figure 20. The above RTL schematics are taken from the Xilinx tool.

ImagePerformanceDWT [2]Haar [3]Kekre’s wavelet [7]DTCWT [13]PCA [24]DWT-PCA
Fused ImageMean44.2532.5332.4145.1453.2255.65
SD40.1436.0734.8251.2437.4453.14
Entropy8.1455.975.910847.216.639.621
MI0.1470.390.55412.120.28323.141

Table 3.

Performance evaluation for different methods.

Figure 20.

Fused image.

Advertisement

5. Conclusion

The proposed architecture has been designed effectively in order to reduce the hardware utilization. In this work, DWT-PCA-IF architecture has been designed to perform the image fusion. In this work, medical images like MRI and CT have been used in the fusion process to obtain more information. The hybrid VLSI architecture provided better fused image compared to previous works. The DWT-PCA-IF architecture was implemented using Verilog code. DWT and PCA method were used to reduce the power and area consumption. The ASIC and FPGA performance were analyzed for different architectures. In 180 nm technology, DWT-PCA-IF architecture achieved 5.145 mm2 area, 298.25 mW power, and 124 ms delay. In Virtex 4, the proposed architecture achieved 3014 LUT, 3987 flip flop, 1968 slices, and 355.14 MHz frequency. From the fused image, 55.658 mean, 53.14 SD, 9.621 entropy, and 3.141 MI value has been evaluated. In the future, different kind of optimization algorithm will be designed to improve the ASIC and FPGA performances.

Advertisement

Acknowledgments

At the outset, I would like to take this as an opportunity to convey my gratitude to Intechopen publishing house and their team to for their consistent support at every step in bringing out this chapter in their book. The process that followed in reviewing this chapter and giving valuable review remarks helped a lot to meet the standards of this book and IntechOpen publishing house has enriched my writing skills. I would like to extend my sincere gratitude to my Ph.D. supervisor Dr. Rajesh K Panakala and Dr.P.Rajesh Kumar for valuable guidance and continuous encouragement in publishing this chapter. I would like to thank my family members for their love and support and the management of our college, PVP Siddhartha Institute of Technology for their constant encouragement to carryout my research work.

References

  1. 1. Mahajan S, Singh A. A comparative analysis of different image fusion techniques. IPASJ International Journal of Computer Science. 2014;2(1):008-015
  2. 2. Hussain DK, Reddy CL, Kumar VA. Implementation of medical image fusion using DWT process on FPGA. International Journal of Computer Applications Technology and Research. 2013;2(6):676-679
  3. 3. Phanindra P, Babu JC, Shree VU. VLSI implementation of medical image fusion using Haar transform. International Journal of Scientific and Engineering Research. 2013;4(9):1437-1442
  4. 4. Yang B, Li S. Pixel-level image fusion with simultaneous orthogonal matching pursuit. Information Fusion. 2012;13(1):10-19
  5. 5. Pavithra C, Bhargavi S. Fusion of two images based on wavelet transform. International Journal of Innovative Research in Science, Engineering and Technology. 2013;2(5):1814-1819
  6. 6. Jose B, Kumar BS. Design of 2-D DWT VLSI architecture for image processing. International Journal of Engineering Research and Technology. 2014;3(4):692-696
  7. 7. B. Kekre H, Sarode T, Dhannawat R. Implementation and comparison of different transform techniques using kekre's wavelet transform for image fusion. International Journal of Computer Applications. 2012;44(10):41-48
  8. 8. Olmos AM, Botella G, Castillo E, Morales DP, Banqueri J, García A. A reconstruction method for electrical capacitance tomography based on image fusion techniques. Digital Signal Processing. 2012;22(6):885-893
  9. 9. Bhatnagar G, Wu QMJ, Liu Z. Directive contrast based multimodal medical image fusion in NSCT domain. IEEE Transactions on Multimedia. 2013;15(5):1014-1024
  10. 10. Fei Y, Wei G, Zongxi S. Medical image fusion based on feature extraction and sparse representation. International Journal of Biomedical Imaging. 2017;2017:1-11
  11. 11. Tank VP, Shah DD, Vyas TV, Chotaliya SB, Manavadaria MS. Image fusion based on wavelet and Curvelet transform. IOSR Journal of VLSI and Signal Processing. 2013;1(5):32-36
  12. 12. Sharmila K, Rajkumar S, Vijayarajan V. Hybrid method for multimodality medical image fusion using discrete wavelet transform and entropy concepts with quantitative analysis. In: Proceedings of International Conference on Communications and Signal Processing. 2013. pp. 489-493
  13. 13. Gurjar R. Hybrid image fusion implemented in DTCWT. International Journal of Engineering Technology and Computer Research. 2014;2(1):688-692
  14. 14. Bhosle DS, Gorde KS. Image registration and wavelet based hybrid image fusion. IOSR Journal of VLSI and Signal Processing. 2014;4(2):1-5
  15. 15. Suraj AA, Francis M, Kavya TS, Nirmal TM. Discrete wavelet transform based image fusion and de-noising in FPGA. Journal of Electrical Systems and Information Technology. 2014;1(1):72-81
  16. 16. Kaur R, Kaur S. An approach for image fusion using PCA and genetic algorithm. International Journal of Computer Applications. 2016;145(6):54-59
  17. 17. Besiris D, Tsagaris V, Fragoulis N, Theoharatos C. An FPGA-based hardware implementation of configurable pixel-level color image fusion. IEEE Transactions on Geoscience and Remote Sensing. 2012;50(2):362
  18. 18. Mishra A, Mahapatra S, Banerjee S. Modified Frei-Chen operator-based infrared and visible sensor image fusion for real-time applications. IEEE Sensors Journal. 2017;17(14):4639-4646
  19. 19. Bavirisetti DP, Dhuli R. Two-scale image fusion of visible and infrared images using saliency detection. Infrared Physics & Technology. 2016;76:52-64
  20. 20. Pemmaraju M, Mashetty SC, Aruva S, Saduvelly M, Edara BB. Implementation of image fusion based on wavelet domain using FPGA. In: Proceedings of International Conference on Trends in Electronics and Informatics. 2017. pp. 500-504
  21. 21. Yang Y, Que Y, Huang S, Lin P. Multimodal sensor medical image fusion based on type-2 fuzzy logic in NSCT domain. IEEE Sensors Journal. 2016;16(10):3735-3745
  22. 22. Bhaskar PC, Munde MV. FPGA implementation of non-subsampled Shearlet transform for image fusion. In: Proceedings of International Conference on Computing, Communication, Control and Automation. 2017. pp. 1-6
  23. 23. Agarwal J, Bedi SS. Implementation of hybrid image fusion technique for feature enhancement in medical diagnosis. Human-Centric Computing and Information Sciences. 2015;5(1):3
  24. 24. Sanjay AR, Soundrapandiyan R, Karuppiah M, Ganapathy R. CT and MRI image fusion based on discrete wavelet transform and Type-2 fuzzy logic. International Journal of Intelligent Engineering and Systems. 2017;10(3):355-362
  25. 25. Surya PRB, Panakala RK, Kumar PR. Hybrid image fusion algorithm using DWT maximum selection rule and PCA. International Journal of Scientific and Engineering Research;8(8):814-820
  26. 26. Surya PRB, Panakala RK, Kumar PR. Qualitative analysis of MRI and enhanced low dose CT scan image fusion. In: Proceedings of International Conference on Advanced Computing and Communication Systems. 2017. pp. 1752-1757

Written By

Surya Prasada Rao Borra, Rajesh K. Panakala and Pullakura Rajesh Kumar

Submitted: 15 August 2019 Reviewed: 21 January 2020 Published: 15 April 2020