如何用matplotlib体素绘制直线3D FEM网格？

如何解决如何用matplotlib体素绘制直线3D FEM网格？

我有一个3D直线网格，它具有不同的尺寸宽度。此网格中的每个体素都有一个关联的值。我想使用python 3.7中的matplotlib绘制此网格定义的体素体积。

以一个示例为例，考虑以下代码定义的体素化球体的表面：

import numpy as np

x_0 = 10
y_0 = 10
z_0 = 10

a = 20
b = 20
c = 20

x_i = np.linspace(-40,40,num=100)
# logspace to simulate differing dimension widths
y_i = np.logspace(0,3,num=100)
z_i = np.logspace(0,num=100)

X,Y,Z = np.meshgrid(x_i,y_i,z_i,indexing='ij')

val = (X-x_0)**2/a**2 + (Y-y_0)**2/b**2 + (Z-z_0)**2/c**2

除每个维度的最后一个索引外，假定每个节点对应于体素的x（0），y（0），z（0），角。

如何获取X，Y，Z和val，并绘制（以3D形式）仅对应val