[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

Re: [Numpy-discussion] Matlab/Numpy index order

2009-05-12 Thread brechmos
t changing the order? Thanks for the help. Ryan May-3 wrote: > > On Tue, May 12, 2009 at 2:51 PM, brechmos wrote: > >> So, in Numpy I have to reshape it so the "slices" are in the first >> dimension. Obviously, I can do a b.transpose( (1,2,0) ) to get it to >

[Numpy-discussion] Matlab/Numpy index order

2009-05-12 Thread brechmos
I am very new to Numpy and relatively new to Python. I have used Matlab for 15+ years now. But, I am starting to lean toward using Numpy for all my work. One thing that I am not understanding is the order of data when read in from a file. Let's say I have a 256x256x150 uint16 dataset (MRI, 150