int iiMRCLoadPCoord(ImodImageFile *inFile, IloadInfo *li, int nx, int ny,
int nz)
{
int i;
b3dUInt16 pcoordxy[2];
b3dInt16 pcoordz;
int offset=1024;
int nread = nz;
MrcHeader *hdr = (MrcHeader *)inFile->header;
int iflag = hdr->nreal;
int nbytes = hdr->nint;
int nextra = hdr->next;
if (!iiMRCcheckPCoord(hdr))
return 0;
if (iflag & TILT_FLAG)
offset += 2;
if (nbytes * nz > nextra) {
nread = nextra / nbytes;
b3dError(stderr, "There are piece coordinates for only %d frames in"
" the extra header\n", nread);
}
li->pcoords = (int *)malloc(sizeof(int) * 3 * nz);
fseek(inFile->fp, offset, SEEK_SET);
for (i = 0; i < nread; i++) {
fread(pcoordxy, 2, 2, inFile->fp);
fread(&pcoordz, 1, 2, inFile->fp);
/* add swapping 10/2/00 */
if (hdr->swapped) {
mrc_swap_shorts((b3dInt16 *)pcoordxy, 2);
mrc_swap_shorts(&pcoordz, 1);
}
if (ferror(inFile->fp)) {
nread = i;
b3dError(stderr, "Error reading piece coordinates from extra"
" header after %d frames\n", i);
break;
}
li->pcoords[(i*3)] = pcoordxy[0];
li->pcoords[(i*3)+1] = pcoordxy[1];
li->pcoords[(i*3)+2] = pcoordz;
offset = nbytes - 6;
if (offset > 0)
fseek(inFile->fp, offset, SEEK_CUR);
}
li->plist = nread;
return(mrc_plist_proc(li, nx, ny, nz));
}
int iiMRCLoadPCoord(ImodImageFile *inFile, struct LoadInfo *li, int nx, int ny,
int nz)
{
int i;
b3dUInt16 pcoordxy[2];
b3dInt16 pcoordz;
int extra_bytes[32];
int extratot = 0;
int offset=1024;
int nread = nz;
MrcHeader *hdr = (MrcHeader *)inFile->header;
int iflag = hdr->nreal;
int nbytes = hdr->nint;
int nextra = hdr->next;
int flag_count;
if(!nextra || !(iflag & MONTAGE_FLAG))
return 0;
b3dHeaderItemBytes(&flag_count, &extra_bytes[0]);
/* DNM 12/10/01: as partial protection against mistaking other entries
for montage information, at least make sure that the total bytes
implied by the bits in the flag equals the nint entry. Also
reject deltavision files */
/* DNM 2/3/02: make sure nreal also does not have bits beyond the flags */
for (i = 0; i < flag_count; i++)
if (iflag & (1 << i))
extratot += extra_bytes[i];
if (extratot != nbytes || hdr->creatid == -16224 ||
iflag >= (1 << flag_count))
return 0;
if (iflag & TILT_FLAG)
offset += 2;
if (nbytes * nz > nextra) {
nread = nextra / nbytes;
b3dError(stderr, "There are piece coordinates for only %d frames in"
" the extra header\n", nread);
}
li->pcoords = (int *)malloc(sizeof(int) * 3 * nz);
fseek(inFile->fp, offset, SEEK_SET);
for (i = 0; i < nread; i++) {
fread(pcoordxy, 2, 2, inFile->fp);
fread(&pcoordz, 1, 2, inFile->fp);
/* add swapping 10/2/00 */
if (hdr->swapped) {
mrc_swap_shorts((b3dInt16 *)pcoordxy, 2);
mrc_swap_shorts(&pcoordz, 1);
}
if (ferror(inFile->fp)) {
nread = i;
b3dError(stderr, "Error reading piece coordinates from extra"
" header after %d frames\n", i);
break;
}
li->pcoords[(i*3)] = pcoordxy[0];
li->pcoords[(i*3)+1] = pcoordxy[1];
li->pcoords[(i*3)+2] = pcoordz;
offset = nbytes - 6;
if (offset > 0)
fseek(inFile->fp, offset, SEEK_CUR);
}
li->plist = nread;
return(mrc_plist_proc(li, nx, ny, nz));
}
