Lesson10——NumPy 迭代陣列

cute_Learner發表於2022-02-14

NumPy 教程目錄

NumPy 迭代陣列

  NumPy 迭代器物件 numpy.nditer 提供了一種靈活訪問一個或者多個陣列元素的方式。

  迭代器最基本的任務的可以完成對陣列元素的訪問。

Example:使用 arange() 函式建立一個 2X3 陣列,並使用 nditer 對它進行迭代。

a = np.arange(6).reshape(2,3)
print('原始陣列是')
print(a)
print('迭代輸出元素')
for x in np.nditer(a):
    print(x,end=' ')
#輸出結果
"""
原始陣列是
[[0 1 2]
 [3 4 5]]
迭代輸出元素
0 1 2 3 4 5 
"""

Example:

a = np.arange(8).reshape(2,2,2)
print('原始陣列是')
print(a)
print('迭代輸出元素')
for x in np.nditer(a):
    print(x,end=' ')

#輸出結果
"""
原始陣列是
[[[0 1]
  [2 3]]

 [[4 5]
  [6 7]]]
迭代輸出元素
0 1 2 3 4 5 6 7 
"""

  以上例項不是使用標準 C 或者 Fortran 順序,選擇的順序是和陣列記憶體佈局一致的,這樣做是為了提升訪問的效率,預設是行序優先(row-major order,或者說是 C-order)。

  這反映了預設情況下只需訪問每個元素,而無需考慮其特定順序。我們可以通過迭代上述陣列的轉置來看到這一點,並與以 C 順序訪問陣列轉置的 copy 方式做對比,如下例項:

Example:

a = np.arange(6).reshape(2,3)
for x in np.nditer(a.T):
    print(x,end=' ')
print()
print(a.T.copy(order ='C'))
for x in np.nditer(a.T.copy(order ='C')):
    print(x,end=' ')
#輸出結果
"""
0 1 2 3 4 5 
[[0 3]
 [1 4]
 [2 5]]
0 3 1 4 2 5 
"""

  從上述例子可以看出, a a.T 的遍歷順序是一樣的,也就是他們在記憶體中的儲存順序也是一樣的,也就是他們在記憶體中的儲存順序也是一樣的,但是出, a.T.copy(order='C') 的遍歷結果是不同的,那是因為它和前兩種的儲存方式是不一樣的,預設是按行訪問。

控制遍歷順序

  • for x in np.nditer(a, order='F'):Fortran order,即是列序優先;
  • for x in np.nditer(a.T, order='C'):C order,即是行序優先;

Example:

a = np.arange(0,60,5) 
a = a.reshape(3,4)  
print ('原始陣列是:') 
print (a) 
print ('\n') 
print ('原始陣列的轉置是:') 
b = a.T 
print (b) 
print ('\n') 
print ('以 C 風格順序排序:') 
c = b.copy(order='C')  
print (c)
for x in np.nditer(c):  
    print (x, end=", " )
print  ('\n') 
print  ('以 F 風格順序排序:')
c = b.copy(order='F')  
print (c)
for x in np.nditer(c):  
    print (x, end=", " )
#輸出結果
"""
原始陣列是:
[[ 0  5 10 15]
 [20 25 30 35]
 [40 45 50 55]]

原始陣列的轉置是:
[[ 0 20 40]
 [ 5 25 45]
 [10 30 50]
 [15 35 55]]

以 C 風格順序排序:
[[ 0 20 40]
 [ 5 25 45]
 [10 30 50]
 [15 35 55]]
0, 20, 40, 5, 25, 45, 10, 30, 50, 15, 35, 55, 

以 F 風格順序排序:
[[ 0 20 40]
 [ 5 25 45]
 [10 30 50]
 [15 35 55]]
0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 
"""

 可以通過顯式設定,來強制 nditer 物件使用某種順序:

Example:

a = np.arange(0,60,5) 
a = a.reshape(3,4)  
print ('原始陣列是:')
print (a)
print ('\n')
print ('以 C 風格順序排序:')
for x in np.nditer(a, order =  'C'):  
    print (x, end="  " )
print ('\n')
print ('以 F 風格順序排序:')
for x in np.nditer(a, order =  'F'):  
    print (x, end=" " )
#輸出結果
"""
原始陣列是:
[[ 0  5 10 15]
 [20 25 30 35]
 [40 45 50 55]]

以 C 風格順序排序:
0  5  10  15  20  25  30  35  40  45  50  55  

以 F 風格順序排序:
0 20 40 5 25 45 10 30 50 15 35 55 
"""

修改陣列中元素的值

  nditer 物件有另一個可選引數 op_flags。 預設情況下,nditer 將視待迭代遍歷的陣列為只讀物件(read-only),為了在遍歷陣列的同時,實現對陣列元素值得修改,必須指定 read-write 或者 write-only 的模式。

Example:

x = np.arange(6).reshape((2,3))
print('The original array is')
print(x)
for tmp in np.nditer(x,op_flags =['readwrite']):
    tmp[...] = 1
print("The modified array is")
print(x)
#output result
"""
The original array is
[[0 1 2]
 [3 4 5]]
The modified array is
[[1 1 1]
 [1 1 1]]
"""

使用外部迴圈

  nditer 類的構造器擁有 flags 引數,它可以接受下列值:

引數 描述
c_index 可以跟蹤 C 順序的索引
f_index 可以跟蹤 Fortran 順序的索引
multi_index 每次迭代可以跟蹤一種索引型別
external_loop 給出的值是具有多個值的一維陣列,而不是零維陣列

  在下面的例項中,迭代器遍歷對應於每列,並組合為一維陣列。

Example:

x = np.arange(9).reshape((3,3))
print("The original array is ")
print(x)
print('The modified array is ')
for x in np.nditer(x,flags=['external_loop'],order='F'):
    print(x)
#output result
"""
The original array is 
[[0 1 2]
 [3 4 5]
 [6 7 8]]
The modified array is 
[0 3 6]
[1 4 7]
[2 5 8]
"""

廣播迭代

  如果兩個陣列是可廣播的,nditer 組合物件能夠同時迭代它們。 假設陣列 $ a $ 的維度為 \(3\times4\),陣列 \(b\) 的維度為 \(1\times4\) ,則使用以下迭代器(陣列 \(b\) 被廣播到 \(a\) 的大小)。

Example:

a = np.arange(0,60,5) 
a = a.reshape(3,4)  
print  ('First array:')
print (a)
print  ('\n')
print ('Second array:')
b = np.array([1,  2,  3,  4], dtype =  int)  
print (b)
print ('\n')
print ('The modified array :')
for x,y in np.nditer([a,b]):  
    print ("%d:%d"  %  (x,y), end=", " )
#output result
"""
First array:
[[ 0  5 10 15]
 [20 25 30 35]
 [40 45 50 55]]

Second array:
[1 2 3 4]

The modified array :
0:1, 5:2, 10:3, 15:4, 20:1, 25:2, 30:3, 35:4, 40:1, 45:2, 50:3, 55:4, 
"""

相關文章