Mayavi:在triangular_mesh中插入面部颜色

Question

我拼凑了下面的代码，用另外一个标量函数指定的颜色绘制一个三角形网格：

#! /usr/bin/env python
import numpy as np
from mayavi import mlab

# Create cone
n = 8
t = np.linspace(-np.pi, np.pi, n)
z = np.exp(1j*t)
x = z.real.copy()
y = z.imag.copy()
z = np.zeros_like(x)
triangles = [(0, i, i+1) for i in range(n)]
x = np.r_[0, x]
y = np.r_[0, y]
z = np.r_[1, z]
t = np.r_[0, t]

# These are the scalar values for each triangle
f = np.mean(t[np.array(triangles)], axis=1)

# Plot it
mesh = mlab.triangular_mesh(x, y, z, triangles,
                            representation='wireframe',
                            opacity=0)
cell_data = mesh.mlab_source.dataset.cell_data
cell_data.scalars = f
cell_data.scalars.name = 'Cell data'
cell_data.update()

mesh2 = mlab.pipeline.set_active_attribute(mesh,
        cell_scalars='Cell data')
mlab.pipeline.surface(mesh2)

mlab.show()

这工作得相当好。然而，与其让每个三角形具有均匀的颜色和三角形之间的锐利过渡，我更愿意在整个曲面上进行平滑插值。

有没有办法做到这一点？

Answer 1

我认为您希望使用点数据而不是单元格数据。对于单元格数据，单个标量值不会本地化到任何点。它被指定给整个面。看起来您只想将t数据分配给顶点。点标量的默认渲染将平滑地跨每个面进行插值。

point_data = mesh.mlab_source.dataset.point_data
point_data.scalars = t
point_data.scalars.name = 'Point data'
point_data.update()

mesh2 = mlab.pipeline.set_active_attribute(mesh,
        point_scalars='Point data')