如何获取矩阵中元素周围的所有邻居值?

时间:2018-03-10 15:02:40

标签: python python-3.x list matrix nearest-neighbor

我需要在python中的矩阵中获取元素周围所有邻居的值。假设我有一个如下所示的矩阵,

  matrix=[[1,2,3,4],
             [5,6,7,8],
             [9,10,11,12]]

对于第一个元素,即matrix[0][0],邻居为[2,5,6]

对于matrix[0][1],邻居为[1,3,5,6,7]

对于matrix[0][2],邻居为[2,4,6,7,8]

对于给定元素,我需要获取这些值列表。

我可以通过比较i = 0时的i,j值,j = 0得到矩阵[0] [1],矩阵[1] [0],矩阵[1] [1]使用开关情况和所以。但它将成为冗长的代码。是否有任何内置功能或任何模块可以使上述任务更简单?

2 个答案:

答案 0 :(得分:3)

如果您不关心效率,请使用import scipy, scipy.ndimage def nb_vals(matrix, indices): matrix = scipy.array(matrix) indices = tuple(scipy.transpose(scipy.atleast_2d(indices))) arr_shape = scipy.shape(matrix) dist = scipy.ones(arr_shape) dist[indices] = 0 dist = scipy.ndimage.distance_transform_cdt(dist, metric='chessboard') nb_indices = scipy.transpose(scipy.nonzero(dist == 1)) return [matrix[tuple(ind)] for ind in nb_indices]

>>> matrix=[[1,2,3,4],
... [5,6,7,8],
... [9,10,11,12]]
>>>
>>> nb_vals(matrix, [1,1])
[1,2,3,5,7,9,10,11]

e.g。

>>> arr_shape = (2,3,4,5)
>>> testMatrix = scipy.array(scipy.random.random(arr_shape)*10, dtype=int)
>>> print(testMatrix)
[[[[7 0 0 1 9]
   [9 5 8 5 8]
   [4 0 0 8 0]
   [1 9 1 3 2]]

  [[9 2 3 3 5]
   [2 3 3 7 9]
   [6 5 6 6 2]
   [9 1 1 0 0]]

  [[8 8 5 7 9]
   [9 0 7 7 6]
   [3 8 7 6 4]
   [8 7 5 5 9]]]

 [[[8 9 2 0 0]
   [8 3 5 5 2]
   [4 0 1 0 3]
   [1 0 9 1 3]]

  [[6 9 2 5 2]
   [2 7 5 5 3]
   [6 7 2 9 5]
   [4 2 7 3 1]]

  [[1 7 7 7 6]
   [5 1 4 1 0]
   [3 9 4 9 7]
   [7 7 6 6 7]]]]

>>> nb_vals(testMatrix, [1,2,2,3])
[3, 7, 9, 6, 6, 2, 1, 0, 0, 7, 7, 6, 7, 6, 4, 5, 5, 9, 5, 5, 3, 2, 9, 5, 7, 3, 1, 4, 1, 0, 4, 7, 6, 6, 7]

这是维度不可知的(适用于输入的任意数量的维度"矩阵")并处理任意数量的索引,您可能希望找到它们周围的邻居。

1

此解决方案在类似图像的二进制数组掩码上使用棋盘式倒角变换,其中0等于掩码上的白色像素,SIGKILL等于黑色像素(背景) )到面具。倒角变换计算所有白色像素与背景的棋盘距离;计算为1的距离的所有像素位置都是邻居,并返回输入数组上这些位置的值。

答案 1 :(得分:1)

您可以在矩阵中构建坐标和相应值的字典,以便进行更简单的查找:

matrix=[[1,2,3,4],
        [5,6,7,8],
        [9,10,11,12]]

def get_neighbors(a, b):
  d = {(i, c):matrix[i][c] for i in range(len(matrix)) for c in range(len(matrix[0]))}
  return filter(None, [d.get(i) for i in
    [(a+1, b+1), (a, b+1), (a+1, b), (a-1, b+1), (a-1, b), (a, b-1), (a+1, b-1)]])

cords = [(0, 0), (0, 1), (0, 2)]
results = [get_neighbors(*i) for i in cords]

输出:

[[6, 2, 5], [7, 3, 6, 1, 5], [8, 4, 7, 2, 6]]