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
Benjamin Root-2 wrote:
>
> On Mon, Sep 13, 2010 at 8:21 AM, mdekauwe <[email protected]> wrote:
>
>>
>> Hi,
>>
>> Well hopefully doing what you suggested correctly...
>>
>> 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))
>> x = np.linspace(-1.74, -1.32, 16)
>> y = np.linspace(15.15, 15.61, 17)
>> px, py = m(x, y)
>> X, Y = np.meshgrid(px, py)
>> m.contourf(X, Y, data2, colors='black')
>> plt.show()
>>
>> Doesn't seem to work either?
>>
>> thanks,
>>
>> Martin
>>
>>
> Martin,
>
> Double-check the shape of the X and the Y arrays. Make sure they have the
> same shape as data2. I would be willing to bet that some of the shapes
> got
> reversed.
>
> Ben Root
>
> ------------------------------------------------------------------------------
> 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
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
--
View this message in context:
http://old.nabble.com/basemap-and-contouring-tp29697864p29698610.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
[email protected]
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
