|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
; ^* F/ c" I* M9 m! o
load tartan;5 ^; I; [7 e0 z& i ^2 }
X1=X;map1=map;
6 s a) l# z. l o2 r) K, |3 B7 zload sinsin;
7 D8 q& z/ x; v4 _9 V2 oX2=X;map2=map;& I- \4 H+ l$ ?0 v0 a0 x
% 打开图像
% P" q/ m7 Y6 P7 hsubplot(121)7 p# {4 j5 v* _0 q- ~) X% w
image(X1),colormap(map1);& H$ n4 ?5 t8 t' u3 @: H
title('图像tartan')
4 |" ]0 C+ C9 ~& A8 {8 x+ ysubplot(122), m6 W! N( Y2 t L/ u: U
image(X2),colormap(map2);, U; o$ j4 N6 d
title('图像sinsin')
' I% k" S9 ~# h* L( x1 g
% r# l( S* O: [! A7 ~4 j+ V9 y+ j6 x& h. O
%%% 图像直接融合
: ~; H* p9 D2 R% lfigure,subplot(131)4 a: }! I% l4 S% \5 x' J0 Q4 a
image((X1+X2)/2),colormap(map2);- y* ]6 w1 e. T/ I& D
% 在空域内图像直接融合
; v O6 C; g5 T5 f O' }+ {$ ^! ztitle('两图像直接相加融合')
) U1 O" a* `$ w& K8 u8 _
8 r3 J7 h3 ]7 y( B# c6 a8 W# m V. Q
% L: ?& a r% O% Y' ^4 v; F%%% 图像直接融合9 E* \: _* ~6 h4 k, a# q
figure,subplot(131)# w, t4 w4 h# z
image((X1+X2)/2),colormap(map2);
+ q1 O2 h' w' M5 p' n! f4 ^7 O% 在空域内图像直接融合
7 I r1 G. O8 Vtitle('两图像直接相加融合')
& H" F$ ?% D* C& D) k% N8 H( k L- g+ r7 d+ [
" y2 Z9 |* j$ r- `, Y8 l%%% 图像小波变换融合
" U4 h7 X C+ c' r. I( M[C1,L1]=wavedec2(X1,2,'sym4');9 @; c7 {, q) Q) d
[C2,L2]=wavedec2(X2,2,'sym4');
+ q/ `( x. q. R5 E/ ^& T% 分别对上述两图像进行分解# u& V4 N- B( Q) M
C=C1+C2;
; C! j0 E) Z) s% 对分解系数进行融合/ C* T# `! y( a* e6 T" i7 ?' n# F
X=waverec2(C,L1,'sym4');2 U0 Q% w5 N8 ]' C
% 融合后的信号进行图像重构
, Q+ K+ o. _- r; \+ |4 bsubplot(133)
) j* ?7 U$ T- v& ^image(X/2),colormap(map2);4 M" Y0 [- X. |* |3 I
% 显示小波变换的图像融合
/ I+ @0 i0 J/ l/ v! Htitle('两图像小波变换融合')7 z& o6 D2 G( e' ?" A
4 m T* }% d6 J& q( O8 q. @+ a6 `! f4 h8 }$ f
Z& W( i( F5 u) _( _" Y |
|