Mayavi组合两个隐式3d表面

我使用Mayavi来绘制隐式3d表面，我想将这些表面组合成一个绘图。然而，当我这样做，我得到类似这样的：Mayavi组合两个隐式3d表面

我使用的代码：

import numpy as np 
from mayavi import mlab 

A=0.24639243776 
B=5.39100472027e-17 
C=1.71555149594 
D=1.72967325617 
E=7.50535440036 
F=-1.17072847143 
G=1.0 

x, y, z = np.ogrid[-3:1:100j, -10:1:100j, 0:3:100j] 

def Fun (x,y,z, A,B,C,D,E,F,G): 
    F1 = (x - A + y - B)/2 + np.sqrt(((x - A - y + B)/2) ** 2 + C * z ** 2) 
    F2 = np.sqrt(x*((F*y)/2 + (E*x)/D) + y*((F*x)/2 + (D*y)/E) + G*z**2) - np.sqrt(D*E) 
    F1 [F1 < 0] = F2 [F1 < 0] 
    return F1 

Fu = Fun (x,y,z,A,B,C,D,E,F,G) 

mlab.contour3d(Fu, contours = [0]) 
mlab.show()

的问题是，无论我如何做到这一点的一个表面上是非常低聚！我曾策划在数学相同，它产生一个完美的情节：

我不指望能得到相同质量的情节一样数学是一个商业产品。但至少我希望两面平滑。

如果Mayavi不是正确的工具，也许你可以推荐一些不同的东西来达到预期的效果。

提前致谢！

来源

2016-11-07 UN4

我想通了。关键是将曲面设置为无，以便修剪。然后可以以相同的方式绘制另一个表面。此外，对于表面根据规模，X，Y被绘制，Z也必须被传递到contour3d

import numpy as np 
from mayavi import mlab 

A=0.24639243776 
B=5.39100472027e-17 
C=1.71555149594 
D=1.72967325617 
E=7.50535440036 
F=-1.17072847143 
G=1.0 

x, y, z = np.mgrid[-3:1:300j, -10:1:300j, 0:3:300j] 

def Fun1 (x,y,z, A,B,C,D,E,F,G): 
    F1 = (x - A + y - B)/2 + np.sqrt(((x - A - y + B)/2) ** 2 + C * z ** 2) 
    F2 = np.sqrt(x*((F*y)/2 + (E*x)/D) + y*((F*x)/2 + (D*y)/E) + G*z**2) - np.sqrt(D*E) 
    F1 [F2 > 0.0] = None 
    return F1 


def Fun2 (x,y,z, A,B,C,D,E,F,G): 
    F1 = (x - A + y - B)/2 + np.sqrt(((x - A - y + B)/2) ** 2 + C * z ** 2) 
    F2 = np.sqrt(x*((F*y)/2 + (E*x)/D) + y*((F*x)/2 + (D*y)/E) + G*z**2) - np.sqrt(D*E) 
    F2 [F1 > 0.0] = None 
    return F2 

Fu1 = Fun1 (x,y,z,A,B,C,D,E,F,G) 
Fu2 = Fun2 (x,y,z,A,B,C,D,E,F,G) 

mlab.contour3d(x,y,z,Fu1, contours = [0]) 
mlab.contour3d(x,y,z,Fu2, contours = [0]) 
mlab.show()

这里是结果：

来源

2016-11-11 09:32:25 UN4

Mayavi组合两个隐式3d表面

回答

相关问题