Py学习  »  Python

数据科学|可视化图解Python科学计算包Numpy

运筹OR帷幄 • 4 年前 • 539 次点击  


↑↑↑↑↑点击上方蓝色字关注我们!







『运筹OR帷幄』转载


作者:梁斌伟


编者按:

你真的会Numpy吗?本文从数组操作、多维处理、公式计算、数据展示三个角度详细介绍了科学计算包Numpy的使用方法,快来看看吧!



NumPy包是python生态系统中数据分析、机器学习和科学计算的主力。它极大地简化了向量和矩阵的操作。Python的一些主要软件包依赖于NumPy作为其基础架构的基础部分(例如scikit-learn、SciPy、pandas和tensorflow)。我们将介绍一些使用NumPy的主要方法,以及在我们为机器学习模型提供服务之前它如何表示不同类型的数据(表格、声音和图像)。

import numpy as np

一、数组操作

1.1创建数组

通过将python列表传递给NumPy对象,使用np.array()创建一个NumPy数组(即ndarray)。Python创建了我们在右边可以看到的数组:

一般情况,我们希望直接使用NumPy作为初始化的数组数据。NumPy为这些情况提供了诸如ones(),zeros()和random.random()类等方法。我们只需要向这些方法传递要生成的元素数量的参数:

一旦我们创建了数组,我们就可以开始以有趣的方式操作它们。

1.2 数组算术

先来创建两个NumPy数组来展示它们的用处。我们称之为data和ones:

这两个数组的加法就像我们+1那么简单(每一行相加)。

可以发现这样的计算不必在循环中实现。这是非常好的抽象处理,可以让你在更高层次上思考问题,而不是陷入如何实现计算的算法里面。

我们还可以像下面一样进行其他运算:

通常情况下,我们希望在数组和单个数字之间执行操作(我们也可以将其称为向量和标量之间的操作)。比如说,我们的数组表示以英里为单位的距离,我们希望将其转换为公里数。我们只需要让数据* 1.6就可以了。

NumPy让每个单元格都会发生相乘叫做广播。

1.3索引数组

我们可以索引、切片NumPy数组。

1.4 数组聚合

NumPy为我们提供非常好用的聚合功能:

除了最小值、最大值和总和之外,还得到非常棒的东西,比如平均值、所有元素相乘的结果、标准差,以及其他很多。

二、多维处理

2.1 创建矩阵

所有的例子都在一个维度上处理向量。NumPy的厉害之处是能够将我们目前所看到的所有内容应用到任意维度上。我们可以传递一个Python列表,让NumPy创建一个矩阵来表示它们:

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

我们也可以使用上面提到的相同方法(ones(),zeros()和random.random()),只要我们给它们一个元组来描述我们正在创建的矩阵的维度:

2.2 矩阵算术

如果两个矩阵的大小相同,我们可以使用算术运算符(+ - * /)来进行矩阵计算。NumPy将这些作为位置操作处理:

只有当不同维度为1时(例如矩阵只有一列或一行),我们才能在不同大小的矩阵上进行这些算术运算,在这种情况下,NumPy将其广播规则用于该操作:

2.3 点积

NumPy为每个矩阵提供了一个dot()方法,我们可以用它来执行与其他矩阵的点积运算:

两个矩阵在它们彼此面对的一侧必须具有相同的尺寸(上图底部红色的数字)。 
您可以将此操作可视化为如下所示:

2.4 矩阵索引

当我们操作矩阵时,索引和切片操作变得更加有用:

冒号表示从哪个位置到哪个位置,留空表示开头或者结尾;逗号表示行和列。NumPy的索引是从0开始的,并且后面的方括号是不包含后面的值(即小于后面的值)。

2.5 矩阵聚合

我们可以像聚合向量一样聚合矩阵:

我们不仅可以聚合矩阵中的所有值,还可以使用axis参数在行或列之间进行聚合:

2.6 转置和重塑

处理矩阵时的一个共同需求是需要旋转矩阵。当我们需要采用两个矩阵的点积并需要对齐它们共享的维度时,通常就是这种情况。NumPy数组有一个方便的属性叫做T来获得矩阵的转置:

在更高级的应用中,您可能会发现自己需要切换某个矩阵的维度。在机器学习应用程序中通常就是这种情况,其中某个模型期望输入的某个形状与您的数据集不同。在这些情况下,NumPy的reshape()法很有用。您只需将矩阵所需的新尺寸传递给它即可。

2.7 更多维度

NumPy可以完成我们在任何维度定义数据。其中心数据结构称为ndarray(N维数组)。

在很多方面,处理新维度只需在NumPy函数的参数中添加逗号:

三、公式计算

我们举一个计算均方误差的例子。实现适用于矩阵和向量的数学公式是考虑NumPy的关键用例。例如均方误差MSE公式,它是监督机器学习、模型处理、回归问题的核心公式:

在NumPy中实现这一点是轻而易举的:

这样做的好处是我们不需要关心predictions和labels是否包含一个或一百个值(只要它们的大小相同)。我们可以通过一个示例逐步执行该代码行中的四个操作:

predictions和labels向量都包含三个值。这意味着n的值为3。在我们执行减法后,我们最终得到如下值:

然后我们可以对矢量中的值进行平方:

然后求和:

结果就是均方误差。

四、数据表示

4.1 表格和电子表格

构建模型所需的所有数据类型(电子表格、图像和音频等),很多都非常适合在n维数组中表示。 
电子表格或值表是二维矩阵。电子表格中的每个工作表都可以是自己的变量。Python中最受欢迎的抽象是pandas数据帧DataFrame,它实际上使用NumPy并在其上构建。

4.2 音频

音频文件是一维样本数组。每个样本都是一个代表音频信号的一小部分的数字。CD质量的音频每秒可能有44,100个样本,每个样本是-65535到65536之间的整数。这意味着如果你有一个10秒的CD质量的WAVE文件,你可以将它加载到长度为10 * 44,100的NumPy数组中 = 441,000个样本。想要提取音频的第一秒?只需将文件加载到我们称之为audio的NumPy数组中,然后获取audio[:44100]。 
以下是一段音频文件:

4.3 图像

图像是高度×宽度的像素矩阵。

如果图像是黑白的(即灰度),则每个像素可以由单个数字表示(通常在0(黑色)和255(白色)之间)。想要裁剪图像的左上角10 x 10像素部分?告诉NumPy让你image[:10,:10]。 
这是一个图像文件的片段:

如果图像是彩色的,则每个像素由三个数字表示(RGB)。在这种情况下,我们需要第三维(因为每个单元格只能包含一个数字)。因此彩色图像由尺寸的ndarray表示:(高x宽x 3)。


『运筹OR帷幄算法社区』知识星球,依托社区25w+专业受众和20+细分领域硕博微信群,特邀腾讯,百度,阿里,华为等公司大咖与您一起聊算法,快来扫码加入!

加入『运筹OR帷幄』知识星球的好处

● 全球Top名校教授|博士和名企研发高管一起交流算法相关学术|研发干货

● 中国你能说出名字的几乎所有大厂|欧美数家大厂(资深)算法工程师入驻

● 依托『运筹OR帷幄』25w+专业受众和20+细分领域硕博微信群的算法技术交流

● 以上所有公司|高校独家内推招聘|实习机会、多家Offer选择指导

● 以面试题|作业题|业界项目为学习资料学习算法干货,从小白变成大咖

● 不定期的线上、线下交流会和聚会,拓展人脉


相关文章推荐

普林斯顿博士后的彪悍项目,用Numpy手写3万行代码、30多个模型,这也许能打造最强机器学习基石?

点击蓝字标题,即可阅读《AI| 『运筹OR帷幄』惊为天人,Numpy手写全部主流机器学习模型,代码超过3万行》


除了Numpy,Pandas也是一个数据处理的高效工具!

点击蓝字标题,即可阅读《数据科学| 『运筹OR帷幄』Pandas数据分析入门》



温馨提示

可以在 公众号后台 回复关键词:“学界”获取大量由我平台编辑精心整理的学习资料,如果觉得有用, 请勿吝啬你的留言和赞哦!~


—— 完 ——




文章申明

Sept. 2019

文章作者:梁斌伟

责任编辑:郭德真

微信编辑:玖蓁

本文转载自公众号 实战统计学

(ID:statrcn)

原文链接:

可视化图解Python科学计算包NumPy




优质公众号推荐

击查看详情

   


Python社区是高质量的Python/Django开发社区
本文地址:http://www.python88.com/topic/40293
 
539 次点击