Create new vtkMRMLModelHierarchyNode() and add models

Fernando · May 30, 2017, 12:37pm

Hi slicers,

This is a simple question:
I have a modelNodeA and modelNodeB.
Using Python, I would like to create a modelHierarchyNode and add my two model nodes inside. Can you please give me a hint?

Best,
Fernando

cpinter · May 30, 2017, 1:09pm

Hi Fernando,

This is how I do it in C++:

  vtkSmartPointer<vtkMRMLModelHierarchyNode> rootModelHierarchyNode = vtkSmartPointer<vtkMRMLModelHierarchyNode>::New();
  scene->AddNode(rootModelHierarchyNode);
  rootModelHierarchyNode->SetName(rootModelHierarchyNodeName.c_str());
  rootModelHierarchyNode->SetSingletonTag(rootModelHierarchyNodeName.c_str());

  if (!rootModelHierarchyNode->GetDisplayNode())
  {
    vtkSmartPointer<vtkMRMLModelDisplayNode> rootModelHierarchyDisplayNode = vtkSmartPointer<vtkMRMLModelDisplayNode>::New();
    scene->AddNode(rootModelHierarchyDisplayNode);
    rootModelHierarchyNode->SetAndObserveDisplayNodeID( rootModelHierarchyDisplayNode->GetID() );
  }

...

  vtkSmartPointer<vtkMRMLModelHierarchyNode> collimatorModelHierarchyNode = vtkSmartPointer<vtkMRMLModelHierarchyNode>::New();
  scene->AddNode(collimatorModelHierarchyNode);
  collimatorModelHierarchyNode->SetModelNodeID(collimatorModelNode->GetID());
  collimatorModelHierarchyNode->SetParentNodeID(rootModelHierarchyNode->GetID());
  collimatorModelHierarchyNode->HideFromEditorsOn();

lassoan · May 30, 2017, 1:36pm

Are you sure you want to build model hierarchy? What is your use case? You may be better off using subject hierarchy instead, as subject hierarchy is:

more robust: model hierarchy tree widget has some issues, for example sometimes items randomly change place within a branch
has much simpler API
has more features: you can define custom actions, you can save the order of nodes in the hierarchy into the scene, and you can group all kinds of related nodes in a branch, not just models
much faster

For example, to create a branch and put a node in it:

shNode = slicer.vtkMRMLSubjectHierarchyNode.GetSubjectHierarchyNode(slicer.mrmlScene)
newParentFolderItem = shNode.CreateFolderItem(shNode.GetSceneItemID(), "my folder")
volumeItem = shNode.GetItemByDataNode(getNode('MRHead'))
shNode.SetItemParent(volumeItem, newParentFolderItem)

Fernando · May 30, 2017, 1:56pm

Hi Csaba and Andras,

Thanks for your quick answers.

@cpinter My Python version:

modelNodeA = getNode('4648J_xscp_model')
modelNodeB = getNode('4648J_RH_lq_model')

rootModelHierarchyNodeName = 'My hierarchy'

rootModelHierarchyNode = slicer.vtkMRMLModelHierarchyNode()
slicer.mrmlScene.AddNode(rootModelHierarchyNode)
rootModelHierarchyNode.SetName(rootModelHierarchyNodeName)
rootModelHierarchyNode.SetSingletonTag(rootModelHierarchyNodeName)  # Not sure what this is

rootModelHierarchyDisplayNode = slicer.vtkMRMLModelDisplayNode()
slicer.mrmlScene.AddNode(rootModelHierarchyDisplayNode)
rootModelHierarchyNode.SetAndObserveDisplayNodeID(rootModelHierarchyDisplayNode.GetID())

modelAHierarchyNode = slicer.vtkMRMLModelHierarchyNode()
slicer.mrmlScene.AddNode(modelAHierarchyNode)
modelAHierarchyNode.SetModelNodeID(modelNodeA.GetID())
modelAHierarchyNode.SetParentNodeID(rootModelHierarchyNode.GetID())
modelAHierarchyNode.HideFromEditorsOn()  # Not sure what this is

modelBHierarchyNode = slicer.vtkMRMLModelHierarchyNode()
slicer.mrmlScene.AddNode(modelBHierarchyNode)
modelBHierarchyNode.SetModelNodeID(modelNodeB.GetID())
modelBHierarchyNode.SetParentNodeID(rootModelHierarchyNode.GetID())
modelBHierarchyNode.HideFromEditorsOn()

What I expect:

What I get:

I had already tried with some code I found (probably yours), and got similar results.

@lassoan I have a lot of models and I want to divide them in groups so that the user can easily change their display properties, using Slicer default modules if possible. I know there are some issues though…
I have never taken the time to look into subject hierarchy. It does looks simpler. Running your snippet makes a folder, from which I cannot modify the children display properties (at least I don’t know how to do it).

cpinter · May 30, 2017, 2:23pm

modelAHierarchyNode.HideFromEditorsOn() # Not sure what this is

This call is supposed to make the leaf hierarchy nodes disappear from the tree (to get from what you get to what you expect). Not sure why it doesn’t work, try SetHideFromEditors(1)

Fernando · May 30, 2017, 2:59pm

SetHideFromEditors(1) worked as expected. I tried to dive into the code to solve the issue I commented before, but I didn’t manage. Do you think it’s hard to fix?

Thank you again

cpinter · May 30, 2017, 4:20pm

You mean the slice intersection issue you linked?

I quickly took a look, and it seems that the checkbox in the Models module that applies color/opacity on the branch drives the vtkMRMLDisplayableHierarchyNode::Expanded flag, which is handled by the 3D displayable manager in vtkMRMLModelDisplayableManager::SetModelDisplayProperty, when getting the hierarchy display node (line 1488: vtkMRMLDisplayNode *hierarchyDisplayNode = this->GetHierarchyDisplayNode(model) ). This call is responsible for getting the right properties. The 2D displayable manager (vtkMRMLModelSliceDisplayableManager) does not use this, and this is what needs to be fixed.

I don’t have time to fix it, but I hope this will help. I’m adding this info on Mantis as well.

Fernando · May 30, 2017, 5:43pm

Yes that’s what I meant. Thank you!

brhoom · May 14, 2020, 7:25pm

I created two model trees using this approach. I can access each tree and its children names. How can I get a model node from a specific child? e.g. to get their display node, polydata , …etc

  self.shNode = slicer.vtkMRMLSubjectHierarchyNode.GetSubjectHierarchyNode(slicer.mrmlScene)
  self.modelTreeID = self.shNode.GetItemChildWithName(self.shNode.GetSceneItemID(), "modelTreeName" ) 
  self.modelTreeChildrenIDs = vtk.vtkIdList()
  self.shNode.GetItemChildren(self.modelTreeID, self.modelTreeChildrenIDs)
  for i in xrange(self.modelTreeChildrenIDs GetNumberOfIds()):
      print("modelChildID   : ", self.modelTreeChildrenIDs.GetId(i)) 
      print("modelChildName : ", self.shNode.GetItemName(self.modelTreeChildrenIDs.GetId(i)))

lassoan · May 15, 2020, 1:41am

Model hierarchy is replaced by the much more generic subject hierarchy in recent Slicer versions. It should be easier to use and has many more capabilities. Script repository contains examples to most common tasks.

brhoom · May 15, 2020, 7:15am

@lassoan thanks for your quick reply.

Model hierarchy is replaced by the much more generic subject hierarchy

I am using subject hierarchy, kindly see the code above.

. Script repository contains examples to most common tasks.

I already checked that. It is not clear to me how to get a “model node” from a child. I can get the child ID e.g. modelTreeChildrenIDs.GetId(i) or the chile name e.g. shNode.GetItemName(self.modelTreeChildrenIDs.GetId(i)) but not the model node.

cpinter · May 15, 2020, 8:35am

You can get the model node using this function

github.com

Slicer/Slicer/blob/master/Libs/MRML/Core/vtkMRMLSubjectHierarchyNode.h#L115



// Get/Set methods
public:
  /// Get the (practically) singleton subject hierarchy node from MRML scene.
  /// Merges subject hierarchy nodes if multiple found, and returns the merged one.
  static vtkMRMLSubjectHierarchyNode* GetSubjectHierarchyNode(vtkMRMLScene* scene);

  /// Get ID of root subject hierarchy item (which can be interpreted as the scene in terms of hierarchy)
  vtkIdType GetSceneItemID();
  /// Get data node for a subject hierarchy item
  vtkMRMLNode* GetItemDataNode(vtkIdType itemID);
  /// Set name for a subject hierarchy item
  void SetItemName(vtkIdType itemID, std::string name);
  /// Get name for a subject hierarchy item
  /// \return Name of the associated data node if any, otherwise the name of the item
  std::string GetItemName(vtkIdType itemID);
  /// Convenience function to set level attribute for a subject hierarchy item
  void SetItemLevel(vtkIdType itemID, std::string level);
  /// Convenience function to get level attribute for a subject hierarchy item
  std::string GetItemLevel(vtkIdType itemID);
  /// Set owner plugin name (role) for a subject hierarchy item

brhoom · May 15, 2020, 10:23am

@cpinter Thanks a lot , this works, Here is the complete code I used:

self.shNode = slicer.vtkMRMLSubjectHierarchyNode.GetSubjectHierarchyNode(slicer.mrmlScene)
self.modelTreeID = self.shNode.GetItemChildWithName(self.shNode.GetSceneItemID(), "modelTreeName" ) 
self.modelTreeChildrenIDs = vtk.vtkIdList()
self.shNode.GetItemChildren(self.modelTreeID, self.modelTreeChildrenIDs)
for i in xrange(self.modelTreeChildrenIDs GetNumberOfIds()):
  print("modelChildID   : ", self.modelTreeChildrenIDs.GetId(i)) 
  print("modelChildName : ", self.shNode.GetItemName(self.modelTreeChildrenIDs.GetId(i)))
  print("modelChildNode  : ", self.shNode.GetItemDataNode(self.modelTreeChildrenIDs.GetId(i)))

Saima · November 7, 2021, 5:29am

Hi,
I tried to get model itself using

shNode.GetItemDataNode(self.modelTreeChildrenIDs.GetId(i))

but I can get the name and id but not the model node.

Is there a bug in the latest version for this 
I get 

modelChildID   :  534
modelChildName :  B_22
modelChildNode  :  None


regards,
Saima Safdar

Fernando · November 7, 2021, 11:13am

Hi @Saima. I think the method is called GetID, not GetId.

Saima · November 8, 2021, 3:34am

Hi,
No the function is GetId() with GetID() its giving error

File “/home/saima/Downloads/Downloads/slicer/BoneModelExtractor/BoneModelExtractor/BoneModelExtractor.py”, line 454, in process
print("modelChildNode : ", shNode.GetItemDataNode(modelTreeChildrenIDs.GetID(i)))
AttributeError: ‘vtkCommonCorePython.vtkIdList’ object has no attribute ‘GetID’

Saima · November 8, 2021, 3:37am

Hi,
fior me

self.shNode.GetItemDataNode(self.modelTreeChildrenIDs.GetId(i)))

is not working its giving none

any help.

Regards,
Saima Safdar

Topic		Replies	Views
The problem between Node and vtkProp3D Development	5	712	January 31, 2019
adding render actors to model node Support	2	256	May 13, 2021
Beam Visualisation in Slicer RT Support	10	665	January 12, 2020
Create a vtkMRMLModelNode from a polyData Development	2	609	October 5, 2019
Model Sequence has no display nodes Support	2	262	March 11, 2019

Create new vtkMRMLModelHierarchyNode() and add models

Related Topics