ERDAS 影像融合方法汇总影像融合是指将多源信道所采集到的关于同一目标的影像数据经过影像处理和计算机技术等,最大限度的提取各自信道中的有利信息,最后综合成高质量的影像,以提高影像信息的利用率、改善计算机解译精度和可靠性、提升原始影像的空间分辨率和光谱分辨率,利于监测。
ERDAS IMAGINE 提供多种影像融合方法,且支持带RPC模型的影像融合处理。
1、(14个波段(蓝/像也为Step4: 利用低通滤波器(5×5)对多光谱影像进行滤波处理,输出多光谱滤波影像。
Step5: 利用相减法对全色锐化影像、全色滤波影像进行相减处理,并将全色锐化影像按照权重融合到多光谱滤波影像各个波段,输出新多光谱影像。
(2)参数说明Input Sensor Type 待融合影像传感器类型,分为Quickbird、IKONOS、Format三种传感器。
输入影像要求:多光谱和全色分辨率比为4:1、同时获取、为同一个传感器,全色为单波段,多光谱为4波段。
Sharpening Filter Center Value3×3锐化卷积窗口中心值,其他值都为-1,默认值根据传感器变化,范围值为11、14、17、20、23、1000,小的中心值会产生更好的锐化效果。
一般来说,大分辨率影像锐化程度要求低,小分辨率影像锐化程度要求高,若全色影像已经经过锐化处理,此处选择1000。
Pan Contribution Weight融合时全色图像所占的比重(权重),范围为0.7-1.3,默认值根据传感器变化,小的锐化值会产生更好的锐化效果。
Create image of subset area根据子区的坐标来定义融合影像范围。
Create image of full area输出所有区域的融合影像,这个范围是全色和多光谱影像的交集(intersect)。
该项勾选时才可设置融合影像成果名称和存放路径。
Null Value 设置输出图像空值的数值。
Mask input Null Values 勾选该项时,可设置输出图像空值。
2、HPF Resolution Merge 高通滤波融合(1)算法说明高通滤波融合采用高通滤波算法(High Pass Filter, HPF)来实现影像融合。
一副影像由不同频率的成分组成,根据图像频谱概念,高的空间频率对应影像中灰度急剧变化的部分,低的频率代表图像中灰度缓慢变化的部分。
高频分量包含影像的空间结构,低频分量包含影像的光谱信息。
高通滤波算法采用高通卷积滤波算子(high pass convolution filter kernel,HPK)来提取高空间分辨率影像中的空间信息,采用像元相加的方法将其加到低空间分辨率的多光谱影像上,实现影像融合目的。
高通滤波算法步骤:Step1: 计算R(比率)。
Step2: 创建高通卷积滤波算子对高分辨率影像进行高通滤波输出高通滤波影像。
Step3: 多光谱影像重采样,分辨率与高通滤波影像一致。
Step4: 将高通滤波影像按照权重融合到多光谱各个波段,输出新多光谱影像。
Step5: 将新多光谱影像进行拉伸处理,使其与原始多光谱影像均值和标准偏差一致。
(2)参数说明R多光谱影像分辨率与高分辨率影像分辨率的比率,自动计算(直接影响后续参数的选择),该值应该大于1。
Kernel Size高通滤波(High Pass Filter, HPF)卷积算子窗口大小。
根据R值来确定、默认即可,用于高分辨率影像处理。
Center Value高通滤波算子窗口中心值,该窗口其他值为-1。
根据R值来确定、默认即可,用于高分辨率影像处理。
Center value /HPK Size options based on HPK3()IHS彩色空间的变换成为IHS变换。
IHS彩色空间中,亮度I反映图像中地物反射的全部能量和图像所包含的空间信息,代表影像的地面分辨率;色度H指组成色彩的主波长,由红绿蓝的比重所决定;饱和度S代表颜色的纯度,H/S代表影像的光谱分辨率。
改进IHS融合采用全色高分辨率影像替换I分量从而使影像同时具有高空间分辨率和高光谱分辨率,融合步骤如下:Step1: 将多光谱影像进行IHS变换,从RGB空间转换至IHS空间。
Step2: 将全色影像与分量I进行直方图匹配,输出新全色影像。
Step3: 将新全色影像作为分量I与分量H、S形成新的IHS图像。
Step4: 将新的IHS图像进行IHS逆变换转换成RGB空间,输出融合后多光谱影像。
改进IHS融合克服了传统IHS融合只能输出三个波段影像的问题。
比如:Quickbird多光谱包含4个波段,可采用432、321组合进行两次迭代融合输出4321融合影像。
(2)参数说明Clip Using Min/Max用多光谱像元的最大最小值来规定重采样后多光谱影像数据的像元值范围,重采样方法选择三次卷积时生效,因为其他两种重采样方法不会超过原来的像元值范围。
波段列表中的字段:Low\High是该波段的波长最低值和最高值,Center是该波段的波长中心值(均值),Range是该波段的波长跨度(High-Low),Overlp是该波段波长范围与全色波长的重叠范围,Beta=Overlap/Pan_Range即重叠范围与全色波长范围的比值。
Resampling Technique重采样方法,提供最邻近、双线性内插和三次卷积,默认双线性内插法。
影像融合过程中重采样方法的选择:根据全色多光谱的相对分辨率来确定。
1) 最近邻法:边缘效应比较强的影像采用最邻近法输出的融合影像仍旧存在。
2) 双线性内插法:全色分辨率较低、多光谱全色分辨率比率太小为2:1甚至更低时(如Landsat全色15m、多光谱30m)采用,效果较好、处理速度较快。
3) 三次卷积法:全色分辨率较高,多光谱全色分辨率必须为4:1甚至更高时(如IKONOS全色1m、多光谱4m)采用,能输出较好的重采样结果。
某些情况下,采用三次卷积法重采样会造成结果错误,这是由于该采样方法会改变多光谱的最大最小像元值范围造成的。
实际应用中,对于大文件数据的影像融合,考虑处理速度和结果的正确性,推荐使用双线性内插法;如果数据量较小,可尝试使用三次卷积法。
Sensor分别选择高分辨影像、多光谱影像的传感器类型。
支持的传感器类型包括:EO-1 ALI、FORMOSAT-2、GeoEye-1、IKONOS、IRS、Landsat 7/8/8 Cirrus、OrbView-3、Pleiades、QuickBird、SkySat、SPOT 4/6/7、WorldView-1/2/3。
Sensor分别选择高分辨影像、多光谱影像的传感器类型。
支持的传感器类型包括:EO-1 ALI、FORMOSAT-2、GeoEye-1、IKONOS、IRS、Landsat 7/8/8 Cirrus、OrbView-3、Pleiades、QuickBird、SkySat、SPOT 4/6/7、WorldView-1/2/3。
Ratio Ceiling设置亮度修正系数的上限。
Batch 支持影像融合批处理操作。
AOI 支持导入AOI文件输出感兴趣区域内的影像融合结果。
Layer Combination Method多光谱组成方法,默认True Color,提供多种选择方法: True Color(真彩色)、False Color IR(假彩色)、Desktop RGB(真彩色)、Desktop BRG(真彩色)、Preference(推荐设置)和Custom(自定义组合)。
Computation Method 计算方法,默认是Single pass – 3 layer RGB(只用所选的三个波段进行输出图像计算),还提供Iterate on multiple RGB combinations(选择多余三个波段进行输出图像计算)输出多波段融合影像。
Add to Iteration list 将当前组合方式添加至迭代列表,当计算方式为Iterate on multiple RGB combinations时可用。
Input to Output Layer Mapping RGB空间到IHS空间转换的波段组合。
Output Files影像融合结果输出路径和名称。
Number of layers影像融合结果波段数。
Data Type 可查看原始影像数据类型、设置输出成果数据类型。
Ignore Zero in Output Statistics勾选时,影像融合结果统计计算时忽略0值。
Ignore Zero in Raster Match勾选时,影像匹配计算时忽略0值。
4、Wavelet Resolution Merge 小波变换融合(1)算法说明小波变换类似于傅里叶变换,不同的是小波变换基于一些小型波,将图像分解成频率域上各个频率上的子图像,来代表原始图像的各个特征分量,因此小波变换融合可根据不同的特征分量采用不同的融合方法来达到最佳的融合效果。
一副图像的小波分解中,绝对值较大的小波高频系数对应亮度急剧变换的点,即图像中对比度变化较大的边缘特征,如边界、亮线及区域轮廓,小波变换融合就是保留此类显著的特征,比如融合前图像A中比图像B中显著,融合后图像中A中的目标被保留,B中的目标被忽略。
小波变换就是基于不同图像中A、B目标的小波变换系数在不同分辨率水平上的主要地位来确保图像中A、B 显著目标都保留。
(2)参数说明Spectral Transform多光谱影像转换为单波段灰度图的方法,提供三种方法:Single Band:选择指定波段,可自定义使用波段数为任意波段(1,2,3,4等)。
IHS:对多光谱进行IHS变换,使用亮度I分量进行融合,仅输出三波段融合成果。
Principal Component:对多光谱进行主成分变换,使用第一主成分进行融合。
Layer Selection 设置用于影像融合的多光谱影像波段。
Resampling Techniques重采样方法,提供Nearest Neighbor(最邻近)和Bilinear Interpolation(双线性内插)两种。
Data Type 可查看原始影像数据类型、设置输出成果数据类型。
Stretch to Unsigned 8 bit 勾选选项时,output无法设置输出类型,输出文件的像元范围拉伸到0~255范围内,并提供两种拉伸方法: Standard Deviation stretch(标准差拉伸)、Min-Max stretch(最大最小值拉伸)。
一般不推荐在融合时进行降位处理。
Ignore Zero in Stats.勾选选项时,影像融合结果统计计算时忽略0值。