【问题标题】:Mayavi: interpolate face colors in triangular_meshMayavi:在 triangular_mesh 中插入面部颜色
【发布时间】:2013-02-19 17:57:38
【问题描述】:

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

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

这工作得相当好。但是,不是每个三角形 在三角形之间具有统一的颜色和清晰的过渡,我会 而是在整个表面上进行平滑插值。

有没有办法做到这一点?

【问题讨论】:

    标签: python enthought mayavi


    【解决方案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')
    

    【讨论】:

    • 这听起来很有希望,但我的标量数据实际上是每个单元格而不是每个点。我是否必须先手动插值才能获取顶点处的值,或者有没有办法为例如指定点数据?单元中心而不是顶点?
    • 在这种情况下,我不清楚您的问题是否恰当。与所有这些面相邻的示例圆锥的峰值处的插值是多少?什么算法可以从中心、面值平滑地插值到那个?您可以考虑将点标量分配给现有顶点,这些顶点是每个相邻单元标量的平均值。不过,这并不一定要在每个面的中心重新创建纯单元标量。
    • 好吧,根据定义,插值从来都不是唯一的(如果这就是你的意思)。关于中心面与峰值顶点,我没有看到问题。在绘制 3d 结构时,我的标量是在 2d 表面上定义的,因此任何 2d 插值算法都应该可以。我错过了什么吗?
    • 2D 表面上的插值仍然需要将值分配给点,而不是面。因此,如果您可以重组数据和/或网格,使您的值位于顶点上,VTK 就会这样做。
    猜你喜欢
    • 2014-12-13
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2013-09-03
    • 1970-01-01
    • 2022-06-13
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多