[Matplotlib-users] basemap map center issue

[Andres L. <and...@gm...>]

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

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()