Convert label map volume into a scalarvolumenode

Laura · July 21, 2018, 11:41am

Good morning,

To convert my label map volume into a scalarvolumenode after correction on Segment editor and “export” the label map, I have some difficulties…
I am able to convert this updated label map into a numpy array (called ‘a’) but then, I don’t succeed in converting it into a vtkMRMLScalarVolumeNode : do you know how can I do this ?
Maybe I have committed some errors before by exporting my label map with my segmentation…

This is my code :

seg = slicer.util.getNode(‘CorrectionSegmentationFoie’)
labelmapVolumeNode = slicer.mrmlScene.AddNewNodeByClass(‘vtkMRMLLabelMapVolumeNode’,‘LabelFoiefinal’)

foiesegfin = slicer.modules.segmentations.logic().ExportAllSegmentsToLabelmapNode(seg, labelmapVolumeNode)
a = slicer.util.array(‘LabelFoiefinal’)
volumeNode = slicer.mrmlScene.AddNewNodeByClass(‘vtkMRMLScalarVolumeNode’,“VolumeFoieSegmente”)
foieseg = slicer.util.updateVolumeFromArray(volumeNode, a)

Thanks a lot in advance !
Laura

P.S : Do you know how to force in my python script to check the “sphere brush” in my 'Segment editor" module and force also my “MasterVolume” to be the label map that I have created from my previous uncorrected volume ?

Fernando · July 22, 2018, 11:02am

Hi Laura,

You can try CreateScalarVolumeFromVolume .

Laura · July 22, 2018, 1:32pm

Hey !

Thanks for this idea, I have tried so much but not this one !
Unfortunately, when I write this on my python script :

seg = slicer.util.getNode('CorrectionSegmentationFoie')
labelmapVolumeNode = slicer.mrmlScene.AddNewNodeByClass('vtkMRMLLabelMapVolumeNode','LabelFoiefinal')
slicer.modules.segmentations.logic().ExportAllSegmentsToLabelmapNode(seg, labelmapVolumeNode)
outputvolumenode = slicer.mrmlScene.AddNewNodeByClass("vtkMRMLScalarVolumeNode", 'FoieSegmenteApresCorrection')
sef = slicer.modules.volumes.logic().CreateScalarVolumeFromVolume(slicer.mrmlScene,labelmapVolumeNode, outputvolumenode)

My 3D Slicer crashes without saying nothing… Just crashes and asks me if I want to reload Slicer without any explanation…
Do I made a mistake ?
Thanks a lot in advance and enjoy your weekend !
Laura

cpinter · July 22, 2018, 3:13pm

I don’t see any problems with your code after a quick look. A few things to check:

What does the log say? You can find it in About menu / Report a problem window
If you show labelmapVolumeNode in the slice view or in 3D volume rendering do you see anything?
If you enter the lines one by one which one causes the crash?

Fernando · July 22, 2018, 4:44pm

I just tried her code and Slicer crashed after running the last command. The last line of my log is not suspicious:

[DEBUG][Qt] 22.07.2018 17:42:34 [] (unknown:0) - Python console user input: outputvolumenode = slicer.mrmlScene.AddNewNodeByClass("vtkMRMLScalarVolumeNode", 'FoieSegmenteApresCorrection')

Also, the label map looks fine.

Fernando · July 22, 2018, 4:52pm

Try swapping the arguments:

sef = slicer.modules.volumes.logic().CreateScalarVolumeFromVolume(slicer.mrmlScene, outputvolumenode, labelmapVolumeNode)

According to the docs, the second argument is your output and the third is the input.

But of course Slicer shouldn’t just crash

cpinter · July 22, 2018, 7:41pm

Thanks, Fernando!
Does this function crash with any data? If so, then it should be an easy fix. Otherwise I can fix it if you give me the labelmap that you get after exporting.

Fernando · July 22, 2018, 7:53pm

Anything should work. I just thresholded MRHead in the Segment Editor, renamed the segmentation node and pasted the code.

cpinter · July 22, 2018, 8:19pm

I fixed the crash, I’ll integrate it in a minute so it shouldn’t crash in tomorrow’s version.

There is also something to fix in VTK too, though, because DeepCopy shouldn’t crash if you give it null pointer, which is what happens here (so my fix in Slicer is only a superficial fix). This is where it starts (and I made sure it doesn’t get this far):

github.com

Slicer/Slicer/blob/master/Modules/Loadable/Volumes/Logic/vtkSlicerVolumesLogic.cxx#L998


  }


// Associate labelmap with the source volume
if (inputVolume->GetID())
  {
  outputVolume->SetNodeReferenceID("AssociatedNodeID", inputVolume->GetID());
  }


// Copy and set image data of the input volume to the label volume
vtkNew<vtkImageData> imageData;
imageData->DeepCopy(inputVolume->GetImageData());
outputVolume->SetAndObserveImageData(imageData.GetPointer());


vtkNew<vtkMatrix4x4> ijkToRas;
inputVolume->GetIJKToRASMatrix(ijkToRas.GetPointer());
outputVolume->SetIJKToRASMatrix(ijkToRas.GetPointer());


outputVolume->SetAndObserveTransformNodeID(inputVolume->GetTransformNodeID());


// Create a display node if the output node does not have one
outputVolume->CreateDefaultDisplayNodes();

Then vtkImageData and vtkDataSet handles the null properly, but then in vtkDataObject this call causes the crash ultimately:

github.com

Kitware/VTK/blob/master/Common/DataModel/vtkDataObject.cxx#L550


    fd->ShallowCopy(src->FieldData);
    this->SetFieldData(fd);
    fd->FastDelete();
  }
}
}


//----------------------------------------------------------------------------
void vtkDataObject::DeepCopy(vtkDataObject *src)
{
vtkFieldData *srcFieldData = src->GetFieldData();


this->InternalDataObjectCopy(src);


if (srcFieldData)
{
  vtkFieldData *newFieldData = vtkFieldData::New();
  newFieldData->DeepCopy(srcFieldData);
  this->SetFieldData(newFieldData);
  newFieldData->FastDelete();
}

Optimally this should be fixed in VTK.

To answer the original question, I think Fernando’s suggestion is the solution, which is to swap the arguments.

lassoan · July 22, 2018, 8:32pm

Unfortunately, DeepCopy has no return value and VTK does not use exceptions, so there is not really a clean way to handle invalid DeepCopy input. Maybe VTK could log an error and return when the input is NULL instead of crashing? There is a chance that the error would be missed, but probably it is still better than crashing.

cpinter · July 22, 2018, 8:35pm

In the DeepCopy functions of vtkImageData and vtkDataSet the null argument was simply handled with an if statement. Indeed, no error was logged. I think adding it would be a nice solution, but then to make it consistent we’d need to change potentially a hundred VTK classes.
Maybe the easiest would be to just fix this one in vtkDataObject.

cpinter · July 22, 2018, 8:36pm

By the way, @Laura, swapping the arguments was the solution to your problem, I could confirm. So the last line should be

sef = slicer.modules.volumes.logic().CreateScalarVolumeFromVolume(slicer.mrmlScene, outputvolumenode, labelmapVolumeNode)

Laura · July 23, 2018, 6:19am

Good morning,

Thank you so so much to everybody, this works perfectly !!!
Have a good day
Laura

Laura · July 23, 2018, 12:13pm

Another little question @cpinter and @Fernando : do you know how to force in my python script to check the “sphere brush” in my 'Segment editor" module and force also my “MasterVolume” to be the label map that I have created from my previous uncorrected volume ?

Thanks again for everything !

Laura

cpinter · July 23, 2018, 1:32pm

Hi Laura,

You can find information about driving Segment Editor from python here:
https://www.slicer.org/wiki/Documentation/Nightly/ScriptRepository#How_to_run_segment_editor_effects_from_a_script

I think if you set the reference volume using the function called SetReferenceImageGeometryParameterFromVolumeNode in the segmentation node, then it should be selected when you enter Segment Editor. Please note that you need to set this as soon as you have both nodes in the scene, and definitely before you enter Segment Editor.

lassoan · July 23, 2018, 7:00pm

You can also explicitly call segmentEditorWidget.setMasterVolumeNode(...) as it is done in almost all segmentation scripting examples.

It is a parameter effect, which you can modify by calling effect.SetParameter(…) as shown in examples, or modifying the segment editor node directly: getNode('SegmentEditor').SetAttribute('BrushSphere','1')

cpinter · July 23, 2018, 7:19pm

That’s right, you’re already using the Segment Editor widget from script anyway, so this will make it simpler.

Laura · July 24, 2018, 7:33am

Thanks a lot for your super help !
I succeed for the brush but I am still looking for the ‘setMasterVolumeNode’ because, when I write :

editorWidget = slicer.modules.segmenteditor.createNewWidgetRepresentation()
editorWidget.setMRMLScene(slicer.app.mrmlScene())
slicer.util.getNode(‘SegmentEditor’).SetAttribute(‘BrushSphere’,‘1’)
editorWidget.setMasterVolumeNode(inputVolCor)
editorWidget.show()

I have this error caused by my last line :

AttributeError: qSlicerScriptedLoadableModuleWidget has no attribute named ‘setMasterVolumeNode’

But I am looking for a solution !
Thanks again !
Laura

lassoan · July 24, 2018, 5:23pm

This line does not seem to be correct. You should not instantiate a new module widget.
You either create a new segment editor widget (as is shown in most examples; I would recommend this):

editorWidget = slicer.qMRMLSegmentEditorWidget()

or access the existing module widget of the SegmentEditor module:

slicer.util.selectModule("SegmentEditor")
editorWidget=slicer.modules.segmenteditor.widgetRepresentation().self().editor

Laura · July 25, 2018, 7:58am

Good morning,

I agree with you but this is the only way I found to create a popup of Segment Editor :

for your first solution :

editorWidget = slicer.qMRMLSegmentEditorWidget()

No segment editor popup appears even when I use the script in most examples with the line ‘.show()’ and I would like to paint and erase as I need so this is not convenient for what I do

for your second solution :

slicer.util.selectModule(“SegmentEditor”)
editorWidget=slicer.modules.segmenteditor.widgetRepresentation().self().editor

The Segment Editor Widget appears of course but not in a popup so I haven’t my extension wizard background so this is not convenient neither

That’s why I have found this one :

editorWidget = slicer.modules.segmenteditor.createNewWidgetRepresentation()

Topic		Replies	Views
How to convert a vtkMRMLScalarVolumeNode to a vtkMRMLSegmentationNode Support segmentation , python , volume , binarylabelmap	7	1008	January 21, 2022
Segmentation to binary label map through python console Support	6	2413	August 21, 2024
How to convert volume to label map programmatically Support	3	885	March 3, 2022
From labelmap to model and from model to labelmap python scripts Support	2	935	February 4, 2021
Get label map node from segmentation node Development segmentation	14	6118	November 26, 2019

Convert label map volume into a scalarvolumenode

Related topics