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 に答える