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

ulppknot

源代码
tic
% function th=Otsu;
5 c! D' {% R4 O/ w$ e; I  a=imread('22.bmp');
  c=a;
   %imshow(a)
   count=imhist(a);%直方图
) x4 y$ o$ s6 `   [m,n]=size(a);
   N=m*n;
   L=256;
   count=count/N;
! I5 B" U$ X. C/ i9 m
for i=1: L%这两个循环是将两头不存在的灰度值去掉，提高运行速度
    if count(i)~=0
        st=i-1;
        break;
3 j- ?1 P$ B; u1 \/ u    end
end
7 E: w2 r- R6 S: s5 Ifor i=L:-1:1
9 d0 u% ]( F$ q" i* l    if count(i)~=0
; {+ ~6 ]; s* i/ d" v7 t3 `* R4 }4 v9 _- D        nd=i-1;
, p' u5 {# R5 F1 T. |, T7 C        break;
    end
/ `% T9 F" M$ Q6 r6 K; M3 I9 Eend
& }9 Z: x; |1 j* O4 q4 R" Ff=count(st+1:nd+1);  %f是每个灰度出现的概率
p=st;   q=nd-st;
" h6 z: I( Y& b, V- e8 pu=0;
for i=1:q
    u=u+f(i)*(p+i-1);  %u是像素的平均值
    ua(i)=u;           %ua（i）是前i个像素的平均灰度值
9 i+ {0 q) a" P$ Q8 `end;
: Z  n: C7 I- b$ B0 {! S
for i=1:q
    w(i)=sum(f(1:i));  %w（i）是前i个像素的累加概率？？
! T# h" C% ^9 Vend;
# @$ c0 A  ?6 O7 T! t- k3 ]
3 Q+ ~  W+ U7 [2 a4 Ad=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);  %可以取出数组的最大值及取最大值的点
th=tp+p;%这一句很重要

for i=1:m%二值化
    for j=1:n
        if a(i,j)>th
+ B: l- b* k# u) q. E. N+ b            a(i,j)=0;
9 q- R! b6 h1 a7 x        else
            a(i,j)=255;
        end
" z9 f) }$ l% ^. f    end
end
. D( S$ d% F% A9 o7 xsubplot(2,1,1),imshow(c);
subplot(2,1,2),imshow(a);
toc%tic，toc可以显示程序运行时间


! K8 C, w: Z2 J: I, t8 ?

SsaaM7

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