Method OnReinit3DInterpolation

Modules/SegmentationUI/src/QmitkSlicesInterpolator.cpp:1206–1267

}

void QmitkSlicesInterpolator::OnReinit3DInterpolation()
{
  //  Step 1. Load from the isContourPlaneGeometry nodes the contourNodes.
  mitk::NodePredicateProperty::Pointer pred =
    mitk::NodePredicateProperty::New("isContourPlaneGeometry", mitk::BoolProperty::New(true));
  mitk::DataStorage::SetOfObjects::ConstPointer contourNodes =
    m_DataStorage->GetDerivations(m_ToolManager->GetWorkingData(0), pred);

  if (contourNodes->Size() != 0)
  {
    if (m_ToolManager->GetWorkingData(0) != nullptr)
    {
      try
      {
        auto labelSetImage = GetData<mitk::MultiLabelSegmentation>(m_ToolManager->GetWorkingData(0));
        if (!labelSetImage->GetTimeGeometry()->IsValidTimePoint(m_TimePoint))
        {
          MITK_ERROR << "Invalid time point requested for interpolation pipeline.";
          return;
        }

        mitk::SurfaceInterpolationController::CPIVector newCPIs;
        //  Adding label and timeStep information for the contourNodes.
        for (auto it = contourNodes->Begin(); it != contourNodes->End(); ++it)
        {
          auto contourNode = it->Value();
          auto labelID = dynamic_cast<mitk::UShortProperty *>(contourNode->GetProperty("labelID"))->GetValue();
          auto timeStep = dynamic_cast<mitk::IntProperty *>(contourNode->GetProperty("timeStep"))->GetValue();

          auto planeGeometry = dynamic_cast<mitk::PlanarFigure *>(contourNode->GetData())->GetPlaneGeometry();
          auto groupID = labelSetImage->GetGroupIndexOfLabel(labelID);
          auto sliceImage = ExtractSliceFromImage(labelSetImage->GetGroupImage(groupID), planeGeometry, timeStep);
          mitk::ImageToContourFilter::Pointer contourExtractor = mitk::ImageToContourFilter::New();
          contourExtractor->SetInput(sliceImage);
          contourExtractor->SetContourValue(labelID);
          contourExtractor->Update();
          mitk::Surface::Pointer contour = contourExtractor->GetOutput();

          if (contour->GetVtkPolyData()->GetNumberOfPoints() == 0)
            continue;

          contour->DisconnectPipeline();
          newCPIs.emplace_back(contour, planeGeometry->Clone(),labelID,timeStep);
        }
        m_SurfaceInterpolator->CompleteReinitialization(newCPIs);
      }
      catch(const std::exception& e)
      {
        MITK_ERROR << "Exception thrown casting toolmanager working data to labelsetImage";
      }
    }
  }
  else
  {
    m_BtnApply3D->setEnabled(false);
    QMessageBox errorInfo;
    errorInfo.setWindowTitle("Reinitialize surface interpolation");
    errorInfo.setIcon(QMessageBox::Information);