% Otsu大津法 公式1
% 公式1 公式2 由同一公式推导，结果差异的原因是计算的精度
clear all; % 清除所有变量 在function里不能使用
% close all; % 关闭所有figure 在function里不能使用
clc;
[filename,pathname]=uigetfile({'*.jpg;*.bmp;*.tif;*.png;*.gif','All Image Files';'*.*','All Files'});
tic;
t1 = clock;
% 读取图像
I = imread([pathname,filename]);
% I = imread('F:\matlab2016\bin\pic\lena.jpg');
figure;
subplot(3,2,1);
imshow(I),title('原图');
% 图像灰度化
numel = numel(size(I));
if numel>2
I_gray = rgb2gray(I);
else
I_gray = I;
end
subplot(3,2,2);
imshow(I_gray),title('灰度图');
% 转为双精度
I_double = double(I_gray);
% 图像大小
[m, n] = size(I_gray);
% 灰度直方图
hist = imhist(I_gray);
% 灰度级
len = length(hist);
% 灰度概率
P = hist/(m*n);
subplot(3,2,3);
imhist(I_gray),title('灰度直方图');
% u 图像平均灰度
u = 0;
for i = 1:len
u = u + i*P(i);
end
% 用于保存每次计算的方差，与下次计算的方差比较大小
g1 = 0;
for t = 1:len
% 前景点比例
w1 = 0;
for tt = 1:t
w1 = w1 + P(tt);
end
% 背景点比例
w2 = 1 - w1;
% u1 前景平均灰度 u2 背景平均灰度
u1 = 0;
u2 = 0;
for i = 1:len
if i<=t
u1 = u1 + i*P(i);
else
u2 = u2 + i*P(i);
end
end
% 计算前景与背景图像方差
% 公式1
g2 = w1*(u1-u)^2 + w2*(u2-u)^2;
% 保留每次计算的方差
g(t) = g2;
% 每次计算方差后，与上一次计算的方差比较。保留最大方差为 finalT
if g2>g1
finalT = t-1;
g1 = g2;
end
end
for i = 1:m;
for j = 1:n
if I_double(i,j)>finalT
bin(i,j) = 255;
else
bin(i,j) = 0;
end
end
end
subplot(3,2,4);
imshow(bin),title({('Ostu阈值分割'),['T=' num2str(finalT)]});
subplot(3,2,5);
plot(1:len,g),title('方差直方图');
disp(['etime程序总运行时间：',num2str(etime(clock,t1))]);
% graythresh Matlab自带Ostu法
% T = graythresh(I_gray);
% bw = im2bw(I_gray,T);
% subplot(3,2,6);
% imshow(bw),title({('Matlab自带Ostu阈值分割'),['T=' num2str(T*255)]});
