3

我有 matlab 数组操作如下：

[M,N]=size(I) ;
J = zeros(2*M,2*N) ;

J(1:2:end,1:2:end) = I ;

J(2:2:end-1,2:2:end-1) = 0.25*I(1:end-1,1:end-1) + 0.25*I(2:end,1:end-1) + 0.25*I(1:end-1,2:end) + 0.25*I(2:end,2:end) ;

J(2:2:end-1,1:2:end) = 0.5*I(1:end-1,:) + 0.5*I(2:end,:) ;
J(1:2:end,2:2:end-1) = 0.5*I(:,1:end-1) + 0.5*I(:,2:end) ;

我正在尝试在 python 中重写相同的操作，我想出了以下内容：

J=numpy.zeros((2*M,2*N))

J[::2,::2] = I ;

J[2:-1:2,2:-1:2] = 0.25*I[1::-1,1::-1] + 0.25*I[2::,1::-1] + 0.25*I[1::-1,2::] + 0.25*I[2::,2::] 

J[2:-1:2,1::2] = 0.5*I[1::-1,] + 0.5*I[2::,]

J[::2,2:-1:2] = 0.5*I[:,1::-1] + 0.5*I[:,2::]

但是python代码给了我不同的结果。

谁能告诉我出了什么问题？

谢谢，

4

2 回答 2

4

逐个浏览这篇文章表明您的范围中有一些错误。我认为您对 python 中的数组有一些误解。

  1. 与 matlab 中数组的第一个元素是 不同array[1]，python 中数组的第一个元素是array[0].
  2. 数组切片语法是array[start:stop:step]，因此要获取从数组中的第五个元素开始到末尾的每隔一个元素array[4::2]。

一块一块地看一遍，你就会发现问题。例如，第二个等式右边的第一个元素应该是：

0.25*I[0:-1, 0:-1]

请注意，您不需要第二个冒号，因为您step是 1，并且如果您想更改步骤，则该步骤最后进行。

于 2013-08-02T04:07:38.337 回答
0

所以这是正确的移植代码：

J[::2,::2] = I ;

J[1:-1:2,2:-1:2] = 0.25*I[0:-1,0:-1] + 0.25*I[1::,0:-1] + 0.25*I[0:-1,1::] + 0.25*I[1::,1::] 

J[1:-1:2,0::2] = 0.5*I[0:-1,] + 0.5*I[1::,]

J[0::2,1:-1:2] = 0.5*I[:,0:-1] + 0.5*I[:,1::]
于 2013-08-07T02:25:31.887 回答