|
From: Andres L. <and...@gm...> - 2010-11-11 10:33:38
|
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() |