misc.c File Reference

Misc functions processing both IMG and TAC structures. More...

#include "tpcclibConfig.h"
#include "tpctacimg.h"

Go to the source code of this file.

Functions
int	tacAllocateWithIMG (TAC *tac, IMG *img, int tacNr)
	Allocate TAC based on data in IMG.
int	imgVb (IMG *img, TAC *btac, double Vb, const int simVb, const int petVolume)
	Correct pixel TACs for vascular blood, or simulate its effect.

Detailed Description

Misc functions processing both IMG and TAC structures.

Definition in file misc.c.

Function Documentation

imgVb()

int imgVb	(	IMG *	img,
		TAC *	btac,
		double	Vb,
		const int	simVb,
		const int	petVolume )

Correct pixel TACs for vascular blood, or simulate its effect.

See also

tacVb

Returns

enum tpcerror (TPCERROR_OK when successful).

Parameters

img	Pointer to IMG data to process.
btac	Pointer to BTAC data to subtract or add; must contain the same sample times as IMG and y values must be in the same units.
Vb	Vb fraction [0,1].
simVb	Switch to either subtract vascular volume (0) or to simulate it (1).
petVolume	Switch to model vascular volume as either 0 : Cpet = (1-Vb)*Ct + Vb*Cb, or 1 : Cpet = Ct + Vb*Cb

Definition at line 56 of file misc.c.

  {
  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);
}

tacAllocateWithIMG()

int tacAllocateWithIMG	(	TAC *	tac,
		IMG *	img,
		int	tacNr )

Allocate TAC based on data in IMG.

See also

tacAllocate, tacAllocateMore, tacDuplicate, tacAllocateWithPAR

Returns

enum tpcerror (TPCERROR_OK when successful).

Author

Vesa Oikonen

Parameters

tac	Pointer to initiated TAC structure; any old contents are deleted. TAC sampleNr and tacNr will be set, although contents will be empty.

Precondition

Structure must be initiated with tacInit.

Parameters

img	Pointer to IMG structure which contains frame times and units.
tacNr	Number of TACs to allocate.

Definition at line 16 of file misc.c.

  {
  /* 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);
}