%Otsu法
clc;
clear?all;
I=imread（‘lena.bmp‘）;
%?I=rgb2gray（I）;
I=double（I）;
figureimshow（uint8（I））title（‘原圖‘）;
[mn]=size（I）;
%Th=Otsu（I）;
Th=OtsuFuzzy（I）;
for?i=1:m
????for?j=1:n
????????if?I（ij）>=Th
????????????I（ij）=255;
????????else
???????????I（ij）=0;
????????end
????end
end

%figureimshow（I）title（‘Otsu法‘）;????
figureimshow（I）title（‘OtsuFuzzy法‘）;????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????399??2012-02-15?09:08??Otsu\demo.m

?????文件?????263222??2006-10-25?18:58??Otsu\lena.bmp

?????文件????????800??2012-02-11?10:54??Otsu\Otsu.m

?????目錄??????????0??2012-02-15?09:16??Otsu

-----------?---------??----------?-----??----

???????????????264421????????????????????4