imgcbv.c File Reference

Subtracts or simulates the contribution of vascular radioactivity to dynamic PET image. More...

Go to the source code of this file.

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

Detailed Description

Subtracts or simulates the contribution of vascular radioactivity to dynamic PET image.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file imgcbv.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 76 of file imgcbv.c.

  {
  int ii, mi, ret;
  char tmp[1024];
  ECAT_HEADERS ehdr;
 
  if(verbose>1) printf("ecatCopyHeaders(%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;
}