Procrustes library有一个例子,它演示了如何通过求解Procrustes问题得到两个矩阵的变换矩阵。这个库看起来很旧,在Python 3中不起作用。
我想知道是否有任何方法可以使用SciPy implementation of the Procrustes problem,并能够解决在library's example中讨论的确切问题。
Another StackOverflow question似乎需要我在这里要求的东西,但我不能让它给予我正确的变换矩阵
,该矩阵将源矩阵
变换为接近
的矩阵。
总之,我希望能够使用SciPy库实现this example。
1条答案
按热度按时间eqqqjvef1#
您可以使用
scipy.linalg.orthogonal_procrustes
。下面是一个演示。请注意,函数generateAB
的存在只是为了生成演示中的数组A和B。计算的关键步骤是将A和B居中,然后调用orthogonal_procrustes
。输出量: