Matlab 优雅地添加行和列

Question

假设我们有以下随机矩阵：

1 2 3 4
5 6 7 8
9 8 7 6
5 4 3 2

我想将其转换为以下内容：

1 0 2 0 3 0 4 0
0 0 0 0 0 0 0 0
5 0 6 0 7 0 8 0
0 0 0 0 0 0 0 0
9 0 8 0 7 0 6 0
0 0 0 0 0 0 0 0
5 0 4 0 3 0 2 0
0 0 0 0 0 0 0 0

由于某种原因，我无法使用mathjax格式，所以它看起来有点糟糕，对此感到抱歉。要点是，我想在当前行和列之间添加零行和列，以便将其大小增加 2 倍。

我想出了以下代码，但是如果我在由于内存限制问题而无法完成的大图像上使用它，它仅适用于非常小的矩阵。

clear all
I=imread('image.png');
I=rgb2gray(I);
B=zeros(2*size(I));
[x, y]=find(-inf<I<inf);
xy=[x,y];
nxy=xy;
%coord change
nxy=2*xy-1;
B(nxy(:,1),nxy(:,2))=I(xy(:,1),xy(:,2)); 

我预计会很快，因为它是完全矢量化的，具有 maltlab 功能，但它失败得很惨。还有其他优雅的方法可以做到这一点吗？

Answer 1

如果你看一下你的索引向量，这类似于 2x2 矩阵的I([1 1 2 2] ,[1 2 1 2] );，这意味着你索引每一行和每一列两次。正确的解决方案是B(1:2:end,1:2:end)=I;索引每隔一行和每第二列。

Answer 2

这也可以通过一个行来完成，假设你的原始矩阵被称为 A ，然后

kron(A,[1,0;0,0])