Question

急需高人将以下MATLAB代码的“每一段”的意思、意义诠释！非常感谢！

clear all;
x=imread('football256.bmp');
a=x;
% a=rgb2gray(x);
imshow(a)
count=imhist(a);
[m,n]=size(a);
N=m*n;
L=256;
count=count/N;

for i=1:L
if count(i)~=0
st=i-1;
break;
end
end
for i=L:-1:1
if count(i)~=0
nd=i-1;
break;
end
end
f=count(st+1:nd+1);
p=st; q=nd-st;
u=0;
for i=1:q
u=u+f(i)*(p+i-1);
ua(i)=u;
end;

for i=1:q
w(i)=sum(f(1:i));
end;

d=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);
th=tp+p;

for i=1:m
for j=1:n
if a(i,j)>th
a(i,j)=0;
else
a(i,j)=255;
end
end
end
imshow(x);figure,imshow(a);

Answer 1

楼主是不是一点matlab经验都没有啊？
x=imread('football256.bmp');%是读图像
a=rgb2gray(x);%将真彩色图像转换成灰度图像
imshow(a)%显示图像
count=imhist(a);%做直方图

%%下面一段是找到256个灰度级中，此图最开始和物简最后的非零值的灰度级及中间所含的灰度级个数。
%%f是掐头去尾留下的图像直方图
for i=1:L
if count(i)~=0
st=i-1;
break;
end
end
for i=L:-1:1
if count(i)~=0
nd=i-1;
break;
end
end
f=count(st+1:nd+1);
p=st; q=nd-st;
%% ua是累加灰罩隐裤度级与相应像素个数乘积矩阵， w是累加灰度级矩阵
for i=1:q
u=u+f(i)*(p+i-1);
ua(i)=u;
end;

for i=1:q
w(i)=sum(f(1:i));
end;

%%这个算式不知道什么意思，就是求最优灰度级，存入th
d=(u*w-ua).^2./(w.*(1-w));
[y,tp]=max(d);
th=tp+p;
%%根据th值，将灰度携销图像变成黑白图像
for i=1:m
for j=1:n
if a(i,j)>th
a(i,j)=0;
else
a(i,j)=255;
end
end
end
imshow(x);figure,imshow(a);

追问

MATLAB不是太会，但是上面的你给的比较详细的那些我都知道，而下面开始的程序，也就是你给的不详细的地方我不太明白，真希望你能把下面的讲解的细致一点，我这是一维otsu法，W应该是前i个像素的累加概率，u是像素的平均值，f是每个灰度出现的概率，我想知道稍微详细一点的，哪怕是稍微，好不好，谢谢！不然我没法完成的—。—！

追答

算法我就不看了，我也不是搞图像处理的，下面注释仅供参考
count=imhist(a);%做直方图
[m,n]=size(a);%m是a的横像素数，n是纵像素数
N=m*n;%N是像素总个数
L=256;%256个灰度级
count=count/N;直方图归一，非说成概率也成吧……

%%下面一段是找到256个灰度级中，此图最开始和最后的非零值的灰度级值及中间所含的灰度级个数。
%%f是掐头去尾留下的图像直方图
for i=1:L %循环，灰度级从1到L遍历搜索，找到最早非零的灰度级的序号。
if count(i)~=0 %判断第i个灰度级是否为0
st=i-1; % st为第i-1个灰度级
break; %如果上述条件为真，即在第i个灰度级上的像素个数非零，就跳出循环，由上句可知1-st间的灰度级上的像素个数都为0，st+1开始非全为零
end
end
for i=L:-1:1%由L到1反向搜索，原理同上，L-nd间的灰度级上的像素都为0，nd-1开始往后非全为零
if count(i)~=0
nd=i-1;
break;
end
end
f=count(st+1:nd+1); %%意思是f是count中st+1--nd+1之间的一段，1--st及nd+2-L都去掉，因为这些灰度值上分布的像素个数为0.
p=st; q=nd-st;
%% 下面的你好像知道了，是具体算法
for i=1:q
u=u+f(i)*(p+i-1);
ua(i)=u;
end;

for i=1:q
w(i)=sum(f(1:i));
end;

%%这个算式不知道什么意思，就是求最优灰度级，存入th
d=(u*w-ua).^2./(w.*(1-w));%%在矩阵中对每个元素进行操作，而不是矩阵乘法
[y,tp]=max(d); %%对d的每一列求最大值
th=tp+p;
%%根据th值，将灰度图像变成黑白图像
for i=1:m %%遍历整个图像，
for j=1:n
if a(i,j)>th %如果第i行第j列的像素灰度级大于th
a(i,j)=0;%那么就令它的灰度级变为0，即变成白色
else
a(i,j)=255;%否则变为黑色
end
end
end
imshow(x);figure,imshow(a); %%出原图像和黑白色两幅图

追问

嗯，非常感谢！有很多我不明白的地方都因为你的回答明白了也想起来了一些好久不用的matlab知识！似乎有一种和你共同学习的感觉...呵呵！对了，在第一个for循环中灰度级是有256个那么在matlab中是不是认为灰度是从0到255（同样是256个）？也就是说“0-st”间灰度的像素个数为0？不知道我的认知对不对？后面那个你不知道的算式是求最大门限也就是最佳阈值。如果可以，愿意和你一起讨论！ 共同学习，共同学习！