遍历数组

Numpy提供了一个**nditer**迭代器对象,它可以配合for循环完成对数组元素的遍历

  1. import numpy as np
  2. a = np.arange(0, 60, 5).reshape(3, 4)
  3. # 使用nditer迭代器,并使用for进行遍历
  4. for x in np.nditer(a):
  5. print(x, end="\t") # 0 5 10 15 20 25 30 35 40 45 50 55

遍历顺序

在内存中,Numpy数组提供了两种存储方式,分别是C(行优先)和F(列优先)。**nditer**迭代器选择了一种与数组内存布局一致的顺序来处理具有特定存储顺序的数组,之所以这样做,是为了提升数据的访问效率。

  • 在默认情况下,当我们遍历数组中元素的时候,不需要考虑数组的存储顺序
  • **a****a.T****a**的转置)的遍历顺序是一样的,意味着它们在内存中的存储顺序是一样的

    指定遍历顺序

    可通过nditer对象的**order**参数来指定数组的遍历顺序。 ```python import numpy as np

a = np.arange(0,60,5) a = a.reshape(3,4) print(a) for x in np.nditer(a, order = ‘C’): print (x,end=”,”) # 0,5,10,15,20,25,30,35,40,45,50,55, for x in np.nditer(a, order = ‘F’): print (x,end=”,”) # 0,20,40,5,25,45,10,30,50,15,35,55,

  1. <a name="kJXhv"></a>
  2. ## 修改数组元素值
  3. `nditer` 对象提供了一个可选参数`op_flags`,它**表示能否在遍历数组时对元素进行修改**。
  4. <a name="UtZ8a"></a>
  5. ### `read-only`
  6. 只读模式,不可修改数组元素值。
  7. <a name="IS2Rm"></a>
  8. ### `read-write`
  9. 读写模式,可以修改数组元素值。
  10. <a name="HuRXM"></a>
  11. ### `write-only`
  12. 只写模式,可以修改数组元素值。
  13. ```python
  14. import numpy as np
  15. a = np.arange(0,60,5).reshape(3,4)
  16. print ("原数组是:",a)
  17. for x in np.nditer(a, op_flags=['readwrite']):
  18. x[...]=2*x
  19. print ('修改后的数组是:',a)

外部循环使用

nditer 对象的构造函数有一个“flags”参数,它可以接受以下参数值(了解即可):

参数值 描述说明
c_index 可以跟踪 C 顺序的索引。
f_index 可以跟踪 Fortran 顺序的索引。
multi_index 每次迭代都会跟踪一种索引类型。
external_loop 返回的遍历结果是具有多个值的一维数组。
  1. import numpy as np
  2. a = np.arange(0,60,5)
  3. a = a.reshape(3,4)
  4. print("原数组",a)
  5. #修改后数组
  6. for x in np.nditer(a, flags = ['external_loop'], order = 'F'):
  7. print(x)

迭代多个数组

如果两个数组都能被广播,那么**nditer**对象就可以同时对它们迭代

  1. import numpy as np
  2. a = np.arange(0, 60, 5).reshape(3, 4)
  3. print(a)
  4. b = np.array([1, 2, 3, 4], dtype=int)
  5. print(b)
  6. # 广播迭代
  7. for x, y in np.nditer([a, b]):
  8. print("%d:%d" % (x, y), end=",")

相关数组操作

数组变维操作

函数名称 函数介绍
reshape 在不改变数组元素的条件下,修改数组的形状
flat 返回的是一个迭代器,可以用 for 循环遍历其中的每一个元素
flatten 以一维数组的形式返回一份数组的副本,对副本的操作不会影响到原数组
ravel 返回一个连续的扁平数组(即展开的一维数组),与 flatten不同,它返回的是数组视图(修改视图会影响原数组)

reshape

略。

flat(遍历每一个元素)

  1. import numpy as np
  2. a = np.arange(9).reshape(3,3)
  3. for row in a:
  4. print (row)
  5. #使用flat属性:
  6. for ele in a.flat:
  7. print (ele,end=",")

flatten(创建副本后操作)

  1. # 语法格式
  2. ndarray.flatten(order='C')
  1. import numpy as np
  2. a = np.arange(8).reshape(2,4)
  3. print (a)
  4. #默认按行C风格展开的数组
  5. print (a.flatten())
  6. #以F风格顺序展开的数组
  7. print (a.flatten(order = 'F'))

ravel(返回视图后操作(将修改原始数据))

  1. # 语法格式
  2. numpy.ravel(a, order='C')
  1. import numpy as np
  2. a = np.arange(8).reshape(2,4)
  3. print ('原数组:')
  4. print (a)
  5. print ('调用 ravel 函数后:')
  6. print (a.ravel())
  7. print ('F 风格顺序调用 ravel 函数之后:')
  8. print (a.ravel(order = 'F'))

数组转置操作

函数名称 说明
transpose 将数组的维度值进行对换,比如二维数组维度(2,4)使用该方法后为(4,2)。
ndarray.T 与 transpose 方法相同
rollaxis 沿着指定的轴向后滚动至规定的位置
swapaxes 对数组的轴进行对换

transpose

  1. numpy.transpose(arr, axes)
  2. 其中:
  3. * arr:要操作的数组;
  4. * axes:可选参数,元组或整数列表,将会按该参数进行转置
  1. import numpy as np
  2. a = np.arange(12).reshape(3,4)
  3. print (a)
  4. print (np.transpose(a))

.Ttranspose的简写)

ndarray.T 的使用方法与其类似。

rollaxis

  1. numpy.rollaxis(arr, axis, start)
  2. 其中:
  3. * arr:要操作的数组;
  4. * axis:沿着哪条轴向后滚动;
  5. * start:默认以0轴开始,可以根据数组维度调整它的值

swapaxes

  1. numpy.swapaxes(arr, axis1, axis2)
  2. 其中:
  3. * arr:要操作的数组
  4. * axis1:要交换的数据轴1
  5. * axis2:要交换的数据轴2
  1. import numpy as np
  2. # 创建了三维的 ndarray
  3. a = np.arange(27).reshape(3, 3, 3)
  4. print(a)
  5. # 对换0轴与2轴
  6. print(np.swapaxes(a, 2, 0))

修改数组维度操作

函数名称 描述说明
broadcast(arr1, arr2) 返回值是数组被广播后的对象
broadcast_to 将数组广播为新的形状
expand_dims 扩展数组的形状
squeeze 从数组的形状中删除一维项

broadcast

  1. import numpy as np
  2. a = np.array([[1], [2], [3]])
  3. b = np.array([4, 5, 6])
  4. # 对b广播a
  5. d = np.broadcast(a, b)
  6. # d它拥有 iterator 属性
  7. r, c = d.iters # 将包含a和b的d对象的迭代器分别传给r和c
  8. print(next(r), next(c))
  9. print(next(r), next(c))
  10. # 使用broadcast将a与b相加
  11. e = np.broadcast(a, b)
  12. f = np.empty(e.shape) # 创建e的数组维度用于创建新的数组
  13. f.flat = [x + y for (x, y) in e] # 遍历e中的迭代器,并传给x和y,最后将x和y相加后传给f的迭代器
  14. print(f)
  15. print(a + b)

broadcast_to

该函数将数组广播到新形状中,它是在原始数组的基础上返回一个数组的只读视图

  • 如果新形状不符合Numpy的广播规则,会抛出ValueError异常。 ```python numpy.broadcast_to(array, shape, subok)

其中:

  • array:要广播的数组
  • shape:要广播的形状,参数类型是tulpe类型
  • subok:可选参数,True表示子类为passed-through,否则返回的数组将被强制为base-class数组(默认) python import numpy as np

a = np.arange(4).reshape(1, 4) print(“原数组”, a) print(‘调用 broadcast_to 函数之后:’) print(np.broadcast_to(a, (4, 4)))

  1. <a name="kTBTY"></a>
  2. ### `expand_dims`(升维)
  3. 在**指定位置插入新的轴**,从而**扩展数组的维度**。
  4. ```python
  5. numpy.expand_dims(arr, axis)
  6. 其中:
  7. * arr:输入的数组
  8. * axis:新轴插入的位置
  1. import numpy as np
  2. x = np.array(([1, 2], [3, 4])) # tuple将转换成数组类型
  3. print('数组 x:')
  4. print(x)
  5. # 在 0 轴处插入新的轴
  6. y = np.expand_dims(x, axis=0) # y被转换成了三维数组
  7. print('数组 y:')
  8. print(y)
  9. print('数组 x 和 y 的形状:')
  10. print(x.shape, y.shape)

squeeze(降维)

删除数组中维度的一项

  1. numpy.squeeze(arr, axis)
  2. 其中:
  3. * arr:输入的数组
  4. * axis:取值为整数或整数的元组,用于指定需要删除的维度所在轴。若为None,则删除维度中所有为1的项
  1. import numpy as np
  2. a = np.arange(9).reshape(1, 3, 3)
  3. print(a)
  4. b = np.squeeze(a)
  5. print(b)
  6. print('数组 a 和 b 的形状:')
  7. print(a.shape, b.shape)

连接和分割数组操作

类型 函数名称 描述说明
连接数组方法 concatenate 沿指定轴连接两个或者多个相同形状的数组
stack 沿着新的轴连接一系列数组
hstack 水平顺序堆叠序列中数组(列方向
vstack 垂直方向堆叠序列中数组(行方向
分割数组方法 split 沿着指定轴,将一个数组分割为多个子数组
hsplit 将一个数组水平分割为多个子数组按列
vsplit 将一个数组垂直分割为多个子数组按行

连接数组操作

concatenate

  1. numpy.concatenate((a1, a2, ...), axis)
  2. 其中:
  3. * a1,a2,...:表示一系列相同类型的数组
  4. * axis:沿着该参数指定的轴连接数组,默认为0
  1. import numpy as np
  2. #创建数组a
  3. a = np.array([[10,20],[30,40]])
  4. print (a)
  5. #创建数组b
  6. b = np.array([[50,60],[70,80]])
  7. print (b)
  8. #沿轴 0 连接两个数组
  9. print (np.concatenate((a,b)))
  10. #沿轴 1 连接两个数组
  11. print (np.concatenate((a,b),axis = 1))
  12. # 输出结果如下:
  13. """
  14. #a
  15. [[10 20]
  16. [30 40]]
  17. #b
  18. [[50 60]
  19. [70 80]]
  20. #axis=0沿着垂直方向
  21. [[10 20]
  22. [30 40]
  23. [50 60]
  24. [70 80]]
  25. #axis=1沿着水平方向
  26. [[10 20 50 60]
  27. [30 40 70 80]]
  28. """

stack

concatenate功能类似。

hstack

略。

vstack

  1. import numpy as np
  2. a = np.array([[1,2],[3,4]])
  3. b = np.array([[5,6],[7,8]])
  4. #垂直堆叠
  5. c = np.vstack((a,b))
  6. print (c)

分割数组操作

split

  1. numpy.split(arr, indices_or_sections, axis)
  2. 其中:
  3. * arr:要分割的数组
  4. * indices_or_sections:若是一个整数,表示用该整数平均切分;若是一个数组,表示沿着轴切分的位置(左开右闭)
  5. * axis:默认为0,表示横向切分;为1时,表示纵向切分
  1. import numpy as np
  2. a = np.arange(6)
  3. # 原数组
  4. print(a)
  5. # 将数组分为二个形状大小相等的子数组
  6. b = np.split(a, 2)
  7. print(b)
  8. # 将数组在一维数组中标明位置分割
  9. c = np.split(a, [3, 4])
  10. print(c)
  11. # 输出结果如下:
  12. """
  13. [0 1 2 3 4 5]
  14. [array([0, 1, 2]), array([3, 4, 5])]
  15. [array([0, 1, 2]), array([3]), array([4, 5])]
  16. """

hsplit

  1. import numpy as np
  2. # arr1数组
  3. arr1 = np.floor(10 * np.random.random((2, 6)))
  4. print(arr1)
  5. # 拆分后数组
  6. print(np.hsplit(arr1, 3))
  7. # 输出结果如下:
  8. """
  9. [[2. 3. 3. 7. 0. 5.]
  10. [9. 5. 2. 9. 5. 9.]]
  11. [array([[2., 3.],[9., 5.]]),
  12. array([[3., 7.],[2., 9.]]),
  13. array([[0., 5.],[5., 9.]])]
  14. """

vsplit

略。