[Chimera-users] molmap and origin index
goddard at sonic.net
Mon Jul 15 10:56:13 PDT 2013
The Matlab calculation of the map is where the problem lies. I don't know what Matlab script you are using so it is hard for me to advise. I would not expect the Matlab script to produce the same voxel size, origin and map size as the Chimera molmap command. Do you tell Matlab the voxel size, origin and map size to calculate? If not, why would you expect them to be the same values as the Chimera molmap command? The molmap command adds padding around the molecule when computing the map and I wouldn't expect the Matlab script to add the same amount of padding. So the Matlab output wouldn't have the same origin grid index.
What you need to know is for each Matlab map grid index i,j,k what is the corresponding position x,y,z in the PDB atom coordinates used to generate the map. Then the matlab script needs to write the MRC xorigin, yorigin, zorigin as the atom x,y,z position for grid index 0,0,0. Alternatively if you cannot change the Matlab script, you can change the "Origin index" in the Chimera volume dialog Coordinates pane. This value should be the i,j,k grid index that corresponds to atom coordinates x,y,z = (0,0,0) and those grid indices can be non-integer.
On Jul 14, 2013, at 7:26 AM, Catalin Buiu wrote:
> Hi Tom,
> thanks for your reply.
> Now the both maps (molmap and Matlab) have the same origin index, the same voxel size and still the Matlab map is shifted away from the other map (and atomic structure). Can you tell me a possible reason for that?
> Thanks again,
> From: Tom Goddard <goddard at sonic.net>
> To: Catalin Buiu <cbuiu at yahoo.com>
> Cc: "chimera-users at cgl.ucsf.edu" <chimera-users at cgl.ucsf.edu>
> Sent: Monday, July 8, 2013 11:46 PM
> Subject: Re: [Chimera-users] molmap and origin index
> 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
> 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.
> 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.
> 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 at cgl.ucsf.edu
> Chimera-users mailing list
> Chimera-users at cgl.ucsf.edu
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Chimera-users