VTK系列57_VTK对几何体网格细化（多分辨率处理）

实例57：网格细化（多分辨率处理）

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkRenderingFreeType);
VTK_MODULE_INIT(vtkInteractionStyle);

#include <vtkSmartPointer.h>
#include <vtkPolyDataReader.h>
#include <vtkPolyData.h>
#include <vtkLinearSubdivisionFilter.h>
#include <vtkLoopSubdivisionFilter.h>
#include <vtkButterflySubdivisionFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkCamera.h>
#include <vtkRenderWindowInteractor.h>

int main()
{
	//读数据
	vtkSmartPointer<vtkPolyDataReader> reader =
		vtkSmartPointer<vtkPolyDataReader>::New();
	reader->SetFileName("C:\\Users\\Administrator\\Desktop\\VTK2\\hellovtk\\vtk_图像处理学习\\第六章_图像处理\\data\\fran_cut.vtk");
	reader->Update();

	vtkSmartPointer<vtkPolyData> orig = reader->GetOutput();
	std::cout << "original" << "-----------------------" << std::endl;
	std::cout << "模型点数为： " << orig->GetNumberOfPoints() << std::endl;
	std::cout << "模型面数为： " << orig->GetNumberOfPolys() << std::endl;
	//线性网格细分滤波器
	vtkSmartPointer<vtkLinearSubdivisionFilter> linear =
		vtkSmartPointer<vtkLinearSubdivisionFilter>::New();
	linear->SetInputData(orig);
	linear->SetNumberOfSubdivisions(4);
	linear->Update();

	vtkSmartPointer<vtkPolyData> linearInfo = linear->GetOutput();
	std::cout << "linear" << "-----------------------" << std::endl;
	std::cout << "模型点数为： " << linearInfo->GetNumberOfPoints() << std::endl;
	std::cout << "模型面数为： " << linearInfo->GetNumberOfPolys() << std::endl;

	//Loop网格细分滤波器
	vtkSmartPointer<vtkLoopSubdivisionFilter> loop =
		vtkSmartPointer<vtkLoopSubdivisionFilter>::New();
	loop->SetInputData(orig);
	loop->SetNumberOfSubdivisions(4);
	loop->Update();

	vtkSmartPointer<vtkPolyData> loopInfo = loop->GetOutput();
	std::cout << "loop" << "-----------------------" << std::endl;
	std::cout << "模型点数为： " << loopInfo->GetNumberOfPoints() << std::endl;
	std::cout << "模型面数为： " << loopInfo->GetNumberOfPolys() << std::endl;

	//butterfly网格细分滤波器
	vtkSmartPointer<vtkButterflySubdivisionFilter> butterfly =
		vtkSmartPointer<vtkButterflySubdivisionFilter>::New();
	butterfly->SetInputData(orig);
	butterfly->SetNumberOfSubdivisions(4);
	butterfly->Update();

	vtkSmartPointer<vtkPolyData> butterflyInfo = butterfly->GetOutput();
	std::cout << "butterfly" << "-----------------------" << std::endl;
	std::cout << "模型点数为： " << butterflyInfo->GetNumberOfPoints() << std::endl;
	std::cout << "模型面数为： " << butterflyInfo->GetNumberOfPolys() << std::endl;
	
	vtkSmartPointer<vtkPolyDataMapper> origMapper =
		vtkSmartPointer<vtkPolyDataMapper>::New();
	origMapper->SetInputData(orig);
	vtkSmartPointer<vtkActor> origActor =
		vtkSmartPointer<vtkActor>::New();
	origActor->SetMapper(origMapper);

	vtkSmartPointer<vtkPolyDataMapper> linearMapper =
		vtkSmartPointer<vtkPolyDataMapper>::New();
	linearMapper->SetInputData(linear->GetOutput());
	vtkSmartPointer<vtkActor> linearActor =
		vtkSmartPointer<vtkActor>::New();
	linearActor->SetMapper(linearMapper);

	vtkSmartPointer<vtkPolyDataMapper> loopMapper =
		vtkSmartPointer<vtkPolyDataMapper>::New();
	loopMapper->SetInputData(loop->GetOutput());
	vtkSmartPointer<vtkActor> loopActor =
		vtkSmartPointer<vtkActor>::New();
	loopActor->SetMapper(loopMapper);

	vtkSmartPointer<vtkPolyDataMapper> butterflyMapper =
		vtkSmartPointer<vtkPolyDataMapper>::New();
	butterflyMapper->SetInputData(butterfly->GetOutput());
	vtkSmartPointer<vtkActor> butterflyActor =
		vtkSmartPointer<vtkActor>::New();
	butterflyActor->SetMapper(butterflyMapper);
	//
	double ltView[4] = { 0, 0, 0.5, 0.5 };
	double rtView[4] = { 0.5, 0, 1, 0.5 };
	double lbView[4] = { 0, 0.5, 0.5, 1 };
	double rbView[4] = { 0.5, 0.5, 1, 1 };

	vtkSmartPointer<vtkRenderer> origRender =
		vtkSmartPointer<vtkRenderer>::New();
	origRender->SetViewport(ltView);
	origRender->AddActor(origActor);
	origRender->SetBackground(1, 0, 0);

	vtkSmartPointer<vtkRenderer> linearRender =
		vtkSmartPointer<vtkRenderer>::New();
	linearRender->SetViewport(rtView);
	linearRender->AddActor(linearActor);
	linearRender->SetBackground(0, 1, 0);

	vtkSmartPointer<vtkRenderer> loopRender =
		vtkSmartPointer<vtkRenderer>::New();
	loopRender->SetViewport(lbView);
	loopRender->AddActor(loopActor);
	loopRender->SetBackground(0, 0, 1);

	vtkSmartPointer<vtkRenderer> butterflyRender =
		vtkSmartPointer<vtkRenderer>::New();
	butterflyRender->SetViewport(rbView);
	butterflyRender->AddActor(butterflyActor);
	butterflyRender->SetBackground(0, 0, 0);
	//
	vtkSmartPointer<vtkRenderWindow> rw =
		vtkSmartPointer<vtkRenderWindow>::New();
	rw->AddRenderer(origRender);
	rw->AddRenderer(linearRender);
	rw->AddRenderer(loopRender);
	rw->AddRenderer(butterflyRender);
	rw->SetSize(640, 640);
	rw->SetWindowName("PolyData Subdivision");

	origRender->GetActiveCamera()->SetPosition(0, -1, 0);
	origRender->GetActiveCamera()->SetFocalPoint(0, 0, 0);
	origRender->GetActiveCamera()->SetViewUp(0, 0, 1);
	origRender->GetActiveCamera()->Azimuth(30);
	origRender->GetActiveCamera()->Elevation(30);
	origRender->ResetCamera();//刷新照相机  
	linearRender->SetActiveCamera(origRender->GetActiveCamera());
	loopRender->SetActiveCamera(origRender->GetActiveCamera());
	butterflyRender->SetActiveCamera(origRender->GetActiveCamera());
	
	vtkSmartPointer<vtkRenderWindowInteractor> rwi =
		vtkSmartPointer<vtkRenderWindowInteractor>::New();
	rwi->SetRenderWindow(rw);
	rwi->Start();
	rw->Render();

	return 0;
}

 

        放大细节 

红色视窗为原始三角网格模型，绿色视窗采用了线性细分算法；蓝色视窗采用Loop细分算法；黑色视窗采用了Butterfly细分算法。可见Loop细分、Butterfly细分能够得到较为光滑的效果；效果优于线性细分算法。

VTK系列目录：

1 VTK基本概念

2 VTK图像处理

3 VTK图形处理

4 VTK体绘制