imgcalc.c File Reference

Simple arithmetic calculations on ECAT files. More...

Go to the source code of this file.

Functions
int	ecat7CopyHeadersNoQuant (const char *file1, const char *file2, int verbose)

Detailed Description

Simple arithmetic calculations on ECAT files.

Previous application name ecatcalc.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Test

Test other image formats than ECAT.

Definition in file imgcalc.c.

Function Documentation

ecat7CopyHeadersNoQuant()

int ecat7CopyHeadersNoQuant	(	const char *	file1,
		const char *	file2,
		int	verbose )

Copy ECAT7 header information from one file to another. Quantitative information is not copied. This function should be kept local to application imgcalc, but this may be useful in other programs, if copied/not copied fields are edited.

Returns

Returns 0 when successful, otherwise >0.

Parameters

file1	ECAT file FROM which header information is copied
file2	ECAT file IN which header information is copied to
verbose	Verbose level

Definition at line 86 of file imgcalc.c.

  {
  int ii, mi, ret;
  char tmp[1024];
  ECAT_HEADERS ehdr;
 
  if(verbose>1) printf("ecatCopyHeadersNoQuant(%s, %s)\n", file1, file2);  
  if(file1==NULL || file2==NULL) return(1);
  
  
  /* Read headers */
  ehdrInitiate(&ehdr);
  ret=ecat7ReadHeaders(file1, &ehdr, verbose-1); if(ret!=0) return(10);
  /* Remove main header information that should NOT be copied */
  strcpy(tmp, "user_process_code");
  ii=iftGet(&ehdr.mh, tmp, 0); iftDeleteItem(&ehdr.mh, ii, 0);
  strcpy(tmp, "ecat_calibration_factor");
  ii=iftGet(&ehdr.mh, tmp, 0); iftDeleteItem(&ehdr.mh, ii, 0);
 
  /* Remove subheader information that should NOT be copied */
  for(mi=0; mi<ehdr.nr; mi++) {
    strcpy(tmp, "scale_factor");
    ii=iftGet(&ehdr.m[mi].sh, tmp, 0); iftDeleteItem(&ehdr.m[mi].sh, ii, 0);
    strcpy(tmp, "image_min");
    ii=iftGet(&ehdr.m[mi].sh, tmp, 0); iftDeleteItem(&ehdr.m[mi].sh, ii, 0);
    strcpy(tmp, "image_max");
    ii=iftGet(&ehdr.m[mi].sh, tmp, 0); iftDeleteItem(&ehdr.m[mi].sh, ii, 0);
    strcpy(tmp, "scan_min");
    ii=iftGet(&ehdr.m[mi].sh, tmp, 0); iftDeleteItem(&ehdr.m[mi].sh, ii, 0);
    strcpy(tmp, "scan_max");
    ii=iftGet(&ehdr.m[mi].sh, tmp, 0); iftDeleteItem(&ehdr.m[mi].sh, ii, 0);
  }
  
  /* Write headers */
  ret=ecat7WriteHeaders(file2, &ehdr, verbose-1);
  ehdrEmpty(&ehdr);
  if(ret!=0) return(10+ret);
 
  return 0;
}