|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
e/ M$ t6 Q7 E1 ^% K# r
load tartan;4 u; w- y6 U* [/ J) {# f) h
X1=X;map1=map;
6 n7 ~) d2 R5 i& T5 iload sinsin;4 t3 Z' p- q. [6 a6 z
X2=X;map2=map;2 F0 ?0 m" z) P2 D5 ~* D
% 打开图像- ^! v1 n/ s: R
subplot(121)" q$ ~& m% A3 p. T( x; D. A q
image(X1),colormap(map1);' ~# \$ I/ ?+ e3 ? }$ \
title('图像tartan')5 d9 B+ L6 N, J" I
subplot(122)
1 _* i8 M/ a8 G2 I5 Yimage(X2),colormap(map2);
# n# L8 F( l( o5 d; k! d4 D0 A. Wtitle('图像sinsin')6 I/ S( ]% {' \0 O0 y1 `# F
- g, n; n) h0 h# ], A
9 ~$ i1 [- T0 q+ D% W2 E; l/ C! \%%% 图像直接融合9 }- g2 N T. _% ^0 \3 {) c
figure,subplot(131)" E8 n8 w/ a4 W' l3 J
image((X1+X2)/2),colormap(map2);
' ]: c% ?5 c. h/ M7 E" W/ K' W% 在空域内图像直接融合- e+ B9 {4 B7 u4 q$ |+ ?
title('两图像直接相加融合')
/ H- }: u5 y) |. ^# t! p, M; m& J G& U" @) N
; a2 U) z* d1 e
%%% 图像直接融合
/ P% \% ~$ C8 L* ]figure,subplot(131)$ q/ [- t+ K+ N& }% \% m
image((X1+X2)/2),colormap(map2);
( O0 X1 H+ a# U% x% o% 在空域内图像直接融合! v5 L4 _* D7 a v+ j
title('两图像直接相加融合')
! L; f2 P: U, G* h1 R# t' V4 `/ v- h- G2 E
, B4 |1 m( {1 i& i
%%% 图像小波变换融合/ I1 {5 k: t8 U1 Y9 S9 B
[C1,L1]=wavedec2(X1,2,'sym4');! R w: s v+ t4 \; s
[C2,L2]=wavedec2(X2,2,'sym4');: H, ?' ^( `" i$ j
% 分别对上述两图像进行分解
. x. A2 b& L7 G& z9 }6 A* v: xC=C1+C2;$ N1 h" B1 y+ U
% 对分解系数进行融合
2 S* G d0 c0 p2 w' CX=waverec2(C,L1,'sym4');
# F4 K$ k3 `& g" \; G& z. C% 融合后的信号进行图像重构9 ?3 T% L2 o; x. c3 u; [
subplot(133)
' E( E6 P9 \! {8 ^0 B' Oimage(X/2),colormap(map2);: t* V3 K g0 u' _3 J
% 显示小波变换的图像融合# b/ {) q% k5 V$ H
title('两图像小波变换融合'). S# s# E3 i2 Z5 X7 d8 b5 T$ U
; G' x5 V1 l9 ]0 {3 k) t) R1 G. q3 Q& F, c5 {6 B. z) e! ?6 a
% ?5 T! E r1 W. H
|
|