Hi Jeff, many thanks that was exactly the problem!
Martin Jeff Whitaker wrote: > > On 9/13/10 8:02 AM, mdekauwe wrote: >> Hi, >> >> Well I tried... >> >> numrows = 17 >> numcols = 16 >> ulat, llat, ulon, llon = 15.61, 15.15, -1.32, -1.74 >> m = Basemap(projection='geos', lon_0=0.0, llcrnrlon=llon, >> llcrnrlat=llat, urcrnrlon=ulon, urcrnrlat=ulat, >> resolution='c') >> >> data2 = np.random.sample((17,16)) >> print 'DATA2 shape:', data2.shape >> x = np.linspace(-1.74, -1.32, 16) >> y = np.linspace(15.15, 15.61, 17) >> print 'x shape:', x.shape >> print 'y shape:', y.shape >> px, py = m(x, y) >> >> >> which gives... >> >> DATA2 shape: (17, 16) >> x shape: (16,) >> y shape: (17,) >> >> So that looks OK. Turns out the error now comes at the line px, py = m(x, >> y). So perhaps I am doing the basemap step incorrectly? >> >> Traceback (most recent call last): >> File "./recontour_plot.py", line 83, in<module> >> px, py = m(x, y) >> File >> "/users/eow/mgdk/sun4u//lib/python/mpl_toolkits/basemap/__init__.py", >> line 823, in __call__ >> return self.projtran(x,y,inverse=inverse) >> File "/users/eow/mgdk/sun4u//lib/python/mpl_toolkits/basemap/proj.py", >> line 241, in __call__ >> outx,outy = self._proj4(x, y, inverse=inverse) >> File >> "/users/eow/mgdk/sun4u//lib/python/mpl_toolkits/basemap/pyproj.py", >> line 193, in __call__ >> _Proj._fwd(self, inx, iny, radians=radians, errcheck=errcheck) >> File "_proj.pyx", line 56, in _proj.Proj._fwd (src/_proj.c:725) >> RuntimeError: Buffer lengths not the same >> >> > Martin: > > You need a meshgrid > > x,y = np.meshgrid(x,y) > > just *before* (not after) > > px,py = m(x,y) > > so that x and y have the same shape as data2. > > -Jeff > > -- > Jeffrey S. Whitaker Phone : (303)497-6313 > Meteorologist FAX : (303)497-6449 > NOAA/OAR/PSD R/PSD1 Email : jeffrey.s.whita...@noaa.gov > 325 Broadway Office : Skaggs Research Cntr 1D-113 > Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg > > > ------------------------------------------------------------------------------ > Start uncovering the many advantages of virtual appliances > and start using them to simplify application deployment and > accelerate your shift to cloud computing > http://p.sf.net/sfu/novell-sfdev2dev > _______________________________________________ > Matplotlib-users mailing list > Matplotlib-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/matplotlib-users > > -- View this message in context: http://old.nabble.com/basemap-and-contouring-tp29697864p29699166.html Sent from the matplotlib - users mailing list archive at Nabble.com. ------------------------------------------------------------------------------ Start uncovering the many advantages of virtual appliances and start using them to simplify application deployment and accelerate your shift to cloud computing http://p.sf.net/sfu/novell-sfdev2dev _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
