Numpy怎么使用最高效,5个技巧送给你!

开发 后端
本文作者将分享 5 个优雅的 Python Numpy 函数,有助于高效、简洁的数据处理。

在 reshape 函数中使用参数-1

Numpy 允许我们根据给定的新形状重塑矩阵,新形状应该和原形状兼容。有意思的是,我们可以将新形状中的一个参数赋值为-1。这仅仅表明它是一个未知的维度,我们希望 Numpy 来算出这个未知的维度应该是多少:Numpy 将通过查看数组的长度和剩余维度来确保它满足上述标准。让我们来看以下例子: 

维度为-1 的不同 reshape 操作图示。 

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

假设我们给定行参数为 1,列参数为-1,那么 Numpy 将计算出 reshape 后的列数为 8。 

  1. a.reshape(1,-1)  
  2. array([[1, 2, 3, 4, 5, 6, 7, 8]]) 

假设我们给定行参数为-1,列参数为 1,那么 Numpy 将计算出 reshape 后的行数为 8。 

  1. a.reshape(-1,1)  
  2. array([[1],  
  3.        [2],  
  4.        [3],  
  5.        [4],  
  6.        [5],  
  7.        [6],  
  8.        [7],  
  9.        [8]]) 

下面的代码也是一样的道理。 

  1. a.reshape(-1,4)  
  2. array([[1, 2, 3, 4],  
  3.        [5, 6, 7, 8]])a.reshape(-1,2)  
  4. array([[1, 2],  
  5.        [3, 4],  
  6.        [5, 6],  
  7.        [7, 8]])a.reshape(2,-1)  
  8. array([[1, 2, 3, 4],  
  9.        [5, 6, 7, 8]])a.reshape(4,-1)  
  10. array([[1, 2],  
  11.        [3, 4],  
  12.        [5, 6],  
  13.        [7, 8]]) 

这也适用于任何更高维度张量的 reshape,但是只有一个维度的参数能赋值为-1。 

  1. a.reshape(2,2,-1)  
  2. array([[[1, 2],  
  3.         [3, 4]],  
  4.        [[5, 6],  
  5.         [7, 8]]])a.reshape(2,-1,1)  
  6. array([[[1],  
  7.         [2],  
  8.         [3],  
  9.         [4]],   
  10.        [[5],  
  11.         [6],  
  12.         [7],  
  13.         [8]]]) 

如果我们尝试 reshape 不兼容的形状或者是给定的未知维度参数多于 1 个,那么将会报错。 

  1. a.reshape(-1,-1)  
  2. ValueError: can only specify one unknown dimensiona.reshape(3,-1)  
  3. ValueError: cannot reshape array of size 8 into shape (3,newaxis) 

总而言之,当试图对一个张量进行 reshape 操作时,新的形状必须包含与旧的形状相同数量的元素,这意味着两个形状的维度乘积必须相等。当使用 -1 参数时,与-1 相对应的维数将是原始数组的维数除以新形状中已给出维数的乘积,以便维持相同数量的元素。

Argpartition:在数组中找到最大的 N 个元素。

Numpy 的 argpartion 函数可以高效地找到 N 个最大值的索引并返回 N 个值。在给出索引后,我们可以根据需要进行值排序。 

  1. array = np.array([10, 7, 4, 3, 2, 2, 5, 9, 0, 4, 6, 0])index = np.argpartition*(array, -5)[-5:]  
  2. index  
  3. array([ 6,  1, 10,  7,  0], dtype=int64)np.sort(array[index])  
  4. array([ 5,  6,  7,  9, 10]) 

Clip:如何使数组中的值保持在一定区间内

在很多数据处理和算法中(比如强化学习中的 PPO),我们需要使得所有的值保持在一个上下限区间内。Numpy 内置的 Clip 函数可以解决这个问题。Numpy clip () 函数用于对数组中的值进行限制。给定一个区间范围,区间范围外的值将被截断到区间的边界上。例如,如果指定的区间是 [-1,1],小于-1 的值将变为-1,而大于 1 的值将变为 1。

Clip 示例:限制数组中的最小值为 2,最大值为 6。 

  1. #Example-1  
  2. array = np.array([10, 7, 4, 3, 2, 2, 5, 9, 0, 4, 6, 0])  
  3. print (np.clip(array,2,6))[6 6 4 3 2 2 5 6 2 4 6 2]#Example-2  
  4. array = np.array([10, -1, 4, -3, 2, 2, 5, 9, 0, 4, 6, 0])  
  5. print (np.clip(array,2,5))[5 2 4 2 2 2 5 5 2 4 5 2] 

Extract:从数组中提取符合条件的元素

我们可以使用 Numpy extract () 函数从数组中提取符合条件的特定元素。

 

  1. arr = np.arange(10)  
  2. arrarray([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])# Define the codition, here we take MOD 3 if zero  
  3. condition = np.mod(arr, 3)==0  
  4. conditionarray([ True, False, False,  True, False, False,  True, False, False,True])np.extract(condition, arr) 
  5. array([0, 3, 6, 9]) 

同样地,如果有需要,我们可以用 AND 和 OR 组合的直接条件,如下所示:

  1. np.extract(((arr > 2) & (arr < 8)), arr)array([3, 4, 5, 6, 7]) 

setdiff1d:如何找到仅在 A 数组中有而 B 数组没有的元素

返回数组中不在另一个数组中的独有元素。这等价于两个数组元素集合的差集。

 

  1. a = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9])  
  2. b = np.array([3,4,7,6,7,8,11,12,14])  
  3. c = np.setdiff1d(a,b)  
  4. carray([1, 2, 5, 9]) 

小结

以上 5 个 Numpy 函数并不经常被社区使用,但是它们非常简洁和优雅。在我看来,我们应该尽可能在出现类似情况时使用这些函数,不仅因为代码量更少,更因为它们是解决复杂问题的绝妙方法。 

 

责任编辑:庞桂玉 来源: 马哥Linux运维
相关推荐

2018-05-04 09:14:09

Git技巧shell命令

2018-06-16 08:35:57

UnixLinux命令

2021-09-13 07:46:06

Kubectl Kubernetes 工具

2021-07-07 09:50:23

NumpyPandasPython

2021-09-12 17:23:57

canvas动画函数

2009-10-27 09:09:06

Eclipse技巧

2017-06-05 12:06:00

2019-02-11 11:16:13

2021-02-14 19:24:45

SpringRegistrar对象

2023-03-29 07:45:58

VS编辑区编程工具

2018-06-20 11:00:06

云应用开发PaaS

2021-12-09 23:24:56

勒索软件攻击网络安全

2019-03-19 13:44:41

Python编程技巧编程语言

2022-02-24 10:05:20

Python编程语言代码

2017-10-23 16:28:33

2019-08-14 09:43:12

开发技能代码

2011-09-25 10:46:18

云计算安全

2024-11-07 15:55:22

PyTorchNumPyPython

2021-02-19 10:59:29

NumpyPandasPython

2022-05-04 12:44:57

Python编程语言
点赞
收藏

51CTO技术栈公众号