hh.rar_最大熵_最大熵准则

%############################################# %2008.1.14 %yuan bao feng %基于最大类间最大距离比准则的图像分割 %图像输入为moon.tif %############################################# %clc; %clear all; %I=imread('moon.tif'); %读入图像 %len=graynum(I); function m=yuanbao1(I) h=imhist(I); %画出灰度直方图 h1=h; len=length(h); %求出所有的可能灰度 [v,r]=size(I); %求出图像的大小 h1=(h1+eps)/(v*r); %算出各灰度点出现的概率 %-------------具体算法------------------------ k=zeros(1,len); for i=1:len if h(i)~=0 k(i)=h(i)/(v*r); else continue; end end for i=1:(len-1) P1=sum(h(1:i))/(v*r); %算出第一类的百分比 P2=sum(h((i+1):len))/(v*r); %算出第二类的百分比 H1=sum((1:i).*k(1:i)); H2=sum(((i+1):len).*k((i+1):len)); m1=(1/P1)*H1; m2=H2*(1/P2); m=P1*m1+P2*m2; S1(i)=P1*((m1-m).^2)+P2*((m2-m).^2); S2(i)=sum((([1:i]-m1*ones(1,length(i))).^2).*k(1:i))+sum((([i+1:len]-m2*ones(1,length(len-i))).^2).*k(i+1:len)); S(i)=S1(i)/S2(i); end m1=max(S); n1=find(S==m1); %-------------结束---------------------------- I1=I(:); for i=1:v*r if I1(i)>=n1 I1(i)=1; else I1(i)=0; end end I2=reshape(I1,v,r); %将向量转化为图像格式 I2=double(I2); I2=mat2gray(I2); imshow(I2);