Landsat 8 OLT遥感数据融合方法及其质量评价

2018-03-06军赵绒绒王婷婷

现代测绘 2018年1期

范军赵绒绒王婷婷

(1.南京市国土资源信息中心,江苏南京210005;2.江苏省地质测绘院,江苏南京211102;3.南京林业大学土木工程学院,江苏南京210037)

0 引言

多种对地观测遥感卫星的发射成功、遥感技术快速发展,多传感器、多时相、多光谱、多空间分辨率遥感影像数据迅猛增加提供了从粗到精、从多光谱到细光谱段的遥感对地观测数据源[1]。在遥感领域中,单一传感器影像数据通常不能提取足够多的信息去满足一些应用需要,它所提取的信息是有限的,不完整的[2]。因此,如何将这些海量遥感数据作为一个整体来综合应用,克服单一遥感影像提供信息不足带来的问题,从而有效利用不同分辨率、不同时相、不同传感器的遥感海量数据信息便成为当前遥感影像处理的一个热门研究课题。

图像融合不是简单的信息叠加,它是要产生新的、包含更多、更有价值的信息的图像,即达到1+1＞2的效果[3]。具有高空间分辨率的全色波段和具有丰富光谱信息的多光谱波段的数据融合,可得到保持多光谱图像光谱特征、空间细节表现能力的融合图像[4]。本文以Landsat 8 OLI全色影像和多光谱影像为数据源(轨道号分别为P122R41影像的B652、B764、B564组合和P125R47影像的B251组合),采用乘积变换(Brovey)、HSV变换、主成分分析(PCA)和光谱锐化法(Gram-Schmidt,GS)4种方法对多光谱影像和全色影像进行融合,运用Matlab软件计算提取融合后图像的均值、标准差、信息熵、平均梯度和相关系数等定量评价指标,结合主观和客观因素评价数据融合效果。

1 数据源与研究方法

1.1 数据源及其预处理

Landsat 8卫星于2013年2月发射,按近极点太阳同步轨道绕地球飞行,轨道高度705 km,倾角98.2°,每98.9 min绕地球一圈,每16 d覆盖地球一遍,降交点时间为当地上午10:00—10:15,卫星数据下行速率为441 Mbps,携带陆地成像仪(Operational Land Imager,OLI)和热红外传感器(Thermal Infrared Sensor,TIRS)[5]。OLI包含ETM+传感器所有波段,为避免大气吸收特征,重新调整了Band5(0.845～0.885μm),排除0.82μm处水汽吸收特征。同时,OLI全色波段Band8范围较窄,增强了植被和无植被特征区分能力;传感器新增两个波段:用于海岸带监测蓝色波段(Band1:0.433～0.453μm)、云检测短波红外波段(Band9:1.360～1.390μm)[6]。本文选择轨道号分别为P122R41、P125R47的Landsat 8 OLI遥感影像为数据源,对其实现几何纠正和大气校正等预处理,完成遥感数据融合,增进对不同融合方法的理解和认识。

1.2 遥感数据融合方法

有效的遥感数据融合算法可尽量保持多光谱影像的光谱信息,提高多光谱影像的空间分辨率,满足实际应用需求。同一种融合方法针对同一影像,观察者感兴趣的部分不同,对融合效果的质量评价也不同;针对不同应用方向,对影像的各项参数要求也不同,需采用不同融合方法[7-9]。因此,本文选择Brovey、HSV、PCA和Gram-Schmidt四种融合方法对Landsat 8 OLI遥感数据进行处理(表1),分析数据质量是否达到预期效果[10]。

表1 不同遥感数据融合方法表达及融合效果

1.3 数据融合质量评价

1.3 .1 主观评价

主观评价建立在人的视觉感官和经验知识基础上,通过人眼对图像的感知能力判断图像质量[11],但其缺点在于不具备客观性,所得到的评价结果可能存在失误。人眼对图像比较明显的感知体现为:图像是否有重影、不同地物分界是否清晰、颜色差别是否很大等。从亮度、清晰度、对比度等方面对融合图像的观察,需在初步主观目视判别基础上进行客观定量评价。

1.3 .2 客观评价

图像融合结果客观定量评价是采用一定量化指标来评价融合效果的,克服了人为主观影响,对图像的评价更加科学、客观[12-14]。常见的客观评价指标主要分为3大类,即影像的亮度信息、空间细节信息、光谱信息量及融合结果的相关类评价参数。本文主要运用Matlab编程计算Landsat OLI 8影像融合前和融合后的均值、标准差、平均梯度、信息熵和相关系数作为客观定量评价指标。

2 结果分析

2.1 B652数据融合及其质量评价

在空间分辨率和清晰度上,利用4种融合方法的融合影像相较原始影像有所改善:融合影像纹理特征更为清晰,保留了尽可能多的原始信息,提高了地物细节特征,更加易于判读。同时,也改善了影像的光谱特征,各类地物光谱畸变较小,保留了原始影像的光谱信息。Brovey融合影像的颜色更鲜艳,绿色更接近于植被真实颜色;HSV融合影像中,裸地部分与植被部分界限十分清晰,但植被部分太暗,难以识别其具体地物;高通滤波融合影像的空间清晰度和PCA影像最接近(图1)。

图1 4种融合算法的融合影像与原始多光谱影像

将4种融合方法提取的定量评价指标列举如下(表2)。

从表2可知:①HSV的亮度最低,其他3种变化RGB值都差不多;②HSV变换后的影像信息量丰富,即经过HSV变换后图像不仅人眼可以识别较多信息,且实际光谱信息也相对较多;③HSV变换明显高于其他任何变换指标,图像细节信息得到明显增强;④ 高通滤波融合算法综合较高,所得的图像与原图像匹配度较高,光谱保真能力较强,而HSV变换较小,对原图像的光谱保真能力较差;⑤Brovey变换、PCA变换和GS变换的效果基本符合主观评价的结果。

表2 不同融合方法提取的定量评价指标

波段组合B652的图像植被类型较丰富,裸地信息明显,可与有作物的耕地区分。HSV变换后图像反差较大,影像层次较多,可分性大,但是HSV融合后第5波段的相关系数较低,植被部分光谱保真较差,会降低后期分类精度,不建议采用。Brovey变换要优于Gram-Schmidt光谱锐化变换和PCA变换,可分性较大,且主观分析图片质量较好,更适合被用作植被解译及信息提取。

2.2 B764数据融合及其质量评价

从影像的空间分辨率、清晰度来看,4种融合算法所得融合影像的空间分辨率、清晰度都较原始多光谱影像有很大提高。从光谱特征看,HSV亮度适中,且建筑物及裸地部分与植被部分界限明显,其他融合影像也均具有较丰富的色彩,但融合后的图像相对较暗,GS相对偏蓝(图2)。

4种融合方法提取的定量评价指标如下(表3)。

图2 4种融合算法的融合影像与原始多光谱影像

表3 不同融合方法提取的定量评价指标

从表3可知:①Brovey、GS和PCA融合虽然目视颜色较丰富,但均值略偏低;② HSV变换信息量丰富,即经过HSV变换后图像具有最多光谱信息;③HSV变换平均梯度明显高于其他任何变换指标,融合后图像细节信息得到明显增强,影像较清晰;④4种融合结果相关系数相差不大,光谱保真力区分不明显;⑤Brovey变换、PCA变换和GS变换在空间信息详细度和影像清晰度上融合效果很接近,HSV融合光谱信息量和清晰度相对较好。

HSV融合后图像反差较大,影像层次较多,建筑物和植被分界明显,可分性较大,且其光谱保真较好:无论从定性还是定量方面,都利于建筑物解译和信息提取。其他3种融合无论是从主观还是客观评价,其效果都比较相似。

2.3 B564数据融合及其质量评价

B564适合用于土壤湿度和植被状况的分析、陆地水体和陆地/水体分界线确定以及山地、丘陵、平原台地等喀斯特地貌景观及各类用地影像特异分异。B564波段合成非标准假彩色图像,红外波段与红色波段合成,水体边界清晰,利于海岸识别;植被有较好显示,但不便于区分具体植被类别。采取5、6、4波段分别赋予红、绿、蓝色合成图像,色彩反差明显,层次丰富,且各类地物的色彩显示规律与常规合成片相似,符合过去常规片的目视解译(图3)。

将4种融合方法提取的定量评价指标列举如下(表4)。

图3 4种融合算法的融合影像与原始多光谱影像

表4 不同融合方法提取的定量评价指标

通过上述4种融合方法,融合影像在空间分辨率和清晰度上均有所改善,融合影像纹理特征变得清晰,保留了尽可能多的原始信息,提高了地物细节特征,更加易于判读。同时,也改善了影像的光谱特征,各类地物光谱畸变较小,保留了原始影像的光谱信息。Brovey融合影像具有较丰富的色彩,HSV融合影像中,水体部分与陆地部分界限十分清晰,但是植被偏暗,裸地部分偏亮,不易分析,难以识别其具体地物。高通滤波融合影像的空间清晰度和PCA影像最接近。

从表4可知:①Brovey变换和GS变换后的均值相对较高;②HSV变换融合信息量丰富,即经过HSV变换后,图像具有最多的光谱信息量;③HSV变换明显高于其他任何变换指标,融合后图像的细节信息得到了较高增强,影像比较清晰;④Brovey算法和高通滤波融合算法的保真能力较高,说明此变换得到的图像与原图像的匹配度较高,失真较小,光谱保真能力较强,而HSV变换和PCA变化结果较小,该算法对原图像的光谱保真能力较差;⑤Brovey变换和GS变换这两种算法无论是在光谱保真上还是影像清晰度上,融合效果都很接近。HSV融合和PCA融合光谱保真效果都很差,但光谱信息量很大。

2.4 B251数据融合及其质量评价

本次融合加入新元素海蓝波段(Band1:0.433～0.453μm),主要应用海岸带观测,将其与波段2和波段5组合应用,对波段组合添加方法,可给地物解译添加新内容。本文利用4种融合方法对B251波段与全色波段进行融合和质量评价。上述4种融合方法的融合影像在空间分辨率和清晰度上比原始影像有所改善,融合影像纹理特征变得清晰,保留了尽可能多的原始信息,提高了地物细节特征,更加易于判读,且保留了原始影像的光谱信息。Brovey融合影像略清晰,但是相对PCA影像和GS影像偏暗。HSV融合影像中,植被部分和海洋部分都偏暗,但分界线很明显。高通滤波融合影像和PCA影像在色度和空间清晰度上比较接近(图4)。

将4种同融合方法提取的定量价指标举如下(表5)。

图4 4种融合算法的融合影像与原始多光谱影像

表5 不同融合方法提取的定量评价指标

从表5可知:①Brovey变换、GS变换和PCA变换其图像亮度也很相近,而HSV影像亮度最低;② 各变换标准差有所降低而信息熵却是增加的,这说明经过各影像变换后光谱信息量有所损失也有所添加,其中HSV变化光谱信息损失得最大;③HSV变换明显高于其他任何变换指标,而Brovey变化其次,说明融合后图像的细节信息得到了较高增强,影像也比较清晰;④Brovey算法综合较高,此变换得到的图像与原图像的匹配度较高,失真较小,光谱保真能力较强,而HSV变换的相对系数最小,该算法对原图像的光谱保真能力较差。⑤4种变化中Brovey算法的融合效果相对较好,HSV算法无论是从色度、光谱保真还是信息量都相对较差,其他两种居于二者之间。

Brovey、GS和PCA变换后图像反差较大,影像层次较多,可分性大,利于后期信息提取及地物分类。再对比这3组融合,Brovey变换后海蓝波段1标准差和信息熵都相对较高,且光谱保真较好,符合本次融合的目的,Brovey变换下波段1的优势使其在海岸带遥感图像中成为首选融合方法。而本次融合HSV变换效果较差。