ana2ecat.c File Reference

Convert Analyze 7.5 image format to ECAT 6.3 or 7. More...

Go to the source code of this file.

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

Detailed Description

Convert Analyze 7.5 image format to ECAT 6.3 or 7.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file ana2ecat.c.

Function Documentation

imgAnalyzeToEcat()

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

Read Analyze database and write ECAT 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	Analyze database name
ecatfile	ECAT image filename, output extension determines the format
flipping	Flipping in Z dimension (usually 1=yes)
scanner_type	Scanner type, 0 if not set
zoom	recon_zoom, 0.0 or 1.0, if not known
verbose	Verbose level; if <=0 then nothing is printed in stdout

Definition at line 355 of file ana2ecat.c.

  {
  IMG img;
  int zi, fi, ret, file_format=0;
 
  if(verbose>0) printf("\nimgAnalyzeToEcat(%s, %s)\n", dbname, ecatfile);
  /* Check the arguments */
  if(dbname==NULL || ecatfile==NULL) return STATUS_FAULT;
 
  /* Read the first frame */
  imgInit(&img); fi=0;
  ret=imgReadAnalyzeFirstFrame(dbname, &img);
  if(ret!=STATUS_OK) {imgEmpty(&img); return ret;}
  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>10) imgInfo(&img);
    if(flipping!=anaFlipping()) imgFlipPlanes(&img);
    /* Write the frame in ECAT format */
    img._fileFormat=0;
    ret=imgWriteFrame(ecatfile, fi, &img, 0); //printf("ret := %d\n", ret);
    if(ret!=STATUS_OK) break;
    if(verbose>2) printf("    frame written.\n");
    /* Try to read the next frame */
    img._fileFormat=file_format;
    ret=imgReadAnalyzeFrame(dbname, fi+1, &img, 0);
  } 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;
}