Method CreateSurface

Modules/AlgorithmsExt/src/mitkLabeledImageToSurfaceFilter.cpp:158–271

}

void mitk::LabeledImageToSurfaceFilter::CreateSurface(int time,
                                                      vtkImageData *vtkimage,
                                                      mitk::Surface *surface,
                                                      mitk::LabeledImageToSurfaceFilter::LabelType label)
{
  vtkImageChangeInformation *indexCoordinatesImageFilter = vtkImageChangeInformation::New();
  indexCoordinatesImageFilter->SetInputData(vtkimage);
  indexCoordinatesImageFilter->SetOutputOrigin(0.0, 0.0, 0.0);

  vtkImageThreshold *threshold = vtkImageThreshold::New();
  threshold->SetInputConnection(indexCoordinatesImageFilter->GetOutputPort());
  // indexCoordinatesImageFilter->Delete();
  threshold->SetInValue(100);
  threshold->SetOutValue(0);
  threshold->ThresholdBetween(label, label);
  threshold->SetOutputScalarTypeToUnsignedChar();
  threshold->ReleaseDataFlagOn();

  vtkImageGaussianSmooth *gaussian = vtkImageGaussianSmooth::New();
  gaussian->SetInputConnection(threshold->GetOutputPort());
  // threshold->Delete();
  gaussian->SetDimensionality(3);
  gaussian->SetRadiusFactor(0.49);
  gaussian->SetStandardDeviation(GetGaussianStandardDeviation());
  gaussian->ReleaseDataFlagOn();
  gaussian->UpdateInformation();
  gaussian->Update();

  // MarchingCube -->create Surface
  vtkMarchingCubes *skinExtractor = vtkMarchingCubes::New();
  skinExtractor->ReleaseDataFlagOn();
  skinExtractor->SetInputConnection(gaussian->GetOutputPort()); // RC++
  indexCoordinatesImageFilter->Delete();
  skinExtractor->SetValue(0, 50);

  vtkPolyData *polydata;
  skinExtractor->Update();
  polydata = skinExtractor->GetOutput();
  polydata->Register(nullptr); // RC++
  skinExtractor->Delete();

  if (m_Smooth)
  {
    vtkSmoothPolyDataFilter *smoother = vtkSmoothPolyDataFilter::New();
    // read poly1 (poly1 can be the original polygon, or the decimated polygon)
    smoother->SetInputData(polydata); // RC++
    smoother->SetNumberOfIterations(m_SmoothIteration);
    smoother->SetRelaxationFactor(m_SmoothRelaxation);
    smoother->SetFeatureAngle(60);
    smoother->FeatureEdgeSmoothingOff();
    smoother->BoundarySmoothingOff();
    smoother->SetConvergence(0);

    polydata->Delete(); // RC--
    smoother->Update();
    polydata = smoother->GetOutput();
    polydata->Register(nullptr); // RC++
    smoother->Delete();