Next Article in Journal
Implementing a Chaotic Cryptosystem by Performing Parallel Computing on Embedded Systems with Multiprocessors
Next Article in Special Issue
M-ary Rank Classifier Combination: A Binary Linear Programming Problem
Previous Article in Journal
Descriptions of Entropy with Fractal Dynamics and Their Applications to the Flow Pressure of Centrifugal Compressor
Previous Article in Special Issue
Efficient Low-PAR Waveform Design Method for Extended Target Estimation Based on Information Theory in Cognitive Radar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Dictionary Construction Based Multimodal Medical Image Fusion Framework

1
School of Instrumentation and Optoelectronic Engineering, Beihang University, Beijing 100083, China
2
School of Physics and Optoelectronic Engineering, Foshan University, Foshan 528000, China
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(3), 267; https://doi.org/10.3390/e21030267
Submission received: 12 January 2019 / Revised: 22 February 2019 / Accepted: 4 March 2019 / Published: 9 March 2019
(This article belongs to the Special Issue Information Theory Applications in Signal Processing)

Abstract

:
Training a good dictionary is the key to a successful image fusion method of sparse representation based models. In this paper, we propose a novel dictionary learning scheme for medical image fusion. First, we reinforce the weak information of images by extracting and adding their multi-layer details to generate the informative patches. Meanwhile, we introduce a simple and effective multi-scale sampling to implement a multi-scale representation of patches while reducing the computational cost. Second, we design a neighborhood energy metric and a multi-scale spatial frequency metric for clustering the image patches with a similar brightness and detail information into each respective patch group. Then, we train the energy sub-dictionary and detail sub-dictionary, respectively by K-SVD. Finally, we combine the sub-dictionaries to construct a final, complete, compact and informative dictionary. As a main contribution, the proposed online dictionary learning can not only obtain an informative as well as compact dictionary, but can also address the defects, such as superfluous patch issues and low computation efficiency, in traditional dictionary learning algorithms. The experimental results show that our algorithm is superior to some state-of-the-art dictionary learning based techniques in both subjective visual effects and objective evaluation criteria.

1. Introduction

In recent decades, various medical imaging techniques have emerged and numerous medical images have been provided. Different modal medical images can only offer a limited characteristic description about human organs due to the diverse imaging mechanism [1]. Medical image fusion technology can integrate the complementary information of multiple single-modal medical images acquired by different imaging sensor modalities and provide a more precise, reliable, and better description of lesions [2]. This technology has witnessed various applications in real products such as medical assistance, as it can provide doctors a comprehensive description on disease tissues, so as to develop appropriate treatment plans and improve diagnostic accuracy.
In past decades, various research efforts have been made to develop an effective medical image fusion technique [3]. In general, according to the characteristics of the medical images, multi-scale transform (MST) based image fusion methods have been the most popular trend due to their simplicity, feasibility, and effectiveness in implementation [4,5,6]. A number of typical MSTs were developed, for example: Stationary Wavelet Transform (SWT) [7], Dual-Tree Complex Wavelet Transform (DTCWT) [8], Curvelet Transform (CVT) [9], Contourlet Transform (CT) [10], Nonsubsampled Contourlet Transform (NSCT) [11,12], etc. The MST-based methods first decompose the source image into a series of sub-images to express the base and detail information with different levels and directions and then fuse the lowpass and high-pass coefficients according to the designed fusion rules respectively. Finally, the fusion result can be obtained by performing the inverse MST, although the MST-based methods can effectively extract some important features from the source image into the fused image. However, these methods cannot adaptively express the content of the image, which heavily limits their fusion performance, meanwhile, each MST tool has its own merits and limitations depending on the context of input images [3]. Furthermore, some useful information may be inevitably lost in the process of decomposition and reconstruction, and the suboptimal results would be generated. Thus choosing an optimal MST is not an obvious and trivial task as it relies on scene contexts and applications.
In order to address the deficiencies of MSTs based methods and obtain an encouraging fusion performance, a plethora of carefully designed medical image fusion algorithms have been proposed, such as neural network (NN) based methods [13,14,15], principal component analysis (PCA) based techniques [16], and mathematical morphology (MM) based methods [17] and so on. However, the NN based method relies too much on a large number of parameters manually set, which is not conducive to the adaptive implementation of the fusion process. The PCA based technique is prone to produce the drawback of spectral distortion. Although the MM based algorithm plays an important role in image fusion, some of the details in source image may be smoothed in the final fusion result, which would affect the fusion performance and hinder the applications of them.
Sparse representation (SR) theory can describe and reconstruct images in a sparse and efficient way by linear combination of sparse coefficients and overcomplete dictionary [18]. In recent years, SR theory has been rapidly developed and successfully applied in many image processing applications, such as image super-resolution reconstruction [19,20], image feature extraction [21], image denoising [22,23], pedestrian re-identification [24,25], image classification [26], and many other fields. At the same time, it has been successfully applied to many image fusion fields and achieved some satisfactory results [27,28,29]. In SR, the overcomplete dictionary is of important significance in medical image representation. It also plays a key role in affecting the quality of fused image. Yang and Li [30] took the first step for utilizing SR theory into the field of multi-focus image fusion, in which DCT was used as a fixed dictionary. Subsequently, they proposed a SR based on Simultaneous Orthogonal Matching Pursuit (SOMP) and successfully applied for multimodal image fusion [31]. Jiang et al. [17] regarded the image as a superposition of two different components of cartoon and texture, respectively, they employed Curvelet and DCT basis as two different dictionaries to express their respective information, and proposed an image fusion method based on SR and morphological component analysis. However, owing to the fact that dictionaries produced by these methods are constructed by the fixed base functions with poor adaptability, this is not an effective way to describe the complex structure information of medical images [18].
Relative to analytical fixed dictionaries, dictionary learning based methods use a small amount of atoms from a trained dictionary instead of a predefined one [27], which can produce state-of-the-art results in many image processing and recognition tasks [2]. Learning from the training sample images to obtain an informative and compact overcomplete dictionary can enhance the adaptability of the dictionary, as well as accurately express the medical image information. In order to produce better fusion quality, many recent dictionary learning based fusion algorithms have been proposed [15]. Meanwhile, some SR based methods that directly train the source images to train the overcomplete dictionary can generate promising fusion results [31]. However, many methods suffered from a general problem that taking all of the patches for dictionary learning, which will unavoidable, would introduce lots of unvalued and redundant information during dictionary training and decrease the fusion performance in medical image fusion. To address this shortcoming in the literature [32], Zhu et al. developed a local density peak-clustering algorithm to refine the process of patch selection, and then established a compact dictionary by K-SVD algorithm. Kim et al. [33] employed a joint clustering technique to classify image patches according to their similar structures, the overcomplete dictionary can be obtained by combing principal components of them. In order to capture the intrinsic features of the image and preserve the hierarchical structure of stationary wavelets, Yin [34] designed a joint dictionary learning method based on all base subbands of stationary wavelets. Liu et al. [35] proposed an adaptive SR method for simultaneous fusion and denoising the images. In this method, they classified a lot of high-quality image patches with gradient information into several categories, and then trained each sub-dictionary of them. Qi et al. [36] proposed an entropy based image fusion method. In this method, the source images were decomposed into low frequency and high frequency images respectively, then the weighted average scheme was utilized to fuse the low-frequency part and the entropy based dictionary learning technique was introduced to fuse the high-frequency part. However, the average scheme for the low-frequency component would inevitably lose some energy of the input image, which would decrease some useful brightness information in fused results. Although comparing with some traditional medical image fusion techniques, the medical image feature can be represented effectively and completely in these methods [32,33], the fusion performance remaining has much room for improvement.
The purpose of the medical image fusion technique is to retain useful information from source images into a fused result as much as possible [37]. To address the shortcomings of low-learning efficiency and weak dictionary expression ability in traditional dictionary learning based algorithms, this paper proposes a novel dictionary learning method based on the brightness and detail clustering for medical image fusion. Our method consists of three steps, the enhancement and multi-scale sampling on images, patch feature classification and dictionary construction, and image fusion. Firstly, we enhance the details of the source images by multi-layer filtering technology, which can significantly improve the information expression ability of dictionary. In the meantime, we conduct a multi-sampling scheme on enhanced images. This operation can significantly increase the richness of patches while decreasing computation load in the dictionary learning stage. Secondly, in order to obtain a fused result with abundant brightness and detail features, i.e., the successful expression of the underlying visual features of the medical image, we developed two feature clustering criteria to classify all patches into two categories including a brightness group and a details group, and construct the energy sub-dictionary and the detail sub-dictionary by K-SVD algorithm. These sub-dictionaries can be directly built by the final compact and informative dictionary. Finally, the sparse coefficient can be calculated by the orthogonal matching pursuit (OMP) algorithm [38], and then reconstructed to the fused images. The main contributions of this paper can be elaborated as follows:
  • We conduct multi-level neighbor distance filtering to enhance the information and take multi-scale sampling to realize the multi-scale expression of images, which can make image patches more informative and flexible, while not increasing the computational complexity in the training stage.
  • Based on the characteristics of the human visual system processing medical images, we develop novel neighborhood energy and a multi-scale spatial frequency to cluster brightness and detail patches, and then train the brightness sub-dictionary and detail sub-dictionary, respectively.
  • A feature discriminative dictionary is constructed by combing the two sub-dictionaries. The final dictionary contains important brightness and detail information, which can effectively describe the useful feature of medical images.
The rest of the paper is organized as follows: The proposed dictionary learning method for medical image fusion is described in detail in Section 2, including reviews the basic theory of SR, dictionary learning and SR subsections. The experiments and results analysis are presented in Section 3. Conclusions and discussion are summarized in Section 4.

2. Proposed Framework

In this section, we present the proposed technique in detail. There are three sub-sections in this part, which include the theory of SR, dictionary learning and image fusion. In the dictionary step, we propose two techniques for images to generate informative training images. Then, the image patches are clustered into brightness and detail groups based on the neighborhood energy and the multi-scale spatial frequency, respectively. The K-SVD algorithm is employed to construct a brightness-based sub-dictionary and a detail-based sub-dictionary, and the final dictionary comsists of the sub-dictionaries. In the image fusion step, a SR model is well established and the sparse coefficients can be calculated by OMP, then generate the final fused image.

2.1. Sparse Representation

For a signal y = ( y 1 , y 2 , y n ) , y n , the basic assumption in SR theory is that y can be approximately represented as a linear combination of a set of base signals { d i } i = 1 m from a redundant dictionary D n × m (n < m), the signal y can be expressed as
y = i = 1 m d i α i = D α
where α = ( α 1 , α 2 , α n ) is the unknown sparse coefficient vector, di is an atom of D . Given a signal y and an overcomplete dictionary D , the process of finding the sparse coefficient α is called sparse representation. The goal of SR is to calculate the sparsest α , which contains the fewest non-zero entries. However, this is an underdetermined problem due to the existence of overcomplete D . Generally, the sparsest α can be obtained by solving the following sparse model:
α ^ = arg min α α 0 ,   s . t . y D α 2 2 < ε
where 0 denotes the l 0 norm that counts the number of non-zero entries, ε > 0 is an error tolerance. Equation (2) is a non-deterministic polynomial-hard (NP-hard) problem which can be solved by the greedy approximation approach OMP [38].

2.2. Proposed Dictionary Learning Approach

For SR based medical image fusion, the ability of information expression of dictionary has a direct, great impact on the fusion results. Numerous studies [32,33,35] indicate that the dictionary obtained by the traditional training methods cannot generate impressive fusion results because of the weak expression ability. In order to construct an informative as well as compact dictionary, we propose a new dictionary learning method illustrated in Figure 1. In the first stage, two effective techniques including detail enhancement and multi-scale sampling are developed for improving the qualities of image patches. In the second stage, a neighborhood energy criterion is introduced to cluster the brightness patches and a multi-scale spatial frequency criterion is proposed to cluster the edge detail patches, and then generate two high-quality training sets. In the last stage, the brightness sub-dictionary and detail sub-dictionary can be constructed by training the two categories of patches using K-SDV, then the overcomplete dictionary is obtained by the combination of the sub-dictionaries. The detailed analyses of these three sub-sections are shown in the rest of the sub-sections.

2.2.1. Detail Enhancement

We employ the way of online learning to obtain the overcomplete dictionary. The advantage is that it can directly utilize useful information from source images and can be beneficial to enhance the expressive ability of the dictionary. However, the general way of dividing the image into a series of overlapping patches is directly performed on the pre-training images, which would lead to some un-conspicuous problems. When the details of certain regions in the source image are relatively weak, some important information in this region could not be well expressed in dictionary training. To address this problem, we develop a detail enhancement technique shown in Figure 2. The high-pass detail information with different levels are first extracted by simple multiple filtering, and then superimposes into the source images, which can be mathematically expressed as:
X ˜ l ( i , j ) = X l ( i , j ) + p = 1 P H X l , p ( i , j )
H X l , p = N D X l
where X l = X 1 ,   X 2 , denote the l-th source image, if only two source images in medical image fusion procedure, l = 1 , 2 , X ˜ l is the detail enhanced version of X l , H X l , p is the p-th level highpass detail information image of X l .
ND denotes the Neighbor Distance filter, some studies [39,40] demonstrated that ND can effectively extract and express the high frequency information of image, so we use ND to extract details in this step, more detailed introduction about ND can be found in [40]. In addition, it should be noted that this paper enhances the source image when training them to generate a strong overcomplete dictionary, rather than enhanced the source images during the fusion process.
As can be seen from Figure 2, by performing the process of detail enhancement, the weak information such as boundary, brightness and texture et al. of the source images are significantly enhanced. For the convenience of comparison, a same region from source image and enhanced image is selected by the red rectangle and enlarged and placed in the lower right of them. The enhanced image can highlight the unobvious but significant information, and express them into patches, which could be beneficial to constructing a strong training set.

2.2.2. Multi-Scale Sampling (MSS)

For dictionary learning, the traditional way of constructing a training set is first dividing the images into a series of overlap patches by a fixed size, and then pull them into a vector respectively. However, the fixed size of image patch usually cannot comprehensive describe image information, especially small target information in medical images. Generally, a large size of image patch is more likely to improve robustness and it can illustrate the shape and location of object, but the description about the details are insufficient, and vice versa. Therefore, choosing a suitable patch size is significant for improving the informative ability of the training set. To integrate the advantages of each patch with different sizes, we present a Multi-scale Sampling (MSS) scheme to obtain a stronger training set. MSS can not only improve multi-scale properties of the training set but also reduce the dictionary training time. The MSS can be obtained by the follows
X ˜ d , l = downsample ( X ˜ l , d )
where X ˜ d , l is the sampling version of X ˜ l by downsampling rate of d ( d = 1 , 2 , ), obviously, X ˜ l = X ˜ 1 , l . As can be seen from Figure 3. A series of different images with various smaller size of source can be generated.
To facilitate the explanation, the white circular areas in each of X ˜ d , l images in Figure 3 are marked by same-sized red rectangles. We can notice that the contour and details for the same region in diverse images with different size are significantly different. By performing MSS, the useful information can be comprehensively utilized from coarse to fine with a different scale, then the informative image patches can be obtained by dividing.
More importantly, MSS can increase the computation efficiency due to the downsampling operation. For example, as illustrated in Figure 3, the size of the input image X ˜ l is 256 × 256, when the downsampling rate d is set to 2, 3, 4 and 5, respectively, we can generate X ˜ 2 , l , X ˜ 3 , l , X ˜ 4 , l and X ˜ 5 , l , the size of each of them is smaller than X ˜ l . Using several images processed by MSS with different sampling rates instead of source images to construct a training set would not decrease dictionary quality, while reducing dictionary learning time. In order to analyze quantitatively, as shown in Table 1, we calculate the number of patches from the upper left to the lower right, which is produced by the patch size 8 × 8 and a step length of two (pixels overlapping is six). As can be seen from Table 1, the summation of patches obtained by the four images { X ˜ d , l } ( d = 2 , 3 , 4 ,   5 ) is much smaller than X ˜ l , which means the process of MSS can significantly reduce the training data, and then decrease the computation complexity in dictionary learning.

2.2.3. Sum of Neighborhood Energy (SNE)

Brightness information and detailed information are two important manifestations for medical image features. Thus, detecting the brightness and detail feature from the source images into the fused image is necessary to produce encouraging fusion results. Considering the energy correlation between pixels in medical images, the sum of neighborhood energy (SNE) is developed to evaluate the brightness of the central pixel, as follows:
e ( i , j ) = m = ( M 1 ) / 2 ( M 1 ) / 2 n = ( N 1 ) / 2 ( N 1 ) / 2 X ˜ d , l ( j + m , j + n ) 2
where e ( i , j ) denote the SNE of X ˜ d , l centered at pixel ( i , j ) ; X ˜ d , l ( i + m , j + n ) is the pixel value of X ˜ d , l at pixel ( i + m , j + n ) ; M × N indicate the size of neighborhood. The MR-T1/MR-T2 images and their corresponding brightness images are displayed in Figure 4.

2.2.4. Multi-Scale Spatial Frequency (MSF)

In medical images, the spatial frequency (SF) of the centered pixel can reflect the difference between the pixel and the surrounding pixel. Thus, SF of an image can be used to express the richness degree of the image detail information. Considering the correlation between pixels, a neighborhood based SF is usually employed, which is defined as:
S P X ˜ d , l ( i , j ) = m = ( M 1 ) / 2 m = ( M 1 ) / 2 n = ( N 1 ) / 2 n = ( N 1 ) / 2 [ X ˜ d , l ( i + m , j + n ) X ( i , j ) ]
where S P X ˜ d , l ( i , j ) is the value of SF of X ˜ d , l centered at ( i , j ) , M × N is the size of neighborhood.
In practical applications, to reduce the computational complexity, the size of the neighborhood is usually set to M = N. However, this technique in Equation (7) is sensitive to the neighborhood size, especially for the dark and low-contrast regions. To address this problem, we propose a novel detail detection technique called Multi-scale spatial frequency (MSF), The MSF of an image is mathematically defined as
M S F X ˜ d , l ( i , j ) = ε 1 | S P X ˜ d , l r 1 ( i , j ) S P X ˜ d , l r 2 ( i , j ) | + ε 2 | S P X ˜ d , l r 1 ( i , j ) S P X ˜ d , l r 3 ( i , j ) |
here
S P X ˜ d , l r ( i , j ) = m = ( r 1 ) / 2 m = ( r 1 ) / 2 n = ( r 1 ) / 2 n = ( r 1 ) / 2 [ X ˜ d , l ( i + m , j + n ) X ˜ d , l ( i , j ) ]
where, r ( r = r 1 , r 2 , r 3 ) is the scale factor, r × r denotes the size of neighborhood, ε 1 and ε 2 are weighting parameter given by the user. S P X ˜ d , l r ( i , j ) is the r-scale SF of X ˜ d , l at pixel ( i , j ) , M S F X ˜ d , l ( i , j ) is the MSF of X ˜ d , l at pixel ( i , j ) .

2.2.5. Clustering and Dictionary Learning

We classify all patches into two groups by taking SNE and MSF as the feature metrics for clustering. The patch clustering and dictionary learning consist of the following five steps.
• Step 1: Patches Collection
Adopt the slide window technique to divide each of X ˜ d , l into a suite of image patches with size ρ × ρ from upper left to lower right, and then overlapping each of them, the step length between two adjacent patches is δ pixels. The patches extraction results are as follows:
X ˜ d , l = X ˜ d , l 1 X ˜ d , l 2 X ˜ d , l k X ˜ d , l K
where denotes merging each of them, X ˜ d , l k is the k-th patch of X ˜ d , l , K is the total number of image patches.
• Step 2: Calculate patch feature by SNE and MSF
Calculate SNE and MSF of each image patch respectively, a series of brightness dominant patches and a number of detail dominant patches are obtained. According to Equation (6), the energy feature of patch X ˜ d , l k can be calculated as following:
E X ˜ d , l k = ( x , y ) X ˜ d , l k e X ˜ d , l k ( x , y )
where e X ˜ d , l k ( x , y ) is the brightness level of the k-th patch in X ˜ d , l at pixel ( x , y ) , E X ˜ d , l k denotes the sum of brightness information of patch X ˜ d , l k . Similarly, the sum of detail information contained in patch X ˜ d , l k can be generated by Equation (8):
C X ˜ d , l k = ( x , y ) X ˜ d , l k c X ˜ d , l k ( x , y )
where c X ˜ d , l k ( x , y ) denotes the detail level of the k-th patch in X ˜ d , l at pixel ( x , y ) , C X ˜ d , l k indicates the sum of the detail information of patch X ˜ d , l k
• Step 3: Construct training sets
Rearrange each patch X ˜ d , l k into a column vector X ¯ d , l k , and then construct a brightness training set by selecting the corresponding image patches containing the best brightness activity level by the following rules:
T e = { T ¯ 1 , T ¯ 2 , T ¯ k , , T ¯ K }
T ¯ k = max { E X ˜ d , 1 k , E X ˜ d , 2 k , E X ˜ d , l k , }
where T e is the brightness dominant training set, which composes the image patches with a better brightness feature, and T ¯ k denotes that the k-th patch from { E X ˜ d , 1 k , E X ˜ d , 2 k , E X ˜ d , l k , } has the largest SNE value among them. Similarly, construct a detail training set by choosing the corresponding image patches containing the richest detail activity level by the following rules:
W c = { W ¯ 1 , W ¯ 1 , W ¯ k , , W ¯ K }
W ¯ k = max { C X ˜ d , 1 k , C X ˜ d , 2 k , C X ˜ d , l k , }
where W c is the detail dominant training set which composes of image patches with better detail feature, and W ¯ k means that the k-th patch from { C X ˜ d , 1 k , C X ˜ d , 2 k , C X ˜ d , l k , } has the largest MSF value is selected.
• Step 4: Sub-dictionaries learning
The brightness sub-dictionary D e can be generated by calculating the objective function in Equation (15) by using K-SVD [41],
{ D e , α e } = arg min D e , α e { T e D e α e F 2 + α e 1 }
where D e and α e are the brightness sub-dictionary and its sparse coefficient, respectively.
To avoid this, a same patch is used in the brightness set as well as the detail set. At the same time, we add the following constraint
E X ˜ d , l k C X ˜ d , l k 0
The detail sub-dictionary D c is obtained by performing K-SVD on detail training set W c as follows
{ D c , α c } = arg min D c , α c { W c D c α c F 2 + α c 1 }
where D c and α c are the detail sub-dictionary and its sparse coefficient, respectively.
• Step 5: Construct the final dictionary
The final informative, compact and discriminative overcomplete dictionary is obtained by the combination of the sub-dictionaries D e and D c in Equation (20)
D = ( D e , D c )
The schematic diagram of the proposed dictionary learning is illustrated in Figure 1. Moreover, the formalized mechanism of the dictionary learning algorithm is described in Algorithms 1. For simplicity, we take two source images as an example to discuss, and it can easily be extended to the case of more than two source images.
Algorithm 1 The proposed dictionary learning algorithm
Inputs: Two group of images { X ˜ d , l } ( l = 1 , 2 )
(1) Extract patches of { X ˜ d , 1 } and { X ˜ d , 2 } from upper left to lower right.
(2) Patch classification based on SNE and MSF in Equations (6) and (8), respectively.
(3) Construct two training sets T e and W c by Equations (13) and (15), respectively.
(4) Obtain the two sub-dictionaries D e and D c by solve Equations (17) and (19), respectively.
(5) Generate the final dictionary D by Equation (20)
Output: The overcomplete dictionary D

2.3. Image Fusion

The whole image fusion schematic diagram is shown in Figure 5. In the sparse coding process, each source image X l is first divided into a series of image patches of size ρ × ρ from upper left to lower right with an overlapping step length of δ pixels, and then rearranging them into column vectors X ¯ l k ( k = 1 , 2 , , K ,   l = 1 , 2 , , n ), where K and n are the number of patches and source images, respectively. X ¯ l k denotes the k-th patch’ column vector of image X l . Calculate the sparse coefficient vectors { α ¯ 1 k , α ¯ 2 k , , α ¯ n k } of source images by using OMP algorithm and the following underdetermined equation can be solved
α ^ = arg min α { y D α 2 2 + α 1 }
where α ^ = { α ¯ 1 k , α ¯ 2 k , , α ¯ n k } , y is the sparse coefficient vectors { X ¯ 1 k , X ¯ 2 k , , X ¯ n k } . By employing the ‘‘max-absolute choosing’’ rule to obtain the fused sparse vector α ¯ F k
α ¯ F k = max { α ¯ 1 k , α ¯ 2 k , , α ¯ n k }
The k-th patch’ column vector of fused image can be generated by combing D and α ¯ F k
F ¯ k = D α ¯ F k
Reshape all of the column vectors F ¯ k into image patches F k and then plug each of F k into its corresponding original position, the final fused image F can be obtained:
F = k = 1 K F k

3. Experiments and Analysis

In this section, we present four subsections to analyze and verify the superiority of the proposed medical image fusion scheme. We first illustrate the detail experimental setups including test methods and parameters setting in Section 3.1. Then, the test images in experiments are displayed and briefly introduced in Section 3.2. We evaluate the fusion results of different methods based on their visual performance and quantitative assessment in Section 3.3 and Section 3.4. The computation efficiency analysis is discussed in Section 3.5. Finally, we extend our method into other type of image fusion applications in Section 3.6.

3.1. Test methods and Parameters Setting

In order to verify the effectiveness of the proposed method, six representative medical image fusion algorithms are selected to compare with the proposed algorithm. The first three compared methods are MST based: DTCWT, CVT and NSCT. The others are three state-of-the-art medical image fusion methods including the adaptive SR (Liu-ASR) proposed by Liu et al. [35], the Kim [33] method, which is based on clustering and principal component analysis methods, and the K-SVD based fusion technique developed by Zhu et al. [32]. The latter three carefully design algorithms and mainly paid their attention to the novel ways for training an informative overcomplete dictionary, they can produce the representative fusion performance for the field of medical image fusion in a few years.
For the parameters setting, the first layer of DTCWT uses ‘LeGall 5-3’ filter, and the other layer employs ‘Qshift-06’ filter. The NSCT method employs the ‘pyrexc’ filter as a pyramid filter, and the ‘vk’ filter as the directional filter, which is decomposed into four layers, and the number of directions from coarse to fine is 4, 8, 8, and 16, respectively. These setting can generate the best fusion performance [3]. For all the MST based methods, they fuse high-frequency coefficients and low-frequency coefficients by “max-absolute” and “weighted averaging” scheme. Furthermore, a 3×3 consistency verification check [42] is adopted to fuse the high-frequency coefficients. The default parameters given by the respective authors are used for the Liu-ASR, Zhu and Kim methods.
For the proposed method, during the training dictionary stage, the size of patch is conventionally set to 8×8, and the overlapping pixels is 6. For the detail enhancement, the number of level P is set to 4 in Equation (3). For the MSS, the downsampling rate is set 3, 4, 5 and 6, respectively in Equation (5). As a consequence, the proposed dictionary learning based method can generate high fusion performance while keep lower computation efficiency. Meanwhile, due to the strong flexibility of multiscale property in MSS, the MSS scheme can also well address the case of small size source images. The size of the neighborhood of the SNE in Equation (6) is set to 11 × 11. For MSF, the three scale factors r 1 , r 2 and r 3 are set to 7, 11 and 15, and the weight coefficients ε 1 and ε 2 are set to 0.33 and 0.67 in Equation (8), respectively. In the sparse representation stage, the reconstruction error is 0.1, and the number of iterations is 50. The setting of the parameters of the proposed algorithm can generate the best fusion results for the experimental images in this paper.

3.2. Test Images

For pixel level multi-modality medical image fusion, most of the test images generally can be obtained at http://www.imagefusion.org and http://www.med.harvard.edu/aanlib/home.html. To demonstrate the effectiveness of the proposed method, we utilize four pairs of real multimodal medical images shown in Figure 6 to evaluate the algorithm in our experiments, the size of all of them is 256 × 256, and all of them are perfectly registered, which means that the objects in each set of images are geometrically aligned. These medical images include MRI, CT, MR-T1, MR-T2, SPECT and PET, the characteristics of each of them are summarized as following.
  • CT has a shorter imaging time and a higher spatial resolution, whereas it provides soft tissue information with low contrast.
  • MRI can clearly display the soft tissue information of the human body, but it is hard to reflect the dynamic information of the metabolic activity in human body.
  • MR-T1 image is sensitive to observe the anatomy, while the MR-T2 image can detect the tissue lesions.
  • SPECT can show the biological activities of cells and molecules, but it is difficult to distinguish human organ tissues due to the low image quality of SPECT.
  • PET can reflect the metabolic activity information of human tissues and organs at the molecular level, but the spatial resolution of PET is relatively low.
In our experiments, all of the tests are implemented in Matlab 2016a using a desktop with a Core (TM) i7-3770 CPU, 3.40GHz, 8 CPUs and 12G RAM.

3.3. Fusion Results by Subjective Visual Effects Analysis

The results of the first set of experiments are shown in Figure 7, and the source MRI and CT images are listed in Figure 7a,b. The fused images obtained by the DCTWT, Curvelet, NSCT, Liu-ASR, Kim, Zhu, and the proposed method are displayed in Figure 7c–i, respectively. To facilitate subjective visual comparisons, local regions enclosed by red colored rectangular boxes in Figure 7 are enlarged and presented in the bottom right corners of their respective images. We can notice that the diverse fusion performance generated by different methods in retaining the brightness information and edge detail the information of source medical images. Despite the MST methods expressing source information with different scales and directions, all of the results obtained by DTCWT, Curvelet and NSCT in Figure 7c–e still produce some distortion, which heavily decreases the qualities of the medical images. Meanwhile, after careful observation, we can find the fusion results generated by Liu, Kim and Zhu can preserve details well from MRI. However, the brightness active level is relatively lower than the source CT image, which indicates the loss of useful information. It can be seen from Figure 7i that the fusion result obtained by the proposed method has the best performance in terms of retaining brightness information and detail information, which indicates that the visual effect of the method is the best.
The fusion results of different methods about “MR-T1/MR-T2” (see Figure 8a,b) are shown in Figure 8c–i. For this experiment, the detail information (shape, edge, texture et. al) in the red rectangle is mainly from the MR-T1 image while the energy information (brightness, contrast et.al) mainly from MR-T2 image. Among these fusion results, it can be seen that the details of the fusion results produced by the DTCWT, Curvelet and NSCT methods are severely damaged, especially the Curvelet method. Although the Liu-ASR, Kim and Zhu methods can relatively effectively protect the edge details of the source image, they do not protect the contrast of the image well. This is very disadvantageous for medical images with high-quality requirements, and is not conducive to subsequent medical image processing and recognition tasks.
By comparison, our fusion result (Figure 8i) can not only effectively protect the edge detail information of the source image, but also maintain the contrast of the source image, which is mainly due to the detail enhancement processing of training set and the clustering techniques to classify brightness and detail groups before dictionary learning. At the same time, little artificial false information is introduced in our fusion result, which means that the visual effect of the proposed method outperforms others methods in this experiment.
To further demonstrate that the proposed method is equally effective for other types of medical image fusion, we test the other two categories of medical image fusion systems include MRI/SPECT and MRI/PET, and the fusion results are shown in Figure 9 and Figure 10. The details of fused images are zoomed in and presented in the bottom left corners of their respective images.
With more careful observation, we can see that the spatial edge details and brightness information in Figure 9i and Figure 10i are more accurate compared with Figure 9c–h and Figure 10c–h generated by the six compared methods. This means the useful information from source images have been successfully transformed into the fused images, that is, our fused results have the best visual features. This is mainly because the novel construction of the overcomplete dictionary in this paper is composed of two parts include the brightness sub-dictionary and the edge detail sub-dictionary, which can fully express the significant features of the medical images. This will be very beneficial to the implementation of medical image fusion in practical medical assistance applications. In conclusion, the proposed algorithm has the best subjective visual effect.

3.4. Fusion Results by Objective Evaluation

To effectively evaluate the quality of fused images, we need take a quantitative evaluation for fused results, the metric can assess the perceptual quality of fused results consistently with visual subjective evaluation is highly desired. For medical image fusion problems, since the reference image (ground truth) does not exist in practice, quantitative evaluation of the quality of the fused image is not an effortless task. Over the last few decades, a number of fusions metrics have been proposed. Unfortunately, due to the various fusion scenarios, none of them are generally considered to be always more reasonable than other metrics. Therefore, it is often necessary to apply multiple indicators for comprehensive evaluation.
In this part, we employ five popular and representative quantitative indices to objectively evaluate all fusion results. These metrics include mutual information (MI) [43], edge information retention operator QAB/F [44], nonlinear correlation information entropy (QNCIE) [45], image fusion metric-based on phase congruency (QP) [46], and Piella’s Metric QS [47]. Among them, MI can reflect the degree of mutual information between the source images and the fused image; QAB/F can effectively evaluate the amount of the edge information transferred from the source images into the fused image; QNCIE is an information-based metric that measures the Correlation between fused image and the source images; QP indicates the degree of retaining significant saliency feature from the source images into the fused image. QS measures the contrast in the neighborhood of pixel.
For all of them, a larger value indicates better fusion results, and the codes of metrics include QNCIE, QP and QS are implemented in an evaluation toolbox by Liu [48]. The quantitative evaluation results of all methods are presented in Table 2 and Table 3. By comparison, we can see that the Zhu’s method will marginally outperform our method for the QP evaluation value in MRI/PET images. But our method can obtain the best performance for other metrics, and the subjective visual effect of our result (Figure 10i) is slightly superior to it (Figure 10h). As a whole, Table 2 and Table 3 clearly indicate our method outperforms the other methods in terms of almost metrics.
Consideration of the subjective visual comparison and objective quantitative evaluation, one can finally conclude that our method can generate visually pleasant fused images contain abundant detail and brightness information in most cases and outperforms some competing methods on visual quality and objective evaluation.

3.5. Computational Efficiency Analysis

As previously mentioned, the proposed dictionary learning method can address the defects of traditional dictionary learning algorithms with some superfluous patches and low computation efficiency. In this subsection, the computational complexity of different fusion methods are analyzed, the results are shown in Table 4. It is well known that SR based fusion method are generally time consuming, especially for the case of learning a dictionary from source images. As can be seen from Table 4, we can see that the first three MST based methods are less time-consuming compared with other four SR based methods.
Nonetheless, compared with the three typical dictionary methods include Liu-ASR, Kim and Zhu methods, the proposed method saves more time. This demonstrates that our method could increase the computational efficiency more than some dictionary learning based fusion methods. Furthermore, the computational efficiency of our method still has much room for improvement, and it can be further improved if we fully optimize the speed of its implementation and by utilizing multithreading and performance boosters like Graphics Processing Unit (GPU) acceleration.

3.6. Extension to Other Type Image Fusion Issues

To exhibit the generalization ability of the proposed method, we extend its application to other types of image fusion including multifocus image fusion, panchromatic-multispectral image fusion and infrared-visible image fusion. The three fusion examples are shown in Figure 11, where the source images are displayed in the first two columns, and the fused results are listed in the last column. We can notice that the important feature include detail, sharp, edge and brightness feature et al. of source images are well preserved in the fused image. Among them, the fusion performance of infrared-visible image fusion is relatively high. This demonstrates that the proposed model can transfer useful information from the source images into the fused result [49]. Meanwhile, fewer undesirable artifacts are introduced in these fusion process, which indicate our method has strong robustness in these applications. Furthermore, compared with MRI and PET imaging, some new functional medical imaging modality, such as photoacoustic computed tomography (PACT) [50], resting-state functional connectivity (RSFC) and functional connectivity photoacoustic tomography (fcPAT) [51] etc., can offer better spatial resolutions with a fast, noninvasive and non-ionizing imaging modality to describe the brain physiology and pathology. Therefore, further study will be required to investigate them for image fusion in the near future.

4. Conclusions

SR based dictionary learning technology has been widely used in the medical image fusion field due to its superior performance. The core problem of these kinds of algorithms is constructing an informative and compact overcomplete dictionary with abundant information. Aiming at the traditional dictionary learning methods lack enough ability to express the source image information. This paper proposes a novel dictionary learning method based on brightness and detail clustering for medical image fusion. The proposed approach consists of three steps. Firstly, multi-layer ND filtering enhances the details of pre-training images, thus the weak information can be reinforced in the training set.
At the same time, we conduct the MSS on images to realize the multi-scale presentation of patches. Secondly, we propose SNE and MSF to classify the brightness and detail information patches groups, and then construct the brightness sub-dictionary and the detail sub-dictionary by K-SVD. The combination of the two sub-dictionaries can generate the final informative and compact dictionary. Finally, a SR model is established, and generates the encouraging fused results. Experimental analysis on some traditional as well the state-of-the art dictionary learning based medical fusion methods on four categories of medical images fusion shows that the proposed method has obvious superiority in both subjective visual effect and quantitative evaluation.

Author Contributions

Conceptualization and methodology, F.Z.; software, Y.C.; writing original draft preparation, F.Z. and X.L.; review and editing, X.L.; visualization, M.Z.; supervision, H.T.

Funding

This research was supported by National Natural Science Foundation of China (No. 61471123).

Acknowledgments

The authors sincerely want to express heartfelt thanks to the dear editor and anonymous reviewers for the valuable comments and suggestions. We thank the supports by National Natural Science Foundation of China (No. 61471123). And we would like to thank Hengjun Zhao, Yu Liu, Zhiqin Zhu, Kim and Zheng Liu for providing the experimental codes for their methods.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, J.; Li, W.; Lu, K.; Xiao, B. An overview of multi-modal medical image fusion. Neurocomputing 2016, 215, 3–20. [Google Scholar] [CrossRef]
  2. Li, H.; He, X.; Tao, D.; Tang, Y.; Wang, R. Joint medical image fusion, denoising and enhancement via discriminative low-rank sparse dictionaries learning. Pattern Recognit. 2018, 79, 130–146. [Google Scholar] [CrossRef]
  3. Li, S.; Yang, B.; Hu, J. Performance comparison of different multi-resolution transforms for image fusion. Inf. Fusion 2011, 12, 74–84. [Google Scholar] [CrossRef]
  4. Li, H.; Li, X.; Yu, Z.; Mao, C. Multifocus image fusion by combining with mixed-order structure tensors and multiscale neighborhood. Inf. Sci. 2016, 349, 25–49. [Google Scholar] [CrossRef]
  5. Goshtasby, A.; Nikolov, S. Image fusion: Advances in the state of the art. Inf. Fusion 2007, 8, 114–118. [Google Scholar] [CrossRef]
  6. Du, J.; Li, W.; Xiao, B.; Nawaz, Q. Medical image fusion by combining parallel features on multi-scale local extrema scheme. Knowl. Based Syst. 2016, 113, 4–12. [Google Scholar] [CrossRef]
  7. Jiang, Q.; Jin, X.; Lee, S.; Yao, S. A novel multi-focus image fusion method based on stationary wavelet transform and local features of fuzzy sets. IEEE Access 2017, 5, 20286–20302. [Google Scholar] [CrossRef]
  8. Lewis, J.; Callaghan, R.; Nikolov, S.; Bull, D.; Canagarajah, N. Pixel and region based image fusion with complex wavelets. Inf. Fusion 2007, 8, 119–130. [Google Scholar] [CrossRef]
  9. Nencini, F.; Garzelli, A.; Baronti, S.; Alparone, L. Remote sensing image fusion using the curvelet transform. Inf. Fusion 2007, 8, 143–156. [Google Scholar] [CrossRef]
  10. Do, M.N.; Vetterli, M. The Contourlet Transform: An Efficient Directional Multiresolution Image Representation. IEEE Trans. Image Process. 2005, 14, 2091–2106. [Google Scholar] [CrossRef]
  11. Li, X.; Li, H.; Yu, Z.; Kong, Y. Multifocus image fusion scheme based on the multiscale curvature in nonsubsampled contourlet transform domain. Opt. Eng. 2015, 54, 073115. [Google Scholar] [CrossRef]
  12. Liu, Y.; Liu, S.; Wang, Z. A general framework for image fusion based on multi-scale transform and sparse representation. Inf. Fusion 2015, 24, 147–164. [Google Scholar] [CrossRef]
  13. Palsson, F.; Sveinsson, J.R.; Ulfarsson, M.O. Multispectral and Hyperspectral Image Fusion Using a 3-D-Convolutional Neural Network. IEEE Geosci. Remote Sens. Lett. 2017, 14, 639–643. [Google Scholar] [CrossRef]
  14. Liu, Y.; Chen, X.; Wang, Z.; Wang, Z.J.; Ward, R.K.; Wang, X. Deep learning for pixel-level image fusion: Recent advances and future prospects. Inf. Fusion 2018, 42, 158–173. [Google Scholar] [CrossRef]
  15. Li, S.; Kang, X.; Fang, L.; Hu, J.; Yin, H. Pixel-level image fusion: A survey of the state of the art. Inf. Fusion 2017, 33, 100–112. [Google Scholar] [CrossRef]
  16. He, C.; Liu, Q.; Li, H.; Wang, H. Multimodal medical image fusion based on IHS and PCA. Procedia Eng. 2010, 7, 280–285. [Google Scholar] [CrossRef] [Green Version]
  17. Jiang, Y.; Wang, M. Image fusion with morphological component analysis. Inf. Fusion 2014, 18, 107–118. [Google Scholar] [CrossRef]
  18. Zhang, Q.; Liu, Y.; Blum, R.; Han, J.; Tao, D. Sparse representation based multi-sensor image fusion for multi-focus and multi-modality images: A review. Inf. Fusion 2018, 40, 57–75. [Google Scholar] [CrossRef]
  19. Shang, L.; Liu, S.; Zhou, Y.; Sun, Z. Modified sparse representation based image super-resolution reconstruction method. Neurocomputing 2017, 228, 37–52. [Google Scholar] [CrossRef]
  20. Gu, S.; Zuo, W.; Xie, Q.; Meng, D.; Feng, X.; Zhang, L. Convolutional Sparse Coding for Image Super-resolution. In Proceedings of the IEEE International Conference on Computer Vision, Santiago, Chile, 13–16 December 2015; pp. 1823–1831. [Google Scholar]
  21. Liu, H.; Liu, Y.; Sun, F. Robust exemplar extraction using structured sparse coding. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 1816–1821. [Google Scholar] [CrossRef]
  22. Xu, S.; Yang, X.; Jiang, S. A fast nonlocally centralized sparse representation algorithm for image denoising. Signal Process. 2017, 131, 99–112. [Google Scholar] [CrossRef]
  23. Mourabit, I.; Rhabi, M.; Hakim, A.; Laghrib, A.; Moreau, E. A new denoising model for multi-frame super-resolution image reconstruction. Signal Process. 2017, 132, 51–65. [Google Scholar] [CrossRef]
  24. Karanam, S.; Li, Y.; Radke, R. Person re-identification with discriminatively trained viewpoint invariant dictionaries. In Proceedings of the IEEE International Conference on Computer Vision, Santiago, Chile, 13–16 December 2015; pp. 4516–4524. [Google Scholar]
  25. An, L.; Chen, X.; Yang, S.; Bhanu, B. Sparse representation matching for person re-identification. Inf. Sci. 2016, 355, 74–89. [Google Scholar] [CrossRef] [Green Version]
  26. Bahrampour, S.; Nasrabadi, N.; Ray, A.; Jenkins, W. Multimodal task-driven dictionary learning for image classification. IEEE Trans. Image Process. 2015, 25, 24–38. [Google Scholar] [CrossRef] [PubMed]
  27. Li, S.; Yin, H.; Fang, L. Group-sparse representation with dictionary learning for medical image denoising and fusion. IEEE Trans. Biomed. Eng. 2012, 59, 3450–3459. [Google Scholar] [CrossRef] [PubMed]
  28. Li, Y.; Sun, Y.; Huang, X.; Qi, G.; Zheng, M.; Zhu, Z. An Image Fusion Method Based on Sparse Representation and Sum Modified-Laplacian in NSCT Domain. Entropy 2018, 20, 522. [Google Scholar] [CrossRef]
  29. Wang, K.; Qi, G.; Zhu, Z.; Chai, Y. A Novel Geometric Dictionary Construction Approach for Sparse Representation Based Image Fusion. Entropy 2017, 19, 306. [Google Scholar] [CrossRef]
  30. Yang, B.; Li, S. Multifocus Image Fusion and Restoration with Sparse Representation. IEEE Trans. Instrum. Meas. 2010, 59, 884–892. [Google Scholar] [CrossRef]
  31. Yang, B.; Li, S. Pixel-level image fusion with simultaneous orthogonal matching pursuit. Inf. Fusion 2012, 13, 10–19. [Google Scholar] [CrossRef]
  32. Zhu, Z.; Chai, Y.; Yin, H.; Li, Y.; Liu, Z. A novel dictionary learning approach for multi-modality medical image fusion. Neurocomputing 2016, 214, 471–482. [Google Scholar] [CrossRef]
  33. Kim, M.; Han, D.K.; Ko, H. Joint patch clustering-based dictionary learning for multimodal image fusion. Inf. Fusion 2016, 27, 198–214. [Google Scholar] [CrossRef]
  34. Yin, H. Sparse representation with learned multiscale dictionary for image fusion. Neurocomputing 2015, 148, 600–610. [Google Scholar] [CrossRef]
  35. Liu, Y.; Wang, Z. Simultaneous image fusion and denoising with adaptive sparse representation. IET Image Process. 2015, 9, 347–357. [Google Scholar] [CrossRef]
  36. Qi, G.; Wang, J.; Zhang, Q.; Zeng, F.; Zhu, Z. An Integrated Dictionary-Learning Entropy-Based Medical Image Fusion Framework. Future Internet 2017, 9, 61. [Google Scholar] [CrossRef]
  37. Zhu, Z.; Yin, H.; Chai, Y.; Li, Y.; Qi, G. A novel multi-modality image fusion method based on image decomposition and sparse representation. Inf. Sci. 2018, 432, 516–529. [Google Scholar] [CrossRef]
  38. Pati, Y.C.; Rezaiifar, R.; Krishnaprasad, P.S. Orthogonal matching pursuit: Recursive function approximation with applications to wavelet decomposition. In Proceedings of the 27th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 1–3 November 1993; pp. 40–44. [Google Scholar]
  39. Li, H.; Liu, X.; Yu, Z.; Zhang, Y. Performance improvement scheme of multifocus image fusion derived by difference images. Signal Process. 2016, 128, 474–493. [Google Scholar] [CrossRef]
  40. Zhao, H.; Shang, Z.; Tang, Y.; Fan, B. Multi-focus image fusion based on the neighbor distance. Pattern Recognit. 2013, 46, 1002–1011. [Google Scholar] [CrossRef]
  41. Aharon, M.; Elad, M.; Bruckstein, A. K-SVD: An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation. IEEE Trans. Signal Process. 2006, 54, 4311–4322. [Google Scholar] [CrossRef]
  42. Li, H.; Manjunath, B.S.; Mitra, S.K. Multisensor image fusion using the wavelet transform. Graph. Models Image Process. 1995, 57, 235–345. [Google Scholar] [CrossRef]
  43. Qu, G.; Zhang, D.; Yan, P. Information measure for performance of image fusion. Electron. Lett. 2002, 38, 313–315. [Google Scholar] [CrossRef]
  44. Xydeas, C.S.; Petrovic, V. Objective image fusion performance measure. Electron. Lett. 2000, 36, 308–309. [Google Scholar] [CrossRef]
  45. Wang, Q.; Shen, Y.; Jin, J. Performance Evaluation of Image Fusion Techniques. In Image Fusion: Algorithms and Applications, 1st ed.; Stathaki, T., Ed.; Elsevier: London, UK, 2008; Volume 19, pp. 469–492. [Google Scholar]
  46. Zhao, J.; Laganiere, R.; Liu, Z. Performance assessment of combinative pixel-level image fusion based on an absolute feature measurement. IJICIC 2007, 3, 1433–1447. [Google Scholar]
  47. Piella, G.; Heijmans, H. A New Quality Metric for Image Fusion. In Proceedings of the International Conference on Image Processing, Barcelona, Spain, 14–17 September 2003; pp. 173–176. [Google Scholar]
  48. Liu, Z.; Blasch, E.; Xue, Z.; Zhao, J.; Laganière, R.; Wu, W. Objective assessment of multiresolution image fusion algorithms for context enhancement in night vision: A comparative study. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 34, 94–108. [Google Scholar] [CrossRef] [PubMed]
  49. Li, H.; Qiu, H.; Yu, Z.; Zhang, Y. Infrared and visible image fusion scheme based on NSCT and low-level visual features. Infrared Phys. Technol. 2016, 76, 174–184. [Google Scholar] [CrossRef]
  50. Zhao, J.; Xia, J.; Maslov, K.I.; Nasiriavanaki, M.; Tsytsarev, V.; Demchenko, A.V.; Wang, L.V. Noninvasive photoacoustic computed tomography of mouse brain metabolism in vivo. NeuroImage 2013, 64, 257–266. [Google Scholar] [Green Version]
  51. Nasiriavanaki, M.; Xia, J.; Wan, H.; Bauer, A.Q.; Culver, J.P.; Wang, L.V. High-resolution photoacoustic tomography of resting-state functional connectivity in the mouse brain. Proc. Natl. Acad. Sci. USA 2013, 111, 21–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Procedure of dictionary learning.
Figure 1. Procedure of dictionary learning.
Entropy 21 00267 g001
Figure 2. Enhancement of training samples.
Figure 2. Enhancement of training samples.
Entropy 21 00267 g002
Figure 3. Sampling.
Figure 3. Sampling.
Entropy 21 00267 g003
Figure 4. Images (left) and their brightness (right) by SNE.
Figure 4. Images (left) and their brightness (right) by SNE.
Entropy 21 00267 g004
Figure 5. Framework of the proposed image fusion.
Figure 5. Framework of the proposed image fusion.
Entropy 21 00267 g005
Figure 6. Pairs of medical source images: (a,b) MRI/CT images, (c,d) MR-T1/MR-T2 images; (e,f) MRI/SPECT images; (g,h) MRI/PET images.
Figure 6. Pairs of medical source images: (a,b) MRI/CT images, (c,d) MR-T1/MR-T2 images; (e,f) MRI/SPECT images; (g,h) MRI/PET images.
Entropy 21 00267 g006
Figure 7. Fused results of MRI/CT medical images by different methods. (a) MRI; (b) CT; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Figure 7. Fused results of MRI/CT medical images by different methods. (a) MRI; (b) CT; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Entropy 21 00267 g007
Figure 8. Fused results of MR-T1/MR-T2 medical images by different methods. (a) MR-T1; (b) MR-T2; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Figure 8. Fused results of MR-T1/MR-T2 medical images by different methods. (a) MR-T1; (b) MR-T2; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Entropy 21 00267 g008
Figure 9. Fused results of MRI/SPECT medical images by different methods. (a) MRI; (b) SPECT; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Figure 9. Fused results of MRI/SPECT medical images by different methods. (a) MRI; (b) SPECT; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Entropy 21 00267 g009aEntropy 21 00267 g009b
Figure 10. Fused results of MRI/PET medical images by different methods. (a) MRI; (b) PET; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Figure 10. Fused results of MRI/PET medical images by different methods. (a) MRI; (b) PET; (c) DCTWT; (d) Curvelet; (e) NSCT; (f) Liu-ASR; (g) Kim; (h) Zhu; (i) Proposed.
Entropy 21 00267 g010
Figure 11. Other types of image fusion examples using the proposed method.
Figure 11. Other types of image fusion examples using the proposed method.
Entropy 21 00267 g011
Table 1. The number of patches by different size images.
Table 1. The number of patches by different size images.
X ˜ l X ˜ 2 , l X ˜ 3 , l X ˜ 4 , l X ˜ 5 , l
Number of Patches1562537211600841529
Summation156256691
Table 2. The quantitative assessment of different methods for MRI/CT and MR-T1/MR-T2 fusion.
Table 2. The quantitative assessment of different methods for MRI/CT and MR-T1/MR-T2 fusion.
ImagesMetricDctwtCurveletNSCTLiu-ASRKimZhuProposed
MRI/CTMI3.38783.15223.55823.78333.66423.66804.9366
QAB/F0.51340.45810.56350.54040.47430.46200.6254
QNICE0.80930.80860.80980.81050.81020.81020.8161
QP0.30940.24600.39480.43530.31680.34690.5421
QS0.72730.68810.75330.77200.76540.76310.7996
MR-T1/MR-T2MI2.97062.86203.05173.29293.18323.28114.7856
QAB/F0.53950.49760.56130.52690.50080.49550.6566
QNICE0.80740.80710.80760.80820.80790.80820.8158
QP0.45130.37090.46930.48960.40890.46140.6735
QS0.72900.68880.76590.81050.80420.80450.8532
The best scores of each metric are highlighted in bold faces.
Table 3. The quantitative assessment of different methods for MRI/SPECT and MRI/PET fusion.
Table 3. The quantitative assessment of different methods for MRI/SPECT and MRI/PET fusion.
ImagesMetricDctwtCurveletNSCTLiu-ASRKimZhuProposed
MRI/SPECTMI2.71042.64302.73472.82512.72892.82353.4636
QAB/F0.66460.65220.66350.63820.57120.61700.6829
QNICE0.80630.80620.80640.80660.80630.80660.8087
QP0.42830.38540.45400.50420.31230.41220.5317
QS0.88900.85130.90970.91430.89500.90530.9186
MRI/PETMI2.49272.46602.56222.70892.59342.75013.3282
QAB/F0.52080.50760.55850.58010.47460.52600.6038
QNICE0.80550.80540.80570.80600.80570.80610.8077
QP0.31980.28750.33660.43830.26610.36410.3519
QS0.76420.71580.82960.82420.80170.82580.8582
The best scores of each metric are highlighted in bold faces.
Table 4. Time consumption of the different fusion methods (Time/s).
Table 4. Time consumption of the different fusion methods (Time/s).
CurveletDCTWTNSCTLiu-ASRKimZhuProposed
MRI/CT16.157.8925.4286.6661.6755.5435.24
MR-T1/MR-215.867.7125.4880.5361.3840.7533.99
MRI/SPECT39.2714.1175.76173.2955.20226.6988.84
MRI/PET39.4414.3879.85178.6659.08227.6782.10

Share and Cite

MDPI and ACS Style

Zhou, F.; Li, X.; Zhou, M.; Chen, Y.; Tan, H. A New Dictionary Construction Based Multimodal Medical Image Fusion Framework. Entropy 2019, 21, 267. https://doi.org/10.3390/e21030267

AMA Style

Zhou F, Li X, Zhou M, Chen Y, Tan H. A New Dictionary Construction Based Multimodal Medical Image Fusion Framework. Entropy. 2019; 21(3):267. https://doi.org/10.3390/e21030267

Chicago/Turabian Style

Zhou, Fuqiang, Xiaosong Li, Mingxuan Zhou, Yuanze Chen, and Haishu Tan. 2019. "A New Dictionary Construction Based Multimodal Medical Image Fusion Framework" Entropy 21, no. 3: 267. https://doi.org/10.3390/e21030267

APA Style

Zhou, F., Li, X., Zhou, M., Chen, Y., & Tan, H. (2019). A New Dictionary Construction Based Multimodal Medical Image Fusion Framework. Entropy, 21(3), 267. https://doi.org/10.3390/e21030267

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop