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.

customticks.py shows how to create custom tick labels for a cylindrical
projection.

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)

hires.py illustrates the use of the high-resolution coastlines, by
drawing the coastlines, political boundaries and rivers of the British
Isles. 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 that accesses
data over http using the dap module. An internet connection is required.

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 and/or an internet connection).