Module gmt.gmt

Introduction

• Contour, boxfill or windbarb maps on various map projections including linear, cylindrical, lambert conformal, orthographic, lambert azimuthal equal area, mercator, north polar stereographic, south polar stereographic, and general stereographic.
• Drawing of symbols, text, lines, curves and polygons on maps.
• Overlaying plots.
• Easy creation of multi-panel figures.
• simple x-y line plots.

Download

• Download module source code. To install, untar and run python setup.py install.
• Sample code and plots.

Requirements

• Python 2.3 or higher.
• GMT 3.4.4 with the windbarb patch.
• ghostscript.
• Numeric.
• Scientific.
• PyX.
• pydods (only needed to run examples).
• To make nicely formatted pdf and html versions of the docstrings, you'll need epydoc.

Example usage

• create a class instance:

    >>> x = GMT()  # see below for optional keyword arguments).
  

• plot a basemap:

    >>> x.basemap() # see below for class attributes that can be set
    >>>             # to affect how the basemap is drawn.
  

• create some GMT 'grdfiles' to plot:

    >>> datgrd = x.wrtgrd(data,lons,lats)
    >>> ugrd = x.wrtgrd(u,lons,lats) # zonal winds.
    >>> vgrd = x.wrtgrd(v,lons,lats) # meridional winds.
  

• plot some contours:

    >>> x.contour(datgrd)  # see below for relevant class attributes.
  

• draw some windbarbs:

    >>> x.windbarbs(ugrd,vgrd) # see below for relevant class attributes.
  

• draw some symbols and text:

    >>> x.drawsymbol(40,254.7,symbol='c',color=(255,0,0))
    >>> xpos,ypos = x.latlontoxy(40,254.7) # convert from lat/lon to x/y.
    >>> x.drawtext('Boulder',xpos,ypos,offset=0.1,justify='CB')
  

• save the grd file for future usage (otherwise it will be lost when close method called):

    >>> savegrd(datgrd,'dat.grd')
  

• close the class instance, cleanup temp files, save postscript:

    >>> x.close(psfilename='test.eps')   
  

• or, instead of saving postscript to file, save in a StringIO instance (useful if the figure is going to be part of a multi-panel plot - see info on PanelPlot class below):

    >>> figpanel = x.close() # figpanel is a StringIO instance 
  

Version: 20041114

Author: Jeffrey Whitaker with contributions from Julien Vienne.

Contact: Jeff Whitaker

Copyright: copyright 2004 by Jeffrey Whitaker.

License: Permission to use, copy, modify, and distribute this software and its documentation for any purpose and without fee is hereby granted, provided that the above copyright notice appear in all copies and that both that copyright notice and this permission notice appear in supporting documentation. THE AUTHOR DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.

Bug: contourfill doesn't work when proj = 'ortho' (use boxfill=True instead).