%改進的bernsen算法
clear?all;clc

local=imread（‘t?（145）.jpg‘）;
f=rgb2gray（local）;%灰度轉化
figure（1）;
imshow（f）;
flag=1;%表明為白底黑字
[ab]=size（f）;
A=a*b;??%計算像素總數
N=3;
M=N;

%計算總體閾值T1
pix=0;
num=zeros（1256）;
for?i=0:255
????for?j=1:a
????????for?m=1:b
????????????if?f（jm）==i
????????????????num（i+1）=num（i+1）+1;???%i為0到255統計個灰度值的個數
????????????end
????????end
????end
end


if（sum（num（1:125））>0.4*A）%如果不是白底黑字的話，要進行轉換
????if（sum（num（1:125））>0.9*A）%排除整體偏暗的情況
????????f=f+20;%稍微調亮
????else
????????f=255-f;%取反
????????flag=0;
????end
end

if（flag==0）
????num=fliplr（num）;
end

for?i=1:256
????pix=pix+num（i）;
????if（pix>=（0.1*A））
????????T1=i;
????????break
????end
end

%初始化參數矩陣
T2=zeros（ab）;
T3=T2;
T4=T3;

for?i=M+1:a-M
????for?j=N+1:b-N???????
????????win=f（i-M:i+Mj-N:j+N）;
????????tmax=max（max（win））;
??