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

ulppknot

源代码
% K, M& T5 ~. w" \( j. ^tic
% function th=Otsu;
! D, Z* E/ t# N1 M- p  a=imread('22.bmp');
$ c* w" o0 G- A8 ]3 ?4 [  c=a;
$ _2 y- t, ?0 U   %imshow(a)
   count=imhist(a);%直方图
   [m,n]=size(a);
   N=m*n;
4 y$ b$ {% n+ k9 ~/ J: H   L=256;
   count=count/N;

for i=1: L%这两个循环是将两头不存在的灰度值去掉，提高运行速度
    if count(i)~=0
        st=i-1;
        break;
; M+ I% x/ y0 @    end
' X3 d9 c; P1 ]5 v: uend
for i=L:-1:1
2 a0 d9 i8 m+ }- s. D! U: S: l" r* l    if count(i)~=0
" T9 ^9 F& F. a$ m  D        nd=i-1;
- i2 b! S: g" X/ \1 J) D        break;
    end
0 n( T( w9 ~. A2 {+ Tend
f=count(st+1:nd+1);  %f是每个灰度出现的概率
p=st;   q=nd-st;
$ C* x# r! }! |( [u=0;
for i=1:q
    u=u+f(i)*(p+i-1);  %u是像素的平均值
  ^, a+ J& N/ Z8 u! J9 g    ua(i)=u;           %ua（i）是前i个像素的平均灰度值
end;

1 q" F- @- f" [, \" ^9 a# S# `" ffor i=1:q
4 W# O5 J  @$ q+ n1 V    w(i)=sum(f(1:i));  %w（i）是前i个像素的累加概率？？
1 T5 r  |' q( g8 e; _( Hend;
/ x" S& \9 Y+ k3 T
d=(u*w-ua).^2./(w.*(1-w));
  J0 H  y. _6 `[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点
th=tp+p;%这一句很重要

for i=1:m%二值化
    for j=1:n
0 W4 k- `; `1 S0 l! a        if a(i,j)>th
) r3 [6 F% V/ \* Q            a(i,j)=0;
        else
3 m& S* M0 L# k6 w, c2 h            a(i,j)=255;
2 }$ j$ ~5 ]6 q  K        end
    end
end
subplot(2,1,1),imshow(c);
; u8 j, `; r% X" Vsubplot(2,1,2),imshow(a);
toc%tic，toc可以显示程序运行时间
8 C$ v* y$ d" g6 i$ ~" T

SsaaM7

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