Difference between revisions of "Documentation/Nightly/Developers/SlicerExecutionModel/Programmatic Invocation"

From Slicer Wiki
Jump to: navigation, search
Line 109: Line 109:
 
== Module entry points ==
 
== Module entry points ==
  
The Tcl example above assumes that the entry point of [http://www.na-mic.org/Slicer/Documentation/Slicer3/html/classModuleDescription.html module description] has been initialized with the correct target. For the shared object modules, description target specifies module entry point as a string. In case if such target has not been initialized properly (i.e., is set to "Unknown"), it has to be initialized manually. Otherwise, module execution will not happen (this can be concluded from the warning ''Module reports that it is a Shared Object Module but does not have a shared object module target'', and the command line will begin with "Unknown" instead of "slicer:x.....").  
+
The Tcl example above assumes that the entry point of [http://www.na-mic.org/Slicer/Documentation/Slicer3/html/classModuleDescription.html module description] has been initialized with the correct target. For the shared object modules, description target specifies module entry point as a string. In case if such target has not been initialized properly (i.e., is set to "Unknown"), it has to be initialized manually. Otherwise, module execution will not happen (this can be concluded from the warning ''Module reports that it is a Shared Object Module but does not have a shared object module target'', and the command line will begin with "Unknown" instead of "slicer:0x.....").  
  
 
Here's how the module entry point can be properly initialized (refer to ProcessGUIEvents() in ''Modules/CommandLineModule/vtkCommandLineModuleGUI.cxx'', Slicer3 trunk). NOTE: this has been tested for the C++ "Linear registration" module.
 
Here's how the module entry point can be properly initialized (refer to ProcessGUIEvents() in ''Modules/CommandLineModule/vtkCommandLineModuleGUI.cxx'', Slicer3 trunk). NOTE: this has been tested for the C++ "Linear registration" module.

Revision as of 21:09, 14 November 2008

Home < 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.


  #
  # 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


Module entry points

The Tcl example above assumes that the entry point of module description has been initialized with the correct target. For the shared object modules, description target specifies module entry point as a string. In case if such target has not been initialized properly (i.e., is set to "Unknown"), it has to be initialized manually. Otherwise, module execution will not happen (this can be concluded from the warning Module reports that it is a Shared Object Module but does not have a shared object module target, and the command line will begin with "Unknown" instead of "slicer:0x.....").

Here's how the module entry point can be properly initialized (refer to ProcessGUIEvents() in Modules/CommandLineModule/vtkCommandLineModuleGUI.cxx, Slicer3 trunk). NOTE: this has been tested for the C++ "Linear registration" module.

  ModuleDescription moduleDesc = moduleNode->GetModuleDescription();

  // Set up the temporary directory for the module (not sure why this is not done in Tcl example)
  moduleLogic->SetTemporaryDirectory(app->GetTemporaryDirectory());

  if(moduleDesc.GetTarget() == "Unknown"){

    // Entry point is unknown
    assert(moduleDesc.GetType() == "SharedObjectModule");
    typedef int (*ModuleEntryPoint)(int argc, char* argv[]);
    itksys::DynamicLoader::LibraryHandle lib =
      itksys::DynamicLoader::OpenLibrary(moduleDesc.GetLocation().c_str());

    if(lib){
      ModuleEntryPoint entryPoint =
        (ModuleEntryPoint) itksys::DynamicLoader::GetSymbolAddress(
          lib, "ModuleEntryPoint");
      if(entryPoint){
        char entryPointAsText[256];
        std::string entryPointAsString;
        sprintf(entryPointAsText, "%p", entryPoint);
        entryPointAsString = std::string("slicer:")+entryPointAsText;

        // Set target in the description
        moduleDesc.SetTarget(entryPointAsString);

        // Important: Update the description in the module node -- description is not a pointer!
        moduleNode->SetModuleDescription(moduleDesc);
      } else {
        std::cerr << "Failed to find entry point. Abort." << std::endl;
        abort();
      }
    } else {
      std::cerr << "Failed to locate module library. Abort." << std::endl;
      abort();
    }
  }

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()->SetCommandLineModuleNode(linearRegistrationNode);

  // Now check and initialize the module entry point, if necessary
  ModuleDescription moduleDesc =
    linearRegistrationNode->GetModuleDescription();
  moduleLogic->SetTemporaryDirectory(app->GetTemporaryDirectory());
  if(moduleDesc.GetTarget() == "Unknown"){
    // Entry point is unknown
    // "Linear registration" is shared object module, at least at this moment
    assert(moduleDesc.GetType() == "SharedObjectModule");
    typedef int (*ModuleEntryPoint)(int argc, char* argv[]);
    itksys::DynamicLoader::LibraryHandle lib =
      itksys::DynamicLoader::OpenLibrary(moduleDesc.GetLocation().c_str());
    if(lib){
      ModuleEntryPoint entryPoint =
        (ModuleEntryPoint) itksys::DynamicLoader::GetSymbolAddress(lib, "ModuleEntryPoint");

      if(entryPoint){
        char entryPointAsText[256];
        std::string entryPointAsString;

        sprintf(entryPointAsText, "%p", entryPoint);
        entryPointAsString = std::string("slicer:")+entryPointAsText;
        moduleDesc.SetTarget(entryPointAsString);
        linearRegistrationNode->SetModuleDescription(moduleDesc);
      } else {
        std::cerr << "Failed to find entry point for Linear registration. Abort." << std::endl;
        abort();
      }
    } else {
      std::cerr << "Failed to locate module library. Abort." << std::endl;
      abort();
    }
  }

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

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

  linearRegistrationNode->Delete();

  // I am not sure next line is necessary, since we have not created the GUI. 
  // GUI is however deleted in the Tcl example above -- let the Slicer experts decide, what is right.
  // linearRegistrationGUI->Delete();
}