imgunits.c File Reference

Setting image calibration unit. More...

#include "libtpcimgio.h"

Go to the source code of this file.

Functions
int	imgUnitId (char *unit)
void	imgUnitFromEcat (IMG *img, int ecat_unit)
void	imgUnitFromEcat7 (IMG *img, ECAT7_mainheader *h)
int	imgUnitToEcat6 (IMG *img)
void	imgUnitToEcat7 (IMG *img, ECAT7_mainheader *h)
char *	imgUnit (int dunit)
int	imgSetUnit (IMG *img, char *unit)

Detailed Description

Setting image calibration unit.

Author

Vesa Oikonen

Definition in file imgunits.c.

Function Documentation

imgSetUnit()

int imgSetUnit	(	IMG *	img,
		char *	unit )

Sets the unit in IMG based on specified unit string.

See also

imgUnitId, imgUnitFromEcat, imgUnitFromEcat7

Returns

Returns 0 if successful.

Definition at line 328 of file imgunits.c.

{
  int new_unit;
 
  if(img==NULL || unit==NULL) return(1);
  new_unit=imgUnitId(unit); if(new_unit<0) return(1-new_unit);
  img->unit=new_unit;
  return(0);
}

Referenced by imgSetUnit().

imgUnit()

char * imgUnit

(

int

dunit

)

Return pointer to string describing the calibrated image data unit

Returns

Pointer to string

Parameters

dunit

Unit id number.

Definition at line 315 of file imgunits.c.

  {
  return(petCunit(dunit));
}

Referenced by clusterTACs(), cunit_check_dft_vs_img(), dftAllocateWithIMG(), imgInfo(), imgUnit(), and sif2dft().

imgUnitFromEcat()

void imgUnitFromEcat	(	IMG *	img,
		int	ecat_unit )

Set IMG calibration unit based on ECAT 6.3 unit.

See also

imgUnitFromEcat7, imgUnitId, imgUnitToEcat6, imgSetUnit

Parameters

img	Pointer to target IMG structure.
ecat_unit	Calibration unit code id.

Definition at line 95 of file imgunits.c.

  {
  switch(ecat_unit) {
    case  0: /* Unknown               */
      img->unit=CUNIT_UNKNOWN; break;
    case  1: /* MBq/mL                */
      img->unit=CUNIT_MBQ_PER_ML; break;
    case  2: /* ECAT counts           */
      img->unit=CUNIT_COUNTS; break;
    case  3: /* uCi/ml                */
      img->unit=CUNIT_UCI_PER_ML; break;
    case  4: /* LMRGlu                */
      img->unit=CUNIT_UNKNOWN; break;
    case  5: /* LMRUGlu umol/min/100g */
      img->unit=CUNIT_UMOL_PER_MIN_PER_100G; break;
    case  6: /* LMRUGlu mg/min/100g   */
      img->unit=CUNIT_MG_PER_MIN_PER_100G; break;
    case  7: /* nCi/mL                */
      img->unit=CUNIT_NCI_PER_ML; break;
    case  8: /* Well counts           */
      img->unit=CUNIT_CPS; break;
    case  9: /* Bq/mL                 */
      img->unit=CUNIT_BQ_PER_ML; break;
    case 10: /* kBq/mL                */
      img->unit=CUNIT_KBQ_PER_ML; break;
    case 11: /* 1/min                 */
      img->unit=CUNIT_PER_MIN; break;
    case 12: /* mL/min/100g           */
      img->unit=CUNIT_ML_PER_DL_PER_MIN; break;
    case 13: /* sec*kBq/mL            */
      img->unit=CUNIT_SEC_KBQ_PER_ML; break;
    case 14: /* sec*nCi/mL            */
      img->unit=CUNIT_UNKNOWN; break;
    case 15: /* 1/sec                 */
      img->unit=CUNIT_PER_SEC; break;
    case 16: /* Unitless              */
      img->unit=CUNIT_UNITLESS; break;
    case 17: /* Hounsfield Unit       */
      img->unit=CUNIT_HU; break;
    case 18: /* Unknown               */
    default:
      img->unit=CUNIT_UNKNOWN; break;
  }
}

Referenced by ecat63ReadAllToImg(), ecat63ReadPlaneToImg(), imgGetEcat63MHeader(), imgReadEcat63Header(), and imgUnitFromEcat().

imgUnitFromEcat7()

void imgUnitFromEcat7	(	IMG *	img,
		ECAT7_mainheader *	h )

Set IMG calibration unit based on ECAT7 main header.

See also

imgUnitFromEcat, imgUnitId, imgUnitToEcat7, imgSetUnit

Parameters

img	Pointer to target IMG structure.
h	Pointer to source ECAT 7 main header structure.

Definition at line 149 of file imgunits.c.

  {
  if(h->calibration_units==0) { /* Not calibrated */
    img->unit=1;
  } else if(h->calibration_units==1) {
    if(h->calibration_units_label==0)
      img->unit=CUNIT_BQ_PER_ML;
    else
      img->unit=imgUnitId(h->data_units);
  } else if(h->calibration_units==2) {
    img->unit=imgUnitId(h->data_units);
  } else {
    img->unit=CUNIT_UNKNOWN;
  }
}

Referenced by imgGetEcat7MHeader(), and imgUnitFromEcat7().

imgUnitId()

int imgUnitId

(

char *

unit

)

Identify the specified unit string as IMG unit.

See also

imgSetUnit, imgUnitId, petCunitId

Returns

unit id number, -1 invalid input, -3 no match found

Parameters

unit	Pointer to string containing the unit to be identified.

Definition at line 14 of file imgunits.c.

  {
  if(unit==NULL) return(-1);
  if(strlen(unit)==0)                             return CUNIT_UNKNOWN;
  else if(strcasecmp(unit, "unknown")==0)         return CUNIT_UNKNOWN;
  else if(strcasecmp(unit, "cnts/sec")==0)        return CUNIT_CPS;
  else if(strcasecmp(unit, "counts/sec")==0)      return CUNIT_CPS;
  else if(strcasecmp(unit, "ECAT counts/sec")==0) return CUNIT_CPS;
  else if(strcasecmp(unit, "cps")==0)             return CUNIT_CPS;
  else if(strcasecmp(unit, "counts")==0)          return CUNIT_COUNTS;
  else if(strcasecmp(unit, "cnts")==0)            return CUNIT_COUNTS;
  else if(strcasecmp(unit, "kBq/cc")==0)          return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBqcc")==0)           return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBq/mL")==0)          return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBqmL")==0)           return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "sec*kBq/cc")==0)      return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "sec*kBq/mL")==0)      return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "integral")==0)        return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "1/sec")==0)           return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "1/s")==0)             return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "s-1")==0)             return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "1/min")==0)           return CUNIT_PER_MIN;
  else if(strcasecmp(unit, "min-1")==0)           return CUNIT_PER_MIN;
  else if(strcasecmp(unit, "mL/mL")==0)           return CUNIT_ML_PER_ML;
  else if(strcasecmp(unit, "mL/cc")==0)           return CUNIT_ML_PER_ML;
  else if(strcasecmp(unit, "mL/dL")==0)           return CUNIT_ML_PER_DL;
  else if(strcasecmp(unit, "mL/100mL")==0)        return CUNIT_ML_PER_DL;
  else if(strcasecmp(unit, "mL/(mL*min)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*mL)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(cc*min)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*cc)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/mL/min")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/min/mL")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/cc/min")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/min/cc")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(dL*min)")==0)     return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*dL)")==0)     return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(100mL*min)")==0)  return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*100mL)")==0)  return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/dL/min")==0)       return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/min/dL")==0)       return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/100mL/min")==0)    return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/min/100mL")==0)    return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "unitless")==0)        return CUNIT_UNITLESS;
  else if(strcasecmp(unit, "Hounsfield Unit")==0) return CUNIT_HU;
  else if(strcasecmp(unit, "HU")==0)              return CUNIT_HU;
  else if(strcasecmp(unit, "nCi/cc")==0)          return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCicc")==0)           return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCi/mL")==0)          return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCimL")==0)           return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "MBq/cc")==0)          return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBqcc")==0)           return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBq/mL")==0)          return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBqmL")==0)           return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "Bq/cc")==0)           return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "Bqcc")==0)            return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "Bq/mL")==0)           return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "BqmL")==0)            return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "uCi/cc")==0)          return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCicc")==0)           return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCi/mL")==0)          return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCimL")==0)           return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "umol/(100g*min)")==0) return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/(min*100g)")==0) return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/100g/min")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/min/100g")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(100g*min)")==0)   return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(min*100g)")==0)   return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/100g/min")==0)     return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/min/100g")==0)     return CUNIT_MG_PER_MIN_PER_100G;
 
  return(-3);
}

Referenced by imgConvertUnit(), imgSetUnit(), imgUnitFromEcat7(), and imgUnitId().

imgUnitToEcat6()

int imgUnitToEcat6

(

IMG *

img

)

Return ECAT 6.3 calibration unit based on IMG unit

See also

imgUnitFromEcat, imgUnitId, imgUnitToEcat7, imgSetUnit

Returns

calibration unit code number.

Parameters

img	Pointer to IMG structure.

Definition at line 175 of file imgunits.c.

  {
  int ecat_unit;
  switch(img->unit) {
    case CUNIT_CPS: ecat_unit=9; break;
    case CUNIT_COUNTS: ecat_unit=2; break;
    case CUNIT_KBQ_PER_ML: ecat_unit=10; break;
    case CUNIT_SEC_KBQ_PER_ML: ecat_unit=13; break;
    case CUNIT_PER_SEC: ecat_unit=15; break;
    case CUNIT_PER_MIN: ecat_unit=11; break;
    case CUNIT_ML_PER_ML: ecat_unit=16; break;
    case CUNIT_ML_PER_DL: ecat_unit=16; break;
    case CUNIT_ML_PER_ML_PER_MIN: ecat_unit=11; break;
    case CUNIT_ML_PER_DL_PER_MIN: ecat_unit=12; break;
    case CUNIT_UNITLESS: ecat_unit=16; break;
    case CUNIT_HU: ecat_unit=17; break;
    case CUNIT_NCI_PER_ML: ecat_unit=7; break;
    case CUNIT_MBQ_PER_ML: ecat_unit=1; break;
    case CUNIT_BQ_PER_ML: ecat_unit=9; break;
    case CUNIT_UCI_PER_ML: ecat_unit=3; break;
    case CUNIT_UMOL_PER_MIN_PER_100G: ecat_unit=5; break;
    case CUNIT_MG_PER_MIN_PER_100G: ecat_unit=6; break;
    default: ecat_unit=0; break;
  }
  return(ecat_unit);
}

Referenced by ecat63AddImg(), ecat63WriteAllImg(), imgSetEcat63MHeader(), imgSetEcat63SHeader(), and imgUnitToEcat6().

imgUnitToEcat7()

void imgUnitToEcat7	(	IMG *	img,
		ECAT7_mainheader *	h )

Set ECAT 7 main header calibration units based on IMG unit.

See also

imgUnitFromEcat7, imgUnitId, imgUnitToEcat6, imgSetUnit

Parameters

img	Pointer to source IMG structure.
h	Pointer to target ECAT 7 main header structure.

Definition at line 209 of file imgunits.c.

  {
  switch(img->unit) {
    case CUNIT_CPS:
    case CUNIT_COUNTS:
      h->calibration_units=0;
      h->calibration_units_label=1;
      strcpy(h->data_units, "ECAT counts/sec");
      break;
    case CUNIT_KBQ_PER_ML:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "kBq/cc");
      break;
    case CUNIT_SEC_KBQ_PER_ML:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "sec*kBq/mL");
      break;
    case CUNIT_PER_SEC:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "1/sec");
      break;
    case CUNIT_PER_MIN:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "1/min");
      break;
    case CUNIT_ML_PER_ML:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "mL/mL");
      break;
    case CUNIT_ML_PER_DL:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "mL/dL");
      break;
    case CUNIT_ML_PER_ML_PER_MIN:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "mL/(mL*min)");
      break;
    case CUNIT_ML_PER_DL_PER_MIN:
      h->calibration_units=2;
      h->calibration_units_label=1;
      strcpy(h->data_units, "mL/(dL*min)");
      break;
    case CUNIT_UNITLESS:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "unitless");
      break;
    case CUNIT_NCI_PER_ML:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "nCi/cc");
      break;
    case CUNIT_MBQ_PER_ML:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "MBq/cc");
      break;
    case CUNIT_BQ_PER_ML:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "Bq/cc");
      break;
    case CUNIT_UCI_PER_ML:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "uCi/cc");
      break;
    case CUNIT_UMOL_PER_MIN_PER_100G:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "umol/min/100g");
      break;
    case CUNIT_MG_PER_MIN_PER_100G:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "mg/min/100g");
      break;
    case CUNIT_HU:
      h->calibration_units=1;
      h->calibration_units_label=1;
      strcpy(h->data_units, "HU");
      break;
    default:
      h->calibration_units=2;
      h->calibration_units_label=0;
      strcpy(h->data_units, "");
      break;
  }
}

Referenced by imgSetEcat7MHeader(), and imgUnitToEcat7().