<div dir="ltr">Dear Elaine and Tom,<div><br></div><div>Thanks so much for the clarification and detailed information with the links of the python scripts.</div><div>The algorithm becomes more refined. I will work though your emails step by step. although there </div><div>is a learning curve.  I understand euler angle can be painful to work with sometimes, thanks for </div><div>suggesting to convert them into a 3 by 3 rotation matrix. I will see how to do it. </div><div></div><div><br></div><div>Best regards,</div><div><br></div><div>Wei</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 2, 2020 at 3:14 PM Tom Goddard <<a href="mailto:goddard@sonic.net">goddard@sonic.net</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">Hi Wei,<div><br></div><div>  If you placed the copies of the subtomogram average, you could then use the Chimera "vop add" command to add them all together (or maybe better to use "vop maximum" so that the overlaps between two copies don't create double density) making a single map.  Then the "vop localCorrelation" command could compute the correlation between this assembled map and the original tomogram and color according to that correlation using the Chimera "Surface Color" tools to color by volume data value (using the local correlation map as the value).</div><div><br></div><div>  Chimera does almost nothing with Euler angles, and there are different Euler angle conventions you would need to worry about.  You would be better off if you had a 3x3 rotation matrix.  There is some Euler angle example of rotating a map and resampling it in the Chimera Python scripts web page but that seems like a hard approach.</div><div><br></div><div><span style="white-space:pre-wrap"> </span><a href="http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/eulermove.py" target="_blank">http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/eulermove.py</a></div><div><br></div><div>  Tom</div><div><br></div><div><div><br><blockquote type="cite"><div>On Jan 2, 2020, at 10:08 AM, Wei Zhang <<a href="mailto:zhangwei@umn.edu" target="_blank">zhangwei@umn.edu</a>> wrote:</div><br><div><div dir="ltr">Elanie,<div><br></div><div>These functions are very very useful. Thank you.</div><div>I believe we can take advantage of these built-in functions in our procedures.</div><div><br></div><div>I have a follow-up question. If I have map (in mrc format), I wish to display it in Chimera at specific x, y, z position and at a certain Eular angle. How can I do it?</div><div>I can think of two ways:</div><div>(1) Modify the header of the map to include the position and orientation of the map, so Chimera can read and display it accordingly</div><div>(2) Have a special text file that include those information and have Chimera display it </div><div><br></div><div>Is there a better way to do it?</div><div><br></div><div>Thanks so much!</div><div><br></div><div>Wei</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 2, 2020 at 11:56 AM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Hi Wei,<br>
Happy 2020!<br>
<br>
Chimera has a “vop localCorrelation” command that given two maps, creates a new map that is the local correlation (using a sliding box, for which the user can specify the box size).  Then the local-correlation map values can be used to color the isosurface of another map using the “Surface Color” tool or “scolor” command, and also show a color key.<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/vop.html#localCorrelation" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/vop.html#localCorrelation</a>><br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/surfcolor/surfcolor.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/surfcolor/surfcolor.html</a>><br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/scolor.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/scolor.html</a>><br>
<br>
ChimeraX also has these features, as commands “volume localCorrelation” and “color sample”, but cannot yet draw the color key.<br>
<<a href="http://rbvi.ucsf.edu/chimerax/docs/user/commands/volume.html#localCorrelation" rel="noreferrer" target="_blank">http://rbvi.ucsf.edu/chimerax/docs/user/commands/volume.html#localCorrelation</a>><br>
<<a href="http://rbvi.ucsf.edu/chimerax/docs/user/commands/color.html#map" rel="noreferrer" target="_blank">http://rbvi.ucsf.edu/chimerax/docs/user/commands/color.html#map</a>><br>
<br>
I hope this helps,<br>
Best,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D.<br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Jan 1, 2020, at 8:40 AM, Wei Zhang <<a href="mailto:zhangwei@umn.edu" target="_blank">zhangwei@umn.edu</a>> wrote:<br>
> <br>
> Hello,<br>
> Happy 2020!<br>
> <br>
> We are doing tomography reconstruction of retrovirus assemblies. Because the raw tomograms are very noisy, they can not be represented nicely using isosurfaces. I see other researchers in the field presented the tomo map in Chimera by superimposing subtomogram-averaged map onto the original tomogram according to the subunit's locations and  relative orientations. The color of the displayed subunits is based on the correlation coefficient between the averaged subunit map and the original density of the tomogram.<br>
> <br>
> One example is Fig.1B in this paper:<br>
> <a href="https://advances.sciencemag.org/content/5/11/eaaw3631" rel="noreferrer" target="_blank">https://advances.sciencemag.org/content/5/11/eaaw3631</a><br>
> <br>
> We do use the program Dynamo for the subtomogram averaging procedure, so our data format is very similar to what was used in this paper. Is this feature is part of Chimera? Could we obtain the extended software? Or is this customer built and we need to contact the authors of the paper? Or shall we work out this procedure by ourselves with the help of Chimera experts?<br>
> <br>
> Thank you,<br>
> Wei<br>
<br>
</blockquote></div>
_______________________________________________<br>Chimera-dev mailing list<br><a href="mailto:Chimera-dev@cgl.ucsf.edu" target="_blank">Chimera-dev@cgl.ucsf.edu</a><br><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev</a><br></div></blockquote></div><br></div></div></blockquote></div>