|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
8 @! d! \2 W: v* Iload tartan;. R" A- `+ _6 b& y! a* g
X1=X;map1=map;2 H' q# u1 f+ K6 o2 h& h9 P
load sinsin;
& Y4 H( X5 q3 L4 x! N/ i( DX2=X;map2=map;
0 P9 u- A% E% q4 e- e% 打开图像: n4 b4 D1 Q6 v: Z3 ~
subplot(121). w' M3 |+ O7 A
image(X1),colormap(map1);
4 a3 d) f5 y" a3 O" W+ Q7 M" ltitle('图像tartan')
) o* H, m) }7 isubplot(122)
& n6 }, D% m2 m9 `1 \' G% oimage(X2),colormap(map2);
- z# |* g0 Q7 k! rtitle('图像sinsin'). `1 f3 M7 Y6 x7 v5 a
, ]3 R8 g6 T$ `
+ x) o4 {) L( C1 F$ w+ l8 n: Y; m3 L%%% 图像直接融合
2 P' V- p, E6 m6 w$ Hfigure,subplot(131)
5 G8 B/ i3 v8 B1 E" ?: limage((X1+X2)/2),colormap(map2);1 w7 U8 z8 s: ]- z% g9 u* n
% 在空域内图像直接融合# o6 z( J, b1 i7 y
title('两图像直接相加融合')% j, ?. Y. U: E# @8 ]) M5 Q
6 n- S/ d1 m- q8 l8 U7 h: O$ s2 U! y6 G
%%% 图像直接融合
# O& S7 _7 W* Ffigure,subplot(131)
3 d6 W% Q' Q. ?2 uimage((X1+X2)/2),colormap(map2);
* w+ L2 s* O! Z9 s% 在空域内图像直接融合
+ o F5 m$ d! o* D# k0 U& Ftitle('两图像直接相加融合')5 D; t# K5 [: b+ n6 P. c
) H( k$ ` ^: j0 @7 ~- }" n
( H7 K3 F3 a1 t/ f9 t& I n7 @3 \# |%%% 图像小波变换融合
& j: R% d9 x+ E2 W[C1,L1]=wavedec2(X1,2,'sym4');
. Q1 y: \/ g/ r[C2,L2]=wavedec2(X2,2,'sym4');
& q4 [4 n7 t# Z' B/ k% 分别对上述两图像进行分解/ J9 `) s ^4 p8 t4 K
C=C1+C2;+ [$ d4 g" p$ j& x6 a0 }5 n
% 对分解系数进行融合
$ ]0 T7 J5 Q8 l2 jX=waverec2(C,L1,'sym4');2 t- ^0 i7 E \: e4 v' h
% 融合后的信号进行图像重构
6 U" R9 c0 k; }& e: ]subplot(133)' p4 W3 u: V9 `4 V
image(X/2),colormap(map2);' z1 N+ c, d/ {
% 显示小波变换的图像融合2 ?: \: C* z) y6 i- J7 Q5 b
title('两图像小波变换融合')
2 M1 e2 I* d. J; a5 T! s: B M7 R. D! R: _! `9 Q6 V4 B) C
) m5 y: U5 s; n. F" }% w4 s5 m9 c: t* @/ o/ P, |: C
|
|