Data handling for projection matrix. More...

#include "libtpcrec.h"

Functions
void	prmatInit (PRMAT *mat)
void	prmatEmpty (PRMAT *mat)
int	prmatAllocate (PRMAT mat, int set, unsigned int rows, unsigned int coords)
unsigned int	prmatGetNV (PRMAT *mat)
unsigned int	prmatGetNB (PRMAT *mat)
unsigned int	prmatGetID (PRMAT *mat)
unsigned int	prmatGetPIX (PRMAT *mat)
unsigned int	prmatGetPixCoord (PRMAT *mat, int row)
unsigned int	prmatGetRays (PRMAT *mat, int row)
unsigned int	prmatGetBinView (PRMAT *mat, int row, int ind)
unsigned int	prmatGetRows (PRMAT *mat)
unsigned int	prmatGetPixels (PRMAT *mat, int row)
unsigned int	prmatGetCoord (PRMAT *mat, int row, int pix)
unsigned int	prmatGetXCoord (PRMAT *mat, int row, int pix)
unsigned int	prmatGetYCoord (PRMAT *mat, int row, int pix)
float	prmatGetFactor (PRMAT *mat, int row, int pix)
float	prmatGetMajor (PRMAT *mat)
float	prmatGetMinor (PRMAT *mat)
float	prmatGetMin (PRMAT *mat)
float	prmatGetMax (PRMAT *mat)
float	prmatGetFactorSum (PRMAT *mat)
float	prmatGetFactorSqrSum (PRMAT *mat, int row)
int	prmatSaveMatrix (PRMAT *mat)
int	prmatReadMatrix (char fname, PRMAT mat)

Variables
int	PRMAT_TEST
	If not 0 drive in test mode.
int	PRMAT_VERBOSE
	If not 0 drive in verbose mode.

Detailed Description

Data handling for projection matrix.

Author: Jarkko Johansson

Date: 2004-04-20

Definition in file prmat.c.

Function Documentation

◆ prmatAllocate()

int prmatAllocate	(	PRMAT *	mat,
		int	set,
		unsigned int	rows,
		unsigned int *	coords )

Allocates memory for PRMAT data. Normally used only in SET-functions.
Allocates memory for look-up table if set = 0 and for projection matrix otherwise.

Precondition: *mat is initialized && rows is positive && *coords list contains as many positive integers as there are rows.

Postcondition: memory is allocated for PRMAT structure mat.

Parameters

mat	pointer to projection matrix for which the allocation is done.
set	tells for which part of the structure allocation is to be done.
rows	number of rows in a matrix.
coords	list of number of entries in each row.

Returns: 0 if ok.

Definition at line 105 of file prmat.c.

{
  unsigned int *iptr, ir, ip;
  if(PRMAT_VERBOSE) printf("PRMAT:prmatAllocate(*mat,%d,*coords) \n",rows);
  
  /*check the arguments*/
  if(mat->status==PRMAT_STATUS_UNINITIALIZED) return(1);
  if(rows<1) return(2);
  
  //If set != 0 allocate for projection matrix.
  if(set){
    /*allocate memory*/
    mat->factor_sqr_sum=calloc(rows,sizeof(float));
    mat->dimr=rows; //for dimr
    mat->dime=calloc(rows,sizeof(int)); 
    mat->_factdata=
     (unsigned short int***)calloc(rows,sizeof(unsigned short int**)); 
    //if not enough memory
    if(mat->_factdata==NULL || mat->dime==NULL) return(4);
    //iterate through the rows
    for(ir=0,iptr=coords; ir<rows; ir++){
      //copy the coords vector in the same loop
      mat->dime[ir]=*iptr; 
      //*iptr is the number of (non-zero) factors in this row
      mat->_factdata[ir]=
       (unsigned short int**)calloc(*iptr,sizeof(unsigned short int*)); 
      //if not enough memory
      if(mat->_factdata[ir]==NULL) return(4); 
      //iterate through the factors in this row
      for(ip=0; ip<(*iptr); ip++){ 
  //NOTE that the leaves are coordinates and factor in this order
  mat->_factdata[ir][ip]=
          (unsigned short int*)calloc(3,sizeof(unsigned short int)); 
  if(mat->_factdata[ir][ip]==NULL) return(4);
      }
      iptr++;
    }
    /*set data pointers; these pointers are used to avoid data loss*/
    mat->fact=mat->_factdata;    
  } else {
    //If set was 0 allocate for look-up table.
    //allocate memory for look up table
    mat->nrp=rows; 
    mat->_linesdata=(unsigned int**)calloc(rows,sizeof(unsigned int*)); 
    //if not enough memory
    if(mat->_linesdata==NULL) return(4);
    //iterate through the rows
    for(ir=0,iptr=coords; ir<rows; ir++){
      //iptr is the number of lors belonging to this pixel
      mat->_linesdata[ir]=(unsigned int*)calloc(*iptr + 2,sizeof(unsigned int)); 
      //if not enough memory
      if(mat->_linesdata[ir]==NULL) return(4); 
      iptr++;
    }
    //set data pointers; these pointers are used to avoid data loss
    mat->lines=mat->_linesdata;    
  }
  return(0);
}

Referenced by prmatReadMatrix(), radonSetLORS(), and radonSetLUT().

◆ prmatEmpty()

void prmatEmpty ( PRMAT * mat )

Frees the memory allocated for mat. All data is cleared.

Parameters

mat	pointer to projection matrix to be emptied.

Definition at line 56 of file prmat.c.

{
  unsigned int ir, ip;
  if(PRMAT_VERBOSE) printf("PRMAT:prmatEmpty() \n");
  
  //if mat is allready empty
  if(mat->status==PRMAT_STATUS_INITIALIZED) return; 
  
  //if look-up table is occupied
  if(mat->status==PRMAT_STATUS_LU_OCCUPIED) { 
    /*free the occupied memory*/
    for(ir=0; ir<mat->nrp; ir++) { //every row
      free((unsigned int*)mat->_linesdata[ir]);
    }
    if(mat->nrp > 0) free(mat->_linesdata);
  }
 
  //if projection matrix is occupied
  if(mat->status==PRMAT_STATUS_PR_OCCUPIED || 
     mat->status==PRMAT_STATUS_BS_OCCUPIED )
  { 
    free((float*)mat->factor_sqr_sum);
    for(ir=0; ir<mat->dimr; ir++){ //every row
      for(ip=0; ip<mat->dime[ir]; ip++){ //every factor
  free((unsigned short int*)mat->_factdata[ir][ip]);
      }
    }
 
    free((int*)mat->dime);
  
    if(mat->dimr>0) free(mat->_factdata); 
    mat->dimr=0; 
  }
  return;
}

Referenced by prmatReadMatrix(), and radonSetLORS().

◆ prmatGetBinView()

unsigned int prmatGetBinView	(	PRMAT *	mat,
		int	row,
		int	ind )

Returns coordinates of a line of response in place (row,ind) in the look-up table.

Precondition: 0<=row<=prmatGetPIX(*mat) && 0<=ind<=prmatGetRays(*mat,row) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	the row index.
ind	the (non-zero) lor index.

Returns: int coordinates (view*binNr + bin) of the referred line of response, some negative value if ERROR.

Definition at line 253 of file prmat.c.

{
  return mat->lines[row][ind+2];
}

Referenced by prmatSaveMatrix().

◆ prmatGetCoord()

unsigned int prmatGetCoord	(	PRMAT *	mat,
		int	row,
		int	pix )

Returns coordinate of a pixel in place (row,pix) in the given projection matrix.

Precondition: 0<=row<=prmatGetLORS(*mat) && 0<=pix<=prmatGetPixels(row,*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	the row index.
pix	the (non-zero) pixel index.

Returns: coordinate of the referred pixel, some negative value if ERROR.

Definition at line 297 of file prmat.c.

{
  return (int)(mat->fact[row][pix][1])*mat->imgDim + (int)(mat->fact[row][pix][0]);
}

◆ prmatGetFactor()

float prmatGetFactor	(	PRMAT *	mat,
		int	row,
		int	pix )

Returns factor (weighting coefficient) for a pixel in place (row,pix) in given projection matrix.

Precondition: 0<=row<=prmatGetLORS(*mat) && 0<=pix<=prmatGetPixels(row,*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	the row index.
pix	the (non-zero) pixel index.

Returns: factor for the referred pixel, some negative value if ERROR.

Definition at line 342 of file prmat.c.

{
  return((float)(mat->scaling_factor)*(mat->fact[row][pix][2]));
}

Referenced by radonBackTransformPRM(), and radonFwdTransformPRM().

◆ prmatGetFactorSqrSum()

float prmatGetFactorSqrSum	(	PRMAT *	mat,
		int	row )

Returns square sum of all factors in given projection matrix in given row.

Precondition: mat is initialized && 0<=row<=prmatGetLORS(*mat).

Parameters

mat	pointer to a projection matrix.
row	row index for accessing factors.

Returns: square sum of all factors in given projection matrix in given row.

Definition at line 409 of file prmat.c.

{
  return(mat->factor_sqr_sum[row]);
}

Referenced by prmatSaveMatrix(), and radonSetLORS().

◆ prmatGetFactorSum()

float prmatGetFactorSum ( PRMAT * mat )

Returns sum of all factors in given projection matrix. For initializing image matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: sum of all factors in given projection matrix.

Definition at line 397 of file prmat.c.

{
  return(mat->factor_sum);
}

◆ prmatGetID()

unsigned int prmatGetID ( PRMAT * mat )

Returns the image dimension in the geometrics.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: the image dimension in the geometrics.

Definition at line 197 of file prmat.c.

{
  return mat->imgDim;
}

Referenced by radonBackTransformPRM(), and radonFwdTransformPRM().

◆ prmatGetMajor()

float prmatGetMajor ( PRMAT * mat )

Returns major semiaxe of the FOV in given projection matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: major semiaxe of the FOV in given projection matrix

Definition at line 353 of file prmat.c.

{
  return(mat->prmatfov[0]);
}

◆ prmatGetMax()

float prmatGetMax ( PRMAT * mat )

Returns greatest value in given projection matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: maximal value in given projection matrix.

Definition at line 386 of file prmat.c.

                              {
  return(mat->max);
}

◆ prmatGetMin()

float prmatGetMin ( PRMAT * mat )

Returns minimal value in given projection matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: minimal value in given projection matrix.

Definition at line 375 of file prmat.c.

{
  return(mat->min);
}

◆ prmatGetMinor()

float prmatGetMinor ( PRMAT * mat )

Returns minor semiaxe of the FOV in given projection matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: minor semiaxe of the FOV in given projection matrix

Definition at line 364 of file prmat.c.

{
  return(mat->prmatfov[1]);
}

◆ prmatGetNB()

unsigned int prmatGetNB ( PRMAT * mat )

Returns the number of bins (distances) in the geometrics.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: the number of bins in the geometrics.

Definition at line 186 of file prmat.c.

{
  return mat->binNr;
}

Referenced by radonBackTransformPRM(), and radonFwdTransformPRM().

◆ prmatGetNV()

unsigned int prmatGetNV ( PRMAT * mat )

Returns the number of views (angles) in the geometrics.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: the number of views in the geometrics.

Definition at line 175 of file prmat.c.

{
  return mat->viewNr;
}

Referenced by radonBackTransformPRM(), and radonFwdTransformPRM().

◆ prmatGetPIX()

unsigned int prmatGetPIX ( PRMAT * mat )

Returns number of rows (pixels) in the look-up table.

Precondition: mat is initialized.

Parameters

mat	pointer to a projection matrix.

Returns: number of pixels in look-up table, some negative value if ERROR.

Definition at line 210 of file prmat.c.

{
  return mat->nrp;
}

Referenced by prmatReadMatrix(), and prmatSaveMatrix().

◆ prmatGetPixCoord()

unsigned int prmatGetPixCoord	(	PRMAT *	mat,
		int	row )

Returns coordinates of a pixel in place 'row' in the look-up table.

Precondition: 0<=row<=prmatGetPIX(*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	index of the pixel whose coordinates are to be returned.

Returns: coordinates of a pixel in the look-up table, some negative value if ERROR.

Definition at line 223 of file prmat.c.

{
  return mat->lines[row][0];
}

Referenced by prmatSaveMatrix().

◆ prmatGetPixels()

unsigned int prmatGetPixels	(	PRMAT *	mat,
		int	row )

Returns number of pixels intersected by the given line of response.

Precondition: 0<=row<=prmatGetLORS(*mat) && mat is initialized.

Parameters

mat	pointer to the projection matrix.
row	index of the line of response for which the number of pixels is returned.

Returns: number of pixels intersected by given line of response, some negative value if ERROR.

Definition at line 282 of file prmat.c.

{
  return(mat->dime[row]);
}

Referenced by prmatReadMatrix(), prmatSaveMatrix(), radonBackTransformPRM(), radonFwdTransformPRM(), radonSetLORS(), and radonSetLUT().

◆ prmatGetRays()

unsigned int prmatGetRays	(	PRMAT *	mat,
		int	row )

Returns the number of lines of response (rays) intersecting the pixel in place 'row'.

Precondition: 0<=row<=prmatGetPIX(*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	index of the pixel for which the number of lors is returned.

Returns: number of lines of response intersecting the pixel, some negative value if ERROR.

Definition at line 237 of file prmat.c.

{
  return mat->lines[row][1];
}

Referenced by prmatReadMatrix(), and prmatSaveMatrix().

◆ prmatGetRows()

unsigned int prmatGetRows ( PRMAT * mat )

Returns number of rows (lines of response) in the given projection matrix.

Precondition: mat is initialized.

Parameters

mat	pointer to the projection matrix.

Returns: number of rows in the given projection matrix, some negative value if ERROR.

Definition at line 267 of file prmat.c.

{
  //NOTE no checking on parameters to accelerate iteration
  return(mat->dimr);
}

Referenced by prmatSaveMatrix(), and radonSetLORS().

◆ prmatGetXCoord()

unsigned int prmatGetXCoord	(	PRMAT *	mat,
		int	row,
		int	pix )

Returns the x-coordinate of a pixel in place (row,pix) in the given projection matrix.

Precondition: 0<=row<=prmatGetLORS(*mat) && 0<=pix<=prmatGetPixels(row,*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	the row index.
pix	the (non-zero) pixel index.

Returns: coordinate of the referred pixel, some negative value if ERROR.

Definition at line 312 of file prmat.c.

{
  return (int)(mat->fact[row][pix][0]);
}

Referenced by prmatSaveMatrix(), radonBackTransformPRM(), radonFwdTransformPRM(), radonSetLORS(), and radonSetLUT().

◆ prmatGetYCoord()

unsigned int prmatGetYCoord	(	PRMAT *	mat,
		int	row,
		int	pix )

Returns the y-coordinate of a pixel in place (row,pix) in the given projection matrix.

Precondition: 0<=row<=prmatGetLORS(*mat) && 0<=pix<=prmatGetPixels(row,*mat) && mat is initialized.

Parameters

mat	pointer to a projection matrix.
row	the row index.
pix	the (non-zero) pixel index.

Returns: coordinate of the referred pixel, some negative value if ERROR.

Definition at line 327 of file prmat.c.

{
  return (int)(mat->fact[row][pix][1]);
}

Referenced by prmatSaveMatrix(), radonBackTransformPRM(), radonFwdTransformPRM(), radonSetLORS(), and radonSetLUT().

◆ prmatInit()

void prmatInit ( PRMAT * mat )

Initializes PRMAT datatype for use. To be utilised before any use of PRMAT type variables.

Parameters

mat	pointer to projection matrix to be initialized.

Definition at line 22 of file prmat.c.

{
  if(PRMAT_VERBOSE) printf("PRMAT:prmatInit() \n"); 
  
  /*init the information of projection matrix*/
  //buffer mat to contain zero
  memset(mat, 0, sizeof(PRMAT)); 
  mat->type=PRMAT_TYPE_NA;
  mat->viewNr=0;
  mat->binNr=0;
  mat->mode=PRMAT_DMODE_01;
  mat->prmatfov= (int*)NULL;
  mat->max=mat->min=0.0;
  mat->scaling_factor=1.0;
  mat->factor_sum=0.0;
  mat->factor_sqr_sum=(float*)NULL;
 
  /*init look-up table*/
  mat->nrp=0;
  mat->lines=mat->_linesdata=(unsigned int**)NULL; //data array
 
  /*init projection data*/
  mat->dimr=0;
  mat->dime=(unsigned int*)NULL;
  mat->fact=mat->_factdata=(unsigned short int***)NULL; //data array
  
  mat->status=PRMAT_STATUS_INITIALIZED;
}

Referenced by radonSetLORS().

◆ prmatReadMatrix()

int prmatReadMatrix	(	char *	fname,
		PRMAT *	mat )

Reads one projection matrix from the given file to the given PRMAT structure mat.

Precondition: mat is initialized.

Postcondition: a projection matrix is read from the file fname.

Parameters

fname	file name to read.
mat	pointer to the projection matrix datatype.

Returns: 0 if OK.

Definition at line 550 of file prmat.c.

{
  unsigned int fov[2], *dimentry=NULL, *nrlor, *iptr, dimtmp, count=0;
  unsigned int row, lor=0, ind, fac;
  unsigned short int  *nrlor_c, *nrlor_l, *usiptr, *usiptr_c;
  unsigned short int x_coordinate=0, y_coordinate=0, value=0;
  float *sqrsum=NULL, *fptr, sqrtmp;
  FILE *fp;
 
  if(PRMAT_VERBOSE) printf("\n PRMAT:prmatReadMatrix(%s,mat) \n",fname);
  if(mat->status==PRMAT_STATUS_UNINITIALIZED) return(1);
  prmatEmpty(mat);
 
  fp=fopen(fname,"rb");
  if(fp==NULL) return(1);
 
  //read first the statical parts of the structure
  if(fread(&(mat->status),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: status: %i \n",mat->status);
  if(fread(&(mat->type),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: scanner type: %i \n",mat->type);
  if(fread(&(mat->viewNr),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(fread(&(mat->binNr),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(fread(&(mat->imgDim),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) 
    printf("PRMAT: scanner geometrics: (%i,%i,%i) \n",
           mat->viewNr,mat->binNr,mat->imgDim);
  if(fread(&(mat->mode),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: discretisation mode: %i \n",mat->mode);
  if(fread(fov,sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  if(fread((fov+1),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
  mat->prmatfov=calloc(2,sizeof(int)); 
  mat->prmatfov[0]=fov[0];
  mat->prmatfov[1]=fov[1];
  if(PRMAT_VERBOSE) 
    printf("PRMAT: FOV: (%i,%i) \n", mat->prmatfov[0],mat->prmatfov[1]);
  // Read min, max, sum and scaling factor.
  if(fread(&(mat->min),sizeof(float),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: min: %f \n",mat->min);
  if(fread(&(mat->max),sizeof(float),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: max: %f \n",mat->max);
  if(fread(&(mat->factor_sum),sizeof(float),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: sum of all factors: %f \n",mat->factor_sum);
  if(fread(&(mat->scaling_factor),sizeof(float),1,fp)==0) {fclose(fp); return(2);}
  if(PRMAT_VERBOSE) printf("PRMAT: scaling factor: %f \n",mat->scaling_factor);  
 
  // If projection matrix has been occupied.
  if(mat->status >= PRMAT_STATUS_BS_OCCUPIED){
    if(PRMAT_VERBOSE) printf("PRMAT: reading projection matrix \n");  
    if(fread(&(mat->dimr),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
    if(PRMAT_VERBOSE) 
      printf("PRMAT: number of rows in the projection matrix: %i \n",mat->dimr);  
    //read square sums
    sqrsum=calloc(mat->dimr,sizeof(float));  
    for(row=0, fptr=sqrsum;row<(mat->dimr); row++, fptr++) {
      if(fread(&sqrtmp,sizeof(float),1,fp)==0) {fclose(fp); return(2);}
      *fptr=sqrtmp;
    }
    //read dime
    dimentry=(unsigned int*)calloc(mat->dimr,sizeof(int));  
    count = 0;
    for(row=0,iptr=dimentry;row<(mat->dimr);row++,iptr++){
      if(fread(&dimtmp,sizeof(int),1,fp)==0) {fclose(fp); return(2);}
      *iptr=dimtmp;
      count += dimtmp;
    }
 
    //ready to allocate; NOTE that dime is also set in function prmatAllocate()    
    prmatAllocate(mat,1,mat->dimr,dimentry);
    if(PRMAT_VERBOSE)
      printf("PRMAT:prmatAllocate(mat) done in prmatReadMatrix(%s,mat)\n",fname);
 
    //now we can set square sums
    for(row=0, fptr=sqrsum;row<(mat->dimr);row++, fptr++){
      mat->factor_sqr_sum[row]=*fptr;
    }
    //read the actual data _factdata
    for(row=0;row<(mat->dimr);row++){
      //read coordinates and values in every row
      for(fac=0;fac<prmatGetPixels(mat,row);fac++){
  if(fread(&x_coordinate,sizeof(unsigned short int),1,fp)==0) {
          fclose(fp); return(2);}
  mat->fact[row][fac][0]=x_coordinate;
  if(fread(&y_coordinate,sizeof(unsigned short int),1,fp)==0) {
          fclose(fp); return(2);}
  mat->fact[row][fac][1]=y_coordinate;
  if(fread(&value,sizeof(unsigned short int),1,fp)==0) {
          fclose(fp); return(2);}
  mat->fact[row][fac][2]=value;
      }
    }
  }// END OF READING PROJECTIONS
 
  // If look-up table has been occupied.
  if(mat->status == PRMAT_STATUS_LU_OCCUPIED){
    if(PRMAT_VERBOSE) printf("PRMAT: reading look-up table \n");  
    if(fread(&(mat->nrp),sizeof(int),1,fp)==0) {fclose(fp); return(2);}
    if(PRMAT_VERBOSE)
      printf("PRMAT: number of pixels inside the FOV: %i \n",mat->nrp);
    //read nrlor
    nrlor_c=(unsigned short int*)calloc(mat->nrp,sizeof(unsigned short int));  
    nrlor_l=(unsigned short int*)calloc(mat->nrp,sizeof(unsigned short int));  
 
    //first coordinates of pixels.
    usiptr = nrlor_c;
    for(row=0; row<mat->nrp; row++, usiptr++){
      // Read next coordinate from file.
      if(fread(usiptr,sizeof(unsigned short int),1,fp)==0) {fclose(fp); return(2);}
    }
    //then number of lors hitting a pixel
    usiptr = nrlor_l;
    for(row=0; row<mat->nrp; row++, usiptr++){
      // Read next lornr from file.
      if(fread(usiptr,sizeof(unsigned short int),1,fp)==0) {fclose(fp); return(2);}
    }
 
    if(PRMAT_VERBOSE){
      // Get number of lors in lut.
      usiptr = nrlor_l;
      count = 0;
      for(row=0; row<mat->nrp; row++, usiptr++){
  count += *usiptr;
      }
      printf("PRMAT: nr of lors in lut:  %i \n",count);
      //printf("  nr of lors in ave:  %i \n",count/mat->nrp);
    }
 
    // Copy list usi format to integer list.
    nrlor = (unsigned int*)calloc(mat->nrp,sizeof(int));
    usiptr = nrlor_l;
    iptr = nrlor;
    for(row=0; row<mat->nrp; row++, iptr++){
      *iptr = (int)*usiptr++;
    }    
 
    // Now we can allocate memory for look-up table.
    prmatAllocate(mat, 0, mat->nrp, nrlor);
    if(PRMAT_VERBOSE)
      printf("PRMAT:prmatAllocate(mat) done in prmatReadMatrix(%s,mat)\n",fname);
    // And fill mat->lines.
    usiptr_c = nrlor_c;
    usiptr = nrlor_l;
    for(row=0;row<mat->nrp;row++ , usiptr++, usiptr_c++){
      mat->lines[row][0] = *usiptr_c;
      mat->lines[row][1] = *usiptr;
    }
 
    // Read the data _linesdata.
    for(row=0;row<prmatGetPIX(mat);row++){
      for(ind=0; ind<prmatGetRays(mat,row); ind++){
  if(fread(&lor,sizeof(unsigned int),1,fp)==0) {fclose(fp); return(2);}
  mat->lines[row][ind+2] = lor;
      }
    }
 
    free((unsigned short int*)nrlor_l);
    free((unsigned short int*)nrlor_c);
  }// END OF READING LOOK-UP TABLE
 
  if(PRMAT_VERBOSE) printf("PRMAT: prmatReadMatrix() done. \n");
  free((float*)sqrsum); free((int*)dimentry);
 
  fclose(fp);
 
  return(0);
}

◆ prmatSaveMatrix()

int prmatSaveMatrix ( PRMAT * mat )

Adds the given projection matrix to the file called prTDXXYY.prm. Where T is the type of the scanner (E=ECAT931 and G=GE Advance), D is the type of discretisation (0='0/1', 1='length of intersection' and 2='exact area'), XX is the major semiaxe of the FOV, and YY the minor semiaxe. If file already exists it is NOT created or rewritten.

Precondition: mat is initialized (contains a projection matrix).

Postcondition: mat is saved in the file prTDXXYY.prm.

Parameters

mat	pointer to the projection matrix to be saved.

Returns: 0 if OK.

Definition at line 428 of file prmat.c.

{
  char matrixfile[FILENAME_MAX], semi_x_c[2], semi_y_c[2], mod[2];
  int semi_x=mat->prmatfov[0], semi_y=mat->prmatfov[1], mode=mat->mode, itmp;
  unsigned short int x_coordinate, y_coordinate, value, p_coord, lors;
  unsigned int lor, row, fac, ind;
  float ftmp;
  FILE *fp;
  
  //Compose the file name to write
  sprintf(semi_x_c,"%i",semi_x);
  sprintf(semi_y_c,"%i",semi_y);
  //decide the type
  if(mat->type==PRMAT_TYPE_ECAT931) strcpy(matrixfile,"prE");
  else 
    if(mat->type==PRMAT_TYPE_GE) strcpy(matrixfile,"prG");
    else strcpy(matrixfile,"prN");
  //decide the mode 
  sprintf(mod,"%i",mode);
  strcat(matrixfile,mod);
  strcat(matrixfile,semi_x_c);
  strcat(matrixfile,semi_y_c);
  strcat(matrixfile,".prm");
  
  if(access(matrixfile, 0) != -1) {
    if(PRMAT_VERBOSE) printf("PRMAT:File %s already exists.\n", matrixfile);
    return(0);
  }
  
  fp=fopen(matrixfile,"wb+");
  if(fp==NULL) return(1);
 
  // Save status.
  fwrite(&(mat->status),sizeof(int),1,fp);  
 
  // Save type of the scanner.
  fwrite(&(mat->type),sizeof(int),1,fp);  
 
  // Save scanner geometrics.
  fwrite(&(mat->viewNr),sizeof(int),1,fp);  
  fwrite(&(mat->binNr),sizeof(int),1,fp);  
  fwrite(&(mat->imgDim),sizeof(int),1,fp);  
 
  // Save mode.
  fwrite(&(mat->mode),sizeof(int),1,fp);  
 
  // Save FOV.
  fwrite(&(mat->prmatfov[0]),sizeof(int),1,fp);
  fwrite(&(mat->prmatfov[1]),sizeof(int),1,fp); 
 
  // Save min, max, sum and scaling factor.
  fwrite(&(mat->min),sizeof(float),1,fp);
  fwrite(&(mat->max),sizeof(float),1,fp);
  fwrite(&(mat->factor_sum),sizeof(float),1,fp);
  fwrite(&(mat->scaling_factor),sizeof(float),1,fp);  
 
  // If projection matrix is set.
  if(mat->status >= PRMAT_STATUS_BS_OCCUPIED){
 
    fwrite(&(mat->dimr),sizeof(int),1,fp);
    //square sums
    for(row=0;row<prmatGetRows(mat);row++){
      ftmp=prmatGetFactorSqrSum(mat,row);
      fwrite(&ftmp,sizeof(float),1,fp);
    }
    //dime
    for(row=0;row<prmatGetRows(mat);row++){
      itmp=prmatGetPixels(mat,row);
      fwrite(&itmp,sizeof(int),1,fp);
    }
    //the actual data _factdata
    for(row=0;row<prmatGetRows(mat);row++){
      //coordinates and values in every row
      for(fac=0;fac<prmatGetPixels(mat,row);fac++){
  x_coordinate=prmatGetXCoord(mat,row,fac);
  fwrite(&x_coordinate,sizeof(unsigned short int),1,fp);
  y_coordinate=prmatGetYCoord(mat,row,fac);
  fwrite(&y_coordinate,sizeof(unsigned short int),1,fp);
  //NOTE that we wan't to save the values in unsigned short ints
  value=mat->fact[row][fac][2];
  fwrite(&value,sizeof(unsigned short int),1,fp);
      }
    }
  }// END OF SAVING PROJECTION MATRIX.
 
  // If look-up table is set.
  if(mat->status == PRMAT_STATUS_LU_OCCUPIED){
    //save first the statical parts of the structure
    fwrite(&(mat->nrp),sizeof(int),1,fp);
    //save lines: first coordinates of pixels inside the fov,
    for(row=0;row<prmatGetPIX(mat);row++){
      p_coord=prmatGetPixCoord(mat,row);
      fwrite(&p_coord,sizeof(unsigned short int),1,fp);
    }
    // then number of lors hitting a pixel.
    for(row=0; row<prmatGetPIX(mat); row++){
      lors=prmatGetRays(mat,row);
      fwrite(&lors,sizeof(unsigned short int),1,fp);
    }
 
    //save lors
    for(row=0;row<prmatGetPIX(mat);row++){
      for(ind=0; ind<prmatGetRays(mat,row); ind++){
  lor = prmatGetBinView(mat,row,ind);
  fwrite(&lor,sizeof(unsigned int),1,fp);
      }
    }
  }// END OF SAVING LOOK-UP TABLE.
 
  fclose(fp);
  return(0);
}

Variable Documentation

◆ PRMAT_TEST

int PRMAT_TEST

If not 0 drive in test mode.

Definition at line 8 of file prmat.c.

◆ PRMAT_VERBOSE

int PRMAT_VERBOSE

If not 0 drive in verbose mode.

Definition at line 9 of file prmat.c.

Referenced by prmatAllocate(), prmatEmpty(), prmatInit(), prmatReadMatrix(), and prmatSaveMatrix().

Functions

Variables

Detailed Description

Function Documentation

◆ prmatAllocate()

◆ prmatEmpty()

◆ prmatGetBinView()

◆ prmatGetCoord()

◆ prmatGetFactor()

◆ prmatGetFactorSqrSum()

◆ prmatGetFactorSum()

◆ prmatGetID()

◆ prmatGetMajor()

◆ prmatGetMax()

◆ prmatGetMin()

◆ prmatGetMinor()

◆ prmatGetNB()

◆ prmatGetNV()

◆ prmatGetPIX()

◆ prmatGetPixCoord()

◆ prmatGetPixels()

◆ prmatGetRays()

◆ prmatGetRows()

◆ prmatGetXCoord()

◆ prmatGetYCoord()

◆ prmatInit()

◆ prmatReadMatrix()

◆ prmatSaveMatrix()

Variable Documentation

◆ PRMAT_TEST

◆ PRMAT_VERBOSE