求解 numpy/scipy 中的最小二乘法？

Question

Lucas Kanade 实施

我正在尝试在不使用任何库的情况下为某个项目实现 Lucas Kanade 算法。我正在关注 this guide from opencv.

我一直在翻译下面的最小二乘方程。

我不确定下面的代码是计算方程式的正确方法吗？

我的代码

inv(X)*np.matrix(T)

其中

X =

T =

如果是，请打赏我。如果没有，请提供演示如何使用 numpy 或 scipy.

I am not sure how to validate the correctness of least square. Otherwise I would not have asked this question

Answer 1

我假设您已经在等式右侧创建了两个矩阵（第一个矩阵 X 和第二个矩阵 T）。一旦你有了这两个矩阵，你就可以做这样的事情

import numpy as np

uvVector = np.dot(np.linalg.inv(X), T)

np.dot 做矩阵乘法。