matplotlib 是python最著名的绘图库,它提供了一整套和matlab类似的命令API,十分适合交互式地行制图。并且也能够方便地将它做为绘图控件,嵌入GUI应用程序中。html
它的文档至关完备,而且Gallery页面中有上百幅缩略图,打开以后都有源程序。所以若是你须要绘制某种类型的图,只须要在这个页面中浏览/复制/粘贴一下,基本上都能搞定。python
在Linux下比较著名的数据图工具还有gnuplot,这个是免费的,Python有一个包能够调用gnuplot,可是语法比较不习惯,并且画图质量不高。express
而 Matplotlib则比较强:Matlab的语法、python语言、latex的画图质量(还可使用内嵌的latex引擎绘制的数学公式)。canvas
Matplotlib.pyplot快速绘图数组
matplotlib其实是一套面向对象的绘图库,它所绘制的图表中的每一个绘图元素,例如线条Line2D、文字Text、刻度等在内存中都有一个对象与之对应。ide
为了方便快速绘图matplotlib经过pyplot模块提供了一套和MATLAB相似的绘图API,将众多绘图对象所构成的复杂结构隐藏在这套API内部。咱们只须要调用pyplot模块所提供的函数就能够实现快速绘图以及设置图表的各类细节。pyplot模块虽然用法简单,但不适合在较大的应用程序中使用。函数
为了将面向对象的绘图库包装成只使用函数的调用接口,pyplot模块的内部保存了当前图表以及当前子图等信息。当前的图表和子图可使用plt.gcf()和plt.gca()得到,分别表示"Get Current Figure"和"Get Current Axes"。在pyplot模块中,许多函数都是对当前的Figure或Axes对象进行处理,好比说:工具
plt.plot()实际上会经过plt.gca()得到当前的Axes对象ax,而后再调用ax.plot()方法实现真正的绘图。oop
能够在Ipython中输入相似"plt.plot??"的命令查看pyplot模块的函数是如何对各类绘图对象进行包装的。
matplotlib所绘制的图表的每一个组成部分都和一个对象对应,咱们能够经过调用这些对象的属性设置方法set_*()或者pyplot模块的属性设置函数setp()设置它们的属性值。
由于matplotlib其实是一套面向对象的绘图库,所以也能够直接获取对象的属性
绘制一幅图须要对许多对象的属性进行配置,例如颜色、字体、线型等等。咱们在绘图时,并无逐一对这些属性进行配置,许多都直接采用了matplotlib的缺省配置。
matplotlib将这些缺省配置保存在一个名为“matplotlibrc”的配置文件中,经过修改配置文件,咱们能够修改图表的缺省样式。配置文件的读入可使用rc_params(),它返回一个配置字典;在matplotlib模块载入时会调用rc_params(),并把获得的配置字典保存到rcParams变量中;matplotlib将使用rcParams字典中的配置进行绘图;用户能够直接修改此字典中的配置,所作的改变会反映到此后建立的绘图元素。
绘制多子图(快速绘图)
Matplotlib 里的经常使用类的包含关系为 Figure -> Axes -> (Line2D, Text, etc.)
一个Figure对象能够包含多个子图(Axes),在matplotlib中用Axes对象表示一个绘图区域,能够理解为子图。
可使用subplot()快速绘制包含多个子图的图表,它的调用形式以下:
subplot(numRows, numCols, plotNum)
subplot将整个绘图区域等分为numRows行* numCols列个子区域,而后按照从左到右,从上到下的顺序对每一个子区域进行编号,左上的子区域的编号为1。若是numRows,numCols和plotNum这三个数都小于10的话,能够把它们缩写为一个整数,例如subplot(323)和subplot(3,2,3)是相同的。subplot在plotNum指定的区域中建立一个轴对象。若是新建立的轴和以前建立的轴重叠的话,以前的轴将被删除。
subplot()返回它所建立的Axes对象,咱们能够将它用变量保存起来,而后用sca()交替让它们成为当前Axes对象,并调用plot()在其中绘图。
绘制多图表(快速绘图)
若是须要同时绘制多幅图表,能够给figure()传递一个整数参数指定Figure对象的序号,若是序号所指定的Figure对象已经存在,将不建立新的对象,而只是让它成为当前的Figure对象。
import numpy as npimport matplotlib.pyplot as pltplt.figure(1) # 建立图表1plt.figure(2) # 建立图表2ax1 = plt.subplot(211) # 在图表2中建立子图1ax2 = plt.subplot(212) # 在图表2中建立子图2x = np.linspace(0, 3, 100)for i in xrange(5): plt.figure(1) #❶ # 选择图表1 plt.plot(x, np.exp(i*x/3)) plt.sca(ax1) #❷ # 选择图表2的子图1 plt.plot(x, np.sin(i*x)) plt.sca(ax2) # 选择图表2的子图2 plt.plot(x, np.cos(i*x))plt.show()
matplotlib的缺省配置文件中所使用的字体没法正确显示中文。为了让图表能正确显示中文,能够有几种解决方案。
在程序中直接指定字体。
在程序开头修改配置字典rcParams。
修改配置文件。
面向对象画图
matplotlib API包含有三层,Artist层处理全部的高层结构,例如处理图表、文字和曲线等的绘制和布局。一般咱们只和Artist打交道,而不须要关心底层的绘制细节。
直接使用Artists建立图表的标准流程以下:
建立Figure对象
用Figure对象建立一个或者多个Axes或者Subplot对象
调用Axies等对象的方法建立各类简单类型的Artists
import matplotlib.pyplot as pltX1 = range(0, 50)Y1 = [num**2 for num in X1] # y = x^2X2 = [0, 1]Y2 = [0, 1] # y = xFig = plt.figure(figsize=(8,4)) # Create a `figure' instanceAx = Fig.add_subplot(111) # Create a `axes' instance in the figureAx.plot(X1, Y1, X2, Y2) # Create a Line2D instance in the axesFig.show()Fig.savefig("test.pdf")
参考:
《Python科学计算》(Numpy视频) matplotlib-绘制精美的图表(快速绘图)(面向对象绘图)(深刻浅出适合系统学习)
什么是 Matplotlib (主要讲面向对象绘图,对于新手可能有点乱)
Matplotlib.pylab快速绘图
matplotlib还提供了一个名为pylab的模块,其中包括了许多NumPy和pyplot模块中经常使用的函数,方便用户快速进行计算和绘图,十分适合在IPython交互式环境中使用。这里使用下面的方式载入pylab模块:
>>> import pylab as pl
1 安装numpy和matplotlib
>>> import numpy
>>> numpy.__version__
>>> import matplotlib
>>> matplotlib.__version__
2 两种经常使用图类型:Line and scatter plots(使用plot()命令), histogram(使用hist()命令)
2.1 折线图&散点图 Line and scatter plots
2.1.1 折线图 Line plots(关联一组x和y值的直线)
import numpy as npimport pylab as plx = [1, 2, 3, 4, 5]# Make an array of x valuesy = [1, 4, 9, 16, 25]# Make an array of y values for each x valuepl.plot(x, y)# use pylab to plot x and ypl.show()# show the plot on the screen
2.1.2 散点图 Scatter plots
把pl.plot(x, y)改为pl.plot(x, y, 'o')便可,下图的蓝色版本
2.2 美化 Making things look pretty
2.2.1 线条颜色 Changing the line color
红色:把pl.plot(x, y, 'o')改为pl.plot(x, y, ’or’)
2.2.2 线条样式 Changing the line style
虚线:plot(x,y, '--')
2.2.3 marker样式 Changing the marker style
蓝色星型markers:plot(x,y, ’b*’)
2.2.4 图和轴标题以及轴坐标限度 Plot and axis titles and limits
import numpy as npimport pylab as plx = [1, 2, 3, 4, 5]# Make an array of x valuesy = [1, 4, 9, 16, 25]# Make an array of y values for each x valuepl.plot(x, y)# use pylab to plot x and ypl.title(’Plot of y vs. x’)# give plot a titlepl.xlabel(’x axis’)# make axis labelspl.ylabel(’y axis’)pl.xlim(0.0, 7.0)# set axis limitspl.ylim(0.0, 30.)pl.show()# show the plot on the screen
2.2.5 在一个坐标系上绘制多个图 Plotting more than one plot on the same set of axes
作法是很直接的,依次做图便可:
import numpy as npimport pylab as plx1 = [1, 2, 3, 4, 5]# Make x, y arrays for each graphy1 = [1, 4, 9, 16, 25]x2 = [1, 2, 4, 6, 8]y2 = [2, 4, 8, 12, 16]pl.plot(x1, y1, ’r’)# use pylab to plot x and ypl.plot(x2, y2, ’g’)pl.title(’Plot of y vs. x’)# give plot a titlepl.xlabel(’x axis’)# make axis labelspl.ylabel(’y axis’)pl.xlim(0.0, 9.0)# set axis limitspl.ylim(0.0, 30.)pl.show()# show the plot on the screen
2.2.6 图例 Figure legends
pl.legend((plot1, plot2), (’label1, label2’), 'best’, numpoints=1)
其中第三个参数表示图例放置的位置:'best’‘upper right’, ‘upper left’, ‘center’, ‘lower left’, ‘lower right’.
若是在当前figure里plot的时候已经指定了label,如plt.plot(x,z,label="$cos(x^2)$"),直接调用plt.legend()就能够了哦。
import numpy as npimport pylab as plx1 = [1, 2, 3, 4, 5]# Make x, y arrays for each graphy1 = [1, 4, 9, 16, 25]x2 = [1, 2, 4, 6, 8]y2 = [2, 4, 8, 12, 16]plot1 = pl.plot(x1, y1, ’r’)# use pylab to plot x and y : Give your plots namesplot2 = pl.plot(x2, y2, ’go’)pl.title(’Plot of y vs. x’)# give plot a titlepl.xlabel(’x axis’)# make axis labelspl.ylabel(’y axis’)pl.xlim(0.0, 9.0)# set axis limitspl.ylim(0.0, 30.)pl.legend([plot1, plot2], (’red line’, ’green circles’), ’best’, numpoints=1)# make legendpl.show()# show the plot on the screen
2.3 直方图 Histograms
import numpy as npimport pylab as pl# make an array of random numbers with a gaussian distribution with# mean = 5.0# rms = 3.0# number of points = 1000data = np.random.normal(5.0, 3.0, 1000)# make a histogram of the data arraypl.hist(data)# make plot labelspl.xlabel(’data’)pl.show()
若是不想要黑色轮廓能够改成pl.hist(data, histtype=’stepfilled’)
2.3.1 自定义直方图bin宽度 Setting the width of the histogram bins manually
增长这两行
bins = np.arange(-5., 16., 1.) #浮点数版本的range
pl.hist(data, bins, histtype=’stepfilled’)
3 同一画板上绘制多幅子图 Plotting more than one axis per canvas
若是须要同时绘制多幅图表的话,能够是给figure传递一个整数参数指定图标的序号,若是所指定
序号的绘图对象已经存在的话,将不建立新的对象,而只是让它成为当前绘图对象。
fig1 = pl.figure(1)
pl.subplot(211)
subplot(211)把绘图区域等分为2行*1列共两个区域, 而后在区域1(上区域)中建立一个轴对象. pl.subplot(212)在区域2(下区域)建立一个轴对象。
You can play around with plotting a variety of layouts. For example, Fig. 11 is created using the following commands:
f1 = pl.figure(1)
pl.subplot(221)
pl.subplot(222)
pl.subplot(212)
当绘图对象中有多个轴的时候,能够经过工具栏中的Configure Subplots按钮,交互式地调节轴之间的间距和轴与边框之间的距离。若是但愿在程序中调节的话,能够调用subplots_adjust函数,它有left, right, bottom, top, wspace, hspace等几个关键字参数,这些参数的值都是0到1之间的小数,它们是以绘图区域的宽高为1进行正规化以后的坐标或者长度。
pl.subplots_adjust(left=0.08, right=0.95, wspace=0.25, hspace=0.45)
4 绘制文件中的数据Plotting data contained in files
4.1 从Ascii文件中读取数据 Reading data from ascii files
读取文件的方法不少,这里只介绍一种简单的方法,更多的能够参考官方文档和NumPy快速处理数据(文件存取)。
numpy的loadtxt方法能够直接读取以下文本数据到numpy二维数组
**********************************************
# fakedata.txt
0 0
1 1
2 4
3 9
4 16
5 25
6 36
7 49
8 64
9 81
0 0
1 1
2 4
3 9
4 16
5 25
6 36
7 49
8 64
9 81
**********************************************
import numpy as npimport pylab as pl# Use numpy to load the data contained in the file# ’fakedata.txt’ into a 2-D array called datadata = np.loadtxt(’fakedata.txt’)# plot the first column as x, and second column as ypl.plot(data[:,0], data[:,1], ’ro’)pl.xlabel(’x’)pl.ylabel(’y’)pl.xlim(0.0, 10.)pl.show()
4.2 写入数据到文件 Writing data to a text file
写文件的方法也不少,这里只介绍一种可用的写入文本文件的方法,更多的能够参考官方文档。
import numpy as np# Let’s make 2 arrays (x, y) which we will write to a file# x is an array containing numbers 0 to 10, with intervals of 1x = np.arange(0.0, 10., 1.)# y is an array containing the values in x, squaredy = x*xprint ’x = ’, xprint ’y = ’, y# Now open a file to write the data to# ’w’ means open for ’writing’file = open(’testdata.txt’, ’w’)# loop over each line you want to write to filefor i in range(len(x)): # make a string for each line you want to write # ’\t’ means ’tab’ # ’\n’ means ’newline’ # ’str()’ means you are converting the quantity in brackets to a string type txt = str(x[i]) + ’\t’ + str(y[i]) + ’ \n’ # write the txt to the file file.write(txt)# Close your filefile.close()
这部分是翻译自:Python Plotting Beginners Guide
对LaTeX数学公式的支持
Matlplotlib对LaTeX有必定的支持,若是记得使用raw字符串语法会很天然:
xlabel(r"$\frac{x^2}{y^4}$")
在matplotlib里面,可使用LaTex的命令来编辑公式,只须要在字符串前面加一个“r”便可
Here is a simple example:
# plain text
plt.title('alpha > beta')
produces “alpha > beta”.
Whereas this:
# math text
plt.title(r'$\alpha > \beta$')
produces "".
这里给你们看一个简单的例子。
import matplotlib.pyplot as plt
x = arange(1,1000,1)
r = -2
c = 5
y = [5*(a**r) for a in x]
fig = plt.figure()
ax = fig.add_subplot(111)
ax.loglog(x,y,label = r"$y = \frac{1}{2\sigma_1^2}, c=5,\sigma_1=-2$")
ax.legend()
ax.set_xlabel(r"x")
ax.set_ylabel(r"y")
程序执行结果如图3所示,这其实是一个power-law的例子,有兴趣的朋友能够继续google之。
再看一个《用Python作科学计算》中的简单例子,下面的两行程序经过调用plot函数在当前的绘图对象中进行绘图:
plt.plot(x,y,label="$sin(x)$",color="red",linewidth=2)
plt.plot(x,z,"b--",label="$cos(x^2)$")
plot函数的调用方式很灵活,第一句将x,y数组传递给plot以后,用关键字参数指定各类属性:
label : 给所绘制的曲线一个名字,此名字在图示(legend)中显示。只要在字符串先后添加"$"符号,matplotlib就会使用其内嵌的latex引擎绘制的数学公式。
color : 指定曲线的颜色
linewidth : 指定曲线的宽度
详细的能够参考matplotlib官方教程:
Writing mathematical expressions
有几个问题:
matplotlib.rcParams属性字典
想要它正常工做,在matplotlibrc配置文件中须要设置text.markup = "tex"。
若是你但愿图表中全部的文字(包括坐标轴刻度标记)都是LaTeX'd,须要在matplotlibrc中设置text.usetex = True。若是你使用LaTeX撰写论文,那么这一点对于使图表和论文中其他部分保持一致是颇有用的。
在matplotlib中使用中文字符串时记住要用unicode格式,例如:u''测试中文显示''