使用vtkDICOMImageReader 读取序列CT图像

/*
使用vtkDICOMImageReader 读取序列CT图像
*/
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkVolume16Reader.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkOutlineFilter.h"
#include "vtkCamera.h"
#include "vtkProperty.h"
#include "vtkPolyDataNormals.h"
#include "vtkContourFilter.h"
#include "vtkMarchingCubes.h"
#include "vtkDICOMImageReader.h"
#include "vtkImageCast.h"
#include "vtkDecimatePro.h"
#include "vtkStripper.h"
#include "vtkImageShrink3D.h"
#include "vtkSmoothPolyDataFilter.h"
#include "vtkTriangleFilter.h"
#include "vtkFeatureEdges.h"
#include "vtkPolyDataWriter.h"
#include "vtkRecursiveDividingCubes.h"
int main ()
{
vtkRenderer *aRenderer = vtkRenderer::New();
vtkRenderWindow *renWin = vtkRenderWindow::New();
renWin->AddRenderer(aRenderer); vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New(); iren->SetRenderWindow(renWin);
vtkDICOMImageReader *reader = vtkDICOMImageReader::New();
reader->SetDataByteOrderToLittleEndian();
reader->SetDirectoryName("E://CT/head10/");
vtkImageShrink3D *shrink=vtkImageShrink3D::New();
shrink->SetInput((vtkDataObject *)reader->GetOutput());
shrink->SetShrinkFactors(4,4,1);
vtkRecursiveDividingCubes *iso=vtkRecursiveDividingCubes::New(); iso->SetInputConnection(shrink->GetOutputPort());
iso->SetValue(300);
iso->SetDistance(0.5);
iso->SetIncrement(2);
/*
vtkPolyDataWriter *wSP=vtkPolyDataWriter::New();
wSP->SetInput(iso->GetOutput());
wSP->SetFileName("E://CT/aaa2a.VTK");
wSP->Write();
wSP->Delete();
//保存为VTK格式
*/
vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New();
skinMapper->SetInput(iso->GetOutput());
skinMapper->ScalarVisibilityOff();
vtkActor *skin = vtkActor::New();
skin->SetMapper(skinMapper);
skin->GetProperty()->SetDiffuseColor(1, 0.49, 0.25);
skin->GetProperty()->SetSpecular(.3);
skin->GetProperty()->SetSpecularPower(20);
vtkCamera *aCamera = vtkCamera::New();
aCamera->SetViewUp (0, 0, -1);
aCamera->SetPosition (0, 1, 0);
aCamera->SetFocalPoint (0, 0, 0);
aCamera->ComputeViewPlaneNormal();
aRenderer->AddActor(skin);
aRenderer->SetActiveCamera(aCamera);
aRenderer->ResetCamera ();
aCamera->Dolly(1.5);
aRenderer->SetBackground(1,1,1);
renWin->SetSize(640, 480);
aRenderer->ResetCameraClippingRange ();
iren->Initialize();
iren->Start();
iso->Delete();
// skinNormals->Delete(); skinMapper->Delete(); skin->Delete(); aCamera->Delete();
iren->Delete();
renWin->Delete(); aRenderer->Delete(); reader->Delete();
return 0;
}。