.. _examples-meteoinfolab-plot_types-plot_3d:

*************************
3-D plots without OpenGL
*************************

Create a 3D axes without OpenGL using ``ax = axes3d(opengl=False)``, then use the plot functions of the 3D axes object (ax)
to plot 3D plots. The functions includes: including ``plot()``, ``scatter()``,
``contour()``, ``contourf()``, ``imshow()``, ``plot_surface()``, ``plot_wireframe()`` and
``plot_layer``.

**Line plots**::

    z = linspace(0, 1, 100)
    x = z * np.sin(20 * z)
    y = z * np.cos(20 * z)

    #Plot
    ax = axes3d(opengl=False)
    ax.plot(x, y, z, '-b')
    title('3D plot example')

.. image:: ../../../_static/plot_3d.png

**Scatter plots**::

    z = linspace(0, 1, 100)
    x = z * np.sin(20 * z)
    y = z * np.cos(20 * z)
    c = x + y

    #Plot
    ax = axes3d(opengl=False)
    points = ax.scatter(x, y, z, c=c)
    colorbar(points,shrink=0.8)
    title('Point 3D plot example')
    
.. image:: ../../../_static/scatter_3d.png

**Wireframe and contour plots**::

    alpha = 0.7
    phi_ext = 2 * pi * 0.5
    N = 25
    x1 = linspace(0, 2*pi, N)
    y1 = linspace(0, 2*pi, N)
    x,y = meshgrid(x1, y1)
    z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
    z = z.T

    #Plot
    ax = axes3d(opengl=False)
    lines = ax.contour(x1, y1, z, 10, offset=-2)
    ax.plot_wireframe(x, y, z, color='b')
    colorbar(lines)
    title('Wireframe and contour 3D plot example')
    
.. image:: ../../../_static/wireframe_contour.png

**Wireframe and contourf plots**::

    alpha = 0.7
    phi_ext = 2 * pi * 0.5
    N = 25
    x1 = linspace(0, 2*pi, N)
    y1 = linspace(0, 2*pi, N)
    x,y = meshgrid(x1, y1)
    z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
    z = z.T

    #Plot
    ax = axes3d(opengl=False)
    lines = ax.contourf(x1, y1, z, 10, offset=-2)
    ax.plot_wireframe(x, y, z, color='b')
    colorbar(lines)
    title('Wireframe and contourf 3D plot example')
    
.. image:: ../../../_static/wireframe_contourf.png

**Surface plots**::

    alpha = 0.7
    phi_ext = 2 * pi * 0.5
    x = linspace(0, 2*pi, 100)
    y = linspace(0, 2*pi, 100)
    x,y = meshgrid(x, y)
    z = 2 + alpha - 2 * cos(y) * cos(x) - alpha * cos(phi_ext - 2 * y)
    z = z.T

    #Plot
    ax = axes3d(opengl=False)
    ls = ax.plot_surface(x, y, z, 20, edge=False)
    colorbar(ls,shrink=0.8)
    title('Surface 3D plot example')
    
.. image:: ../../../_static/surface_1.png

**Image plots**::

    fn = 'D:/Temp/nc/air_clm.nc'
    f = addfile(fn)
    ps = f['aveair'][0,:,:,'120']
    yy = linspace(0, 1., ps.shape[0])
    ps.setdimvalue(0, yy)

    #Map layer
    layer = shaperead('D:/Temp/map/110m-land.shp')

    #Plot
    ax = axes3d(opengl=False)
    ax.plot_layer(layer, color='c', edgecolor='b')
    ls = ax.imshow(ps, 10, offset=120, zdir='x', alpha=0.8)
    colorbar(ls)
    zlim(0, 1)
    xlim(0, 180)
    title('3D imshow x direction example')
    
.. image:: ../../../_static/imshow_xdir_3d.png

**Trajectory plots**::

    #Open trajectory data and get trajectory layer
    fn = 'D:/Temp/HYSPLIT/traj_20131211_00'
    f = addfile_hytraj(fn)
    lon = f['lon'][:]
    lat = f['lat'][:]
    alt = f['height'][:]
    pres = f['PRESSURE'][:]

    #Relief data
    fn = 'D:/Temp/nc/elev.0.25-deg.nc'
    f = addfile(fn)
    elev = f['data'][0,::8,::8]
    elev = elev[:,'0:180']
    elev[elev<0] = 0

    #Plot
    ax = axes3d(opengl=False)
    surf(elev, 20, cmap='MPL_terrain', edge=False)
    geoshow('coastline', facecolor='g')
    traj = plot3(lon, lat, alt, mvalues=pres, linewidth=2)
    scatter3(lon[:,0], lat[:,0], alt[:,0], marker='o', fill=False,
        size=14, edgecolor='gray')
    zlim(0, 10000)
    xlim(0, 180)
    colorbar(traj, shrink=0.8, aspect=30)
    title('3D trajectory example')
    
.. image:: ../../../_static/trajectory_3d_1.png