SurfaceMesh generating in FemurCut

[Han_Lin]

Accoring to suggestions about how to compute an SDF:

LevelSetModel construction change

Hello,I have read the FemurCut example ,in which you first compute the SDF of a leg nii and then modify their voxels.But Could you give me some advice or examples on how to modify a SDF of a SurfaceMesh(such as .stl)? Is there a filter available in the current library? Thanks a lot!

I tried to building one by my own nii pixel data.
That’s how i converts a NIfTI 3D voxel image to a triangulated mesh

vtkSmartPointer<vtkDiscreteMarchingCubes>discrete = vtkSmartPointer<vtkDiscreteMarchingCubes>::New();
discrete->SetInputConnection(niiReader->GetOutputPort());
iscrete->Update();

vtkSmartPointer<vtkWindowedSincPolyDataFilter> smoother = vtkSmartPointer<vtkWindowedSincPolyDataFilter>::New();
smoother->SetInputConnection(discrete->GetOutputPort());
smoother->SetNumberOfIterations(100);

Without generating .obj data, i try to directly convert mesh from vtkpolydata to surfacemesh.
so i displace MeshIO::read with imstk::GeometryUtils::copyToSurfaceMesh like this

std::shared_ptr<imstk::SurfaceMesh> spinalSurfMesh = imstk::GeometryUtils::copyToSurfaceMesh(smoother->GetOutput());

Results shows that part of surface becomes shattered and i can’t get properly collision(tool will disappear) but it goes well when i generate a obj file and reading it by MeshIO::read.

What’s different bewteen two operators or can any suggestions provide?

[Han_Lin]

imstk::GeometryUtils::copyToSurfaceMesh is effective indeed.

Problem also caused by normals-flipped when generating surfacemesh

The collision working good after do SurfaceMesh->flipNormals().

[Harald_Scheirich]

Great to see, does that solve your problem ?

[Han_Lin]

it did. Thanks for your attention