|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
. B) u- u. {/ t+ H/ k& Q" d
load tartan;
% C8 r: E4 S+ D1 M1 s, EX1=X;map1=map;
8 a5 Q% `+ l4 p- @% b/ F1 B3 J& zload sinsin;/ a6 o \& |" d9 j; k
X2=X;map2=map;
# {: i6 z1 e6 K1 i% y% 打开图像& f4 A# X) F4 h* s" j9 m
subplot(121)
0 D. I8 f% ` Himage(X1),colormap(map1);( @! t6 N; u8 B+ a, R# M) Z8 s9 F
title('图像tartan')" O# L6 l& W. |8 u6 k; J; d
subplot(122)
: s R# _/ ]' R6 L8 kimage(X2),colormap(map2);
5 [; o8 w6 e* ^/ R; D% f! p) dtitle('图像sinsin')
' C! F2 {0 Z2 M' J
3 N4 e# |5 `' f
% w' ], h, N* b" \7 k, \8 j9 b%%% 图像直接融合; C9 F9 A: L# K: [! D7 d
figure,subplot(131)+ G' G& |+ k# X/ C9 |5 E5 d4 L
image((X1+X2)/2),colormap(map2);
1 p2 L6 z4 |$ K; ^7 @/ b# ?' I% 在空域内图像直接融合
# F; ]: i& B& |' A0 N0 otitle('两图像直接相加融合')
* m# {9 }" l0 W% \" v; ]% f% R3 D H; O0 A4 C! B
7 {! W' B/ z- C9 b
%%% 图像直接融合
{6 a! p% E$ [- q2 V* y: c( @# afigure,subplot(131): L- ]# Z9 q5 R* ]. d; s! \5 N
image((X1+X2)/2),colormap(map2);
- p/ T2 O& {' r$ y* N% R3 Q) y, w% 在空域内图像直接融合1 t4 T) t4 m8 [
title('两图像直接相加融合')- ]& Z" \: n G5 d& i# w( s
, S: C# j0 h* d; ^
7 b9 k' N' C) I/ K3 n! E
%%% 图像小波变换融合 L& G0 U4 A+ f3 n# X$ {2 d! X
[C1,L1]=wavedec2(X1,2,'sym4');
4 ?1 H3 Z# n# Q* a, x[C2,L2]=wavedec2(X2,2,'sym4');
2 Q) u; Z/ ^2 u8 r* O% 分别对上述两图像进行分解3 P( ?1 }! P$ q6 P1 b7 S
C=C1+C2;6 m" B \* e6 e0 p' M6 i
% 对分解系数进行融合) F% O o, U% X
X=waverec2(C,L1,'sym4');7 t8 M6 p( @# |6 R9 t
% 融合后的信号进行图像重构
2 r" e$ X2 }) o3 T5 Vsubplot(133)
9 K+ A% @3 g# wimage(X/2),colormap(map2);1 B( F8 z0 q7 n6 i, i" D5 U4 ^2 B
% 显示小波变换的图像融合
# }6 Y8 p0 @& m& K5 Ptitle('两图像小波变换融合')+ `! r, l' I- \; D% }
: c( t6 L9 s$ }" l0 v
* E; L. \$ { o' X( B5 _( y% L6 v; d: s" ]0 ?
|
|