独立着色和旋转3D VTK图示符

Question

我尝试使用pythonwrapper for VTK在一个场景中显示多个glyph，每个glyph都有自己的颜色和旋转。不幸的是，vtk只考虑了旋转（使用vtkTensorGlyph）。相反，当我使用vtkGlyph3D时，只考虑颜色。在

下面是一段带有vtkTensorGlyph的现成代码。每个立方体应该有一个随机的颜色，但是它们都是相同的颜色。我反复阅读了VTK的文档，但没有找到任何解决方案。在

#!/usr/bin/env python
# -*- coding: utf-8 -*-

import vtk
import scipy.linalg as sc
import random as ra
import numpy as np
import itertools

points = vtk.vtk.vtkPoints()  # where to locate each glyph in the scene
tensors = vtk.vtkDoubleArray() # rotation for each glyph
tensors.SetNumberOfComponents(9)
colors = vtk.vtkUnsignedCharArray() # should be the color for each glyph
colors.SetNumberOfComponents(3)

# let's make 10 cubes in the scene
for i in range(0, 50, 5): 
    points.InsertNextPoint(i, i, i) # position of a glyph
    colors.InsertNextTuple3(ra.randint(0, 255), ra.randint(0, 255), ra.randint(0, 255) ) # pick random color
    rot = list(itertools.chain(*np.reshape(sc.orth(np.random.rand(3, 3)).transpose(), (1, 9)).tolist())) # random rotation matrix (row major)
    tensors.InsertNextTuple9(*rot)

polydata = vtk.vtkPolyData() # create the polydatas
polydata.SetPoints(points) 
polydata.GetPointData().SetTensors(tensors)
polydata.GetPointData().SetScalars(colors)

cubeSource = vtk.vtkCubeSource()
cubeSource.Update()

glyphTensor = vtk.vtkTensorGlyph()
glyphTensor.SetColorModeToScalars() # is it really work ?

try:
    glyphTensor.SetInput(polydata)
except AttributeError:
    glyphTensor.SetInputData(polydata)

glyphTensor.SetSourceConnection(cubeSource.GetOutputPort())
glyphTensor.ColorGlyphsOn() # should not color all cubes independently ?
glyphTensor.ThreeGlyphsOff()
glyphTensor.ExtractEigenvaluesOff()
glyphTensor.Update()

# next is usual vtk code
mapper = vtk.vtkPolyDataMapper()
mapper.SetInputConnection(glyphTensor.GetOutputPort())

actor = vtk.vtkActor()
actor.SetMapper(mapper)

ren = vtk.vtkRenderer()
ren.SetBackground(0.2, 0.5, 0.3)
ren.AddActor(actor)

renwin = vtk.vtkRenderWindow()
renwin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor()
iren.SetInteractorStyle(vtk.vtkInteractorStyleTrackballCamera())
iren.SetRenderWindow(renwin)

renwin.Render()
iren.Initialize()
renwin.Render()
iren.Start()

Answer 1

您可以按以下方式为vtktensorglyph涂上不同的颜色。我会用不同颜色的张量字形来展示我的C++代码的一些片段。我希望您可以将类似的东西集成到您的python代码中：

vtkSmartPointer<vtkFloatArray> colors = vtkSmartPointer<vtkFloatArray>::New();
vtkSmartPointer<vtkPolyData> data = vtkSmartPointer<vtkPolyData>::New();

//Add for each single object the following data to colors: 
float color[3] = {r, g, b};
colors->InsertNextTupleValue(color);

//Then proceed here after adding all colors
colors->SetName("colors");
colors->SetNumberOfComponents(3);

//Add all necessary information to the vtkPolyData element, including the color information
data->GetPointData()->SetScalars(colors);

...

tensorGlyph->SetInputData(data);
tensorGlyph->ColorGlyphsOn();
tensorGlyph->SetColorModeToScalars();

虽然VTK开发人员可能打算使用RGB元组来定义颜色，但只需使用第一个参数“r”。在[0,1]的间隔内改变“r”将导致不同的颜色。在

使用VTK 6.1.0

Answer 2

我终于想办法做点什么了。 如果我很清楚，字形需要一个曲面，尤其是构成该曲面的每个点上的法线。每个法线给出每个字形的方向。 vtk源代码中Examples/VisualizationAlgorithms/Python/spikeF.py处的示例说明了我的意思。 下面是一个单独着色和旋转字形的小例子。在

#!/usr/bin/env python
# -*- coding: utf-8 -*-

from vtk import *

# Let's create a surface with 3 points unconnected.
points = vtkPoints()
points.InsertNextPoint(1, 0, 0)
points.InsertNextPoint(0, 0, 0)
points.InsertNextPoint(0, 1, 0)

polydata = vtkPolyData()
polydata.SetPoints(points)

colors = vtkUnsignedCharArray()
colors.SetNumberOfComponents(3)
colors.SetNumberOfTuples(polydata.GetNumberOfPoints())

colors.InsertTuple3(0, 147, 25, 98)
colors.InsertTuple3(1, 32, 84, 247)
colors.InsertTuple3(2, 198, 214, 36)

polydata.GetPointData().SetScalars(colors)

# Now, let's control normal on each point composing or 'fake' surface
# We should say, let's give a direction to each point, a normal in strange for a point.
pointNormalsArray = vtkDoubleArray()
pointNormalsArray.SetNumberOfComponents(3)
pointNormalsArray.SetNumberOfTuples(polydata.GetNumberOfPoints())

pN1 = [1.0, 0.0, 0.0]
pN2 = [0.0, 1.0, 0.0]
pN3 = [0.0, 0.0, 1.0]

pointNormalsArray.SetTuple(0, pN1)
pointNormalsArray.SetTuple(1, pN2)
pointNormalsArray.SetTuple(2, pN3)

polydata.GetPointData().SetNormals(pointNormalsArray)

sources = vtkConeSource()
sources.SetResolution(6)
sources.Update()

glyph = vtkGlyph3D()
try:
    glyph.SetInput(polydata)
except AttributeError:
    glyph.SetInputData(polydata)

glyph.SetSourceConnection(sources.GetOutputPort())
glyph.SetColorModeToColorByScalar()
glyph.SetVectorModeToUseNormal()
glyph.ScalingOff()
glyph.Update()

mapper = vtkPolyDataMapper()
mapper.SetInputConnection(glyph.GetOutputPort())

actor = vtkActor()
actor.SetMapper(mapper)

ren = vtkRenderer()
ren.SetBackground(0.2, 0.5, 0.3)
ren.AddActor(actor)

renwin = vtk.vtkRenderWindow()
renwin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor()
iren.SetInteractorStyle(vtk.vtkInteractorStyleTrackballCamera())
iren.SetRenderWindow(renwin)

renwin.Render()
iren.Initialize()
renwin.Render()
iren.Start()

顺便说一句，还有一个更简单的解决方案，即创建一个源，以及所需的任何映射器和参与者。事实上，给每个演员一个颜色和方向是很容易的。问题是，对于很多演员（我的机器上大约有5000场）来说，表演对互动是很不利的。在