Re: [Numpy-discussion] 3D interpolation of large array

2009-09-14 Thread Gael Varoquaux
On Sun, Sep 13, 2009 at 11:55:43PM -0400, David Warde-Farley wrote: > On 13-Sep-09, at 8:51 PM, brechmos wrote: > > I have a large dataset (e.g., 70 x 500 x 500) and want to > > interpolate points > > (for example to double the size). What it seems I want is: > > []: newx,newy,newz=mgrid[1:70

Re: [Numpy-discussion] 3D interpolation of large array

2009-09-13 Thread David Warde-Farley
On 13-Sep-09, at 8:51 PM, brechmos wrote: > > I have a large dataset (e.g., 70 x 500 x 500) and want to > interpolate points > (for example to double the size). What it seems I want is: > > []: newx,newy,newz=mgrid[1:70:0.5,1:500:0.5,1:500:0.5] > []: coords = array([newz, newy, newx]) > []: dou

[Numpy-discussion] 3D interpolation of large array

2009-09-13 Thread brechmos
I have a large dataset (e.g., 70 x 500 x 500) and want to interpolate points (for example to double the size). What it seems I want is: []: newx,newy,newz=mgrid[1:70:0.5,1:500:0.5,1:500:0.5] []: coords = array([newz, newy, newx]) []: dout = np.map_coordinates(d, coords) The problem is that mgri