misc.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include "tpctacimg.h"
/*****************************************************************************/
 
/*****************************************************************************/ 
int tacAllocateWithIMG(
  TAC *tac,
  IMG *img,
  int tacNr
) {
  /* Check the input */
  if(tac==NULL || img==NULL || tacNr<1) return(TPCERROR_FAIL);
  if(img->dimt<1) return(TPCERROR_NO_DATA);
 
  /* Allocate TAC */
  int ret=tacAllocate(tac, img->dimt, tacNr);
  if(ret!=TPCERROR_OK) return(ret);
  tac->sampleNr=img->dimt; tac->tacNr=tacNr;
 
  /* Copy frame times */
  tac->isframe=1;
  for(int i=0; i<img->dimt; i++) {
    tac->x1[i]=img->x1[i];
    tac->x2[i]=img->x2[i];
    tac->x[i]=0.5*(tac->x1[i]+tac->x2[i]);
  }
 
  /* Copy data units */
  tac->cunit=img->cunit;
  tac->tunit=img->tunit;
 
  return(TPCERROR_OK);
}
/*****************************************************************************/
 
/*****************************************************************************/
int imgVb(
  IMG *img,
  TAC *btac,
  double Vb,
  const int simVb,
  const int petVolume
) {
  if(img==NULL || btac==NULL) return(TPCERROR_FAIL);
  if(img->dimt<1 || btac->tacNr<1) return(TPCERROR_NO_DATA);
  if(img->dimt>btac->sampleNr) return(TPCERROR_INCOMPATIBLE_DATA);
  if(!(Vb>=0.0) || !(Vb<=1.0)) return TPCERROR_INVALID_VALUE;
  if(Vb==1.0 && petVolume==0 && simVb==0) // combination would lead to divide by zero
    return(TPCERROR_INVALID_VALUE);
 
  if(simVb==0) {
    for(int zi=0; zi<img->dimz; zi++)
      for(int yi=0; yi<img->dimy; yi++)
        for(int xi=0; xi<img->dimx; xi++)
          for(int ti=0; ti<img->dimt; ti++) {
            img->m[zi][yi][xi][ti]-=Vb*btac->c[0].y[ti];
            if(petVolume==0) img->m[zi][yi][xi][ti]/=(1.0-Vb);
          }
  } else {
    for(int zi=0; zi<img->dimz; zi++)
      for(int yi=0; yi<img->dimy; yi++)
        for(int xi=0; xi<img->dimx; xi++)
          for(int ti=0; ti<img->dimt; ti++) {
            if(petVolume==0) img->m[zi][yi][xi][ti]*=(1.0-Vb);
            img->m[zi][yi][xi][ti]+=Vb*btac->c[0].y[ti];
          }
  }
 
  return(TPCERROR_OK);
}
/*****************************************************************************/
 
/*****************************************************************************/