I've been trying to help a friend who wants to plot directional data on a "Wulff net" <http://en.wikipedia.org/wiki/Pole_figure#Geometry_in_the_pole_figure>, which is a stereographic projection plot. She wants to plot points specified by latitude and longitude in degrees. We hoped to be able to use the basemap toolkit's "stere" plot, centred at lat_0=0, lon_0=0, with limits set at +/-90deg lat and +/-180deg lon, but we keep getting tracebacks and I wondered whether this is possible, based on a comment from Jeff Whitaker <http://www.nabble.com/-basemap--stereographic-projection-bounding-boxes-tf2170166.html#a6000978> which implies that basemap's stereographic projection code can't handle these default limits. Is this the case? Would it be asking too much to request a small sample of generating some Wulff-net axes? Plotting points on these seems simple enough. We started with the polarmaps.py example and the code below is as close as we could get. Any suggestions would be welcome.
Gary Ruben
···
--
from matplotlib.toolkits.basemap import Basemap
from pylab import *
# loop over projections, one for each panel of the figure.
fig = figure(figsize=(4,4))
# setup map projection
m = Basemap(projection='stere',lat_0=0.,lon_0=0.,llcrnrlat=-50.,llcrnrlon=-120., urcrnrlat=90., urcrnrlon=90.)
ax = fig.add_subplot(1,1,1)
# draw parallels and meridians.
m.drawparallels(arange(-180.,180.,10.))
m.drawmeridians(arange(-90.,90.,10.))
# draw boundary around map region.
m.drawmapboundary()
show()