I've been playing around with stereographic projections, and it appears that the bounding-box for the 'stere' projection isn't computed. Being blissfully unaware of all the complexities involved, I thought I'd send this email to see if there's an easy way to deal with the situation.
Basically, I'm trying to plot on an equatorial stereographic projection created like so: m=Basemap(projection='stere', lat_ts = 0.0, lat_0 = 0, lon_0 = 90.0) This creates a projection that seems to work well, but the problem is that the [ll|ur]crn* attributes aren't set to anything useful, meaning that lots of the nice basemap goodies aren't working (drawmeridians, etc.) (Note: I don't really understand what the lat_ts is doing -- the docstring says it is the "natural origin". I can't find reference to this at http://en.wikipedia.org/wiki/Stereographic_projection . I guess this is related to the discussion of "scale error at the extremities" at http://www.remotesensing.org/geotiff/proj_list/oblique_stereographic.html , and thus I assume that if lat_ts is set to lat_0, the scale factor will be 1 at the point of tangency. ) So, the question is, can some default bounding box be computed for equatorial (and possibly oblique) stereographic projections in any reasonable way? It would be nice to make plots like those found on the wikipedia page. Cheers! Andrew ------------------------------------------------------------------------- Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
