国土资源遥感, 2018, 30(2): 100-106 doi: 10.6046/gtzyyg.2018.02.14

技术方法

一种国产高分卫星遥感影像变分融合方法

尹峰,1, 孟祥超,2, 梁鹏1

1.湖北省国土资源研究院,武汉 430071

2. 武汉大学资源与环境科学学院,武汉 430079

A variational fusion method for remote sensing images of China’s domestic high-resolution satellites

YIN Feng,1, MENG Xiangchao,2, LIANG Peng1

1. Hubei Institute of Land and Resources, Wuhan 430071, China

2.School of Resources and Environmental Sciences, Wuhan University, Wuhan 430079, China

通讯作者: 孟祥超(1989-),男,博士研究生,主要从事影像融合与遥感图像处理研究。Email:mengxc@whu.edu.cn

第一联系人:

第一作者: 尹 峰(1982-),男,硕士,主要从事遥感应用、国土资源调查等工作。Email: 89642740@qq.com

收稿日期: 2016-10-13   修回日期: 2016-12-20   网络出版日期: 2018-06-15

基金资助: 湖北省国土资源科研专项项目“基于多源高分遥感数据的多层次土地监测监管关键技术与应用研究”.  编号: ETZ2016A16

Received: 2016-10-13   Revised: 2016-12-20   Online: 2018-06-15

Fund supported: .  编号: ETZ2016A16

摘要

由于现有全色/多光谱融合方法对国产高分卫星遥感影像数据特点考虑不足,提出一种针对国产高分卫星遥感影像的变分融合方法。该方法充分考虑融合影像波段间光谱关系的保持,构建基于光谱梯度的三维光谱高保真模型,并针对国产高分卫星全色影像存在模糊降质的数据特点,发展顾及模糊降质的空间增强模型。在此基础上,结合影像先验知识建立融合目标函数,最后采用梯度下降法优化求解得到融合影像。通过高分一号(GF-1)和高分二号(GF-2)影像数据对提出的融合方法进行实验验证,并与典型GS,PRACS和ATWT-M3等融合方法分别从定性和定量2方面进行比较分析。实验结果表明,该融合方法充分考虑了国产高分卫星影像数据特点,在对照实验的几种方法中得到了最优的融合结果,可在有效提升多光谱影像空间分辨率的同时,很大程度上保持原有的光谱信息。

关键词: 融合 ; 多光谱影像 ; 全色影像 ; 变分 ; 遥感

Abstract

The current panchromatic (PAN) / multispectral (MS) fusion methods do not comprehensively consider the characteristics of the remote sensing images from China’s domestic high-resolution satellites. Therefore, this paper proposes a variational fusion method for China’s domestic high-resolution images. On the one hand, the three-dimensional spectral high-fidelity model based on the spectral gradient is proposed by comprehensive consideration of the relations between the spectral bands. On the other hand, according to the existing blurring characteristics of the PAN image acquired by China’s domestic high-resolution satellites, the spatial enhancement model in consideration of the blurring degradation is developed. Finally, the fusion energy function is constructed by combining the prior knowledge of the remote sensing images, and it is solved by the classical gradient decent methods to obtain the fused image. The proposed method was tested and verified by the Gaofen-1 (GF-1) and Gaofen-2 (GF-2) satellite datasets. In addition, the popular GS, PRACS, and ATWT-M3 methods were applied for comparison from both qualitative and quantitative aspects. The experimental results show that the proposed variational high-fidelity PAN/MS fusion method comprehensively considers the characteristics of China’s domestic satellites, and hence it can maximally preserve the spectral information while effectively improve the spatial resolution of the MS images, thus achieving the best fused results.

Keywords: fusion ; multispectral image ; panchromatic image ; variational model ; remote sensing

PDF (4022KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

尹峰, 孟祥超, 梁鹏. 一种国产高分卫星遥感影像变分融合方法. 国土资源遥感[J], 2018, 30(2): 100-106 doi:10.6046/gtzyyg.2018.02.14

YIN Feng, MENG Xiangchao, LIANG Peng. A variational fusion method for remote sensing images of China’s domestic high-resolution satellites. REMOTE SENSING FOR LAND & RESOURCES[J], 2018, 30(2): 100-106 doi:10.6046/gtzyyg.2018.02.14

0 引言

近年来,国产高分卫星得到了快速发展,然而,相比于QuickBird和IKONOS等国外卫星影像,部分国产高分卫星(如高分一号(GF-1)、高分二号(GF-2)和吉林一号等)各波段影像相对较模糊。因此,有针对性地发展国产高分卫星遥感影像全色/多光谱波段的融合方法,以提高影像数据的整体质量十分必要。

全色/多光谱影像融合方法起源于20世纪80年代[1,2]。1986年SPOT-1卫星同时提供全色和多光谱影像以来,融合方法得到了近30 a的快速发展。一般而言,融合方法可归为3类[3,4]: 成分替换类融合方法、多分辨率分析融合方法和基于模型的融合方法。其中,成分替换类方法是最简单的也是最流行的融合方法,该类方法已被广泛应用到了ENVI和ERDAS等专业遥感软件中。该类方法首先基于光谱变换得到亮度分量,然后通过全色影像替换亮度分量的方式对多光谱影像进行空间信息增强,经典方法有主成分分析融合方法[5]、Gram-Schimidt(GS)融合方法[6]、Intensity-Hue-Saturation (IHS) 融合方法[5]等。多分辨率分析融合方法基于小波变换或拉普拉斯金字塔等工具提取全色影像的高空间结构信息,并采用一定的注入模型将提取的空间结构信息注入到多光谱影像中得到高空间分辨率融合影像[7],如多孔小波融合方法[8]、拉普拉斯金字塔融合方法[9]和Contourlet 小波融合方法[10]等。针对成分替换类融合方法和多分辨率分析融合方法,Tu 等[11]进一步将其扩展到同一个融合框架,很大程度上促进了全色/多光谱融合方法的发展。

尽管已提出了大量成分替换类融合方法和多分辨率分析融合方法,然而这些方法都是根据某种假设进行简单正向求解得到融合影像,缺少强有力的数学理论基础和严密的逻辑关系,因此基于模型的影像融合方法 [12,13,14,15] 的开发得到了广泛关注。该类方法将融合影像的求解过程看成病态逆问题,基于影像观测模型建立能量函数,通过优化求解得到融合影像。其中,基于变分的融合方法[16,17]最具代表性,其建立的能量函数总体可分为3项: 光谱保真项、空间增强项和先验项,三者之中光谱保真项和空间增强项最为关键。然而现有方法中这2项对影像数据特点考虑不足,主要表现为光谱保真项仅简单考虑融合影像各波段与多光谱观测影像对应波段之间一对一的空间降质关系,对影像波段间光谱关系考虑不足; 空间增强项则未顾及国产卫星全色影像存在的模糊降质问题。

因此,针对上述问题,本文提出一种基于变分的国产高分卫星全色/多光谱融合方法。该方法充分考虑国产高分卫星影像特点,基于光谱梯度的三维光谱高保真项和顾及全色影像模糊降质的改进空间增强项,对国产高分卫星遥感影像在有效提升多光谱影像空间分辨率的同时,最大程度地保持其光谱信息。

1 理论与方法

充分考虑融合影像与全色、多光谱影像之间的关系,以及国产高分卫星影像特点,提出一种基于变分的全色/多光谱高保真融合方法。假设融合影像为X=[X1,X2,...,XB]T,其中B为波段数,原始多光谱影像为Y=[Y1,Y2,...,YB]T,全色影像为Z,则融合模型表示为

E(X)=fspectral(Y,X)+fspatial(Z,X)+fprior(X), (1)

式中: fspectral(Y,X)为光谱保真项,建立融合影像X与多光谱影像Y之间的关系; fspatial(Z,X)为空间增强项,建立融合影像X与全色影像Z之间的关系; fprior(X)为先验项。

1.1 改进的三维光谱高保真项

传统光谱保真项基于多光谱影像模型[1, 15, 18],建立融合影像各波段与多光谱影像相应波段之间一对一的空间降质关系,表示为

fspectral(Yb,Xb)=‖Yb-AXb2, (2)

式中:A表示模糊降采样过程; b表示波段序号。然而,该保真项对多光谱影像波段与波段之间的光谱关系考虑不足,故提出基于光谱梯度的三维光谱高保真项,即

fspectral(Yb,Xb)=‖Y'b-AX'b2, (3)

式中: Y'b= [Yb,Yb-Y1,...,Yb-YB]T; X'b= [Xb,Xb-X1,...,Xb-XB]T。该三维光谱高保真项不仅建立了融合影像与多光谱影像间的空间降质关系,同时基于光谱维梯度顾及了融合影像各波段间的光谱关系保持,可进一步提升融合模型的光谱保真能力。

1.2 改进的空间增强项

传统空间增强项假设理想高空间分辨率融合影像和全色影像具有相似的空间结构信息,然而,通过比较国产卫星(如GF-1,GF-2和吉林一号等)全色影像与国外卫星(QuickBird和IKONOS等)全色影像发现,国产高分卫星全色影像相对较模糊。基于梯度结构信息,提出顾及全色影像模糊降质的改进空间增强项,合理建立理想高空间分辨率融合影像与全色影像之间的关系,即

fspatial(Z,Xb)=‖▽Z-▽SXb2, (4)

式中: ▽=[▽H,▽V]T表示梯度算子,H表示水平方向,V表示垂直方向; S表示模糊算子。此外,考虑到多光谱影像各波段和全色影像灰度范围可能存在的不一致,在空间增强项中引入矩匹配操作,即

fspatial(Z,Xb)=‖▽Z-f(▽SXb) 2, (5)

式中f(·)表示矩匹配函数[19]

1.3 融合模型与优化求解

基于改进的三维光谱高保真项和顾及全色影像模糊降质的空间增强项,并结合经典拉普拉斯先验建立融合目标函数,即

E(Xb)= λ12Y'b-AX'b2+ 12||▽Z-f(▽SXb)+ λ22ΔXb2, (6)

式中λ1λ2为正则化参数。

针对目标函数,通过梯度下降法求解得到融合影像。

对式(6)求导,即

式中: ΔHΔV分别为水平方向和垂直方向上的拉普拉斯运算; (·)std为标准差运算。

通过连续的迭代逼近运算得到融合影像,迭代公式为

Xb,n+1= Xb,n-tb,n▽E(Xb,n), (8)

式中: n为迭代次数; tb,n为第b波段的迭代步长,通过对目标函数进行二阶泰勒级数展开得到。迭代终止条件为

Xn+1-Xn2Xn2d, (9)

式中d为预设迭代终止阈值,本文设置为10-7

2 实验与分析

选用GF-1和GF-2国产卫星影像进行融合实验。其中,GF-1影像数据获取地点为捷克某地,获取时间为2013年4月28日,全色影像空间分辨率为2 m,多光谱影像空间分辨率为8 m。GF-2影像数据获取地点为伊朗某地,获取时间为2014年9月4日,全色影像空间分辨率为0.81 m,多光谱影像空间分辨率为3.24 m。为了进行全面验证,所选用实验数据包含了植被、水体、建筑物和农田等多种地表覆盖类型,并基于模拟实验和真实实验,从定性和定量2方面分别对融合方法进行评价分析,其中定量评价采用全色/多光谱融合中4个最常用的评价指标,包括相关系数 (correlation coefficient,CC)[15]、峰值信噪比(peak signal-to-noise ratio,PSNR)[20]、相对全局误差 (relative dimensionless global error in synthesis,ERGAS)[21]和光谱角(spectral angle mapper,SAM)[15]。此外,提出方法与典型的GS融合方法、局部自适应成份替换融合方法(partial replacement adaptive component substitution,PRACS) [22]和ATWT-M3融合方法[8]进行比较分析。为了进一步验证改进的三维光谱高保真项在融合影像光谱保持上的优势,实验结果与基于传统光谱保真项的融合结果也进行了对比。本文提出方法中模型参数根据大量实验测试以人工经验设定,除特别说明外,参数设置为: λ1=10,λ2=0.001,空间增强项模糊核大小根据全色与多光谱空间分辨率比率设定为(2r+1)×(2r+1),其中空间分辨率比率r=4,方差为0.5。

2.1 模拟实验

根据Wald等[23]提出的获取参考影像的方法,首先,按全色和多光谱影像空间分辨率比率对原始影像数据进行空间降质; 然后,原始多光谱影像作为参考影像对融合结果进行评价。模拟实验采用GF-1和GF-2卫星影像数据,实验结果分别如图1图2所示。从定性评价方面,GF-1实验结果中GS方法具有较好的空间结构信息,但存在较大的光谱畸变; PRACS和ATWT-M3的融合结果光谱保持较好,但空间结构信息较模糊。相比于其他方法,本文方法在2种光谱保真项模式下均能取得较好的融合效果,在色彩上更加接近参考影像,同时空间结构更加清晰。与GF-1融合结果相比,GF-2模拟实验结果视觉上GS方法融合效果在光谱上有所改善; 其他方法展现出与GF-1类似的实验结果。进一步开展定量评价研究,定量评价结果如表1表2所示。

图1

图1   GF-1模拟实验结果

Fig.1   Fusion results of the GF-1 simulated experiment


图2

图2   GF-2模拟实验结果

Fig.2   Fusion results of the GF-2 simulated experiment


表1   GF-1模拟实验定量评价

Tab.1  Quantitative evaluation of the fusion results in the GF-1 simulated experiment

评价指标融合方法
GS方法PRACS方法ATWT-M3方法基于传统光谱
保真项融合方法
本文方法
CC0.841 70.857 90.860 70.884 80.894 0
PSNR24.246 525.101 325.112 125.394 725.610 3
ERGAS6.340 55.838 75.891 95.580 75.462 1
SAM4.885 54.316 04.574 74.529 54.204 1

新窗口打开| 下载CSV


表2   GF-2模拟实验定量评价

Tab.2  Quantitative evaluation of the fusion results in the GF-2 simulated experiment

评价指标融合方法
GS方法PRACS方法ATWT-M3方法基于传统光谱
保真项融合方法
本文方法
CC0.979 40.976 60.973 30.990 10.990 8
PSNR39.291 939.194 338.721 343.675 043.800 5
ERGAS0.672 10.707 40.779 70.423 60.416 2
SAM0.547 50.530 20.645 70.459 10.437 9

新窗口打开| 下载CSV


表1表2可以发现,相比于其他方法,本文方法定量评价结果最好,其中最优定量评价结果以加粗表示,次之以斜体表示。相比于传统光谱保真项,本文方法基于三维光谱保真项在光谱保持上也更有优势。其中,在CC,PSNR,ERGAS和SAM这4个定量评价指标中,SAM优势较为明显,这是因为改进的三维光谱保真项可更好地保持波段间光谱关系,具有更佳的光谱保持能力。

2.2 真实实验

基于原始GF-1和GF-2影像数据进行真实实验。由于在该实验中没有参考影像,根据Wald准则[23],将融合影像重采样到与原始多光谱影像相同的空间分辨率进行定量评价。实验结果分别如图3图4所示,表3表4分别为定量评价结果。通过与模拟实验结果对比发现,真实实验在定性和定量上总体展示了类似的结果。但是,GF-2真实实验中,GS融合结果在目视上具有较好的空间结构信息,这是因为在该实验数据中,多光谱影像白色建筑物在重采样过程中存在较为明显的膨胀现象,导致多光谱影像与全色影像存在一定的地物不匹配现象,GS方法通过全部成分替换的方式更有利于空间信息的增强,而这也造成了其融合影像的光谱信息存在一定的损失,可从表3表4定量评价指标上明显看出。在进行实验对比的各种方法中,本文提出的改进高保真融合方法在光谱保持和空间信息增强方面总体表现最好,相比于传统光谱保真项,在所有定量评价指标上均有所提升,其中以SAM提升最为明显,分别提升了19%和40%,具有明显的优势。

图3

图3   GF-1真实实验结果

Fig.3   Fusion results of the GF-1 real experiment


图4

图4   GF-2真实实验结果

Fig.4   Fusion results of the GF-2 real experiment


表3   GF-1真实实验定量评价

Tab.3  Quantitative evaluation of the fusion results in the GF-1 real experiment

评价指标融合方法
GS方法PRACS方法ATWT-M3方法基于传统光谱
保真项融合方法
本文方法
CC0.798 00.996 40.994 20.997 40.998 0
PSNR25.913 044.590 542.419 045.775 146.892 2
ERGAS4.239 40.542 20.697 50.449 60.403 2
SAM6.218 90.426 40.581 90.406 00.328 7

新窗口打开| 下载CSV


表4   GF-2真实实验定量评价

Tab.4  Quantitative evaluation of the fusion results in the GF-2 real experiment

评价指标融合方法
GS方法PRACS方法ATWT-M3方法基于传统光谱
保真项融合方法
本文方法
CC0.948 40.988 10.991 40.992 00.994 0
PSNR30.268 138.633 939.046 639.441 140.442 4
ERGAS1.153 50.491 70.492 10.445 90.391 6
SAM0.349 90.290 20.355 00.418 20.250 8

新窗口打开| 下载CSV


3 结论

本文针对现有全色/多光谱融合方法对国产高分卫星遥感影像数据特点考虑不足的缺点,提出了一种针对国产高分卫星全色/多光谱遥感影像的高保真变分融合方法。

1)在模型构建中,充分考虑融合影像波段间关系保持,发展了基于光谱梯度的三维光谱高保真模型,并考虑国产卫星全色影像模糊降质的问题,进一步发展了顾及全色降质的空间增强模型。

2)实验结果表明,本文方法相比于对照实验的其他方法,具有良好的光谱信息保持和空间信息增强能力,针对国产卫星影像可得到最优的融合结果。

3)不足之处在于现有方法求解效率相对较低,后期研究将引入并行计算等加速策略,进一步提升模型求解效率。

参考文献

Shen H F, Meng X C, Zhang L P .

An integrated framework for the spatio-temporal-spectral fusion of remote sensing images

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016,54(12):7135-7148.

DOI:10.1109/TGRS.2016.2596290      URL     [本文引用: 2]

Remote sensing satellite sensors feature a tradeoff between the spatial, temporal, and spectral resolutions. In this paper, we propose an integrated framework for the spatio-temporal-spectral fusion of remote sensing images. There are two main advantages of the proposed integrated fusion framework: it can accomplish different kinds of fusion tasks, such as multiview spatial fusion, spatio-spectral fusion, and spatio-temporal fusion, based on a single unified model, and it can achieve the integrated fusion of multisource observations to obtain high spatio-temporal-spectral resolution images, without limitations on the number of remote sensing sensors. The proposed integrated fusion framework was comprehensively tested and verified in a variety of image fusion experiments. In the experiments, a number of different remote sensing satellites were utilized, including IKONOS, the Enhanced Thematic Mapper Plus (ETM+), the Moderate Resolution Imaging Spectroradiometer (MODIS), the Hyperspectral Digital Imagery Collection Experiment (HYDICE), and Syste me Pour l' Observation de la Terre-5 (SPOT-5). The experimental results confirm the effectiveness of the proposed method.

张良培, 沈焕锋 .

遥感数据融合的进展与前瞻

[J]. 遥感学报, 2016,20(5):1050-1061.

[本文引用: 1]

Zhang L P, Shen H F .

Progress and future of remote sensing data fusion

[J]. Journal of Remote Sensing, 2016,20(5):1050-1061.

[本文引用: 1]

Aiazzi B, Alparone L, Baronti S , et al. Twenty-five years of pansharpening:A critical review and new developments[M] //Chen C H.Signal and Image Processing for Remote Sensing. 2nd edition.Boca Raton,FL: CRC Press, 2012: 533-548.

[本文引用: 1]

Vivone G, Alparone L, Chanussot J , et al.

A critical comparison among pansharpening algorithms

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014,53(5):2565-2586.

DOI:10.1109/IGARSS.2014.6946389      URL     [本文引用: 1]

In this paper state-of-the-art and advanced methods for multispectral pansharpening are reviewed and evaluated on two very high resolution datasets acquired by IKONOS-2 (four bands) and WorldView-2 (eight bands). The experimental analysis allows us to highlight the performances of the two main pansharpening approaches (i.e. component substitution and multiresolution analysis).

Chavez Jr P S, Sides S C, Anderson J A .

Comparison of three different methods to merge multiresolution and multispectral data: Landsat TM and SPOT panchromatic

[J]. Photogrammetric Engineering and Remote Sensing, 1991,57(3):295-303.

DOI:10.1163/157006707X222768      URL     [本文引用: 2]

The merging of multisensor image data is becoming a widely used procedure because of the complementary nature of various data sets. Ideally, the method used to merge data sets with high-spatial and high-spectral resolution should not distort the spectral characteristics of the high-spectral resolution data. This paper compares the results of three different methods used to merge the information contents of the Landsat Thermatic Mapper (TM) and Satellite Pour l'Observation de la Terre (SPOT) panchromatic data. The comparison is based on spectral characteristics and is made using statistical, visual, and graphical analyses of the results

Laben C A, Brower B V .

Process for enhancing the spatial resolution of multispectral imagery using pan-sharpening:United States, 6011875

[P].2000 -01-04.

[本文引用: 1]

Meng X C, Li J, Shen H F , et al.

Pansharpening with a guided filter based on three-layer decomposition

[J]. Sensors, 2016,16(7):1068.

DOI:10.3390/s16071068      URL     PMID:4970115      [本文引用: 1]

State-of-the-art pansharpening methods generally inject the spatial structures of a high spatial resolution (HR) panchromatic (PAN) image into the corresponding low spatial resolution (LR) multispectral (MS) image by an injection model. In this paper, a novel pansharpening method with an edge-preserving guided filter based on three-layer decomposition is proposed. In the proposed method, the PAN image is decomposed into three layers: A strong edge layer, a detail layer, and a low-frequency layer. The edge layer and detail layer are then injected into the MS image by a proportional injection model. In addition, two new quantitative evaluation indices, including the modified correlation coefficient (MCC) and the modified universal image quality index (MUIQI) are developed. The proposed method was tested and verified by IKONOS, QuickBird, and Gaofen (GF)-1 satellite images, and it was compared with several of state-of-the-art pansharpening methods from both qualitative and quantitative aspects. The experimental results confirm the superiority of the proposed method.

Ranchin T, Wald L .

Fusion of high spatial and spectral resolution images:The ARSIS concept and its implementation

[J]. Photogrammetric Engineering and Remote Sensing, 2000,66(1):49-61.

DOI:10.1002/2014WR016042      URL     [本文引用: 2]

In various applications of remote sensing, when high spatial resolution is required in addition with classification results, sensor fusion is a solution. From a set of images with different spatial and spectral resolutions, the aim is to synthesize images with the highest spatial resolution available in the set and with an appropriate spectral content. Several sensor fusion methods exist; most of them improve the spatial resolution but with a poor quality of the spectral content of the resulting image. Based on a multiresolution modeling of the information, the ARSIS concept (from its French name "Am lioration de la R solution Spatiale par Injection de Structures") was designed in the aim of improving the spatial resolution together with a high-quality in the spectral content of the synthesized images. The general case of application of this concept is described. A quantitative comparison of all presented methods is achieved for a SPOT image. Another example of the fusion of SPOT XS (20 m) and KVR-1000 (2 m) images is given. Practical information for the implementation of the wavelet transform, the multiresolution analysis, and the ARSIS concept by practitioners is given with particular relevance to SPOT and Landsat imagery.

Alparone L, Aiazzi B .

MTF-tailored multiscale fusion of high-resolution MS and Pan imagery

[J]. Photogrammetric Engineering and Remote Sensing, 2006,72(5):591-596.

DOI:10.14358/PERS.72.5.591      URL     [本文引用: 1]

This work presents a multiresolution framework for merging a multispectral image having an arbitrary number of bands with a higher-resolution panchromatic observation. The fusion method relies on the generalized Laplacian pyramid (GLP), which is a multiscale, oversampled structure. The goal is to selectively perform injection of spatial frequencies from an image to another with the constraint of thoroughly retaining the spectral information of the coarser data. The novel idea is that a model of the modulation transfer functions (MTF) of the multispectral scanner is exploited to design the GLP reduction filter. Thus, the interband structure model (IBSM), which is calculated at the coarser scale, where both MS and PAN data are available, can be extended to the finer scale, without the drawback of the poor enhancement occurring when MTFs are assumed to be ideal filters. Experiments carried out on QuickBird data demonstrate that a superior spatial enhancement, besides the spectral quality typical of injection methods, is achieved by means of the MTF-adjusted fusion.

李文静, 温文鹏, 王清和 .

基于Contourlet变换的遥感图像融合方法研究

[J]. 国土资源遥感, 2015,27(2):44-50.doi: 10.6046/gtzyyg.2015.02.07.

URL     Magsci     [本文引用: 1]

<p>为了充分利用多源遥感图像的影像信息,针对不同分辨率的遥感图像进行融合算法研究。通过对基于小波变换(warelet transform,WT)与IHS变换的改进算法研究,提出了基于轮廓波变换(Contourlet transform,CT)与IHS变换的改进算法: 结合传统IHS彩色空间变换,将经IHS变换获得的多光谱图像亮度分量与原全色图像分别进行CT; 然后对得到的低频分量采用自适应融合规则、高频分量采用基于区域相似度的阈值控制规则分别进行融合; 最后对融合后的高频和低频分量进行Contourlet逆变换,得到最终的融合图像。对比实验结果表明: 本文提出的方法能够在有效保留光谱信息的同时,纳入全色图像丰富的空间细节信息。融合之后的结果图像与原多光谱图像具有更高的相关系数和更小的光谱畸变度,并且信息熵和标准差较传统WT及CT更优,具有一定的实用性。</p>

Li W J, Wen W P, Wang Q H .

A study of remote sensing image fusion method based on Contourlet transform

[J]. Remote Sensing for Land and Resources, 2015,27(2):44-50.doi: 10.6046/gtzyyg.2015.02.07.

Magsci     [本文引用: 1]

Tu T M, Su S C, Shyu H C , et al.

A new look at IHS-like image fusion methods

[J]. Information Fusion, 2001,2(3):177-186.

DOI:10.1016/S1566-2535(01)00036-7      URL     [本文引用: 1]

The intensity-hue-saturation (IHS) method, principal component analysis (PCA), Brovey transform (BT) and wavelet transform (WT) are the contemporary image fusion methods in remote sensing community. However, they often face color distortion problems with fused images. In other words, they are sensitive to the characteristics of the analyzed area. To investigate this color distortion problem, this work presents a relatively detailed study indicating that the color distortion problem arises from the change of the saturation during the fusion process. Meanwhile, PCA, BT, and WT are evaluated and found to be IHS-like image merging techniques. Experimental results for distinct image fusion methods are also demonstrated in this paper.

Meng X C, Shen H F, Li H F, et al.

Improving the spatial resolution of hyperspectral image using panchromatic and multispectral images:An integrated method

[C]//Proceedings of the 7th Workshop on Hyperspectral Image and Signal Processing:Evolution in Remote Sensing (WHISPERS).Tokyo,Japan:IEEE, 2015: 1-4.

[本文引用: 1]

孟祥超, 沈焕锋, 张洪艳 , .

基于梯度一致性约束的多光谱/全色影像最大后验融合方法

[J]. 光谱学与光谱分析, 2014,34(5):1332-1337.

URL     [本文引用: 1]

多光谱/全色影像融合可以得到高空间分辨率的多光谱影像, 在影像解译和分类等方面具有十分重要的意义。 提出一种基于梯度一致性约束的遥感影像融合方法。 该方法在最大后验概率框架下, 通过梯度一致性约束建立理想高空间分辨率多光谱影像和全色影像之间的关系, 并结合多光谱影像观测模型和Huber-Markov影像先验, 构建融合目标函数, 最后采用梯度下降法求解得到融合影像。 本文方法在目标函数中引入了梯度一致性约束, 克服了现有的同类方法受限于波段数量的缺陷, 同时在求解中自适应确定每个波段的迭代步长, 充分顾及了各波段的光谱特性, 从而既确保了融合影像的光谱信息保真度, 也提高了融合影像的空间信息融入度。 通过IKONOS和WorldView-2影像对该方法进行了验证, 并和GS, AIHS和AMBF等融合方法从定性和定量两方面进行了比较分析。 实验结果表明, 相比于其他方法, 该方法可以在更好保持光谱信息的同时增强影像的空间分辨率, 具有更广泛的适用范围和更佳的融合效果。

Meng X C, Shen H F, Zhang H Y , et al.

Maximum a posteriori fusion method based on gradient consistency constraint for multispectral/panchromatic remote sensing images

[J]. Spectroscopy and Spectral Analysis, 2014,34(5):1332-1337.

[本文引用: 1]

Meng X C, Shen H F, Zhang L P, et al.

A unified framework for spatio-temporal-spectral fusion of remote sensing images

[C]//Proceedings of IEEE International Geoscience and Remote Sensing Symposium (IGARSS).Milan,Italy:IEEE, 2015: 2584-2587.

[本文引用: 1]

Zhang L P, Shen H F, Gong W , et al.

Adjustable model-based fusion method for multispectral and panchromatic images

[J]. IEEE Transactions on Systems,Man,and Cybernetics,Part B (Cybernetics), 2012,42(6):1693-1704.

DOI:10.1109/TSMCB.2012.2198810      URL     PMID:22736648      [本文引用: 4]

In this paper, an adjustable model-based image fusion method for multispectral (MS) and panchromatic (PAN) images is developed. The relationships of the desired high spatial resolution (HR) MS images to the observed low-spatial-resolution MS images and HR PAN image are formulated with image observation models. The maximum a posteriori framework is employed to describe the inverse problem of image fusion. By choosing particular probability density functions, the fused HR MS images are solved using a gradient descent algorithm. In particular, two functions are defined to adaptively determine most regularization parameters using the partially fused results at each iteration, retaining one parameter to adjust the tradeoff between the enhancement of spatial information and the maintenance of spectral information. The proposed method has been tested using QuickBird and IKONOS images and compared to several known fusion methods using quantitative evaluation indices. The experimental results verify the efficacy of this method.

Ballester C, Caselles V, Igual L , et al.

A variational model for P+XS image fusion

[J]. International Journal of Computer Vision, 2006,69(1):43-58.

DOI:10.1007/s11263-006-6852-x      URL     [本文引用: 1]

We propose an algorithm to increase the resolution of multispectral satellite images knowing the panchromatic image at high resolution and the spectral channels at lower resolution. Our algorithm is based on the assumption that, to a large extent, the geometry of the spectral channels is contained in the topographic map of its panchromatic image. This assumption, together with the relation of the panchromatic image to the spectral channels, and the expression of the low-resolution pixel in terms of the high-resolution pixels given by some convolution kernel followed by subsampling, constitute the elements for constructing an energy functional (with several variants) whose minima will give the reconstructed spectral images at higher resolution. We discuss the validity of the above approach and describe our numerical procedure. Finally, some experiments on a set of multispectral satellite images are displayed.

Palsson F, Sveinsson J R, Ulfarsson M O .

A new pansharpening algorithm based on total variation

[J]. IEEE Geoscience and Remote Sensing Letters, 2014,11(1):318-322.

DOI:10.1109/LGRS.2013.2257669      URL     [本文引用: 1]

In this letter, we present a new method for the pansharpening of multispectral satellite imagery. Pansharpening is the process of synthesizing a high spatial resolution multispectral image from a low spatial resolution multispectral image and a high-resolution panchromatic (PAN) image. The method uses total variation to regularize an ill-posed problem dictated by a widely used explicit image formation model. This model is based on the assumptions that a linear combination of the bands of the pansharpened image gives the PAN image and that a decimation of the pansharpened image gives the original multispectral image. Experimental results are based on two real datasets and the quantitative quality of the pansharpened images is evaluated using a number of spatial and spectral metrics, some of which have been recently proposed and do not need a reference image. The proposed method compares favorably to other well-known methods for pansharpening and produces images of excellent spatial and spectral quality.

Shen H F, Zhang L P, Huang B , et al.

A MAP approach for joint motion estimation,segmentation, and super resolution

[J]. IEEE Transactions on Image Processing, 2007,16(2):479-490.

DOI:10.1109/TIP.2006.888334      URL     PMID:17269640      [本文引用: 1]

Super resolution image reconstruction allows the recovery of a high-resolution (HR) image from several low-resolution images that are noisy, blurred, and down sampled. In this paper, we present a joint formulation for a complex super-resolution problem in which the scenes contain multiple independently moving objects. This formulation is built upon the maximum a posteriori (MAP) framework, which judiciously combines motion estimation, , and super resolution together. A cyclic coordinate descent optimization procedure is used to solve the MAP formulation, in which the motion fields, fields, and HR images are found in an alternate manner given the two others, respectively. Specifically, the gradient-based methods are employed to solve the HR image and motion fields, and an iterated conditional mode optimization method to obtain the fields. The proposed algorithm has been tested using a synthetic image sequence, the "Mobile and Calendar" sequence, and the original "Motorcycle and Car" sequence. The experiment results and error analyses verify the efficacy of this algorithm.

Li H F, Zhang L P, Shen H F , et al.

A variational gradient-based fusion method for visible and SWIR imagery

[J]. Photogrammetric Engineering and Remote Sensing, 2012,78(9):947-958.

DOI:10.14358/PERS.78.9.947      URL     [本文引用: 1]

Abstract This paper presents a new variational gradient-based fusion method for visible and short-wave infrared (SWIR) imagery. The proposed method enables spatial enhancement and dehazing of visible imagery. Integrating gradients from SWIR imagery into visible imagery produces a single image with true color and sharp gradients. A constraint based on band correlation is included to improve the enhancement and implement dehazing. The band correlation is according to the quantitative relationship between the wavelength and the atmospheric effect caused by Rayleigh scattering. In this study, both clear and hazy Landsat ETM[H11001] images are used in the experiments. By visual assessment, the gradient of the fused image is more salient than that of the original image, and the true color is well preserved. With the inclusion of the band correlation constraint, the proposed fusion method yields almost haze-free results. Quantitatively, the Metric Q of the fused images is significantly higher than that of the original images; the largest increase of the Metric Q in the experimental results is from 0.0114 to 0.0611. Moreover, for the results of the proposed method, the Metric Q increase in the visible bands declines from blue band to red band.

Wang Z, Bovik A C .

Mean squared error:Love it or leave it? A new look at signal fidelity measures

[J]. IEEE Signal Processing Magazine, 2009,26(1):98-117.

DOI:10.1109/MSP.2008.930649      URL     [本文引用: 1]

In this article, we have reviewed the reasons why we (collectively) want to love or leave the venerable (but perhaps hoary) MSE. We have also reviewed emerging alternative signal fidelity measures and discussed their potential application to a wide variety of problems. The message we are trying to send here is not that one should abandon use of the MSE nor to blindly switch to any other particular signal fidelity measure. Rather, we hope to make the point that there are powerful, easy-to-use, and easy-to-understand alternatives that might be deployed depending on the application environment and needs. While we expect (and indeed, hope) that the MSE will continue to be widely used as a signal fidelity measure, it is our greater desire to see more advanced signal fidelity measures being used, especially in applications where perceptual criteria might be relevant. Ideally, the performance of a new signal processing algorithm might be compared to other algorithms using several fidelity criteria. Lastly, we hope that we have given further motivation to the community to consider recent advanced signal fidelity measures as design criteria for optimizing signal processing algorithms and systems. It is in this direction that we believe that the greatest benefit eventually lies.

Wald L.

Quality of high resolution synconfproced images: Is there a simple criterion?

[C]//Proceedings of the 3rd Conference “Fusion of Earth Data:Merging Point Measurements,Raster Maps and Remotely Sensed Images”. Sophia Antipolis,France:SEE/URISCA, 2000: 99-103.

[本文引用: 1]

Choi J, Yu K, Kim Y .

A new adaptive component-substitution-based satellite image fusion by using partial replacement

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011,49(1):295-309.

DOI:10.1109/TGRS.2010.2051674      URL     [本文引用: 1]

Preservation of spectral information and enhancement of spatial resolution are regarded as important issues in remote sensing satellite image fusion. In previous research, various algorithms have been proposed. Although they have been successful, there are still some margins of spatial and spectral quality that can be improved. In addition, a new method that can be used for various types of sensors is required. In this paper, a new adaptive fusion method based on component substitution is proposed to merge a high-spatial-resolution panchromatic (PAN) image with a multispectral image. This method generates high-/low-resolution synthetic component images by partial replacement and uses statistical ratio-based high-frequency injection. Various remote sensing satellite images, such as IKONOS-2, QuickBird, LANDSAT ETM+, and SPOT-5, were employed in the evaluation. Experiments showed that this approach can resolve spectral distortion problems and successfully conserve the spatial information of a PAN image. Thus, the fused image obtained from the proposed method gave higher fusion quality than the images from some other methods. In addition, the proposed method worked efficiently with the different sensors considered in the evaluation.

Wald L, Ranchin T, Mangolini M .

Fusion of satellite images of different spatial resolutions:Assessing the quality of resulting images

[J]. Photogrammetric Engineering And Remote Sensing, 1997,63(6):691-699.

DOI:10.1016/S0924-2716(97)00008-7      URL     [本文引用: 2]

Methods have been proposed to produce multispectral images with enhanced spatial resolution using one or more images of the same scene of better spatial resolution. Assuming that the main concern of the user is the quality of the transformation of the multispectral content when increasing the spatial resolution, this paper defines the properties of such enhanced multispectral images. It then proposes both a formal approach and some criteria to provide a quantitative assessment of the spectral quality of these products. Five sets of criteria are defined. They measure the performance of a method to synthesize the radiometry in a single spectral band as well as the multispectral information when increasing the spatial resolution. The influence of the type of landscape present in the scene upon the assessment of the quality is underlined, as well as its dependence with scale. The whole approach is illustrated by the case of a SPOT image and three different standard methods to enhance the spatial resolution

/

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编:100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发