IMPEM_BEGIN_INTERNAL_NAMESPACE
void MRCHeader::FromDensityHeader(const DensityHeader &h)
{
std::string empty;
nz=h.get_nz(); ny=h.get_ny(); nx=h.get_nx(); // map size
// mode
if(h.get_data_type()==0) // data type not initialized
mode = 2;
if(h.get_data_type()==1)
mode=0; // 8-bits
else if(h.get_data_type()==2)
mode=1;// 16-bits
else if(h.get_data_type()==5)
mode=2;// 32-bits
// number of first columns in map (default = 0)
nxstart=h.nxstart ; nystart=h.nystart ; nzstart=h.nzstart;
mx=h.mx ; my=h.my ; mz=h.mz; // Number of intervals along each dimension
xlen=h.xlen ; ylen=h.ylen ; zlen=h.zlen; // Cell dimensions (angstroms)
alpha=h.alpha ; beta=h.beta ; gamma=h.gamma; // Cell angles (degrees)
/* Axes corresponding to columns (mapc), rows (mapr) and sections (maps)
(1,2,3 for x,y,z) */
mapc=h.mapc;
mapr=h.mapr;
maps=h.maps;
/* Minimum, maximum and mean density value */
dmin=h.dmin;
dmax=h.dmax;
dmean=h.dmean;
ispg=h.ispg; // Sapce group number 0 or 1 (default 0)
nsymbt=h.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++)
user[i]=h.user[i];
strcpy(map,"MAP\0"); // character string 'MAP ' to identify file type
// Origin used for transforms
xorigin=h.get_xorigin() ; yorigin=h.get_yorigin() ; zorigin=h.get_zorigin();
// machine stamp (0x11110000 bigendian, 0x44440000 little)
machinestamp=h.machinestamp;
rms=h.rms; // RMS deviation of map from mean density
nlabl=h.nlabl; // Number of labels being used
// Copy comments
for(int i=0;i<nlabl;i++)
strcpy(labels[i],h.comments[i]);
// Fill empty coments with null character
const char *c="\0";
empty.resize(IMP_MRC_LABEL_SIZE,*c);
for(int i=nlabl;i<IMP_MRC_NUM_LABELS;i++)
strcpy(labels[i],empty.c_str());
}
//! Function to transfer the (compatible) information content from DensityHeader
//! to ImageHeader
void DensityHeader_to_ImageHeader(const DensityHeader& dh,ImageHeader& h) {
// map size
h.set_number_of_slices(dh.get_nz());
h.set_number_of_rows(dh.get_ny());
h.set_number_of_columns(dh.get_nx());
// mode
switch (dh.get_data_type()) {
case 1:
h.set_image_type(ImageHeader::VOL_BYTE);
break;
case 2:
h.set_image_type(ImageHeader::VOL_INT);
break;
case 5:
h.set_image_type(ImageHeader::VOL_IMPEM);
break;
}
// Reversed image?
if(algebra::get_is_big_endian()) {
h.set_reversed(dh.lswap!=1);
} else {
h.set_reversed(dh.lswap==1);
}
// Origins and pixel size
h.set_origin(algebra::Vector3D(dh.get_xorigin(),
dh.get_yorigin(),
dh.get_zorigin()) );
h.set_object_pixel_size(dh.get_spacing());
// Statistical values. There is no field in DensityHeader to guarantee that
// they are computed or correct, so they are ignored
h.set_fImami(0.);
h.set_fFmax(0.);
h.set_fFmin(0.);
h.set_fAv(0.);
h.set_fSig(0.);
// These fields don't have an equivalent in DensityHeader
h.set_Scale(1.);
h.set_fAngle1(0.);
h.set_Weight(1.);
h.set_Flip(0.);
h.set_fNrec(0.);
h.set_fNlabel(0.);
h.set_fIhist(0.);
h.set_fLabrec(0.);
h.set_fIangle(0.);
h.set_fLabbyt(0.);
h.set_fLenbyt(0.);
h.set_fFlag(0.); // Ignore all the Euler angles for the Spider Image
// Set time and date
h.set_date();
h.set_time();
}
