%圖象融合（基于小波變化與局部能量）
clear;
clear;
Y=imread（‘clockA-t.jpg‘）;?%裝入原圖像
J?=?imnoise（Y?‘gaussian‘）;
%對圖形的一半加高斯噪聲
%?Y=rgb2gray（Y）;
%??Y=double（Y）/255;
%?figure;
subplot（221）
imshow（Y）;
%?Y=double（Y）/255;
X1=double（Y）/255;
%?a（::1）=Y（1:2561:1281）;
%?a（::2）=Y（1:2561:1282）;
%?a（::3）=Y（1:2561:1283）;
%?a=Y（1:2561:128）;
%?
%?J?=?imnoise（a?‘gaussian‘）;
%?%?Y（1:2561:1281）=J（::1）;
%?%?Y（1:2561:1282）=J（::2）;
%?%?Y（1:2561:1283）=J（::3）;
%?Y（1:2561:128）=J;
%?%?[Xmap]=rgb2ind（Y256）;%將真彩色圖象變成索引色圖象
%?[Xmap]=gray2ind（Y256）;
%?X1=X;map1=map;
%?save?‘picture1‘?X?map;
%?load?picture1;
%?subplot（222）;
%?%?figure;
%?imshow（X1）;
%?%?colormap（map1）;
%?title（‘圖片1‘）;
%?axis?square?%畫出picture1
Y1=imread（‘clockB-t.jpg‘）;?%裝入原圖像
%對圖形的另一半加高斯噪聲
%?Y1=rgb2gray（Y1）;
%?Y1=double（Y1）/255;
X2=double（Y1）/255;
%?b（::1）=Y1（1:256129:2561）;
%?b（::2）=Y1（1:256129:2562）;
%?b（::3）=Y1（1:256129:2563）;
%?b=Y1（1:256129:256）;
%?
%?
%?b?=?imnoise（b?‘gaussian‘）;
%?Y1（1:256129:256）=b;
%?
%?%?Y1（1:256129:2561）=b（::1）;
%?%?Y1（1:256129:2562）=b（::2）;
%?%?Y1（1:256129:2563）=b（::3）;
%?%?[Xmap]=rgb2ind（Y1256）;
%?[Xmap]=gray2ind（Y1256）;
%?X2=X;map2=map;
%?save?‘picture2‘?X?map;
%?load?picture2;
%?%?figure;
%?subplot（223）;
%?imshow（X2）;
%?%?colormap（map2）;
%?title（‘圖片2‘）;
%?axis?square?%畫出wbarb?圖像
subplot（222）;
imshow（X2）;
%對上述二圖像進(jìn)行一級小波分解
%?X1=double（X1）/255;
%?X2=double（X2）/255;
[c1s1]=wavedec2（X11‘sym4‘）;
[c2s2]=wavedec2（X21‘sym4‘）;


%?cA31=appcoef2（c1s1‘sym4‘3）;
%?cH31=detcoef2（‘h‘c1s13）;
%?cV31=detcoef2（‘v‘c1s13）;
%?cD31=detcoef2（‘d‘c1s13）;
%?
%?cH21=detcoef2（‘h‘c1s12）;????
%?cV21=detcoef2（‘v‘c1s12）;
%?cD21=detcoef2（‘d‘c1s12）;
cA11=appcoef2（c1s1‘sym4‘1）;
cH11=detcoef2（‘h‘c1s11）;
cV11=detcoef2（‘v‘c1s11）;
cD11=detcoef2（‘d‘c1s11）;


%?cA32=appcoef2（c2s2‘sym4‘3）;
%?cH32=detcoef2（‘h‘c2s23）;
%?cV32=detcoef2（‘v‘c2s23）;
%?cD32=detcoef2（‘d‘c2s23）;
%?cH22=detcoef2（‘h‘c2s22）;
%?cV22=detcoef2（‘v‘c2s22）;
%?cD22=detcoef2（‘d‘c2s22）;
cA12=appcoef2（c2s2‘sym4‘1）;
cH12=detcoef2（‘h‘c2s21）;
cV12=detcoef2（‘v‘c2s21）;
cD12=detcoef2（‘d‘c2s21）;

R=100;
%將低頻系數(shù)乘1/R
cA11=cA11*（1/R）;
cA12=cA12*（1/R）;
%用新的低頻系數(shù)和原來的高頻系數(shù)一起重構(gòu)新圖象
A1=waverec2（c1s1‘sym4‘）;
B1=waverec2（c2s2‘sym4‘）;
subplot（223）;
imshow（A1）;
subplot（224）;
imshow（B1）;
%用3*3滑動窗口求每個象素的能量
A11=nlfilter（A1[3?3]@myfunction3）;
B11=nlfilter（B1[3?3]@myfunction3）;
[mn]=size（A11）;
for?i=1:m
????for?j=1:n
????????if?abs（A11（ij））>abs（B11（ij））
????????????c（ij）=A1（ij）;
????????else
????????????c（ij）=B1（ij）;
????????end
????end
end
????????
????




%?cA31=cA31.^2;
%?cA32=cA32.^2;
%?cA33=cA31.*cA32;
%?
%?Ea=nlfilter（cA31[3?3]@myfunction1）;
%?Eb=nlfilter（cA32[3?3]@myfunction1）;
%?E=nlfilter（cA33[3?3]@myfunction1）;
%?T=1.5;
%?for?i=1:numel（Ea）
%?M（i）=2*E（i）/（Ea（i）+Eb（i））;
%?Wmin（i）=0.5-0.5*（（1-M（i））/（1-T））;
%?
%?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-09-11?14:07??matlab下小波圖像融合代碼實(shí)現(xiàn)\
?????目錄???????????0??2007-12-20?10:06??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\
?????文件????????7786??2007-04-02?11:27??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\clockA-t.jpg
?????文件????????7271??2007-04-02?11:27??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\clockB-t.jpg
?????文件????????3926??2009-05-27?22:32??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\exercise1.asv
?????文件????????3929??2009-05-28?20:56??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\exercise1.m
?????文件????????2907??2009-05-28?18:50??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\exercise2.asv
?????文件????????2907??2009-05-28?13:48??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\exercise2.m
?????文件?????????180??2007-04-05?11:27??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\myfunction1.m
?????文件?????????140??2007-05-11?15:45??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\myfunction2.m
?????文件?????????125??2007-04-05?11:26??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\myfunction3.m
?????文件?????????282??2007-04-05?11:45??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\myfunction4.m
?????文件?????????134??2011-09-11?14:07??matlab下小波圖像融合代碼實(shí)現(xiàn)\基于matlab的小波圖像融合（多種算法）\說明.txt