idimask.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpcidi.h"
/*****************************************************************************/
 
/*****************************************************************************/
int imgMaskPixelTACs(
  IMG *img,
  IMG *mask,
  double thrs,
  DFT *dft,
  int verbose
) {
  if(verbose>0) printf("%s()\n", __func__);
 
  int xi, yi, zi, fi, ri, ret;
 
  if(img->status<IMG_STATUS_OCCUPIED) return(1);
  if(mask->status<IMG_STATUS_OCCUPIED) return(2);
  if(mask->dimz!=img->dimz) return(3);
  if(mask->dimy!=img->dimy) return(4);
  if(mask->dimx!=img->dimx) return(5);
  if(img->dimt<1 || mask->dimt<1) return(6);
  if(dft==NULL) return(7);
 
  /* Delete any previous TAC data; do it first so that caller can check
     the nr of extracted TACs inside DFT */
  dftEmpty(dft);
 
  /* Calculate the nr of voxels to extract */
  long long nr=0;
  for(zi=0; zi<mask->dimz; zi++)
    for(yi=0; yi<mask->dimy; yi++)
      for(xi=0; xi<mask->dimx; xi++)
        if(mask->m[zi][yi][xi][0]>thrs) nr++;
  if(verbose>1) printf("mask_pixel_nr := %lld\n", nr);
  /* Return, if no voxels were marked in mask image */
  if(nr==0) return(0);
 
  /* Allocate memory for pixel TACs */
  ret=dftAllocateWithIMG(dft, nr, img);
  if(ret!=0) {
    if(verbose>0) fprintf(stderr, "Error: cannot allocate memory for %lld pixel TACs.\n", nr);
    return(10);
  }
  strcpy(dft->studynr, img->studyNr);
 
  /* Save the pixel TACs in DFT */
  ri=0;
  for(zi=0; zi<mask->dimz; zi++)
    for(yi=0; yi<mask->dimy; yi++)
      for(xi=0; xi<mask->dimx; xi++) if(mask->m[zi][yi][xi][0]>thrs)
  {
    sprintf(dft->voi[ri].voiname, "%d", 1+xi);
    sprintf(dft->voi[ri].hemisphere, "%d", 1+yi);
    sprintf(dft->voi[ri].place, "%d", 1+zi);
    sprintf(dft->voi[ri].name, "%d %d %d", 1+xi, 1+yi, 1+zi);
    for(fi=0; fi<img->dimt; fi++) dft->voi[ri].y[fi]=img->m[zi][yi][xi][fi];
    dft->voi[ri].size=mask->m[zi][yi][xi][0];
    ri++;
  }
  dft->voiNr=ri;
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/