Region Growth Over Contour on STL Surface

christeen · January 9, 2025, 6:39am

I used the vtkPolygonalSurfacePointPlacer and vtkPolygonalSurfaceContourLineInterpolator to draw a contour over the surface of an actor using vtkContourWidget. Now, I want to extract the upper portion of the contour on the STLactor and apply a region growing effect similar to the following:

into Desired outcome as, :

I am looking for guidance on how to proceed with this development, as I am new to this area. Your support in helping me achieve this would be greatly appreciated. I’m currently using VTK 8.2 and C++ for this task.

lassoan · January 9, 2025, 2:11pm

Thise feature is already implemented in 3D Slicer, which is a VTK-based open-source application, therefore you can take the full implementation from there:

Source code for curve cut tool:

github.com

Slicer/SlicerSurfaceToolbox/blob/e8b8f70930883adb6f4a227ad9d7339d20120f2c/DynamicModeler/Logic/vtkSlicerDynamicModelerCurveCutTool.cxx#L208-L315


      
          double bounds_World[6] = {0.0};
          inputModelNode->GetRASBounds(bounds_World);
          
          this->CleanFilter->SetInputData(inputModelNode->GetPolyData());
          if (inputModelNode->GetParentTransformNode())
            {
            inputModelNode->GetParentTransformNode()->GetTransformToWorld(this->InputModelToWorldTransform);
            }
          else
            {
            this->InputModelToWorldTransform->Identity();
            }
          
          this->SelectionFilter->SetLoop(curveNode->GetCurvePointsWorld());
          
          vtkMRMLMarkupsFiducialNode* insidePointNode = vtkMRMLMarkupsFiducialNode::SafeDownCast(
            surfaceEditorNode->GetNthNodeReference(CURVE_CUT_INPUT_INSIDE_POINT_REFERENCE_ROLE, 0));
          if (insidePointNode && insidePointNode->GetNumberOfControlPoints() > 0)
            {
            double insidePointPositionWorld[3] = { 0.0, 0.0, 0.0 };

This file has been truncated. show original

Your experimentation with VTK is extremely valuable learning experience, but when it comes to building the software that you ship to your customers (researchers, dentists, etc.) then I would strongly recommend to move to a more full featured VTK-based medical application platform, such as 3D Slicer. This platform already provides all other essential features for a dental implant planning software, such as convenient DICOM browser, automatic AI dental segmentation tools, panoramic image reconstruction, curve and thick planar reformat views, integration with several surgical navigation systems, tool trackers, real-time imaging devices, etc. All these are available for completely free, without any restrictions, for any use, including using it as a basis for commercial, FDA-approved products. By building on a platform you save millions of dollars worth of software development and maintenance, so you can focus your efforts on things that make a difference (unique algorithms, ideas, workflows), not on redeveloping yet another DICOM browser or 2D/3D image viewer.

There is even a full dental implant planning suite implemented in 3D Slicer that you might potentially leverage if you have some VTK-based tools to deliver:

This application is not open-source, so if you are interested in starting from this then you would need to contact its developers (see their names for example on this recent 3D Slicer Project Week page on automatic dental segmentation.

christeen · January 10, 2025, 6:34am

Thank you so much sir…

christeen · January 23, 2025, 10:19am

I have a requirement as per continuation of Region growth…Now I want to make region growth like as like the following image

and I would be like to clip the some portions exactly as per picture given below…

How can able to consider these operations using VTK 8.2 with cpp…Currently I done region growing using vtkCleanPolyData(input polydata as STLactor and loop as contourwidget) to clipping procedure for region growing…Now I want to thicken the region grown area and clip some of the portions as accurate…Which is the best choice for that…vtkpolydata / mesh or any other methods…

Currently I have done process by making layers by using vtkCellArray and vtkPoints with the help of vtkAppendPolyData…I’ll give that output…

I have to improve the requirement…I am looking for guidance on how to proceed with this development, as I am new to this area. Your support in helping me achieve this would be greatly appreciated. I’m currently using VTK 8.2 and C++ for this task.

lassoan · January 23, 2025, 3:58pm

There are many questions, hard to decode them, and not directly related to the original question. It would be easier if you posted each new question as separate topic, with detailed description and illustrations.