[Numpy-discussion] 3D interpolation of large array

[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 mgrid can't do such a large set of values (at least on
my setup! MacBook Pro with 4Gig RAM).  The error from the mgrid command is
"ValueError: dimensions too large."

My version of Numpy is numpy-1.3.0n1-py2.5-macosx-10.3 (from the Enthought
package).

Is there any other way that I could do this?
-- 
View this message in context: http://www.nabble.com/3D-interpolation-of-large-array-tp25428880p25428880.html
Sent from the Numpy-discussion mailing list archive at Nabble.com.