output labels or the data of "elem(:,5)" from surf2mesh about iso2mesh HOT 2 CLOSED

@tmedani, sorry for the delay.

the reason you got 0 labels in some of the regions was because your seeds were not generated properly.

because in this case, 3 of the surfaces are enclosed by the outer layer, therefore, raysurf.m may not return the seed that is specific for a given region.

raysurf.m works by shooting a ray across a closed surface, finding a pair of intersection points, and use the mid-point as the seed. This only works when a closed surface is not enclosed/enclosing another surface - because the mid-point may be inside the surface being enclosed.

to do this properly, you need to use something like below:

%% define seeds along the electrode axis
[t,baryu,baryv,faceidx]=raytrace(orig,v0,allnode,allface);
t=sort(t(faceidx));
t=(t(1:end-1)+t(2:end))*0.5;
seedlen=length(t);
seeds=repmat(orig(:)',seedlen,1)+repmat(v0(:)',seedlen,1).*repmat(t(:)-1,1,3);

where the allnode/allface are the combined surfaces of the 3 layers. The raytrace call returns all intersection points with distance from the ray origin (orig) stored in t. So, the mid-point of every segment between adjacent intersection points defines a seed that belongs to a specific region.

Please give above code a try and let me know if you need me to explain this further.

from iso2mesh.

Dear @fangq
Thank you for your answer and no worry for the delay ;)

I tried your solution, it works fine, thanks.

However, since I have three surfaces (in this example), I tried just to use the 3 upper points,
in this case, the label 0 is still here.
So I took the whole 6 points (length(t)) and the outputs are perfect.

I changed the last line to :
seeds=repmat(orig(:)',seedlen,1)+repmat(v0(:)',seedlen,1).*repmat(t(:),1,3);

I removed the (-1)

Thank you again for your help.

Kind regards
Takfarinas

from iso2mesh.