tpcnifti.h

Go to the documentation of this file.

#ifndef _TPCNIFTI_H_
#define _TPCNIFTI_H_
/*****************************************************************************/
 
/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
/*****************************************************************************/
#include "tpcextensions.h"
#include "tpcift.h"
#include "tpcfileutil.h"
/*****************************************************************************/
 
/*****************************************************************************/
#define NIFTI1_HEADER_SIZE 348
#define NIFTI1_HEADER_EXTENDER_SIZE 4
#define NIFTI2_HEADER_SIZE 540
#define NIFTI2_HEADER_EXTENDER_SIZE 4
#define ANALYZE_HEADER_SIZE 348
/*****************************************************************************/
#define NIFTI_UNITS_UNKNOWN 0
#define NIFTI_UNITS_METER   1
#define NIFTI_UNITS_MM      2
#define NIFTI_UNITS_MICRON  4
#define NIFTI_UNITS_SEC     8
#define NIFTI_UNITS_MSEC   16
#define NIFTI_UNITS_USEC   24
#define NIFTI_UNITS_HERTZ   32
#define NIFTI_UNITS_PPM     40
#define NIFTI_UNITS_RADS    48
/*****************************************************************************/
#define ANALYZE_DT_NONE 0
#define ANALYZE_DT_UNKNOWN 0
#define ANALYZE_DT_BINARY 1
#define ANALYZE_DT_UNSIGNED_CHAR 2
#define ANALYZE_DT_SIGNED_SHORT 4
#define ANALYZE_DT_SIGNED_INT 8
#define ANALYZE_DT_FLOAT 16
#define ANALYZE_DT_COMPLEX 32
#define ANALYZE_DT_DOUBLE 64
#define ANALYZE_DT_RGB 128
#define ANALYZE_DT_ALL 255

#define NIFTI_DT_NONE ANALYZE_DT_NONE
#define NIFTI_DT_UNKNOWN ANALYZE_DT_UNKNOWN
#define NIFTI_DT_BINARY ANALYZE_DT_BINARY
#define NIFTI_DT_UNSIGNED_CHAR ANALYZE_DT_UNSIGNED_CHAR
#define NIFTI_DT_SIGNED_SHORT ANALYZE_DT_SIGNED_SHORT
#define NIFTI_DT_SIGNED_INT ANALYZE_DT_SIGNED_INT
#define NIFTI_DT_FLOAT ANALYZE_DT_FLOAT
#define NIFTI_DT_COMPLEX ANALYZE_DT_COMPLEX
#define NIFTI_DT_DOUBLE ANALYZE_DT_DOUBLE
#define NIFTI_DT_RGB ANALYZE_DT_RGB
#define NIFTI_DT_ALL ANALYZE_DT_ALL
#define NIFTI_DT_SIGNED_CHAR 256
#define NIFTI_DT_UNSIGNED_SHORT 512
#define NIFTI_DT_UNSIGNED_INT 768
#define NIFTI_DT_LONG_LONG 1024
#define NIFTI_DT_UNSIGNED_LONG_LONG 1280
#define NIFTI_DT_LONG_DOUBLE 1536
#define NIFTI_DT_DOUBLE_PAIR 1792
#define NIFTI_DT_LONG_DOUBLE_PAIR 2048
#define NIFTI_DT_RGBA 2304
/*****************************************************************************/
#define NIFTI_INTENT_NONE 0
#define NIFTI_INTENT_CORREL 2
#define NIFTI_INTENT_TTEST 3
#define NIFTI_INTENT_FTEST 4 
#define NIFTI_INTENT_ZSCORE 5
#define NIFTI_INTENT_CHISQ 6
#define NIFTI_INTENT_BETA 7
#define NIFTI_INTENT_BINOM 8
#define NIFTI_INTENT_GAMMA 9
#define NIFTI_INTENT_POISSON 10
#define NIFTI_INTENT_NORMAL 11
#define NIFTI_INTENT_FTEST_NONC 12
#define NIFTI_INTENT_CHISQ_NONC 13
#define NIFTI_INTENT_LOGISTIC 14
#define NIFTI_INTENT_LAPLACE 15
#define NIFTI_INTENT_UNIFORM 16
#define NIFTI_INTENT_TTEST_NONC 17
#define NIFTI_INTENT_WEIBULL 18
#define NIFTI_INTENT_CHI 19
#define NIFTI_INTENT_INVGAUSS 20
#define NIFTI_INTENT_EXTVAL 21
#define NIFTI_INTENT_PVAL 22
#define NIFTI_INTENT_LOGPVAL 23 
#define NIFTI_INTENT_LOG10PVAL 24
#define NIFTI_INTENT_ESTIMATE 1001
#define NIFTI_INTENT_LABEL 1002
#define NIFTI_INTENT_NEURONAME 1003
#define NIFTI_INTENT_GENMATRIX 1004
#define NIFTI_INTENT_SYMMATRIX 1005
#define NIFTI_INTENT_DISPVECT 1006
#define NIFTI_INTENT_VECTOR 1007
#define NIFTI_INTENT_POINTSET 1008
#define NIFTI_INTENT_TRIANGLE 1009
#define NIFTI_INTENT_QUATERNION 1010
#define NIFTI_INTENT_DIMLESS 1011
/*****************************************************************************/
#define NIFTI_XFORM_UNKNOWN 0
#define NIFTI_XFORM_SCANNER_ANAT 1
#define NIFTI_XFORM_ALIGNED_ANAT 2
#define NIFTI_XFORM_TALAIRACH 3
#define NIFTI_XFORM_MNI_152 4
/*****************************************************************************/
 
/*****************************************************************************/
typedef struct {
  int sizeof_hdr;
  char data_type[10];
  char db_name[18];
  int extents;
  short int session_error;
  char regular;
  char dim_info;

  short int dim[8];
  float intent_p1;
  float intent_p2;
  float intent_p3;
  short int intent_code;
  short int datatype;
  short int bitpix;
  short int slice_start;
  float pixdim[8];
  float vox_offset;
  float scl_slope;
  float scl_inter;
  short int slice_end;
  char slice_code;
  char xyzt_units;
  float cal_max;
  float cal_min;
  float slice_duration;
  float toffset;
  int glmax;
  int glmin;

  char descrip[80];
  char aux_file[24];
  short int qform_code;
  short int sform_code;
  float quatern_b;
  float quatern_c;
  float quatern_d;
  float qoffset_x;
  float qoffset_y;
  float qoffset_z;
  float srow_x[4];
  float srow_y[4];
  float srow_z[4];
  char intent_name[16];
  char magic[4];
} NIFTI_1_HEADER;

typedef struct {
  int sizeof_hdr;
  char magic[8];
  int16_t datatype;
  int16_t bitpix;
  int64_t dim[8];
  double intent_p1;
  double intent_p2;
  double intent_p3;
  double pixdim[8];
  int64_t vox_offset;
  double scl_slope;
  double scl_inter;
  double cal_max;
  double cal_min;
  double slice_duration;
  double toffset;
  int64_t slice_start;
  int64_t slice_end;
  char descrip[80];
  char aux_file[24];
  int qform_code;
  int sform_code;
  double quatern_b;
  double quatern_c;
  double quatern_d;
  double qoffset_x;
  double qoffset_y;
  double qoffset_z;
  double srow_x[4];
  double srow_y[4];
  double srow_z[4];
  int slice_code;
  int xyzt_units;
  int intent_code;
  char intent_name[16];
  char dim_info;
  char unused_str[15];
} NIFTI_2_HEADER;

typedef struct {
  char extension[4];
} NIFTI_EXTENDER;

typedef struct {
  NIFTI_1_HEADER h1;
  NIFTI_2_HEADER h2; 
  NIFTI_EXTENDER e;
  int n;
  int byte_order;
} NIFTI_DSR;
/*****************************************************************************/
 
/*****************************************************************************/
typedef struct {
  int sizeof_hdr;
  char data_type[10];
  char db_name[18];
  int extents;
  short int session_error;
  char regular;
  char hkey_un0;
  short int dim[8];
  short int unused8;
  short int unused9;
  short int unused10;
  short int unused11;
  short int unused12;
  short int unused13;
  short int unused14;
  short int datatype;
  short int bitpix;
  short int dim_un0;
  float pixdim[8];
  float vox_offset;
  float funused1;
  float funused2;
  float funused3;
  float cal_max;
  float cal_min;
  float compressed;
  float verified;
  int glmax;
  int glmin;
  char descrip[80];
  char aux_file[24];
  char orient;
  char originator[10];
  char generated[10];
  char scannum[10];
  char patient_id[10];
  char exp_date[10];
  char exp_time[10];
  char hist_un0[3];
  int views;
  int vols_added;
  int start_field;
  int field_skip;
  int omax;
  int omin;
  int smax;
  int smin;
} ANALYZE_HEADER;

typedef struct {
  ANALYZE_HEADER h;
  int byte_order;
} ANALYZE_DSR;
/*****************************************************************************/
 
/*****************************************************************************/
/* niftiname */
extern void niftiBasename(char *filename);
/*****************************************************************************/
 
/*****************************************************************************/
/* niftiio */
extern int niftiExists(
  const char *filename, char *hdrfile, char *imgfile, char *siffile, NIFTI_DSR *header,
  TPCSTATUS *status
);
extern int niftiReadHeader(const char *filename, NIFTI_DSR *dsr, int verbose);
extern int niftiWriteHeader(const char *filename, NIFTI_DSR *dsr, int verbose);
/*****************************************************************************/
 
/*****************************************************************************/
/* analyzeio */
extern int anaExists(
  const char *filename, char *hdrfile, char *imgfile, char *siffile, ANALYZE_DSR *header,
  TPCSTATUS *status
);
extern int anaReadHeader(const char *filename, ANALYZE_DSR *dsr, int verbose);
/*****************************************************************************/
 
/*****************************************************************************/
#endif /* TPCNIFTI */