Next Article in Journal
Deep Learning-Based Change Detection in Remote Sensing Images: A Review
Next Article in Special Issue
Stepwise Fusion of Hyperspectral, Multispectral and Panchromatic Images with Spectral Grouping Strategy: A Comparative Study Using GF5 and GF1 Images
Previous Article in Journal
The Impact of Vegetation on the Visibility of Archaeological Features in Airborne Laser Scanning Datasets from Different Acquisition Dates
Previous Article in Special Issue
A New Adaptive Remote Sensing Extraction Algorithm for Complex Muddy Coast Waterline
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multispectral and SAR Image Fusion Based on Laplacian Pyramid and Sparse Representation

1
School of Resource and Environmental Sciences, Wuhan University, Wuhan 430079, China
2
Collaborative Innovation Center of Geospatial, Wuhan University, Wuhan 430079, China
3
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(4), 870; https://doi.org/10.3390/rs14040870
Submission received: 14 December 2021 / Revised: 7 February 2022 / Accepted: 8 February 2022 / Published: 11 February 2022

Abstract

:
Complementary information from multi-sensors can be combined to improve the availability and reliability of stand-alone data. Typically, multispectral (MS) images contain plentiful spectral information of the Earth’s surface that is beneficial for identifying land cover types, while synthetic aperture radar (SAR) images can provide abundant information on the texture and structure of target objects. Therefore, this paper presents a fusion framework to integrate the information from MS and SAR images based on the Laplacian pyramid (LP) and sparse representation (SR) theory. LP is performed to decompose both the multispectral and SAR images into high-frequency components and low-frequency components, so that different processing strategies can be applied to multi-scale information. Low-frequency components are merged based on SR theory, whereas high-frequency components are combined based on a certain activity-level measurement, identifying salient features. Finally, LP reconstruction is performed to obtain the integrated image. We conduct experiments on several datasets to verify the effectiveness of the proposed method. Both visual interpretation and statistical analyses demonstrate that the proposed method strikes a satisfactory balance between spectral information preservation and the enhancement of spatial and textual characteristics. In addition, a further discussion regarding the adjustability property of the proposed method shows its flexibility for further application scenarios.

1. Introduction

Multi-sensor image fusion aims to combine complementary information from different sources to obtain improved results with greater quality and reliability [1,2]. Fusion of multispectral (MS) and synthetic aperture radar (SAR) images, as one important branch of multi-sensor image fusion, has been attracting increasing attention. Multispectral image reflects abundant information of spectral signature of the ground objects and thus is useful for distinguishing different land cover types, whereas it heavily depends on solar illumination and weather conditions. On the other hand, SAR operates as an all-weather, all-time earth observation system, capable of penetrating clouds, rain, smoke, and fog. Moreover, SAR image usually has good contrast and provides sufficient textural and structural information of observed objects. Therefore, the integration of multispectral and SAR images can effectively utilize the complementary information and achieve the full potential of both datasets, which can improve the accuracy and efficiency of remote sensing image interpretation and information extraction, as well as further applications [3,4,5,6].
Due to the differences in imaging mechanism and spectral characteristics, multi-source heterogeneous image fusion (e.g., optical and SAR image fusion) has always been a sophisticated and challenging topic. Therefore, the literature that focuses on pixel-level multispectral and SAR image fusion algorithms is limited and usually borrows the techniques utilized in the field of pansharpening. To our knowledge, most existing methods can be grouped into four categories: component substitution (CS)-based methods, multi-resolution analysis (MRA)-based methods, model-based methods, and hybrid methods [7]. In recent years, the deep learning-based method has gradually come into people’s vision [8,9,10]. It is immature at the current stage due to the lack of massive training datasets and the adaptability to various sensors [11,12,13], which is beyond the scope of this paper. The earliest methods are basically based on component substitution, which is the substitution of a component obtained from a certain transformation of MS bands with the high spatial resolution (HR) SAR image. Representative CS-based methods include intensity-hue-saturation (IHS) [14,15], principal component analysis (PCA) [16,17], Gram–Schmidt (GS) [18], and Brovey [19]. CS-based methods can remarkably enhance the spatial features of MS images but usually cause more severe spectral distortion than in the pansharpening task. This is mainly because the assumption that the substituted component and high spatial resolution image have a strong correlation may be unreasonable in multispectral and SAR image fusion. The MRA-based methods rely on the injection of high-frequency information that is extracted from the high spatial resolution SAR image into the MS image via MRA tools. Popular MRA-based methods comprise wavelet transform (WT) [20], high-pass filter (HPF) [21], à trous wavelet transform (AWT) [22] and Laplacian pyramid (LP) [23]. Since only the spatial contents of MS images are changed, the spectral information is well preserved. Thus, MRA-based methods generally perform better than CS-based methods in the conservation of spectral properties. However, MRA-based methods unavoidably lose spatial characteristics contained in the low frequency of source images [24]. In addition, MRA-based approaches are sensitive to aliasing and misregistration, which is vulnerable to blocking artifacts and blur effects [25]. Model-based methods are also extended for the integration of MS and SAR images, which are originally proposed for the pansharpening tasks [26,27]. These methods may cause less registration errors but usually lead to a higher computational complexity.
According to the aforementioned introduction of the three groups of methods, they all have their own advantages and disadvantages. Naturally, hybrid methods that incorporate the merits of three groups of methods were proposed. Alparone et al. [28] combined the generalized IHS transform and à trous wavelet decomposition to inject panchromatic and SAR features into MS images. Chibani [29] conjointly used the modified Brovey transform (MBT) and à trous wavelet decomposition to integrate SAR features into MS images. Hong et al. [30] developed a fusion method based on IHS transform and wavelet decomposition to fuse moderate spatial resolution MS imagery and high spatial resolution SAR image. Yin [31] proposed a fusion algorithm based on support value transform (SVT) and sparse representation for the purpose of target recognition. Shao et al. [32] combined the IHS transform and gradient transfer fusion (GTF) algorithm to maintain the spatial details from images of both MS and SAR. Generally, hybrid methods have become more and more popular because of their comprehensive abilities of dealing with different problems [7].
Since most of the previous studies simply exploited the elementary methods developed many years ago, complementary information contained in MS and SAR images can not be utilized simultaneously. Therefore, it is necessary to explore more advanced tools such as sparse representation (SR) to excavate the underlying information. Sparse representation theory has been successfully applied in various fields of computer vision and image processing [33,34,35,36]. Sparse representation theory assumes that natural signals, such as images, can be represented or approximately represented by a linear combination of a small number of atoms, which are columns of an over-complete dictionary [37]. Li et al. [38] first brought the SR theory into multi-focus image fusion in 2009. Subsequently, ever-growing interests can be seen in SR-based image fusion [39,40], especially in the remote sensing image fusion [41] and multi-modal image fusion [42,43]. The reason why SR theory works effectively in image fusion tasks is that sparse coefficients can well describe the salient features and structure information contained in images, which is immensely beneficial for image fusion [38]. Nevertheless, stand-alone SR-based image fusion approaches have some drawbacks [24]. Firstly, the sparse coding stage in SR-based image fusion usually consumes a considerable computational burden, and the running time will grow sharply as the size of source images increases. Secondly, the sliding window technique employed in sparse representation may cause smoothness and missing of details, especially when the overlapped regions between adjacent patches are large.
Therefore, in this paper, we proposes a fusion framework combining the merits of MRA-based and SR-based methods to integrate multispectral and SAR images. Specifically, Laplacian pyramid (LP) is chosen as the multi-scale transform tool in our study since images can be decomposed into uncorrelated high-frequency channels without losing details such as edges and textures [44]. Considering the first weakness, the running time can be sharply reduced, as the sparse coding operation is only implemented in small-sized low-pass bands. For the second demerit, most spatial information is separated by the LP and preserved in the high-frequency components via effective fusion rule. Furthermore, this paper also explores the flexibility and adjustability of the proposed fusion framework for potential applications.
The remainder of this paper is organized as follows. Section 2 describes the proposed multispectral and SAR image fusion framework in detail. The experimental results and analyses are provided in Section 3. Discussion about the adjustability of fusion results is included in Section 4. Finally, Section 5 concludes this paper.

2. Proposed Fusion Framework

The conventional Laplacian-pyramid-based image fusion methods generally employ the “averaging” rule to merge the low-frequency components, which usually leads to the loss of contrasts and details. On the other hand, SR-based image fusion approaches are limited by the high computation expenses and may be smoothed due to the sliding window technique. In order to overcome the disadvantages of stand-alone methods, we combine the Laplacian pyramid and sparse presentation theory to integrate multispectral and SAR images via the proposed fusion method.
In this paper, we adopt a band-wise fusion strategy for each channel of multispectral images, that is to say each channel follows the same fusion procedures to be integrated with SAR image. (1) LP generation. First, multispectral and SAR image LPs are constructed. Both the multispectral and SAR images are decomposed into high-frequency components and low-frequency components. (2) High-frequency components fusion. For high-frequency components, an appropriate fusion rule is employed to maximally reserve the high-frequency features, such as edges, point targets, and lines. (3) Low-frequency components fusion. Subsequently, sparse representation theory is applied to low-frequency components fusion. (4) Reconstruction. Finally, a reconstruction process is implemented using fused high-frequency and low-frequency components according to the inverse Laplacian pyramid transform. The proposed fusion method is outlined in Figure 1.

2.1. LP Generation and Reconstruction

The Laplacian pyramid can be seen as a hierarchical image structure that represents an image as a sequence of bandpass sub-images at different spatial scales. This versatile data structure makes LP useful for image fusion [45]. Some multi-scale transform tools, such as wavelet transform, are critically subsampled and thus unavoidably create artifacts due to misregistration. Compared with them, LP is oversampled, without detail decimation, and therefore it is considerably appropriate for our study.
After the original multispectral image I M S and despeckled SAR image I S A R are aligned geometrically, LP decomposition is conducted to both MS and SAR images, and then the corresponding LPs are generated. The implementation details can refer to [46]. The N-level LPs of MS and SAR images are denoted as { L l M S } l = 0 N and { L l S A R } l = 0 N , respectively. High-frequency information including textures and edges are maintained in a series of detail images { L l M S } l = 0 N 1 and { L l S A R } l = 0 N 1 at different spatial scales, while low-frequency information is stored in the approximation images L N M S and L N S A R at a much smaller scale compared to the original images.
Once the fused high-frequency { L l F } l = 0 N 1 and low-frequency L N F components are obtained, the reconstruction process is performed following inverse Laplacian pyramid transform, and then the fused image L N F is obtained. The reconstruction process can be expressed as:
I F = L 0 F + EXPAND ( L 1 F + EXPAND ( L 2 F + + EXPAND ( L N F ) ) )
where the function EXPAND consists of up-sampling and interpolating operations [46].

2.2. High-Frequency COMPONENTS Fusion

High-frequency components are also known as detail images in LP, which contain abundant details such as textures and edges. The absolute value of each pixel in the detail image is set as activity level measurement, which represents the sharpness or edges of an image. Accordingly, a popular “max-absolute” fusion rule [47] is exploited for merging the detail images { L l M S } l = 0 N 1 and { L l S A R } l = 0 N 1 . The fused high-frequency components { L l F } l = 0 N 1 are obtained according to Equation (2):
L l F ( i , j ) = { L l M S ( i , j ) ,   if   | L l M S ( i , j ) | | L l S A R ( i , j ) | L l S A R ( i , j ) ,   otherwise
where l = 0 , 1 , , N 1 , 0 i < R l , 0 j < C l , R l , and C l are the row and column number of L l M S , L l S A R , and L l F .

2.3. Low-Frequency Component Fusion

The low-frequency components in LP are also called the approximation images, which approximately represent the source image in a smaller scale. The traditional multi-resolution analysis (MRA) fusion methods often adopt a simple average rule when the low-frequency information is fused. Since most of the energy in an image is concentrated in the low-frequency information, the averaging operation inevitably reduces the contrast of the fusion results [24]. Consequently, in the proposed low-frequency component fusion, sparse representation theory is utilized to maintain the information that may be lost in the LP-based fusion framework. The SR-based fusion method can extract the local structure information contained in the low-frequency components thanks to the powerful image feature expression ability of the over-complete dictionary [38]. Therefore, the energy of original images is maintained as much as possible, and the contrast of the fused image is improved. The flowchart of the low-frequency component fusion, based on sparse representation, is illustrated in Figure 2.
For convenience, the low-frequency component L N M S of MS image Laplacian pyramid is denoted as A, the low-frequency component L N S A R of SAR image Laplacian pyramid is denoted as B, and the fused low-frequency component L N F of fusion image Laplacian pyramid is denoted as F. The size of image A, B, and F is M × N . The SR-based low-frequency components fusion in this paper takes the following steps:
(1)
Patch generation. To make full use of the local information of source images, the sliding window technique is applied to divide the source images A and B into image patches p k A , p k B   ( k = 1 , 2 , , K   and   K = ( M n + 1 ) × ( N n + 1 ) ) of size n × n , starting from the top-left to the bottom-right with a fixed step length s .
(2)
Vectorization. Then, image patches p k A , p k B are rearranged to n × 1 vectors v k A , v k B in a column-wise way. Each vector is normalized to zero-mean via subtracting the mean value according to the following Equation (3), and the mean values are stored for subsequent reconstruction process [47],
{ v ^ k A = v k A v ¯ k A 1 v ^ k B = v k B v ¯ k B 1
where v ¯ k A and v ¯ k B are the mean values of vectors v k A and v k B , respectively, and 1 is a n × 1 vector of all ones.
(3)
Sparse coding. Calculate the sparse coefficients of vectors v ^ k A and v ^ k B according to Equation (4) using the simultaneous orthogonal matching pursuit (SOMP) algorithm [48]. The SOMP algorithm is employed here for its high computing efficiency and suitability for image fusion,
{ α k A = arg min α α 0   s . t . v ^ k A D α 2 2 ε α k B = arg min α α 0   s . t . v ^ k B D α 2 2 ε
where α 0 denotes the number of nonzero elements in α , D refers to a pre-defined dictionary, and ε is the error tolerance.
(4)
Coefficient fusion. The activity level measurement and fusion rule are two important issues in image fusion tasks [47]. In this paper, the absolute value of the sparse coefficient is chosen to describe the activity level, and the popular max-absolute rule is selected as the fusion rule to combine the corresponding sparse coefficients. The detailed fusion process can be described by Equation (5):
α k F ( t ) = { α k A ( t ) ,   if   | α k A ( t ) | | α k B ( t ) | α k B ( t ) ,   otherwise
where α k A ( t ) and α k B ( t ) are the t-th element in α k A and α k B , respectively.
(5)
Vector reconstruction. The fused sparse vector v ^ k F is obtained via the fused sparse coefficient α k F multiplied by the same dictionary used in Step (3). The local mean subtracted in Step (2) is added back, and the final fused vector v k F is obtained.
v ^ k F = D α k F
v ¯ k F = 1 2 ( v ¯ k A + v ¯ k B )
v k F = v ^ k F + v ¯ k F 1
(6)
Final reconstruction. Every fused sparse vector v k F is reshaped to a n × n patch and placed in the corresponding position in the fused image F. As the patches may be overlapped, the same pixel in the source image may appear in multiple patches. In other words, one position in F may relate to multiple patches. Therefore, each pixel’s value in the fused image F is the average value of the corresponding elements in all related patches. Finally, the fused low-frequency component L N F is obtained.

3. Experiments

3.1. Experiment Settings

3.1.1. Data Description

In order to verify the validity and effectiveness of the proposed fusion framework, we selected three pairs of airborne multispectral images and X-band intensity images of airborne SAR AeS-1 [49,50] for experiments, covering the Trudering region of Munich in Germany. The original multispectral and SAR images have the same spatial resolution of 1.5 m and are co-registered. In this paper, multispectral images were spatially degraded to 6 m in order to verify the spatial enhancement of various algorithms. An advanced SAR despeckling algorithm [51] was selected to reduce the impact of speckle noise. Subsequently, multispectral images were resampled to the same size of SAR images.
As shown in Figure 3, area 1 is featured with agricultural fields and field ridges, area 2 is dominated by buildings and grasslands, and area 3 contains more diversified land cover types. The selected areas feature abundant land cover types and ground features, with different shapes, sizes, and orientations that are appropriate for evaluation of the fusion methods in terms of spectral preservation and spatial enhancement.

3.1.2. Evaluation Metrics

Both visual inspection and quantitative analysis are included to evaluate the performance of all the fusion methods. Several popular evaluation indicators are employed to quantitatively analyze the fusion results in terms of spectral preservation, spatial enhancement, and the amount of information. Correlation coefficient (CC) and spectral angle mapper (SAM) [52] are selected to evaluate the spectral fidelity between the spatially degraded HR multispectral and fused image. Spatial correlation coefficient (SCC) [53], structural similarity (SSIM) [54], and average gradient (AG) [50] are chosen as indicators measuring the spatial enhancement of various fusion approaches. Feature mutual information (FMI) [55] is used to assess the amount of information transferred from the source images. Detailed descriptions can be found in the literatures listed above. All the metrics are calculated as average values across multi bands, except that SAM is averaged over all the pixels.
In order to synthetically compare the statistical results, we also define a comprehensive index (CI) considering spectral preservation, spatial enhancement, and information incorporation simultaneously, which is formulated as:
CI = i ω i X i i ω i
where X i denotes the i-th evaluation metric, and ω i represents the corresponding weight of the evaluation metric.

3.1.3. Comparison Methods

Five typical fusion methods covering CS-based, MRA-based, and hybrid methods were chosen for comparison with the proposed Laplacian pyramid and sparse representation-based (LPSR) method, including IHS fusion method [14], BTH fusion method [56], Laplacian pyramid (LP)-based fusion method [23], AWLPR fusion method [57], and IHS+Wavelet (IW) fusion method [30]. In addition, we also give the fusion results employing a competitive pansharpening method (GFTD) proposed by [58] with the best parameters authors suggested. The IW fusion method used the ‘Db4’ filter. The patch size in the sparse coding stage of the proposed method was 8 × 8, and the sliding window’s step size was set to 1. The DCT dictionary was chosen for generality in this work, and the error tolerance ε in SR-based low-frequency components fusion was set to 0.01. The source code of the proposed method can be downloaded at http://sendimage.whu.edu.cn/send-resource-download/ (accessed on 27 January 2022).

3.2. Experimental Results

Figure 4 shows the fusion results of different methods in area 1, which features agricultural lands. The integration of MS and SAR images can effectively enhance the transitions and boundaries between buildings, roads, and cultivated lands and retain the colors at the same time. Figure 4a,b shows the HR SAR image and up-sampled LR multispectral image, which are the input images for various fusion methods. The HR multispectral image is provided for spectral evaluation (but not the ground truth), as shown in Figure 4j. Compared with Figure 4b, Figure 4c–h demonstrates that all the fusion methods improve the sharpness of MS image and effectively enhance the spatial and texture features. As can be seen from Figure 4c,d, the spatial enhancement of IHS and BTH fusion is obvious, but the colors also appear unnatural. The spectral distortions of agricultural land in the upper right corner and the building in the lower left are serious. At the same time, serious noises can be seen in the red circle marked in the upper right corner in Figure 4c,d. Compared with the original HR multispectral image, the fusion results of AWLPR and IW fusion (see Figure 4f,g) perform well in spectral preservation but integrate limited spatial and texture features. Besides, obvious artifacts can be seen near the boundaries of ground features. Figure 4h (GFTD fusion method) exhibits appealing visual performance, while it suffers from aberrant colors in the areas marked with yellow arrows. Figure 4e,i shows that the LP-based and proposed methods strike the balance between spectral preservation and spatial enhancement. The colors of these two methods are comparable to the HR multispectral image, and the buildings and boundaries are markedly enhanced. Moreover, the noises and artifacts are effectively avoided. Compared with the LP-based method, the fusion result of the LPSR method looks sharper in the region marked with a green rectangle, and overall, the image is clearer, with stronger contrast.
The fact that the ground truth image of multispectral and SAR image fusion does not exist makes the quantitative evaluation challenging and less reliable. Therefore, we just give the qualitative comparison for supplements. In order to compare the spectral fidelity of different fusion results, the first property of Wald’s protocol [59] is checked via calculating the spectral indexes between the spatially degraded HR multispectral image and the degraded fusion images (reduced resolution). For testing the second and third property of Wald’s protocol, we compared the differences of original HR multispectral and fused images (full resolution) in the spectral aspect. Compared with the pansharpening field, estimation of the similarity between the fused and HR SAR image in spatial details is much more frequently used and important. In order to analyze the spatial quality, Zhou’s procedure [53] is adopted to measure the spatial similarity between SAR and the fused images. Additionally, a comprehensive metric based on information theory [55] is used to calculate the amount of gradient information conducted from both the input multispectral and SAR images to the fused images.
Table 1 displays the quantitative evaluation results of area 1. Overall, the statistical indicators are in agreement with the visual inspection. In terms of two spectral fidelity evaluation indicators, CC and SAM, the competitive pansharpening fusion methods GFTD and AWLPR have the strongest spectral preservation abilities, no matter if in the reduced or full spatial resolution. LP-based, IW, and the proposed LPSR fusion methods are a little bit weaker than the GFTD and AWLPR methods. IHS and BTH fusion methods have the worst results in all spectral fidelity indexes. Moreover, the CC of IHS method is almost zero, which indicates that the fusion result has no correlation with the multispectral image because of losing too many spectral characteristics. With regard to the spatial enhancement, IHS fusion method performs the best in SCC and SSIM indexes, which demonstrates that IHS fusion method integrates more spatial and texture features so that the fused images are closer to the HR SAR image. The SCC and SSIM of the proposed method are lower than the IHS method but better than the other methods. The BTH and GFTD method are the worst in terms of SCC and SSIM indicators because the intensity and histogram-matched SAR image have low contrast and limited details due to the low correlation between MS and SAR images. This phenomenon is extraordinarily common when simply borrowing the method from the pansharpening field. The LP-based, AWLPR, and IW methods lose edges because some structures and contents contained in the low frequency of SAR are discarded. The AGs of all the methods are higher than that of the original LR multispectral image (5.0008), which proves that the multispectral and SAR image fusion can enhance the spatial and texture features of the original MS image. FMI indicates that IHS and our method integrate more gradient information into the fused image from the source images.
Figure 5 exhibits the fusion results in area 2. Area 2 is dominated by buildings and roads. Consequently, the goal of image fusion is to improve the spatial resolution and effectively enhance the interpretability of the image. As can be seen from Figure 5, all the methods enhance the contrast of MS image and make the image more interpretable, which verifies the potentials of MS and SAR image fusion in urban areas. Figure 5c,d indicates that the IHS and BTH fusion methods can significantly improve the spatial resolution of the images. However, the fused images of IHS and BTH methods lose too many spectral characteristics of the multispectral image, which affects the visual interpretation. Figure 5e–h suggests that LP-based, AWLPR fusion, IW fusion, and GFTD methods can better maintain the spectral contents of MS image, whereas the fused images are blurred. In addition, the results of IW methods introduce visible artifacts in the boundaries. Again, GFTD method suffers from weird spectral distortion. In contrast, Figure 5i is the best result intuitively, which prevents the considerable disturbance of spectral contents and severe artifacts. Compared with Figure 5e, the proposed LPSR method shows better maintenance of spatial details and makes the image more interpretable than the LP-based method. Experimental results indicate that the proposed method achieves a satisfactory trade-off between spectral maintenance and spatial enhancement, which is the most suitable one in urban areas.
Table 2 displays the quantitative evaluation results of area 2, from which similar conclusions to area 1 can be drawn. AWLPR and GFTD methods perform the best in terms of indicators measuring spectral quality, followed by LP-based, IW, and the proposed LPSR methods. Again, the IHS method performed the worst in spectral preservation because it introduce excessive SAR features, including speckle noise. As to the indexes measuring spatial enhancement, IHS is the highest in SCC and SSIM. The proposed method is just inferior to the IHS method in SCC and SSIM and notably outperforms LP-based, AWLPR, IW, and GFTD methods. The SCC and SSIM indicators demonstrate that the proposed LPSR fusion method also has a satisfactory spatial enhancing ability. The AGs of all the methods are higher than that of the original LR multispectral image (13.6610), which proves the advantage of doing such fusion. Similarly, IHS has the highest FMI values, followed by our method, which are superior to the other fusion methods.
Figure 6 displays the fusion results in area 3. As shown in Figure 6, all the methods improve the spatial resolution of the LR multispectral image to some extent, indicating that these methods more or less integrate the spatial details of SAR image into the fused images. Figure 6c,d shows that the results of IHS and BTH fusion are similar, which is to say the results look sharper but suffer from severe color changes. Compared to the original HR multispectral image, the roads become black and the color of vegetation changes greatly. IHS and BTH fusion methods put too strong of an emphasis on spatial detail enhancement, resulting in the lack of necessary interpretability of the fusion results. Figure 6e–g shows that LP-based, AWLPR, and IW fusion methods maintain the spectral characteristics better. For example, the colors of grasslands and roads appear more natural and are closer to the colors in HR multispectral image. This is because LP-based, AWLPR, and IW fusion methods only incorporate the high-frequency information of SAR image and leave the low-frequency components of MS image unchanged, resulting in small spectral distortion. Compared with Figure 6e, obvious artifacts can be seen near the linear objects such as roads in Figure 6g. This phenomenon is common in WT-based fusion methods. By contrast, LP-based fusion methods abandon the decimation operator in the spatial frequency domain without losing the spatial connectivity of the edges and textures, and thus effectively preserve the spatial information of the original images [44]. AWLPR and GFTD method (Figure 6f,h) obtains the spectrally closest result compared to the HR multispectral image, while the fused image looks blurry and with low contrast. As shown in Figure 6i, the proposed LPSR method not only maintains a satisfactory spectral fidelity, but also integrates abundant spatial and texture features.
Table 3 shows the quantitative evaluation results of area 3. IHS method performs the worst in terms of the two spectral evaluation indexes CC and SAM. The other methods perform much better in spectral assessment compared to his. Among them, AWLPR and GFTD achieve the best results in spectral fidelity. Indicators SCC and SSIM show that the proposed method has a better spatial enhancement ability than BTH, LP-based, AWLPR, IW, and GFTD methods and is only worse than the IHS method. The AGs of all the methods are higher compared to that of the original LR multispectral image (7.7424), indicating that the fused images are clearer and contain more abundant spatial and texture information. Similar to the aforementioned two areas, the FMI values of IHS and the proposed method are much higher than other methods.
The evaluation metrics for CI are correlation coefficient on the full spatial resolution (CCfull), structural similarity (SSIM), and feature mutual information (FMI). We selected them because they measure the spectral preservation, spatial enhancement, and information incorporation abilities. Moreover, the value ranges of them are similar. The weights are equal for them because we considered the abilities of spectral preservation, spatial enhancement and information incorporation equally in this paper. Figure 7 illustrates the CI values of all the methods in three areas, from which we can get a more intuitive impression on the quantitative performance. The proposed LPSR fusion method obtains the highest CI values in all the three areas, which indicates that our method achieves the best performances from a comprehensive angle considering spectral, spatial, and information characteristics. In addition, the order of CI values highly accords with visual inspection results, which validates the effectiveness of the proposed index.
To verify the stability of the results disregarding the variability of input images, we conducted experiments on a larger set of images, which consists of 116 pairs of multispectral and SAR images of size 100 × 100. Table 4 exhibits the average performance and standard deviation of each method across all the images. CS-based methods (IHS and BTH) behave the worst in terms of spectral information preservation (indexes CC and SAM), while MRA-based method AWLPR achieves the best results in spectral fidelity indexes CC and SAM. Regarding the spatial enhancement performance, CS-based method IHS introduces the most spatial and texture features in terms of SCC and SSIM, while MRA-based methods (LP and AWLPR) behave poorly in maintaining spatial characteristics. All methods improve the average gradients of the LR multispectral images (9.4163), which demonstrates that the MS-SAR fusion can powerfully incorporates SAR features into MS images and improve the information richness of remotely sensed images. The extraordinarily high value of AG from BTH may owe to the outliers in heterogeneous areas (such as in Figure 4d). The proposed method strikes a satisfactory trade-off between spectral and spatial preservation, which can be validated by the CI value in Figure 8.
Based on the consideration of visual inspection and statistical analysis, the proposed LPSR fusion method can significantly enhance the spatial and texture characteristics of multispectral images, while preserving the spectral signature of ground features at the same time. Compared with the CS-based (IHS and BTH), MRA-based (LP-based and AWLPR), hybrid (IW), and GFTD methods, the LPSR method achieves the best balance between spectral characteristics preservation and spatial and texture features enhancement.

4. Discussion

4.1. Adjustment Capability

In practical usage, different application scenarios may have entirely different requirements on what the multispectral and SAR fusion images should be like. The flexible adjustability will make the fusion method more powerful and meet different needs in practical use. This section explores the adjustability of fusion results via tuning the decomposition level of the Laplacian pyramid to make the fusion results strike different balances between spectral preservation or spatial and texture enhancement. Therefore, the decomposition levels 1–5 are chosen for comparison experiment, and the corresponding fusion results are denoted as LPSR-1 to LPSR-5. In addition, the representative CS-based method (IHS fusion) and MRA-based method (LP-based fusion) are included as baselines of spatial enhancement and spectral preservation. Relevant experimental results are shown in Figure 9, Figure 10 and Figure 11.
As can be seen from Figure 9, Figure 10 and Figure 11, the fused image incorporates more characteristics from SAR image. Therefore, the spatial and texture enhancement is more significant when the decomposition level arises from 1 to 5. Nevertheless, the color changes more considerably when the decomposition level increases. When the decomposition level is 1 or 2, the fused image is visually similar to the LP-based fusion image in spectral consistency, while the spatial improvement is close to the IHS fusion method when the decomposition level reaches 4 or 5. Therefore, different trade-offs can be achieved between spectral fidelity and spatial improvement by adjusting the decomposition level of the Laplacian pyramid. From the statistical results in Figure 12, we can draw a similar conclusion. When the decomposition level equals to 1, the proposed LPSR fusion method achieves a good performance in spectral index CC, which is close to the LP-based method. While the decomposition level increases from 1 to 5, the spatial evaluation indexes (SCC and SSIM) of the proposed LPSR method approach the IHS fusion method, which indicates that the fusion result can be converted from the mode emphasizing spectral preservation to focusing on spatial and texture details injection. The verified adjustability makes the proposed fusion method more applicable and can satisfy different application requirements.

4.2. Time Complexity

As mentioned before, the combined use of LP and sparse representation can overcome the drawback of the stand-alone SR-based fusion method because the most time-consuming stage, which is to say the sparse coding operation, is only performed in small-sized low-pass bands. In this part, we are going to discuss the relationship between time complexity and Laplacian pyramid level. All the fusion methods in this part are implemented in MATLAB 2020b on a computer with an Intel Core i9-10900K 3.7 GHz CPU and 128 GB RAM.
As can be seen in Table 5, the stand-alone SR-based method consumes a high computing burden. Our proposed method with only one decomposition level can also sharply reduce the running time to some extent. As the decomposition level increases, the low-pass component has a smaller size, leading to a lighter sparse coding operation. Therefore, the time complexity is also reduced. The running time of proposed method with decomposition level 5 is comparable to those of the traditional methods. This part confirms that the computational efficiency can be improved when sparse representation is combined with LP.

5. Conclusions

This paper proposes a fusion method combining the Laplacian pyramid and sparse representation. Experimental results demonstrate that the proposed method strikes the balance between spectral preservation and spatial enhancement in terms of visual inspection and statistical comparison. The fusion framework integrates the merits of multi-scale transform and sparse representation theory, leading to the prevention of artifacts and inclination to smooth details. Additionally, the capability of achieving different balances between spectral consistency and detail improvement extends the applicability and practicability of the proposed method. Evaluation of pixel-level multispectral and SAR image fusion is challenging, considering that data fusion is more or less data-driven and application-oriented. Therefore, future work will concern more about mapping and classification tasks with and without fusion strategy, in order to validate the advantages of products derived from multispectral and SAR image fusion over datasets generated by stand-alone data sources.

Author Contributions

Conceptualization, H.S. and H.Z.; methodology, H.Z.; software, H.Z.; validation, H.Z.; formal analysis, H.Z.; investigation, H.Z.; resources, H.Z. and H.S.; data curation, H.Z.; writing—original draft preparation, H.Z.; writing—review and editing, H.Z., Q.Y., H.S., and X.G.; visualization, H.Z.; supervision, H.S. and Q.Y.; project administration, H.S.; funding acquisition, H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by National Natural Science Foundation of China (NSFC), grant number 41971303.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to thank Olaf Hellwich (Berlin University of Technology) for providing the datasets for experiments. We are also grateful to the editor and the anonymous reviewers for their detailed review, valuable comments, and constructive suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pohl, C.; van Genderen, J.L. Review article multisensor image fusion in remote sensing: Concepts, methods and applications. Int. J. Remote Sens. 1998, 19, 823–854. [Google Scholar] [CrossRef] [Green Version]
  2. Schmitt, M.; Zhu, X.X. Data fusion and remote sensing: An ever-growing relationship. IEEE Geosci. Remote Sens. Mag. 2016, 4, 6–23. [Google Scholar] [CrossRef]
  3. Pastorino, M.; Montaldo, A.; Fronda, L.; Hedhli, I.; Moser, G.; Serpico, S.B.; Zerubia, J. Multisensor and Multiresolution Remote Sensing Image Classification through a Causal Hierarchical Markov Framework and Decision Tree Ensembles. Remote Sens. 2021, 13, 849. [Google Scholar] [CrossRef]
  4. Tong, X.; Luo, X.; Liu, S.; Xie, H.; Chao, W.; Liu, S.; Liu, S.; Makhinov, A.; Makhinova, A.; Jiang, Y. An approach for flood monitoring by the combined use of Landsat 8 optical imagery and COSMO-SkyMed radar imagery. J. Photogramm. Remote Sens. 2018, 136, 144–153. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Zhang, H.; Lin, H. Improving the impervious surface estimation with combined use of optical and SAR remote sensing images. Remote Sens. Environ. 2014, 141, 155–167. [Google Scholar] [CrossRef]
  6. Bao, N.; Li, W.; Gu, X.; Liu, Y. Biomass Estimation for Semiarid Vegetation and Mine Rehabilitation Using Worldview-3 and Sentinel-1 SAR Imagery. Remote Sens. 2019, 11, 2855. [Google Scholar] [CrossRef] [Green Version]
  7. Samadhan, C.K.; Priti, P.R. Pixel Level Fusion Techniques for SAR and Optical Images: A Review. Inf. Fusion 2020, 59, 13–19. [Google Scholar]
  8. Nomura, R.; Oki, K. Downscaling of MODIS NDVI by Using a Convolutional Neural Network-Based Model with Higher Resolution SAR Data. Remote Sens. 2021, 13, 732. [Google Scholar] [CrossRef]
  9. Kong, Y.; Hong, F.; Leung, H.; Peng, X. A Fusion Method of Optical Image and SAR Image Based on Dense-UGAN and Gram–Schmidt Transformation. Remote Sens. 2021, 13, 4274. [Google Scholar] [CrossRef]
  10. Jiang, M.; Shen, H.; Li, J.; Zhang, L. An Integrated Framework for the Heterogeneous Spatio-Spectral-Temporal Fusion of Remote Sensing Images. arXiv 2021, arXiv:2109.00400. [Google Scholar]
  11. Lin, J.; Zhao, L.; Wang, Q.; Ward, R.; Wang, Z.J. DT-LET: Deep transfer learning by exploring where to transfer. Neurocomputing 2020, 390, 99–107. [Google Scholar] [CrossRef] [Green Version]
  12. Lin, J.; Mou, L.; Yu, T.; Zhu, X.; Wang, Z.J. Dual Adversarial Network for Unsupervised Ground/Satellite-To-Aerial Scene Adaptation. In Proceedings of the 28th ACM International Conference on Multimedia, Seattle, WA, USA, 12–16 October 2020. [Google Scholar]
  13. Lin, J.; Yu, T.; Mou, L.; Zhu, X.; Wang, Z.J. Unifying top–down views by task-specific domain adaptation. IEEE Trans. Geosci. Remote Sens. 2020, 59, 4689–4702. [Google Scholar] [CrossRef]
  14. Harris, J.R.; Murray, R.; Hirose, T. IHS transform for the integration of radar imagery with other remotely sensed data. Photogramm. Eng. Remote Sens. 1990, 56, 1631–1641. [Google Scholar]
  15. Yonghong, J.; Deren, L. Comparison of IHS transformation for integrating SAR and TM images. J. Remote. Sens. Land Res. 1997, 9, 34–39. [Google Scholar]
  16. Shettigara, V.K. A generalized component substitution technique for spatial enhancement of multispectral images using a higher resolution data set. Photogramm. Eng. Remote Sens. 1992, 58, 561–567. [Google Scholar]
  17. Yonghong, J. Fusion of landsat TM and SAR images based on principal component analysis. Remote Sens. Technol. Appl. 2012, 13, 46–49. [Google Scholar]
  18. Laben, C.; Brower, B. Process for Enhancing the Spatial Resolution of Multispectral Imagery Using Pan-Sharpening. U.S. Patent Application 6,011,875, 4 January 2000. [Google Scholar]
  19. Gillespie, A.R.; Kahle, A.B.; Walker, R.E. Color enhancement of highly correlated images. I. Decorrelation and HSI contrast stretches. Remote Sens. Environ. 1986, 20, 209–235. [Google Scholar] [CrossRef]
  20. Garzelli, A. Wavelet-based fusion of optical and SAR image data over urban area. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2002, 34, 59–62. [Google Scholar]
  21. Tateishi, R.; Wikantika, K.; Munadi, K.; Aslam, M. Study on the spectral quality preservation derived from multisensor image fusion techniques between JERS-1 SAR and Landsat TM data. In Proceedings of the 2003 IEEE International Geoscience and Remote Sensing Symposium, IGARSS ’03, Toulouse, France, 21–25 July 2003. [Google Scholar]
  22. Abdikan, S. Exploring image fusion of ALOS/PALSAR data and LANDSAT data to differentiate forest area. Geocarto Int. 2018, 33, 21–37. [Google Scholar] [CrossRef]
  23. Burt, P.J.; Adelson, E.H. Merging images through pattern decomposition. Proc. SPIE Appl. Digit. Image Process. 1985, 575, 173–181. [Google Scholar]
  24. 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]
  25. Cheng, J.; Liu, H.; Liu, T.; Wang, F.; Li, H. Remote sensing image fusion via wavelet transform and sparse representation. ISPRS J. Photogramm. Remote Sens. 2015, 104, 158–173. [Google Scholar] [CrossRef]
  26. Zhang, W.; Yu, L. SAR and Landsat ETM+ image fusion using variational model. In Proceedings of the 2010 International Conference on Computer and Communication Technologies in Agriculture Engineering, Chengdu, China, 12–13 June 2010. [Google Scholar]
  27. Huang, B.; Li, Y.; Han, X.; Cui, Y.; Li, W.; Li, R. Cloud removal from optical satellite imagery with SAR imagery using sparse representation. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1046–1050. [Google Scholar] [CrossRef]
  28. Alparone, L.; Baronti, S.; Garzelli, A.; Nencini, F. Landsat ETM+ and SAR image fusion based on generalized intensity modulation. IEEE Trans. Geosci. Remote Sens. 2004, 42, 2832–2839. [Google Scholar] [CrossRef]
  29. Chibani, Y. Additive integration of SAR features into multispectral SPOT images by means of the à trous wavelet decomposition. J. Photogramm. Remote Sens. 2006, 60, 306–314. [Google Scholar] [CrossRef]
  30. Hong, G.; Zhang, Y.; Mercer, B. A wavelet and IHS integration method to fuse high resolution SAR with moderate resolution multispectral images. Photogramm. Eng. Remote Sens. 2009, 75, 1213–1223. [Google Scholar] [CrossRef]
  31. Yin, Z. Fusion algorithm of optical images and SAR with SVT and sparse representation. Int. J. Smart Sens. Intell. Syst. 2015, 8, 1123–1141. [Google Scholar]
  32. Shao, Z.; Wu, W.; Guo, S. IHS-GTF: A Fusion Method for Optical and Synthetic Aperture Radar Data. Remote Sens. 2020, 12, 2796. [Google Scholar] [CrossRef]
  33. Wright, J.; Yang, A.Y.; Ganesh, A.; Sastry, S.S.; Ma, Y. Robust face recognition via sparse representation. IEEE Trans. Pattern Anal. Mach. Intell. 2008, 31, 210–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Elad, M.; Aharon, M. Image denoising via sparse and redundant representations over learned dictionaries. IEEE Trans. Image Process. 2006, 15, 3736–3745. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, J.; Wright, J.; Huang, T.S.; Ma, Y. Image super-resolution via sparse representation. IEEE Trans. Image Process. 2010, 19, 2861–2873. [Google Scholar] [CrossRef] [PubMed]
  36. Olshausen, B.A.; Field, D.J. Emergence of simple-cell receptive field properties by learning a sparse code for natural images. Nature 1996, 381, 607–609. [Google Scholar] [CrossRef] [PubMed]
  37. Bruckstein, A.M.; Donoho, D.L.; Elad, M. From sparse solutions of systems of equations to sparse modeling of signals and images. SIAM Rev. 2009, 51, 34–81. [Google Scholar] [CrossRef] [Green Version]
  38. Yang, B.; Li, S. Multifocus image fusion and restoration with sparse representation. IEEE Trans. Instrum. Meas. 2009, 59, 884–892. [Google Scholar] [CrossRef]
  39. Yu, N.; Qiu, T.; Bi, F.; Wang, A. Image features extraction and fusion based on joint sparse representation. IEEE J. Sel. Top. Signal Process. 2011, 5, 1074–1082. [Google Scholar] [CrossRef]
  40. 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]
  41. Zhu, X.X.; Bamler, R. A sparse image fusion algorithm with application to pan-sharpening. IEEE Trans. Geosci. Remote Sens. 2012, 51, 2827–2836. [Google Scholar] [CrossRef]
  42. 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]
  43. 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] [Green Version]
  44. Aiazzi, B.; Alparone, L.; Baronti, S.; Garzelli, A. Context-driven fusion of high spatial and spectral resolution images based on oversampled multiresolution analysis. IEEE Trans. Geosci. Remote Sens. 2002, 40, 2300–2312. [Google Scholar] [CrossRef]
  45. Wang, W.; Chang, F. A Multi-focus Image Fusion Method Based on Laplacian Pyramid. J. Comput. 2011, 6, 2559–2566. [Google Scholar] [CrossRef]
  46. Burt, P.; Adelson, E. The Laplacian pyramid as a compact image code. IEEE Trans. Commun. 1983, 31, 532–540. [Google Scholar] [CrossRef]
  47. Yang, B.; Li, S. Pixel-level image fusion with simultaneous orthogonal matching pursuit. Inf. Fusion 2012, 13, 10–19. [Google Scholar] [CrossRef]
  48. Tropp, J.A.; Gilbert, A.C.; Strauss, M.J. Algorithms for simultaneous sparse approximation. Part I: Greedy pursuit. Signal Process. 2006, 86, 572–588. [Google Scholar] [CrossRef]
  49. Hellwich, O.; Reigber, A.; Lehmann, H. Sensor and data fusion contest: Test imagery to compare and combine airborne SAR and optical sensors for mapping. In Proceedings of the IEEE International Geosciences and Remote Sensing Symposium and the 24th Canadian Symposium on Remote Sensing, Toronto, ON, Canada, 24–28 June 2002. [Google Scholar]
  50. Zhang, J.; Yang, J.; Zhao, Z.; Li, H.; Zhang, Y. Block-regression based fusion of optical and SAR imagery for feature enhancement. Int. J. Remote Sens. 2010, 31, 2325–2345. [Google Scholar] [CrossRef]
  51. Ma, X.; Shen, H.; Zhao, X.; Zhang, L. SAR image despeckling by the use of variational methods with adaptive nonlocal functionals. IEEE Trans. Geosci. Remote Sens. 2016, 54, 3421–3435. [Google Scholar] [CrossRef]
  52. Yuhas, R.H.; Goetz, A.F.H.; Boardman, J.W. Discrimination Among Semi-Arid Landscape Endmembers Using the Spectral Angle Mapper (SAM) Algorithm. In JPL, Summaries of the Third Annual JPL Airborne Geoscience Workshop; JPL Publication: Washington, DC, USA; NASA: Washington, DC, USA, 1992; pp. 147–149. [Google Scholar]
  53. Zhou, J.; Civco, D.L.; Silander, J. A wavelet transform method to merge Landsat TM and SPOT panchromatic data. Int. J. Remote Sens. 1998, 19, 743–757. [Google Scholar] [CrossRef]
  54. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [Green Version]
  55. Haghighat, M.B.A.; Aghagolzadeh, A.; Seyedarabi, H. A non-reference image fusion metric based on mutual information of image features. Comput. Electr. Eng. 2011, 37, 744–756. [Google Scholar] [CrossRef]
  56. Lolli, S.; Alparone, L.; Garzelli, A.; Vivone, G. Haze correction for contrast-based multispectral pansharpening. IEEE Geosci. Remote Sens. Lett. 2017, 14, 2255–2259. [Google Scholar] [CrossRef]
  57. Alparone, L.; Garzelli, A.; Vivone, G. Intersensor statistical matching for pansharpening: Theoretical issues and practical solutions. IEEE Trans. Geosci. Remote Sens. 2017, 55, 4682–4695. [Google Scholar] [CrossRef]
  58. Meng, X.; Li, J.; Shen, H.; Zhang, L.; Zhang, H. Pansharpening with a guided filter based on three-layer decomposition. Sensors 2016, 16, 1068. [Google Scholar] [CrossRef] [PubMed]
  59. Wald, L.; Ranchin, T.; Mangolini, M. Fusion of satellite images of different spatial resolutions: Assessing the quality of resulting images. Photogramm. Eng. Remote Sens. 1997, 63, 691–699. [Google Scholar]
Figure 1. Schematic diagram of proposed fusion method.
Figure 1. Schematic diagram of proposed fusion method.
Remotesensing 14 00870 g001
Figure 2. The flowchart of SR-based low-frequency component fusion.
Figure 2. The flowchart of SR-based low-frequency component fusion.
Remotesensing 14 00870 g002
Figure 3. Multispectral and SAR images used for fusion experiments, located in the Trudering region of Munich in Germany. (a) Airborne multispectral images; (b) airborne SAR intensity images.
Figure 3. Multispectral and SAR images used for fusion experiments, located in the Trudering region of Munich in Germany. (a) Airborne multispectral images; (b) airborne SAR intensity images.
Remotesensing 14 00870 g003aRemotesensing 14 00870 g003b
Figure 4. Fusion results of different methods in area 1. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Figure 4. Fusion results of different methods in area 1. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Remotesensing 14 00870 g004
Figure 5. Fusion results of different methods in area 2. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Figure 5. Fusion results of different methods in area 2. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Remotesensing 14 00870 g005
Figure 6. Fusion results of different methods in area 3. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Figure 6. Fusion results of different methods in area 3. (a) HR SAR image; (b) up-sampled LR multispectral image; (c) IHS fusion; (d) BTH fusion; (e) LP-based fusion; (f) AWLPR fusion; (g) IW fusion; (h) GFTD fusion; (i) proposed LPSR fusion; (j) original HR multispectral image.
Remotesensing 14 00870 g006
Figure 7. Bar chart of the comprehensive index.
Figure 7. Bar chart of the comprehensive index.
Remotesensing 14 00870 g007
Figure 8. The box plot and mean values of the comprehensive index on a large set of images.
Figure 8. The box plot and mean values of the comprehensive index on a large set of images.
Remotesensing 14 00870 g008
Figure 9. Fusion results of adjustability experiment in area 1. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Figure 9. Fusion results of adjustability experiment in area 1. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Remotesensing 14 00870 g009
Figure 10. Fusion results of adjustability experiment in area 2. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Figure 10. Fusion results of adjustability experiment in area 2. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Remotesensing 14 00870 g010
Figure 11. Fusion results of adjustability experiment in area 3. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Figure 11. Fusion results of adjustability experiment in area 3. (a) IHS fusion; (b) LPSR-5; (c) LPSR-4; (d) LPSR-3; (e) LPSR-2; (f) LPSR-1; (g) LP-based fusion; (h) original HR multispectral image.
Remotesensing 14 00870 g011
Figure 12. Statistical results of adjustability experiment: (ac) display the CC f u l l , SCC, SSIM values of area 1; (df) for area 2; (gi) for area 3.
Figure 12. Statistical results of adjustability experiment: (ac) display the CC f u l l , SCC, SSIM values of area 1; (df) for area 2; (gi) for area 3.
Remotesensing 14 00870 g012
Table 1. Quantitative evaluation results of different fusion methods in area 1.
Table 1. Quantitative evaluation results of different fusion methods in area 1.
IndexesIHSBTHLPAWLPRIWGFTDLPSRIdeal
CC 10.04180.42100.96200.97530.96080.97510.92351
0.02470.35660.86630.89170.85980.89640.82381
SAM 24.76682.98801.36851.00671.20481.02491.97900
7.53315.60515.11604.36804.77994.38355.56350
SCC0.97460.64360.95040.95060.95310.87140.95601
SSIM0.89980.68620.70230.66320.66560.65720.73401
AG20.45219.740321.045820.693621.091016.139220.8458+∞
FMI0.51610.40610.47940.46780.48030.46520.50961
1,2 The first row of indicator CC/SAM shows the spectral fidelity in the reduced spatial resolution, in order to check the first property of Wald’s protocol. The second row corresponds to the results in the full spatial resolution, which verifies the second and third property of Wald’s protocol.
Table 2. Quantitative evaluation results of different fusion methods in area 2.
Table 2. Quantitative evaluation results of different fusion methods in area 2.
IndexesIHSBTHLPAWLPRIWGFTDLPSRIdeal
CC 1−0.19510.03040.86790.92880.85820.90440.74011
−0.1832−0.01420.68460.80290.65980.72670.56031
SAM 26.46273.32833.89412.63143.71142.77484.63090
14.26319.631110.40759.052510.26109.475210.89050
SCC0.98420.57340.84980.54530.88210.85710.91501
SSIM0.96810.74280.40460.14290.34660.32420.57771
AG24.641533.344125.755219.411824.735126.254423.4259+∞
FMI0.51300.37850.40330.39740.40710.42170.45251
1,2 The first row of indicator CC/SAM shows the spectral fidelity in the reduced spatial resolution, in order to check the first property of Wald’s protocol. The second row corresponds to the results in the full spatial resolution, which verifies the second and third property of Wald’s protocol.
Table 3. Quantitative evaluation results of different fusion methods in area 3.
Table 3. Quantitative evaluation results of different fusion methods in area 3.
IndexesIHSBTHLPAWLPRIWGFTDLPSRIdeal
CC 1−0.19420.28460.87960.93350.84110.92120.76901
−0.19230.20080.67120.79260.61510.74810.56681
SAM 24.12202.19202.36371.57172.03411.65253.04980
7.59905.66086.34454.97565.88445.19186.99010
SCC0.99720.91990.91380.61630.94580.91390.96291
SSIM0.95830.79200.48850.21150.43760.34680.65341
AG16.914917.225517.839411.624717.177813.9111116.9217+∞
FMI0.51640.45360.43470.42230.43260.43690.48731
1,2 The first row of indicator CC/SAM shows the spectral fidelity in the reduced spatial resolution, in order to check the first property of Wald’s protocol. The second row corresponds to the results in the full spatial resolution, which verifies the second and third property of Wald’s protocol.
Table 4. Quantitative analysis on a large set of images.
Table 4. Quantitative analysis on a large set of images.
MethodsCC 1SAM 2SCCSSIMAGFMI
IHS−0.1263 ± 0.23007.7138 ± 5.94780.9551 ± 0.05630.8416 ± 0.115018.5431 ± 8.50310.5029 ± 0.0314
BTH0.1577 ± 0.25316.8539 ± 4.29260.6217 ± 0.24890.5792 ± 0.190228.9911 ± 12.51800.3723 ± 0.0569
LP0.6617 ± 0.11316.1318 ± 4.18530.8661 ± 0.05970.5488 ± 0.131920.1218 ± 8.74190.4266 ± 0.0342
AWLPR0.6933 ± 0.11205.4492 ± 3.50530.8731 ± 0.05810.5031 ± 0.130821.7485 ± 9.71970.4326 ± 0.0228
IW0.6235 ± 0.13385.9936 ± 4.18490.8851 ± 0.05280.5375 ± 0.120519.3256 ± 8.55870.4315 ± 0.0355
GFTD0.6707 ± 0.13565.6733 ± 3.80640.7607 ± 0.15820.4773 ± 0.136321.2761 ± 12.00010.4188 ± 0.0317
LPSR0.5696 ± 0.13206.6137 ± 4.53350.9070 ± 0.04930.6399 ± 0.113418.9086 ± 8.14010.4699 ± 0.0315
Ideal1011+∞1
1,2 The values of indicator CC/SAM correspond to the results in the full spatial resolution.
Table 5. Running time (in seconds) of the proposed method with different LP levels.
Table 5. Running time (in seconds) of the proposed method with different LP levels.
SRLPSR-1LPSR-2LPSR-3LPSR-4LPSR-5GFTDIWAWLPRLPBTHIHS
Time/s9308.662226.38536.63119.5923.362.950.28960.11370.20470.08300.13620.0234
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, H.; Shen, H.; Yuan, Q.; Guan, X. Multispectral and SAR Image Fusion Based on Laplacian Pyramid and Sparse Representation. Remote Sens. 2022, 14, 870. https://doi.org/10.3390/rs14040870

AMA Style

Zhang H, Shen H, Yuan Q, Guan X. Multispectral and SAR Image Fusion Based on Laplacian Pyramid and Sparse Representation. Remote Sensing. 2022; 14(4):870. https://doi.org/10.3390/rs14040870

Chicago/Turabian Style

Zhang, Hai, Huanfeng Shen, Qiangqiang Yuan, and Xiaobin Guan. 2022. "Multispectral and SAR Image Fusion Based on Laplacian Pyramid and Sparse Representation" Remote Sensing 14, no. 4: 870. https://doi.org/10.3390/rs14040870

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