leijianfangcha.m

这是计算各种物体类间距离的一种方法

%*******************************最大类间方差法__edit_by_wdb,20050106**************

function [II,th]=leijianfangcha(I);
%I=rgb2gray(I);    
II=I;
count=imhist(I);
L=256;                 %定义灰度级数


[m,n]=size(I);
N=m*n-sum(sum(find(I==0),1));                 %总象素数
count=count/N;       %各象素概率

%找出概率不为0的最小灰度(st)
for i=2:L
    if count(i)~=0;
        st=i-1;
        break;
    end;
end;

%找出概率不为0的最大灰度值(nd)
for i=L:-1:1
    if count(i)~=0;
        nd=i-1;
        break
    end;
end;

f=count(st+1:nd+1);
p=st; q=nd-st;              %q为非0灰度值区间长度

%计算图像平均灰度
u=0;
for i=1:q
    u=u+f(i)*(p+i-1);
    ua(i)=u;
end;

%计算不同门限k时，灰度为1～k的象素点的概率
for i=1:q
    w(i)=sum(f(1:i));
end;

%计算类间方差
d=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);
th=tp+p;

%二值化图像
for i=1:m
    for j=1:n;
        if II(i,j)<=th;
            II(i,j)=0;
        else
            II(i,j)=255;
        end;
    end;
end;