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
|
