[Matplotlib-users] basemap, boundinglat, and spstere

[Ken M. <ma...@gm...>]

Hi,

I'm having trouble with Basemap South Pole Stereographic
projections. The required "boundinglat" doesn't seem to work properly,
and I hope someone can help me figure out the correct syntax.

I have the following information about the data:

extent: -3,333,500 to 3,333,500 m
true scale latitude: -71
south pole stereographic projection.

Given this, I try the following:

m = Basemap(resolution='c',projection='spstere',
lat_ts=-71,lat_0=0.,lon_0=0.)

And it crashes reporting that I need boundinglat. I can't give
width/height which are accepted by the 'stere' projection.

I use http://www.pgc.umn.edu/tools/conversion and determine that
3333500 is 60 degrees south. If then add "boundinglat=-60" to the
Basemap call, the data appears nowhere near the correct location.

If I add "boundinglat=-89.999999" it looks almost perfect. If I change
that to "-90" or add a few more 9's, it crashes.

Am I determining the bounding lat correctly from the provided
information? Should I be calling Basemap differently? Any help will be
much appreciated.

Thanks,

  Ken Mankoff