flat2img.c File Reference

Construct an ECAT PET file from binary flat file. More...

Go to the source code of this file.

Functions
void	zero_hot_spots (float *fdata, int dim)
int	readMatrixInformation (char *fname, int *planeNr, int *frameNr, int *dim1, int *dim2)

Detailed Description

Construct an ECAT PET file from binary flat file.

Application name was previously flo2ecat.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file flat2img.c.

Function Documentation

readMatrixInformation()

int readMatrixInformation	(	char *	fname,
		int *	planeNr,
		int *	frameNr,
		int *	dim1,
		int *	dim2 )

Read ASCII file containing information on matrix dimensions

Returns

Returns 0 when successful, otherwise <> 0.

Parameters

fname	Matrix information filename
planeNr	matrix plane number (zdim)
frameNr	matrix time frame number (tdim)
dim1	matrix column number (xdim)
dim2	matrix row number (ydim)

Definition at line 468 of file flat2img.c.

  {
  FILE *fp;
  int n;
  
  if(fname==NULL || planeNr==NULL || frameNr==NULL || dim1==NULL || dim2==NULL)
    return(1);
  fp=fopen(fname, "r"); if(fp==NULL) return(2);
  n=fscanf(fp, "%d %d %d %d", planeNr, frameNr, dim1, dim2);
  fclose(fp); if(n!=4) return(3);
  if(*planeNr<1 || *frameNr<1 || *dim1<1 || *dim2<1) return(4);
  return(0);
}

zero_hot_spots()

void zero_hot_spots	(	float *	fdata,
		int	dim )

Calculate the distance between the pixel and the center of matrix. If the distance is greater than radius of FOV, pixel value is set to zero.

Parameters

fdata	Pointer to array containing one 2D image matrix
dim	Dimensions (x and y) of the image matrix

Definition at line 441 of file flat2img.c.

  {
  float   *fptr;
  float    r, d;
  int      i;
 
  r=(float)dim/2.0; if(r>25.0) r-=5.0;
  fptr=fdata;
  for(i=0; i<dim*dim; i++, fptr++ ) {
    d=hypotf((float)(i/dim-dim/2), (float)(i-dim*(i/dim)-dim/2));
/*
    d=sqrt((i/dim-dim/2)*(i/dim-dim/2)+
           (i-dim*(i/dim)-dim/2)*(i-dim*(i/dim)-dim/2));
*/
    if(d>r) *fptr=0.0;
  }
}