基于小波变换的图像融合方法及其系统技术方案

一种基于小波变换的图像融合方法，包括：离散小波变换、高频图像处理、低频图像处理以及分解图像融合。本发明专利技术方法能够对源图像的高、低频图像进行分别处理，可保留更多的有用信息，能够得到质量更好的融合图像。

【技术实现步骤摘要】
【专利摘要】一种基于小波变换的图像融合方法，包括：离散小波变换、高频图像处理、低频图像处理以及分解图像融合。本专利技术方法能够对源图像的高、低频图像进行分别处理，可保留更多的有用信息，能够得到质量更好的融合图像。【专利说明】基于小波变换的图像融合方法及其系统
本专利技术属于图像融合
，尤其涉及一种基于小波变换的图像融合方法及其 系统。
技术介绍
人们出于不同的目的有时会对于同一个场景拍摄多张图片。图像融合技术可以用 来对于同一场景的多张图片进行处理，融合成为一张质量更好的图片，以达到通过一张图 片获取这一场景更多更有效信息的目的。图像融合技术可以用在生活中的各个领域，比如 遥感、多焦距融合、医疗图像融合等等，与人们的生活息息相关。 图像融合技术已经存在了许多年，现有的技术比如HIS、PCA融合技术已经可以满 足日常需求。但是，这两种方法都有各自的不足之处，在有些场景下不能够显示出足够多的 重要细节部分，比如两种医疗图像PET与MRI的融合时不能同时保留两种图像的优势细节部 分，为了实现更好的融合效果，融合技术还是在不断改进以满足更高的要求。
技术实现思路
基于此，针对上述技术问题，提供一种基于小波变换的图像融合方法及其系统。 为解决上述技术问题，本专利技术采用如下技术方案： -种基于小波变换的图像融合方法，包括： 离散小波变换:通过离散小波变换算法将同一场景的多张图像分解成多张高频图 像以及多张低频图像； 高频图像处理:在所述多张高频图像中，取对应相同空间位置的多个像素点中像 素绝对值最大的像素点的值作为高频融合图像在该空间位置的像素点的值，获得高频融合 图像在每个空间位置的像素点的值，从而获得高频融合图像； 低频图像处理:对多张低频图像中对应相同空间位置的多个像素点进行加权平均 计算，将加权平均值作为低频融合图像在该空间位置的像素点的值，获得低频融合图像在 每个空间位置的像素点的值，从而获得低频融合图像； 分解图像融合:通过逆小波变换算法将所述高频融合图像以及低频图像融合成一 张图像。 所述离散小波变换步骤通过matlab的DWT函数进行1-3级小波变换。 所述高频图像处理步骤还包括:在所述每个空间位置的像素点的周围进行3*3的 窗口搜索，若某一像素点周围的8个点都取自同一张高频图像，则将该像素点的值校正为该 高频图像中相应像素点的值。 本方案还涉及一种基于小波变换的图像融合系统，包括： 离散小波变换单元，用于通过离散小波变换算法将同一场景的多张图像分解成多 张高频图像以及多张低频图像； 高频图像处理单元，用于在所述多张高频图像中，取对应相同空间位置的多个像 素点中像素绝对值最大的像素点的值作为高频融合图像在该空间位置的像素点的值，获得 高频融合图像在每个空间位置的像素点的值，从而获得高频融合图像； 低频图像处理单元，用于对多张低频图像中对应相同空间位置的多个像素点进行 加权平均计算，将加权平均值作为低频融合图像在该空间位置的像素点的值，获得低频融 合图像在每个空间位置的像素点的值，从而获得低频融合图像； 分解图像融合单元，用于通过逆小波变换算法将所述高频融合图像以及低频图像 融合成一张图像。 所述离散小波变换单元通过matlab的DWT函数进行1-3级小波变换。 所述高频图像处理单元还用于在所述每个空间位置的像素点的周围进行3*3的窗 口搜索，若某一像素点周围的8个点都取自同一张高频图像，则将该像素点的值校正为该高 频图像中相应像素点的值。 本专利技术方法能够对源图像的高、低频图像进行分别处理，可保留更多的有用信息， 能够得到质量更好的融合图像。【附图说明】 下面结合附图和【具体实施方式】本专利技术进行详细说明： 图1为本专利技术的一种基于小波变换的图像融合方法的流程图； 图2为本专利技术的一种基于小波变换的图像融合系统的结构示意图。【具体实施方式】 如图1所示，一种基于小波变换的图像融合方法，包括： S110、离散小波变换:通过离散小波变换算法将同一场景的多张图像分解成多张 高频图像以及多张低频图像。 具体地，本实施例通过matlab的DWT函数进行1-3级小波变换，可取得较好的效果。 每增加一级小波变换可以分解出更过的高频图像，以及得到更少冗余信息的低频 图像，在一定程度上可以融合出质量更好的图像;但是如果级数过多，那么低频图像将会丢 失过多的图像信息，合成后的图像可能会损失原图的信息。所以，最优的级数在1-3级之间。 S120、高频图像处理： 由于高频图像包含的是图像的边缘等信息，所以应该予以保留。 在多张高频图像中，取对应相同空间位置的多个像素点中像素绝对值最大的像素 点的值作为高频融合图像在该空间位置的像素点的值，从而获得高频融合图像在每个空间 位置的像素点的值，进而获得高频融合图像。 为了提高融合图像的质量，本步骤还可在每个空间位置的像素点的周围进行3*3 的窗口搜索，若某一像素点周围的8个点都取自同一张高频图像，则将该像素点的值校正为 该高频图像中相应像素点的值。 若8个点不是取自同一张高频图像，则保留原值。 S130、低频图像处理： 对多张低频图像中对应相同空间位置的多个像素点进行加权平均计算，将加权平 均值作为低频融合图像在该空间位置的像素点的值，从而获得低频融合图像在每个空间位 置的像素点的值，进而获得低频融合图像。由于低频段图像有很多是冗余信息，故这里采用加权平均源图片的每个像素点， 以此来得到融合图像的像素值。现以两张图X、Y的融合为例进行说明：每个低频段的图像系数使用公式（1)计算，D代表系数值，W代表权值，Z是融合后的 图像，P为像素点： Dz(p) =wx(p)Dx(p)+wy(p)Dy(p) (1); 计算融合等级A，假设AX〈AY，W是加权平均窗口，使用公式(2)计算A，其中，w(s，t) 代表窗口的权值，并且每个窗口权值的和为l，s和t代表水平和竖直方向的下标，这里的融 合等级A使用窗内每一个系数的高斯区域方差来表示，这意味着窗内越靠近p点的系数的权 值越大，m，η，k，1为高斯窗口的系数： Ai(p)= 2ses,terw(s，t) |Di(m+s，n+t，k，l) I， （2); 计算用来得到权值的度量值MXY，使用公式(3)，参数由公式(1)、（2)计算得到：牌 若MXY大于或者等于一个阈值α，那么权值用公式(4)计算：(4)<·· 若MXY小于阈值α，那么WY = I，WX = 0。 S140、分解图像融合:通过逆小波变换算法将高频融合图像以及低频图像融合成 一张图像。 本专利技术方法能够对源图像的高、低频图像进行分别处理，可保留更多的有用信息， 能够得到质量更好的融合图像。 如图2所示，本方案还涉及一种基于小波变换的图像融合系统，包括： 离散小波变换单元110,用于通过离散小波变换算法将同一场景的多张图像分解 成多张高频图像以及多张低频图像。 具体地，本实施例通过matlab的DWT函数进行1-3级小波变换，可取得较好的效果。 每增加一级小波变换可以分解出更过的高频图像，以及得到更少冗余信息的低频 图像，在一定程度上可以融合出质量更好的图像;但是如果级数过多，那么低频图像将会丢 失过多的图像信息，合成后的图像可能会损失原图的信息。所以，最优的级数本文档来自技高网...

【技术保护点】
一种基于小波变换的图像融合方法，其特征在于，包括：离散小波变换：通过离散小波变换算法将同一场景的多张图像分解成多张高频图像以及多张低频图像；高频图像处理：在所述多张高频图像中，取对应相同空间位置的多个像素点中像素绝对值最大的像素点的值作为高频融合图像在该空间位置的像素点的值，获得高频融合图像在每个空间位置的像素点的值，从而获得高频融合图像；低频图像处理：对多张低频图像中对应相同空间位置的多个像素点进行加权平均计算，将加权平均值作为低频融合图像在该空间位置的像素点的值，获得低频融合图像在每个空间位置的像素点的值，从而获得低频融合图像；分解图像融合：通过逆小波变换算法将所述高频融合图像以及低频图像融合成一张图像。

【技术特征摘要】

【专利技术属性】
技术研发人员：宋高飞，张伟，张彩红，谯帅，刘洋，
申请(专利权)人：博康智能网络科技股份有限公司，
类型：发明
国别省市：上海;31