如何在基于截视图的python脚本中显示vtkUnstructuredGrid？
提问于 2019-02-08 20:45:21
我在我的Ubuntu18.04系统上安装了段视图5.6，我想编写一个python脚本来显示vtkUnstructuredGrid。
import numpy as np
from paraview.simple import *
import paraview.vtk as vtk
from paraview.vtk.numpy_interface import dataset_adapter as dsa
import paraview.vtk.util.numpy_support as vnp
node = np.array(
[[0.0, 0.0, 0.0],
[1.0, 0.0, 0.0],
[1.0, 1.0, 0.0],
[0.0, 1.0, 0.0]], dtype=np.float)
cell = np.array([[3, 1, 2, 0], [3, 3, 0, 2]], dtype=np.int)
NC = cell.shape[0]
points = vtk.vtkPoints()
points.SetData(vnp.numpy_to_vtk(node))
cells = vtk.vtkCellArray()
cells.SetCells(NC, vnp.numpy_to_vtkIdTypeArray(cell))
uGrid =vtk.vtkUnstructuredGrid()
uGrid.SetPoints(points)
uGrid.SetCells(vtk.VTK_TRIANGLE, cells)
# how to put uGrid into the following codes
view = GetActiveViewOrCreate('RenderView')
dispaly = Show()
render = Render()
Interact()在python脚本中，我无法从互联网上找到任何这样的例子。所以我需要你的帮助，非常感谢。
更新：
我尝试编写一个源代码类，如下所示：
import numpy as np
from paraview.simple import *
import vtk
import vtk.util.numpy_support as vnp
from vtkmodules.util.vtkAlgorithm import VTKPythonAlgorithmBase
from vtkmodules.numpy_interface import dataset_adapter as dsa
from paraview.util.vtkAlgorithm import smproxy, smproperty, smdomain
@smproxy.source(name="MeshSource", label="triangle mesh!")
class MeshSource(VTKPythonAlgorithmBase):
def __init__(self):
print("Initialize the source!")
VTKPythonAlgorithmBase.__init__(self,
nInputPorts=0,
nOutputPorts=1,
outputType='vtkUnstructuredGrid')
node = np.array(
[[0.0, 0.0, 0.0],
[1.0, 0.0, 0.0],
[1.0, 1.0, 0.0],
[0.0, 1.0, 0.0]], dtype=np.float)
cell = np.array([[3, 1, 2, 0], [3, 3, 0, 2]], dtype=np.int)
NN = node.shape[0]
NC = cell.shape[0]
points = vtk.vtkPoints()
points.SetData(vnp.numpy_to_vtk(node))
cells = vtk.vtkCellArray()
cells.SetCells(NC, vnp.numpy_to_vtkIdTypeArray(cell))
self.mesh = vtk.vtkUnstructuredGrid()
self.mesh.SetPoints(points)
self.mesh.SetCells(vtk.VTK_TRIANGLE, cells)
rho = vnp.numpy_to_vtk(np.zeros(NN))
rho.SetName('rho_A')
self.mesh.GetPointData().AddArray(rho)
self.Port = 0
def RequestData(self, request, inInfo, outInfo):
print("Request the data!")
output = vtk.vtkUnstructuredGrid.GetData(outInfo)
optput.ShallowCopy(self.mesh)
return 1
def UpdatePointData(self, rho):
print("Update the point data!")
rho = vnp.numpy_to_vtk(rho)
rho.SetName('rho_A')
self.mesh.GetPointData().AddArray(rho)
self.Modified()
source = MeshSource()
view = GetActiveViewOrCreate('RenderView')
display = Show(source, view)
Interact()但我有一些错误：
Traceback (most recent call last):
File "test_triangle.py", line 55, in <module>
dispaly = Show(source, view)
File "/home/why/local/lib/python3.6/site-packages/paraview/simple.py", line 482, in Show
rep = controller.Show(proxy, proxy.Port, view)
File "/home/why/local/lib/python3.6/site-packages/paraview/servermanager.py", line 158, in __ConvertArgumentsAndCall
retVal = func(*newArgs)
TypeError: Show argument 1: method requires a vtkSMSourceProxy, a vtkPythonAlgorithm was provided.我一定是错过了什么。
回答 1
Stack Overflow用户
回答已采纳
发布于 2019-02-11 07:24:21
需要理解的关键是，在ParaView中有两个级别的Python可用。较低级别是您可以使用VTK创建或筛选数据的地方。较高级别允许您控制ParaView的操作，例如显示数据、设置显示属性等。您所缺少的是脚本中两个级别之间的桥梁。
在您最初的示例中，您正在VTK中创建一个非结构化网格。要将它带到ParaView可以使用它的地方，请添加以下内容：
# how to put uGrid into the following codes
view = GetActiveViewOrCreate('RenderView')
# create a trivial producer to bridge between the VTK object and ParaView
tp = TrivialProducer()
tp.GetClientSideObject().SetOutput(uGrid)
dispaly = Show(tp)这将创建一个ParaView代理(TrivialProducer)，它是一个名为vtkTrivialProducer的简单VTK数据源的代理。它所做的就是获取一个数据集，并将其传递给请求它的任何下游过滤器。(注意:这只在您在内置服务器模式下运行时起作用，因为GetClientSideObject()是相当常见的)。
定义单元格的方法也有问题。单元格定义中的第一个条目需要定义单元格的点数。所以把这一行改为
cell = np.array([[3, 1, 2, 0], [3, 3, 0, 2]], dtype=np.int)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/54603267
复制相关文章
腾讯云开发者
