void MRCHeader::ToDensityHeader(DensityHeader &h)
{
std::string empty;
h.update_map_dimensions(nx,ny,nz);
h.update_cell_dimensions();
// mode
if(mode==0)
h.set_data_type(1);
else if(mode==1)
h.set_data_type(2);
else if(mode==2)
h.set_data_type(5);
// number of first columns in map (default = 0)
h.nxstart=nxstart ; h.nystart=nystart ; h.nzstart=nzstart;
h.mx=mx ; h.my=my ; h.mz=mz; // Number of intervals along each dimension
h.xlen=xlen;h.ylen=ylen;h.zlen=zlen;//Cell dimentions(amgstroms)
h.alpha=alpha ; h.beta=beta ; h.gamma=gamma; // Cell angles (degrees)
/* Axes corresponding to columns (mapc), rows (mapr) and sections (maps)
(1,2,3 for x,y,z) */
h.mapc=mapc;
h.mapr=mapr;
h.maps=maps;
/* Minimum, maximum and mean density value */
h.dmin=dmin;
h.dmax=dmax;
h.dmean=dmean;
h.ispg=ispg; // Sapce group number 0 or 1 (default 0)
h.nsymbt=nsymbt; // Number of bytes used for symmetry data (0 or 80)
// extra space used for anything - 0 by default
for(int i=0;i<IMP_MRC_USER;i++)
h.user[i]=user[i];
strcpy(h.map,"MAP\0"); // character string 'MAP ' to identify file type
// Origin used for transforms
h.set_xorigin(xorigin) ; h.set_yorigin(yorigin) ; h.set_zorigin(zorigin);
// machine stamp (0x11110000 bigendian, 0x44440000 little)
h.machinestamp=machinestamp;
h.rms=rms; // RMS deviation of map from mean density
h.nlabl=nlabl; // Number of labels being used
// Copy comments
for(int i=0;i<h.nlabl;i++) {
//to make sure there is not memory leak
std::string temp;
temp.copy(labels[i],DensityHeader::COMMENT_FIELD_SINGLE_SIZE,0);
strcpy(h.comments[i],temp.c_str());
}
// Fill empty coments with null character
const char *c="\0";
empty.resize(IMP_MRC_LABEL_SIZE,*c);
for(int i=h.nlabl;i<IMP_MRC_NUM_LABELS;i++)
strcpy(h.comments[i],empty.c_str());
}
