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() ------------------------------------------------------------------------------ Centralized Desktop Delivery: Dell and VMware Reference Architecture Simplifying enterprise desktop deployment and management using Dell EqualLogic storage and VMware View: A highly scalable, end-to-end client virtualization framework. Read more! http://p.sf.net/sfu/dell-eql-dev2dev _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
