Tim Gray wrote:

Hi, I've got a question for everyone. I'm trying to some code from IDL that uses the interpolate function. Basically I've got a big 2d array of data with time as one index and radius on the other index. Each value in the array is "data". I basically need to interpolate this data to other radii.

Here's what IDL's interpolate does (to 2d arrays). It takes an array such as: x =[[ 0, 1, 2, 3,] [ 4, 5, 6, 7,] [ 8, 9,10,11,] [12,13,14,15,]]

and given two vectors which define the desired indices to interpolate the data to, it spits out a new array. So for an input of: interpolate(x,[.5,1.5,2.5],[.5,1.5,2.5]), you get a 3x3 array of: [[ 2.5, 3, 4.5] [ 6.5, 7.5, 8.5] [10.5, 11.5, 12.5]]

I'm sure something like this is in Scipy already, I'm just having problems finding it. Any suggestions?

Look at scipy.interpolate.interpolate.interp2d . Yes, there are 2 "interpolate"s in there. I don't know why interp2d isn't exposed like interp1d is. Pearu, is this an oversight? Or is interp2d not robust enough, yet? -- Robert Kern rkern@ucsd.edu "In the fields of hell where the grass grows high Are the graves of dreams allowed to die." -- Richard Harter