I have two VTK mesh models in 3D Slicer as shown in the figure. One of them is the brain model of the MRI data which is rendered using the Marching Cubes algorithm while the other model in color is the model generated within the intraoperative trepanation boundary which is the exposed brain surface during the surgery. As you see, the intraoperative model is slightly above the brain surface. This is because the intraoperative points were acquired from the skull surface and not brain surface. Now I just want the points of the intersection of intraoperative model and the brain surface something like shown in the second figure. Could anyone let me know how I can achieve this in 3D Slicer? Thanks!
You can import both meshes into a Segmentation node, then use Logical operators effect in Segment Editor to compute intersection. You can use Margin effect on the intraoperative model if you want to increase the amount of overlap between the segments.
Thanks for your reply. Actually I tried using segment editor with the intersect option in Logical operator however I cannot get the image as shown in my question post. I think I missed my first figure in my previous post. The image below is what I have currently but I want to compute the intersection of the both model and just get the surface surrounding it.
What is the “both model”?
What do you mean by “surface surrounding it”? Segments are represented as surfaces as well and you can export segment to surface model using Export/Import section in Segmentations module.
Maybe you could attach a few more screenshots of what you would like to have as an end result (from different views, with different display settings) because from the ones above it is not yet clear.
Dear Andras Lasso,
Sorry for not making myself clear. Let me try to explain the problem in a simpler manner:
Well as you see in the left side of the below figure, I have an intraoperative mesh model which is orange in color and brain mesh model which is showcased in green. After using the Intersect + Smoothing operations of Segment editor, I was successfully able to extract the intersection of the models giving me a new intraoperative model as shown in the right figure. As you see, the new intraoperative model goes till the tip of the cone starting from the surface of the brain model. However, I would like to extract only the vertices of the intraoperative model which is at the surface of the brain model.
Secondly, after I have extracted vertices at the surface of the brain model, I would like to find out the MRI intensities at those vertices so that I could interpolate those intensities on a 2D image grid.
You may use Mask volume effect (provided by SegmentEditorExtraEffects extension) to set MRI intensities to -1000 outside the resection region, then use “Probe volume with model” module on the masked volume with the full brain surface model to get image intensities of the brain surface. You can set up display of probe volume output in Modules module to hide very low intensities (so the regions that were set to -1000 will be cut away).