Re: [Matplotlib-users] basemap map center issue

[Jeff W. <js...@fa...>]

On 11/11/10 3:33 AM, Andres Luhamaa wrote:
> Hi,
> I have a problem creating map in "aeqd" projection. The problem seems 
> to be that the map center, given by arguments lon_0 and lat_0 is 
> taking into account only full part of the coordinate. If I want to get 
> map center at coordinates lon_0=25.5 and lat_0=58.5, I get center in 
> lon_0=25 and lat_0=58, instead. I run numpy and matplotlib dev 
> versions (a few weeks old) and basemap 1.0
>
> Best regards,
> Andres Luhamaa

Andres:  I can't reproduce this with GTKAgg or any other backend, with 
basemap svn.  I don't see anything that has changed since 1.0 that would 
have affected this behavior though.

-Jeff
>
> I use following example:
>
> import pylab as plt
> from mpl_toolkits.basemap import Basemap
> lat_0=58.5
> lon_0=25.5
> width=400000
> height=400000
>
> m=Basemap(projection='aeqd',lat_0=lat_0, 
> lon_0=lon_0,width=width,height=height,resolution='i')
>
> m.drawcoastlines()
> m.drawmeridians(range(20,30),labels=[0,0,1,0])
> m.drawparallels(range(45,65),labels=[1,0,0,0])
>
> ## what I expect to be map center
> X,Y=m(lon_0,lat_0)
> m.plot([X],[Y],'ko')
> plt.text(X,Y,'expected')
>
> ## actual map center
> X,Y=m(int(lon_0),int(lat_0))
> m.plot([X],[Y],'ro')
> plt.text(X,Y,'actual')
>
> plt.show()