我最近注意到NumPy ndarays的Python打印功能不一致。例如,它水平打印水平的一维数组:
import numpy as np
A1=np.array([1,2,3])
print(A1)
#--> [1 2 3]
但是一维水平阵列,其垂直有多余的括号:
A2=np.array([[1],[2],[3]])
print(A2)
#--> [[1]
# [2]
# [3]]
水平的一维垂直阵列:
A3=np.array([[1,2,3]])
print(A3)
#--> [[1 2 3]]
和2D数组:
B=np.array([[11,12,13],[21,22,23],[31,32,32]])
print(B)
# --> [[11 12 13]
# [21 22 23]
# [31 32 32]]
第一个尺寸现在是垂直的。对于更高尺寸的尺寸,情况更糟,因为它们都是垂直打印的:
C=np.array([[[111,112],[121,122]],[[211,212],[221,222]]])
print(C)
#--> [[[111 112]
# [121 122]]
#
# [[211 212]
# [221 222]]]
我认为,一致的行为是水平打印偶数尺寸,垂直打印奇数尺寸。使用Unicode字符,可以很好地格式化它。我想知道是否可以创建一个函数来打印上面的数组:
A1 --> [1 2 3]
A2 --> ┌┌─┐┌─┐┌─┐┐
│ 1 2 3 │
└└─┘└─┘└─┘┘
A3 --> ┌┌─┐┐ # \u250c\u2500\u2510
│ 1 │ # \u2502
│ 2 │
│ 3 │
└└─┘┘ # \u2514\u2500\u2518
B --> ┌┌──┐┌──┐┌──┐┐
│ 11 21 31 │
│ 12 22 32 │
│ 13 23 33 │
└└──┘└──┘└──┘┘
C --> ┌┌─────────┐┌─────────┐┐
│ [111 112] [211 212] │
│ [121 122] [221 222] │
└└─────────┘└─────────┘┘
我发现这个gist可以处理不同数量的数字。我试图为递归函数创建原型以实现上述概念:
def npprint(A):
assert isinstance(A, np.ndarray), "input of npprint must be array like"
if A.ndim==1 :
print(A)
else:
for i in range(A.shape[1]):
npprint(A[:,i])
这有点适用于A1
,A2
,A3
和B
,但不适用于C
。如果您能帮助我知道npprint
应该如何实现上述任意维度numpy ndarrays的输出,将不胜感激?
P.S.1。。在Jupyter环境中,可以在Markdown中使用LaTeX \mathtools
\underbracket
和\overbracket
。 Sympy的漂亮打印功能也是一个很好的起点。它可以使用ASCII,Unicode,LaTeX ...
P.S.2。。有人告诉我,ndarray的打印方式确实存在一致性。但是恕我直言,这是一种有线的和非直觉的。具有灵活漂亮的打印功能可以帮助以各种形式显示ndarray。
P.S.3。。Sympy家伙已经考虑了我在这里提到的两点。它们的Matrix模块非常一致(A1
和A2
是一致的),并且它们还具有pprint
函数,该函数执行相同的操作,我希望在这里npprint。
P.S.4。。对于那些遵循这个想法的人,我在this Jupyter Notebook
中集成了所有内容。答案 0 :(得分:7)
这对我来说是一个很大的启示,了解numpy数组与我想到的MATLAB矩阵或多维数学数组完全不同。它们是同质且统一的嵌套Python列表。我还了解到,numpy数组的第一个维度是水平括号中最深/最里面的一对方括号,水平打印,然后从那里垂直打印第二个维度,第三个维度以空格隔开...
无论如何,我想拥有一个ppring
函数(受Sympy的命名约定启发)会很有帮助。所以我要在这里放一个非常糟糕的实现,希望它能启发其他高级Python使用者提出更好的解决方案:
def pprint(A):
if A.ndim==1:
print(A)
else:
w = max([len(str(s)) for s in A])
print(u'\u250c'+u'\u2500'*w+u'\u2510')
for AA in A:
print(' ', end='')
print('[', end='')
for i,AAA in enumerate(AA[:-1]):
w1=max([len(str(s)) for s in A[:,i]])
print(str(AAA)+' '*(w1-len(str(AAA))+1),end='')
w1=max([len(str(s)) for s in A[:,-1]])
print(str(AA[-1])+' '*(w1-len(str(AA[-1]))),end='')
print(']')
print(u'\u2514'+u'\u2500'*w+u'\u2518')
结果对于1D和2D数组是可以接受的:
B1=np.array([[111,122,133],[21,22,23],[31,32,33]])
pprint(B1)
#┌─────────────┐
# [111 122 133]
# [21 22 23 ]
# [31 32 33 ]
#└─────────────┘
这确实是一个非常糟糕的代码,它仅适用于整数。希望其他人会提出更好的解决方案。
P.S.1。 Eric Wieser已经为IPython / Jupiter实现了一个非常漂亮的HTML原型,可以看到here:
您可以关注有关numpy邮件列表here的讨论。
P.S.2。我也发布了这个想法here on Reddit。
P.S.3 我花了一些时间将代码扩展到3D维度数组:
def ndtotext(A, w=None, h=None):
if A.ndim==1:
if w == None :
return str(A)
else:
s= '['
for i,AA in enumerate(A[:-1]):
s += str(AA)+' '*(max(w[i],len(str(AA)))-len(str(AA))+1)
s += str(A[-1])+' '*(max(w[-1],len(str(A[-1])))-len(str(A[-1]))) +'] '
elif A.ndim==2:
w1 = [max([len(str(s)) for s in A[:,i]]) for i in range(A.shape[1])]
w0 = sum(w1)+len(w1)+1
s= u'\u250c'+u'\u2500'*w0+u'\u2510' +'\n'
for AA in A:
s += ' ' + ndtotext(AA, w=w1) +'\n'
s += u'\u2514'+u'\u2500'*w0+u'\u2518'
elif A.ndim==3:
h=A.shape[1]
s1=u'\u250c' +'\n' + (u'\u2502'+'\n')*h + u'\u2514'+'\n'
s2=u'\u2510' +'\n' + (u'\u2502'+'\n')*h + u'\u2518'+'\n'
strings=[ndtotext(a)+'\n' for a in A]
strings.append(s2)
strings.insert(0,s1)
s='\n'.join(''.join(pair) for pair in zip(*map(str.splitlines, strings)))
return s
,例如:
shape = 4, 3, 6
B2=np.arange(np.prod(shape)).reshape(shape)
print(B2)
print(ndtotext(B2))
[[[ 0 1 2 3 4 5]
[ 6 7 8 9 10 11]
[12 13 14 15 16 17]]
[[18 19 20 21 22 23]
[24 25 26 27 28 29]
[30 31 32 33 34 35]]
[[36 37 38 39 40 41]
[42 43 44 45 46 47]
[48 49 50 51 52 53]]
[[54 55 56 57 58 59]
[60 61 62 63 64 65]
[66 67 68 69 70 71]]]
┌┌───────────────────┐┌───────────────────┐┌───────────────────┐┌───────────────────┐┐
│ [0 1 2 3 4 5 ] [18 19 20 21 22 23] [36 37 38 39 40 41] [54 55 56 57 58 59] │
│ [6 7 8 9 10 11] [24 25 26 27 28 29] [42 43 44 45 46 47] [60 61 62 63 64 65] │
│ [12 13 14 15 16 17] [30 31 32 33 34 35] [48 49 50 51 52 53] [66 67 68 69 70 71] │
└└───────────────────┘└───────────────────┘└───────────────────┘└───────────────────┘┘
答案 1 :(得分:0)
在每种情况下,最终尺寸的每个实例都打印在一行上。这里没有任何矛盾之处。
尝试以下各种形式:
a = np.random.rand(5, 4, 3)
print(a)
更改a
中的维数(例如,添加更多用逗号分隔的整数)。您会发现,每次打印a
时,打印对象中的每一行都将具有k
值,其中k
是a
形状中的最后一个整数。