fengg.m

一个用matlab编写的 关于笔画提取的仿真试验程序

aa=imread('ej.bmp');
%aa=rgb2gray(A);
imshow(aa);
[m,n]=size(aa);
sigma=zeros([1,256]);
double mean1=0;
double mean2=0;
amap1=aa;
%amap2=edge(amap1,'prewitt');
%figure(3),imshow(amap2);
amap=double(amap1);
int th;
th=0;
for th=0:255
    int thh;
    thh=th+1;
    double grayv1;
    grayv1=0;
    double grayv2;
    grayv2=0;
    double sum1;
    sum1=0;
    double sum2;
    sum2=0;
    for i=1:m
        for j=1:n
            if amap(i,j)<th
                sum1=sum1+1;
                grayv1=grayv1+amap(i,j);
            else 
                sum2=sum2+1;
                grayv2=grayv2+amap(i,j);
            end
        end
    end
    if sum1>0
        mean1=grayv1/sum1;
    else 
        mean1=0;
    end
    if sum2>0
        mean2=grayv2/sum2;
    else mean2=0;
    end
    sigma(1,thh)=sum1*sum2*(mean1-mean2)*(mean1-mean2);
end
double maxsigma;
maxsigma=sigma(1,1);
int thresh;
thresh=0;
for i=0:255
    int thh;
    thh=i+1;
    if maxsigma<sigma(1,thh)
        maxsigma=sigma(1,thh);
        thresh=i;
    end
end

yy=thresh;
for i=1:m
    for j=1:n
       if  amap(i,j) > yy
          amap1(i,j)=255;
      else 
          amap1(i,j)=0;
      end
  end
end

figure(2),imshow(amap1);