Method RequestData

Filters/General/vtkPointsMatchingTransformFilter.cxx:95–136

------------------------------------------------------------------------------

//------------------------------------------------------------------------------
int vtkPointsMatchingTransformFilter::RequestData(vtkInformation* vtkNotUsed(request),
  vtkInformationVector** inputVector, vtkInformationVector* outputVector)
{
  vtkSmartPointer<vtkPointSet> input = vtkPointSet::GetData(inputVector[0]);
  vtkSmartPointer<vtkPointSet> output = vtkPointSet::GetData(outputVector);
  if (!input)
  {
    vtkErrorMacro(<< "Invalid or missing input");
    return 0;
  }

  // Build transform matrix
  vtkNew<vtkMatrix4x4> srcInv;
  srcInv->DeepCopy(this->SourceMatrix);
  if (std::abs(srcInv->Determinant()) < 10e-4)
  {
    vtkWarningMacro("Source matrix is not invertible. Source points are likely coplanar.");
    output->ShallowCopy(input);
    return 1;
  }
  srcInv->Invert();
  vtkNew<vtkMatrix4x4> transformMatrix;
  vtkMatrix4x4::Multiply4x4(this->TargetMatrix, srcInv, transformMatrix);

  if (this->RigidTransform)
  {
    vtkSmartPointer<vtkMatrix3x3> rotation = ::ExtractRotationFromMatrix4x4(transformMatrix);
    rotation = ::PolarDecomposition(rotation);
    ::SetRotationInMatrix4x4(rotation, transformMatrix);
  }

  // Apply the transform
  vtkNew<vtkTransform> transform;
  transform->SetMatrix(transformMatrix);
  vtkNew<vtkTransformFilter> transformFilter;
  transformFilter->SetInputData(input);
  transformFilter->SetTransform(transform);
  transformFilter->Update();

  output->ShallowCopy(transformFilter->GetOutput());
  return 1;
}

//------------------------------------------------------------------------------
void vtkPointsMatchingTransformFilter::SetSourcePoint1(double x, double y, double z)