from
OTSU Binarization
by Shirley
The image "cameraman.tif" will be binarize with the classic OTSU algortithm.
|
| myotsu.m |
close all
clear all
clc
I=imread('D:\Program Files\MATLAB\R2011a\toolbox\images\imdemos\cameraman.tif');
figure(1),imshow(I);
figure(2),imhist(I);
n=imhist(I);
N=sum(n);
max=0;
for i=1:256
P(i)=n(i)/N;
end
for T=2:255
w0=sum(P(1:T));
w1=sum(P(T+1:256));
u0=dot([0:T-1],P(1:T))/w0;
u1=dot([T:255],P(T+1:256))/w1;
sigma=w0*w1*((u1-u0)^2);
if sigma>max
max=sigma;
threshold=T-1;
end
end
bw=im2bw(I,threshold/255);
figure(3),imshow(bw);
|
|
Contact us