matplotlib 是python最著名的绘图库,它提供了一整套和matlab相似的命令API,十分适合交互式地进行制图。而且也可以方便地将它作为绘图控件,嵌入GUI应用程序中。
它的文档相当完备,并且 Gallery页面 中有上百幅缩略图,打开之后都有源程序。因此如果你需要绘制某种类型的图,只需要在这个页面中浏览/复制/粘贴一下,基本上都能搞定。
本章节作为matplotlib的入门介绍,将较为深入地挖掘几个例子,从中理解和学习matplotlib绘图的一些基本概念。
5.1 快速绘图
matplotlib的pyplot子库提供了和matlab类似的绘图API,方便用户快速绘制2D图表。让我们先来看一个简单的例子:
# -*- coding: utf-8 -*-import numpy as npimport ma as pltx = np.linspace(0, 10, 1000)y = np.sin(x)z = np.cos(x**2(figsize=(8,4(x,y,label="$sin(x)$",color="red",linewidth=2(x,z,"b--",label="$cos(x^2)$"("Time(s)"("Volt"("PyPlot First Example"()
图5.1 调用pyplot库快速将数据绘制成曲线图
matplotlib中的快速绘图的函数库可以通过如下语句载入:
import ma as pltpylab模块
matplotlib还提供了名为pylab的模块,其中包括了许多numpy和pyplot中常用的函数,方便用户快速进行计算和绘图,可以用于IPython中的快速交互式使用。
接下来调用figure创建一个绘图对象,并且使它成为当前的绘图对象。
(figsize=(8,4))也可以不创建绘图对象直接调用接下来的plot函数直接绘图,matplotlib会为我们自动创建一个绘图对象。如果需要同时绘制多幅图表的话,可以是给figure传递一个整数参数指定图标的序号,如果所指定序号的绘图对象已经存在的话,将不创建新的对象,而只是让它成为当前绘图对象。
通过figsize参数可以指定绘图对象的宽度和高度,单位为英寸;dpi参数指定绘图对象的分辨率,即每英寸多少个像素,缺省值为80。因此本例中所创建的图表窗口的宽度为8*80 = 640像素。
但是用工具栏中的保存按钮保存下来的png图像的大小是800*400像素。这是因为保存图表用的函数savefig使用不同的DPI配置,savefig函数也有一个dpi参数,如果不设置的话,将使用matplotlib配置文件中的配置,此配置可以通过如下语句进行查看,关于配置文件将在后面的章节进行介绍:
>>> import matplotlib>>> ma[";]100下面的两行程序通过调用plot函数在当前的绘图对象中进行绘图:
(x,y,label="$sin(x)$",color="red",linewidth=2(x,z,"b--",label="$cos(x^2)$")plot函数的调用方式很灵活,第一句将x,y数组传递给plot之后,用关键字参数指定各种属性:
label : 给所绘制的曲线一个名字,此名字在图示(legend)中显示。只要在字符串前后添加"$"符号,matplotlib就会使用其内嵌的latex引擎绘制的数学公式。
color : 指定曲线的颜色
linewidth : 指定曲线的宽度
第二句直接通过第三个参数"b--"指定曲线的颜色和线型,这个参数称为格式化参数,它能够通过一些易记的符号快速指定曲线的样式。其中b表示蓝色,"--"表示线型为虚线。在IPython中输入 "?" 可以查看格式化字符串的详细配置。
接下来通过一系列函数设置绘图对象的各个属性:
("Time(s)"("Volt"("PyPlot First Example"()xlabel : 设置X轴的文字
ylabel : 设置Y轴的文字
title : 设置图表的标题
ylim : 设置Y轴的范围
legend : 显示图示
最后调用()显示出我们创建的所有绘图对象。
5.1.1 配置属性
matplotlib所绘制的图的每个组成部分都对应有一个对象,我们可以通过调用这些对象的属性设置方法set_*或者pyplot的属性设置函数setp设置其属性值。例如plot函数返回一个 ma 对象的列表,下面的例子显示如何设置Line2D对象的属性:
>>> import numpy as np>>> import ma as plt>>> x = np.arange(0, 5, 0.1)>>> line, = (x, x*x) # plot返回一个列表,通过line,获取其第一个元素>>> # 调用Line2D对象的set_*方法设置属性值>>> line.set_antialiased(False)>>> # 同时绘制sin和cos两条曲线,lines是一个有两个Line2D对象的列表>>> lines = (x, np.sin(x), x, np.cos(x)) #>>> # 调用setp函数同时配置多个Line2D对象的多个属性值>>> (lines, color="r", linewidth=2.0)这段例子中,通过调用Line2D对象line的set_antialiased方法,关闭对象的反锯齿效果。或者通过调用函数配置多个Line2D对象的颜色和线宽属性。
同样我们可以通过调用Line2D对象的get_*方法,或者函数获取对象的属性值:
>>> line.get_linewid;>> (lines[0], "color") # 返回color属性'r'>>> (lines[1]) # 输出全部属性alpha = 1.0animated = Falseantialiased or aa = Trueaxes = Axe;0.775x0.8)... ...注意getp函数只能对一个对象进行操作,它有两种用法:
指定属性名:返回对象的指定属性的值
不指定属性名:打印出对象的所有属性和其值
matplotlib的整个图表为一个Figure对象,此对象在调用函数时返回,我们也可以通过函数获取当前的绘图对象:
>>> f = ()>>> (f)alpha = 1.0animated = False...Figure对象有一个axes属性,其值为AxesSubplot对象的列表,每个AxesSubplot对象代表图表中的一个子图,前面所绘制的图表只包含一个子图,当前子图也可以通过获得:
>>> (f, "axes")[<ma object at 0x05CDD170>]>>> ()<ma object at 0x05CDD170>用可以发现AxesSubplot对象有很多属性,例如它的lines属性为此子图所包括的 Line2D 对象列表:
>>> alllines = ((), "lines")>>> alllines<a list of 3 Line2D objects>>>> alllines[0] == line # 其中的第一条曲线就是最开始绘制的那条曲线True通过这种方法我们可以很容易地查看对象的属性和它们之间的包含关系,找到需要配置的属性。
5.2 绘制多轴图
一个绘图对象(figure)可以包含多个轴(axis),在Matplotlib中用轴表示一个绘图区域,可以将其理解为子图。上面的第一个例子中,绘图对象只包括一个轴,因此只显示了一个轴(子图)。我们可以使用subplot函数快速绘制有多个轴的图表。subplot函数的调用形式如下:
subplot(numRows, numCols, plotNum)subplot将整个绘图区域等分为numRows行 * numCols列个子区域,然后按照从左到右,从上到下的顺序对每个子区域进行编号,左上的子区域的编号为1。如果numRows,numCols和plotNum这三个数都小于10的话,可以把它们缩写为一个整数,例如subplot(323)和subplot(3,2,3)是相同的。subplot在plotNum指定的区域中创建一个轴对象。如果新创建的轴和之前创建的轴重叠的话,之前的轴将被删除。
下面的程序创建3行2列共6个轴,通过axisbg参数给每个轴设置不同的背景颜色。
for idx, color in enumerate("rgbyck"): (320+idx+1, axisbg=color)()
图5.2 用subplot函数将Figure分为六个子图区域
如果希望某个轴占据整个行或者列的话,可以如下调用subplot:
(221) # 第一行的左图(222) # 第一行的右图(212) # 第二整行()
图5.3 将Figure分为三个子图区域
当绘图对象中有多个轴的时候,可以通过工具栏中的Configure Subplots按钮,交互式地调节轴之间的间距和轴与边框之间的距离。如果希望在程序中调节的话,可以调用subplots_adjust函数,它有left, right, bottom, top, wspace, hspace等几个关键字参数,这些参数的值都是0到1之间的小数,它们是以绘图区域的宽高为1进行正规化之后的坐标或者长度。
5.3 配置文件
一幅图有许多需要配置的属性,例如颜色、字体、线型等等。我们在绘图时,并没有一一对这些属性进行配置,许多都直接采用了Matplotlib的缺省配置。Matplotlib将缺省配置保存在一个文件中,通过更改这个文件,我们可以修改这些属性的缺省值。
Matplotlib 使用配置文件 matplotlibrc 时的搜索顺序如下:
当前路径 : 程序的当前路径
用户配置路径 : 通常为 HOME/.matplotlib/,可以通过环境变量MATPLOTLIBRC修改
系统配置路径 : 保存在 matplotlib的安装目录下的 mpl-data 下
通过下面的语句可以获取用户配置路径:
>>> import matplotlib>>> ma()'C:\\Documents and Settings\\zhang\\.matplotlib'通过下面的语句可以获得目前使用的配置文件的路径:
>>> import matplotlib>>> ma()'C:\\Python26\\lib\\site-packages\\matplotlib\\mpl-data\\matplotlibrc'由于在当前路径和用户配置路径中都没有找到位置文件,因此最后使用的是系统配置路径下的配置文件。如果你将matplotlibrc复制一份到脚本的当前目录下:
>>> import os>>> os.getcwd()'C:\\zhang\\doc'复制配置文件之后再运行:
>>> ma()'C:\\zhang\\doc\\matplotlibrc'如果你用文本编辑器打开此配置文件的话,你会发现它实际上是定义了一个字典。为了对众多的配置进行区分,关键字可以用点分开。
配置文件的读入可以使用 rc_params 函数,它返回一个配置字典:
>>> ma(){'agg.;: 0, 'axes.axisbelow': False, 'axes.edgecolor': 'k', 'axes.facecolor': 'w', ... ...在matplotlib模块载入的时候会调用rc_params,并把得到的配置字典保存到rcParams变量中:
>>> ma{'agg.;: 0,'axes.axisbelow': False,... ...matplotlib将使用rcParams中的配置进行绘图。用户可以直接修改此字典中的配置,所做的改变会反映到此后所绘制的图中。例如下面的脚本所绘制的线将带有圆形的点标识符:
>>> ma["lines.marker"] = "o">>> import pylab>>> ([1,2,3])>>> ()为了方便配置,可以使用rc函数,下面的例子同时配置点标识符、线宽和颜色:
>>> ma("lines", marker="x", linewidth=2, color="red")如果希望恢复到缺省的配置(matplotlib载入时从配置文件读入的配置)的话,可以调用 rcdefaults 函数。
>>> madefaults()如果手工修改了配置文件,希望重新从配置文件载入最新的配置的话,可以调用:
>>> ma.update( ma() )5.4 Artist对象
matplotlib API包含有三层:
backend_ba : 图表的绘制领域
backend_ba : 知道如何在FigureCanvas上如何绘图
ar : 知道如何使用Renderer在FigureCanvas上绘图
FigureCanvas和Renderer需要处理底层的绘图操作,例如使用wxPython在界面上绘图,或者使用PostScript绘制PDF。Artist则处理所有的高层结构,例如处理图表、文字和曲线等的绘制和布局。通常我们只和Artist打交道,而不需要关心底层的绘制细节。
Artists分为简单类型和容器类型两种。简单类型的Artists为标准的绘图元件,例如Line2D、 rectangle、 Text、AxesImage 等等。而容器类型则可以包含许多简单类型的Artists,使它们组织成一个整体,例如Axis、 Axes、Figure等。
直接使用Artists创建图表的标准流程如下:
创建Figure对象
用Figure对象创建一个或者多个Axes或者Subplot对象
调用Axies等对象的方法创建各种简单类型的Artists
下面首先调用辅助函数创建Figure对象,然后调用Figure对象的add_axes方法在其中创建一个Axes对象,add_axes的参数是一个形如[left, bottom, width, height]的列表,这些数值分别指定所创建的Axes对象相对于fig的位置和大小,取值范围都在0到1之间:
>>> import ma as plt>>> fig = ()>>> ax = ([0.15, 0.1, 0.7, 0.3])然后我们调用ax的plot方法绘图,创建一条曲线,并且返回此曲线对象(Line2D)。
>>> line, = ax.plot([1,2,3],[1,2,1])>>> ax.lines[<ma object at 0x0637A3D0>]>>> line<ma object at 0x0637A3D0>ax.lines是一个为包含ax的所有曲线的列表,后续的ax.plot调用会往此列表中添加新的曲线。如果想删除某条曲线的话,直接从此列表中删除即可。
Axes对象还包括许多其它的Artists对象,例如我们可以通过调用set_xlabel设置其X轴上的标题:
>>> ax.set_xlabel("time")如果我们查看set_xlabel的源代码的话,会发现它是通过调用下面的语句实现的:
(xlabel)如果我们一直跟踪下去,会发现Axes的xaxis属性是一个XAxis对象:
>>> ax.xaxis<ma object at 0x06343230>XAxis的label属性是一个Text对象:
>>> ax.xaxis.label<ma object at 0x06343290>而Text对象的_text属性为我们设置的值:
>>> ax.xaxis.label._text'time'这些对象都是Artists,因此也可以调用它们的属性获取函数来获得相应的属性:
>>> ax.xaxis.label.get_text()'time'5.4.1 Artist的属性
图表中的每个元素都用一个matplotlib的Artist对象表示,而每个Artist对象都有一大堆属性控制其显示效果。例如Figure对象和Axes对象都有patch属性作为其背景,它的值是一个Rectangle对象。通过设置此它的一些属性可以修改Figrue图表的背景颜色或者透明度等属性,下面的例子将图表的背景颜色设置为绿色:
>>> fig = ()>>> ()>>> ("g")>>> ()patch的color属性通过set_color函数进行设置,属性修改之后并不会立即反映到图表的显示上,还需要调用()函数才能够更新显示。
下面是Artist对象都具有的一些属性:
alpha : 透明度,值在0到1之间,0为完全透明,1为完全不透明
animated : 布尔值,在绘制动画效果时使用
axes : 此Artist对象所在的Axes对象,可能为None
clip_box : 对象的裁剪框
clip_on : 是否裁剪
clip_path : 裁剪的路径
contains : 判断指定点是否在对象上的函数
figure : 所在的Figure对象,可能为None
label : 文本标签
picker : 控制Artist对象选取
transform : 控制偏移旋转
visible : 是否可见
zorder : 控制绘图顺序
Artist对象的所有属性都通过相应的 get_* 和 set_* 函数进行读写,例如下面的语句将alpha属性设置为当前值的一半:
>>> *())如果你想用一条语句设置多个属性的话,可以使用set函数:
>>> (alpha=0.5, zorder=2)使用前面介绍的 ma.getp 函数可以方便地输出Artist对象的所有属性名和值。
>>> ) aa = True alpha = 1.0 animated = False antialiased or aa = True ... ...5.4.2 Figure容器
现在我们知道如何观察和修改已知的某个Artist对象的属性,接下来要解决如何找到指定的Artist对象。前面我们介绍过Artist对象有容器类型和简单类型两种,这一节让我们来详细看看容器类型的内容。
最大的Artist容器是ma,它包括组成图表的所有元素。图表的背景是一个Rectangle对象,用Figure.patch属性表示。当你通过调用add_subplot或者add_axes方法往图表中添加轴(子图时),这些子图都将添加到Figure.axes属性中,同时这两个方法也返回添加进axes属性的对象,注意返回值的类型有所不同,实际上AxesSubplot是Axes的子类。
>>> fig = ()>>> ax1 = (211)>>> ax2 = ([0.1, 0.1, 0.7, 0.3])>>> ax1<ma object at 0x056BCA90>>>> ax2<ma object at 0x056BC910>>>> [<ma object at 0x056BCA90>,<ma object at 0x056BC910>]为了支持pylab中的gca()等函数,Figure对象内部保存有当前轴的信息,因此不建议直接对Figure.axes属性进行列表操作,而应该使用add_subplot, add_axes, delaxes等方法进行添加和删除操作。但是使用for循环对axes中的每个元素进行操作是没有问题的,下面的语句打开所有子图的栅格。
>>> for ax in : ax.grid(True)Figure对象可以拥有自己的文字、线条以及图像等简单类型的Artist。缺省的坐标系统为像素点,但是可以通过设置Artist对象的transform属性修改坐标系的转换方式。最常用的Figure对象的坐标系是以左下角为坐标原点(0,0),右上角为坐标(1,1)。下面的程序创建并添加两条直线到fig中:
>>> from ma import Line2D>>> fig = ()>>> line1 = Line2D([0,1],[0,1], transform=, figure=fig, color="r")>>> line2 = Line2D([0,1],[1,0], transform=, figure=fig, color="g")>>> ([line1, line2])>>> ()
图5.4 在Figure对象中手工绘制直线
注意为了让所创建的Line2D对象使用fig的坐标,我们将赋给Line2D对象的transform属性;为了让Line2D对象知道它是在fig对象中,我们还设置其figure属性为fig;最后还需要将创建的两个Line2D对象添加到属性中去。
Figure对象有如下属性包含其它的Artist对象:
axes : Axes对象列表
patch : 作为背景的Rectangle对象
images : FigureImage对象列表,用来显示图片
legends : Legend对象列表
lines : Line2D对象列表
patches : patch对象列表
texts : Text对象列表,用来显示文字
5.4.3 Axes容器
Axes容器是整个matplotlib库的核心,它包含了组成图表的众多Artist对象,并且有许多方法函数帮助我们创建、修改这些对象。和Figure一样,它有一个patch属性作为背景,当它是笛卡尔坐标时,patch属性是一个Rectangle对象,而当它是极坐标时,patch属性则是Circle对象。例如下面的语句设置Axes对象的背景颜色为绿色:
>>> fig = ()>>> ax = (111)>>> ax.("green")当你调用Axes的绘图方法(例如plot),它将创建一组Line2D对象,并将所有的关键字参数传递给这些Line2D对象,并将它们添加进Axes.lines属性中,最后返回所创建的Line2D对象列表:
>>> x, y = np.random.rand(2, 100)>>> line, = ax.plot(x, y, "-", color="blue", linewidth=2)>>> line<ma object at 0x03007030>>>> ax.lines[<ma object at 0x03007030>]注意plot返回的是一个Line2D对象的列表,因为我们可以传递多组X,Y轴的数据,一次绘制多条曲线。
与plot方法类似,绘制直方图的方法bar和绘制柱状统计图的方法hist将创建一个Patch对象的列表,每个元素实际上都是Patch的子类Rectangle,并且将所创建的Patch对象都添加进Axes.patches属性中:
>>> ax = (111)>>> n, bins, rects = ax.hi(1000), 50, facecolor="blue")>>> rects<a list of 50 Patch objects>>>> rects[0]<ma object at 0x05BC2350>>>> ax.patches[0]<ma object at 0x05BC2350>一般我们不会直接对Axes.lines或者Axes.patches属性进行操作,而是调用add_line或者add_patch等方法,这些方法帮助我们完成许多属性设置工作:
>>> fig = ()>>> ax = (111)>>> rect = ma((1,1), width=5, height=12)>>> print rect.get_axes() # rect的axes属性为空None>>> rect.get_transform() # rect的transform属性为缺省值BboxTransformTo(Bbox(array([[ 1., 1.], [ 6., 13.]])))>>> ax.add_patch(rect) # 将rect添加进ax<ma object at 0x05C34E50>>>> rect.get_axes() # 于是rect的axes属性就是ax<ma object at 0x05C09CB0>>>> # rect的transform属性和ax的transData相同>>> rect.get_transform()... # 太长,省略>>> ax.transData... # 太长,省略>>> ax.get_xlim() # ax的X轴范围为0到1,无法显示完整的rec, 1.0)>>> ax.da() # 数据的范围和rect的大小一致, 1.0, 5.0, 12.0)>>> ax.autoscale_view() # 自动调整坐标轴范围>>> ax.get_xlim() # 于是X轴可以完整显示rect, 6.0)>>> ()
通过上面的例子我们可以看出,add_patch方法帮助我们设置了rect的axes和transform属性。
下面详细列出Axes包含各种Artist对象的属性:
artists : Artist对象列表
patch : 作为Axes背景的Patch对象,可以是Rectangle或者Circle
collections : Collection对象列表
images : AxesImage对象列表
legends : Legend对象列表
lines : Line2D对象列表
patches : Patch对象列表
texts : Text对象列表
xaxis : XAxis对象
yaxis : YAxis对象
下面列出Axes的创建Artist对象的方法:
Axes的方法 | 所创建的对象 | 添加进的列表 |
annotate | Annotate | texts |
bars | Rectangle | patches |
errorbar | Line2D, Rectangle | lines,patches |
fill | Polygon | patches |
hist | Rectangle | patches |
imshow | AxesImage | images |
legend | Legend | legends |
plot | Line2D | lines |
scatter | PolygonCollection | Collections |
text | Text | texts |
下面以绘制散列图(scatter)为例,验证一下:
>>> fig = ()>>> ax = (111)>>> t = ax.scatter(20), np.random.rand(20))>>> t # 返回值为CircleCollection对象<ma object at 0x06004230>>>> ax.collections # 返回的对象已经添加进了collections列表中[<ma object at 0x06004230>]>>> ()>>> t.get_sizes() # 获得Collection的点数20
图5.5 用scatter函数绘制散列图
5.4.4 Axis容器
Axis容器包括坐标轴上的刻度线、刻度文本、坐标网格以及坐标轴标题等内容。刻度包括主刻度和副刻度,分别通过Axis.get_major_ticks和Axis.get_minor_ticks方法获得。每个刻度线都是一个XTick或者YTick对象,它包括实际的刻度线和刻度文本。为了方便访问刻度线和文本,Axis对象提供了get_ticklabels和get_ticklines方法分别直接获得刻度线和刻度文本:
>>> ([1,2,3],[4,5,6])[<ma object at 0x0AD3B670>]>>> ()>>> axis = ().xaxis>>> axis.get_ticklocs() # 获得刻度的位置列表array([ 1. , 1.5, 2. , 2.5, 3. ])>>> axis.get_ticklabels() # 获得刻度标签列表<a list of 5 Text major ticklabel objects>>>> [x.get_text() for x in axis.get_ticklabels()] # 获得刻度的文本字符串[u'1.0', u'1.5', u'2.0', u'2.5', u'3.0']>>> axis.get_ticklines() # 获得主刻度线列表,图的上下刻度线共10条<a list of 10 Line2D ticklines objects>>>> axis.get_ticklines(minor=True) # 获得副刻度线列表<a list of 0 Line2D ticklines objects>
获得刻度线或者刻度标签之后,可以设置其各种属性,下面设置刻度线为绿色粗线,文本为红色并且旋转45度:
>>> for label in axis.get_ticklabels():... label.set_color("red")... label.set_rotation(45)... label.set_fontsize(16)...>>> for line in axis.get_ticklines():... line.set_color("green")... line.set_markersize(25)... line.set_markeredgewidth(3)
最终的结果图如下:
图5.6 手工配置X轴的刻度线和刻度文本的样式
上面的例子中,获得的副刻度线列表为空,这是因为用于计算副刻度的对象缺省为NullLocator,它不产生任何刻度线;而计算主刻度的对象为AutoLocator,它会根据当前的缩放等配置自动计算刻度的位置:
>>> axis.get_minor_locator() # 计算副刻度的对象<ma instance at 0x0A014300>>>> axis.get_major_locator() # 计算主刻度的对象<ma instance at 0x09281B20>我们可以使用程序为Axis对象设置不同的Locator对象,用来手工设置刻度的位置;设置Formatter对象用来控制刻度文本的显示。下面的程序设置X轴的主刻度为pi/4,副刻度为pi/20,并且主刻度上的文本以pi为单位:
# -*- coding: utf-8 -*-import ma as plfrom ma import MultipleLocator, FuncFormatterimport numpy as npx = np.arange(0, 4*np.pi, 0.01)y = np.sin(x(figsize=(8,4))(x, y)ax = ()def pi_formatter(x, pos): """ 比较罗嗦地将数值转换为以pi/4为单位的刻度文本 """ m = np.round(x / )) n = 4 if m%2==0: m, n = m/2, n/2 if m%2==0: m, n = m/2, n/2 if m == 0: return "0" if m == 1 and n == 1: return "$\pi$" if n == 1: return r"$%d \pi$" % m if m == 1: return r"$\frac{\pi}{%d}$" % n return r"$\frac{%d \pi}{%d}$" % (m,n)# 设置两个坐标轴的范围(0, np.max(x))# 设置图的底边距(bottom = 0.15() #开启网格# 主刻度为pi( MultipleLocator) )# 主刻度文本用pi_formatter函数计算ax.xaxis.set_major_formatter( FuncFormatter( pi_formatter ) )# 副刻度为pi( MultipleLocator) )# 设置刻度文本的大小for tick in ax.xaxis.get_major_ticks(): (16)()关于刻度的定位和文本格式的东西都在ma中定义,程序中使用到如下两个类:
MultipleLocator : 以指定值的整数倍为刻度放置刻度线
FuncFormatter : 使用指定的函数计算刻度文本,他会传递给所指定的函数两个参数:刻度值和刻度序号,程序中通过比较笨的办法计算出刻度值所对应的刻度文本
此外还有很多预定义的Locator和Formatter类,详细内容请参考相应的API文档。