如何使用python实现以两图像的像元值为基础通过线性计算得到新的像元值？

Question

就是假定两个不同时期图像的像元值之间存在线性关系，线性回归得出关系式。具体改用什么库做呢，我在numpy中似乎没有看到可以的方法。。。