Read Me

Basemap example scripts

simpletest.py is the simple example given in the Basemap docstrings.

simpletest_oo.py is the same example, but using the "OO" interface (does
not import pylab).

test.py exercises all the available map projections by plotting the etopo 
20 minute bathymetry/topography data using imshow.

contour_demo.py demonstrates the use of filled contours with map projections.

plotmap.py is the example on the matplotlib 'screenshots' page (included
in test.py) which shows the ETOPO topography as an image on a Lambert
Conformal projection (using imshow).

plotmap_oo.py is a version of plotmap.py that does not import pylab - it uses
the matplotlib OO interface instead.

plotmap_masked.py is a version of plotmap.py that shows how to mask out
portions of the plot (in this case the oceans) using pcolor and imshow.

nytolondon.py shows how to draw a great circle on a map (NY to London)

ireland.py draws maps of ireland with the crude, low and intermediate
resolution boundary datasets. Also does high res coastlines if they
are installed. Shows how to pickle and re-load a
Basemap class instance (useful for high-res coastlines when you are plotting
the same small map region many times in different scripts).

quiver_demo.py shows how to plot wind vectors on a map.

randompoints.py demonstrates the use of scatter to plot randomly distributed
points on the earth.

ortho_demo.py demonstrates the orthographic projection - you can interactively
enter the reference lat/lon (the viewpoint of the satellite). Also shows
how to plot the default 5 minute land-sea mask, specifying the colors
to be used for land and ocean areas.

geos_demo.py demonstrates the geostationary projection - you can interactively
enter the reference longitude and height above the equator of the
satellite.

geos_demo_2.py demonstrates how to read in a geostationary satellite image
from a jpeg file, then plot only a portion of the full earth (contributed
by Scott Sinclair, requires PIL).
fcstmaps.py is a sample multi-panel plot.  Care is taken to preserve the aspect ratio of the map in each panel, and a common title and colorbar is created.  
Requires the opendap module (http://opendap.oceanografia.org) and an 
active internet connection to fetch the data to be plotted.

wiki_example.py is the example from the MatplotlibCookbook scipy wiki page
(http://www.scipy.org/wikis/topical_software/MatplotlibCookbook/wikipage_view).

fillstates.py is an example that shows how to draw 
and fill polygons (in this case U.S. States) from
and ESRI shapefile. 

hurrtrack.py plots hurricane tracks from shapefile data 
(from nationalatlas.gov).

testgdal.py is an example that shows how to plot raster geospatial data
using the gdal module (http://gdal.maptools.org).  Requires two data
files that must be downloaded manually, see docstrings for URLs.

panelplot.py shows how to make multi-panel plots, taking special care to 
make sure that map aspect ratios are preserved (so you don't get squished
and/or stretched maps).

polarmaps.py illustrates the use of the special-case versions of the 
stereographic, lambert azimuthal equal area, azimuthal equidistant 
projections (prefixed by 'np' and 'sp').

warpimage.py shows how to interpolate (warp) an image from one
map projection to another. Requires PIL, and an image from
http://www.space-graphics.com/earth_topo-bathy.htm, 

garp.py makes a 'World According to Garp' map - an azimuthal equidistant
projection centered on a specified location.  Straight lines from that
location to any point on the map have true scale and show the shortest
route.  Shows how to use the 'width' and 'height' keywords
to specify the map projection region (instead of specifying
the lat/lon of the upper right and lower left corners).

More examples of using the 'width' and 'height' keywords (instead of
'llcrnrlat','llcrnrlon','urcrnrlat','urcrnrlon') are given in setwh.py.

plot_tissot.py - tissot's indicatrix, shows the distorting effects of
different map projections. Exercises the shapefile methods.

pnganim.py shows how to make an animation.

ccsm_popgrid.py shows how to plot data on 2d grids which are not 
rectilinear in map projection or geographic coordinates.  The NCAR
Parallel Ocean Program grid is used as an example.

show_colormaps.py plots the extra colormaps included with basemap.

plotprecip.py use nonlinear precip colormap included with basemap
to make a rainfall plot.

run_all.py is a driver script that runs all the examples except fcstmaps.py,
testgdal.py, geos_demo_2.py, warpimage.py, and pnganim.py (which
rely on external dependencies).