[Chimera-users] transform coordinates of multiple volumes
meng at cgl.ucsf.edu
Fri Apr 27 10:26:40 PDT 2012
On Apr 27, 2012, at 9:26 AM, Ludovic Renault wrote:
> Hi Elaine,
> I have 120 or so 3d volumes (EM maps). They are all equivalent and are centered at the same position.
> I would like to move each one of these 120 volumes by a specific value in terms of euler angles and shift (in pixels I assume). I can use the 'trAnsform coordinates' tool but then if I save the map, it is not saved with the new coordinates.
> 1- How to (Or can I) save the new coordinates for each volume?
> 2- Can I use command lines to automate the process for example by giving the command a text file with the specific coordinates for each volume?
> I did do the rotation with spider software, but the shifting doesn't work because the volume(s) go outside of the box. So if all the above is not doable, is it possible to increase the volumes boxe size in chimera, then I could use the shift function in spider?
> I hope I'm clear enough.
It is recommended to send Chimera questions to the chimera-users list (CC'd) -- I'm not the top expert on volumes, so others may have better answers!
I think the shift in Transform Coordinates is not in pixels but in units of length.
How the volume grid units map to units of length can be viewed by using Features...Coordinates in the Volume Viewer menu.
Chimera has commands for translation (move) and rotation (turn) in X,Y,Z coordinate systems, but I don't know of a command to apply Euler rotations.
About saving a transformed map:
The map formats don't include transformation information, so there is no way to save EXACTLY the same data, except transformed. However, you can resample the map onto a new grid, and then save the new map. The values at the new grid points are determined by interpolation, so there will be some (generally minor) data loss.
The resampling of one map onto the grid of another can be done with command "vop resample", for example:
vop #1 resample onGrid #0
To save the new resampled volume to a file, you could use "File... Save Map as" in the Volume Viewer dialog or the command "volume" with its "save" option.
However, I see that "vop resample" is most pertinent for saving a map that has been fitted to another. In your case you are simply transforming the map, so it is unclear on what grid you could resample. Maybe there is some way to make a "dummy" grid in the right place.
Although "vop cover" extends a map, it assumes periodicity of that map, whereas I assume you want zeroes in the expanded parts. Perhaps expanding a map with zeroes could be done using "vop add" to add a dummy zero map to your data set of interest.
Sorry I couldn't give a more complete solution,
Elaine C. Meng, Ph.D.
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco
More information about the Chimera-users