我在使用 MATLAB 中的 find 函数时遇到一些问题。我有一个由零和一组成的矩阵(表示结构元素的几何形状),当矩阵元素 = 1 时存在 Material ,当矩阵元素 = 0 时不存在 Material 。矩阵可能具有一般形式如下所示(它会随着几何形状的变化而更新,但这并不重要)。
Geometry = [0 0 0 0 0 0 0 0 0 0;
0 0 1 0 1 0 1 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 1 0 0;
0 0 0 0 0 0 0 0 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 1 1 1 0 1 0 0;
0 0 0 0 0 0 0 0 0 0;]
我试图找到不连续连接的行和列(即行和列在行或列的外部范围之间不都等于 1 的位置),然后更新它们,使它们全部连接的。 IE。上面的矩阵变为:
Geometry = [0 0 0 0 0 0 0 0 0 0;
0 0 1 1 1 1 1 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 0 0 0 0 1 0 0;
0 0 1 1 1 1 1 1 0 0;
0 0 0 0 0 0 0 0 0 0;]
我遇到的问题是我希望能够找到每行(和列)中等于 1 的第一个和最后一个元素的索引,然后将用于更新几何矩阵。
理想情况下,我想用向量表示这些元素,因此遍历各列,找到第一个元素的行号等于 1 并将其存储在名为 rowfirst 的向量中。
即:
rowfirst = zeros(1,numcols)
for i = 1:numcols % Going across the columns
rowfirst(i) = find(Geometry(i,1) == 1, 1,'first')
% Store values in vector called rowfirst
end
并对列重复此操作并找到每行中的最后一个元素。
出于某种原因,我无法将值正确存储在向量中,有人知道我哪里出错了吗?
提前致谢。如果不清楚,请告诉我,因为我可能没有很好地解释问题。
0) bwmorph(Geometry,'close')
在一行中完成所有操作。如果孔可能更大,请尝试 bwmorph(Geometry,'close',Inf)
。
关于你的尝试:
1) 应该是 Geometry(i,:)
而不是 Geometry(i,1)
。
2) 你真正的问题是空矩阵。实际上,如果第 i 行中没有 1,您希望 rowfirst(i)
是什么?
我是一名优秀的程序员,十分优秀!