【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】

原创

Matlab科研 2021-11-08 11:04:15 ©著作权

©著作权归作者所有：来自51CTO博客作者Matlab科研的原创作品，请联系作者获取转载授权，否则将追究法律责任

一、小波变换图像融合技术简介

1 案例背景

图像融合，指通过对同一目标或同一场景用不同的传感器(或用同一传感器采用不同的方式)进行图像采集得到多幅图像，对这些图像进行合成得到单幅合成图像，而该合成图像是单传感器无法采集得到的。图像融合所输出的合成图像往往能够保持多幅原始图像中的关键信息，进而为对目标或场景进行更精确、更全面的分析和判断提供条件。图像融合属于数据融合范畴，是数据融合的子集，兼具数据融合和图像可视化的优点。因此，图像融合能够在一定程度上提高传感器系统的有效性和信息的使用效率，进而提高待分析目标的分辨率，抑制不同传感器所产生的噪声，改善图像处理的效果。

图像融合最早是以数据融合理论为基础的，通过计算像素算术平均的方式得到合成图像。该方法忽略了像素间的相互关系，往往会产生融合图像的对比度差、可视化效果不理想等问题。因此，为了提高目标检测的分辨率，抑制不同传感器的检测噪声，本案例选择了一种基于小波变换的图像数据融合方法，首先通过小波变换将图像分解到高频、低频，然后分别进行融合处理，最后再逆变换到图像矩阵"。在融合过程中，为了尽可能保持多源图像的特征，在小波分解的高频域内，选择图像邻域平均绝对值较大的系数作为融合小波重要系数；在小波分解的低频域内，选择对多源图像的低频系数进行加权平均作为融合小波近似系数。在反变换过程中，利用重要小波系数和近似小波系数作为输入进行小波反变换。在融合图像输出后，对其做进一步的处理。实验结果表明，基于小波变换的图像数据融合方法运行效率高，具有良好的融合效果，并可用于广泛的研究领域，具有一定的使用价值。

根据融合的作用对象，图像融合一般可以分为3个层次：像素级图像融合、特征级图像融合和决策级图像融合。其中，像素级融合是作用于图像像素点最底层的融合，本章所研究的图像融合是像素级图像融合。

1.2理论基础

传统的直接像素算术平均进行图像融合的方法往往会造成融合结果对比度降低、可视化效果不理想等问题，为此研究人员提出了基于金字塔的图像融合方法，其中包括拉普拉斯金字塔、梯度金字塔等多分辨率融合方法。20世纪80年代中期发展起来的小波变换技术为图像融合提供了新的工具，小波分解的紧致性、对称性和正交性使其相对于金字塔分解具有更好的图像融合性能。此外，小波变换具有“数学显微镜”聚焦的功能，能实现时间域和频率域的步调统一，能对频率域进行正交分解，因此小波变换在图像处理中具有非常广泛的应用，已经被运用到图像处理的几乎所有分支，如图像融合、边缘检测、图像压缩、图像分割等领域。

假设对一维连续小波w…(t)和连续小波变换W，(a，b)进行离散化，其中，a表示尺度参数， b表示平移参数，在离散化过程中分别取a=a{和b=bf，其中， jeZ， a， >1，则对应的离散小波函数如下：

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_数据

离散化的小波变换系数如下：

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_小波变换_02

小波重构公式如下：

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_小波变换_03

式中，C为常数且与数据信号无关。根据对连续函数进行离散化逼近的步骤，如果选择的a，和b，越小，则生成的网格节点就越密集，所计算的离散小波函数wj x® 和离散小波系数Cj就越多，进而数据信号重构的精确度也越高。

由于数字图像是二维矩阵，所以需要将一维信号的小波变换推广到二维信号。假设(x)是一个一维的尺度函数，p(x)是相应的小波函数，那么可以得到一个二维小波变换的基础函数：

v’(x，y)=0(x)v(y)w(x，y)=y(x)(y)v(x，y)=V(x)w(y)

由于数字图像是二维矩阵，一般假设图像矩阵的大小为NxN，且N=2"(n为非负整数)，所以经一层小波变换后，原始图像便分解为4个分辨率为原来尺寸一的子带区域，如图11-1所示，分别包含了相应频带的小波系数，这一过程相当于在水平方向和垂直方向上进行隔点采样。

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_高通滤波_04

进行下一层小波变换时，变换数据集中在LL子带上。说明了图像小波变换的数学原型。

(1)LL频带保持了原始图像的内容信息，图像的能量集中于此频带：

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_高通滤波_05

对图像进行小波变换的原理就是通过低通滤波器和高通滤波器对图像进行卷积滤波，再进行二取一的下抽样。因此，图像通过一层小波变换可以被分解为1个低频子带和3个高频子带。其中，低频子带LL；通过对图像水平方向和垂直方向均进行低通滤波得到；高频子带HL通过对图像水平方向高通滤波和垂直方向低通滤波得到；高频子带LH通过对图像水平方向低通滤波和垂直方向高通滤波得到：高频子带HH通过对图像水平方向高通滤波和垂直方向高通滤波得到。各子带的分辨率为原始图像的，同理，对图像进行二层小波变换时只对低频子带LL进行，可以将LLi子带分解为LL 2、LH 2、HL 2和HH，各子带的分辨率为原始图像的一。，以此类推可得到三层及更高层的小波变换结果。所以，进行

一层小波变换得到4个子带，进行二层小波变换得到7个子带，进行x层分解就得到3·x+1个子带。如图11-2所示为三层小波变换后的系数分布。

【图像融合】基于matlab GUI小波变换可见光与红外光图像融合【含Matlab源码 701期】_小波变换_06

二、部分源代码

function varargout = main(varargin)
% MAIN MATLAB code for main.fig
%      MAIN, by itself, creates a new MAIN or raises the existing
%      singleton*.
%
%      H = MAIN returns the handle to a new MAIN or the handle to
%      the existing singleton*.
%
%      MAIN('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in MAIN.M with the given input arguments.
%
%      MAIN('Property','Value',...) creates a new MAIN or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before main_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to main_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help main

% Last Modified by GUIDE v2.5 14-Dec-2017 19:28:26

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @main_OpeningFcn, ...
                   'gui_OutputFcn',  @main_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before main is made visible.
function main_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to main (see VARARGIN)

% Choose default command line output for main
handles.output = hObject;
handles.imfusion1=[];
handles.imfusion2=[];
axis(handles.axes1,'off');
axis(handles.axes2,'off');
axis(handles.axes3,'off')
movegui(handles.figure1,'center');
% Update handles structure
guidata(hObject, handles);
uiwait(handles.figure1);

% UIWAIT makes main wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = main_OutputFcn(hObject, eventdata, handles) 
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;


% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
img =handles.imfusion1;
% img = rgb2gray(img);
% img = img(20:160,140:280);
nn=10;
img=imresize(img,[58*nn 77*nn]);


img2 =handles.imfusion2;
% img2 = rgb2gray(img2);
% img2 = img2(20:160,140:280);
img2=imresize(img2,[58*nn 77*nn]);
[RcA1,RcH1,RcV1,RcD1] = dwt2(img,'sym4');%LL1,HL1,LH1,HH1
[VcA1,VcH1,VcV1,VcD1] = dwt2(img2,'sym4');%LL1,HL1,LH1,HH1
[RcA2,RcH2,RcV2,RcD2] = dwt2(RcA1,'sym4');%LL1,HL1,LH1,HH1
[VcA2,VcH2,VcV2,VcD2] = dwt2(VcA1,'sym4');%LL1,HL1,LH1,HH1
%%%%%% 2级小波融合
%%%%%% 策略1：  低频平均；高频模值较大者
[M1,N1,P1] = size(RcA1);
[M2,N2,P2] = size(RcA2);
%%%%%%1级
cA10 = (RcA1 + VcA1)*0.5;
for i=1:M1
    for j=1:N1
        for k=1:3
        if abs(RcH1(i,j,k)) > abs(VcH1(i,j,k))
              cH10(i,j,k) = RcH1(i,j,k);      
        else
              cH10(i,j,k) = VcH1(i,j,k); 
        end
        if abs(RcV1(i,j,k)) > abs(VcV1(i,j,k))
              cV10(i,j,k) = RcV1(i,j,k);      
        else
              cV10(i,j,k) = VcV1(i,j,k); 
        end
        if abs(RcD1(i,j,k)) > abs(VcD1(i,j,k))
              cD10(i,j,k) = RcD1(i,j,k);      
        else
              cD10(i,j,k) = VcD1(i,j,k); 
        end        
        end
    end
end
%%%%%%%2级
cA20 = (RcA2 + VcA2)*0.5;
for i=1:M2
    for j=1:N2
        for k=1:3
        if abs(RcH2(i,j,k)) > abs(VcH2(i,j,k))
              cH20(i,j,k) = RcH2(i,j,k);      
        else
              cH20(i,j,k) = VcH2(i,j,k); 
        end
        if abs(RcV2(i,j,k)) > abs(VcV2(i,j,k))
              cV20(i,j,k) = RcV2(i,j,k);      
       
              cD20(i,j,k) = RcD2(i,j,k);      
        else
              cD20(i,j,k) = VcD2(i,j,k); 
        end        
        end
    end
end
%%%%%%反变换
img_fuse2 = idwt2(cA20,cH20,cV20,cD20,'sym4');%LL1,HL1,LH1,HH1
img_fuse2 = img_fuse2(1:M1,1:N1,1:3);
img_fuse0 = idwt2(img_fuse2,cH10,cV10,cD10,'sym4');
for k=1:3
%     a=max(max(img_fuse0(:,:,k)));
%     b=min(min(img_fuse0(:,:,k)));
%     img_fuse0(:,:,k)=1/(b-a)*img_fuse0(:,:,k)-(1/(b-a)*a);
 
imshow(img_fuse0);
guidata(hObject, handles);


% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
[FileName, FilePath]=uigetfile('*.tif;*.jpg;*.png;*.img;*.gif;*.bmp;','请选择图像数据');
str=[FilePath FileName];
image_1=imread(str);
handles.imfusion1=image_1;
axes(handles.axes1);
imshow(image_1);
[FileName, FilePath]=uigetfile('*.tif;*.jpg;*.png;*.img;*.gif;*.bmp;','请选择图像数据');
str=[FilePath FileName];
image_1=imread(str);
handles.imfusion2=image_1;
axes(handles.axes2);
imshow(image_1);
guidata(hObject, handles);