如何在numpy中旋转索引数组？

Question

我有一个网格，这是一个形状(522, 476)的numpy数组。

然后我有另一个numpy数组，它是形状(2, 3866)的各个xy点（网格的索引）的轨迹。

我使用np.rot90(grid)成功旋转网格。我的问题是如何以相同的方式旋转轨迹，使得各个xy点继续与网格对齐。

Answer 1

如果你总是只想旋转一次90度并始终保持相同的方向（所以rot90没有任何其他参数）你可以使用这个公式：

idx2 = np.array([[n.shape[0]-1-x[1], x[0]] for x in idx])

假设idx是你的索引数组，而n是旋转的数组。它只是使用单个旋转对元素的作用的知识，即将第一维度切换到第二维度，并使第二维度从第一维度的末尾开始计算。

Answer 2

您可以定义旋转功能：

def rotate(origin, point, angle):
    """
    Rotate a point counterclockwise by a given angle around a given origin.

    The angle should be given in radians.
    """
    ox, oy = origin
    px, py = point

    qx = ox + math.cos(angle) * (px - ox) - math.sin(angle) * (py - oy)
    qy = oy + math.sin(angle) * (px - ox) + math.cos(angle) * (py - oy)
    return qx, qy

然后将此函数应用于轨迹的所有（X，Y）点。

origin = tuple(0, 0)
newTrajectory = []
for i in range(0:len(trajectory[0])):
    p = tuple(trajectory[i][0], trajectory[i][1])
    newP = rotate(origin, p, math.pi/2)
    row = [newP[0], newP[1]]
    newTrajectory.append(row)

最好