Question

matlab A=I((X(k)-1):(X(k)+1),(Y(k)-1):(Y(k)+1));

Answer 1

语句的功能是对I矩阵进行截取，X(k)-1):(X(k)+1)这三行，和(Y(k)-1):(Y(k)+1)这三列，生成3X3矩阵

追问

　　帮我看看整段代码吧！
　　%I:待分割的图像，rmin ,rmax:瞳孔半径的最大最小值 function [ci,cp,out]=thresh(I,rmin,rmax);
　　I=im2double(I); rows=size(I,1); cols=size(I,2); [X,Y]=find(Irmin)&(Y(k)>rmin)&(X(k)<=(rows-rmin))&(Y(k)<(cols-rmin))
　　A=I((X(k)-1):(X(k)+1),(Y(k)-1):(Y(k)+1));　M=min(min(A));　if I(X(k),Y(k))~=MX(k)=NaN;Y(k)=NaN;　　end

追答

你给的程序有错误，我稍微改了下，感觉大题的意思是寻找亮度小于0.5的区域中的亮度最低点的坐标，你ci cp没用到，不知道做什么用的。：
[ci,cp,out]=thresh(I,rmin,rmax);
I=im2double(I);
rows=size(I,1);
cols=size(I,2);
[X,Y]=find(Irmin)&(Y(k)>rmin)&(X(k)<=(rows-rmin))&(Y(k)<(cols-rmin))%如果点在瞳孔范围之内
A=I((X(k)-1):(X(k)+1),(Y(k)-1):(Y(k)+1));
M=min(min(A));%M取最小值
end
if I(X(k),Y(k))~=M%如果X(k),Y(k)点不是最小值
X(k)=NaN;%赋空
Y(k)=NaN;　　
end
end
建议贴程序认真点。格式自己调好了，再贴上来。不然浪费大家时间。

追问

不好意思！格式问题是字数导致的，上面的问题我有一个最难懂的地方就是，怎样确定点在瞳孔范围内，其中的坐标系的原点、还有X(k)>rmin的含义？ 我要用c怎么实现？

追答

1.(X(k)>rmin)&(Y(k)>rmin)&(X(k)<=(rows-rmin))&(Y(k)<(cols-rmin))这句是判断点是否在瞳孔范围之内的。按我的理解，例如，总大小是10x10的时候，而中间4X4是瞳孔范围的时候，那rows=cols=10，而rmin=3，也就是x在3-7范围内认为是瞳孔范围内的
2.一般程序认为图像坐标系的原点是在左上角；
3.见1
4.matlab转c需要较多的for，架构是基本一样的，比如
for k=1：s
end
变成
for（k=1;k<=s;k++)
{
}