x.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include "tpctacimg.h"
/*****************************************************************************/
 
/*****************************************************************************/
int tacimgXMatch(
  TAC *tac,
  IMG *img
) {
  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); 
}
/*****************************************************************************/
 
/*****************************************************************************/
int tacimgXCopy(
  TAC *tac,
  IMG *img
) {
  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); 
}
/*****************************************************************************/
 
/*****************************************************************************/