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

ulppknot

源代码
3 ]6 X; s5 @: p  o' `% mtic
% function th=Otsu;
0 r: T. B8 ^. V8 n5 [+ H  a=imread('22.bmp');
  c=a;
   %imshow(a)
' s) F. W6 i/ Q! W2 z9 O9 U5 Y   count=imhist(a);%直方图
& ^5 j, s& @) j. d# G6 L# K   [m,n]=size(a);
   N=m*n;
   L=256;
) S' u" r9 G- v9 W5 G8 y   count=count/N;

for i=1: L%这两个循环是将两头不存在的灰度值去掉，提高运行速度
4 X8 M2 L7 l  n/ C& T6 ]5 T    if count(i)~=0
        st=i-1;
        break;
    end
end
% ?* v! p6 e8 A; B2 B% ~; q+ Dfor i=L:-1:1
    if count(i)~=0
        nd=i-1;
; }$ I. ~0 c9 h6 o        break;
9 V0 d5 m, s* ~+ U: J4 w    end
+ c% r4 R. C: C2 |5 Dend
f=count(st+1:nd+1);  %f是每个灰度出现的概率
/ ~6 O; J0 |  `8 o/ m0 np=st;   q=nd-st;
  S" F$ c% n9 K* _2 Ou=0;
- b: g' r. ^6 X: M& ~( Hfor i=1:q
    u=u+f(i)*(p+i-1);  %u是像素的平均值
    ua(i)=u;           %ua（i）是前i个像素的平均灰度值
! k( r5 m  n' }& ~end;

for i=1:q
    w(i)=sum(f(1:i));  %w（i）是前i个像素的累加概率？？
3 q# R/ M% M- j, ?end;

+ [! i9 [( w+ C. wd=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点
th=tp+p;%这一句很重要
8 _5 f; B; P" [
for i=1:m%二值化
    for j=1:n
( k" ~/ v' U. ?        if a(i,j)>th
# s: W4 ^1 W4 m            a(i,j)=0;
        else
            a(i,j)=255;
        end
    end
end
subplot(2,1,1),imshow(c);
subplot(2,1,2),imshow(a);
8 ~# ]9 ^; ~* h- b3 T7 K( Mtoc%tic，toc可以显示程序运行时间

* I% H4 _8 K* l  e5 D2 j
6 a" T0 W  Q! v* ^; b" l
! @. C! V) }  T' Z

SsaaM7

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