![](https://secure.gravatar.com/avatar/b951736e99a6f8d9129537b5fea565fa.jpg?s=120&d=mm&r=g)
My code is below. I am not projecting first. I tried just now, but I don't think I know exactly how I would do that, could you please elaborate? I can imagine x,y = m(x,y) ... but what about the lons, lats? They are not the same size. # Set up a basemap and get a figure ## This is just a convenience function to return an m instance and a fig handle fig,m = mp.get_base1(region='NPOLE') ## create interpolator print "creating interpolators" Z = mlab.griddata(x,y,z,lons,lats) #where x,y are in lon,lat for my data and #where lons,lats are from [-180,180], [50,90] regularly spaced at 0.5 degrees #transform to nx x ny regularly spaced native projection grid dx = 2.*np.pi*m.rmajor/len(lons) nx = int((m.xmax-m.xmin)/dx)+1; ny = int((m.ymax-m.ymin)/dx)+1 imdat = m.transform_scalar(Z,lons,lats,nx,ny) m.imshow(imdat) -- View this message in context: http://www.nabble.com/2d-interpolation%2C-non-regular-lat-lon-grid-tp2490968... Sent from the Scipy-User mailing list archive at Nabble.com.