Thread: [Matplotlib-users] Round plot with 'npstere'-projection

Dear all,

using a polarstereograhic projection with Basemap (

m = Basemap(projection = 'npstere', resolution = 'c', lat_0 = 90., lon_0 = 0., 
area_thresh = 1000., boundinglat = 40.

)

one get a square map. Is there any way to make the map round, not square?

Many thanks,

On 1/16/12 11:08 AM, Serhiy Hrechanyy wrote:
>>> m = Basemap(projection = 'npstere', resolution = 'c', lat_0 = 90.,
>>> lon_0 = 0.,
>>> area_thresh = 1000., boundinglat = 40.
>>>
>>> )
>>>
>>> one get a square map. Is there any way to make the map round, not 
>>> square?
>> No.
>
> Thanks, Jeff,
>
> could you suggest any other projection to do the job (round polar plot).
>
> 'ortho' makes a round map, but goes almost till Equator (I would 
> rather prefer 40°)
>
> Do 'npaeqd' and 'nplaea' make only "squares" as well?
Serhiy:

Yes, they all make squares.  Why is it so important the plot be round?

-Jeff

> Yes, they all make squares. Why is it so important the plot be round?

Well, it is definitely not a matter of life or death. Just think that 
round polar plots look a bit nicer.

I haven't tried your workaround solution yet (must do something urgent 
first). But thanks a lot for your help and of course for the  Basemap. 
Very helpful tool even with square maps))!!

Serhiy

On 1/17/12 3:50 AM, Serhiy Hrechanyy wrote:
>> Yes, they all make squares. Why is it so important the plot be round?
>
> Well, it is definitely not a matter of life or death. Just think that 
> round polar plots look a bit nicer.
>
> I haven't tried your workaround solution yet (must do something urgent 
> first). But thanks a lot for your help and of course for the  Basemap. 
> Very helpful tool even with square maps))!!
>
> Serhiy
>
Serhiy:  I've added this feature to my github fork 
(https://github.com/jswhit/basemap).  You just specify 'round=True' when 
creating the pole-centered basemap.   Meridian labels are not working 
yet though.

-Jeff

On 1/16/12 11:08 AM, Serhiy Hrechanyy wrote:
>>> m = Basemap(projection = 'npstere', resolution = 'c', lat_0 = 90.,
>>> lon_0 = 0.,
>>> area_thresh = 1000., boundinglat = 40.
>>>
>>> )
>>>
>>> one get a square map. Is there any way to make the map round, not 
>>> square?
>> No.
>
> Thanks, Jeff,
>
> could you suggest any other projection to do the job (round polar plot).
>
> 'ortho' makes a round map, but goes almost till Equator (I would 
> rather prefer 40°)
>
> Do 'npaeqd' and 'nplaea' make only "squares" as well?
Serhiy:  Here's a hacky workaround that does what you want.  Run it in 
the matplotlib examples directory (it uses data files in there).

-Jeff