task2.cpp

本程序完成任务并行等值面提取功能

#include "programme.h"
#include "vtkVolume16Reader.h"
#include "vtkContourFilter.h"
#include "vtkImageData.h"
#include "vtkOutputPort.h"
#include "vtkPolyData.h"
#include "vtkPolyDataNormals.h"
#include "vtkStripper.h"

void task2()
{  
  vtkVolume16Reader *v16 = vtkVolume16Reader::New();
    v16->SetDataDimensions(64,64);
    v16->SetDataByteOrderToLittleEndian();
	v16->SetFilePrefix ("D:/headsq/quarter");
    v16->SetImageRange(1, 93);
    v16->SetDataSpacing (3.2, 3.2, 1.5);

	 vtkContourFilter* contour = vtkContourFilter::New();
  contour->SetInput(v16->GetOutput());
  contour->SetNumberOfContours(1);
  contour->SetValue(0, 1700);
  vtkPolyDataNormals *boneNormals = vtkPolyDataNormals::New();
       boneNormals->SetInput(contour->GetOutput());
       boneNormals->SetFeatureAngle(60.0);
       vtkStripper *boneStripper = vtkStripper::New();
       boneStripper->SetInput(boneNormals->GetOutput());



  // Output port
  vtkOutputPort* op = vtkOutputPort::New();
  op->SetInput(contour->GetOutput());
  op->SetTag(11);

  // Process requests
  op->WaitForUpdate();

  // Cleanup
  op->Delete();
  v16->Delete();
  contour->Delete();
  boneNormals->Delete();
  boneStripper->Delete();

}