【图像融合】主成分结合小波离散变换PCA-DWT三张图像融合【含Matlab源码 2199期】

?博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
??座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

?一、小波变换彩色图像融合简介

1 基于小波的图像融合
1.1 小波的分解和重构
小波变换是一种能够用来检测信号局部特征的数学工具。当然也可以将二维信号分解成不同分辨率的子带信号。由于图像为二维信号, 可以作以下小波分解:

其中, f (x, y) 为源图像, C0, H, G为一维小波滤波器, h, v, d分别代表水平、垂直和对角分量, H′, G′表示H, G的转置矩阵。

小波变换属于可逆变换, 可以通过以下方式实现信号的重构, 从而实现信号的恢复:

1.2 小波图像融合算法
该方法的基本思想是将待融合图像分解到不同频率上, 然后在不同频率上根据一定的规则完成信息的融合, 最后通过小波逆变换恢复图像。

基本的融合步骤描述如下:

(1) 配准待融合的源图像, 保证相关像素对齐。

(2) 将待融合图像各自完成小波分解, 如果进行层分解, 将会得到个高频子带和1个低频子带。

(3) 据不同融合规则完成高频和低频的系数融合, 低频采用平均法融合, 高频采用绝对值取大完成融合。

(4) 对步骤 (3) 得到的融合图像进行小波逆变换, 得到最终融合结果图像。

2 图像块清晰度的度量
2.1 主成分分析
主成分分析属于一种线性分析方法之一, 其被成功应用于人脸图像识别领域。它的主题思想是寻找最小均方意义下的逼近原始数据的最佳投影[。而对于数据的变异性和主成分分析中特征值的关系给出了证明。

数据矩阵C (标准化矩阵) 分解成正交矩阵U和V, V满足VTV=VVT=I:

由 (5) 式知, 描述变异性的分量ν1也可由特征值λ解释。

根据多聚焦图像的特征, 模糊区域类似于对该区域清晰图像实施了低通平滑滤波, 消除了图像中的部分高频信息, 所以模糊区域的高频局部变异较小, 而清晰区域高频部分变异较大。因此, 多聚焦图像相应像素块内高频信息变异大小可以作为清晰与否的判断依据, 即特征值越大, 图像越清晰。

2.2 图像块按清晰度分类
本文结合分块的思想依据, 将源图像进行区域的划分, 对每一区域块使用PCA分析法, 把特征值作为每一区域块的清晰与否的判定依据。

假设H1和H2为待融合的多聚焦图像P1和P2的高通滤波图像。邻域内像素一般情况下具有较强关联性, 因此, 为提高辨识的准确性, 进一步对H1和H2进行低通滤波, 得到A1和A2。假定Α1k, Α2k代表A1, A2的第k个子块图像, 根据PCA原理计算出AT1kA1k, AT2kA2k的第一主分量对应的特征值λk1, λk2, 则由下式得到描述图像清晰属性的标识矩阵:

其中, T为阀值, Lk1, Lk2和是描述源图像P1和P2第k块清晰属性的标识, 如果为1表示此块是清晰, 为0表示此块模糊, 为0.5表示此块是属于两者间的边界区域。

3 图像融合方法
小波方法需要将整幅图像进行小波变换。考虑到多聚焦图像只有部分模糊的情况下, 这种变换就存在一定的冗余操作, 对于对应区域图像块都是清晰的情况, 此种变换就显得多余了, 而且会增加源图像和重构图像之间的偏差。因此考虑图像融合分成如下两部分来完成:

对于确定的清晰部分直接作为结果图像的对应部分, 即:

其中, Rk表示第k块结果图像块;、A1k, A2k为源图像P1和P2的第k块;Lk1、Lk2对应源图像第k块的标识元素。

对于介于清晰与模糊之间的区域采用小波法完成融合:

其中, F (A1k, A2k) 表示对子块A1k、A2k进行小波融合。

?二、部分源代码

clc; clear all;close all;
% Give number of input images
n=input(‘Number of images for fusion=’);
for i=1:1:n;
[imagefile1 , pathname]= uigetfile(‘.jpg;.bmp;’,‘Open file Eye image’);
name=[pathname,imagefile1];
if imagefile1 ~= 0
a{i} = double(imread([name]));
[row,col,dim]=size(a{i});
if (dim>2)
a{i} = double(rgb2gray(imread([name])));
va{i}=cov(a{i},a{i});
end
end;
[ca{i},ch{i},cv{i},cd{i}] = dwt2(a{i},‘db3’);
vadwt{i}=cov(ca{i},ca{i});
vhdwt{i}=cov(ch{i},ch{i});
vvdwt{i}=cov(cv{i},cv{i});
vddwt{i}=cov(cd{i},cd{i});
avdwt{i}=(vadwt{i}+vhdwt{i}+vvdwt{i}+vddwt{i})/n;
end
% Evaluation of principal components for the relevant coefficients
[fca m1]=fuse_pcaany(ca,n);
[fch m2]=fuse_pcaany(ch,n);
[fcv m3]=fuse_pcaany(cv,n);
[fcd m4]=fuse_pcaany(cd,n);
% principal component averaging
for i=1:1:n;
pc(i)=(m1(i,1)+m2(i,1)+m3(i,1)+m4(i,1))/4;
end
[z3 s3]=size(a{1});
out=zeros(z3,s3);
% Fusion based on average of principal components
for i=1:1:n
dwtpcav = pc(i)*a{i};
out=out+dwtpcav;
end
figure,
for i=1:1:n
subplot(1,n+1,i)
imshow(a{i},[]);
title([‘原图’,num2str(i)])
end

?三、运行结果

?四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]向昌成,杨兴江,李清茂.PCA和小波在多聚焦图像融合中的应用研究[J].阿坝师范高等专科学校学报. 2011,28(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合