如何在mayavi中正確顯示包含非立方體素的體積

Question

我正在使用mayavi（3.3.2）顯示體積等值面。

通常，我的體積沒有立方體素;例如，採樣網格在X和Y方向上可能爲1mm x 1mm，但在Z方向上爲1.4mm。

如何使用mayavi的mlab.contour3d或mlab.pipeline.iso_surface以正確的空間比例顯示這樣的體積？我真的不想將這些卷重新採樣到立方格。

陳述問題的另一種方式：我可以做些什麼來獲得下面的代碼來顯示一個球體而不是展平的elipsoid（將volume與其預期的1：1：2縱橫比體素作爲給定值，並且無需重新生成或重新採樣該卷）。

import numpy as np 
from enthought.mayavi import mlab 

def sqr(x): return x*x 

s=64 
x,y,z = np.ogrid[0:s,0:s,0:s/2] 

volume = np.sqrt(sqr(x-s/2)+sqr(y-s/2)+sqr(2*z-s/2)) 

isos = mlab.contour3d(volume,contours=[5,15,25],transparent=True) 
mlab.show() 

我猜，就必須有在底層VTK圖形管道獲得的某種方式（其變換等），並插入相應的各向異性縮放（如果更直接地通過沒有這樣做的一些方法mlab API）。

Answer 1

爲此，從輸入數據中明確創建scalar_field對象是最容易的。

實際上，我經常這樣做，因爲我們喜歡在地質學中深入研究（正向下）。這意味着你需要在Z方向上有負增量。如果它僅僅是各種mlab函數的參數，那將是非常好的，但它仍然不是很難做到。

from mayavi import mlab 
import numpy as np 

s=64 
x,y,z = np.ogrid[0:s,0:s,0:s/2] 

data = np.sqrt((x-s/2)**2 + (y-s/2)**2 + (2*z-s/2)**2) 

grid = mlab.pipeline.scalar_field(data) 
grid.spacing = [1.0, 1.0, 2.0] 

contours = mlab.pipeline.contour_surface(grid, 
         contours=[5,15,25], transparent=True) 
mlab.show()