VTK使用入门(四)

　　#include "vtkDICOMImageReader.h"
　　#include "vtkRenderer.h"
　　#include "vtkRenderWindow.h"
　　#include "vtkRenderWindowInteractor.h"
　　//#include "vtkVolume16Reader.h"
　　#include "vtkPolyDataMapper.h"
　　#include "vtkActor.h"
　　#include "vtkProperty.h"
　　#include "vtkOutlineFilter.h"
　　#include "vtkCamera.h"
　　#include "vtkPolyDataNormals.h"
　　#include "vtkStripper.h"
　　#include "vtkMarchingCubes.h" int main(int argc, char **argv)
　　{
　　if (argc SetDirectoryName("d:/headsq/");
　　v16->SetFilePrefix(argv[1]);
　　v16->SetDataSpacing(1.0,1.0,1.0);
　　v16->Update();
　　//=====================重建皮肤=============================
　　vtkMarchingCubes *skinExtractor = vtkMarchingCubes::New();
　　skinExtractor->SetInputConnection(v16->GetOutputPort ());
　　skinExtractor->SetValue(0,500);
　　vtkPolyDataNormals *skinNormals = vtkPolyDataNormals::New();
　　skinNormals->SetInputConnection(skinExtractor->GetOu tputPort());
　　skinNormals->SetFeatureAngle(60.0);
　　vtkStripper *skinStripper = vtkStripper::New();
　　skinStripper->SetInputConnection(skinNormals->GetOut putPort());
　　vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New();
　　skinMapper->SetInputConnection(skinNormals->GetOutpu tPort());
　　skinMapper->ScalarVisibilityOff();
　　vtkActor *skin = vtkActor::New();
　　skin->SetMapper(skinMapper);
　　skin->GetProperty()->SetDiffuseColor(1, .49, .25);
　　skin->GetProperty()->SetSpecular( .3);
　　skin->GetProperty()->SetSpecularPower(20);
　　skin->GetProperty()->SetOpacity(1.0);
　　//=======================重建骨骼================================
　　vtkMarchingCubes *boneExtraractor = vtkMarchingCubes::New();
　　boneExtraractor->SetInputConnection(v16->GetOutputPo rt());
　　boneExtraractor->SetValue(0,1150);
　　vtkPolyDataNormals *boneNormals = vtkPolyDataNormals::New();
　　boneNormals->SetInputConnection(boneExtraractor->Get OutputPort());
　　boneNormals->SetFeatureAngle(60.0);
　　vtkStripper *boneStripper = vtkStripper::New();
　　boneStripper->SetInputConnection(boneStripper->GetOu tputPort());
　　vtkPolyDataMapper *boneMapper = vtkPolyDataMapper::New();
　　boneMapper->SetInputConnection(boneMapper->GetOutput Port());
　　boneMapper->ScalarVisibilityOff();
　　vtkActor *bone = vtkActor::New();
　　bone->SetMapper(boneMapper);
　　bone->GetProperty()->SetDiffuseColor(1,1, .9412);
　　//=============================建立轮廓============================
　　vtkOutlineFilter *outlineData = vtkOutlineFilter::New();
　　outlineData->SetInputConnection(v16->GetOutputPort() );
　　vtkPolyDataMapper *mapOutline = vtkPolyDataMapper::New();
　　mapOutline->SetInputConnection(outlineData->GetOutpu tPort());
　　vtkActor *outline = vtkActor::New();
　　outline->SetMapper(mapOutline);
　　outline->GetProperty()->SetColor(0,0,0);
　　//========================建立相机===============================
　　vtkCamera *aCamera = vtkCamera::New();
　　aCamera->SetViewUp(0,0,-1);
　　aCamera->SetPosition(0,1,0);
　　aCamera->SetFocalPoint(0,0,0);
　　aCamera->ComputeViewPlaneNormal();
　　//=======================准备演示=================================
　　vtkRenderer *aRenderer = vtkRenderer::New();
　　vtkRenderWindow *renWin = vtkRenderWindow::New();
　　renWin->AddRenderer(aRenderer);
　　vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
　　iren->SetRenderWindow(renWin);
　　aRenderer->AddActor(outline);
　　aRenderer->AddActor(skin);
　　aRenderer->AddActor(bone);
　　aRenderer->SetActiveCamera(aCamera);
　　aRenderer->ResetCamera();
　　aCamera->Dolly(1.5);
　　aRenderer->SetBackground(1,1,1);
　　aRenderer->ResetCameraClippingRange();
　　renWin->SetSize(640,480);
　　//=========================开始演示==============================
　　iren->Initialize();
　　iren->Start();
　　//==========================演示结束=============================
　　v16->Delete();
　　skinExtractor->Delete();
　　skinNormals->Delete();
　　skinStripper->Delete();
　　skinMapper->Delete();
　　skin->Delete();
　　boneExtraractor->Delete();
　　boneNormals->Delete();
　　boneStripper->Delete();
　　boneMapper->Delete();
　　bone->Delete();
　　outlineData->Delete();
　　mapOutline->Delete();
　　outline->Delete();
　　aCamera->Delete();
　　aRenderer->Delete();
　　renWin->Delete();
　　iren->Delete();
　　return 0;
　　}
　　＝＝＝＝＝＝＝＝＝＝＝＝＝程序来源于《Visual C++数字图象处理技术详解》＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
　　以上是在VS2005环境下VTK对医学图像成功3V重建的代码，基于MC算法，当然，前提是你配置好了VTK。
　　
　　
　　【不知道等上传后会不会出现代码丢失。。。。】
　　1，头文件的用法，声明，需要进一步搞明白；
　　2，程序的模块化后，一步一步，依次，读取-》重建皮肤-》重建骨骼-》建立轮廓-》建立相机-》准备演示-》开始演示-》演示结束~~~ 其实大致就包括信息读取，处理，新建窗口，窗口显示；
　　3，需要搞明白其中的算法；
　　4，还有很多~~~~~~~~~
　　这个只是VTK初探，慢慢学吧，呵呵。