<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Oleksiy,<div class=""><br class=""></div><div class=""> I’m not sure the cause of your MATLAB positioning problem, but here is an explanation of the Chimera matrices that may be helpful.</div><div class=""><br class=""></div><div class=""> First the Chimera command matrixget gives the position in the Chimera scene including any mouse rotation and translation you did to view the models. So you would have to consider positioning matrices for both maps. Even when you first open the maps and don’t move them, then fit one in the other, both maps probably have non-zero translations since the initial centering of the view on the first map does a translation. So using matrixget is not a good way to get the position of one map relative to another, unless you are going to invert one of the matrices and multiply times the other. Instead use the “measure rotation” command (e.g. measure rotation #0 #1) which reports a matrix in the reply log giving the motion of map #1 relative to the coordinates of map #2.</div><div class=""><br class=""></div><div class=""> Another thing to be aware of is that these 3x4 matrices contain a rotation (the first 3 columns) and a translation (4th column) and the rotation is applied first, followed by the translation. The rotation is about center x,y,z = (0,0,0) but that is not necessarily grid point i=j=k=0 because the MRC file can have an origin offset. If you look in the Chimera volume viewer dialog coordinates panel (volume viewer menu Features / Coordinates) the origin grid index is given as “Origin index”. That is the grid point the rotation is centered on.</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Tom</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Apr 24, 2017, at 9:30 AM, Oleksiy Kovtun <<a href="mailto:okovtun@mrc-lmb.cam.ac.uk" class="">okovtun@mrc-lmb.cam.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><br class="">
<br class="">
<br class="">
<br class=""><p class="">Dear all, <br class="">
I need to merge two maps in MATLAB using<br class="">
alignment data from chimera. Inputs: two maps of the same box size,
origin. I use matrixget - to report translation and rotation<br class="">
matrices for the aligned map . Translation operation works well when
transferred into MATLAB (briefly, translation matrix gets converted in
eulers and pased on to tom_rotate command in MATLAB, TOM package). However
the x y z translation coordinates (the fourth column top to bottom in
matrixget - output) are way off the real translation when using tom_shift
command in MATLAB, TOM.<br class="">
I found out that if I open the rotated map<br class="">
and fit it again, it gives close to correct translation<br class="">
matrix.<br class="">
<br class="">
What's going wrong here? Does anyone have a<br class="">
solution how to overlay 3D arrays in MATLAB using chimera-reported
transformation matrices?<br class="">
<br class="">
Many<br class="">
thanks,<br class="">
Oleksiy <br class="">
<br class="">
<br class="">
<br class="">
Dr Oleksiy Kovtun<br class="">
<br class="">
Research Fellow<br class="">
<br class="">
Briggs group<br class="">
<br class="">
MRC Laboratory of Molecular Biology<br class="">
<br class="">
Francis Crick Avenue<br class="">
<br class="">
Cambridge<br class="">
<br class="">
CB2 0QH<br class="">
<br class="">
UK<br class="">
<br class="">
<br class="">
<br class="">
Tel: + 44 1223 267551<br class="">
<br class="">
Email: <a href="mailto:okovtun@mrc-lmb.cam.ac.uk" class="">okovtun@mrc-lmb.cam.ac.uk</a><br class="">
<br class="">
<br class="">
<br class="">
Dr Oleksiy Kovtun<br class="">
Research Fellow<br class="">
Briggs group<br class="">
MRC Laboratory of Molecular Biology<br class="">
Francis Crick Avenue<br class="">
Cambridge<br class="">
CB2 0QH<br class="">
UK<br class="">
<br class="">
Tel: + 44 1223 267551<br class="">
Email: <a href="mailto:okovtun@mrc-lmb.cam.ac.uk" class="">okovtun@mrc-lmb.cam.ac.uk</a></p><span id="cid:3A600793-90D6-41D3-AF80-5FE0A07DD3DB@cgl.ucsf.edu"><untitled-[2].html></span>_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></blockquote></div><br class=""></div></body></html>