.. _mplot3d-tricontourf3d_demo: mplot3d example code: tricontourf3d_demo.py =========================================== .. plot:: /home/tcaswell/src/p/matplotlib/doc/mpl_examples/mplot3d/tricontourf3d_demo.py :: """ ================================= Triangular 3D filled contour plot ================================= Filled contour plots of unstructured triangular grids. The data used is the same as in the second plot of trisurf3d_demo2. tricontour3d_demo shows the unfilled version of this example. """ import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D import matplotlib.tri as tri import numpy as np # First create the x, y, z coordinates of the points. n_angles = 48 n_radii = 8 min_radius = 0.25 # Create the mesh in polar coordinates and compute x, y, z. radii = np.linspace(min_radius, 0.95, n_radii) angles = np.linspace(0, 2*np.pi, n_angles, endpoint=False) angles = np.repeat(angles[..., np.newaxis], n_radii, axis=1) angles[:, 1::2] += np.pi/n_angles x = (radii*np.cos(angles)).flatten() y = (radii*np.sin(angles)).flatten() z = (np.cos(radii)*np.cos(angles*3.0)).flatten() # Create a custom triangulation. triang = tri.Triangulation(x, y) # Mask off unwanted triangles. xmid = x[triang.triangles].mean(axis=1) ymid = y[triang.triangles].mean(axis=1) mask = np.where(xmid*xmid + ymid*ymid < min_radius*min_radius, 1, 0) triang.set_mask(mask) fig = plt.figure() ax = fig.gca(projection='3d') ax.tricontourf(triang, z, cmap=plt.cm.CMRmap) # Customize the view angle so it's easier to understand the plot. ax.view_init(elev=45.) plt.show() Keywords: python, matplotlib, pylab, example, codex (see :ref:`how-to-search-examples`)