Dear community,
I use this segment of code frequently to make 3D models using the electrical impedance tomography (EIT) package EIDORS. I am producing a multi-layered head model by subtracting the overlapping tetrehedra of STL files of skull, brain, etc. from meshed head models using the MATLAB functions find() and any().
How can I speed up this function given that the (number of) overlapping elements for each layer depend(s) on the geometries of each individual model?
fwd_mdlvent = stl_read(STLfilenames{v});
elems_with_nodes_inside_vent_border = [];
node_row_value = intriangulation(fwd_mdlvent.nodes,fwd_mdlvent.elems,fwd_mdlHeadWithElec.nodes);
node_row_value_if_inside_border = find(node_row_value);
for k = 1:size(node_row_value_if_inside_border,1)
elems_with_this_node_value = find(any(fwd_mdlHeadWithElec.elems==node_row_value_if_inside_border(k),2));
elems_with_nodes_inside_vent_border = [elems_with_nodes_inside_vent_border; elems_with_this_node_value];
elems_with_nodes_inside_vent_border = unique(elems_with_nodes_inside_vent_border);
elems_with_nodes_inside_brain_border = setdiff(elems_with_nodes_inside_brain_border,elems_with_nodes_inside_vent_border);
save(filename,'headImageCases');
The result is a nested struct with a forward model of the layered-head with shell elements (points) separated by material index.
I hope this explains my problem, though I know it is highly specific. Thank you in advance and I look forward to discussing further with you.
-Jared
**PS - I recognize that it may be helpful for me to provide the full code for this process in the case that my explanations are unclear, but the models and required files make this quite complicated so I hoped that someone may be able to help with just the provided snippet. If including everything will help you to address my problem then please leave a comment and I will take the time to upload/explain the required files and full code.
