Как вставить точки, извлеченные из нескольких файлов .vtp, в одну полидату - PullRequest
0 голосов
/ 29 марта 2019

Я сохранил в полиданные 10 точки из среза 10, а в полиданные11 точки из среза11

slice10 = 'Slices\Slice10\Slice10_0_0.vtp'
slice11 = 'Slices\Slice11\Slice11_0_0.vtp'

readerSlice10 = vtk.vtkXMLPolyDataReader()
readerSlice10.SetFileName(slice10)
readerSlice10.Update()

readerSlice11 = vtk.vtkXMLPolyDataReader()
readerSlice11.SetFileName(slice11)
readerSlice11.Update()

polydata10 = vtk.vtkPolyData()
polydata10.SetPoints(readerSlice10.GetOutput().GetPoints())

polydata11 = vtk.vtkPolyData()
polydata11.SetPoints(readerSlice11.GetOutput().GetPoints())

Теперь я хотел, чтобы у меня были единичные полиданные со всеми точками из срезов 10 и срезов 11, как я могу достичьэто?

1 Ответ

0 голосов
/ 10 апреля 2019

Это может быть то, что вы ищете.

import vtk

# Generate two sets of points for this example.
points0 = vtk.vtkPoints()
points1 = vtk.vtkPoints()

points0.InsertNextPoint(1., 0., 0.)
points0.InsertNextPoint(1., 1., 0.)

points1.InsertNextPoint(1., 0., 1.)
points1.InsertNextPoint(1., 1., 1.)

polydata10 = vtk.vtkPolyData()
polydata11 = vtk.vtkPolyData()
polydata10.SetPoints(points0)
polydata11.SetPoints(points1)
#-------------------------------------

# Create a set of points that joints the two sets of points from polydata10 and polydata11.
pointsJoined = vtk.vtkPoints()
for i in range(polydata10.GetNumberOfPoints()):
    pointsJoined.InsertNextPoint(polydata10.GetPoint(i))

for i in range(polydata11.GetNumberOfPoints()):
    pointsJoined.InsertNextPoint(polydata11.GetPoint(i))

# Initialize a polydata object and set the joint point set.
polydata = vtk.vtkPolyData()
polydata.SetPoints(pointsJoined)

# Create verticies so the points can be visualized
verts = vtk.vtkCellArray()
for i in range(polydata.GetNumberOfPoints()):
    verts.InsertNextCell(1) # Create a 1 dimensional cell
    verts.InsertCellPoint(i) # Append point i to the verts array
polydata.SetVerts(verts)

# Setup and run the visualization
ren = vtk.vtkRenderer() #: vtk.vtkRenderer, The renderer for the visualization
renWin = vtk.vtkRenderWindow() #: vtk.vtkRenderWindow, The render window
renWin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor() #: vtk.vtkRenderWindowInteractor, The render window interactor
iren.SetRenderWindow(renWin)

mapper = vtk.vtkPolyDataMapper()
mapper.SetInputData(polydata)
actor = vtk.vtkActor()
actor.GetProperty().SetPointSize(5) # Increase the size of the points
actor.SetMapper(mapper)

ren.AddActor(actor)
iren.Start()
...