close all,clear;
I1=imread('medA.gif');
I2=imread('medB.gif');
I1=I1(:,:,1);%%%%%%%%%
I2=I2(:,:,1);
figure(1),imshow(I1)
figure(2),imshow(I2)
[c1,l1]=wavedec2(I1,2,'sym4');
[c2,l2]=wavedec2(I2,2,'sym4');
%c=max(c1,c2);
c=(c1+c2)./2;
%c=(c1+c2);
%X=mean(c);
%计算一副图像的熵
%随机生成图像
%A=floor(rand(8,8).*255);
A=I1;
[M,N]=size(A);
temp=zeros(1,256);
%对图像的灰度值在[0,255]上做统计
for m=1:M;
for n=1:N;
if A(m,n)==0;
i=1;
else
i=A(m,n);
end
temp(i)=temp(i)+1;
end
end
temp=temp./(M*N);
%由熵的定义做计算
result=0;
for i=1:length(temp)
if temp(i)==0;
result=result;
else
result=result-temp(i)*log2(temp(i));
end
end
result
xx=waverec2(c,l1,'sym4');
xx=xx./260;%%%%%%%%%%%%
figure(3),imshow(xx);
评论0