[SciPy-User] [SciPy-user] 2d interpolation, non-regular lat/lon grid - help with delauney/natgrid??

John [H2O] washakie at gmail.com
Sat Sep 26 14:55:57 EDT 2009




Robert Kern-2 wrote:
> 
> 
> 
> Just transform reg_lon,reg_lat to the projection space an interpolate
> using that.
> 
> Z0 = interp(*m(reg_lon, reg_lat))
> 
> 

I get the following error now:

Traceback (most recent call last):
  File "./irregular_interp.py", line 541, in <module>
    main()
  File "./irregular_to_regulargrid.py", line 78, in main
    triangulate=method_options[1])
  File "./irregular_to_regulargrid.py", line 474, in grid_points
    newx,newy,Z =
regrid(x,y,z,m,dres=dres,method=method,triangulate=triangulate)
  File "./irregular_to_regulargrid.py", line 336, in regrid
    Z0 = interp(*m(newx,newy)) 
  File "/dist/site-packages/mpl_toolkits/basemap/__init__.py", line 823, in
__call__
    return self.projtran(x,y,inverse=inverse)
  File "/dist/site-packages/mpl_toolkits/basemap/proj.py", line 241, in
__call__
    outx,outy = self._proj4(x, y, inverse=inverse)
  File "/dist/site-packages/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:876)
RuntimeError: Buffer lengths not the same


-- 
View this message in context: http://www.nabble.com/2d-interpolation%2C-non-regular-lat-lon-grid-tp24909685p25627893.html
Sent from the Scipy-User mailing list archive at Nabble.com.




More information about the SciPy-User mailing list