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;
}