数据可视化 指的是通过可视化表示来探索数据,它与数据挖掘 紧密相关,而数据挖掘指的是使用代码来探索数据集的规律和关联。数据集可以是用一行代码就能表示的小型数字列表,也可以是数以吉字节的数据。
漂亮地呈现数据关乎的并非仅仅是漂亮的图片。以引人注目的简洁方式呈现数据,让观看者能够明白其含义,发现数据集中原本未意识到的规律和意义。

所幸即便没有超级计算机,也能够可视化复杂的数据。鉴于Python的高效性,使用它在笔记本电脑上就能快速地探索由数百万个数据点组成的数据集。数据点并非必须 是数字,利用本书前半部分介绍的基本知识,也可以对非数字数据进行分析。

在基因研究、天气研究、政治经济分析等众多领域,大家都使用Python来完成数据密集型工作。数据科学家使用Python编写了一系列令人印象深刻的可视化和分析工 具,其中很多也可供你使用。最流行的工具之一是matplotlib,它是一个数学绘图库,我们将使用它来制作简单的图表,如折线图和散点图。然后,我们将基于随机漫 步概念生成一个更有趣的数据集——根据一系列随机决策生成的图表。
我们还将使用Pygal包,它专注于生成适合在数字设备上显示的图表。通过使用Pygal,可在用户与图表交互时突出元素以及调整其大小,还可轻松地调整整个图表的尺 寸,使其适合在微型智能手表或巨型显示器上显示。我们将使用Pygal以各种方式探索掷骰子的结果。

15.1安装matplotlib
首先,需要安装matplotlib,我们将使用它来制作开始的几个图表。如果你还未使用过pip,请参阅12.2.1节。

15.1.1在Linux系统中安装matplotlib
如果你使用的是系统自带的Python版本,可使用系统的包管理器来安装matplotlib,为此只需执行一行命令:

$  sudo  apt-get  install  python3-matplotlib

如果你使用的是Python 2.7,请执行如下命令:

$  sudo  apt-get  install  python-matplotlib

如果你安装了较新的Python版本,就必须安装matplotlib依赖的一些库:

$  sudo  apt-get  install  python3.5-dev  python3.5-tk  tk-dev
$  sudo  apt-get  install  libfreetype6-dev  g++

再使用pip来安装matplotlib:

15.1.2在OS X系统中安装matplotlib
Apple的标准Python安装自带了matplotlib。要检查系统是否安装了matplotlib,可打开一个终端会话并尝试导入matplotlib。如果系统没有自带matplotlib,且你的Python是使用Homebrew安装的,则可以像下面这样安装matplotlib:

$  pip  install  --user  matplotlib

注意 安装包时可能需要使用pip3 ,而不是pip 。另外,如果这个命令不管用,你可能需要删除标志–user 。

15.1.3在Windows系统中安装matplotlib
在Windows系统中,首先需要安装Visual Studio。为此,请访问https://dev.windows.com/ ,单击Downloads,再查找Visual Studio Community——一组免费的Windows开发工具。请下载并运行该安装程序。

接下来,需要下载matplotlib安装程序。为此,请访问https://pypi.python.org/pypi/matplotlib/ ,并查找与你使用的Python版本匹配的wheel文件(扩展名为.whl的文件)。例如,如果你使用的是32位的Python 3.5,则需要下载matplotlib-1.4.3-cp35-none-win32.whl。
注意 如果找不到与你安装的Python版本匹配的文件,请去http://www.lfd.uci.edu/-gohlke/pythonlibs/#matplotlib 看看,这个网站发布安装程序的时间通常比matplotlib官网早些。
将这个.whl文件复制到你的项目文件夹,打开一个命令窗口,并切换到该项目文件夹,再使用pip来安装matplotlib:

>  cd  python_work
python_work>  python  -m  pip  install  --user  matplotlib-1.4.3-cp35-none-win32.whl

15.1.4测试matplotlib
安装必要的包后,对安装进行测试。为此,首先使用命令python 或python3 启动一个终端会话,再尝试导入matplotlib:

$python3
>>>  import  matplotlib
>>>

如果没有出现任何错误消息,就说明你的系统安装了matplotlib,可以接着阅读下一节。
注意 如果你在安装过程中遇到了麻烦,请参阅附录C。如果依然无济于事,请向他人寻求帮助。对于你遇到的问题,只要向经验丰富的Python程序员提供少量的信息,他们很可能很快就能帮你解决。

15.1.5matplotlib画廊
要查看使用matplotlib可制作的各种图表,请访问http://matplotlib.org/ 的示例画廊。单击画廊中的图表,就可查看用于生成图表的代码。

15.2绘制简单的折线图
下面来使用matplotlib绘制一个简单的折线图,再对其进行定制,以实现信息更丰富的数据可视化。我们将使用平方数序列1、4、9、16和25来绘制这个图表。 只需向matplotlib提供如下数字,matplotlib就能完成其他的工作:
mpl_squares.py

import  matplotlib.pyplot  as  plt
squares  =  [1,  4,  9,  16,  25]
plt.plot(squares)
plt.show()

我们首先导入了模块pyplot ,并给它指定了别名plt ,以免反复输入pyplot 。在线示例大都这样做,因此这里也这样做。模块pyplot 包含很多用于生成图表的函数。
我们创建了一个列表,在其中存储了前述平方数,再将这个列表传递给函数plot() ,这个函数尝试根据这些数字绘制出有意义的图形。plt.show() 打开matplotlib查看器,并显示绘制的图形,如图15-1所示。查看器让你能够缩放和导航图形,另外,单击磁盘图标可将图形保存起来。

图15-1 使用matplotlib可制作的最简单的图表

15.2.1修改标签文字和线条粗细
图15-1所示的图形表明数字是越来越大的,但标签文字太小,线条太细。所幸matplotlib让你能够调整可视化的各个方面。

下面通过一些定制来改善这个图形的可读性,如下所示:
mpl_squares.py

import  matplotlib.pyplot  as  plt
squares  =  [1,  4,  9,  16,  25]
plt.plot(squares,  linewidth=5)# 设置图表标题,并给坐标轴加上标签
plt.title("Square  Numbers",  fontsize=24)
plt.xlabel("Value",  fontsize=14)
plt.ylabel("Square  of  Value",  fontsize=14)
# 设置刻度标记的大小
plt.tick_params(axis='both',  labelsize=14)
plt.show()

参数linewidth (见❶)决定了plot() 绘制的线条的粗细。函数title() (见❷)给图表指定标题。在上述代码中,出现了多次的参数fontsize 指定了图表中文字的大小。
函数xlabel() 和ylabel() 让你能够为每条轴设置标题(见❸);而函数tick_params() 设置刻度的样式(见❹),其中指定的实参将影响x 轴和y 轴上的刻度
(axes=‘both’ ),并将刻度标记的字号设置为14(labelsize=14 )。
最终的图表阅读起来容易得多了,如图15-2所示:标签文字更大,线条也更粗。

            图15-2 现在图表阅读起来容易得多

15.2.2校正图形
图形更容易阅读后,我们发现没有正确地绘制数据:折线图的终点指出4.0的平方为25!下面来修复这个问题。
当你向plot() 提供一系列数字时,它假设第一个数据点对应的 x 坐标值为0,但我们的第一个点对应的 x 值为1。为改变这种默认行为,我们可以给plot() 同时提供输入值和输出值:
mpl_squares.py

import  matplotlib.pyplot  as  plt
input_values  =  [1,  2,  3,  4,  5]
squares  =  [1,  4,  9,  16,  25]
plt.plot(input_values,  squares,  linewidth=5)# 设置图表标题并给坐标轴加上标签
plt.title("Square  Numbers",  fontsize=24)
plt.xlabel("Value",  fontsize=14)
plt.ylabel("Square  of  Value",  fontsize=14)
# 设置刻度标记的大小
plt.tick_params(axis='both',  labelsize=14)
plt.show()

现在plot() 将正确地绘制数据,因为我们同时提供了输入值和输出值,它无需对输出值的生成方式作出假设。最终的图形是正确的,如图15-3所示。

图15-3 根据数据正确地绘制了图形

使用plot() 时可指定各种实参,还可使用众多函数对图形进行定制。本章后面处理更有趣的数据集时,将继续探索这些定制函数。

15.2.3使用scatter() 绘制散点图并设置其样式
有时候,需要绘制散点图并设置各个数据点的样式。例如,你可能想以一种颜色显示较小的值,而用另一种颜色显示较大的值。绘制大型数据集时,你还可以对每个点都设置同样的样式,再使用不同的样式选项重新绘制某些点,以突出它们。
要绘制单个点,可使用函数scatter() ,并向它传递一对 x 和 y 坐标,它将在指定位置绘制一个点:
scatter_squares.py

import  matplotlib.pyplot  as  pltplt.scatter(2,  4)
plt.show()

下面来设置输出的样式,使其更有趣:添加标题,给轴加上标签,并确保所有文本都大到能够看清:

import  matplotlib.pyplot  as  plt
plt.scatter(2,  4,  s=200)# 设置图表标题并给坐标轴加上标签plt.title("Square  Numbers",  fontsize=24) plt.xlabel("Value",  fontsize=14) plt.ylabel("Square  of  Value",  fontsize=14)# 设置刻度标记的大小
plt.tick_params(axis='both',  which='major',  labelsize=14)
plt.show()

在❶处,我们调用了scatter() ,并使用实参s 设置了绘制图形时使用的点的尺寸。如果此时运行scatter_squares.py,将在图表中央看到一个点,如图15-4所示。

图15-4 绘制单个点

15.2.4使用scatter() 绘制一系列点
要绘制一系列的点,可向scatter() 传递两个分别包含x 值和y 值的列表,如下所示:
scatter_squares.py

import  matplotlib.pyplot  as  plt
x_values  =  [1,  2,  3,  4,  5]
y_values  =  [1,  4,  9,  16,  25]
plt.scatter(x_values,  y_values,  s=100)
# 设置图表标题并给坐标轴指定标签
plt.title("Square  Numbers",  fontsize=24)
plt.xlabel("Value",  fontsize=14)
plt.ylabel("Square  of  Value",  fontsize=14)# 设置刻度标记的大小
plt.tick_params(axis='both',  which='major',  labelsize=14)
plt.show()

列表x_values 包含要计算其平方值的数字,而列表y_values 包含前述每个数字的平方值。将这些列表传递给scatter() 时,matplotlib依次从每个列表中读取一个值来绘制一个点。要绘制的点的坐标分别为 (1, 1)、(2, 4)、(3, 9)、(4, 16)和(5, 25),最终的结果如图15-5所示。

图15-5 由多个点组成的散点图

15.2.5自动计算数据
手工计算列表要包含的值可能效率低下,需要绘制的点很多时尤其如此。可以不必手工计算包含点坐标的列表,而让Python循环来替我们完成这种计算。下面是绘制1000个点的代 码:
scatter_squares.py


import  matplotlib.pyplot  as  pltx_values  =  list(range(1,  1001))
y_values  =  [x**2  for  x  in  x_values]
plt.scatter(x_values,  y_values,  s=40)
# 设置图表标题并给坐标轴加上标签
plt.title("Square  Numbers",  fontsize=24)
plt.xlabel("Value",  fontsize=14)
plt.ylabel("Square  of  Value",  fontsize=14)
# 设置每个坐标轴的取值范围
plt.axis([0,  1100,  0,  1100000])
plt.show()

我们首先创建了一个包含 x 值的列表,其中包含数字1~1000(见❶)。接下来是一个生成 y 值的列表解析,它遍历 x 值(for x in x_values ),计算其平方值(x**2 ),并将结果存储到列表y_values 中。然后,将输入列表和输出列表传递给scatter() (见❷)。

由于这个数据集较大,我们将点设置得较小,并使用函数axis() 指定了每个坐标轴的取值范围(见❸)。函数axis() 要求提供四个值:x 和 y 坐标轴的最小值和最大值。在这里,我们将 x 坐标轴的取值范围设置为0~1100,并将 y 坐标轴的取值范围设置为0~1 100 000。结果如图15-6所示。

图15-6 Python绘制1000个点与绘制5个点一样容易

15.2.6删除数据点的轮廓
matplotlib允许你给散点图中的各个点指定颜色。默认为蓝色点和黑色轮廓,在散点图包含的数据点不多时效果很好。但绘制很多点时,黑色轮廓可能会粘连在一起。要删除数据点的轮廓,可在调用scatter() 时传递实参edgecolor=‘none’ :

plt.scatter(x_values,  y_values,  edgecolor='none',  s=40)

将相应调用修改为上述代码后,如果再运行scatter_squares.py,在图表中看到的将是蓝色实心点。

15.2.7自定义颜色
要修改数据点的颜色,可向scatter() 传递参数c ,并将其设置为要使用的颜色的名称,如下所示:

plt.scatter(x_values,  y_values,  c='red',  edgecolor='none',  s=40)

你还可以使用RGB颜色模式自定义颜色。要指定自定义颜色,可传递参数c ,并将其设置为一个元组,其中包含三个0~1之间的小数值,它们分别表示红色、绿色和蓝色分量。例 如,下面的代码行创建一个由淡蓝色点组成的散点图:

plt.scatter(x_values,  y_values,  c=(0,  0,  0.8),  edgecolor='none',  s=40)

值越接近0,指定的颜色越深,值越接近1,指定的颜色越浅。

15.2.8使用颜色映射
颜色映射(colormap)是一系列颜色,它们从起始颜色渐变到结束颜色。在可视化中,颜色映射用于突出数据的规律,例如,你可能用较浅的颜色来显示较小的值,并使用较深的颜色来显示较大的值。
模块pyplot 内置了一组颜色映射。要使用这些颜色映射,你需要告诉pyplot 该如何设置数据集中每个点的颜色。下面演示了如何根据每个点的 y 值来设置其颜色:
scatter_squares.py

import  matplotlib.pyplot  as  pltx_values  =  list(range(1001))
y_values  =  [x**2  for  x  in  x_values]plt.scatter(x_values,  y_values,  c=y_values,  cmap=plt.cm.Blues, edgecolor='none',  s=40)# 设置图表标题并给坐标轴加上标签
plt.title("Square  Numbers",  fontsize=24)
plt.xlabel("Value",  fontsize=14)
plt.ylabel("Square  of  Value",  fontsize=14)
# 设置每个坐标轴的取值范围
plt.axis([0,  1100,  0,  1100000])
plt.show()


图15-7 使用颜色映射Blues 的图表

注意 要了解pyplot 中所有的颜色映射,请访问http://matplotlib.org/ ,单击Examples,向下滚动到Color Examples,再单击colormaps_reference。

15.2.9自动保存图表
要让程序自动将图表保存到文件中,可将对plt.show() 的调用替换为对plt.savefig() 的调用:

plt.savefig('squares_plot.png',  bbox_inches='tight')

第一个实参指定要以什么样的文件名保存图表,这个文件将存储到scatter_squares.py所在的目录中;第二个实参指定将图表多余的空白区域裁剪掉。如果要保留图表周围多余的空 白区域,可省略这个实参。

15.3随机漫步
在本节中,我们将使用Python来生成随机漫步数据,再使用matplotlib以引人瞩目的方式将这些数据呈现出来。随机漫步 是这样行走得到的路径:每次行走都完全是随机的,没有明确的方向,结果是由一系列随机决策决定的。你可以这样认为,随机漫步就是蚂蚁在晕头转向的情况下,每次都沿随机的方向前行所经过的路径。

在自然界、物理学、生物学、化学和经济领域,随机漫步都有其实际用途。例如,漂浮在水滴上的花粉因不断受到水分子的挤压而在水面上移动。水滴中的分子运动是随机的, 因此花粉在水面上的运动路径犹如随机漫步。我们稍后将编写的代码模拟了现实世界的很多情形。

15.3.1创建RandomWalk() 类
为模拟随机漫步,我们将创建一个名为RandomWalk 的类,它随机地选择前进方向。这个类需要三个属性,其中一个是存储随机漫步次数的变量,其他两个是列表,分别存储随机漫步经过的每个点的 x 和 y 坐标。
RandomWalk 类只包含两个方法:init() 和fill_walk() ,其中后者计算随机漫步经过的所有点。下面先来看看__init__() ,如下所示:
random_walk.py

from  random  import  choice
class  RandomWalk():"""一个生成随机漫步数据的类"""
def  __init__(self,  num_points=5000):""" 初 始 化 随 机 漫 步 的 属 性 """self.num_points  =  num_points#  所有随机漫步都始于(0,  0)self.x_values = [0]self.y_values = [0]

为做出随机决策,我们将所有可能的选择都存储在一个列表中,并在每次做决策时都使用choice() 来决定使用哪种选择(见❶)。接下来,我们将随机漫步包含的默认点数设置为5000,这大到足以生成有趣的模式,同时又足够小,可确保能够快速地模拟随机漫步(见❷)。然后,在❸处,我们创建了两个用于存储x 和y 值的列表,并让每次漫步都从点(0, 0)出发。

15.3.2选择方向
我们将使用fill_walk() 来生成漫步包含的点,并决定每次漫步的方向,如下所示。请将这个方法添加到random_walk.py中:
random_walk.py

from  random  import  choice
class  RandomWalk():"""一个生成随机漫步数据的类"""def  __init__(self,  num_points=5000):""" 初 始 化 随 机 漫 步 的 属 性 """self.num_points  =  num_points#  所有随机漫步都始于(0,  0)self.x_values = [0]self.y_values = [0]def  fill_walk(self):"""计算随机漫步包含的所有点"""# 不断漫步,直到列表达到指定的长度while  len(self.x_values)  <  self.num_points:# 决定前进方向以及沿这个方向前进的距离x_direction  =  choice([1,  -1])x_distance  =  choice([0,  1,  2,  3,  4])x_step  =  x_direction  *  x_distancey_direction = choice([1, -1])y_distance = choice([0, 1, 2, 3, 4])y_step = y_direction * y_distance# 拒绝原地踏步if x_step == 0 and y_step == 0:continue# 计算下一个点的x和y值next_x = self.x_values[-1] + x_stepnext_y = self.y_values[-1] + y_stepself.x_values.append(next_x)self.y_values.append(next_y)

在❶处,我们建立了一个循环,这个循环不断运行,直到漫步包含所需数量的点。这个方法的主要部分告诉Python如何模拟四种漫步决定:向右走还是向左走?沿指定的方向走多 远?向上走还是向下走?沿选定的方向走多远?

我们使用choice([1, -1]) 给x_direction 选择一个值,结果要么是表示向右走的1,要么是表示向左走的-1(见❷)。接下来,choice([0, 1, 2, 3, 4]) 随机地选择一个0~4之间的整数,告诉Python 沿指定的方向走多远(x_distance )。(通过包含0,我们不仅能够沿两个轴移动,还能够沿y 轴移动。)
在❸和❹处,我们将移动方向乘以移动距离,以确定沿 x 和 y 轴移动的距离。如果x_step 为正,将向右移动,为负将向左移动,而为零将垂直移动;如果y_step 为正,就意味着向上移动,为负意味着向下移动,而为零意味着水平移动。如果x_step 和y_step 都为零,则意味着原地踏步,我们拒绝这样的情况,接着执行下一次循环(见❺)。

为获取漫步中下一个点的 x 值,我们将x_step 与x_values 中的最后一个值相加(见❻),对于 y 值也做相同的处理。获得下一个点的 x 值和 y 值后,我们将它们分别附加到列表x_values 和y_values 的末尾。

15.3.3绘制随机漫步图
下面的代码将随机漫步的所有点都绘制出来:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk#  创建一个RandomWalk实例,并将其包含的点都绘制出来
rw  =  RandomWalk()
rw.fill_walk()
plt.scatter(rw.x_values,  rw.y_values,  s=15)
plt.show()

我们首先导入了模块pyplot 和RandomWalk 类,然后创建了一个RandomWalk 实例,并将其存储到rw 中(见❶),再调用fill_walk() 。在❷处,我们将随机漫步包含的 x 和 y 值传递给scatter() ,并选择了合适的点尺寸。图15-8显示了包含5000个点的随机漫步图(本节的示意图未包含matplotlib查看器部分,但你运行rw_visual.py时,依然会看到)。

图15-8 包含5000个点的随机漫步

15.3.4模拟多次随机漫步
每次随机漫步都不同,因此探索可能生成的各种模式很有趣。要在不多次运行程序的情况下使用前面的代码模拟多次随机漫步,一种办法是将这些代码放在一个while 循环中,如下所示:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
while  True:
#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk()rw.fill_walk()plt.scatter(rw.x_values,  rw.y_values,  s=15)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if  keep_running  ==  'n':break


这些代码模拟一次随机漫步,在matplotlib查看器中显示结果,再在不关闭查看器的情况下暂停。如果你关闭查看器,程序将询问你是否要再模拟一次随机漫步。如果你输入y ,可模拟多次随机漫步:这些随机漫步都在起点附近进行,大多沿特定方向偏离起点,漫步点分布不均匀等。要结束程序,请输入n 。
注意 如果你使用的是Python 2.7,别忘了将❶处的input() 替换为raw_input() 。

15.3.5设置随机漫步图的样式
在本节中,我们将定制图表,以突出每次漫步的重要特征,并让分散注意力的元素不那么显眼。为此,我们确定要突出的元素,如漫步的起点、终点和经过的路径。接下来确定要使其不那么显眼的元素,如刻度标记和标签。最终的结果是简单的可视化表示,清楚地指出了每次漫步经过的路径。

15.3.6给点着色
我们将使用颜色映射来指出漫步中各点的先后顺序,并删除每个点的黑色轮廓,让它们的颜色更明显。为根据漫步中各点的先后顺序进行着色,我们传递参数c ,并将其设置为一个列表,其中包含各点的先后顺序。由于这些点是按顺序绘制的,因此给参数c 指定的列表只需包含数字1~5000,如下所示:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
# while  True:
# #  创建一个RandomWalk实例,并将其包含的点都绘制出来
#     rw  =  RandomWalk()
#     rw.fill_walk()
#     plt.scatter(rw.x_values,  rw.y_values,  s=15)
#     plt.show()
#
#     keep_running  =  input("Make  another  walk?  (y/n):  ")
#     if  keep_running  ==  'n':
#         breakwhile  True:#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk()rw.fill_walk()point_numbers  =  list(range(rw.num_points))plt.scatter(rw.x_values,  rw.y_values,  c=point_numbers,  cmap=plt.cm.Blues, edgecolor='none',  s=15)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if keep_running == 'n':break

在❶处,我们使用了range() 生成了一个数字列表,其中包含的数字个数与漫步包含的点数相同。接下来,我们将这个列表存储在point_numbers 中,以便后面使用它来设 置每个漫步点的颜色。我们将参数c 设置为point_numbers ,指定使用颜色映射Blues ,并传递实参edgecolor=none 以删除每个点周围的轮廓。最终的随机漫步图从浅蓝色渐变为深蓝色,如图15-9所示。

图15-9 使用颜色映射Blues 着色的随机漫步图

15.3.7重新绘制起点和终点
除了给随机漫步的各个点着色,以指出它们的先后顺序外,如果还能呈现随机漫步的起点和终点就更好了。为此,可在绘制随机漫步图后重新绘制起点和终点。我们让起点和终点变得更大,并显示为不同的颜色,以突出它们,如下所示:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
# while  True:
# #  创建一个RandomWalk实例,并将其包含的点都绘制出来
#     rw  =  RandomWalk()
#     rw.fill_walk()
#     plt.scatter(rw.x_values,  rw.y_values,  s=15)
#     plt.show()
#
#     keep_running  =  input("Make  another  walk?  (y/n):  ")
#     if  keep_running  ==  'n':
#         breakwhile  True:#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk()rw.fill_walk()point_numbers  =  list(range(rw.num_points))plt.scatter(rw.x_values,  rw.y_values,  c=point_numbers,  cmap=plt.cm.Blues, edgecolor='none',  s=15)# 突出起点和终点plt.scatter(0, 0, c='green', edgecolors='none', s=100)plt.scatter(rw.x_values[-1], rw.y_values[-1], c='red', edgecolors='none', s=100)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if keep_running == 'n':break


为突出起点,我们使用绿色绘制点(0, 0),并使其比其他点大(s=100 )。为突出终点,我们在漫步包含的最后一个 x 和 y 值处绘制一个点,将其颜色设置为红色,并将尺寸设置为100。请务必将这些代码放在调用plt.show() 的代码前面,确保在其他点的上面绘制起点和终点。
如果你现在运行这些代码,将能准确地知道每次随机漫步的起点和终点(如果起点和终点不明显,请调整它们的颜色和大小,直到明显为止)。

15.3.8隐藏坐标轴
下面来隐藏这个图表中的坐标轴,以免我们注意的是坐标轴而不是随机漫步路径。要隐藏坐标轴,可使用如下代码:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
# while  True:
# #  创建一个RandomWalk实例,并将其包含的点都绘制出来
#     rw  =  RandomWalk()
#     rw.fill_walk()
#     plt.scatter(rw.x_values,  rw.y_values,  s=15)
#     plt.show()
#
#     keep_running  =  input("Make  another  walk?  (y/n):  ")
#     if  keep_running  ==  'n':
#         breakwhile  True:#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk()rw.fill_walk()point_numbers  =  list(range(rw.num_points))plt.scatter(rw.x_values,  rw.y_values,  c=point_numbers,  cmap=plt.cm.Blues, edgecolor='none',  s=15)# 突出起点和终点plt.scatter(0, 0, c='green', edgecolors='none', s=100)plt.scatter(rw.x_values[-1], rw.y_values[-1], c='red', edgecolors='none', s=100)# 隐藏坐标轴plt.axes().get_xaxis().set_visible(False)plt.axes().get_yaxis().set_visible(False)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if keep_running == 'n':break


为修改坐标轴,使用了函数plt.axes() (见❶)来将每条坐标轴的可见性都设置为False 。随着你越来越多地进行数据可视化,经常会看到这种串接方法的方式。如果你现在运行rw_visual.py,将看到一系列图形,但看不到坐标轴。

15.3.9增加点数
下面来增加点数,以提供更多的数据。为此,我们在创建RandomWalk 实例时增大num_points 的值,并在绘图时调整每个点的大小,如下所示:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
# while  True:
# #  创建一个RandomWalk实例,并将其包含的点都绘制出来
#     rw  =  RandomWalk()
#     rw.fill_walk()
#     plt.scatter(rw.x_values,  rw.y_values,  s=15)
#     plt.show()
#
#     keep_running  =  input("Make  another  walk?  (y/n):  ")
#     if  keep_running  ==  'n':
#         breakwhile  True:#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk(50000)rw.fill_walk()# 绘制点并将图形显示出来point_numbers = list(range(rw.num_points))plt.scatter(rw.x_values, rw.y_values, c=point_numbers, cmap=plt.cm.Blues, edgecolor='none', s=1)point_numbers  =  list(range(rw.num_points))plt.scatter(rw.x_values,  rw.y_values,  c=point_numbers,  cmap=plt.cm.Blues, edgecolor='none',  s=15)# 突出起点和终点plt.scatter(0, 0, c='green', edgecolors='none', s=100)plt.scatter(rw.x_values[-1], rw.y_values[-1], c='red', edgecolors='none', s=100)# 隐藏坐标轴plt.axes().get_xaxis().set_visible(False)plt.axes().get_yaxis().set_visible(False)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if keep_running == 'n':break


图15-10 包含50 000个点的随机漫步

这个示例模拟了一次包含50000个点的随机漫步(以模拟现实情况),并将每个点的大小都设置为1。最终的随机漫步图更纤细,犹如云朵,如图15-10所示。正如你看到的,我们使用简单的散点图制作出了一件艺术品!
请尝试修改上述代码,看看将漫步包含的点数增加到多少后,程序的运行速度变得极其缓慢或绘制出的图形变得很难看。

15.3.10调整尺寸以适合屏幕
图表适合屏幕大小时,更能有效地将数据中的规律呈现出来。为让绘图窗口更适合屏幕大小,可像下面这样调整matplotlib输出的尺寸:
rw_visual.py

import  matplotlib.pyplot  as  plt
from  random_walk  import  RandomWalk
# 只要程序处于活动状态,就不断地模拟随机漫步
# while  True:
# #  创建一个RandomWalk实例,并将其包含的点都绘制出来
#     rw  =  RandomWalk()
#     rw.fill_walk()
#     plt.scatter(rw.x_values,  rw.y_values,  s=15)
#     plt.show()
#
#     keep_running  =  input("Make  another  walk?  (y/n):  ")
#     if  keep_running  ==  'n':
#         breakwhile  True:#  创建一个RandomWalk实例,并将其包含的点都绘制出来rw  =  RandomWalk()rw.fill_walk()# 设置绘图窗口的尺寸plt.figure(figsize=(10, 6))# 绘制点并将图形显示出来point_numbers = list(range(rw.num_points))plt.scatter(rw.x_values, rw.y_values, c=point_numbers, cmap=plt.cm.Blues, edgecolor='none', s=1)point_numbers  =  list(range(rw.num_points))plt.scatter(rw.x_values,  rw.y_values,  c=point_numbers,  cmap=plt.cm.Blues, edgecolor='none',  s=15)# 突出起点和终点plt.scatter(0, 0, c='green', edgecolors='none', s=100)plt.scatter(rw.x_values[-1], rw.y_values[-1], c='red', edgecolors='none', s=100)# 隐藏坐标轴plt.axes().get_xaxis().set_visible(False)plt.axes().get_yaxis().set_visible(False)plt.show()keep_running  =  input("Make  another  walk?  (y/n):  ")if keep_running == 'n':break

函数figure() 用于指定图表的宽度、高度、分辨率和背景色。你需要给形参figsize 指定一个元组,向matplotlib指出绘图窗口的尺寸,单位为英寸。
Python假定屏幕分辨率为80像素/英寸,如果上述代码指定的图表尺寸不合适,可根据需要调整其中的数字。如果你知道自己的系统的分辨率,可使用形参dpi 向figure() 传递该分辨率,以有效地利用可用的屏幕空间,如下所示:

plt.figure(dpi=128,  figsize=(10,  6))

15.4使用Pygal模拟掷骰子
在本节中,我们将使用Python可视化包Pygal来生成可缩放的矢量图形文件。对于需要在尺寸不同的屏幕上显示的图表,这很有用,因为它们将自动缩放,以适合观看者的屏幕。 如果你打算以在线方式使用图表,请考虑使用Pygal来生成它们,这样它们在任何设备上显示时都会很美观。
在这个项目中,我们将对掷骰子的结果进行分析。掷6面的常规骰子时,可能出现的结果为1~6点,且出现每种结果的可能性相同。然而,如果同时掷两个骰子,某些点数出现的可能性将比其他点数大。为确定哪些点数出现的可能性最大,我们将生成一个表示掷骰子结果的数据集,并根据结果绘制出一个图形。
在数学领域,常常利用掷骰子来解释各种数据分析,但它在赌场和其他博弈场景中也得到了实际应用,在游戏《大富翁》以及众多角色扮演游戏中亦如此。

15.4.1安装Pygal
请使用pip 来安装Pygal(如果还未使用过pip ,请参阅12.2.1节)。

在Linux和OS X系统中,应执行的命令类似于下面这样:

pip  install  --user  pygal

在Windows系统中,命令类似于下面这样:

python  -m  pip  install  --user  pygal

要在python安装路径下执行

注意 你可能需要使用命令pip3 而不是pip ,如果这还是不管用,你可能需要删除标志–user 。

15.4.2Pygal画廊
要了解使用Pygal可创建什么样的图表,请查看图表类型画廊:访问http://www.pygal.org/ ,单击Documentation,再单击Chart types。每个示例都包含源代码,让你知道这些图表是如何生成的。
15.4.3创建Die 类
下面的类模拟掷一个骰子:
die.py

from  random  import  randint
class  Die():"""表示一个骰子的类"""
def  __init__(self,  num_sides=6):""" 骰 子 默 认 为 6 面 """self.num_sides  =  num_sidesdef  roll(self):""""返回一个位于1和骰子面数之间的随机值"""return  randint(1,  self.num_sides)

方法__init__() 接受一个可选参数。创建这个类的实例时,如果没有指定任何实参,面数默认为6;如果指定了实参,这个值将用于设置骰子的面数(见❶)。骰子是根据面数命名的,6面的骰子名为D6,8面的骰子名为D8,以此类推。
方法roll() 使用函数randint() 来返回一个1和面数之间的随机数(见❷)。这个函数可能返回起始值1、终止值num_sides 或这两个值之间的任何整数。

15.4.4掷骰子

from  die  import  Die
# 创建一个D6
die  =  Die()# 掷几次骰子,并将结果存储在一个列表中
results  =  []
for  roll_num  in  range(100):result  =  die.roll()results.append(result)print(results)

使用这个类来创建图表前,先来掷D6骰子,将结果打印出来,并检查结果是否合理:
die_visual.py

在❶处,我们创建了一个Die 实例,其面数为默认值6。在❷处,我们掷骰子100次,并将每次的结果都存储在列表results 中。下面是一个示例结果集:

[1, 2, 1, 4, 4, 2, 5, 5, 4, 4, 2, 5, 1, 5, 2, 6, 1, 1, 1, 2, 1, 4, 4, 1, 3, 1, 1, 1, 2, 2, 2, 4, 5, 1, 1, 4, 3, 1, 3, 3, 2, 2, 5, 2, 6, 5, 6, 2, 3, 2, 6, 4, 6, 1, 1, 6, 5, 5, 2, 1, 1, 3, 3, 3, 3, 4, 6, 4, 3, 5, 3, 6, 5, 4, 4, 1, 1, 4, 2, 4, 1, 3, 1, 6, 2, 1, 5, 4, 4, 4, 6, 4, 3, 6, 3, 2, 2, 4, 6, 6]

通过快速扫描这些结果可知,Die 类看起来没有问题。我们见到了值1和6,这表明返回了最大和最小的可能值;我们没有见到0或7,这表明结果都在正确的范围内。我们还看到了1~6的所有数字,这表明所有可能的结果都出现了。

15.4.5分析结果
为分析掷一个D6骰子的结果,我们计算每个点数出现的次数:
die_visual.py

from  die  import  Die
# 创建一个D6
die  =  Die()# 掷几次骰子,并将结果存储在一个列表中
results  =  []
for  roll_num  in  range(100):result  =  die.roll()results.append(result)# 分析结果frequencies  =  []for  value  in  range(1,  die.num_sides+1):frequency  =  results.count(value)frequencies.append(frequency)print(frequencies)

由于我们将使用Pygal来进行分析,而不是将结果打印出来,因此可以将模拟掷骰子的次数增加到1000(见❶)。为分析结果,我们创建了空列表frequencies ,用于存储每种点数出现的次数。在❷处,我们遍历可能的点数(这里为1~6),计算每种点数在results 中出现了多少次(见❸),并将这个值附加到列表frequencies 的末尾(见❹)。接下来,我们在可视化之前将这个列表打印出来:

[1, 0, 0, 0, 0, 0]
[1, 0, 0, 1, 0, 0]
[1, 0, 0, 1, 1, 0]
[1, 1, 0, 1, 1, 0]
[1, 1, 0, 2, 1, 0]
[1, 2, 0, 2, 1, 0]
[1, 2, 0, 2, 2, 0]
[1, 3, 0, 2, 2, 0]
[1, 4, 0, 2, 2, 0]
[1, 4, 0, 3, 2, 0]
[1, 4, 0, 4, 2, 0]
[1, 5, 0, 4, 2, 0]
[1, 5, 0, 4, 3, 0]
[1, 5, 0, 4, 3, 1]
[2, 5, 0, 4, 3, 1]
[2, 6, 0, 4, 3, 1]
[3, 6, 0, 4, 3, 1]
[3, 6, 1, 4, 3, 1]
[3, 6, 2, 4, 3, 1]
[3, 6, 2, 4, 4, 1]
[3, 7, 2, 4, 4, 1]
[3, 7, 3, 4, 4, 1]
[3, 7, 3, 5, 4, 1]
[3, 8, 3, 5, 4, 1]
[3, 8, 3, 6, 4, 1]
[3, 8, 3, 7, 4, 1]
[4, 8, 3, 7, 4, 1]
[4, 8, 3, 7, 5, 1]
[4, 9, 3, 7, 5, 1]
[4, 9, 3, 8, 5, 1]
[4, 9, 3, 8, 5, 2]
[4, 9, 3, 8, 5, 3]
[4, 10, 3, 8, 5, 3]
[4, 10, 3, 8, 6, 3]
[4, 11, 3, 8, 6, 3]
[4, 11, 4, 8, 6, 3]
[4, 11, 4, 9, 6, 3]
[4, 11, 4, 9, 6, 4]
[4, 11, 4, 10, 6, 4]
[4, 12, 4, 10, 6, 4]
[4, 12, 5, 10, 6, 4]
[4, 12, 5, 11, 6, 4]
[4, 12, 5, 12, 6, 4]
[4, 12, 5, 12, 6, 5]
[4, 12, 5, 13, 6, 5]
[4, 13, 5, 13, 6, 5]
[5, 13, 5, 13, 6, 5]
[5, 13, 6, 13, 6, 5]
[5, 13, 6, 13, 7, 5]
[6, 13, 6, 13, 7, 5]
[7, 13, 6, 13, 7, 5]
[7, 13, 7, 13, 7, 5]
[7, 13, 7, 14, 7, 5]
[7, 13, 7, 14, 7, 6]
[7, 13, 7, 14, 7, 7]
[7, 13, 7, 14, 7, 8]
[7, 13, 7, 14, 8, 8]
[7, 13, 7, 14, 9, 8]
[7, 13, 7, 15, 9, 8]
[7, 13, 7, 15, 9, 9]
[7, 14, 7, 15, 9, 9]
[7, 14, 7, 15, 9, 10]
[7, 14, 8, 15, 9, 10]
[7, 14, 8, 15, 9, 11]
[7, 14, 8, 16, 9, 11]
[7, 14, 9, 16, 9, 11]
[7, 14, 9, 16, 10, 11]
[7, 14, 9, 16, 10, 12]
[7, 14, 9, 16, 11, 12]
[7, 15, 9, 16, 11, 12]
[7, 15, 9, 16, 12, 12]
[8, 15, 9, 16, 12, 12]
[8, 15, 9, 17, 12, 12]
[9, 15, 9, 17, 12, 12]
[9, 15, 9, 17, 13, 12]
[9, 16, 9, 17, 13, 12]
[10, 16, 9, 17, 13, 12]
[10, 16, 9, 17, 14, 12]
[10, 16, 10, 17, 14, 12]
[10, 17, 10, 17, 14, 12]
[10, 17, 10, 17, 14, 13]
[10, 17, 11, 17, 14, 13]
[10, 17, 11, 17, 14, 14]
[10, 17, 11, 17, 15, 14]
[10, 17, 11, 17, 16, 14]
[10, 17, 11, 17, 16, 15]
[10, 17, 11, 17, 16, 16]
[10, 17, 11, 17, 16, 17]
[10, 17, 12, 17, 16, 17]
[11, 17, 12, 17, 16, 17]
[11, 17, 12, 18, 16, 17]
[11, 18, 12, 18, 16, 17]
[12, 18, 12, 18, 16, 17]
[12, 18, 12, 18, 17, 17]
[12, 19, 12, 18, 17, 17]
[12, 19, 12, 18, 17, 18]
[12, 20, 12, 18, 17, 18]
[12, 20, 12, 18, 17, 19]
[12, 20, 12, 18, 17, 20]
[12, 20, 12, 18, 17, 21]
[12, 20, 13, 18, 17, 21]
[12, 20, 13, 18, 18, 21]
[12, 20, 13, 18, 18, 22]
[13, 20, 13, 18, 18, 22]
[13, 20, 14, 18, 18, 22]
[13, 20, 14, 19, 18, 22]
[13, 20, 15, 19, 18, 22]
[14, 20, 15, 19, 18, 22]
[14, 20, 15, 20, 18, 22]
[14, 20, 15, 20, 18, 23]
[14, 21, 15, 20, 18, 23]
[14, 21, 16, 20, 18, 23]
[14, 21, 16, 21, 18, 23]
[14, 21, 16, 21, 18, 24]
[14, 21, 16, 21, 18, 25]
[14, 22, 16, 21, 18, 25]
[14, 22, 17, 21, 18, 25]
[14, 22, 17, 22, 18, 25]
[14, 22, 17, 22, 18, 26]
[14, 22, 18, 22, 18, 26]
[14, 22, 19, 22, 18, 26]
[14, 22, 20, 22, 18, 26]
[14, 23, 20, 22, 18, 26]
[14, 23, 21, 22, 18, 26]
[15, 23, 21, 22, 18, 26]
[15, 24, 21, 22, 18, 26]
[15, 25, 21, 22, 18, 26]
[15, 26, 21, 22, 18, 26]
[15, 26, 21, 23, 18, 26]
[15, 26, 21, 23, 18, 27]
[15, 26, 22, 23, 18, 27]
[16, 26, 22, 23, 18, 27]
[16, 26, 22, 24, 18, 27]
[16, 26, 23, 24, 18, 27]
[16, 26, 23, 25, 18, 27]
[16, 27, 23, 25, 18, 27]
[16, 27, 24, 25, 18, 27]
[16, 27, 24, 25, 18, 28]
[17, 27, 24, 25, 18, 28]
[17, 27, 25, 25, 18, 28]
[17, 27, 26, 25, 18, 28]
[18, 27, 26, 25, 18, 28]
[18, 27, 26, 26, 18, 28]
[18, 27, 26, 26, 18, 29]
[18, 27, 26, 26, 18, 30]
[19, 27, 26, 26, 18, 30]
[19, 27, 27, 26, 18, 30]
[20, 27, 27, 26, 18, 30]
[21, 27, 27, 26, 18, 30]
[21, 27, 28, 26, 18, 30]
[21, 27, 29, 26, 18, 30]
[21, 27, 29, 26, 18, 31]
[22, 27, 29, 26, 18, 31]
[22, 27, 29, 26, 18, 32]
[22, 27, 29, 26, 19, 32]
[23, 27, 29, 26, 19, 32]
[23, 27, 29, 26, 19, 33]
[23, 27, 30, 26, 19, 33]
[23, 27, 31, 26, 19, 33]
[23, 27, 31, 27, 19, 33]
[23, 27, 31, 27, 20, 33]
[23, 27, 31, 27, 21, 33]
[24, 27, 31, 27, 21, 33]
[24, 27, 32, 27, 21, 33]
[25, 27, 32, 27, 21, 33]
[25, 27, 32, 27, 22, 33]
[25, 27, 33, 27, 22, 33]
[25, 27, 33, 27, 23, 33]
[25, 27, 33, 28, 23, 33]
[26, 27, 33, 28, 23, 33]
[26, 27, 34, 28, 23, 33]
[26, 27, 34, 28, 24, 33]
[26, 27, 34, 28, 25, 33]
[26, 28, 34, 28, 25, 33]
[26, 28, 34, 29, 25, 33]
[26, 29, 34, 29, 25, 33]
[26, 29, 34, 30, 25, 33]
[26, 29, 34, 31, 25, 33]
[26, 30, 34, 31, 25, 33]
[26, 30, 34, 31, 25, 34]
[26, 31, 34, 31, 25, 34]
[26, 31, 34, 31, 26, 34]
[26, 31, 35, 31, 26, 34]
[27, 31, 35, 31, 26, 34]
[27, 32, 35, 31, 26, 34]
[27, 32, 35, 32, 26, 34]
[27, 32, 36, 32, 26, 34]
[27, 32, 36, 32, 26, 35]
[28, 32, 36, 32, 26, 35]
[28, 33, 36, 32, 26, 35]
[28, 33, 36, 32, 27, 35]
[28, 33, 36, 32, 28, 35]
[28, 33, 36, 32, 29, 35]
[29, 33, 36, 32, 29, 35]
[30, 33, 36, 32, 29, 35]
[30, 33, 37, 32, 29, 35]
[30, 33, 37, 32, 29, 36]
[30, 33, 38, 32, 29, 36]
[30, 33, 38, 33, 29, 36]
[30, 34, 38, 33, 29, 36]
[31, 34, 38, 33, 29, 36]
[31, 35, 38, 33, 29, 36]
[31, 35, 38, 33, 29, 37]
[31, 35, 38, 33, 29, 38]
[31, 35, 38, 34, 29, 38]
[31, 35, 39, 34, 29, 38]
[31, 35, 40, 34, 29, 38]
[31, 35, 40, 34, 29, 39]
[31, 35, 40, 35, 29, 39]
[31, 35, 40, 36, 29, 39]
[31, 35, 41, 36, 29, 39]
[32, 35, 41, 36, 29, 39]
[32, 35, 41, 36, 29, 40]
[32, 36, 41, 36, 29, 40]
[32, 37, 41, 36, 29, 40]
[33, 37, 41, 36, 29, 40]
[33, 37, 41, 36, 30, 40]
[33, 38, 41, 36, 30, 40]
[33, 38, 42, 36, 30, 40]
[34, 38, 42, 36, 30, 40]
[34, 38, 42, 36, 31, 40]
[34, 39, 42, 36, 31, 40]
[34, 39, 42, 36, 32, 40]
[34, 39, 43, 36, 32, 40]
[34, 39, 43, 36, 33, 40]
[34, 39, 43, 36, 33, 41]
[34, 39, 44, 36, 33, 41]
[34, 39, 45, 36, 33, 41]
[34, 39, 45, 36, 33, 42]
[34, 39, 45, 36, 33, 43]
[34, 40, 45, 36, 33, 43]
[34, 40, 45, 36, 33, 44]
[34, 40, 46, 36, 33, 44]
[34, 40, 47, 36, 33, 44]
[34, 40, 48, 36, 33, 44]
[34, 41, 48, 36, 33, 44]
[34, 42, 48, 36, 33, 44]
[34, 42, 48, 37, 33, 44]
[34, 42, 49, 37, 33, 44]
[35, 42, 49, 37, 33, 44]
[35, 43, 49, 37, 33, 44]
[35, 43, 50, 37, 33, 44]
[35, 43, 50, 37, 34, 44]
[36, 43, 50, 37, 34, 44]
[36, 43, 50, 37, 34, 45]
[36, 43, 50, 38, 34, 45]
[36, 44, 50, 38, 34, 45]
[36, 45, 50, 38, 34, 45]
[36, 45, 51, 38, 34, 45]
[37, 45, 51, 38, 34, 45]
[38, 45, 51, 38, 34, 45]
[38, 45, 52, 38, 34, 45]
[38, 45, 52, 38, 35, 45]
[38, 45, 52, 39, 35, 45]
[39, 45, 52, 39, 35, 45]
[40, 45, 52, 39, 35, 45]
[41, 45, 52, 39, 35, 45]
[41, 45, 52, 39, 36, 45]
[41, 45, 52, 40, 36, 45]
[42, 45, 52, 40, 36, 45]
[43, 45, 52, 40, 36, 45]
[43, 46, 52, 40, 36, 45]
[43, 46, 52, 40, 37, 45]
[43, 46, 52, 40, 37, 46]
[43, 46, 53, 40, 37, 46]
[43, 46, 53, 40, 38, 46]
[43, 47, 53, 40, 38, 46]
[43, 47, 54, 40, 38, 46]
[43, 47, 54, 40, 39, 46]
[43, 47, 54, 40, 40, 46]
[43, 47, 55, 40, 40, 46]
[44, 47, 55, 40, 40, 46]
[44, 47, 55, 41, 40, 46]
[44, 47, 55, 42, 40, 46]
[44, 47, 55, 42, 40, 47]
[45, 47, 55, 42, 40, 47]
[46, 47, 55, 42, 40, 47]
[46, 48, 55, 42, 40, 47]
[46, 49, 55, 42, 40, 47]
[46, 50, 55, 42, 40, 47]
[46, 50, 55, 42, 40, 48]
[46, 50, 56, 42, 40, 48]
[47, 50, 56, 42, 40, 48]
[47, 50, 56, 42, 40, 49]
[47, 51, 56, 42, 40, 49]
[47, 51, 56, 43, 40, 49]
[47, 51, 57, 43, 40, 49]
[47, 51, 57, 43, 40, 50]
[47, 51, 57, 44, 40, 50]
[47, 51, 57, 44, 41, 50]
[47, 51, 58, 44, 41, 50]
[48, 51, 58, 44, 41, 50]
[48, 51, 58, 44, 42, 50]
[49, 51, 58, 44, 42, 50]
[49, 52, 58, 44, 42, 50]
[49, 53, 58, 44, 42, 50]
[50, 53, 58, 44, 42, 50]
[50, 53, 58, 44, 42, 51]
[50, 53, 59, 44, 42, 51]
[51, 53, 59, 44, 42, 51]
[51, 53, 59, 45, 42, 51]
[51, 53, 59, 45, 42, 52]
[51, 54, 59, 45, 42, 52]
[51, 55, 59, 45, 42, 52]
[51, 56, 59, 45, 42, 52]
[52, 56, 59, 45, 42, 52]
[52, 56, 59, 45, 43, 52]
[52, 57, 59, 45, 43, 52]
[52, 57, 59, 45, 44, 52]
[52, 58, 59, 45, 44, 52]
[53, 58, 59, 45, 44, 52]
[53, 58, 59, 45, 45, 52]
[54, 58, 59, 45, 45, 52]
[54, 58, 59, 45, 45, 53]
[54, 58, 59, 45, 46, 53]
[54, 58, 59, 46, 46, 53]
[54, 58, 59, 47, 46, 53]
[54, 58, 59, 48, 46, 53]
[55, 58, 59, 48, 46, 53]
[55, 58, 59, 48, 46, 54]
[56, 58, 59, 48, 46, 54]
[56, 58, 59, 48, 47, 54]
[57, 58, 59, 48, 47, 54]
[57, 58, 59, 48, 47, 55]
[57, 58, 59, 48, 48, 55]
[57, 58, 59, 48, 49, 55]
[57, 58, 59, 49, 49, 55]
[57, 59, 59, 49, 49, 55]
[57, 59, 59, 50, 49, 55]
[57, 59, 59, 51, 49, 55]
[57, 59, 59, 51, 49, 56]
[57, 59, 59, 51, 50, 56]
[57, 59, 59, 51, 51, 56]
[57, 59, 59, 51, 51, 57]
[57, 59, 59, 51, 51, 58]
[57, 59, 60, 51, 51, 58]
[57, 60, 60, 51, 51, 58]
[58, 60, 60, 51, 51, 58]
[58, 61, 60, 51, 51, 58]
[58, 61, 60, 52, 51, 58]
[59, 61, 60, 52, 51, 58]
[59, 61, 60, 52, 52, 58]
[59, 61, 60, 53, 52, 58]
[59, 61, 60, 53, 52, 59]
[59, 61, 60, 54, 52, 59]
[59, 61, 61, 54, 52, 59]
[60, 61, 61, 54, 52, 59]
[61, 61, 61, 54, 52, 59]
[61, 61, 62, 54, 52, 59]
[61, 61, 63, 54, 52, 59]
[61, 62, 63, 54, 52, 59]
[61, 62, 63, 54, 53, 59]
[62, 62, 63, 54, 53, 59]
[62, 62, 63, 54, 54, 59]
[62, 62, 64, 54, 54, 59]
[62, 62, 64, 54, 55, 59]
[62, 62, 64, 54, 55, 60]
[62, 62, 65, 54, 55, 60]
[62, 62, 66, 54, 55, 60]
[62, 62, 66, 54, 56, 60]
[62, 62, 66, 54, 57, 60]
[62, 62, 66, 55, 57, 60]
[62, 62, 66, 55, 58, 60]
[62, 62, 66, 55, 59, 60]
[62, 62, 66, 55, 60, 60]
[63, 62, 66, 55, 60, 60]
[63, 62, 66, 55, 60, 61]
[63, 63, 66, 55, 60, 61]
[63, 63, 67, 55, 60, 61]
[63, 63, 67, 55, 61, 61]
[64, 63, 67, 55, 61, 61]
[64, 63, 67, 56, 61, 61]
[64, 63, 68, 56, 61, 61]
[64, 63, 68, 56, 61, 62]
[64, 63, 69, 56, 61, 62]
[64, 64, 69, 56, 61, 62]
[65, 64, 69, 56, 61, 62]
[65, 64, 70, 56, 61, 62]
[65, 64, 71, 56, 61, 62]
[65, 64, 71, 56, 62, 62]
[65, 64, 71, 56, 62, 63]
[65, 64, 72, 56, 62, 63]
[65, 64, 72, 56, 62, 64]
[65, 64, 72, 57, 62, 64]
[65, 65, 72, 57, 62, 64]
[65, 65, 72, 58, 62, 64]
[65, 65, 72, 59, 62, 64]
[66, 65, 72, 59, 62, 64]
[67, 65, 72, 59, 62, 64]
[67, 65, 72, 59, 62, 65]
[67, 65, 72, 60, 62, 65]
[67, 65, 72, 61, 62, 65]
[67, 65, 73, 61, 62, 65]
[67, 65, 73, 61, 62, 66]
[67, 65, 73, 61, 62, 67]
[67, 65, 74, 61, 62, 67]
[67, 65, 74, 62, 62, 67]
[67, 65, 75, 62, 62, 67]
[67, 65, 75, 62, 63, 67]
[67, 66, 75, 62, 63, 67]
[67, 66, 75, 62, 63, 68]
[67, 66, 75, 62, 63, 69]
[67, 66, 76, 62, 63, 69]
[67, 66, 76, 63, 63, 69]
[67, 67, 76, 63, 63, 69]
[68, 67, 76, 63, 63, 69]
[68, 67, 76, 64, 63, 69]
[68, 67, 76, 65, 63, 69]
[68, 67, 76, 65, 63, 70]
[68, 68, 76, 65, 63, 70]
[68, 68, 76, 65, 64, 70]
[68, 69, 76, 65, 64, 70]
[68, 69, 76, 65, 65, 70]
[69, 69, 76, 65, 65, 70]
[69, 69, 76, 66, 65, 70]
[70, 69, 76, 66, 65, 70]
[70, 69, 77, 66, 65, 70]
[70, 69, 77, 66, 66, 70]
[70, 69, 77, 67, 66, 70]
[71, 69, 77, 67, 66, 70]
[71, 69, 77, 68, 66, 70]
[71, 70, 77, 68, 66, 70]
[72, 70, 77, 68, 66, 70]
[72, 71, 77, 68, 66, 70]
[72, 71, 77, 68, 67, 70]
[72, 71, 77, 69, 67, 70]
[72, 71, 78, 69, 67, 70]
[72, 72, 78, 69, 67, 70]
[72, 72, 78, 69, 68, 70]
[73, 72, 78, 69, 68, 70]
[73, 72, 78, 69, 69, 70]
[73, 72, 78, 70, 69, 70]
[73, 72, 78, 70, 69, 71]
[73, 73, 78, 70, 69, 71]
[74, 73, 78, 70, 69, 71]
[74, 73, 78, 70, 69, 72]
[74, 73, 78, 70, 70, 72]
[74, 73, 78, 70, 70, 73]
[74, 73, 78, 70, 70, 74]
[74, 73, 78, 70, 70, 75]
[74, 73, 79, 70, 70, 75]
[74, 73, 79, 71, 70, 75]
[74, 73, 79, 71, 70, 76]
[74, 73, 79, 72, 70, 76]
[74, 73, 79, 73, 70, 76]
[74, 74, 79, 73, 70, 76]
[75, 74, 79, 73, 70, 76]
[75, 74, 79, 73, 71, 76]
[75, 74, 80, 73, 71, 76]
[75, 74, 80, 73, 71, 77]
[75, 74, 80, 73, 72, 77]
[75, 75, 80, 73, 72, 77]
[75, 76, 80, 73, 72, 77]
[75, 76, 80, 74, 72, 77]
[76, 76, 80, 74, 72, 77]
[76, 76, 80, 74, 73, 77]
[76, 76, 80, 74, 73, 78]
[76, 77, 80, 74, 73, 78]
[77, 77, 80, 74, 73, 78]
[77, 77, 80, 75, 73, 78]
[77, 77, 80, 75, 74, 78]
[77, 77, 80, 75, 74, 79]
[77, 77, 81, 75, 74, 79]
[77, 77, 82, 75, 74, 79]
[78, 77, 82, 75, 74, 79]
[79, 77, 82, 75, 74, 79]
[80, 77, 82, 75, 74, 79]
[80, 77, 82, 75, 74, 80]
[80, 78, 82, 75, 74, 80]
[80, 78, 82, 75, 75, 80]
[80, 78, 82, 75, 76, 80]
[80, 78, 82, 76, 76, 80]
[80, 78, 82, 76, 77, 80]
[80, 78, 82, 76, 78, 80]
[80, 78, 82, 76, 79, 80]
[80, 78, 82, 76, 79, 81]
[80, 78, 83, 76, 79, 81]
[80, 78, 83, 77, 79, 81]
[80, 78, 83, 78, 79, 81]
[80, 78, 83, 78, 79, 82]
[81, 78, 83, 78, 79, 82]
[81, 79, 83, 78, 79, 82]
[81, 79, 84, 78, 79, 82]
[81, 79, 84, 78, 79, 83]
[81, 79, 85, 78, 79, 83]
[81, 80, 85, 78, 79, 83]
[81, 80, 86, 78, 79, 83]
[81, 80, 86, 79, 79, 83]
[81, 80, 86, 79, 79, 84]
[81, 80, 86, 80, 79, 84]
[81, 81, 86, 80, 79, 84]
[81, 81, 86, 80, 80, 84]
[81, 81, 87, 80, 80, 84]
[81, 81, 88, 80, 80, 84]
[81, 81, 88, 80, 80, 85]
[82, 81, 88, 80, 80, 85]
[82, 81, 88, 81, 80, 85]
[82, 82, 88, 81, 80, 85]
[82, 83, 88, 81, 80, 85]
[82, 83, 89, 81, 80, 85]
[82, 83, 89, 82, 80, 85]
[83, 83, 89, 82, 80, 85]
[83, 83, 89, 82, 80, 86]
[83, 84, 89, 82, 80, 86]
[84, 84, 89, 82, 80, 86]
[84, 85, 89, 82, 80, 86]
[84, 85, 90, 82, 80, 86]
[85, 85, 90, 82, 80, 86]
[85, 86, 90, 82, 80, 86]
[85, 87, 90, 82, 80, 86]
[85, 88, 90, 82, 80, 86]
[86, 88, 90, 82, 80, 86]
[86, 88, 91, 82, 80, 86]
[86, 88, 92, 82, 80, 86]
[86, 88, 93, 82, 80, 86]
[86, 88, 93, 83, 80, 86]
[86, 88, 93, 83, 81, 86]
[86, 88, 93, 83, 82, 86]
[86, 89, 93, 83, 82, 86]
[87, 89, 93, 83, 82, 86]
[87, 89, 93, 83, 82, 87]
[87, 89, 93, 83, 83, 87]
[87, 90, 93, 83, 83, 87]
[87, 90, 93, 83, 83, 88]
[87, 90, 93, 83, 83, 89]
[87, 90, 93, 83, 84, 89]
[87, 90, 93, 83, 84, 90]
[87, 91, 93, 83, 84, 90]
[87, 91, 93, 83, 85, 90]
[88, 91, 93, 83, 85, 90]
[88, 91, 93, 84, 85, 90]
[88, 91, 94, 84, 85, 90]
[88, 92, 94, 84, 85, 90]
[89, 92, 94, 84, 85, 90]
[90, 92, 94, 84, 85, 90]
[91, 92, 94, 84, 85, 90]
[91, 92, 94, 85, 85, 90]
[91, 93, 94, 85, 85, 90]
[92, 93, 94, 85, 85, 90]
[92, 93, 94, 85, 85, 91]
[93, 93, 94, 85, 85, 91]
[93, 93, 94, 85, 85, 92]
[93, 93, 95, 85, 85, 92]
[94, 93, 95, 85, 85, 92]
[94, 93, 96, 85, 85, 92]
[95, 93, 96, 85, 85, 92]
[95, 93, 96, 85, 86, 92]
[95, 93, 96, 85, 86, 93]
[95, 93, 96, 85, 86, 94]
[95, 93, 96, 85, 87, 94]
[96, 93, 96, 85, 87, 94]
[96, 93, 96, 85, 88, 94]
[96, 93, 96, 85, 88, 95]
[96, 93, 96, 85, 88, 96]
[96, 93, 96, 85, 89, 96]
[97, 93, 96, 85, 89, 96]
[98, 93, 96, 85, 89, 96]
[98, 94, 96, 85, 89, 96]
[98, 94, 96, 86, 89, 96]
[98, 94, 96, 86, 89, 97]
[98, 94, 96, 86, 89, 98]
[98, 94, 97, 86, 89, 98]
[98, 94, 97, 86, 89, 99]
[98, 94, 97, 86, 90, 99]
[98, 95, 97, 86, 90, 99]
[98, 95, 97, 87, 90, 99]
[98, 95, 97, 87, 91, 99]
[98, 95, 98, 87, 91, 99]
[98, 96, 98, 87, 91, 99]
[98, 96, 98, 88, 91, 99]
[98, 96, 98, 88, 91, 100]
[98, 96, 98, 88, 92, 100]
[98, 97, 98, 88, 92, 100]
[98, 97, 98, 88, 93, 100]
[98, 97, 98, 88, 93, 101]
[99, 97, 98, 88, 93, 101]
[100, 97, 98, 88, 93, 101]
[100, 97, 99, 88, 93, 101]
[100, 97, 100, 88, 93, 101]
[100, 97, 101, 88, 93, 101]
[100, 97, 101, 89, 93, 101]
[100, 97, 102, 89, 93, 101]
[101, 97, 102, 89, 93, 101]
[101, 97, 102, 89, 94, 101]
[102, 97, 102, 89, 94, 101]
[102, 97, 102, 90, 94, 101]
[102, 97, 102, 90, 94, 102]
[103, 97, 102, 90, 94, 102]
[103, 97, 102, 91, 94, 102]
[103, 97, 102, 91, 95, 102]
[103, 98, 102, 91, 95, 102]
[104, 98, 102, 91, 95, 102]
[104, 98, 102, 91, 96, 102]
[104, 98, 102, 91, 96, 103]
[104, 98, 103, 91, 96, 103]
[104, 98, 103, 91, 96, 104]
[104, 98, 103, 92, 96, 104]
[104, 98, 103, 92, 96, 105]
[104, 98, 103, 93, 96, 105]
[104, 98, 103, 93, 97, 105]
[104, 98, 103, 93, 97, 106]
[105, 98, 103, 93, 97, 106]
[105, 98, 103, 94, 97, 106]
[105, 98, 103, 94, 97, 107]
[105, 98, 103, 95, 97, 107]
[105, 98, 103, 96, 97, 107]
[105, 98, 103, 97, 97, 107]
[105, 98, 103, 97, 97, 108]
[105, 98, 103, 97, 98, 108]
[106, 98, 103, 97, 98, 108]
[106, 98, 104, 97, 98, 108]
[107, 98, 104, 97, 98, 108]
[107, 98, 104, 97, 99, 108]
[107, 98, 104, 98, 99, 108]
[107, 98, 104, 99, 99, 108]
[107, 98, 104, 99, 100, 108]
[107, 98, 104, 100, 100, 108]
[107, 98, 104, 101, 100, 108]
[107, 99, 104, 101, 100, 108]
[108, 99, 104, 101, 100, 108]
[108, 99, 104, 102, 100, 108]
[108, 100, 104, 102, 100, 108]
[108, 100, 104, 102, 100, 109]
[108, 100, 104, 102, 101, 109]
[108, 100, 104, 102, 101, 110]
[109, 100, 104, 102, 101, 110]
[109, 100, 104, 102, 101, 111]
[109, 100, 104, 103, 101, 111]
[109, 100, 104, 104, 101, 111]
[109, 100, 104, 104, 101, 112]
[109, 100, 105, 104, 101, 112]
[110, 100, 105, 104, 101, 112]
[110, 100, 105, 104, 101, 113]
[110, 100, 105, 105, 101, 113]
[110, 101, 105, 105, 101, 113]
[110, 101, 105, 105, 102, 113]
[110, 101, 105, 105, 102, 114]
[110, 101, 105, 105, 103, 114]
[110, 101, 105, 105, 103, 115]
[110, 101, 105, 105, 103, 116]
[110, 102, 105, 105, 103, 116]
[110, 102, 105, 106, 103, 116]
[110, 102, 105, 107, 103, 116]
[110, 102, 105, 107, 104, 116]
[110, 102, 105, 107, 105, 116]
[111, 102, 105, 107, 105, 116]
[111, 102, 105, 107, 105, 117]
[111, 102, 105, 107, 105, 118]
[111, 102, 105, 108, 105, 118]
[111, 102, 105, 108, 105, 119]
[111, 102, 105, 108, 106, 119]
[111, 102, 105, 108, 106, 120]
[111, 102, 105, 108, 106, 121]
[111, 103, 105, 108, 106, 121]
[111, 103, 105, 108, 106, 122]
[112, 103, 105, 108, 106, 122]
[112, 103, 105, 109, 106, 122]
[112, 103, 105, 109, 106, 123]
[112, 103, 106, 109, 106, 123]
[113, 103, 106, 109, 106, 123]
[113, 103, 106, 109, 107, 123]
[114, 103, 106, 109, 107, 123]
[114, 103, 106, 109, 107, 124]
[115, 103, 106, 109, 107, 124]
[115, 103, 106, 109, 107, 125]
[115, 103, 106, 109, 107, 126]
[115, 103, 106, 109, 108, 126]
[116, 103, 106, 109, 108, 126]
[117, 103, 106, 109, 108, 126]
[117, 103, 107, 109, 108, 126]
[117, 103, 107, 109, 108, 127]
[117, 104, 107, 109, 108, 127]
[117, 104, 107, 109, 108, 128]
[117, 105, 107, 109, 108, 128]
[117, 105, 108, 109, 108, 128]
[118, 105, 108, 109, 108, 128]
[118, 106, 108, 109, 108, 128]
[119, 106, 108, 109, 108, 128]
[120, 106, 108, 109, 108, 128]
[121, 106, 108, 109, 108, 128]
[121, 106, 108, 109, 108, 129]
[121, 106, 108, 109, 108, 130]
[121, 106, 108, 110, 108, 130]
[122, 106, 108, 110, 108, 130]
[123, 106, 108, 110, 108, 130]
[123, 107, 108, 110, 108, 130]
[123, 107, 108, 110, 109, 130]
[123, 107, 109, 110, 109, 130]
[123, 107, 110, 110, 109, 130]
[123, 108, 110, 110, 109, 130]
[124, 108, 110, 110, 109, 130]
[124, 109, 110, 110, 109, 130]
[124, 109, 110, 110, 110, 130]
[124, 110, 110, 110, 110, 130]
[125, 110, 110, 110, 110, 130]
[126, 110, 110, 110, 110, 130]
[127, 110, 110, 110, 110, 130]
[127, 110, 110, 110, 111, 130]
[127, 110, 110, 110, 111, 131]
[127, 110, 110, 110, 112, 131]
[128, 110, 110, 110, 112, 131]
[129, 110, 110, 110, 112, 131]
[129, 110, 110, 110, 113, 131]
[129, 110, 110, 111, 113, 131]
[129, 110, 110, 112, 113, 131]
[129, 110, 110, 113, 113, 131]
[129, 110, 111, 113, 113, 131]
[129, 110, 112, 113, 113, 131]
[129, 110, 112, 113, 114, 131]
[129, 110, 112, 113, 115, 131]
[129, 110, 112, 113, 116, 131]
[129, 110, 112, 113, 117, 131]
[129, 110, 112, 114, 117, 131]
[129, 110, 112, 114, 117, 132]
[129, 111, 112, 114, 117, 132]
[129, 111, 113, 114, 117, 132]
[130, 111, 113, 114, 117, 132]
[130, 111, 113, 115, 117, 132]
[131, 111, 113, 115, 117, 132]
[131, 111, 113, 116, 117, 132]
[131, 111, 113, 116, 118, 132]
[131, 111, 113, 117, 118, 132]
[131, 111, 113, 118, 118, 132]
[132, 111, 113, 118, 118, 132]
[132, 111, 113, 118, 119, 132]
[132, 111, 113, 118, 119, 133]
[133, 111, 113, 118, 119, 133]
[133, 111, 113, 119, 119, 133]
[133, 111, 113, 120, 119, 133]
[133, 112, 113, 120, 119, 133]
[133, 113, 113, 120, 119, 133]
[133, 113, 113, 121, 119, 133]
[133, 113, 114, 121, 119, 133]
[133, 113, 114, 121, 120, 133]
[134, 113, 114, 121, 120, 133]
[134, 114, 114, 121, 120, 133]
[134, 114, 115, 121, 120, 133]
[135, 114, 115, 121, 120, 133]
[135, 114, 115, 121, 121, 133]
[136, 114, 115, 121, 121, 133]
[136, 114, 115, 122, 121, 133]
[136, 115, 115, 122, 121, 133]
[137, 115, 115, 122, 121, 133]
[137, 116, 115, 122, 121, 133]
[137, 117, 115, 122, 121, 133]
[137, 117, 115, 122, 122, 133]
[137, 117, 116, 122, 122, 133]
[137, 118, 116, 122, 122, 133]
[137, 118, 116, 122, 123, 133]
[138, 118, 116, 122, 123, 133]
[139, 118, 116, 122, 123, 133]
[140, 118, 116, 122, 123, 133]
[141, 118, 116, 122, 123, 133]
[141, 119, 116, 122, 123, 133]
[141, 119, 116, 122, 123, 134]
[141, 119, 116, 122, 124, 134]
[141, 120, 116, 122, 124, 134]
[141, 120, 116, 122, 125, 134]
[141, 121, 116, 122, 125, 134]
[141, 121, 116, 122, 125, 135]
[141, 122, 116, 122, 125, 135]
[141, 122, 117, 122, 125, 135]
[141, 122, 117, 122, 126, 135]
[141, 122, 118, 122, 126, 135]
[141, 123, 118, 122, 126, 135]
[141, 123, 118, 122, 126, 136]
[141, 124, 118, 122, 126, 136]
[141, 124, 118, 122, 126, 137]
[141, 124, 118, 123, 126, 137]
[141, 125, 118, 123, 126, 137]
[141, 125, 118, 123, 126, 138]
[141, 125, 118, 123, 126, 139]
[141, 125, 118, 123, 126, 140]
[141, 125, 118, 123, 127, 140]
[141, 125, 119, 123, 127, 140]
[142, 125, 119, 123, 127, 140]
[142, 125, 119, 123, 128, 140]
[142, 125, 119, 123, 128, 141]
[142, 125, 119, 124, 128, 141]
[143, 125, 119, 124, 128, 141]
[143, 126, 119, 124, 128, 141]
[143, 127, 119, 124, 128, 141]
[144, 127, 119, 124, 128, 141]
[144, 127, 119, 124, 128, 142]
[145, 127, 119, 124, 128, 142]
[145, 127, 119, 125, 128, 142]
[145, 128, 119, 125, 128, 142]
[145, 128, 119, 125, 129, 142]
[145, 129, 119, 125, 129, 142]
[146, 129, 119, 125, 129, 142]
[146, 129, 119, 125, 129, 143]
[146, 129, 119, 125, 130, 143]
[146, 129, 119, 125, 130, 144]
[147, 129, 119, 125, 130, 144]
[147, 129, 119, 125, 130, 145]
[147, 129, 120, 125, 130, 145]
[147, 129, 120, 125, 131, 145]
[147, 129, 120, 126, 131, 145]
[147, 129, 121, 126, 131, 145]
[147, 130, 121, 126, 131, 145]
[147, 131, 121, 126, 131, 145]
[147, 131, 121, 126, 131, 146]
[147, 131, 121, 127, 131, 146]
[147, 131, 121, 127, 131, 147]
[147, 131, 122, 127, 131, 147]
[148, 131, 122, 127, 131, 147]
[149, 131, 122, 127, 131, 147]
[149, 132, 122, 127, 131, 147]
[149, 132, 123, 127, 131, 147]
[149, 132, 124, 127, 131, 147]
[149, 133, 124, 127, 131, 147]
[149, 133, 124, 128, 131, 147]
[149, 134, 124, 128, 131, 147]
[149, 134, 124, 129, 131, 147]
[149, 134, 124, 130, 131, 147]
[149, 134, 124, 131, 131, 147]
[149, 134, 124, 131, 132, 147]
[149, 134, 124, 131, 132, 148]
[149, 134, 124, 131, 132, 149]
[149, 134, 124, 131, 132, 150]
[149, 134, 124, 131, 133, 150]
[149, 134, 124, 131, 134, 150]
[150, 134, 124, 131, 134, 150]
[150, 134, 124, 131, 134, 151]
[150, 135, 124, 131, 134, 151]
[150, 136, 124, 131, 134, 151]
[150, 136, 124, 131, 135, 151]
[150, 136, 124, 131, 135, 152]
[150, 136, 124, 131, 135, 153]
[150, 136, 124, 131, 136, 153]
[150, 136, 125, 131, 136, 153]
[150, 137, 125, 131, 136, 153]
[150, 137, 125, 131, 137, 153]
[150, 138, 125, 131, 137, 153]
[150, 138, 125, 131, 138, 153]
[150, 138, 125, 131, 139, 153]
[150, 138, 125, 131, 139, 154]
[150, 138, 125, 131, 140, 154]
[150, 138, 125, 131, 140, 155]
[151, 138, 125, 131, 140, 155]
[151, 138, 125, 131, 141, 155]
[151, 138, 125, 131, 142, 155]
[151, 138, 126, 131, 142, 155]
[151, 138, 126, 131, 143, 155]
[151, 139, 126, 131, 143, 155]
[151, 139, 127, 131, 143, 155]
[151, 139, 127, 131, 143, 156]
[151, 140, 127, 131, 143, 156]
[151, 140, 127, 132, 143, 156]
[151, 140, 128, 132, 143, 156]
[151, 140, 128, 132, 143, 157]
[152, 140, 128, 132, 143, 157]
[152, 141, 128, 132, 143, 157]
[152, 141, 128, 132, 144, 157]
[152, 141, 129, 132, 144, 157]
[152, 141, 129, 132, 144, 158]
[152, 141, 129, 132, 144, 159]
[153, 141, 129, 132, 144, 159]
[153, 141, 129, 132, 145, 159]
[153, 141, 129, 133, 145, 159]
[153, 141, 130, 133, 145, 159]
[153, 141, 131, 133, 145, 159]
[153, 141, 131, 133, 146, 159]
[154, 141, 131, 133, 146, 159]
[154, 142, 131, 133, 146, 159]
[155, 142, 131, 133, 146, 159]
[155, 142, 131, 134, 146, 159]
[155, 142, 132, 134, 146, 159]
[155, 142, 132, 135, 146, 159]
[155, 142, 132, 135, 147, 159]
[155, 142, 133, 135, 147, 159]
[155, 142, 134, 135, 147, 159]
[155, 142, 134, 135, 147, 160]
[155, 143, 134, 135, 147, 160]
[155, 143, 135, 135, 147, 160]
[155, 143, 135, 136, 147, 160]
[156, 143, 135, 136, 147, 160]
[156, 144, 135, 136, 147, 160]
[156, 144, 136, 136, 147, 160]
[156, 144, 136, 136, 148, 160]
[156, 144, 137, 136, 148, 160]
[156, 144, 137, 136, 148, 161]
[156, 145, 137, 136, 148, 161]
[156, 145, 138, 136, 148, 161]
[157, 145, 138, 136, 148, 161]
[158, 145, 138, 136, 148, 161]
[159, 145, 138, 136, 148, 161]
[159, 145, 139, 136, 148, 161]
[159, 145, 139, 136, 148, 162]
[159, 145, 139, 137, 148, 162]
[159, 145, 139, 138, 148, 162]
[159, 145, 139, 138, 148, 163]
[159, 145, 139, 139, 148, 163]
[159, 145, 139, 140, 148, 163]
[159, 145, 140, 140, 148, 163]
[159, 145, 140, 140, 149, 163]
[159, 145, 140, 140, 150, 163]
[159, 145, 140, 141, 150, 163]
[159, 146, 140, 141, 150, 163]
[159, 146, 141, 141, 150, 163]
[159, 146, 141, 141, 151, 163]
[160, 146, 141, 141, 151, 163]
[160, 146, 142, 141, 151, 163]
[160, 146, 142, 141, 152, 163]
[160, 146, 142, 142, 152, 163]
[160, 146, 142, 142, 152, 164]
[160, 146, 142, 142, 153, 164]
[160, 146, 142, 142, 153, 165]
[160, 146, 142, 142, 154, 165]
[160, 146, 142, 142, 155, 165]
[160, 146, 143, 142, 155, 165]
[160, 146, 143, 143, 155, 165]
[160, 146, 143, 144, 155, 165]
[160, 146, 144, 144, 155, 165]
[160, 146, 144, 145, 155, 165]
[160, 146, 145, 145, 155, 165]
[160, 146, 145, 146, 155, 165]
[161, 146, 145, 146, 155, 165]
[161, 147, 145, 146, 155, 165]
[161, 147, 146, 146, 155, 165]
[161, 147, 146, 146, 155, 166]
[161, 148, 146, 146, 155, 166]
[162, 148, 146, 146, 155, 166]
[162, 149, 146, 146, 155, 166]
[162, 149, 146, 146, 155, 167]
[163, 149, 146, 146, 155, 167]
[163, 150, 146, 146, 155, 167]
[164, 150, 146, 146, 155, 167]
[165, 150, 146, 146, 155, 167]
[165, 150, 147, 146, 155, 167]
[165, 150, 148, 146, 155, 167]
[165, 150, 148, 146, 155, 168]
[165, 150, 148, 146, 156, 168]
[165, 151, 148, 146, 156, 168]
[165, 151, 149, 146, 156, 168]
[165, 151, 150, 146, 156, 168]
[165, 151, 150, 146, 156, 169]
[165, 151, 151, 146, 156, 169]
[166, 151, 151, 146, 156, 169]
[166, 151, 152, 146, 156, 169]
[166, 151, 153, 146, 156, 169]
[166, 151, 153, 146, 157, 169]
[167, 151, 153, 146, 157, 169]
[167, 151, 153, 147, 157, 169]
[167, 151, 153, 147, 158, 169]
[167, 151, 153, 148, 158, 169]
[167, 151, 154, 148, 158, 169]
[168, 151, 154, 148, 158, 169]
[168, 151, 155, 148, 158, 169]
[168, 151, 155, 149, 158, 169]
[169, 151, 155, 149, 158, 169]
[169, 151, 155, 149, 159, 169]
[170, 151, 155, 149, 159, 169]
[170, 152, 155, 149, 159, 169]
[170, 153, 155, 149, 159, 169]
[170, 153, 155, 149, 159, 170]
[170, 154, 155, 149, 159, 170]
[170, 154, 155, 149, 160, 170]
[170, 154, 155, 149, 161, 170]
[171, 154, 155, 149, 161, 170]
[171, 155, 155, 149, 161, 170]
[171, 155, 155, 150, 161, 170]
[171, 156, 155, 150, 161, 170]
[171, 157, 155, 150, 161, 170]
[171, 157, 156, 150, 161, 170]
[171, 157, 157, 150, 161, 170]
[171, 157, 157, 150, 162, 170]
[171, 157, 157, 150, 162, 171]
[171, 157, 158, 150, 162, 171]
[172, 157, 158, 150, 162, 171]
[172, 157, 159, 150, 162, 171]
[172, 158, 159, 150, 162, 171]
[172, 158, 159, 150, 162, 172]
[172, 158, 159, 151, 162, 172]
[172, 158, 160, 151, 162, 172]
[172, 158, 160, 151, 162, 173]
[173, 158, 160, 151, 162, 173]
[173, 158, 160, 151, 163, 173]
[173, 158, 160, 151, 164, 173]
[173, 158, 160, 151, 164, 174]
[173, 159, 160, 151, 164, 174]
[173, 159, 161, 151, 164, 174]
[173, 159, 162, 151, 164, 174]
[174, 159, 162, 151, 164, 174]
[174, 159, 162, 151, 165, 174]
[174, 159, 162, 152, 165, 174]
[174, 159, 163, 152, 165, 174]
[174, 159, 163, 153, 165, 174]
[174, 159, 163, 153, 166, 174]
[174, 159, 163, 153, 166, 175]
[175, 159, 163, 153, 166, 175]
[176, 159, 163, 153, 166, 175]
[176, 159, 163, 154, 166, 175]
[176, 159, 163, 154, 167, 175]
[176, 159, 163, 155, 167, 175]
[176, 159, 163, 156, 167, 175]
[176, 160, 163, 156, 167, 175]
[176, 160, 163, 156, 168, 175]
[176, 160, 163, 156, 169, 175]
[176, 160, 163, 157, 169, 175]
[176, 160, 163, 157, 169, 176]

结果看起来是合理的:我们看到了6个值——掷D6骰子时可能出现的每个点数对应一个;我们还发现,没有任何点数出现的频率比其他点数高很多。下面来可视化这些结果。

15.4.6绘制直方图
有了频率列表后,我们就可以绘制一个表示结果的直方图。直方图 是一种条形图,指出了各种结果出现的频率。创建这种直方图的代码如下:
die_visual.py

from  die  import  Die
import  pygal
# 创建一个D6
die  =  Die()# 掷几次骰子,并将结果存储在一个列表中
results  =  []
for  roll_num  in  range(1001):result  =  die.roll()results.append(result)# 分析结果frequencies  =  []for  value  in  range(1,  die.num_sides+1):frequency  =  results.count(value)frequencies.append(frequency)# 对结果进行可视化hist = pygal.Bar()hist.title = "Results  of  rolling  one  D6  1000  times."hist.x_labels = ['1', '2', '3', '4', '5', '6']hist.x_title = "Result"hist.y_title = "Frequency  of  Result"hist.add('D6', frequencies)hist.render_to_file('die_visual.svg')

为创建条形图,我们创建了一个pygal.Bar() 实例,并将其存储在hist 中(见❶)。接下来,我们设置hist 的属性title (用于标示直方图的字符串),将掷D6骰子的可能结果用作 x 轴的标签(见❷),并给每个轴都添加了标题。在❸处,我们使用add() 将一系列值添加到图表中(向它传递要给添加的值指定的标签,还有一个列表,其中包含将出现在图表中的值)。最后,我们将这个图表渲染为一个SVG文件,这种文件的扩展名必须为.svg。
要查看生成的直方图,最简单的方式是使用Web浏览器。为此,在任何Web浏览器中新建一个标签页,再在其中打开文件die_visual.svg(它位于die_visual.py所在的文件夹中)。你 将看到一个类似于图15-11所示的图表(为方便印刷,我稍微修改了这个图表;默认情况下,Pygal生成的图表的背景比你在图15-11中看到的要暗)。

图15-11 使用Pygal创建的简单条形图

注意,Pygal让这个图表具有交互性:如果你将鼠标指向该图表中的任何条形,将看到与之相关联的数据。在同一个图表中绘制多个数据集时,这项功能显得特别有用。

15.4.7同时掷两个骰子
同时掷两个骰子时,得到的点数更多,结果分布情况也不同。下面来修改前面的代码,创建两个D6骰子,以模拟同时掷两个骰子的情况。每次掷两个骰子时,我们都将两个骰子的点数相加,并将结果存储在results 中。请复制die_visual.py并将其保存为dice_visual.py,再做如下修改:
dice_visual.py

import  pygal
from  die  import  Die # 创建两个D6骰子
die_1  =  Die()
die_2  =  Die()
# 掷骰子多次,并将结果存储到一个列表中
results  =  []
for  roll_num  in  range(1000):result  =  die_1.roll()  +  die_2.roll()results.append(result)# 分析结果frequencies = []max_result = die_1.num_sides + die_2.num_sidesfor value in range(2, max_result + 1): frequency = results.count(value)frequencies.append(frequency# 可视化结果hist = pygal.Bar()hist.title = "Results  of  rolling  two  D6  dice  1000  times."hist.x_labels = ['2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12']hist.x_title = "Result"hist.y_title = "Frequency  of  Result"hist.add('D6  +  D6', frequencies)hist.render_to_file('dice_visual.svg')

创建两个Die 实例后,我们掷骰子多次,并计算每次的总点数(见❶)。可能出现的最大点数12为两个骰子的最大可能点数之和,我们将这个值存储在了max_result 中(见
❷)。可能出现的最小总点数2为两个骰子的最小可能点数之和。分析结果时,我们计算2到max_result 的各种点数出现的次数(见❸)。我们原本可以使用range(2, 13)
,但这只适用于两个D6骰子。模拟现实世界的情形时,最好编写可轻松地模拟各种情形的代码。前面的代码让我们能够模拟掷任何两个骰子的情形,而不管这些骰子有多少面。
创建图表时,我们修改了标题、x 轴标签和数据系列(见❹)。(如果列表x_labels 比这里所示的长得多,那么编写一个循环来自动生成它将更合适。)运行这些代码后,在浏览器中刷新显示图表的标签页,你将看到如图15-12所示的图表。

图15-12 模拟同时掷两个6面骰子1000次的结果

这个图表显示了掷两个D6骰子时得到的大致结果。正如你看到的,总点数为2或12的可能性最小,而总点数为7的可能性最大,这是因为在6种情况下得到的总点数都为7。这6种 情况如下:1和6、2和5、3和4、4和3、5和2、6和1。

15.4.8同时掷两个面数不同的骰子
下面来创建一个6面骰子和一个10面骰子,看看同时掷这两个骰子50 000次的结果如何:
different_dice.py

import  pygal
from  die  import  Die # 创建两个D6骰子
die_1  =  Die()
die_2  =  Die()
# 掷骰子多次,并将结果存储到一个列表中
results  =  []
for  roll_num  in  range(50000):result  =  die_1.roll()  +  die_2.roll()results.append(result)# 分析结果
frequencies = []
max_result = die_1.num_sides + die_2.num_sides
for value in range(2, max_result + 1):frequency = results.count(value)frequencies.append(frequency)# 可视化结果
hist  =  pygal.Bar()
hist.title  =  "Results  of  rolling  a  D6  and  a  D10  50,000  times."
hist.x_labels  =  ['2',  '3',  '4',  '5',  '6',  '7',  '8',  '9',  '10',  '11',  '12','13',  '14',  '15',  '16']hist.x_title  =  "Result"
hist.y_title  =  "Frequency  of  Result"
hist.add('D6  +  D10',  frequencies)
hist.render_to_file('dice_visual.svg')

为创建D10骰子,我们在创建第二个Die 实例时传递了实参10 (见❶)。我们还修改了第一个循环,以模拟掷骰子50 000次而不是1000次。可能出现的最小总点数依然是2,但现在可能出现的最大总点数为16,因此我们相应地调整了标题、x 轴标签和数据系列标签(见❷)。
图15-13显示了最终的图表。可能性最大的点数不是一个,而是5个,这是因为导致出现最小点数和最大点数的组合都只有一种(1和1以及6和10),但面数较小的骰子限制了得到

中间点数的组合数:得到总点数7、8、9、10和11的组合数都是6种。因此,这些总点数是最常见的结果,它们出现的可能性相同。

图15-13 同时掷6面骰子和10面骰子50 000次的结果

通过使用Pygal来模拟掷骰子的结果,能够非常自由地探索这种现象。只需几分钟,就可以掷各种骰子很多次。

15.5小结
在本章中,你学习了:如何生成数据集以及如何对其进行可视化;如何使用matplotlib创建简单的图表,以及如何使用散点图来探索随机漫步过程;如何使用Pygal来创建直方图, 以及如何使用直方图来探索同时掷两个面数不同的骰子的结果。

使用代码生成数据集是一种有趣而强大的方式,可用于模拟和探索现实世界的各种情形。完成后面的数据可视化项目时,请注意可使用代码模拟哪些情形。请研究新闻媒体中的可视化,看看其中是否有图表是以你在这些项目中学到的类似方式生成的。
在第16章中,我们将从网上下载数据,并继续使用matplotlib和Pygal来探索这些数据。

第 15 章 生成数据相关推荐

  1. 《Python编程:从入门到实践》读书笔记:第15章 生成数据

    目录 第15章 生成数据 15.1 安装Matplotlib 15.2 绘制简单的折线图 15.2.1 修改标签文字和线条粗细 15.2.2 校正图形 15.2.3 使用内置样式 15.2.4 使用s ...

  2. 《从Excel到Python——数据分析进阶指南》一第1章 生成数据表

    本节书摘来自异步社区<从Excel到Python--数据分析进阶指南>一书中的第1章,第1.1节,作者 王彦平(蓝鲸),更多章节内容可以访问云栖社区"异步社区"公众号查 ...

  3. 《Python数据分析与挖掘实战》第15章 ——电商产品评论数据情感分析(LED)

    文章目录 1.挖掘背景与目标 2.2 数据探索与预处理 2.1 数据筛选 2.2 数据去重 2.3 删除前缀评分 2.4 jieba分词 3 基于LDA 模型的主题分析 4.权重 5.如何在主题空间比 ...

  4. SAP BW4 410(数据抽取、信息对象、建模等)学习1-5章

    SAP BW4 410学习1-5章 1 Introduction to SAP HANA 1.1 描述SAPHANA的演变和数据布局 1.1.1 Introduction to SAP HANA 1. ...

  5. 《C++程序设计原理与实践》笔记 第15章 绘制函数图和数据图

    本章讨论绘制函数图和数据图的基本机制.关键例子是绘制一元函数图像,以及展示从文件中读取的值. 15.1 引言 我们的主要目标不是输出的美观性,而是理解如何生成这样的图形输出以及所使用的编程技术.你会发 ...

  6. 《从Excel到R 数据分析进阶指南》一第1章 生成数据表1.1 导入数据表

    本节书摘来自异步社区<从Excel到R 数据分析进阶指南>一书中的第1章,第1.1节,作者 王彦平(蓝鲸),更多章节内容可以访问云栖社区"异步社区"公众号查看 第1章 ...

  7. 【2023】Kotlin教程 第二篇 面向对象与函数式编程 第15章 数据容器——数组和集合 15.2 集合概述

    [2023]Kotlin教程 文章目录 [2023]Kotlin教程 第二篇 面向对象与函数式编程 第15章 数据容器--数组和集合 15.2 集合概述 第二篇 面向对象与函数式编程 第15章 数据容 ...

  8. 云计算与大数据第15章 分布式大数据处理平台Hadoop习题带答案

    第15章 分布式大数据处理平台Hadoop习题 15.1 选择题 1.分布式系统的特点不包括以下的(  D  ). A. 分布性     B. 高可用性        C. 可扩展性     D.串行 ...

  9. 《大数据之路:阿里巴巴大数据实践》-第3篇 数据管理篇 -第15章 数据质量

    <大数据之路:阿里巴巴大数据实践>系列丛书  第1章 总述 第1篇 数据技术篇  第2章 日志釆集  第3章 数据同步  第4章 离线数据开发  第5章 实时技术  第6章 数据服务  第 ...

最新文章

  1. 清华构建新一代数据集NICO,定义图像分类新标准
  2. 手机存储sqlite还是mysql_手机存储 和sd 路径
  3. Python基本语法_函数_参数的多类型传值
  4. mysql 存储过程 无限分类_mysql存储过程实现无限分类
  5. js函数 every some map ()=a+b;
  6. Java黑皮书课后题第6章:*6.20(计算一个字符串中字母的个数)编写一个方法,使用下面的方法体计算字符串中的字母个数。编写一个测试程序,提示用户输入字符串,然后显示字符串中的字母个数
  7. 类的6个默认成员函数:构造函数、析构函数、拷贝构造函数、重载运算符、三/五法则
  8. fastJson toJSONString注意点
  9. 将一副图片编译进uboot
  10. RestTemplate与Feign使用对比
  11. php验证码图片看不清更换代码,php如何实现验证码看不清换一张的效果
  12. 计量数据分析数据库-计量分析资料大全空间计量分析资料大全
  13. word2019每页设置不同页眉
  14. 高登学苑-解密携程4.9分的秘诀学习笔记
  15. 210 Ruby 方法、类和模块【Rails后端开发训练营】
  16. Python闭包理解
  17. 聚类算法(4)--Hierarchical clustering层次聚类
  18. PDF怎么编辑注释? 5 个免费好用的 PDF 注释器
  19. 15.React-router6的编程式导航
  20. 项目管理笔记-第十三章 项目相关方管理

热门文章

  1. idea打jar包方式
  2. 软件设计师 软考 真题练习 (四)
  3. python基因差异分析_玉米RNA-seq测序数据差异基因分析
  4. 教你用OCR文字识别工具快速识别图片上的文字
  5. Uncaught SyntaxError: Cannot use import statement outside a module
  6. Java时间类型转换详细代码示例
  7. 线程间通信的几种方式
  8. 青灯教育python视频教程_PHD的生涯
  9. 关于密码和商用密码应用
  10. 蓝桥杯物联网竞赛基础图文教程——时钟选择