ird.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpcimgio.h"
/*****************************************************************************/
 
/*****************************************************************************/
int string_to_xyzf(
  const char *str,
  IMG_PIXEL *v
) {
  if(v==NULL) return 1;
  v->x=v->y=v->z=v->f=0;
 
  char *cptr, tmp[256];
  strncpy(tmp, str, 255); tmp[255]=(char)0;
  cptr=strtok(tmp, " ,;:()|-"); if(cptr==NULL) return 1;
  v->x=atoi(cptr); if(v->x<1) return 1;  
  cptr=strtok(NULL, " ,;:()|-"); if(cptr==NULL) return 2;
  v->y=atoi(cptr); if(v->y<1) return 1;
  cptr=strtok(NULL, " ,;:()|-"); if(cptr==NULL) return 3;
  v->z=atoi(cptr); if(v->z<1) return 1;
  cptr=strtok(NULL, " ,;:()|-"); if(cptr==NULL) return 0;
  v->f=atoi(cptr);
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int irdReorder(
  IMG_RANGE *img_range
) {
  int i;
 
  /* Check that input is ok */
  if(img_range==NULL) return 1;
  /* Change the order if necessary */
  if(img_range->x1<0 || img_range->x2<0) return 2;
  if(img_range->x2<img_range->x1) {
    i=img_range->x1; img_range->x1=img_range->x2; img_range->x2=i;}
  if(img_range->y1<0 || img_range->y2<0) return 3;
  if(img_range->y2<img_range->y1) {
    i=img_range->y1; img_range->y1=img_range->y2; img_range->y2=i;}
  if(img_range->z1<0 || img_range->z2<0) return 4;
  if(img_range->z2<img_range->z1) {
    i=img_range->z1; img_range->z1=img_range->z2; img_range->z2=i;}
  if(img_range->f1<0 || img_range->f2<0) return 5;
  if(img_range->f2<img_range->f1) {
    i=img_range->f1; img_range->f1=img_range->f2; img_range->f2=i;}
 
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int irdRead(
  char *irdfile,
  IMG_RANGE *img_range,     
  char *status
) {
  int ret, ii;
  IFT ift;
  char key[256];
  IMG_PIXEL v;
 
  /* Check that input is ok */
  if(irdfile==NULL || strnlen(irdfile, 2)<1 || img_range==NULL) {
    if(status!=NULL) strcpy(status, "program error");
    return 1;
  }
  /* Read IDF as IFT file */
  iftInit(&ift); ret=iftRead(&ift, irdfile, 1, 0); if(ret) {
    if(status!=NULL) strcpy(status, ift.status);
    iftEmpty(&ift); return 2;
  }
  /* Try to find keys 'corner1' and 'corner2' */
  strcpy(key, "corner1"); ii=iftGet(&ift, key, 0);
  if(ii>=0) {
    ret=string_to_xyzf(ift.item[ii].value, &v);
    if(ret==0) {
      img_range->x1=v.x; img_range->y1=v.y; img_range->z1=v.z; img_range->f1=v.f;
      strcpy(key, "corner2"); ii=iftGet(&ift, key, 0);
      if(ii>=0) {
        ret=string_to_xyzf(ift.item[ii].value, &v);
        img_range->x2=v.x; img_range->y2=v.y; img_range->z2=v.z; img_range->f2=v.f;
        if(ret==0) {
          irdReorder(img_range);
          if(status!=NULL) strcpy(status, "ok");
          iftEmpty(&ift); return 0;
        }
      }
    }
  }
  /* We are here only if keys were not found */
  /* Lets not care about keys at all */
  for(ii=0, ret=0; ii<ift.keyNr; ii++) {
    if(ret==0 && string_to_xyzf(ift.item[ii].value, &v)==0)
    {
      img_range->x1=v.x; img_range->y1=v.y; img_range->z1=v.z; img_range->f1=v.f;
      ret++; continue;
    }
    if(ret==1 && string_to_xyzf(ift.item[ii].value, &v)==0)
    {
      img_range->x2=v.x; img_range->y2=v.y; img_range->z2=v.z; img_range->f2=v.f;
      ret++; break;
    }
  }
  if(ret<2) {
    if(status!=NULL) strcpy(status, "volume definitions not found");
    iftEmpty(&ift); return 2;
  }
 
  irdReorder(img_range);
  if(status!=NULL) strcpy(status, "ok");
  iftEmpty(&ift);
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int irdCheck(
  IMG_RANGE *r,
  IMG *img
) {
  if(r==NULL || img==NULL) return(1);
  if(img->dimx<1 || img->dimy<1 || img->dimz<1 || img->dimt<1) return(2);
 
  if(r->x1<1 || r->x1>img->dimx) return(11);
  if(r->x2<1 || r->x2>img->dimx) return(12);
 
  if(r->y1<1 || r->y1>img->dimy) return(21);
  if(r->y2<1 || r->y2>img->dimy) return(22);
 
  if(r->z1<1 || r->z1>img->dimz) return(31);
  if(r->z2<1 || r->z2>img->dimz) return(32);
 
  /* If time frame range is not set, then set it now */
  if(r->f1<1 && r->f2<1) {
    r->f1=1; r->f2=img->dimt;
    return(0);
  }
  /* else check as usual */
  if(r->f1<1 || r->f1>img->dimt) return(41);
  if(r->f2<1 || r->f2>img->dimt) return(42);
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/