Graph Plotting in Python | Set 1Graph Plotting in Python | Set 2
Matplotlib is a pretty extensive library which supports Animations of graphs as well. The animation tools center around the matplotlib.animation base class, which provides a framework around which the animation functionality is built. The main interfaces are TimedAnimation and FuncAnimation and out of the two, FuncAnimation is the most convenient one to use.
Installation:
Numpy: You can install numpy module using following pip command:
pip install numpy
FFMPEG: It is required only for saving the animation as a video. The executable can be downloaded from here.
Implementation:Python
# importing required modulesimportmatplotlib.pyplotaspltimportmatplotlib.animationasanimationimportnumpyasnp# create a figure, axis and plot elementfig=plt.figure()ax=plt.axes(xlim=(-50,50),ylim=(-50,50))line,=ax.plot([],[],lw=2)# initialization functiondefinit():# creating an empty plot/frameline.set_data([],[])returnline,# lists to store x and y axis pointsxdata,ydata=[],[]# animation functiondefanimate(i):# t is a parametert=0.1*i# x, y values to be plottedx=t*np.sin(t)y=t*np.cos(t)# appending new points to x, y axes points listxdata.append(x)ydata.append(y)# set/update the x and y axes dataline.set_data(xdata,ydata)# return line objectreturnline,# setting a title for the plotplt.title('A growing coil!')# hiding the axis detailsplt.axis('off')# call the animator anim=animation.FuncAnimation(fig,animate,init_func=init,frames=500,interval=20,blit=True)# save the animation as mp4 video fileanim.save('animated_coil.mp4',writer='ffmpeg',fps=30)# show the plotplt.show()
Here is how the output animation looks like:
Now, let us try to understand the code in pieces:
Here, we first create a figure, i.e a top level container for all our subplots.
Then we create an axes element ax which acts as a subplot. The range/limit for x and y axis are also defined while creating the axes element.
Finally, we create the plot element, named as line . Initially, the x and y axis points have been defined as empty lists and line-width (lw) has been set as 2.
def init():
line.set_data([], [])
return line,
Now, we declare a initialization function, init . This function is called by animator to create the first frame.
def animate(i):
# t is a parameter
t = 0.1*i
# x, y values to be plotted
x = t*np.sin(t)
y = t*np.cos(t)
# appending new points to x, y axes points list
xdata.append(x)
ydata.append(y)
# set/update the x and y axes data
line.set_data(xdata, ydata)
# return line object
return line,
This is the most important function of above program. animate() function is called again and again by the animator to create each frame. The number of times this function will be called is determined by number of frames, which is passed as frames argument to animator.
animate() function takes the index of ith frame as argument.
t = 0.1*i
Here, we cleverly use the index of current frame as a parameter!
x = t*np.sin(t)
y = t*np.cos(t)
Now, since we have the parameter t, we can easily plot any parametric equation. For example, here, we are plotting a spiral using its parametric equation.
line.set_data(xdata, ydata)
return line,
Finally, we use set_data() function to set x and y data and then return plot object, line .
Now, we create the FuncAnimation object, anim . It takes various arguments explained below:
fig : figure to be plotted.
animate : the function to be called repeatedly for each frame.init_func : function used to draw a clear frame. It is called once before the first frame.
frames : number of frames. (Note: frames can also be an iterable or generator.)
interval : duration between frames ( in milliseconds)
blit : setting blit=True means that only those parts will be drawn, which have changed.
Now, we save the animator object as a video file using save() function. You will need a movie writer for saving the animation video. In this example, we have used FFMPEG movie writer. So, writer is set as 'ffmpeg'.
fps stands for frame per second.
Example 2
This example shows how one can make a rotating curve by applying some simple mathematics!
Python
# importing required modulesimportmatplotlib.pyplotaspltimportmatplotlib.animationasanimationimportnumpyasnp# create a figure, axis and plot elementfig=plt.figure()ax=plt.axes(xlim=(-25,25),ylim=(-25,25))line,=ax.plot([],[],lw=2)# initialization functiondefinit():# creating an empty plot/frameline.set_data([],[])returnline,# set of points for a star (could be any curve)p=np.arange(0,4*np.pi,0.1)x=12*np.cos(p)+8*np.cos(1.5*p)y=12*np.sin(p)-8*np.sin(1.5*p)# animation functiondefanimate(i):# t is a parametert=0.1*i# x, y values to be plottedX=x*np.cos(t)-y*np.sin(t)Y=y*np.cos(t)+x*np.sin(t)# set/update the x and y axes dataline.set_data(X,Y)# return line objectreturnline,# setting a title for the plotplt.title('A rotating star!')# hiding the axis detailsplt.axis('off')# call the animator anim=animation.FuncAnimation(fig,animate,init_func=init,frames=100,interval=100,blit=True)# save the animation as mp4 video fileanim.save('basic_animation.mp4',writer='ffmpeg',fps=10)# show the plotplt.show()
Here is how the output of above program looks like:
Here, we have used some simple mathematics to rotate a given curve.
The star shape is obtained by putting k = 2.5 and 0<t<4*pi in the parametric equation given below:
x= [a-b]\cos (t) + b\cos [t(\frac{a}{b}-1)]
y= [a-b]\sin (t) - b\sin [t(\frac{a}{b}-1)] , k= \frac{a}{b}
The same has been applied here:
p = np.arange(0, 4*np.pi, 0.1)
x = 12*np.cos(p) + 8*np.cos(1.5*p)
y = 12*np.sin(p) - 8*np.sin(1.5*p)
Now, in each frame, we rotate the star curve using concept of rotation in complex numbers. Let x, y be two ordinates. Then after rotation by angle theta, the new ordinates are:
{x}' = x\cos \theta - y\sin \theta
{y}' = x\sin \theta + y\cos \theta
The same has been applied here:
X = x*np.cos(t) - y*np.sin(t)
Y = y*np.cos(t) + x*np.sin(t)
All in all, animations are a great tool to create amazing stuff and many more things can be created using them.
So, this was how animated plots can be generated and saved using Matplotlib.