[python][科学计算][numpy]使用指南

栏目: Python · 发布时间: 5年前

内容简介:最后一次更新日期: 2019/4/13点击下方链接可前往各小节使用指南1 - 属性 (数据类型,形状,维数,元素数,元素大小,字节数,顺序)

最后一次更新日期: 2019/4/13

NumPy 是一个 Python 包。 它代表 “Numeric Python”。 它是一个由多维数组对象(ndarray)和用于处理数组的例程集合组成的库。 使用NumPy,开发人员可以执行以下操作:

  • 数组的算数和逻辑运算。
  • 傅立叶变换和用于图形操作的例程。
  • 与线性代数有关的操作。 NumPy 拥有线性代数和随机数生成的内置函数。

使用前先导入模块: import numpy as np

点击下方链接可前往各小节

使用指南1 - 属性 (数据类型,形状,维数,元素数,元素大小,字节数,顺序)

使用指南2 - 创建 (从已有数据创建,快速填充)

使用指南3 - 运算 (数组/集合/位/字符串/统计/线性代数运算,常量,广播)

使用指南4 - 查找 (索引,遍历,抽样,视图)

使用指南5 - 变更 (更新,扩增,删除)

使用指南6 - 重构 (类型转换,重塑,排序,去重,拆分)

使用指南7 - 读写 (保存,读取)

使用指南8 - 其他 (矩阵类型,张量运算,傅里叶变换,图像处理)

一. 属性

1. 数据类型

ndarray.dtype

numpy常用数据类型

dtype type 类型名 说明
np.dtype('bool') np.bool_ 布尔类型 True or False
np.dtype('int8' | 'i1' | 'b') np.int8 8位整数(字节) -2^7 to 2^7-1
np.dtype('int16' | 'i2') np.int16 16位整数 -2^15 to 2^15-1
np.dtype('int32' | 'i4' | 'int' | 'i') np.int32 32位整数 -2^31 to 2^31-1
np.dtype('int64' | 'i8') np.int64 64位整数 -2^63 to 2^63-1
np.dtype('uint8' | 'u1' | 'B') np.uint8 8位无符号整数 0 to 2^8-1
np.dtype('uint16' | 'u2') np.uint16 16位无符号整数 0 to 2^16-1
np.dtype('uint32' | 'u4' | 'uint') np.uint32 32位无符号整数 0 to 2^32-1
np.dtype('uint64' | 'u8') np.uint64 64位无符号整数 0 to 2^64-1
np.dtype('float16' | 'f2') np.float16 半精度浮点数 1符号位+5指数位+10尾数位
np.dtype('float32' | 'f4' | 'f') np.float32 单精度浮点数 1符号位+8指数位+23尾数位
np.dtype('float64' | 'f8' | 'float' | 'd') np.float64 双精度浮点数 1符号位+11指数位+52尾数位
np.dtype('complex64' | 'c8') np.complex64 64位复数 双32位浮点数(实部+虚部)
np.dtype('complex128' | 'complex' | 'c16') np.complex128 128位复数 双64位浮点数(实部+虚部)
np.dtype('object' | 'O') object 对象 可用于存储引用类型
np.dtype('string_' | 'S' | 'S1' | 'S2' ...) np.bytes_ 定长字符串 需要声明长度
np.dtype('unicode' | 'U' | 'U1' | 'U2' ...) np.str_ 定长Unicode字符串 需要声明长度
np.dtype('datetime64' | 'M') np.datetime64 日期时间 可指定日期单位

注:同一种dtype可以通过多种字符串标识去声明的,见 | 分割的多项。 python的基本数据类型可被numpy识别,转换为对应类别下的默认类型, int 对应 np.int32float 对应 np.float64complex 对应 np.complex128

dtype与type的转换

获取dtype的type: dtype.type 通过type创建dtype: dtype(type)

dtype与type可以直接进行比较

np.dtype('int')==np.int32

dtype的字节顺序

np.dtype('<i8')

字节顺序是通过对数据类型预先设定"<"或">"来决定的。"<"意味着小端法(最小值存储在最小的地址,即低位组放在最前面)。">"意味着大端法(最重要的字节存储在最小的地址,即高位组放在最前面)。一般情况下采用默认设置即可。

2. 元素个数

ndarray.size

3. 元素的子节大小

ndarray.itemsize

由数据类型决定,每8位为1字节。

4. 总字节数

ndarray.nbytes

由size和itemsize计算得到。

5. 维数

ndarray.ndim

ndarray的每一个维度被描述为轴(axis),维数即轴的个数,轴的编号从0开始。

6. 形状

ndarray.shape

tuple类型,每一位对应到每个轴,例如0轴方向长为10个元素,1轴方向长为5个元素的数组形状为 (10,5)

7. 内存布局

numpy有两种数据存储的方式,一种是C-order,即行序优先,另一种是Fortran-order,即列序优先,未显式指定时默认C-order。由于内部的优化机制,多数情况下两者在使用时没有明显区别,一般不需要更改默认设置。 但在使用诸如迭代器之类的方法时,可通过 order 参数指定元素访问顺序,其中'C'指代C-order,'F'指代F-order。

二. 创建

1. 从已有数据创建

(1). 从list创建

#一维数组
a1=np.array([1,2])
#二维数组
a2=np.array([[1.,2.],[3.,4.]])
#三维数组
a3=np.array([[[1,2],[3,4]],[[5,6],[7,8]]],dtype='float')
#以此类推
复制代码

从嵌套的序列结构创建多维数组时,最外面的一层对应数组第一个轴,以此类推。 第二个参数 dtype 可以显式声明数组的数据类型,可传入 dtypetype (python基础数据类型或numpy数据类型)或与 dtype 对应的字符串标识,不声明的情况下会根据传入数据自动采用最合适的数据类型。

(2). 从tuple创建

a=np.array((1,2))
复制代码

list 是一样, listtuple 的嵌套同理。

(3). 从dict创建(不合适)

In [138]: np.array({"x1":1,"x2":2})
Out[138]: array({'x1': 1, 'x2': 2}, dtype=object)
复制代码

无法正常转换,整个 dict 会作为一个对象存入数组,可以尝试用pandas库去处理。

(4). 从其他类数组结构中创建,如PIL的图片类型

from PIL import Image
image= Image.open("D:\\test.jpg")
a=np.asarray(image)
复制代码

np.asarray 在多数情况下与 np.array 等效,区别在于 np.asarray 会避免没有必要的重复创建,当数据源同样是 ndarraydtype 无变化时,不会返回新的数组。

2. 快速填充

(1). 0/1填充

#填充0
a1=np.zeros((2,2))
#填充1
a2=np.ones((2,2))
#声明类型
a2=np.ones((2,2),dtype='int')
复制代码

第一个参数 shape 为数组形状,必须赋值; 默认数据类型为 float64 ,可通过 dtype 参数指定类型。

(2). 对角矩阵

In [151]: np.eye(3,3,0)
Out[151]: 
array([[ 1.,  0.,  0.],
       [ 0.,  1.,  0.],
       [ 0.,  0.,  1.]])

In [150]: np.eye(3,3,1)
Out[150]: 
array([[ 0.,  1.,  0.],
       [ 0.,  0.,  1.],
       [ 0.,  0.,  0.]])
复制代码

第一个参数 N 为0轴长度,必须赋值; 第二个参数 M 为1轴长度,不赋值时与 N 值一致; 第一个参数 k 为对角线偏移量,默认0,正数向上偏移,负数向下偏移。 对角矩阵固定为二维数组。

(3).单位矩阵

In [5]: np.identity(1)
Out[5]: array([[ 1.]])

In [6]: np.identity(2)
Out[6]: 
array([[ 1.,  0.],
       [ 0.,  1.]])
复制代码

第一个参数 n 为统一的轴长度,必须赋值; 单位矩阵固定为二维数组。

(4). 指定值填充

a=np.full((2,2),'a')
复制代码

第一个参数 shape 为数组形状,必须赋值; 第二个参数 fill_value 为填充值,必须赋值。

(5). 空值填充

a=np.empty((2,2))
复制代码

第一个参数 shape 为数组形状,必须赋值。 创建一个数组,但不初始化其中的值。

(6). 参考已有数组形状创建

a=np.ones((2,2))
a2=np.zeros_like(a)
a3=np.full_like(a,'a')
复制代码

(7). 等差数列

#方法一
a1=np.arange(0,100,1)
#方法二
a2=np.linspace(0,99,100)
复制代码

方法一,类似 range 函数,默认 int 类型, 三个参数分别为:开始,结束,步长(区间前闭后开); 方法二,线性空间,默认 float 类型, 和 range 不一样,结束值默认是包含于区间的, 且第三个参数不是步长而是元素个数。

(8). 随机数

浮点数随机数

#[0,1]浮点随机数
a1=np.random.rand(5,5)
#标准正太分布随机数
a2=np.random.randn(5,5)
复制代码

需要注意的是,上面的方法不通过 shape 而通过全部参数来定义数组形状。

整数随机数

#产生指定范围的整数随机数
a=np.random.randint(0,10,(4,3))
复制代码

第一个参数 low 表示区间下限,必须赋值; 第二个参数 high 表示区间上限,未赋值时会将 low 作为上限,0作为下限; 第三个参数 size 表示数组形状,未赋值时函数会返回标量值。

正态分布随机数

a=np.random.normal(100,10,(4,3))
复制代码

第一个参数 loc 表示分布的平均值; 第二个参数 scale 表示分布的标准偏差; 第三个参数 size 表示数组形状,未赋值时函数会返回标量值。

均匀分布随机数

a=np.random.uniform(0,10,(4,3))
复制代码

第一个参数 low 表示区间下限,必须赋值; 第二个参数 high 表示区间上限,未赋值时会将 low 作为上限,0作为下限; 第三个参数 size 表示数组形状,未赋值时函数会返回标量值。

泊松分布随机数

a=np.random.poisson(1.0,(4,3))
复制代码

第一个参数 lam 是lambda系数; 第二个参数 size 表示数组形状,未赋值时函数会返回标量值。

(9). 网格数据

In [26]: X,Y=np.mgrid[1:2:2j,1:3:3j]

In [27]: X
Out[27]: 
array([[ 1.,  1.,  1.],
       [ 2.,  2.,  2.]])

In [28]: Y
Out[28]: 
array([[ 1.,  2.,  3.],
       [ 1.,  2.,  3.]])

In [33]: x=np.array([1,2,3])

In [34]: y=np.array([1,2])

In [35]: X,Y=np.meshgrid(x,y)

In [36]: X
Out[36]: 
array([[1, 2, 3],
       [1, 2, 3]])

In [37]: Y
Out[37]: 
array([[1, 1, 1],
       [2, 2, 2]])
复制代码

np.mgrid 通过索引选取的方式直接获得网格数据, j 表示步数且包含停止值,去掉 j 该项代表步长且不包含停止值,第一个索引位声明 X 在轴0方向上的增长,第二个索引位声明 Y 在轴1方向上的增长,然后将 XY 广播至相同大小,存在更多索引位时以此类推。 np.meshgrid 在已有数据基础上构造网格数据,跟 mgrid 不一样,第一个参数是表示第二个轴方向上增长的向量,第二个参数对应第一个轴,第三个参数对应第三个轴,第四个参数对应第四个轴,之后以此类推。 网格数据常用于绘图。

(10). 复制

a1=np.zeros((2,2))
a2=np.copy(a1)
复制代码

三. 运算

1. 运算符

numpy对python中的运算符作了重载,可通过同样的用法实现数组运算。

数组与标量值的运算

In [191]: a=np.arange(0, 4)

In [192]: a
Out[192]: array([0, 1, 2, 3])

In [193]: a+1
Out[193]: array([1, 2, 3, 4])

In [193]: a*2
Out[194]: array([0, 2, 4, 6])
复制代码

数组与数组的运算

In [199]: a1=np.arange(0, 4);a2=np.arange(4, 8)

In [200]: a1,a2
Out[200]: (array([0, 1, 2, 3]), array([4, 5, 6, 7]))

In [201]: a1+a2
Out[201]: array([ 4,  6,  8, 10])

In [202]: a2**a1
Out[202]: array([  1,   5,  36, 343], dtype=int32)
复制代码

2. 标量值函数

标量值函数会对数组中每一个元素进行同样的计算。

一元函数

函数 作用 说明
np.abs 绝对值 计算浮点数/整数/复数的绝对值
np.fabs 绝对值 计算浮点数/整数的绝对值,速度更快(?)
np.sqrt 平方根 x^0.5
np.square 平方 x^2
np.log 自然对数 -
np.log2 2为底的对数 -
np.log10 10为底的对数
np.log1p x+1的自然对数 用于数值过小时保证计算的有效性
np.ceil 向上取整 -
np.floor 向下取整 -
np.rint 舍入取整 -
np.around 舍入指定位数 第二个参数 decimals 为舍入位数
np.exp 自然指数 e^x
np.sign 符号值 三种值:1(正)、0(0)、-1(负)
np.modf 拆分小数和整数部分 以两个独立的数组方式返回
np.isnan 判断是否为NaN 返回bool型数组
np.isfinite 判断是否是有穷 值非inf,非NaN;返回bool型数组
np.isinf 判断是否是有穷 值为inf或-inf;返回bool型数组
np.sin,np.sinh 正弦,双曲正弦 -
np.cos,np.cosh 余弦,双曲余弦 -
np.tan,np.tanh 正切,双曲正切 -
np.arcsin,np.arcsinh 反正弦,反双曲正弦 -
np.arccos,np.arccosh 反余弦,反双曲余弦 -
np.arctan,np.arctanh 反正切,反双曲正切 -
np.logical_not 逻辑非 -

多元函数

函数 作用 说明
np.add(a1,a2) 相加 a1+a2
np.sub(a1,a2) 相减 a1-a2
np.multiply(a1,a2) 相乘 a1*a2
np.divide(a1,a2) 相除 a1/a2
np.power(a1,a2) 乘幂 a1**a2
np.floor_divide(a1,a2) 整除 a1//a2
np.mod(a1,a2) 取模 a1%a2
np.maxinum(a1,a2,a3) 最大值 逐个元素进行比较,返回全部最大值的数组
np.fmax(a1,a2,a3) 最大值(忽略NaN) 逐个元素进行比较,返回全部最大值的数组
np.mininum(a1,a2,a3) 最小值 逐个元素进行比较,返回全部最小值的数组
np.fmin(a1,a2,a3) 最小值(忽略NaN) 逐个元素进行比较,返回全部最小值的数组
np.copysign(a1,a2) 复制符号 将a2的符号复制到a1中
np.greater(a1,a2) 大于 a1>a2
np.greater_equal(a1,a2) 大于等于 a1>=a2
np.less(a1,a2) 小于 a1<a2
np.less_equal(a1,a2) 小于等于 a1<=a2
np.equal(a1,a2) 等于 a1==a2
np.not_equal(a1,a2) 不等于 a1!=a2
np.logical_and(a1,a2) 逻辑与 -
np.logical_or(a1,a2) 逻辑或 -
np.logical_xor(a1,a2) 逻辑异或 -

3. 聚合函数

聚合函数会减少数组的维数,通常可以指定一个轴方向 axis 进行聚合,结果数组会减少一个维度,不指定方向时会在所有轴方向上聚合,结果为一个标量值。 大多数既可以静态调用,也可以直接调用ndarray对象的方法。

函数 作用 说明
np.sum 求和 -
np.mean 平均值 -
np.max 最大值 -
np.min 最小值 -
np.prod 连乘 -
np.any 至少一个为True 返回True/False
np.all 全部为True 返回True/False

np.maxnp.min 有对应的 np.argmaxnp.argmin 的方法用于返回索引,详见查找章节。

(以下是部分示例)

np.sum

In [313]: a=np.array([[1,3],[4,2]])

In [314]: a
Out[314]: 
array([[1, 3],
       [4, 2]])

In [315]: a.sum()
Out[315]: 10

In [316]: a.sum(axis=0)
Out[316]: array([5, 5])

In [317]: a.sum(axis=1)
Out[317]: array([4, 6])
复制代码

np.all

In [322]: a=np.array([[True,False],[True,True]])

In [323]: a
Out[323]: 
array([[ True, False],
       [ True,  True]], dtype=bool)

In [324]: a.all()
Out[324]: False

In [325]: a.all(axis=0)
Out[325]: array([ True, False], dtype=bool)
复制代码

4. 复合统计函数

函数 作用 说明
np.cumsum 累加 -
np.cumprod 累乘 -
np.std 标准差 ((a-a.mean())**2).sum()/a.size
np.var 方差 np.sqrt(((a-a.mean())**2).sum()/a.size)
np.average 加权平均数 第三个参数 weights 为权重;ndarray无对应方法
np.bincount 分箱计数 只支持整数,分箱区间根据最大最小值自动生成,间隔为1
np.histogram 直方图统计 第二个参数 bins 指定分箱方式,比 np.bincount 更灵活

(以下是部分示例)

np.cumsum

In [317]: a=np.array([[1,3],[4,2]])

In [319]: a.cumsum()
Out[319]: array([ 1,  4,  8, 10], dtype=int32)

In [320]: a.cumsum(axis=0)
Out[320]: 
array([[1, 3],
       [5, 5]], dtype=int32)
复制代码

np.average

In [331]: a=np.array([[1,3],[4,2]])

In [332]: w=np.array([[0.4,0.1],[0.2,0.3]])

In [333]: np.average(a)
Out[333]: 2.5

In [334]: np.average(a,weights=w)
Out[334]: 2.1000000000000001
复制代码

5. 字符串函数

函数 作用 说明
np.char.add 字符串相加 逐个元素执行字符串相加
np.char.multiply 字符串重复 第二个参数 i 为重复次数
np.char.center 字符串居中 第二个参数 width 为长度,第三个参数 fillchar 为填充字符
np.char.capitalize 首字母大写 -
np.char.title 单词首字母大写 -
np.char.lower 转换为小写 -
np.char.upper 转换为大写 -
np.char.split 字符串分割 第二个参数 sep 为分隔符,返回list数组
np.char.splitlines 行分割 以换行符分割,返回list数组
np.char.strip 移除头尾指定字符 第二个参数 chars 为需要移除的字符
np.char.join 以指定分隔符拼接字符 第一个参数 sep 为分隔符
np.char.replace 替换字符串 第二个参数 old 为旧字符串,第三个参数 new 为新字符串
np.char.decode 解码 对每个元素调用str.decode
np.char.encode 编码 对每个元素调用str.encode

(以下是部分示例)

np.char.add

In [301]: a1=np.array(['a','b']);a2=np.array(['c','d'])

In [302]: np.char.add(a1,a2)
Out[302]: 
array(['ac', 'bd'],
      dtype='<U2')
复制代码

np.char.multiply

In [303]: a=np.array(['a','b'])

In [304]: np.char.multiply(a,3)
Out[304]: 
array(['aaa', 'bbb'],
      dtype='<U3')
复制代码

np.char.center

In [305]: a=np.array(['a','b'])

In [306]: np.char.center(a,10,'*')
Out[306]: 
array(['****a*****', '****b*****'],
      dtype='<U10')
复制代码

np.char.split

In [307]: a=np.array(['a,b','c,d'])

In [308]: np.char.split(a,',')
Out[308]: array([list(['a', 'b']), list(['c', 'd'])], dtype=object)
复制代码

np.char.join

In [309]: a=np.array(['ab','cd'])

In [310]: np.char.join(',',a)
Out[310]: 
array(['a,b', 'c,d'],
      dtype='<U3')

In [311]: a=np.array(['a,b','c,d'])

In [312]: np.char.join(',',np.char.split(a,','))
Out[312]: 
array(['a,b', 'c,d'],
      dtype='<U3')
复制代码

注意,该方法无法实现多维数组的聚合计算。在数组元素为字符串时,会对字符串中的每个元素进行拼接;在数组元素为字符串序列时,会对序列中的字符串进行拼接。 与split互为逆运算。

6. 线性代数运算

函数 作用 说明
np.dot(a1,a2) 点乘 多维数组下会将a1的最后一个轴和a2的倒数第二个轴作为向量的维度,可视作向量的栈
np.vdot(a1,a2) 向量点乘 高维数组会被展开计算
np.inner(a1,a2) 向量内积 多维数组会将最后一个轴作为向量的维度
np.matmul(a1,a2) 矩阵乘积 多维数组下会将最后两个轴作为矩阵的维度,可视作元素是矩阵的数组
np.linalg.det(a) 行列式 行列式描述的是矩阵所表示的线性变换对“体积”的影响
np.linalg.solve(A,b) 求解线性方程组 求解线性方程组Ax = b,A为系数矩阵(方阵),b为常数矩阵
np.linalg.lstsq(A,b) 求解线性方程组 求线性方程组Ax = b的最小二乘解,A为系数矩阵,b为常数矩阵
np.linalg.inv(a) 逆矩阵 AB=BA=E,E为单位矩阵,则B为A的逆矩阵
np.linalg.pinv(a) 广义逆矩阵 可以输入非方阵
np.linalg.eig(a) 特征值和特征向量 返回两个数组
np.linalg.qr(a) 正交分解 -
np.linalg.svd(a) 奇异值分解 -
ndarray.T 转置 对一维数组转置是无效的

(以下是部分示例)

np.dot, np.vdot, np.inner, np.matmul

In [339]: a=np.array([1,2]);b=np.array([[3,4],[5,6]])

In [340]: np.dot(a,a), np.vdot(a,a), np.inner(a,a), np.matmul(a,a)
Out[340]: (5, 5, 5, 5)

In [341]: np.dot(a,b), np.dot(b,a)
Out[341]: (array([13, 16]), array([11, 17]))

In [342]: np.vdot(a,b)
Traceback (most recent call last):

  File "<ipython-input-358-f2388a21d848>", line 1, in <module>
    np.vdot(a,b)

ValueError: cannot reshape array of size 4 into shape (2,)

In [343]: np.inner(a,b), np.inner(b,a)
Out[343]: (array([11, 17]), array([11, 17]))

In [344]: np.matmul(a,b), np.matmul(b,a)
Out[344]: (array([13, 16]), array([11, 17]))

In [345]: np.dot(b,b)
Out[345]: 
array([[29, 36],
       [45, 56]])

In [346]: np.vdot(b,b)
Out[346]: 86

In [347]: np.inner(b,b)
Out[347]: 
array([[25, 39],
       [39, 61]])

In [348]: np.matmul(b,b)
Out[348]: 
array([[29, 36],
       [45, 56]])
复制代码

这四个方法执行的运算都基于向量内积,非常相似,但在具体行为上有区别,很容易混淆。

dot 是将数组看作向量的集合,第一个数组的最后一个轴和第二个数组的倒数第二个轴作为向量的轴,抽取每个向量两两匹配进行点积运算,因此这两个轴需要长度一致。例如, a.shape=(2,3,4)b.shape=(5,4,6) ,则 np.dot(a,b).shape=(2,3,5,6) 。在二维情况下,实现的计算即是矩阵乘法。

vdot 的规则比较简单,会将数组展开为向量再计算点积,要求数组的 size 一致。

innerdot 类似,但规则更简单,两个数组均以最后一个轴作为向量的轴,即最后一个轴长度要保持一致。例如, a.shape=(2,3,4)b.shape=(5,6,4) ,则 np.dot(a,b).shape=(2,3,5,6)

matmul 的计算针对矩阵,和 dot 一样在二维情况下表示矩阵乘法,二维以上视作元素为矩阵的数组,参与运算的两个数组均以最后两个轴作为矩阵的轴,逐个元素进行矩阵乘法并向其他轴方向上广播。例如, a.shape=(1,3,4)b.shape=(5,4,6) ,则 np.dot(a,b).shape=(5,3,6) ,如果 a.shape=(2,3,4) ,计算就会报错,因为不满足广播的规则。

7. 集合运算

函数 作用 说明
np.intersect1d(x, y) 交集 xy 的公共元素
np.union1d(x, y) 并集 xy 的所有元素
np.setdiff1d(x, y) 集合差 x 中存在, y 中不存在的元素
np.setxor1d(x, y) 集合异或 xy 的独占元素

以上方法适用于一维数组。

(以下是部分示例)

In [810]: x=np.array([1,3,4])

In [811]: y=np.array([2,3,5])

In [812]: np.intersect1d(x,y)
Out[812]: array([3])

In [813]: np.union1d(x,y)
Out[813]: array([1, 2, 3, 4, 5])

In [814]: np.setdiff1d(x,y)
Out[814]: array([1, 4])

In [815]: np.setxor1d(x,y)
Out[815]: array([1, 2, 4, 5])
复制代码

8. 位运算

函数 作用 说明
np.invert 按位取反 等效于 ~ 运算符
np.bitwise_and 按位与 等效于 & 运算符
np.bitwise_or 按位或 等效于 \| 运算符
np.bitwise_xor 按位异或 等效于 ^ 运算符
np.left_shift 位左移 等效于 << 运算符
np.right_shift 位右移 等效于 >> 运算符

(以下是部分示例)

In [858]: a=np.array([2,3,5,8,13])

In [859]: b=np.array([3,4,7,11,18])

In [860]: np.invert(a)
Out[860]: array([ -3,  -4,  -6,  -9, -14], dtype=int32)

In [862]: bin(a[4]),bin(~a[4])
Out[862]: ('0b1101', '-0b1110')

In [863]: np.bitwise_and(a,b)
Out[863]: array([2, 0, 5, 8, 0], dtype=int32)

bin(a[3]),bin(b[3]),bin(a[3]&b[3])
Out[865]: ('0b1000', '0b1011', '0b1000')
复制代码

np.invert 对于有符号整数,取对应二进制数的补码,然后 +1。二进制数形式的最高位为0表示正数,最高位为 1 表示负数。

9. 广播

numpy在进行不同形状的数组之间的计算时,会自动沿长度不足且长度为1的轴方向进行广播。当维数不一致时,会向前补齐,这要求后面的轴长度相同或是有一方长度等于1,即从后至前进行轴长比对,例如形状(a,b,c)与(b,c)可计算,(a,b,c)与(b,1)可计算,(a,b,c)与(a,b,1)可计算,(a,b,c)与(a,b)不可计算,(a,b,c)与(a,c,b)不可计算。

In [414]: a1=np.zeros((2,3,4))
         ...: a2=np.random.randint(0,10,(3,4))
         ...: a3=np.random.randint(0,10,(2,3))
         ...: a4=np.random.randint(0,10,(2,3,1))

In [418]: a1+a2
Out[418]: 
array([[[ 3.,  6.,  0.,  6.],
        [ 9.,  6.,  3.,  4.],
        [ 2.,  3.,  1.,  5.]],

       [[ 3.,  6.,  0.,  6.],
        [ 9.,  6.,  3.,  4.],
        [ 2.,  3.,  1.,  5.]]])

In [419]: a1+a3
Traceback (most recent call last):

  File "<ipython-input-419-d778f9717621>", line 1, in <module>
    a1+a3

ValueError: operands could not be broadcast together with shapes (2,3,4) (2,3) 

In [420]: a1+a4
Out[420]: 
array([[[ 0.,  0.,  0.,  0.],
        [ 9.,  9.,  9.,  9.],
        [ 1.,  1.,  1.,  1.]],

       [[ 4.,  4.,  4.,  4.],
        [ 0.,  0.,  0.,  0.],
        [ 2.,  2.,  2.,  2.]]])
复制代码

10. 常量

属性名 常量名
np.e 自然指数 2.718281828459045
np.pi 圆周率 3.141592653589793
np.euler_gamma 欧拉常数 0.5772156649015329
np.inf 无穷大 -
np.nan 非数值 -

四. 查找

1. 索引

通过 ndarray[index1,index2...] 的形式指定索引。

(1). 定位元素

In [452]: a=np.array([[1,2],[3,4]])

In [453]: a[1,1]
Out[453]: 4
复制代码

(2). 数据切片

In [454]: a[:,1]
Out[454]: array([2, 4])

In [455]: a[0:1,:]
Out[455]: array([[1, 2]])
复制代码

通过 low:high 的形式限定索引范围,区间前闭后开,上下限留空表示不限制。 标量值索引会使对应维度消失,范围索引则不会。

(3). 倒序索引

In [457]: a[-1,:-1]
Out[457]: array([3])
复制代码

-n 表示倒数第n个。

(4). 按步长选取

In [458]: a[0:2:2]
Out[458]: array([[1, 2]])
复制代码

通过 low:high:step 的形式限定索引范围和步长。

(5). 序列反转

In [459]: a[::-1]
Out[459]: 
array([[3, 4],
       [1, 2]])
复制代码

等效于 a[-1:-len(a)-1:-1]

(6). 布尔索引

In [460]: a>2
Out[460]: 
array([[False, False],
       [ True,  True]], dtype=bool)

In [461]: (a>2)&(a<4)
Out[461]: 
array([[False, False],
       [ True, False]], dtype=bool)

In [462]: a[(a>2)&(a<4)]
Out[462]: array([3])
复制代码

ndarray 应用逻辑运算符会得到布尔索引,布尔索引标识了每个元素符合逻辑判断条件的情况,使用该索引筛选数组将得到所有满足条件的元素构成的一维数组。

(7). 数组索引

In [463]: a[[0,1],[1]]
Out[463]: array([2, 4])

In [464]: a[[0,1],[1,0]]
Out[464]: array([2, 3])
复制代码

将对应轴上需要选取的索引以数组形式传入,当在多个轴上传入数组索引时,索引会被一一对应。支持list,tuple,ndarray,支持倒序索引。

(8). 综合示例

获取第二列大于2的行的第一列

In [465]: a[a[:,1]>2,0]
Out[465]: array([3])
复制代码

获取二维数组的顶角元素

In [466]: a[[0,-1]][:,[0,-1]]
Out[466]: 
array([[1, 2],
       [3, 4]])

In [467]: a[[0,0,-1,-1],[0,-1,0,-1]].reshape((2,2))
Out[467]: 
array([[1, 2],
       [3, 4]])
复制代码

此处使用的 reshape 重塑方法会在重构章节讲到。

(9). 获取布尔索引指示的位置

In [473]: a>2
Out[473]: 
array([[False, False],
       [ True,  True]], dtype=bool)

In [474]: np.where(a>2)
Out[474]: (array([1, 1], dtype=int64), array([0, 1], dtype=int64))

In [476]: a[np.where(a>2)]
Out[476]: array([3, 4])

In [477]: np.where(a>2,a,-a)
Out[477]: 
array([[-1, -2],
       [ 3,  4]])
复制代码

只输入数组时,内部会调用 nonzero 方法,返回值为一个 tuple ,包含每个轴上的索引值序列。 通过第二个参数 x 和第三个参数 y ,可以根据是否满足条件对元素进行不同的计算,返回同样形状的数组。

(10). 特定值查找

指定值序列

In [655]: a=np.arange(5)

In [656]: np.isin(a,[3,4])
Out[656]: array([False, False, False,  True,  True], dtype=bool)
复制代码

第二个参数 test_elements 为指定值序列,返回布尔索引。

排序 搜索

In [683]: a=np.array([1,3,5,7,9])

In [684]: np.searchsorted(a,[3,6])
Out[684]: array([1, 3], dtype=int64)
复制代码

第一个参数 a 为升序一维数组,非升序可通过第四个参数指定排序 sorter=np.argsort(a) ,第二个参数 v 为需要被插入 a 的值或值序列,第三个参数 side 为判断方式, 'left' 表示 a[i-1] < v <= a[i]'right' 表示 a[i-1] <= v < a[i] 。该方法将逐个搜索v的每个元素在升序数组a中合适的插入位置,返回数组索引。

空值

In [661]: a=np.array([np.nan,1,2,3])

In [662]: np.isnan(a)
Out[662]: array([ True, False, False, False], dtype=bool)
复制代码

返回布尔索引。

(11). 返回索引的方法

函数 作用 说明
np.argsort 返回排序后索引 -
np.argmax 返回最大值索引 -
np.argmin 返回最小值索引 -
np.argpartition 返回分区索引 第二个参数 kth 指定用于分区的元素索引
np.argwhere 返回符合条件的值索引 np.where 类似,但返回值的形式不太适合索引
np.where 返回符合条件的值索引 -
np.isin 返回判断元素是否在指定列表中的布尔索引 第二个参数 test_elements 为指定的值序列
np.isnan 返回判断元素是否nan值的布尔索引 -
np.searchsorted 返回待插入值在升序序列中的插入位置 第一个参数 a 为升序一维数组,第二个参数 v 为待插入值
np.digitize 返回分箱后的索引 第二个参数 bins 指定分箱方式

2. 遍历

(1). 序列结构遍历

a=np.array([[1,2],[3,4]])

for buf1 in a:
    for buf2 in buf1:
        print(buf2)
复制代码

(2). 索引遍历

for i in range(a.shape[0]):
    for j in range(a.shape[1]):
        print(a[i,j])
复制代码

以上两种方法在遍历元素时,维数是多少就要嵌套多少层循环,效率不算高,但可以在每一层循环中嵌入额外的计算。

(3). 快速迭代器遍历

for item in np.nditer(a):
    print(item)
复制代码

默认选择元素的顺序是和数组内存布局一致的,而不是使用标准C或者Fortran顺序。这是为了使用效率而设计的,这反映了默认情况下只需访问每个元素,而无需考虑其特定顺序。可通过 order 参数来指定特定的顺序来访问数组。

(4). 平铺迭代器遍历

for item in a.flat:
    print(item)
复制代码

将数组转换为1-D的迭代器。

3. 抽样

(1). 无放回抽样

In [506]: a=np.array([[1,2],[3,4]])

In [507]: idx=np.random.choice(np.arange(a.shape[0]),size=1,replace=False)

In [508]: a[idx]
Out[508]: array([[3, 4]])
复制代码

(2). 有放回抽样

In [511]: idx=np.random.randint(0,a.shape[0],size=3)

In [512]: a[idx]
Out[512]: 
array([[1, 2],
       [3, 4],
       [1, 2]])
复制代码

也可使用 np.random.choice(np.arange(a.shape[0]),size=3) 生成索引。 以上都是在单个轴方向上抽样,如果想在多个轴方向上抽样,由于抽样后必然会破坏数组结构,建议先将用于抽样的轴展开,比如用 reshape 方法,见下面的示例,关于 reshape 方法的说明在重构章节。

(3). 多轴方向抽样

In [524]: a=np.array([[[1,2],[3,4]],[[5,6],[7,8]]])

In [525]: a2=a.reshape((2,4))

In [526]: idx=np.random.randint(0,a2.shape[0],size=2)

In [527]: a2[:,idx]
Out[527]: 
array([[1, 2],
       [5, 6]])
复制代码

4. 视图

在对数组进行了索引切片后,返回的通常是原数组的一个视图,不会完整的拷贝数据,因此在这种情况下进行更新操作会影响到原数组和所有视图。 可通过 ndarray.view() 方法获得一个完整视图,等效于 ndarray[:,:,...] 。 想将视图转换为拷贝可使用 ndarray.copy() 方法。

In [542]: a=np.array([[1,2],[3,4]])

In [543]: a2=a[1,:]

In [544]: a2
Out[544]: array([3, 4])

In [545]: a2[:]=0

In [546]: a
Out[546]: 
array([[1, 2],
       [0, 0]])
复制代码

五. 变更

1. 更新

(1). 更新整个数组

In [555]: a=np.array([[1,2],[3,4]])

In [556]: a=a+1

In [557]: a
Out[557]: 
array([[2, 3],
       [4, 5]])

In [558]: a[:,:]=a-1

In [559]: a
Out[559]: 
array([[1, 2],
       [3, 4]])
复制代码

两种方式都能更新整个数组,第一种将计算得到的新数组的引用重新赋给了a,第二种根据计算得到的新数组更新了原数组中相应位置的值。

(2). 更新指定位置

In [569]: a[a>2]+=1

In [570]: a
Out[570]: 
array([[1, 2],
       [4, 5]])

In [571]: a[a>2]=a[a>2]+1

In [572]: a
Out[572]: 
array([[1, 2],
       [5, 6]])

In [573]: a[0,:]=0

In [574]: a
Out[574]: 
array([[0, 0],
       [5, 6]])

In [575]: a[a>2]=[3,4]

In [576]: a
Out[576]: 
array([[0, 0],
       [3, 4]])
复制代码

值数组形状需要与筛选后的原数组一致或遵循广播的规则。

(3). 定值填充

In [9]: a.fill(1)

In [10]: a
Out[10]: 
array([[1, 1],
       [1, 1]])
复制代码

当填充值数据类型与数组数据类型不一致时,会尝试转换,失败时才会报错。

2. 扩增

(1). 插入

In [577]: a=np.array([[1,2],[3,4]])

In [578]: np.insert(a,1,[5,6],axis=0)
Out[578]: 
array([[1, 2],
       [5, 6],
       [3, 4]])
复制代码

第二个参数 obj 是插入的位置索引,第三个参数 values 是待插入的值,需要与指定轴方向上的切片形状一致或满足广播规则,第四个参数 axis 是指定的轴。不影响原数组,返回的是一个拷贝。

(2). 附加

In [578]: np.append(a,[[5,6]],axis=0)
Out[578]: 
array([[1, 2],
       [3, 4],
       [5, 6]])
复制代码

第二个参数 values 是待插入的值,需要与指定轴方向上的切片形状一致或满足广播规则,第三个参数 axis 是指定的轴。只能将新数据附加到数组末尾。不影响原数组,返回的是一个拷贝。

(3). 堆叠

In [589]: np.c_[a,a]
Out[589]: 
array([[1, 2, 1, 2],
       [3, 4, 3, 4]])

In [590]: np.column_stack((a,a))
Out[590]: 
array([[1, 2, 1, 2],
       [3, 4, 3, 4]])

In [591]: np.concatenate((a,a),axis=1)
Out[591]: 
array([[1, 2, 1, 2],
       [3, 4, 3, 4]])

In [592]: np.r_[a,a]
Out[592]: 
array([[1, 2],
       [3, 4],
       [1, 2],
       [3, 4]])

In [593]: np.row_stack((a,a))
Out[593]: 
array([[1, 2],
       [3, 4],
       [1, 2],
       [3, 4]])

In [594]: np.concatenate((a,a),axis=0)
Out[594]: 
array([[1, 2],
       [3, 4],
       [1, 2],
       [3, 4]])

In [595]: np.stack((a,a),axis=0)
Out[595]: 
array([[[1, 2],
        [3, 4]],

       [[1, 2],
        [3, 4]]])
复制代码

np.c_np.column_stack 是沿轴1进行堆叠,其他轴长度需要相同或满足广播规则,等效于 np.concatenate(axis=1)np.r_np.row_stack 是沿轴0进行堆叠,其他轴长度需要相同或满足广播规则,等效于 np.concatenate(axis=0)np.stack 是沿新轴进行堆叠,所有轴长度需要相同或满足广播规则。

(4). 重复

In [93]: a.repeat(3,axis=1)
Out[93]: 
array([[1, 1, 1, 2, 2, 2],
       [3, 3, 3, 4, 4, 4]])

In [96]: a.reshape((2,2,1)).repeat(3,axis=2)
Out[96]: 
array([[[1, 1, 1],
        [2, 2, 2]],

       [[3, 3, 3],
        [4, 4, 4]]])

In [98]: a.repeat(3).reshape((2,2,3))
Out[98]: 
array([[[1, 1, 1],
        [2, 2, 2]],

       [[3, 3, 3],
        [4, 4, 4]]])

In [99]: np.tile(a,2)
Out[99]: 
array([[1, 2, 1, 2],
       [3, 4, 3, 4]])

In [100]: np.tile(a,[2,2])
Out[100]: 
array([[1, 2, 1, 2],
       [3, 4, 3, 4],
       [1, 2, 1, 2],
       [3, 4, 3, 4]])

In [101]: np.tile(a.ravel(),2)
Out[101]: array([1, 2, 3, 4, 1, 2, 3, 4])
复制代码

repeat 方法将数组中的元素重复,可通过 axis 参数指定轴方向,默认会将数组展开后在唯一的轴方向上重复元素。可配合 ndarray.reshape 在新轴上复制元素。 tile 方法将数组重复,注意,重复的是整个数组,不是单个元素,得到的结果中同元素不一定是紧挨着的。

3. 删除

(1). 索引筛选

In [616]: a[~(a[0]==1),:]
Out[616]: array([[3, 4]])
复制代码

通过索引筛选可得到删除指定内容的数组。

(2). 删除方法

In [617]: np.delete(a,[0,1],axis=0)
Out[617]: array([], shape=(0, 2), dtype=int32)
复制代码

通过相应方法获得删除指定索引位置内容的数组。第二个参数 obj 为索引位置,第三个参数 axis 为指定轴。

六. 重构

1. 类型转换

In [619]: a=np.array([[1,2],[3,4]])

In [620]: a.dtype
Out[620]: dtype('int32')

In [621]: a=a.astype('float64')

In [622]: a.dtype
Out[622]: dtype('float64')

In [623]: a=np.int32(a)

In [624]: a.dtype
Out[624]: dtype('int32')
复制代码

使用 ndarray.astype 方法或是使用数据类型同名方法都可以转换类型,关于numpy支持的数据类型可以查看属性章节。转换类型后返回一个新数组。

2. 重塑

(1). 改变形状

In [625]: a=np.array([[1,2],[3,4]])

In [626]: a.reshape((1,4))
Out[626]: array([[1, 2, 3, 4]])

In [627]: a.reshape((-1,4))
Out[627]: array([[1, 2, 3, 4]])
复制代码

使用tuple类型的参数声明新的形状。允许有一个新轴的大小为-1,表示自动计算。 改变前后元素数 size 需要保持一致。元素在轴上的排列是从最后一个轴开始往前面的轴方向上堆叠,见如下图示,可通过 order 参数指定其他排序方式。轴的相对位置不会改变,所以一些复杂的变形可能需要结合 transposeswapaxes 此类轴交换方法使用。

[python][科学计算][numpy]使用指南
[python][科学计算][numpy]使用指南

(2). 平铺

In [640]: a.ravel()
Out[640]: array([1, 2, 3, 4])

In [641]: a.flatten()
Out[641]: array([1, 2, 3, 4])
复制代码

将数组平铺为向量,等效于reshape((-1,)),可通过 order 参数指定其他排序方式。

(3). 转置

In [643]: a.T
Out[643]: 
array([[1, 3],
       [2, 4]])
复制代码

数组为一维时转置无效,为二维时即矩阵的转置,多于二维时交换第一个和最后一个轴。

(4). 轴交换

In [646]: a.swapaxes(0,1)
Out[646]: 
array([[1, 3],
       [2, 4]])

In [647]: a.transpose([1,0])
Out[647]: 
array([[1, 3],
       [2, 4]])
复制代码

swapaxes 一次只能指定两个轴进行交换, transpose 可以重新为所有轴排序。

3. 排序

(1). 直接排序

In [734]: a=np.array([[2,3],[1,4]])

In [735]: np.sort(a,axis=None)
Out[735]: array([1, 2, 3, 4])

In [736]: a.sort(axis=0)

In [737]: a
Out[737]: 
array([[1, 3],
       [2, 4]])
复制代码

ndarray.sort 会直接在原数组上排序,可通过第一个参数 axis 指定排序的轴,会将沿着该轴方向的每个向量单独排序,默认 -1 ,除沿最后一个轴外,指定其他轴都会在排序时生成数据的临时副本,因此沿最后一个轴排序最快。 等效方法 np.sort ,返回的是排序后的副本,还可指定 axis=None ,会将数组展开再排序。 当数组的维度具备实际含义时,直接排序会打乱数据结构,得到不被期望的结果,这种情况下需要使用间接排序。

(2). 间接排序

In [740]: a=np.array([[2,3,5],[1,1,4],[1,2,3]])

In [741]: a
Out[741]: 
array([[2, 3, 5],
       [1, 1, 4],
       [1, 2, 3]])

In [742]: idx1=np.argsort(a[:,0])

In [743]: a[idx1]
Out[743]: 
array([[1, 1, 4],
       [1, 2, 3],
       [2, 3, 5]])

In [744]: idx2=np.lexsort((a[:,0],a[:,2]))

In [745]: a[idx2]
Out[745]: 
array([[1, 2, 3],
       [1, 1, 4],
       [2, 3, 5]])
复制代码

argsort 可用于单键间接排序, lexsort 可用于多键间接排序。

(3). 随机排序

In [763]: a=np.arange(12).reshape((3,4))

In [764]: a
Out[764]: 
array([[ 0,  1,  2,  3],
       [ 4,  5,  6,  7],
       [ 8,  9, 10, 11]])

In [765]: np.random.shuffle(a)

In [766]: a
Out[766]: 
array([[ 4,  5,  6,  7],
       [ 0,  1,  2,  3],
       [ 8,  9, 10, 11]])

In [768]: idx=np.random.permutation(a.shape[1])

In [769]: a[:,idx]
Out[769]: 
array([[ 4,  5,  7,  6],
       [ 0,  1,  3,  2],
       [ 8,  9, 11, 10]])
复制代码

方法一 np.random.shuffle 只能沿第一个轴进行随机排序,方法二是通过 np.random.permutation ,该方法也只能沿第一个轴随机排序,但在输入参数 xint 类型时,会对np.arange(x)进行随机排序,可以快速生成乱序索引,再通过索引查找得到乱序数组。

(4). 分区排序

In [799]: a=np.array([2,3,1,5,4,0,8])

In [800]: np.partition(a,1)
Out[800]: array([0, 1, 3, 5, 4, 2, 8])

In [801]: np.partition(a,4)
Out[801]: array([0, 1, 2, 3, 4, 5, 8])
复制代码

分区排序是一种不完整的排序,用于不需要获取完整排序序列的情况下。该方法只保证 kth 指定位置的元素是正确排序的,其他小于该元素的元素前移,大于的后移。可用于快速找出第k大或第k小的元素。 也可通过 ndarray.partition 调用, axis 参数指定轴方向,还有对应的 np.argpartition 用于获取分区后的索引。

4. 去重

In [775]: a=np.array([3,2,2,3,1,1,4])

In [776]: np.unique(a,return_index=True,return_inverse=True,return_counts=True,axis=None)
Out[776]: 
(array([1, 2, 3, 4]),
 array([4, 1, 0, 6], dtype=int64),
 array([2, 1, 1, 2, 0, 0, 3], dtype=int64),
 array([2, 2, 2, 1], dtype=int64))

In [777]: a=np.array([[1,3],[2,4],[1,3]])

In [778]: np.unique(a,axis=0)
Out[778]: 
array([[1, 3],
       [2, 4]])
复制代码

axis 指定去重的轴,默认 None 会将数组展开后再去重。 返回值一为去重后的有序值列表; 返回值二为唯一值在原数组中的索引,仅在 return_index=True 时提供; 返回值三为根据唯一值重建原数组的索引,仅在 return_inverse=True 时提供; 返回值四为唯一值的出现计数,仅在 return_counts=True 时提供。

5. 拆分

(1). 索引拆分

In [780]: a=np.arange(9).reshape((3,3))

In [781]: a[:2,:],a[2:,:]
Out[781]: 
(array([[0, 1, 2],
        [3, 4, 5]]), array([[6, 7, 8]]))

In [782]: a[a>2],a[~(a>2)]
Out[782]: (array([3, 4, 5, 6, 7, 8]), array([0, 1, 2]))

In [790]: idx=np.random.permutation(a.shape[0])

In [791]: sp_idx=int(a.shape[0]*0.8)

In [792]: a[idx[:sp_idx]],a[idx[sp_idx:]]
Out[792]: 
(array([[3, 4, 5],
        [6, 7, 8]]), array([[0, 1, 2]]))
复制代码

最灵活的方式,复杂的拆分需要写较多的代码,可使用其他方法配合生成用于拆分的索引,比如使用 np.digitize 进行分箱。

(2). 拆分方法

In [802]: a=np.arange(8)

In [803]: np.split(a,2)
Out[803]: [array([0, 1, 2, 3]), array([4, 5, 6, 7])]

In [805]: np.split(a,[2,5])
Out[805]: [array([0, 1]), array([2, 3, 4]), array([5, 6, 7])]

In [806]: a=a.reshape((2,4))

In [807]: a
Out[807]: 
array([[0, 1, 2, 3],
       [4, 5, 6, 7]])

In [809]: np.split(a,[1,3],axis=1)
Out[809]: 
[array([[0],
        [4]]), array([[1, 2],
        [5, 6]]), array([[3],
        [7]])]
复制代码

第二个参数 indices_or_sections 指定分割方式, int 类型表示等分数量,一维数组类型表示用于分割的索引值,例如 [2,5] 表示分割为 a[:2],a[2:5],a[5:] 。 第三个参数 axis 可以指定轴方向。

七. 读写

numpy 可使用专用的 npynpz 格式或常见的 txt 格式存储 ndarray 的数据。

1. 保存

In [803]: a=np.array([1,2,3,4])

In [804]: b=np.array([[5,6],[7,8]])

In [805]: np.save('D:\\out.npy',a)

In [806]: np.savez('D:\\out.npz',a,b=b)

In [807]: np.savetxt('D:\\out.txt',b,fmt="%d", delimiter=",")
复制代码

拓展名可以省略,会自动补全。无法自动创建文件夹。 save 用于保存单个数组为 npy 格式文件。 savez 用于保存多个数组为 npz 格式文件,没有使用关键字参数传递的数组会自动命名为 arr_0,arr_1,...savetxt 用于保存单个数组为 txt 格式文件,参数 fmt 指定保存时的字符串转换,参数 delimiter 指定分隔符,注意在读取时也需要指定分隔符。分隔符的设置对一维数组无效,二维以上的数组不适合用该方法保存。

2. 加载

In [835]: np.load('D:\\out.npy')
Out[835]: array([1, 2, 3, 4])

In [836]: npz=np.load('D:\\out.npz')

In [837]: npz['arr_0']
Out[837]: array([1, 2, 3, 4])

In [838]: npz['b']
Out[838]: 
array([[5, 6],
       [7, 8]])

In [841]: np.loadtxt('D:\\out.txt',dtype='int',delimiter=',')
Out[841]: 
array([[5, 6],
       [7, 8]])
复制代码

np.load 读取 npz 格式文件会得到一个 NpzFile 对象,之后通过保存时设置的名称进行 [] 索引可以取得每一个数组。

八. 其他

1. 矩阵类型

numpy提供了一个专用的矩阵对象 matrix ,是基于 ndarray 作了进一步的封装得到的,能够更加快捷地进行一些矩阵相关的运算,但相比 ndarray 没有性能上的优势且维数限制在二维,并不推荐使用。

(1). 创建

In [79]: m1=np.matrix([[1,2],[3,4]])

In [81]: m1
Out[81]: 
matrix([[1, 2],
        [3, 4]])
复制代码

创建方式与 ndarray 类似。

(2). 与 ndarray 的相互转换

In [82]: a=np.array([[5,6],[7,8]])

In [83]: m2=np.matrix([[5,6],[7,8]])

In [84]: np.asmatrix(a)
Out[84]: 
matrix([[5, 6],
        [7, 8]])

In [85]: np.matrix(a)
Out[85]: 
matrix([[5, 6],
        [7, 8]])

In [88]: np.asarray(m2)
Out[88]: 
array([[5, 6],
       [7, 8]])

In [89]: np.array(m2)
Out[89]: 
array([[5, 6],
       [7, 8]])

In [90]: m2.base
Out[90]: 
array([[5, 6],
       [7, 8]])
复制代码

(3). 矩阵运算

In [94]: m1*m2
Out[94]: 
matrix([[19, 22],
        [43, 50]])

In [95]: m1.I
Out[95]: 
matrix([[-2. ,  1. ],
        [ 1.5, -0.5]])

In [96]: m1.T
Out[96]: 
matrix([[1, 3],
        [2, 4]])
复制代码

运算符 * 用在 matrix 上表示矩阵乘法,等效于 np.dot(m1,m2) ,要实现元素相乘需要使用 np.multiply(m1,m2)matrix.T 表示转置矩阵, matrix.I 表示逆矩阵。 matrix 可以使用大部分 ndarray 的方法,比如 maxsumsort 等。

2. 张量运算

张量是向量、矩阵这类概念的推广,标量是0阶张量,向量是1阶张量,矩阵是2阶张量。 numpy提供了广义的张量点积运算 np.tensordot

In [2]: a=np.arange(1,9).reshape((2,2,2))

In [3]: b=np.arange(1,5).reshape((2,2))

In [4]: a
Out[4]: 
array([[[1, 2],
        [3, 4]],

       [[5, 6],
        [7, 8]]])

In [5]: b
Out[5]: 
array([[1, 2],
       [3, 4]])

In [8]: np.tensordot(a,b,axes=1)
Out[8]: 
array([[[ 7, 10],
        [15, 22]],

       [[23, 34],
        [31, 46]]])

In [9]: np.tensordot(a,b,axes=(-1,0))
Out[9]: 
array([[[ 7, 10],
        [15, 22]],

       [[23, 34],
        [31, 46]]])

In [10]: np.tensordot(a,b,axes=2)
Out[10]: array([30, 70])

In [12]: np.tensordot(a,b,axes=([-2,-1],[0,1]))
Out[12]: array([30, 70])

In [13]: np.dot(a,b)
Out[13]: 
array([[[ 7, 10],
        [15, 22]],

       [[23, 34],
        [31, 46]]])

In [14]: np.tensordot(a,b,axes=(-1,-2))
Out[14]: 
array([[[ 7, 10],
        [15, 22]],

       [[23, 34],
        [31, 46]]])

In [15]: np.inner(a,b)
Out[15]: 
array([[[ 5, 11],
        [11, 25]],

       [[17, 39],
        [23, 53]]])

In [16]: np.tensordot(a,b,axes=(-1,-1))
Out[16]: 
array([[[ 5, 11],
        [11, 25]],

       [[17, 39],
        [23, 53]]])
复制代码

前两个参数 ab 为参与运算的两个张量。 第三个参数 axes 用于指定收缩的轴,完整格式形如 ([a_axis1,a_axis2,...],[b_axis1,b_axis2,...]) ,两个序列分别指定 ab 的轴,轴方向上的元素会按照被指定的顺序对应、相乘并相加;可使用 (a_axis,b_axis) 的形式仅指定一个轴; 可使用 int 类型快速指定 a 中最后N个轴和 b 中前N个轴用于收缩,即 0 等效于 ([],[]) ,对应张量积运算, 1 等效于 ([-1],[0]) ,对应张量内积运算, 2 等效于 ([-2,-1],[0,1]) ,对应张量双收缩运算, axes 的默认值为 2np.dotnp.inner 这类运算可视作该函数表示的几个特例, np.dot(a1,a2) 等效于 np.tensordot(a1,a2,axes=(-1,-2))np.inner(a1,a2) 等效于 np.tensordot(a1,a2,axes=(-1,-1))

3. 傅里叶变换

(1) 频率序列

In [12]: np.fft.fftfreq(10,1.0)
Out[12]: array([ 0. ,  0.1,  0.2,  0.3,  0.4, -0.5, -0.4, -0.3, -0.2, -0.1])
复制代码

返回离散傅里叶变换采样频率,第一个参数 n 为窗口长度, int 类型,第二个参数 d 为采样间距(采样率的倒数),默认为1.0,返回值单位与采样间距单位相对应。 返回值序列的计算方式: 如果 n 是偶数, f = [0, 1, ..., n/2-1, -n/2, ..., -1] / (d*n) ; 如果 n 是奇数, f = [0, 1, ..., (n-1)/2, -(n-1)/2, ..., -1] / (d*n)

(2) 快速傅里叶变换

In [13]: x=np.cos(np.linspace(0,2*np.pi,30))

In [14]: y=np.fft.fft(x)

In [15]: x2=np.fft.ifft(y)

In [16]: np.abs(x2-x).max()
Out[16]: 3.8864883384594504e-16
复制代码

np.fftnp.ifft 互为逆运算,用于一维快速傅里叶变换,经 np.fft 变换后的序列可通过 np.ifft 近似还原为原序列。 第二个参数 n 指定输出的变换轴长度,超长裁剪,不足补0; 第三个参数 axis 指定用于变换的轴,默认最后一个轴。

(3) 移频

In [21]: x=np.linspace(0,2*np.pi,8)

In [22]: y=np.fft.fft(x)

In [23]: y2=np.fft.fftshift(y)

In [24]: y3=np.fft.ifftshift(y2)

In [32]: np.abs(y3-y).max()
Out[32]: 0.0
复制代码

np.fftshiftnp.ifftshift 互为逆运算,用于将傅里叶变换输出中的直流分量移动到频谱的中央。第二个参数 axis 可指定用于转移的轴。

4. 图像处理

图像数据的存储方式是类似于数组的,可借助 PIL 库读取图片,再将图像数据转成 ndarray 进行计算处理。 以下提供一些使用 numpy 配合 PIL 处理图片数据的方法:

(1) 图片的创建、读取、缩放、保存

In [2]: from PIL import Image
   ...: image1= Image.open("D:\\test.jpg")

In [3]: image1.size
Out[3]: (1015, 610)

In [4]: image2=image1.resize((500,300))

In [5]: image2.save("D:\\test2.jpg")

In [6]: image0=Image.new('RGB',image2.size)

In [7]: image2
Out[7]: 
复制代码
[python][科学计算][numpy]使用指南

Image.open 用于打开一张图片, mode 可以设置读取模式,最常用的是 'L' 灰度图和 'RGB' 彩色图,一般会自动匹配不需要设置。 图片对象的 resize 方法可以缩放图片,大小参数以 tuple 类型 (width,height) 格式传入。 图片对象的 save 方法可以保存图片,通过保存路径中的文件拓展名或是 format 参数指定保存文件类型, quality 参数可以设置保存图像的质量,取值1(最差)~ 95(最佳),默认75。 Image.new 用于创建一个新的图片, mode 参数指定模式, size 指定大小。 在 IPython 中,直接输入图片变量名就可以显示图片。

(2) 图片与数组之间的转换

In [7]: a=np.asarray(image2)

In [8]: a.shape
Out[8]: (300, 500, 3)

In [9]: image3=Image.fromarray(a,'RGB')
复制代码

np.asarraynp.array 可以将图片转换为 ndarraynp.asarray 返回图片数据的 ndarray 类型视图,不能更改。根据类型的不同,得到的数组形状也不一样,常见的两种,灰度图转换得到形状为 (height,width) 的数组,彩色图转换得到形状为 (height,width,channel) 的数组, channel 即颜色通道, RGB 模式下 channel=3 ,分别对应红绿蓝。 Image.fromarray 可以将 ndarray 转换为图片,可通过 mode 参数指定模式,默认会根据数组形状自动匹配,当指定某个模式时,数组形状也需要匹配。 注意,数组转图片需要是 uint8 数据类型,取值 0~255 ,如不符合要进行另外的转换。

(3) 像素点绘制

In [10]: a0=np.array(image0)

In [11]: a0+=np.random.randint(0,256,a.shape,dtype='uint8')

In [12]: image4=Image.fromarray(a0)

In [13]: image4.putpixel((0,0),(255,255,255))

In [14]: image4 
Out[14]: 
复制代码
[python][科学计算][numpy]使用指南

图片对象的 putpixel 方法可以添加单个像素点,第一个参数 xy(x,y) 的形式声明添加像素点的位置,第二个参数 value 指定像素点的值,例如, L 灰度图模式下为标量值, RGB 彩色图模式下为 (r,g,b) 形式的 tuple ,取值均在 0~255 之间。该方法一次只能绘制一个像素点,效率低,在需要批量绘制时建议转换为 ndarray 处理。 上面的示例中使用 ndarray 的方法为新图片的每一个像素点添加了随机色彩。

(4) 灰度图和彩色图之间的转换

In [72]: a2=np.asarray(image2,dtype='float')

In [73]: a5=(11*a2[:,:,0]+16*a2[:,:,1]+5*a2[:,:,2])/32

In [74]: image5=Image.fromarray(np.uint8(a5))

In [75]: a6=a5.repeat(3).reshape(a5.shape+(3,))

In [76]: image6=Image.fromarray(np.uint8(a6))

In [77]: image5
Out[77]: 
复制代码
[python][科学计算][numpy]使用指南
彩色图转灰度图, L=11*R+16*G+5*B 只是一种可行的公式,也有其他公式可用。 灰度图转彩色图,较为简单,即将灰度值拷贝到 RGB

3个通道上,转换后颜色还是灰色,因为灰度图不具备色彩方面的信息,即使先将彩色图转灰度图,再转换回彩色图,色彩信息同样会丢失。

(5) 图片的翻转、旋转、裁剪

In [84]:a2=np.asarray(image2)

In [85]:Image.fromarray(a2[::-1,::-1,:])
Out[85]: 
复制代码
[python][科学计算][numpy]使用指南
In [90]:Image.fromarray(a2.transpose([1,0,2]))
Out[90]: 
复制代码
[python][科学计算][numpy]使用指南
In [91]:Image.fromarray(a2[:150,:,:])
Out[91]: 
复制代码
[python][科学计算][numpy]使用指南

以上示例展示了 左右翻转+上下翻转,左右翻转+逆时针转90°,截取上半部分 三种情况,借助倒序索引和轴交换的组合可以得到90°倍数旋转和上下左右翻转的所有组合情形。精确的旋转需要使用矩阵运算,此处不作展开。


以上就是本文的全部内容,希望本文的内容对大家的学习或者工作能带来一定的帮助,也希望大家多多支持 码农网

查看所有标签

猜你喜欢:

本站部分资源来源于网络,本站转载出于传递更多信息之目的,版权归原作者或者来源机构所有,如转载稿涉及版权问题,请联系我们

Dreamweaver基础教程

Dreamweaver基础教程

李振华、季小武、季小武、李振华 / 清华大学 / 2005-6 / 23.00元

本书通过实例的方式介绍了Macromedia公司的Dreamweaver MX 2004的使用方法和技巧。 全书由14章组成,第1章和第2章介绍了软件的应用领域、知识结构、界面组成等;第3章到第12章是本书的重点部分,通过实例制作介绍了站点的建立,表格、文本及样式的创建,链接、图像、行为的使用,层、表单、框架的创建和使用以及动画、多媒体的制作等主要知识点;第13章和第14章介绍了插件技......一起来看看 《Dreamweaver基础教程》 这本书的介绍吧!

XML、JSON 在线转换
XML、JSON 在线转换

在线XML、JSON转换工具

RGB HSV 转换
RGB HSV 转换

RGB HSV 互转工具