Python ValueError：操作数无法与形状一起广播

Question

我正在执行 SVD，当我尝试运行我的代码时，我收到以下错误：

ValueError: 操作数不能与形状一起广播 (375, 375) (375, 500)

我正在使用尺寸为 (500, 375) 的图像

这是我的代码：

from PIL import Image
from Image import new
from numpy import *
import numpy as np
from scipy.linalg import svd

im = Image.open("lake.tif")
pix = im.load()
im.show()
r, g, b = im.split()
R = np.array(r.getdata())
R.shape = (500, 375)
Ur, Sr, VrT = svd(R.T, full_matrices=False)
R1 = Ur * diag(Sr) * VrT

Answer 1

您正在执行组件明智的产品。将这些东西制作成矩阵或使用：

 R1 = np.dot(Ur, np.dot(diag(SR), VrT))

或使用类似的东西：

Ur, Sr, VrT = map(np.asmatrix, (Ur, diag(Sr), Vrt))
R1 = Ur * Sr * VrT

如果你做很多矩阵产品（比如这一行），这会更干净，否则数组通常是更可取的，因为它们是基本类型。如果您愿意，当然也可以单独调用 np.asmatrix。