nii2ecat.c File Reference

Convert NIfTI-1 PET image to ECAT 7 format. More...

Go to the source code of this file.

Functions
int	imgNiftiToEcat (char *dbname, char *ecatfile, int scanner_type, float zoom, int verbose)

Detailed Description

Convert NIfTI-1 PET image to ECAT 7 format.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file nii2ecat.c.

Function Documentation

imgNiftiToEcat()

int imgNiftiToEcat	(	char *	dbname,
		char *	ecatfile,
		int	scanner_type,
		float	zoom,
		int	verbose )

Read NIfTI-1 image database and write ECAT7 image frame-by-frame

Returns

Returns errstatus, which is STATUS_OK (0) when call was successful, and >0 in case of an error.

Parameters

dbname	NIfTI database name
ecatfile	ECAT image filename, output extension must be .v
scanner_type	Scanner type, 0 if not set
zoom	recon_zoom, 0.0 or 1.0, if not known
verbose	Verbose level; if zero, then nothing is printed to stderr or stdout

Definition at line 305 of file nii2ecat.c.

  {
  IMG img;
  int zi, fi, ret, file_format=0;
 
  if(verbose>0) printf("\nimgNiftiToEcat(%s, %s)\n", dbname, ecatfile);
  /* Check the arguments */
  if(dbname==NULL || ecatfile==NULL) return STATUS_FAULT;
 
  /* Read the first frame */
  imgInit(&img); fi=0;
  if(verbose>1) printf("  reading the first frame\n");
  ret=imgReadNiftiFirstFrame(dbname, &img, verbose-1);
  if(ret!=STATUS_OK) {imgEmpty(&img); return ret;}
  /* ... and then the rest of the frames */
  do {
    fi++;
    if(verbose>1) printf("  frame %d\n", fi);
    /* Set header information */
    if(fi==0)  // must be saved for reading next frame
      file_format=img._fileFormat;
    if(scanner_type>0) imgSetScanner(&img, scanner_type);
    if(zoom>0.0) img.zoom=zoom;
    for(zi=0; zi<img.dimz; zi++) img.planeNumber[zi]=zi+1;
    if(verbose>15) imgInfo(&img);
    /* Write the frame in ECAT 7 format */
    img._fileFormat=IMG_E7;
    ret=imgWriteFrame(ecatfile, fi, &img, 0); //printf("ret := %d\n", ret);
    if(ret!=STATUS_OK) break;
    if(verbose>1) printf("  frame written.\n");
    /* Try to read the next frame */
    img._fileFormat=file_format;
    ret=imgReadNiftiFrame(dbname, fi+1, &img, 0, verbose-1);
  } while(ret==0);
  imgEmpty(&img);
  if(verbose>0 && ret==STATUS_NOMATRIX) {
    fprintf(stdout, "  %d frame(s) processed.\n", fi);
  }
  if(ret!=STATUS_OK && ret!=STATUS_NOMATRIX) {
    remove(ecatfile); return ret;
  }
  /*if(verbose>2)
    fprintf(stdout, "  Image saved in ECAT format in %s\n", ecatfile);*/
  return STATUS_OK;
}