Hi Catalin,
The Chimera molmap command creates a density map for a molecule using a grid that has its 3 axes aligned with the x, y, and z axes of the molecular data (atom coordinates) with padding 3 times the requested map resolution. This is described in the molmap documentation
http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html
The origin of the grid equals the x,y,z coordinates of the grid point at grid index 0,0,0. That origin position (floating point values) is saved in an MRC 2000 file in the
xorigin, yorigin, zorigin header fields.
http://www2.mrc-lmb.cam.ac.uk/image2000.html
It cannot be saved in the nxstart,nystart,nzstart header fields since those are integer values and can't accomodate the exact floating point origin. Here's the Chimera Python source code that reads MRC files if you want more details of how that is handled.
http://plato.cgl.ucsf.edu/trac/chimera/browser/trunk/libs/VolumeData/mrc/mrc_format.py
Tom
On Jul 8, 2013, at 1:08 PM, Catalin Buiu wrote:
Hello all,
I have the following problem: when I use molmap to generate
a map from an atomic structure, the result is fine: both align.
When I use another (Matlab) program to generate the map, the two,
when opened in Chimera do not align anymore.
My question is how molmap generates an origin index such that
both the map and the structure do align and how Chimera reads (and converts)
the xyz origin data from the MRC file?
Thank you.
Catalin Buiu
University of Bucharest, Romania
Chimera-users mailing list
Chimera-users@cgl.ucsf.edu
http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users