.. _pylab_examples-tricontour_vs_griddata:

pylab_examples example code: tricontour_vs_griddata.py
======================================================



.. plot:: /media/TOSHI/temp/numpy_scipy_matplotlib/matplotlib/matplotlib-1.1.1~rc2/doc/mpl_examples/pylab_examples/tricontour_vs_griddata.py

::

    """
    Comparison of griddata and tricontour for an unstructured triangular grid.
    """
    import matplotlib.pyplot as plt
    import matplotlib.tri as tri
    import numpy as np
    from numpy.random import uniform, seed
    from matplotlib.mlab import griddata
    import time
    
    seed(0)
    npts = 200
    ngridx = 100
    ngridy = 200
    x = uniform(-2,2,npts)
    y = uniform(-2,2,npts)
    z = x*np.exp(-x**2-y**2)
    
    # griddata and contour.
    start = time.clock()
    plt.subplot(211)
    xi = np.linspace(-2.1,2.1,ngridx)
    yi = np.linspace(-2.1,2.1,ngridy)
    zi = griddata(x,y,z,xi,yi,interp='linear')
    plt.contour(xi,yi,zi,15,linewidths=0.5,colors='k')
    plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
    plt.colorbar() # draw colorbar
    plt.plot(x, y, 'ko', ms=3)
    plt.xlim(-2,2)
    plt.ylim(-2,2)
    plt.title('griddata and contour (%d points, %d grid points)' % (npts, ngridx*ngridy))
    print 'griddata and contour seconds:', time.clock() - start
    
    # tricontour.
    start = time.clock()
    plt.subplot(212)
    triang = tri.Triangulation(x, y)
    plt.tricontour(x, y, z, 15, linewidths=0.5, colors='k')
    plt.tricontourf(x, y, z, 15, cmap=plt.cm.jet)
    plt.colorbar()
    plt.plot(x, y, 'ko', ms=3)
    plt.xlim(-2,2)
    plt.ylim(-2,2)
    plt.title('tricontour (%d points)' % npts)
    print 'tricontour seconds:', time.clock() - start
    
    plt.show()
    

Keywords: python, matplotlib, pylab, example, codex (see :ref:`how-to-search-examples`)