Ostu算法的Matlab源代码以及程序解析

ulppknot

源代码
* ~+ S4 W: `9 ^! c6 i- n1 B5 U9 btic
' n8 N- v. D1 X- W+ e8 L % function th=Otsu;
: t& Z, g6 p" d/ N* i: C  a=imread('22.bmp');
. k: P7 w# M  ~8 X0 m( S/ R3 A  c=a;
   %imshow(a)
   count=imhist(a);%直方图
   [m,n]=size(a);
   N=m*n;
   L=256;
   count=count/N;
  Q2 m/ |, i0 X  f
for i=1: L%这两个循环是将两头不存在的灰度值去掉，提高运行速度
4 |+ k- Z4 ]# D/ s" e5 E6 T5 m0 e, F, `    if count(i)~=0
        st=i-1;
        break;
    end
end
for i=L:-1:1
    if count(i)~=0
4 y- ~# m' z& w. \: ^- g        nd=i-1;
        break;
$ s) O5 g9 e" x4 b/ x    end
end
f=count(st+1:nd+1);  %f是每个灰度出现的概率
p=st;   q=nd-st;
; `' W: e; X3 I4 C9 A, y4 H/ E) Eu=0;
for i=1:q
    u=u+f(i)*(p+i-1);  %u是像素的平均值
    ua(i)=u;           %ua（i）是前i个像素的平均灰度值
( {/ P7 @/ b" z4 O8 }' x' Qend;

' }$ b$ x+ X1 G  ?9 O  a7 yfor i=1:q
/ `- z) ?) ^& r& j. Q    w(i)=sum(f(1:i));  %w（i）是前i个像素的累加概率？？
end;
7 ~) r1 k/ w7 t$ Z" G; c  g
d=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点
7 o/ ~, j9 L: O# C% Tth=tp+p;%这一句很重要
4 M/ X$ \) z2 r7 f/ Z
for i=1:m%二值化
3 S" L) e1 r2 K    for j=1:n
        if a(i,j)>th
6 \* F7 N* i" n7 h1 @: C            a(i,j)=0;
        else
            a(i,j)=255;
0 }) j& [. o! A. S        end
, }, t, |3 @) m! D    end
% @0 a# `# e  O3 mend
8 |# F  b) i4 q' l0 ^3 Xsubplot(2,1,1),imshow(c);
7 s7 y7 D, v) [9 ]% {subplot(2,1,2),imshow(a);
6 T5 b1 b  T9 Y3 V( @5 Y. ztoc%tic，toc可以显示程序运行时间

1 O7 \: @* a4 u) j* c& P

2 E9 i: S* a) X3 u2 q3 m

SsaaM7

Ostu算法的Matlab源代码以及程序解析