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

ulppknot

源代码
tic
% function th=Otsu;
4 g( i  n9 K; J7 F4 ~  a=imread('22.bmp');
2 Z- ?& P! A& }" y  c=a;
# x- k4 n# A2 b. z   %imshow(a)
) a. t- F$ ~; a; O, G* W   count=imhist(a);%直方图
   [m,n]=size(a);
   N=m*n;
   L=256;
   count=count/N;
3 t. h9 D+ T+ r. V( @+ _
for i=1: L%这两个循环是将两头不存在的灰度值去掉，提高运行速度
+ e& X+ f" h7 A$ D% l. p) t    if count(i)~=0
( y  o% q4 ?8 ~8 D% B, E        st=i-1;
0 E8 }  G2 {* k        break;
    end
, B' a3 z  o* ~7 fend
for i=L:-1:1
    if count(i)~=0
        nd=i-1;
; x& l- n# E' T' @( f! ?        break;
    end
- \9 m7 {: @- Nend
: P* z. D+ [1 c! E! N; ^7 |; K3 m) _f=count(st+1:nd+1);  %f是每个灰度出现的概率
p=st;   q=nd-st;
u=0;
+ `4 q1 P. r7 _( ufor i=1:q
    u=u+f(i)*(p+i-1);  %u是像素的平均值
( A" Q- M( B$ P- d1 M) ]% }% L8 Y    ua(i)=u;           %ua（i）是前i个像素的平均灰度值
+ w& A: m3 x8 T" w* Nend;

for i=1:q
    w(i)=sum(f(1:i));  %w（i）是前i个像素的累加概率？？
' u/ H) f& m9 Q9 q" R, vend;

1 b; ~  d* R* g4 C# Q. H3 _d=(u*w-ua).^2./(w.*(1-w));
$ V/ ~- N' @. ^2 q[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点
$ Y) T# T( a# b+ r' P$ m$ Qth=tp+p;%这一句很重要

  H' S# \; a6 K$ v8 u* I1 Gfor i=1:m%二值化
3 f2 g, Z6 ]1 l    for j=1:n
, o6 G4 }6 k3 d  a. v, |        if a(i,j)>th
            a(i,j)=0;
0 V& K9 N' i: R        else
            a(i,j)=255;
        end
    end
8 B4 T/ Z8 l& Bend
subplot(2,1,1),imshow(c);
subplot(2,1,2),imshow(a);
toc%tic，toc可以显示程序运行时间
6 c* q' j1 X: e* }/ Z

SsaaM7

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