我拼凑了 下面的代码用一个三角形网格绘制颜色指定的颜色 额外的标量函数:
#! /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()
这种方法运作得相当好。但是,而不是每个三角形 我会用均匀的颜色和三角形之间的锐利过渡 而是在整个表面上进行平滑插值。
有办法吗?
答案 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')