mitk::Gizmo::HandleType mitk::GizmoInteractor::PickFrom2D(const InteractionPositionEvent* positionEvent)
{
BaseRenderer* renderer = positionEvent->GetSender();
auto mapper = GetDataNode()->GetMapper(BaseRenderer::Standard2D);
auto gizmo_mapper = dynamic_cast<GizmoMapper2D*>(mapper); auto& picker = m_Picker[renderer];
if ( picker == nullptr )
{
picker = vtkSmartPointer<vtkCellPicker>::New();
picker->SetTolerance(0.005);
if ( gizmo_mapper )
{ // doing this each time is bizarre
picker->AddPickList(gizmo_mapper->GetVtkProp(renderer));
picker->PickFromListOn();
}
}
auto displayPosition = positionEvent->GetPointerPositionOnScreen();
picker->Pick( displayPosition[0], displayPosition[1], 0,
positionEvent->GetSender()->GetVtkRenderer() );
vtkIdType pickedPointID = picker->GetPointId();
if (pickedPointID == -1)
{
return Gizmo::NoHandle;
}
vtkPolyData* polydata = gizmo_mapper->GetVtkPolyData(renderer);
if (polydata && polydata->GetPointData() && polydata->GetPointData()->GetScalars())
{
double dataValue = polydata->GetPointData()->GetScalars()->GetTuple1(pickedPointID);
return m_Gizmo->GetHandleFromPointDataValue(dataValue);
}
return Gizmo::NoHandle;
}
void mitk::CreateDistanceImageFromSurfaceFilter::PreprocessContourPoints()
{
unsigned int numberOfInputs = this->GetNumberOfIndexedInputs();
if (numberOfInputs == 0)
{
MITK_ERROR << "mitk::CreateDistanceImageFromSurfaceFilter: No input available. Please set an input!" << std::endl;
itkExceptionMacro("mitk::CreateDistanceImageFromSurfaceFilter: No input available. Please set an input!");
return;
}
// First of all we have to extract the nomals and the surface points.
// Duplicated points can be eliminated
vtkSmartPointer<vtkPolyData> polyData;
vtkSmartPointer<vtkDoubleArray> currentCellNormals;
vtkSmartPointer<vtkCellArray> existingPolys;
vtkSmartPointer<vtkPoints> existingPoints;
double p[3];
PointType currentPoint;
PointType normal;
for (unsigned int i = 0; i < numberOfInputs; i++)
{
auto currentSurface = this->GetInput(i);
polyData = currentSurface->GetVtkPolyData();
if (polyData->GetNumberOfPolys() == 0)
{
MITK_INFO << "mitk::CreateDistanceImageFromSurfaceFilter: No input-polygons available. Please be sure the input "
"surface consists of polygons!"
<< std::endl;
}
currentCellNormals = vtkDoubleArray::SafeDownCast(polyData->GetCellData()->GetNormals());
existingPolys = polyData->GetPolys();
existingPoints = polyData->GetPoints();
existingPolys->InitTraversal();
vtkIdType *cell(nullptr);
vtkIdType cellSize(0);
for (existingPolys->InitTraversal(); existingPolys->GetNextCell(cellSize, cell);)
{
for (vtkIdType j = 0; j < cellSize; j++)
{
existingPoints->GetPoint(cell[j], p);
currentPoint.copy_in(p);
int count = std::count(m_Centers.begin(), m_Centers.end(), currentPoint);
if (count == 0)
{
double currentNormal[3];
currentCellNormals->GetTuple(cell[j], currentNormal);
normal.copy_in(currentNormal);
m_Normals.push_back(normal);
m_Centers.push_back(currentPoint);
}
} // end for all points
} // end for all cells
} // end for all outputs
}
