<div dir="ltr"><div>mplot3d is good for real simple 3d scenes, but it falls apart when you start composing things like this. You might want to consider using something like glumpy or mayavi for better results. There are some discussions in the past on the problems inherent in creating spheres in mplot3d. There are some interesting solutions, but they get fairly complex pretty quickly.</div><div><br></div><div>Cheers!<br></div><div>Ben Root<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 28, 2018 at 8:16 AM, Bala subramanian <span dir="ltr"><<a href="mailto:bala.biophysics@gmail.com" target="_blank">bala.biophysics@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><p style="margin-bottom:0in;line-height:100%">Dear mpl friends,</p>
<p style="margin-bottom:0in;line-height:100%">

</p>
<p style="margin-bottom:0in;line-height:100%">My objective is to
create a molecular surface using matplotlib. I used the snippet
pasted below to create a series of spheres (attached png), drawn
using the atom attributes (their positions and radius). From here, I
would like to create a smooth surface which can depict the shape of
the molecule surface (curves,cavities etc). Based on examples in net, I get a mixed feeling that either triangulation or
making 3D grids are some directions to proceed. If anyone has some
suggestions or better protocols for such task with some example data, it would be
of great help to me.</p>
<p style="margin-bottom:0in;line-height:100%">

</p>
<p style="margin-bottom:0in;line-height:100%">fig =
plt.figure(figsize=(6,6),dpi=1<wbr>40)</p>
<p style="margin-bottom:0in;line-height:100%">ax1 =
fig.add_subplot(111, projection='3d')</p>
<p style="margin-bottom:0in;line-height:100%">

</p>
<p style="margin-bottom:0in;line-height:100%"># Define spherical
coordinates</p>
<p style="margin-bottom:0in;line-height:100%">u, v =
np.mgrid[0:2*np.pi:20j, 0:np.pi:10j]</p>
<p style="margin-bottom:0in;line-height:100%">S =
np.cos(u)*np.sin(v)</p>
<p style="margin-bottom:0in;line-height:100%">C =
np.sin(u)*np.sin(v)</p>
<p style="margin-bottom:0in;line-height:100%">T = np.cos(v)</p>

<p style="margin-bottom:0in;line-height:100%">

</p>
<p style="margin-bottom:0in;line-height:100%">for P,R in
zip(positions,radii):</p><p style="margin-bottom:0in;line-height:100%">         ax1.plot_wireframe(R*S+P[0],R*<wbr>C+P[1],R*T+P[2],color='r',alph<wbr>a=0.4)</p>

<br clear="all">


        
        <span></span>
        
        


<p style="margin-bottom:0in;line-height:100%">Thanks in advance,</p>
<p style="margin-bottom:0in;line-height:100%">Bala</p><span class="gmail-HOEnZb"><font color="#888888"><br clear="all"><br>-- <br><div class="gmail-m_313166897354034870gmail_signature">C. Balasubramanian<br></div>
</font></span></div>
<br>______________________________<wbr>_________________<br>
Matplotlib-users mailing list<br>
<a href="mailto:Matplotlib-users@python.org">Matplotlib-users@python.org</a><br>
<a href="https://mail.python.org/mailman/listinfo/matplotlib-users" rel="noreferrer" target="_blank">https://mail.python.org/<wbr>mailman/listinfo/matplotlib-<wbr>users</a><br>
<br></blockquote></div><br></div>