3

我将以下代码拼凑在一起，以绘制具有附加标量函数指定颜色的三角形网格：

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

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

有没有办法做到这一点？

4

1 回答 1

6

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