Hi all,
I want to export a volumetric model along with HU values for FEA analysis.
Based on several very helpful forum posts, I was able to create a volumetric mesh (single material) with the SegmentMesher module. I then used Probe Volume with Model to get a map of the densities across the model.
I have 2 questions.
1) export of volumetric mesh (single material)
Unfortunately my FEA program of choice (Artisynth) is unable to read the vtk FE mesh file. I get an error that it is of type “Polydata” and unable to read that (vtk ascii should be supported).
To double check that it isn’t an error with the FE program, I also tried using meshio in the Slicer Python console to try to read the vtk mesh and convert to .inp (which Artisynth can also read):
import meshio
mesh = meshio.read(“T:/Slicer_Tests/Material Mapping/Test_single_material_skull.vtk”)
meshio.write(“T:/Slicer_Tests/Material Mapping/test2.inp”, mesh)
However, when trying to read the .vtk file, I also get the error:
Only VTK ‘UNSTRUCTURED_GRID’, ‘STRUCTURED_POINTS’, ‘STRUCTURED_GRID’, ‘RECTILINEAR_GRID’ supported (not POLYDATA).
And when I open the .vtk file itself it says “DATASET POLYDATA”.
My question is whether this is a bug in the .vtk export (I did select “unstructured grid”, not “polydata”).
I also get the error vtkMRMLModelDisplayNode (00000218420FB2D0): Can not use color node scalar range since model display node does not have a color node. during the export, and it looks like in the type of export, it quickly flashes to “polydata” instead of “unstructured grid” before the pop up window closes.
2) mapping material properties
Related to question 1, when I export the output from Probe Volume with Model, I get a polydata .vtk file. Looking at the file it is unclear to me how it is structured, e.g. where the coordinates are and where the HU values are.
Thank you very much!
Eva