Calculating distance using Maurer distance map

Support
14

Here you go:

import vtkSlicerSegmentComparisonModuleLogicPython
polHist = vtkSlicerSegmentComparisonModuleLogicPython.vtkPolyDataDistanceHistogramFilter()
segNode = slicer.util.getNode('vtkMRMLSegmentationNode1')
poly1 = segNode.GetSegmentation().GetSegment('Segment_1').GetRepresentation(slicer.vtkSegmentationConverter.GetSegmentationClosedSurfaceRepresentationName())
poly2 = segNode.GetSegmentation().GetSegment('Segment_2').GetRepresentation(slicer.vtkSegmentationConverter.GetSegmentationClosedSurfaceRepresentationName())
polHist.SetInputReferencePolyData(poly1)
polHist.SetInputComparePolyData(poly2)
polHist.Update()
histTable = polHist.GetOutputHistogram()
tableNode = slicer.vtkMRMLTableNode()
tableNode.SetAndObserveTable(histTable)
slicer.mrmlScene.AddNode(tableNode)

Then you can save the table to csv with Save data dialog, or oyu can make a plot with it from Tables module.

Oh, and you need to install the SlicerRT extension first to have the Segment Comparison module if you haven’t.