基于HIS变换的遥感影像融合系统
项目介绍
本系统是一款基于MATLAB开发的遥感影像处理工具,专门用于实现多光谱影像(MS)与全色影像(PAN)的融合。在卫星遥感领域,由于传感器物理特性的限制,通常难以同时获得高空间分辨率和高光谱分辨率的影像。本系统通过经典的HIS(强度、色调、饱和度)色彩变换算法,将多光谱影像的光谱信息与全色影像的结构细节进行有效结合。该技术能够生成既具有全色影像高空间分辨率细节,又保留了多光谱影像颜色特征的融合影像,广泛应用于资源调查、城市规划及环境监测等领域。
功能特性
- 遥感数据仿真模拟:内置数据生成模块,通过对标准图像进行下采样和对比度增强,模拟真实遥感场景下的多光谱波段与全色波段。
- 线性对比度拉伸:采用2%至98%的分位数截断拉伸技术,剔除影像表现中的异常极值,增强原始数据的视觉层次感。
- 亚像素级图像重采样:利用双三次插值算法(Bicubic Interpolation)将低分辨率的多光谱影像放大至全色影像的尺寸,确保空间位置的精确对齐。
- IHS色彩空间解耦:通过数学模型将RGB影像转换为IHS空间,将代表能量的强度分量(I)与代表色彩的色度(H)、饱和度(S)分量彻底分离。
- 全色影像直方图匹配:为了减小融合过程中的光谱扭曲,系统能够将全色影像的亮度分布自动调节至与多光谱强度分量一致。
- 分量替换融合:使用调整后的全色影像替换原始强度分量,结合原始色彩分量逆转回RGB空间,实现影像增强。
- 多维度可视化评估:系统自动生成包含原始图、处理中间过程及最终融合结果的对比矩阵,便于直观评估算法效果。
实现逻辑与算法细节
系统的运行流程严格遵循遥感影像融合的标准物理模型,具体逻辑如下:
首先,系统通过对输入图像进行处理,模拟产生一组配准后的遥感数据。全色影像模拟高空间分辨率传感器,经过锐化处理以突出边缘;多光谱影像通过下采样模拟低分辨率传感器。
在预处理阶段,系统对每一路波段执行线性拉伸。通过计算像素强度的累积分布函数,定位2%和98%的增益边界,将影像动态范围归一化到标准化的[0, 1]区间。随后,多光谱影像被重采样,使其行列数与全色影像完全一致。
核心变换算法基于三角形IHS模型。强度分量I是通过计算RGB三通道的平均值获得的。饱和度S基于像素点中最小色彩分量相对于总亮度的占比计算。色调H则利用几何投影公式,根据像素在色彩空间中的偏转角度确定,并归一化至[0, 1]区间以方便计算。
在融合关键步骤中,系统调用直方图匹配算法,以多光谱影像的I分量为参考,对全色影像进行非线性灰度映射。这一步是防止光谱偏移的核心,它确保了全色影像在替换I分量时,其平均亮度和对比度与原图保持一致。
最后,系统执行逆向色彩变换。该变换过程根据色调H所属的扇区(0-120度、120-240度、240-360度)分别应用不同的三角函数公式,计算得出新的R、G、B分量,并对超出范围的像素值进行边界截断处理,得到最终的融合图像。
使用方法
- 准备环境:确保计算机已安装MATLAB R2016a或更高版本,并配置好图像处理工具箱(Image Processing Toolbox)。
- 运行系统:打开MATLAB软件,将工作目录切换至本项目文件夹。
- 执行程序:在命令行窗口输入入口函数名称并回车。
- 查看结果:系统将自动弹出图形化界面,展示从原始数据、预处理状态、中间特征量到最终融合图的全流程视图。
系统要求
软件平台:MATLAB 2016b 或以上版本。
必要工具箱:Image Processing Toolbox(用于执行灰度调整、重采样及直方图匹配)。
硬件要求:建议内存4GB以上,以支持高分辨率遥感矩阵的运算。