VTK使用矢量数据弯曲几何体-白红宇

VTK使用矢量数据弯曲几何体

阅读量：5292 次

发布时间：2019-06-14

本文共 5624 字，大约阅读时间需要 18 分钟。

　　 is a filter that modifies point coordinates by moving points along vector times the scale factor. Useful for showing flow profiles or mechanical deformation. The filter passes both its point data and cell data to its output.　

　　沿法向膨胀

#!/usr/bin/env pythonimport vtkinputPolyData = vtk.vtkPolyData()sphereSource = vtk.vtkSphereSource()sphereSource.SetPhiResolution(15)sphereSource.SetThetaResolution(15)sphereSource.Update()inputPolyData = sphereSource.GetOutput()# merge duplicate points, and/or remove unused points and/or remove degenerate cellsclean = vtk.vtkCleanPolyData()clean.SetInputData(inputPolyData)# Generate normalsnormals = vtk.vtkPolyDataNormals()normals.SetInputConnection(clean.GetOutputPort())normals.SplittingOff()# vtkWarpVector is a filter that modifies point coordinates by moving points along vector times the scale factor# Warp using the normals (deform geometry with vector data)warp = vtk.vtkWarpVector() warp.SetInputConnection(normals.GetOutputPort())# Set the input data arrays that this algorithm will process# The fieldAssociation refers to which field in the data object the array is storedwarp.SetInputArrayToProcess(0, 0, 0,vtk.vtkDataObject.FIELD_ASSOCIATION_POINTS, vtk.vtkDataSetAttributes.NORMALS)# Specify value to scale displacementwarp.SetScaleFactor(1.0) # Visualize the original and warped modelscolors = vtk.vtkNamedColors()mapper = vtk.vtkPolyDataMapper()mapper.SetInputConnection(warp.GetOutputPort())warpedActor = vtk.vtkActor()warpedActor.SetMapper(mapper)warpedActor.GetProperty().SetColor(colors.GetColor3d("Flesh"))originalMapper = vtk.vtkPolyDataMapper()originalMapper.SetInputConnection(normals.GetOutputPort())originalActor = vtk.vtkActor()originalActor.SetMapper(originalMapper)originalActor.GetProperty().SetInterpolationToFlat() # Set the shading interpolation method for an objectoriginalActor.GetProperty().SetColor(colors.GetColor3d("Flesh"))renderWindow =vtk.vtkRenderWindow()renderWindow.SetSize(640, 480)# Create a camera for all rendererscamera = vtk.vtkCamera()# Define viewport ranges: (xmin, ymin, xmax, ymax)leftViewport = [0.0, 0.0, 0.5, 1.0]rightViewport =[0.5, 0.0, 1.0, 1.0]# Setup both renderersleftRenderer = vtk.vtkRenderer()leftRenderer.SetViewport(leftViewport)leftRenderer.SetBackground(colors.GetColor3d("Burlywood"))leftRenderer.SetActiveCamera(camera)rightRenderer = vtk.vtkRenderer()rightRenderer.SetViewport(rightViewport)rightRenderer.SetBackground(colors.GetColor3d("CornFlower"))rightRenderer.SetActiveCamera(camera)leftRenderer.AddActor(originalActor)rightRenderer.AddActor(warpedActor)rightRenderer.ResetCamera()renderWindow.AddRenderer(rightRenderer)renderWindow.AddRenderer(leftRenderer)style = vtk.vtkInteractorStyleTrackballCamera()interactor = vtk.vtkRenderWindowInteractor()interactor.SetRenderWindow(renderWindow)interactor.SetInteractorStyle(style)renderWindow.Render()interactor.Start()

View Code

　　膨胀导入的STL模型

#!/usr/bin/env pythonimport vtkinputPolyData = vtk.vtkPolyData()readerSTL = vtk.vtkSTLReader()readerSTL.SetFileName("Suzanne.stl")readerSTL.Update()inputPolyData = readerSTL.GetOutput()# merge duplicate points, and/or remove unused points and/or remove degenerate cellsclean = vtk.vtkCleanPolyData()clean.SetInputData(inputPolyData)# Generate normalsnormals = vtk.vtkPolyDataNormals()normals.SetInputConnection(clean.GetOutputPort())normals.SplittingOff() # Turn off the splitting of sharp edges# vtkWarpVector is a filter that modifies point coordinates by moving points along vector times the scale factor# Warp using the normals (deform geometry with vector data)warp = vtk.vtkWarpVector() warp.SetInputConnection(normals.GetOutputPort())# Set the input data arrays that this algorithm will process# The fieldAssociation refers to which field in the data object the array is storedwarp.SetInputArrayToProcess(0, 0, 0,vtk.vtkDataObject.FIELD_ASSOCIATION_POINTS, vtk.vtkDataSetAttributes.NORMALS)# Specify value to scale displacementwarp.SetScaleFactor(0.3) # Visualize the original and warped modelscolors = vtk.vtkNamedColors()mapper = vtk.vtkPolyDataMapper()mapper.SetInputConnection(warp.GetOutputPort())warpedActor = vtk.vtkActor()warpedActor.SetMapper(mapper)warpedActor.GetProperty().SetColor(colors.GetColor3d("Flesh"))originalMapper = vtk.vtkPolyDataMapper()originalMapper.SetInputConnection(normals.GetOutputPort())originalActor = vtk.vtkActor()originalActor.SetMapper(originalMapper)originalActor.GetProperty().SetInterpolationToFlat() # Set the shading interpolation method for an objectoriginalActor.GetProperty().SetColor(colors.GetColor3d("Flesh"))renderWindow =vtk.vtkRenderWindow()renderWindow.SetSize(640, 480)# Create a camera for all rendererscamera = vtk.vtkCamera()# Define viewport ranges: (xmin, ymin, xmax, ymax)leftViewport = [0.0, 0.0, 0.5, 1.0]rightViewport =[0.5, 0.0, 1.0, 1.0]# Setup both renderersleftRenderer = vtk.vtkRenderer()leftRenderer.SetViewport(leftViewport)leftRenderer.SetBackground(colors.GetColor3d("Burlywood"))leftRenderer.SetActiveCamera(camera)rightRenderer = vtk.vtkRenderer()rightRenderer.SetViewport(rightViewport)rightRenderer.SetBackground(colors.GetColor3d("CornFlower"))rightRenderer.SetActiveCamera(camera)leftRenderer.AddActor(originalActor)rightRenderer.AddActor(warpedActor)rightRenderer.ResetCamera()renderWindow.AddRenderer(rightRenderer)renderWindow.AddRenderer(leftRenderer)style = vtk.vtkInteractorStyleTrackballCamera()interactor = vtk.vtkRenderWindowInteractor()interactor.SetRenderWindow(renderWindow)interactor.SetInteractorStyle(style)renderWindow.Render()interactor.Start()

View Code