ecatio.c

Go to the documentation of this file.

 
/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
/*****************************************************************************/
#include "tpcecat.h"
/*****************************************************************************/
 
/*****************************************************************************/
int ecatReadBlock(
  const char *filename,
  FILE *fp,
  const unsigned int blocknumber, 
  unsigned char *data
) {
  if(filename==NULL && fp==NULL) return(TPCERROR_FAIL);
  if(data==NULL) return(TPCERROR_FAIL);
 
  if(fp!=NULL) { // read from given file pointer
 
    /* If block number was given, then seek the block */
    if(blocknumber>0) {
      if(fseeko(fp, (off_t)(blocknumber-1)*ECATBLKSIZE, SEEK_SET) != 0) 
        return(TPCERROR_NO_DATA);
    }
    /* Read the block */
    if(fread(data, ECATBLKSIZE, 1, fp)<1) return(TPCERROR_CANNOT_READ);
 
  } else { // read from given file name
 
    /* Open file for read */
    FILE *localfp=fopen(filename, "rb");
    if(localfp==NULL) return(TPCERROR_CANNOT_OPEN);
 
    /* If block number was given, then seek the block */
    if(blocknumber>0) {
      if(fseeko(localfp, (off_t)(blocknumber-1)*ECATBLKSIZE, SEEK_SET) != 0)
        {fclose(localfp); return(TPCERROR_NO_DATA);}
    }
    /* Read the block */
    if(fread(data, ECATBLKSIZE, 1, localfp)<1) {fclose(localfp); return(TPCERROR_CANNOT_READ);}
    fclose(localfp);
 
  }
 
  /* That's it then */
  return(TPCERROR_OK);
}
/*****************************************************************************/
 
/*****************************************************************************/
int ecatVerifyMagic(
  const char *filename,
  FILE *fp
) {
 
  /* Read the first data block (512 bytes) */
  unsigned char buf[ECATBLKSIZE];
  int ret=ecatReadBlock(filename, fp, 0, buf);
  if(fp!=NULL) rewind(fp);
  if(ret!=TPCERROR_OK) return(0);
 
  /* Check if the start matches the ECAT 7 magic numbers */
  if(strncmp((char*)buf, "MATRIX72v", 9)==0 || strncmp((char*)buf, "MATRIX7011", 10)==0) return(7);
 
  /* ECAT 6 does not have defined magic numbers, but check anyway */
  if(strncmp((char*)buf, "MATRIX6", 7)==0 || strncmp((char*)buf, "ECAT6", 5)==0) return(6);
 
  /* Check that we do NOT have DICOM magic number */
  if(strncmp((char*)buf+128, "DICM", 4)==0) return(0);
 
  /* Check that we do NOT have NIfTI magic number */
  if(strncmp((char*)buf+344, "ni1", 3)==0 || strncmp((char*)buf, "n+1", 3)==0) return(0);
 
  /* Check that we do NOT have Interfile magic number */
  if(strncmp((char*)buf, "!INTERFILE", 10)==0) return(0);
 
  /* There is no indication that this IS NOT an ECAT 6 file */
  return(6);
}
/*****************************************************************************/
 
/*****************************************************************************/