tpcimage.h

Go to the documentation of this file.

 
#ifndef _TPCIMAGE_H_
#define _TPCIMAGE_H_
/*****************************************************************************/
 
/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
/*****************************************************************************/
#include "tpcextensions.h"
#include "tpcift.h"
#include "tpcisotope.h"
#include "tpctac.h" // SIF is TAC
/*****************************************************************************/
#include "tpcecat.h"
#include "tpcdcm.h"
#include "tpcnifti.h"
/*****************************************************************************/
 
/*****************************************************************************/
typedef enum {
  IMG_FORMAT_UNKNOWN,     
  IMG_FORMAT_DICOM,       
  IMG_FORMAT_E63,         
  IMG_FORMAT_E7,          
  IMG_FORMAT_E7_2D,       
  IMG_FORMAT_POLARMAP,    
  IMG_FORMAT_ANA,         
  IMG_FORMAT_ANA_L,       
  IMG_FORMAT_INTERFILE,   
  IMG_FORMAT_NIFTI_1D,    
  IMG_FORMAT_NIFTI_1S,    
  IMG_FORMAT_NIFTI_2D,    
  IMG_FORMAT_NIFTI_2S,    
  IMG_FORMAT_MICROPET,    
  IMG_FORMAT_FLAT,        
  IMG_FORMAT_LAST         
} imgformat;
/*****************************************************************************/
 
/*****************************************************************************/
typedef enum {
  IMG_MODALITY_UNKNOWN, 
  IMG_MODALITY_PET,     
  IMG_MODALITY_SPECT,   
  IMG_MODALITY_CT,      
  IMG_MODALITY_MRI,     
  IMG_MODALITY_LAST     
} imgmodality;
/*****************************************************************************/
 
/*****************************************************************************/
typedef enum {
  IMG_CONTENT_UNKNOWN,   
  IMG_CONTENT_IMAGE,     
  IMG_CONTENT_RAW,       
  IMG_CONTENT_ATTN,      
  IMG_CONTENT_POLARMAP,  
  IMG_CONTENT_LAST       
} imgcontent;
/*****************************************************************************/
 
/*****************************************************************************/
typedef struct {
 
  char studyNr[MAX_STUDYNR_LEN+1];
 
  isotope isot;
 
  decaycorrection decayCorrection;
 
  char scanStart[20];
 
  imgcontent content;
 
  imgmodality modality;
 
  imgformat format;
  imgformat oformat;
 
  /* Coordination */
 
  unsigned short int dimt;
  unsigned short int dimx;
  unsigned short int dimy;
  unsigned short int dimz;
 
  float sizex;
  float sizey;
  float sizez;
 
  float gapx;
  float gapy;
  float gapz;
 
  short int xform[2];
  float quatern[6];
  float srow[12];
  float iop[6];
  float ipp[3];
  float mt[12];
 
  /* Pixel data */
 
 
  float *_t;
  float **_x;
  float ***_y;
  float ****_z;
 
  /* Pointers for data to be used */
  float ****m;
  float *p;
 
  /*
   *  Frame times
   */
  float *x1;
  float *x2;
  float *x;
 
  /*
   *  Frame weights
   */
  weights weighting;
  float *weight;
  float *prompts;
  float *randoms;
 
  /*
   *  Data units
   */
  unit cunit;
  unit tunit;
 
  IFT ih;
  IFT oh;
 
} IMG;
/*****************************************************************************/
 
/*****************************************************************************/
/* image */
extern void imgInit(IMG *image);
extern void imgFree(IMG *image);
extern int imgAllocate(
  IMG *img, const unsigned int dimz, const unsigned int dimy,
  const unsigned int dimx, const unsigned int dimt, TPCSTATUS *status
);
extern char *imgContentDescr(imgcontent c);
extern char *imgModalityDescr(imgmodality c);
extern int imgHasData(IMG *img);
extern int imgHasTimes(IMG *img);
extern int imgHasCounts(IMG *img);
extern int imgHasWeights(IMG *img);
extern void imgContents(IMG *img, FILE *fp);
extern unsigned long long imgNaNs(IMG *img, int fix);
extern int imgMinMax(IMG *img, float *minvalue, float *maxvalue);
extern int imgXRange(IMG *img, double *xmin, double *xmax);
/*****************************************************************************/
 
/*****************************************************************************/
/* imageheader */
extern int imgCopyHeader(IMG *img1, IMG *img2);
extern int imgFillOHeader(IMG *img, TPCSTATUS *status);
/*****************************************************************************/
 
/*****************************************************************************/
/* imageio */
extern char *imgFormatDescr(imgformat c);
extern char *imgDefaultExtension(imgformat c);
extern int imgRead(IMG *img, const char *fname, TPCSTATUS *status);
extern int imgWrite(IMG *img, const char *fname, TPCSTATUS *status);
extern imgformat imgFormatFromFName(const char *fname, TPCSTATUS *status);
extern imgformat imgFormatIdentify(const char *s);
/*****************************************************************************/
 
/*****************************************************************************/
/* imagedcm */
extern int imgReadDICOM(IMG *img, const char *fname, TPCSTATUS *status);
extern int imgWriteDICOM(IMG *img, const char *fname, TPCSTATUS *status);
/*****************************************************************************/
 
/*****************************************************************************/
/* imagenii */
extern int niftiCreateFNames(
  const char *filename, char *hdrfile, char *imgfile, char *siffile, int fileformat
);
extern int imgReadNifti(IMG *img, const char *fname, TPCSTATUS *status);
extern int imgWriteNifti(IMG *img, const char *fname, TPCSTATUS *status);
extern int imgSetNiftiHeader(IMG *img, NIFTI_DSR *dsr, int verbose);
extern int imgGetNiftiHeader(IMG *img, NIFTI_DSR *dsr, int verbose);
/*****************************************************************************/
 
/*****************************************************************************/
/* imagesif */
extern int imgToSIF(IMG *img, TAC *sif, int copy_header, int copy_frames, int copy_counts, int verbose);
extern int imgFromSIF(IMG *img, TAC *sif, int copy_header, int copy_frames, int copy_counts, int verbose);
/*****************************************************************************/
 
/*****************************************************************************/
/* imagecomp */
extern int imgCompareMatrixSize(IMG *d1, IMG *d2);
extern int imgCompareUnit(IMG *d1, IMG *d2, TPCSTATUS *status);
extern int imgCompareConc(
  IMG *d1, IMG *d2, const float test_abs, const float test_rel, TPCSTATUS *status);
extern int imgCompareTimes(
  IMG *d1, IMG *d2, const float test_abs, const float test_rel, TPCSTATUS *status);
/*****************************************************************************/
 
/*****************************************************************************/
#endif /* TPCIMAGE */