How to I find the values associated with a vertex position on an isosurface?

Hello,

I am interested in assessing the vector field values at the vertices of an isosurface as using 'find' does not seem to work the way I would expect.

I used 'isosurface(xx,yy,zz,omega,ISOVALUE)' to generate the coherent vortex structure from measured velocity data that was extracted from 4D MRI. In order to make the field amenable to vector operations, I interpolated the array data onto a regular, 3D grid. When I perform this step, I am able to successfully visualise the structures I expect to appear in the flow.

However, I have not been able to find a good way to find the point in the 3D grid of vorticity magnitude values 'omega' that correspond with the vertex positions on the isosurface. I tried to use 'find()' in a for loop like so:

S = isosurface(xx,yy,zz,omega,ISOVALUE);

v = S.vertices;

for i = 1:length(v)

idx(i) = find(xx == v(i,1) & yy == v(i,2) & zz == v(i,3);

end

but when I try this, the method is unable to identify any values for yy that equal any values in the corresponding column of the three column vertex array.

I am puzzled by this chiefly because my assumption was that the vertices used to create the isosurface had to come directly from the values of the position matricies input into the original isosurface function, but that doesn't seem to have been the case here.

Any help with this, or alternative methods to obtaining the values/indices of omega that correspond with the vertex values of the isosurface are greatly appreciated. Thank you in advance.

Cheers,

John

[KSSV]

Already you have the required vertices of the surface in S. Why you want to find them again?

[John Kearns]

Because the vertex values does not store the value of the field from which they are obtained. That is, I want the values of omega at that point.

[KSSV]

Show us your full code.