.. _news-meteoinfo_1.4.6:


******************************************
MeteoInfo 1.4.6 was released (2017-7-31)
******************************************

Several 3D plot function were added in Axes3D class, 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()
    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()
    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()
    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()
    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()
    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(bbox=True)
    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)
    tlayer = f.trajlayer()
    stlayer = f.trajsplayer()

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

    #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()
    ls = ax.plot_surface(elev, 20, cmap='MPL_terrain', edge=False)
    ax.plot_layer(layer, edgecolor='g')
    ax.plot_layer(tlayer)
    ax.plot_layer(stlayer, fill=False)
    zlim(0, 10000)
    xlim(0, 180)
    colorbar(ls)
    title('3D trajectory example')
    
.. image:: ../_static/trajectory_3d_1.png