python vtk 样条插值_用python中的vtk插值点

import vtk

import numpy as np

src=np.random.rand(32,3)

vtkPoints=vtk.vtkPoints()

for i in range(0,src.shape[0]):

x=src[i,0]

y=src[i,1]

z=src[i,2]

p=[x,y,z]

vtkPoints.InsertNextPoint(p)

poly = vtk.vtkPolyData()

poly.SetPoints(vtkPoints)

#Triangulate the grid points.

#Create a grid of points to interpolate over

gridPoints =vtk.vtkPoints()

for x in range (0,src.shape[0]):

for y in range (0,src.shape[1]):

gridPoints.InsertNextPoint(x,y,0)

#Triangulate the grid points.

vtkDelaunay=vtk.vtkDelaunay2D()

vtkDelaunay.SetInputData(poly)

vtkDelaunay.Update()

#Create a dataset from the grid points

gridPolyData =vtk.vtkPolyData()

gridPolyData.SetPoints(gridPoints)

#Perform the interpolation

probeFilter =vtk.vtkProbeFilter()

probeFilter.SetSourceConnection(vtkDelaunay.GetOutputPort())

probeFilter.SetInputData(gridPolyData)

probeFilter.Update()

gridWarpScalar =vtk.vtkWarpScalar()

gridWarpScalar.SetInputConnection(probeFilter.GetOutputPort())

gridWarpScalar.Update()

#contour filter:

surf = vtk.vtkSurfaceReconstructionFilter()

surf.SetInputData(poly)

surf.GetOutput().GetPointData().SetScalars(weights)

surf.Update()

contourFilter = vtk.vtkContourFilter()

contourFilter.SetInputConnection(surf.GetOutputPort())

contourFilter.SetValue(0, 0.0)

contourFilter.Update()

# visualization

Mapper =vtk.vtkPolyDataMapper()

Mapper.SetInputConnection(contourFilter.GetOutputPort())

Mapper.ScalarVisibilityOff()

Actor=vtk.vtkActor()

Actor.SetMapper(Mapper)

Actor.GetProperty().SetColor(1.0, 0, 0)

Actor.GetProperty().SetPointSize(3)

gridMapper =vtk.vtkDataSetMapper()

gridMapper.SetInputConnection(vtkDelaunay.GetOutputPort())

gridMapper.ScalarVisibilityOff()

gridActor =vtk.vtkActor()

gridActor.SetMapper(gridMapper);

gridActor.GetProperty().SetColor(0, 0, 1.0)

gridActor.GetProperty().SetPointSize(3)

renderer =vtk.vtkRenderer()

renderWindow =vtk.vtkRenderWindow()

renderWindow.AddRenderer(renderer)

renderWindowInteractor =vtk.vtkRenderWindowInteractor()

renderWindowInteractor.SetRenderWindow(renderWindow)

style =vtk.vtkInteractorStyleTrackballCamera()

renderWindowInteractor.SetInteractorStyle(style)

renderer.AddActor(Actor)

renderer.AddActor(gridActor)

renderWindow.Render()

renderWindowInteractor.Start()