二维OTSU（二维大津法），matlab实现，（复制粘贴即可），二维熵

我发现CSDN搜索故意屏蔽调免费、公开的代码，特别推荐的  都是收费，收币的代码或项目。我分享的都是开源的，大家最好关注，收藏我，否则在CSDN 搜索中，即使你属于关键词准确，也并不能找到我。

欢迎 收藏 ，关注，加点赞。 

就喜欢开源，免费。

好了，废话不说直接上代码。

二维Otsu matlab 代码 

% 二维Otsu
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(2,3,1);
imshow(I),title('原图');numel = numel(size(I));
% 图像灰度化
if numel>2I_gray = rgb2gray(I);
elseI_gray = I;
end
subplot(2,3,2);
imshow(I_gray),title('灰度图');% 图像大小
[m,n] = size(I_gray);% NI=filter2(fspecial('average',3),I_gray); % 第一行  第一列  最后一行  最后一列  不满足八邻域  补0计算NI = I_gray;% 3×3模板
for i = 2:m-1for j = 2:n-1pixsum=sum(sum(I_gray(i-1:i+1,j-1:j+1)));pix=round(pixsum/9);NI(i,j)=pix;end
end
subplot(2,3,3);
imshow(NI),title('邻域平均灰度图');Hist=zeros(256,256);
for i = 1:mfor j = 1:npix1=I_gray(i,j);pix2=NI(i,j);Hist(pix1+1,pix2+1)=Hist(pix1+1,pix2+1)+1;           end
end
subplot(2,3,4);
mesh(double(Hist));
title('二维灰度直方图');
xlim([0 256]);% X轴范围
ylim([0 256]);% Y轴范围
xlabel('灰度级');
ylabel('邻域平均灰度级');
zlabel('像素数');% 灰度级
len = length(Hist);
% 灰度概率
p = Hist/(m*n);ux = 0;
uy = 0;
for t = 1:lenfor s = 1:lenux = ux + t * p(s,t);uy = uy + s * p(s,t);end
endg1 = 0;
for t = 1:lenfor s = 1:lenw1 = sum(sum(p(1:t,1:s)));w2 = sum(sum(p(t+1:len,s+1:len)));u1x = 0;u1y = 0;u2x = 0;u2y = 0;for i = 1:tfor j = 1:su1x = u1x + i * p(i,j);u1y = u1y + j * p(i,j);endendfor i = t+1:lenfor j = s+1:lenu2x = u2x + i * p(i,j);u2y = u2y + j * p(i,j);endendg2=w1*((u1x-ux)^2+(u1y-uy)^2)+w2*((u2x-ux)^2+(u2y-uy)^2);if g2>g1T = t-1;S = s-1;g1 = g2;endend
endfor x = 1:mfor y = 1:nif I_gray(x,y)>TI_gray(x,y)=255;elseI_gray(x,y)=0;endend
endsubplot(2,3,5);
imshow(I_gray),title({('二维Ostu阈值分割'),['T=' num2str(T)],['S=' num2str(S)]});
disp(['etime程序总运行时间：',num2str(etime(clock,t1))]);