ENH: add createDentalImplantCylindersFiducialList, createCylindersFromFiducialListAndMandibleReconstruction, onDentalImplantPlaneMoved, updateFibulaDentalImplantCylinders functions; add fibulaDentalImplantCylinders to boolean operations

Author: mauigna06

BoneReconstructionPlanner/BRPLib/helperFunctions.py

@@ -101,15 +101,15 @@ def getAverageNormalFromModel(model):
 
   return averageNormal
 
-def getAverageNormalFromModelPoint(model,point):
+def getAverageNormalFromModelPoint(model,point,radius=3):
   normalsOfModel = slicer.util.arrayFromModelPointData(model, 'Normals')
 
   modelMesh = model.GetMesh()
   pointID = modelMesh.FindPoint(point)
   normalAtPointID = normalsOfModel[pointID]
 
   cylinder = vtk.vtkCylinder()
-  cylinder.SetRadius(3)
+  cylinder.SetRadius(radius)
   cylinder.SetCenter(point)
   cylinder.SetAxis(normalAtPointID)
 
@@ -175,19 +175,21 @@ def createBox(X, Y, Z, name):
   miterBox.SetAndObservePolyData(triangleFilter.GetOutput())
   return miterBox
 
-def createCylinder(R,name):
+def createCylinder(name,R,H=50):
   cylinder = slicer.mrmlScene.CreateNodeByClass('vtkMRMLModelNode')
   cylinder.SetName(slicer.mrmlScene.GetUniqueNameByString(name))
   slicer.mrmlScene.AddNode(cylinder)
   cylinder.CreateDefaultDisplayNodes()
   lineSource = vtk.vtkLineSource()
-  lineSource.SetPoint1(0, 0, 25)
-  lineSource.SetPoint2(0, 0, -25)
+  lineSource.SetPoint1(0, 0, H/2)
+  lineSource.SetPoint2(0, 0, -H/2)
   tubeFilter = vtk.vtkTubeFilter()
   tubeFilter.SetInputConnection(lineSource.GetOutputPort())
   tubeFilter.SetRadius(R)
   tubeFilter.SetNumberOfSides(50)
   tubeFilter.CappingOn()
   tubeFilter.Update()
   cylinder.SetAndObservePolyData(tubeFilter.GetOutput())
+  cylinder.SetAttribute('radius',str(R))
+  cylinder.SetAttribute('height',str(H))
   return cylinder