Question

我有一个形状为(Height, Width)的NumPy数组。现在我想计算这个数组中存在的所有元素的值。该值是元素位置[x,y]的函数。

假设我想将所有元素的值分配为(x**2+y**2)/2。我想在不使用for循环的情况下执行此操作。有没有办法做到这一点？

Answer 1

必须是这样的：

numpy.fromfunction(lambda i, j: (i**2+j**2)/2, (3, 3), dtype=int)

更多信息：http://docs.scipy.org/doc/numpy/reference/generated/numpy.fromfunction.html#numpy.fromfunction

Answer 2

您还可以查看meshgrid，mgrid和/或indices：

>>> H, W = 4,5
>>> x, y = np.indices([H, W])
>>> m
array([[  0. ,   0.5,   2. ,   4.5,   8. ],
       [  0.5,   1. ,   2.5,   5. ,   8.5],
       [  2. ,   2.5,   4. ,   6.5,  10. ],
       [  4.5,   5. ,   6.5,   9. ,  12.5]])

这是有效的，因为x和y是具有相应x和y坐标的数组：

>>> x
array([[0, 0, 0, 0, 0],
       [1, 1, 1, 1, 1],
       [2, 2, 2, 2, 2],
       [3, 3, 3, 3, 3]])
>>> y
array([[0, 1, 2, 3, 4],
       [0, 1, 2, 3, 4],
       [0, 1, 2, 3, 4],
       [0, 1, 2, 3, 4]])

meshgrid和mgrid允许更精细的控制，例如

>>> x, y = np.meshgrid(np.linspace(0, 1, 5), np.linspace(0, 10, 3))
>>> x
array([[ 0.  ,  0.25,  0.5 ,  0.75,  1.  ],
       [ 0.  ,  0.25,  0.5 ,  0.75,  1.  ],
       [ 0.  ,  0.25,  0.5 ,  0.75,  1.  ]])
>>> y
array([[  0.,   0.,   0.,   0.,   0.],
       [  5.,   5.,   5.,   5.,   5.],
       [ 10.,  10.,  10.,  10.,  10.]])

NumPy地图计算取决于索引

2 个答案: