How To Join/connect/group Multiple Objects In Mayavi2

I would like to combine multiple Mayavi objects into a single 'grouped' object so that I control all of their properties together. For example I created the following bi-convex len

Solution 1:

I don't know of a way to combine sources in the way you are looking for. I think in fact that is probably impossible since under the hood the BuiltinSurface object has specific vtk sources that are not what you want. It should however be possible to simply use a different source that gives what you want. In this case you could generate a biconvex lens with mlab.mesh:

a,c,h=3,1,.2

phi,theta = np.mgrid[0:2*np.pi:np.pi/250, 0:2*np.pi:np.pi/250]
x=a*np.cos(theta)*np.sin(phi)
y=a*np.sin(theta)*np.sin(phi)
z=c*np.cos(phi)+(h*(-1)**(np.cos(phi)<0))

mlab.mesh(x,y,z,color=(1,1,1)
mlab.show()

One minor difference is that this surface is smooth. This is the nature of sampling a single surface --i.e., this result is a direct consequence of what your question asks to do. If this is an important feature of your figure, I would suggest an entirely different approach: wrap the 3 sources in a class and have the event handler update the relevant attributes on all three.

Solution 2:

Based on the code by @aestrivex, here is one way of getting the desired output (lens with sharp edges). Note that this is not a solution for connecting multiple Mayavi objects.

import numpy as np
from mayavi import mlab

# Control parameters# r is the semi-diameter of the lens # c controls the center thickness of the lens# h controls the curvature of the surfaces (lesser the value more the curvature)

r, c, h = 3, .75, .9 

delta_phi = np.pi/250.0# phi == azimuth  (0 <= phi <= pi)
delta_theta = np.pi/100.0# theta == zenith (0 <= theta <= pi)
phi, theta = np.mgrid[0:2.0*np.pi + delta_phi:delta_phi,0:np.pi + delta_theta:delta_theta]
# The Exact threshold values for masking tz, txy will change depending upon the# sampling of theta. txy is always slightly less than tz. tz should be around 0.3
tz, txy = 0.279, 0.275
x = r*np.sin(theta)*np.cos(phi)*(np.abs(np.cos(theta)) > txy)
y = r*np.sin(theta)*np.sin(phi)*(np.abs(np.cos(theta)) > txy)
z = c*np.cos(theta)*(h**(-1)*( np.abs(np.cos(theta)) > tz))

mlab.mesh(x,y,z,color=(1,1,1))

mlab.show()

And here is the output: