我需要访问 3D 中某个点周围的多个索引。
例如,对于点 (x1
,y1
,z1
) 我需要获取其 3x3x3 邻域的所有索引,以便 ( x1
,y1
,z1
) 居中。对于大小为 3 的邻域,我使用
[x,y,z] = meshgrid(-1:1,-1:1,-1:1);
x_neighbors = bsxfun(@plus,x,x1);
y_neighbors = bsxfun(@plus,y,y1);
z_neighbors = bsxfun(@plus,z,z1);
在这里,我通过添加距 ( x1
,y1
,z1
) 到 3x3x3 框中的任意点。
这给了我 (x1
,y1
,z1
) 3x3x3 邻域的坐标。然后我需要将它们转换为线性索引以便我可以访问它们:
lin_ind = sub2ind(size(volume),y_neighbors,x_neighbors,z_neighbors);
我所做的是昂贵的。
我的问题是,如何避免sub2ind
。如果inx
是(x1
,y1
,z1
)的线性索引,
inx = sub2ind(size(volume),y1,x1,z1);
如何通过对inx
进行加减或任何其他简单操作来找到线性索引的 3x3x3 邻域?
最佳答案
只要知道 3D 数组的维度,就可以计算 3x3x3 邻域中所有元素的线性偏移。为了说明这一点,请考虑一个 4x5 矩阵的二维示例。线性索引如下所示:
1 5 9 13 17
2 6 10 14 18
3 7 11 15 19
4 8 12 16 20
10 的 3x3 邻域是 [5 6 7 9 10 11 13 14 15]
。 15 的 3x3 邻域是 [10 11 12 14 15 16 18 19 20]
。如果我们减去中心元素的索引,在这两种情况下我们都会得到 [-5 -4 -3 -1 0 1 3 4 5]
。更一般地,对于 M
xN
矩阵,我们将有 [-M-1 -M -M+1 -1 0 1 M-1 M M+1 ]
,或 [(-M+[-1 0 1]) -1 0 1 (M+[-1 0 1])]
。
推广到三个维度,如果数组是M
xN
xP
,从中心元素的线性索引偏移量将是[(-M*N+[-M-1 -M -M+1 -1 0 1 M-1 M M+1]) [-M-1 -M -M+1 -1 0 1 M-1 M M+1] (M*N+[-M-1 -M -M+1 -1 0 1 M-1 M M+1])]
。如果您愿意,可以将其重新整形为 3x3x3。
请注意,这种索引不能很好地处理边缘;如果你想在数组的边缘找到一个元素的邻居,你应该先在数组的所有边上填充(从而改变 M
、N
和 P
).
关于matlab - 避免 sub2ind 和 ind2sub,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/16144291/