Method OnMousePressed

Modules/Segmentation/src/Interactions/mitkRegionGrowingTool.cpp:287–373

}

void mitk::RegionGrowingTool::OnMousePressed(StateMachineAction*, InteractionEvent* interactionEvent)
{
  auto* positionEvent = dynamic_cast<mitk::InteractionPositionEvent*>(interactionEvent);
  if (nullptr == positionEvent)
  {
    return;
  }

  m_LastEventSender = positionEvent->GetSender();
  m_LastEventSlice = m_LastEventSender->GetSlice();
  m_LastScreenPosition = Point2I(positionEvent->GetPointerPositionOnScreen());

  // ReferenceSlice is from the underlying image, WorkingSlice from the active segmentation (can be empty)
  m_ReferenceSlice = FeedbackContourTool::GetAffectedReferenceSlice(positionEvent);
  m_WorkingSlice = FeedbackContourTool::GetAffectedWorkingSlice(positionEvent);

  if (m_WorkingSlice.IsNull())
  {
    // can't do anything without a working slice (i.e. a possibly empty segmentation)
    return;
  }

  // Determine if the user clicked inside or outside of the working slice (i.e. the whole volume)
  mitk::BaseGeometry::Pointer workingSliceGeometry;
  workingSliceGeometry = m_WorkingSlice->GetGeometry();
  workingSliceGeometry->WorldToIndex(positionEvent->GetPositionInWorld(), m_SeedPoint);
  itk::Index<2> indexInWorkingSlice2D;
  indexInWorkingSlice2D[0] = m_SeedPoint[0];
  indexInWorkingSlice2D[1] = m_SeedPoint[1];

  if (!workingSliceGeometry->IsIndexInside(m_SeedPoint))
  {
    MITK_DEBUG << "OnMousePressed: point " << positionEvent->GetPositionInWorld() << " (index coordinates "
               << m_SeedPoint << ") is not inside working slice";
    return;
  }

  mitk::BaseGeometry::Pointer referenceSliceGeometry;
  referenceSliceGeometry = m_ReferenceSlice->GetGeometry();
  itk::Index<3> indexInReferenceSlice;
  itk::Index<2> indexInReferenceSlice2D;
  referenceSliceGeometry->WorldToIndex(positionEvent->GetPositionInWorld(), indexInReferenceSlice);
  indexInReferenceSlice2D[0] = indexInReferenceSlice[0];
  indexInReferenceSlice2D[1] = indexInReferenceSlice[1];

  // Get seed neighborhood
  ScalarType averageValue(0);
  AccessFixedDimensionByItk_3(m_ReferenceSlice, GetNeighborhoodAverage, 2, indexInReferenceSlice2D, &averageValue, 1);
  m_SeedValue = averageValue;
  MITK_DEBUG << "Seed value is " << m_SeedValue;

  // Calculate initial thresholds
  AccessFixedDimensionByItk(m_ReferenceSlice, CalculateInitialThresholds, 2);
  m_Thresholds[0] = m_InitialThresholds[0];
  m_Thresholds[1] = m_InitialThresholds[1];

  // Perform region growing
  mitk::Image::Pointer resultImage = mitk::Image::New();