Hi folks,<br><br>This is probably a very simple question, but it has me stumped...<br><br>I have an integer 2D array containing 3rd dimesion indicies that I'd like to use to index values in a 3D array.<br><br>Basically, I want the equivalent of:<br>

<br>> output = np.zeros((ny,nx))<br>><br>> for i in xrange(ny):<br>>     for j in xrange(nx):<br>>         z = grid[i,j]<br>>         output[j,i] = bigVolume[j,i,z]<br><br>Where grid is my 2D array of indicies and bigVolume is my 3D array.<br>

<br>I've read the numpy-user and numpybook sections on indexing with an
integer ndarray, but I'm still not quite able to wrap my head around
how it should work.  I'm sure I'm missing something obvious (I haven't been using numpy particularly long).  <br><br>If it helps anyone visualize it, I'm essentially trying to extract attributes out of a seismic volume along the surface of a horizion.<br>

<br>Thanks!<br>-Joe<br><br><br><br><br>