ecalibr.c File Reference

Calibrate reconstructed ECAT 931 image from counts to kBq/mL using specific file containing calibration coefficients. More...

Go to the source code of this file.

Functions
int	selectEcat931Calibrationfile (char *foldername, struct tm *stm, char *fname, int verbose)
int	readEcat931Calibrationfile (char *fname, double *coef, int verbose)

Detailed Description

Calibrate reconstructed ECAT 931 image from counts to kBq/mL using specific file containing calibration coefficients.

Updated cal4img version 1.3 (1997-09-14) which could only be compiled and run in Sun UNIX workstations.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file ecalibr.c.

Function Documentation

readEcat931Calibrationfile()

int readEcat931Calibrationfile	(	char *	fname,
		double *	coef,
		int	verbose )

Read ECAT 931 calibration coefficients from specific file.

Returns

0 if successful.

Parameters

fname	Filename which contains calibration factors for each image plane
coef	Pointer to double array of size 15, into which calibration coefficients are written.
verbose	Verbose level; if zero, then nothing is printed to stderr or stdout

Definition at line 453 of file ecalibr.c.

  {
  if(verbose>0) printf("readEcat931Calibrationfile('%s', ...)\n", fname);
  if(fname==NULL || coef==NULL) return(1);
 
  /* Read file into IFT struct */
  IFT ift; iftInit(&ift);
  if(iftRead(&ift, fname, 0, 0)!=0) return(2);
  if(ift.keyNr<15) {iftEmpty(&ift); return(3);}
 
  /* Jump over the first line which should contain the date, and then
     read plane factors from the next 15 lines */
  double v;
  for(int i=1; i<=15; i++) {
    if(verbose>1) printf("\t%d\t%s\n", i, ift.item[i].value);
    if(atof_with_check(ift.item[i].value, &v)) {iftEmpty(&ift); return(4);}
    if(v<=0.0) {iftEmpty(&ift); return(5);}
    coef[i-1]=v;
  }
  iftEmpty(&ift);
  return(0);
}

selectEcat931Calibrationfile()

int selectEcat931Calibrationfile	(	char *	foldername,
		struct tm *	stm,
		char *	fname,
		int	verbose )

Select correct calibration file from given folder.

Returns

0 if successful.

Parameters

foldername	Folder where calibration files are located
stm	Pointer to scan start time; calibration file is selected based on it.
fname	Pointer to allocated string where filename will be written
verbose	Verbose level; if zero, then nothing is printed to stderr or stdout

Definition at line 378 of file ecalibr.c.

  {
  if(verbose>0) printf("selectEcat931Calibrationfile('%s', ...)\n", foldername);
  if(foldername==NULL || stm==NULL || fname==NULL) return(1);
  strcpy(fname, "");
 
  char buf[256];
 
  /* Check whether foldername is a directory */
  struct stat fst;
  stat(foldername, &fst);
  if(!S_ISDIR(fst.st_mode)) { /* it is not */
    if(verbose>1) printf("  %s is not a directory\n", foldername);
    return(1);
  }
 
  /* Open the directory */
  DIR *dp;
  dp=opendir(foldername); if(dp==NULL) {
    if(verbose>1) printf("Error: cannot open directory %s\n", foldername);
    return(2);
  }
 
  /* Go throught the directory */
  struct dirent *de;
  char year[10], month[10], day[10], cfname[FILENAME_MAX];
  struct tm ftm;
  double tdif, tdif_min=-1.0E30;
  cfname[0]=(char)0;
  while((de=readdir(dp))!=NULL) {
    if(verbose>2) printf("d_name='%s'\n", de->d_name);
    if(de->d_name[0]=='.') continue; /* Ignore hidden and 'system' dirs */
    if(strlen(de->d_name)<6) continue;
    /* Identify the date from filename which should be in format YYMMDD */
    strlcpy(year, de->d_name, 3);
    strlcpy(month, de->d_name+2, 3);
    strlcpy(day, de->d_name+4, 3);
    sprintf(buf, "%s.%s.%s 16:00:00", day, month, year);
    /* Convert to struct tm */
    if(get_datetime(buf, &ftm, 0)) continue;
    if(verbose>3) printf("  %s\n", buf);
    /* Calculate time difference between calibration and scan time */
    tdif=tmDifference(&ftm, stm);
    if(verbose>3) printf("  sec_difference := %.0f\n", tdif);
    /* If later than scan time, then forget this */
    if(tdif>0.0) continue;
    /* If smaller difference than before, then save this */
    if(tdif>tdif_min) {tdif_min=tdif; strcpy(cfname, de->d_name);}
  }
  closedir(dp);
  if(tdif_min<-1.0E10 || !cfname[0]) return(3);
  if(verbose>1) printf("selected_file := %s\n", cfname);
 
  /* Combine path and name */
  if(foldername[strlen(foldername)-1]!='/')
    sprintf(fname, "%s/%s", foldername, cfname);
  else
    sprintf(fname, "%s%s", foldername, cfname);
 
  return(0);
}