Isodose model is only 1 file for all of them

Support
1

Dear Slicer Users,
I have Slicer 5.6.0 now with Radiotherapy module. When I generate isodoses from RTdose dicom file, only 1 model for all isodoses is created (looks like sum of all).
In previous version of Slicer it was possible to see each isodose separately and also to export isodoses to models separately.
I am not sure if I have a chance to do something differently or there is no longer possibility to generate isodoses separately.

Thank you!

Lukas

2

It has been a conscious decision to change the output of the Isodose module to a single model, where the isodose levels are indicated as scalars in that model. If you need separate models, you can use VTK to separate the surfaces for each scalar value. You can use a function like this:

#----------------------------------------------------------------------------
def getPolyDataFromScalarValue(modelPolyData, targetScalarValue, scalarName):
  """
  Get polydata based on the value of the scalar array.
  :param modelPolyData: input polydata
  :param scalarValue: points with this scalar value will be kept
  :param scalarName: name of the scalar array
  :return vtkPolyData: output polydata
  """
  # Select points in polydata
  selectedPointIds = vtk.vtkIdTypeArray()
  numPoints = modelPolyData.GetNumberOfPoints()
  modelPointValues = modelPolyData.GetPointData().GetArray(scalarName)
  for pointId in range(numPoints):
    scalarValue = modelPointValues.GetTuple(pointId)[0]
    if scalarValue == targetScalarValue:
      loc = selectedPointIds.InsertNextValue(pointId)
  # Extract selected points
  selectionNode = vtk.vtkSelectionNode()
  selectionNode.SetFieldType(vtk.vtkSelectionNode.POINT)
  selectionNode.SetContentType(vtk.vtkSelectionNode.INDICES)
  selectionNode.SetSelectionList(selectedPointIds)
  selectionNode.GetProperties().Set(vtk.vtkSelectionNode.CONTAINING_CELLS(), 1)
  selection = vtk.vtkSelection()
  selection.AddNode(selectionNode)
  extractSelection = vtk.vtkExtractSelection()
  extractSelection.SetInputData(0, modelPolyData)
  extractSelection.SetInputData(1, selection)
  extractSelection.Update()
  convertToPolydata = vtk.vtkDataSetSurfaceFilter()
  convertToPolydata.SetInputConnection(extractSelection.GetOutputPort())
  convertToPolydata.Update()
  return convertToPolydata.GetOutput()
3

Thanks a lot for this! I will try.

Best regards

Lukas