vtktubefilter，根据场数据创建半径不同的圆柱体

最新推荐文章于 2025-05-02 14:00:00 发布

原创最新推荐文章于 2025-05-02 14:00:00 发布 · 734 阅读

4 ·

CC 4.0 BY-SA版权

文章标签：

#开发语言

本文介绍了一个使用VTK（Visualization Toolkit）创建管状几何的例子。通过管状滤波器（vtkTubeFilter），将简单的线条转换为立体的管状结构，并展示了如何根据自定义数组调整管的半径。此教程还涵盖了从设置源数据到渲染管线配置的完整过程。

网上关于vtk的东西好少，自己写了点

先看vtkuserguide：

Lines are often difficult to see and create useful images from. In this example we wrap the lines with

a tube filter. The tube filter is configured to vary the radius of the tube inversely proportional to the

velocity magnitude (i.e., a flux-preserving relationship if the flow field is incompressible). The Set

VaryRadiusToVaryRadiusByVector() enables this. You can also vary the radius by scalar value (Set

VaryRadiusToVaryRadiusByScalar()) or turn off the variable radius

(SetVaryRadiusToVaryRadiusOff()). Note that the tube filter has to be told which array to use when

scaling its radius. In this case, the array with the name “vectors” was selected using SetInputArrayToProcess().

可以使用标量或者数组的形式设定半径，因为标量影响太大，所以我这里选择使用数组

（SetInputArrayToProcess()使用

）


int main() {

	vtkSmartPointer<vtkLineSource> lineSource = vtkSmartPointer<vtkLineSource>::New();
	vtkSmartPointer<vtkPoints> points =
		vtkSmartPointer<vtkPoints>::New();
	points->InsertPoint(0, 1, 0, 0);
	points->InsertPoint(1, 2, 0, 0);
	points->InsertPoint(2, 3, 0, 0);
	lineSource->SetPoints(points);
	lineSource->Update();
	vtkSmartPointer<vtkDoubleArray> tubeRadius = vtkSmartPointer<vtkDoubleArray>::New();
	unsigned int n = lineSource->GetOutput()->GetNumberOfPoints();
	tubeRadius->SetNumberOfTuples(3);
	tubeRadius->SetName("TubeRadius");
	tubeRadius->SetTuple1(0, 0.5);
	tubeRadius->SetTuple1(1, 1.0);
	tubeRadius->SetTuple1(2, 1.0);
	lineSource->GetOutput()->GetPointData()->AddArray(tubeRadius);

	vtkSmartPointer<vtkTubeFilter> tubeFilter =
		vtkSmartPointer<vtkTubeFilter>::New();
	tubeFilter->SetInputData(lineSource->GetOutput());
	tubeFilter->SetNumberOfSides(20);
	tubeFilter->SetVaryRadiusToVaryRadiusByVector();
	tubeFilter->SetInputArrayToProcess(1, 0, 0, vtkDataObject::FIELD_ASSOCIATION_POINTS, "TubeRadius");
	tubeFilter->CappingOn();
	tubeFilter->Update();
	
	vtkSmartPointer<vtkPolyDataWriter> writer = vtkSmartPointer<vtkPolyDataWriter>::New();
	writer->SetInputData(tubeFilter->GetOutput());
	writer->SetFileName("g:/AppendCylinder.vtk");
	writer->Update();

	vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
	mapper->SetInputData(tubeFilter->GetOutput());

	vtkSmartPointer<vtkActor> Actor = vtkSmartPointer<vtkActor>::New();
	Actor->SetMapper(mapper);
	mapper->Update();

	vtkSmartPointer<vtkRenderer> renderer =
        vtkSmartPointer<vtkRenderer>::New();
    vtkSmartPointer<vtkRenderWindow> renderWindow =
        vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);
    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    renderWindowInteractor->SetRenderWindow(renderWindow);
    renderer->AddActor(Actor);
    renderer->SetBackground(.2, .3, .4);
    renderWindow->Render();
    renderWindowInteractor->Start();

    return EXIT_SUCCESS;
}