使用vtkDICOMImageReader 读取序列CT图像.docxVIP

/*使用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.hint 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();aRend