如何在 Matlab 中轻松创建动态二维向量矩阵

Question

给定我的代码：

G=zeros(height,width); %zeros or whatever

for y = 1 : height
    for x = 1 : width
         magnitude = sqrt(Gx(y,x)^2 + Gy(y,x)^2);
         gradient_direction = atan(Gy(y,x)/Gx(y,x));
         G(y,x) = [magnitude gradient_direction];
    end
end

我不断得到这个（如果我不使用 zeros()）：

下标分配维度不匹配。

或者这个：

Assignment 的非单件 rhs 维度比非单件多 下标

Answer 1

虽然@atru 回答有效，但如果有帮助的话，我想建议一种更快、更整洁的矢量化方式。这里的操作可以很容易地转化为向量化操作：

G=cat(3,hypot(Gy,Gx),atan(Gy./Gx));

Answer 2

通过使用G(y,x,:) = [magnitude, gradient_direction];，您试图将两个值分配给一个为具有索引(y,x) 的单个值保留的位置。解决此问题的一种方法是改用 3 维数组 G，

G=zeros(height,width,2);

for y = 1 : height
    for x = 1 : width
         magnitude = sqrt(Gx(y,x)^2 + Gy(y,x)^2);
         gradient_direction = atan(Gy(y,x)/Gx(y,x));
         G(y,x,:) = [magnitude, gradient_direction];
    end
end

现在在每个点 G(y,x)，您都可以存储这两个值并访问它们，例如 G(1,2,1) 用于 magnitude 在(1,2) 位置和G(1,2,2) 用于gradient_direction。这假设 Gx 和 Gy 都是大小为高 x 宽的数组。

---

需要注意的重要一点是，沿第三维的 G 切片也将是 3D 数组，即 mag_dir = G(3,2,:) 的大小为 [1 1 2] 而不是 [1 2]。这可能会在某些应用程序中导致错误，示例包括尝试将mag_dir 与另一个向量（没有额外维度）和线性代数运算连接。

要解决此问题，请使用 reshape 将维度显式更改为目标维度。对于这里的向量，它将是reshape(mag_dir, 1, 2)。同样适用于像more_md = G(1,:,:) 这样的二维切片——这将需要例如more_md = reshape(more_md,2,5)。