Mayavi:在triangular_mesh中插入面部颜色

时间:2013-02-19 17:57:38

标签: python enthought mayavi

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

#! /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()

这种方法运作得相当好。但是,而不是每个三角形 我会用均匀的颜色和三角形之间的锐利过渡 而是在整个表面上进行平滑插值。

有办法吗?

1 个答案:

答案 0 :(得分:5)

我认为您想使用点数据而不是单元格数据。对于单元格数据,单个标量值不会本地化到任何点。它被分配给整个面部。看起来您只想将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')