22

所以我发现了这个：

转换 MATLAB 代码时，可能需要先将矩阵整形为线性序列，执行一些索引操作，然后再整形。由于 reshape （通常）在同一存储上生成视图，因此应该可以相当有效地执行此操作。

请注意，Numpy 中 reshape 使用的扫描顺序默认为“C”顺序，而 MATLAB 使用 Fortran 顺序。如果您只是简单地转换为线性序列并返回，这无关紧要。但是，如果您要从依赖于扫描顺序的 MATLAB 代码转换整形，那么这个 MATLAB 代码：

z = reshape(x,3,4);

应该成为

z = x.reshape(3,4,order='F').copy()

在 Numpy 中。

mafs当我在 MATLAB 中执行以下操作时，我有一个名为 的多维 16*2 数组：

mafs2 = reshape(mafs,[4,4,2])

我得到的东西与我在 python 中所做的不同：

mafs2 = reshape(mafs,(4,4,2))

甚至

mafs2 = mafs.reshape((4,4,2),order='F').copy()

对此有什么帮助吗？谢谢你们。

4

2 回答 2

32

例子：

MATLAB：

>> mafs = [(1:16)' (17:32)']
mafs =
     1    17
     2    18
     3    19
     4    20
     5    21
     6    22
     7    23
     8    24
     9    25
    10    26
    11    27
    12    28
    13    29
    14    30
    15    31
    16    32

>> reshape(mafs,[4 4 2])
ans(:,:,1) =
     1     5     9    13
     2     6    10    14
     3     7    11    15
     4     8    12    16
ans(:,:,2) =
    17    21    25    29
    18    22    26    30
    19    23    27    31
    20    24    28    32

Python：

>>> import numpy as np
>>> mafs = np.c_[np.arange(1,17), np.arange(17,33)]
>>> mafs.shape
(16, 2)
>>> mafs[:,0]
array([ 1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16])
>>> mafs[:,1]
array([17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32])

>>> r = np.reshape(mafs, (4,4,2), order="F")
>>> r.shape
(4, 4, 2)
>>> r[:,:,0]
array([[ 1,  5,  9, 13],
       [ 2,  6, 10, 14],
       [ 3,  7, 11, 15],
       [ 4,  8, 12, 16]])
>>> r[:,:,1]
array([[17, 21, 25, 29],
       [18, 22, 26, 30],
       [19, 23, 27, 31],
       [20, 24, 28, 32]])
于 2012-08-09T22:54:32.350 回答
0

我自己也有类似的问题，因为我也在尝试从 MATLAB 过渡到 Python。我终于能够将一个以深度、行、列、格式给出的 numpy 矩阵转换为单张列向量（每张图像）。

在 MATLAB 中，我会做类似的事情：

output = reshape(imStack,[row*col,depth])

在 Python 中，这似乎转化为：

import numpy as np
output=np.transpose(imStack)
output=output.reshape((row*col, depth), order='F')
于 2020-03-10T14:49:12.800 回答