Documentation/Nightly/Developers/SlicerExecutionModel/Programmatic Invocation

Goal

A large amount of functionality is available in command line modules for processing images and models. Some applications would like to use several modules together in a workflow or as a background processing task.

Implementation Example

This example shows calling the Model Maker (from the Editor). This example is in Tcl, but the same API can be used from C++ or Python. The idea is to query for the Model Maker module (and bail out if it is not loaded). Then to create a MRMLNode with the parameters you want for the model maker. These would be the same parameters you would set interactively. The model maker logic and gui classes are then told about the node and told to start processing. The processing then goes on in the background.

Note: check the command line module xml file for the to know the correct symbols and data types for parameters (see the <name> field for the correct parameter name string).

  #
  # create a model using the command line module
  # based on the current editor parameters
  #

  #
  # get the image data for the label layer
  #
  set sliceLogic [lindex [vtkSlicerSliceLogic ListInstances] 0]
  set layerLogic [$sliceLogic GetLabelLayer]
  set volumeNode [$layerLogic GetVolumeNode]
  if { $volumeNode == "" } {
    errorDialog "Cannot make model - no volume node for $layerLogic in $sliceLogic."
    return
  }

  #
  # find the Model Maker
  # - call Enter to be sure GUI has been built
  # - find the GUI - for non-tcl use, try vtkSlicerApplication::GetModuleGUIByName
  set modelMaker ""
  foreach gui [vtkCommandLineModuleGUI ListInstances] {
    if { [$gui GetGUIName] == "Model Maker" } {
      set modelMaker $gui
    }
  }

  if { $modelMaker == "" } {
    errorDialog "Cannot make model: no Model Maker Module found."
  }

  $modelMaker Enter

  #
  # set up the model maker node
  #
  set moduleNode [$::slicer3::MRMLScene CreateNodeByClass "vtkMRMLCommandLineModuleNode"]
  $::slicer3::MRMLScene AddNode $moduleNode
  $moduleNode SetModuleDescription "Model Maker"

  $moduleNode SetParameterAsString "Name" "Quick Model"
  $moduleNode SetParameterAsString "FilterType" "Sinc"
  $moduleNode SetParameterAsBool "GenerateAll" "0"
  $moduleNode SetParameterAsString "Labels" [EditorGetPaintLabel]
  $moduleNode SetParameterAsBool "JointSmooth" 1
  $moduleNode SetParameterAsBool "SplitNormals" 1
  $moduleNode SetParameterAsBool "PointNormals" 1
  $moduleNode SetParameterAsBool "SkipUnNamed" 1
  $moduleNode SetParameterAsInt "Start" -1
  $moduleNode SetParameterAsInt "End" -1
  if { [[$o(smooth) GetWidget] GetSelectedState] } {
    $moduleNode SetParameterAsDouble "Decimate" 0.25
    $moduleNode SetParameterAsDouble "Smooth" 10
  } else {
    $moduleNode SetParameterAsDouble "Decimate" 0
    $moduleNode SetParameterAsDouble "Smooth" 0
  }

  $moduleNode SetParameterAsString "InputVolume" [$volumeNode GetID]

  #
  # output 
  # - make a new hierarchy node if needed
  #
  set outHierarchy [[$::slicer3::MRMLScene GetNodesByClassByName "vtkMRMLModelHierarchyNode" "Editor Models"] GetItemAsObject 0]
  if { $outHierarchy == "" } {
    set outHierarchy [$::slicer3::MRMLScene CreateNodeByClass "vtkMRMLModelHierarchyNode"]
    $outHierarchy SetScene $::slicer3::MRMLScene
    $outHierarchy SetName "Editor Models"
    $::slicer3::MRMLScene AddNode $outHierarchy
  }

  $moduleNode SetParameterAsString "ModelSceneFile" [$outHierarchy GetID]



  # 
  # run the task (in the background)
  # - use the GUI to provide progress feedback
  # - use the GUI's Logic to invoke the task
  # - model will show up when the processing is finished
  #
  $modelMaker SetCommandLineModuleNode $moduleNode
  [$modelMaker GetLogic] SetCommandLineModuleNode $moduleNode
  [$modelMaker GetLogic] Apply $moduleNode

  $this statusText "Model Making Started..."

  #
  # clean up our references
  #
  $moduleNode Delete
  $outHierarchy Delete
  $modelMaker Enter

When you are calling a shared object module, you need to add another call to initialize the pointer to the module entry point:

  ...
  $modelMaker SetCommandLineModuleNode $moduleNode
  [$modelMaker GetLogic] SetCommandLineModuleNode $moduleNode
  [$modelMaker GetLogic] LazyEvaluateModuleTarget $moduleNode
  [$modelMaker GetLogic] Apply $moduleNode
  ...

C++ example

Update the CMakeLists.txt to include

${CommandLineModule_SOURCE_DIR}
${CommandLineModule_BINARY_DIR}

The following code shows an example of executing Linear registration CLM from the C++ code. This was done in the context of ChangeTracker Slicer GUI module.

// Use LinearRegistration module to rigidly align input scans
void DoITKGlobalRegistration(vtkSlicerApplication *app){
  // Init some useful references
  vtkCommandLineModuleGUI *linearRegistrationGUI = NULL;
  vtkMRMLCommandLineModuleNode *linearRegistrationNode = NULL;
  vtkMRMLScene *scene = this->ChangeTrackerNode->GetScene();
  vtkMRMLChangeTrackerNode *ctNode = this->ChangeTrackerNode;
  vtkMRMLScalarVolumeNode *outputVolumeNode = NULL;
  vtkMRMLTransformNode *lrTransform = NULL;

  linearRegistrationGUI = static_cast<vtkCommandLineModuleGUI*>(app->GetModuleGUIByName("Linear registration"));
  if(!linearRegistrationGUI){
    std::cerr << "Cannot find LinearRegistration module. Aborting." << std::endl;
    abort();
  }

  linearRegistrationGUI->Enter();

  linearRegistrationNode = static_cast<vtkMRMLCommandLineModuleNode*>(scene->CreateNodeByClass("vtkMRMLCommandLineModuleNode"));
  if(!linearRegistrationNode){
    std::cerr << "Cannot create LinearRegistration node. Aborting." << std::endl;
    abort;
  }

  // Add node to the scene
  scene->AddNode(linearRegistrationNode);
  // This will actually search and find the initialized description for the node
  linearRegistrationNode->SetModuleDescription("Linear registration");

  // Create output volume node
  vtkMRMLScalarVolumeNode *scan1Node = static_cast<vtkMRMLScalarVolumeNode*>(scene->GetNodeByID(ctNode->GetScan1_Ref()));
  outputVolumeNode = this->CreateVolumeNode(scan1Node, "TG_scan2_Global");
  ctNode->SetScan2_GlobalRef(outputVolumeNode->GetID());

  linearRegistrationNode->SetParameterAsString("FixedImageFileName",
    ctNode->GetScan1_Ref());

  linearRegistrationNode->SetParameterAsString("MovingImageFileName",
    ctNode->GetScan2_Ref());

  linearRegistrationNode->SetParameterAsString("ResampledImageFileName",
    ctNode->GetScan2_GlobalRef());

  vtkCommandLineModuleLogic *moduleLogic = linearRegistrationGUI->GetLogic();

  linearRegistrationGUI->SetCommandLineModuleNode(linearRegistrationNode);
  linearRegistrationGUI->GetLogic()->LazyEvaluateModuleTarget(linearRegistrationNode);
  linearRegistrationGUI->GetLogic()->SetCommandLineModuleNode(linearRegistrationNode);
  linearRegistrationGUI->GetLogic()->Apply(linearRegistrationNode);

  std::cout << "Linear registration module started..." << std::endl;

  linearRegistrationNode->Delete();
}

Python Example

The following example calls the Substract Images command line module directly from Python (presented by Luca Antiga at the AHM 2009).

import Slicer
from Slicer import slicer
volume1 = slicer.MRMLScene.GetNodeByID(“vtkMRMLScalarVolumeNode1”)
volume2 = slicer.MRMLScene.GetNodeByID(“vtkMRMLScalarVolumeNode2”)
plugin = Slicer.Plugin(“Subtract Images”)
plugin.Execute(volume1,volume2)

The next snippet calls the Grayscale Model Maker on a provided vtkMRMLScalarVolumeNode named outVolume and saves the generated gray scale model to a provided vtkMRMLModelNode named outGeometry.

import Slicer
from Slicer import slicer
# get or create the outVolume and outGeometry nodes
plugin = Slicer.Plugin("Grayscale Model Maker")
plugin.Execute(outVolume,outGeometry,threshold=10,smooth=15,decimate=0.25,splitnormals=True,pointnormals=True)

In both cases the called command line module will run hidden.