x.c File Reference

Process time frames in image and TACs. More...

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

Go to the source code of this file.

Functions
int	tacimgXMatch (TAC *tac, IMG *img)
int	tacimgXCopy (TAC *tac, IMG *img)

Detailed Description

Process time frames in image and TACs.

Definition in file x.c.

Function Documentation

tacimgXCopy()

int tacimgXCopy	(	TAC *	tac,
		IMG *	img )

Copy x values (time frames) from IMG structure into TAC structure. If IMG has more frames than has been allocated for TAC (_sampleNr), then additional frames are ignored. If IMG has less frames than TAC, then error is returned.

Returns

Returns tpcerror; 0 if successful and >0 in case of an error.

See also

tacXCopy, tacIsX, tacSetX, tacimgXMatch, imgXRange

Parameters

tac	Pointer to TAC structure.
img	Pointer to IMG structure.

Definition at line 74 of file x.c.

  {
  if(img==NULL || tac==NULL) return(TPCERROR_NO_DATA);
  if(tac->sampleNr>img->dimt || img->dimt<1) return(TPCERROR_INCOMPATIBLE_DATA);
 
  int sNr=tac->_sampleNr; if(img->dimt<sNr) sNr=img->dimt;
 
  for(int i=0; i<sNr; i++) {
    tac->x1[i]=img->x1[i];
    tac->x2[i]=img->x2[i];
    tac->x[i]=0.5*(tac->x1[i]+tac->x2[i]);
  }
  tac->isframe=1;
  tac->tunit=img->tunit;
  return(TPCERROR_OK); 
}

tacimgXMatch()

int tacimgXMatch	(	TAC *	tac,
		IMG *	img )

Check whether sample (frame) times are the same (or very close to) in TAC and IMG structure structures. If number of samples/frames is different, then only common number of samples are compared.

Returns

Returns 0 in case of no match, 1 if times do match, and <1 in case of an error.

See also

imgCompareTimes, tacXMatch, tacXUnitConvert, tacimgXCopy

Parameters

tac	Pointer to TAC structure.
img	Pointer to IMG structure.

Definition at line 17 of file x.c.

  {
  if(img==NULL || tac==NULL) return(-1);
 
  double acceptedDif=2.2; // sec
 
  /* Get the last common sample nr */
  int sNr=tac->sampleNr; if(img->dimt<sNr) sNr=img->dimt;
  if(sNr==0) return(0);
 
  /* Get the x of the last common sample */
  double tLastImg=0.5*(img->x1[sNr-1]+img->x2[sNr-1]);
  double tLastTac=tac->x[sNr-1]; if(tac->isframe) tLastTac=0.5*(tac->x1[sNr-1]+tac->x2[sNr-1]);
  if(acceptedDif>tLastImg) acceptedDif=0.1*tLastImg;
 
  /* Need to convert x unit? */
  double tscale=nan("");
  if(tac->tunit==UNIT_SEC && img->tunit==UNIT_SEC) tscale=1.;
  else if(tac->tunit==UNIT_SEC && img->tunit==UNIT_MIN) tscale=1.0/60;
  else if(tac->tunit==UNIT_MIN && img->tunit==UNIT_SEC) tscale=60.;
  else if(tac->tunit==UNIT_MIN && img->tunit==UNIT_MIN) tscale=1.;
  else {
    if(tLastImg>30.*tLastTac) tscale=60.;
    else if(fabs(tLastImg-tLastTac)<30.) tscale=1.;
  }
  if(!isfinite(tscale)) return(0);
 
  /* Compare sample times frame-by-frame */
  int n=0;
  if(tac->isframe==0) { // check frame mid times
    for(int i=0; i<sNr; i++) {
      double f=fabs(tac->x[i]*tscale - 0.5*(img->x1[i]+img->x2[i]));
      if(f>acceptedDif) n++;
    }
  } else {
    for(int i=0; i<sNr; i++) {
      double f=fabs(tac->x1[i]*tscale - img->x1[i]);
      if(f>acceptedDif) {n++; continue;}
      f=fabs(tac->x2[i]*tscale - img->x2[i]);
      if(f>acceptedDif) n++;
    }
  }
  if(n==0) return(1); else return(0); 
}