Lambert azimuthal equal-area projection for full globe

[rwogburn]

The Lambert azimuthal equal-area projection can show the full globe, as in the example image in the Wikipedia article. However, with Basemap(projection='laea',...) I'm only able to plot a small portion of the globe. Once the width or height are > 15 million or so, it gives an error "ValueError: width and/or height too large for this projection, try smaller values".

Is there a good reason why basemap's laea projection doesn't do the full sphere? I'm using release 1.0.7_1 from Macports.

[micahcochran]

_choosecorners() or _choosecornersllur() are the only functions that raise that error, which are here in the code.

Could you write a simple example demonstrating the problem?

[rwogburn]

from mpl_toolkits.basemap import Basemap
# Size on a side in km -- 5000 works, 6000 gives error
skm = 6000
s = skm * 1000 * 3.141592
map = Basemap(projection='laea',lat_0=-90,lon_0=0,resolution='l',width=s,height=s)