forked from NIFTI-Imaging/nifti_clib
-
Notifications
You must be signed in to change notification settings - Fork 0
/
nifti2_io.h
743 lines (583 loc) · 34.7 KB
/
nifti2_io.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
/** \file nifti2_io.h
\brief Data structures for using nifti2_io API.
- Written by Bob Cox, SSCC NIMH
- Revisions by Rick Reynolds, SSCC NIMH
*/
#ifndef _NIFTI2_IO_HEADER_
#define _NIFTI2_IO_HEADER_
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <limits.h>
#include <ctype.h>
#include <inttypes.h>
#ifndef DONT_INCLUDE_ANALYZE_STRUCT
#define DONT_INCLUDE_ANALYZE_STRUCT /*** not needed herein ***/
#endif
#include "nifti1.h" /*** NIFTI-1 header specification ***/
#include "nifti2.h" /*** NIFTI-2 header specification ***/
#include <znzlib.h>
/*=================*/
#ifdef __cplusplus
extern "C" {
#endif
/*=================*/
/*****===================================================================*****/
/***** File nifti2_io.h == Declarations for nifti2_io.c *****/
/*****...................................................................*****/
/***** This code is a modification of nifti1_io.h. *****/
/*****...................................................................*****/
/***** This code is released to the public domain. *****/
/*****...................................................................*****/
/***** Author: Robert W Cox, SSCC/DIRP/NIMH/NIH/DHHS/USA/EARTH *****/
/***** Date: August 2003 *****/
/*****...................................................................*****/
/***** Neither the National Institutes of Health (NIH), nor any of its *****/
/***** employees imply any warranty of usefulness of this software for *****/
/***** any purpose, and do not assume any liability for damages, *****/
/***** incidental or otherwise, caused by any use of this document. *****/
/*****===================================================================*****/
/* ......................................................................
Modified by: Mark Jenkinson (FMRIB Centre, University of Oxford, UK)
Date: July/August 2004
Mainly adding low-level IO and changing things to allow gzipped files
to be read and written
Full backwards compatability should have been maintained
......................................................................
Modified by: Rick Reynolds (SSCC/DIRP/NIMH, National Institutes of Health)
Date: December 2004
Modified and added many routines for I/O, particularly involving
extensions and nifti_brick_list.
......................................................................
Modified by: Rick Reynolds (SSCC/DIRP/NIMH, National Institutes of Health)
Date: August 2013
Converted to be based on nifti_2_header.
** NOT BACKWARD COMPATABLE **
These routines will read/write both NIFTI-1 and NIFTI-2 image files,
but modification to the _calling_ routies is necessary, since:
a. the main nifti_image type has changed (to nifti2_image)
b. some image field types have been altered (to have larger size)
c. some routines have been changed to apply to multiple NIFTI types
*/
/********************** Some sample data structures **************************/
typedef struct { /** 4x4 matrix struct **/
float m[4][4] ;
} mat44 ;
typedef struct { /** 3x3 matrix struct **/
float m[3][3] ;
} mat33 ;
typedef struct { /** 4x4 matrix struct (double) **/
double m[4][4] ;
} nifti_dmat44 ;
typedef struct { /** 3x3 matrix struct (double) **/
double m[3][3] ;
} nifti_dmat33 ;
/*...........................................................................*/
/*! \enum analyze_75_orient_code
* \brief Old-style analyze75 orientation
* codes.
*/
typedef enum _analyze75_orient_code {
a75_transverse_unflipped = 0,
a75_coronal_unflipped = 1,
a75_sagittal_unflipped = 2,
a75_transverse_flipped = 3,
a75_coronal_flipped = 4,
a75_sagittal_flipped = 5,
a75_orient_unknown = 6
} analyze_75_orient_code;
/*! \struct nifti_image
\brief High level data structure for open nifti datasets in the
nifti2_io API. Note that this structure is not part of the
nifti2 format definition; it is used to implement one API
for reading/writing datasets in the nifti1 or nifti2 formats.
Field types changed for NIFTI-2 (note: ALL floats to doubles):
nx, ny, ..., nw, dim, nvox,
dx, dy, ..., dw, pixdim,
scl_slope, scl_inter, cal_min, cal_max,
slice_start, slice_end, slice_duration,
quatern_b,c,d, qoffset_x,y,z, qfac,
qto_xyz,ijk, sto_xyz,ijk,
toffset, intent_p1,2,3, iname_offset
*/
typedef struct { /*!< Image storage struct **/
int64_t ndim ; /*!< last dimension greater than 1 (1..7) */
int64_t nx ; /*!< dimensions of grid array */
int64_t ny ; /*!< dimensions of grid array */
int64_t nz ; /*!< dimensions of grid array */
int64_t nt ; /*!< dimensions of grid array */
int64_t nu ; /*!< dimensions of grid array */
int64_t nv ; /*!< dimensions of grid array */
int64_t nw ; /*!< dimensions of grid array */
int64_t dim[8] ; /*!< dim[0]=ndim, dim[1]=nx, etc. */
int64_t nvox ; /*!< number of voxels = nx*ny*nz*...*nw */
int nbyper ; /*!< bytes per voxel, matches datatype */
int datatype ; /*!< type of data in voxels: DT_* code */
double dx ; /*!< grid spacings */
double dy ; /*!< grid spacings */
double dz ; /*!< grid spacings */
double dt ; /*!< grid spacings */
double du ; /*!< grid spacings */
double dv ; /*!< grid spacings */
double dw ; /*!< grid spacings */
double pixdim[8] ; /*!< pixdim[1]=dx, etc. */
double scl_slope ; /*!< scaling parameter - slope */
double scl_inter ; /*!< scaling parameter - intercept */
double cal_min ; /*!< calibration parameter, minimum */
double cal_max ; /*!< calibration parameter, maximum */
int qform_code ; /*!< codes for (x,y,z) space meaning */
int sform_code ; /*!< codes for (x,y,z) space meaning */
int freq_dim ; /*!< indexes (1,2,3, or 0) for MRI */
int phase_dim ; /*!< directions in dim[]/pixdim[] */
int slice_dim ; /*!< directions in dim[]/pixdim[] */
int slice_code ; /*!< code for slice timing pattern */
int64_t slice_start ; /*!< index for start of slices */
int64_t slice_end ; /*!< index for end of slices */
double slice_duration ; /*!< time between individual slices */
/*! quaternion transform parameters
[when writing a dataset, these are used for qform, NOT qto_xyz] */
double quatern_b , quatern_c , quatern_d ,
qoffset_x , qoffset_y , qoffset_z ,
qfac ;
nifti_dmat44 qto_xyz ; /*!< qform: transform (i,j,k) to (x,y,z) */
nifti_dmat44 qto_ijk ; /*!< qform: transform (x,y,z) to (i,j,k) */
nifti_dmat44 sto_xyz ; /*!< sform: transform (i,j,k) to (x,y,z) */
nifti_dmat44 sto_ijk ; /*!< sform: transform (x,y,z) to (i,j,k) */
double toffset ; /*!< time coordinate offset */
int xyz_units ; /*!< dx,dy,dz units: NIFTI_UNITS_* code */
int time_units ; /*!< dt units: NIFTI_UNITS_* code */
int nifti_type ; /*!< see NIFTI_FTYPE_* codes, below:
0==ANALYZE,
1==NIFTI-1 (1 file),
2==NIFTI-1 (2 files),
3==NIFTI-ASCII (1 file)
4==NIFTI-2 (1 file),
5==NIFTI-2 (2 files) */
int intent_code ; /*!< statistic type (or something) */
double intent_p1 ; /*!< intent parameters */
double intent_p2 ; /*!< intent parameters */
double intent_p3 ; /*!< intent parameters */
char intent_name[16] ; /*!< optional description of intent data */
char descrip[80] ; /*!< optional text to describe dataset */
char aux_file[24] ; /*!< auxiliary filename */
char *fname ; /*!< header filename (.hdr or .nii) */
char *iname ; /*!< image filename (.img or .nii) */
int64_t iname_offset ; /*!< offset into iname where data starts */
int swapsize ; /*!< swap unit in image data (might be 0) */
int byteorder ; /*!< byte order on disk (MSB_ or LSB_FIRST) */
void *data ; /*!< pointer to data: nbyper*nvox bytes */
int num_ext ; /*!< number of extensions in ext_list */
nifti1_extension * ext_list ; /*!< array of extension structs (with data) */
analyze_75_orient_code analyze75_orient; /*!< for old analyze files, orient */
} nifti_image ;
/* allow clarity */
typedef nifti_image nifti2_image;
typedef struct {
int ndim ; /*!< last dimension greater than 1 (1..7) */
int nx ; /*!< dimensions of grid array */
int ny ; /*!< dimensions of grid array */
int nz ; /*!< dimensions of grid array */
int nt ; /*!< dimensions of grid array */
int nu ; /*!< dimensions of grid array */
int nv ; /*!< dimensions of grid array */
int nw ; /*!< dimensions of grid array */
int dim[8] ; /*!< dim[0]=ndim, dim[1]=nx, etc. */
int64_t nvox ; /*!< number of voxels = nx*ny*nz*...*nw */
int nbyper ; /*!< bytes per voxel, matches datatype */
int datatype ; /*!< type of data in voxels: DT_* code */
float dx ; /*!< grid spacings */
float dy ; /*!< grid spacings */
float dz ; /*!< grid spacings */
float dt ; /*!< grid spacings */
float du ; /*!< grid spacings */
float dv ; /*!< grid spacings */
float dw ; /*!< grid spacings */
float pixdim[8] ; /*!< pixdim[1]=dx, etc. */
float scl_slope ; /*!< scaling parameter - slope */
float scl_inter ; /*!< scaling parameter - intercept */
float cal_min ; /*!< calibration parameter, minimum */
float cal_max ; /*!< calibration parameter, maximum */
int qform_code ; /*!< codes for (x,y,z) space meaning */
int sform_code ; /*!< codes for (x,y,z) space meaning */
int freq_dim ; /*!< indexes (1,2,3, or 0) for MRI */
int phase_dim ; /*!< directions in dim[]/pixdim[] */
int slice_dim ; /*!< directions in dim[]/pixdim[] */
int slice_code ; /*!< code for slice timing pattern */
int slice_start ; /*!< index for start of slices */
int slice_end ; /*!< index for end of slices */
float slice_duration ; /*!< time between individual slices */
/*! quaternion transform parameters
[when writing a dataset, these are used for qform, NOT qto_xyz] */
float quatern_b , quatern_c , quatern_d ,
qoffset_x , qoffset_y , qoffset_z ,
qfac ;
mat44 qto_xyz ; /*!< qform: transform (i,j,k) to (x,y,z) */
mat44 qto_ijk ; /*!< qform: transform (x,y,z) to (i,j,k) */
mat44 sto_xyz ; /*!< sform: transform (i,j,k) to (x,y,z) */
mat44 sto_ijk ; /*!< sform: transform (x,y,z) to (i,j,k) */
float toffset ; /*!< time coordinate offset */
int xyz_units ; /*!< dx,dy,dz units: NIFTI_UNITS_* code */
int time_units ; /*!< dt units: NIFTI_UNITS_* code */
int nifti_type ; /*!< 0==ANALYZE, 1==NIFTI-1 (1 file),
2==NIFTI-1 (2 files),
3==NIFTI-ASCII (1 file) */
int intent_code ; /*!< statistic type (or something) */
float intent_p1 ; /*!< intent parameters */
float intent_p2 ; /*!< intent parameters */
float intent_p3 ; /*!< intent parameters */
char intent_name[16] ; /*!< optional description of intent data */
char descrip[80] ; /*!< optional text to describe dataset */
char aux_file[24] ; /*!< auxiliary filename */
char *fname ; /*!< header filename (.hdr or .nii) */
char *iname ; /*!< image filename (.img or .nii) */
int iname_offset ; /*!< offset into iname where data starts */
int swapsize ; /*!< swap unit in image data (might be 0) */
int byteorder ; /*!< byte order on disk (MSB_ or LSB_FIRST) */
void *data ; /*!< pointer to data: nbyper*nvox bytes */
int num_ext ; /*!< number of extensions in ext_list */
nifti1_extension * ext_list ; /*!< array of extension structs (with data) */
analyze_75_orient_code analyze75_orient; /*!< for old analyze files, orient */
} nifti1_image ;
/* struct for return from nifti_image_read_bricks() */
typedef struct {
int64_t nbricks; /* the number of allocated pointers in 'bricks' */
int64_t bsize; /* the length of each data block, in bytes */
void ** bricks; /* array of pointers to data blocks */
} nifti_brick_list;
/*****************************************************************************/
/*------------------ NIfTI version of ANALYZE 7.5 structure -----------------*/
/* (based on fsliolib/dbh.h, but updated for version 7.5) */
typedef struct {
/* header info fields - describes the header overlap with NIfTI */
/* ------------------ */
int sizeof_hdr; /* 0 + 4 same */
char data_type[10]; /* 4 + 10 same */
char db_name[18]; /* 14 + 18 same */
int extents; /* 32 + 4 same */
short int session_error; /* 36 + 2 same */
char regular; /* 38 + 1 same */
char hkey_un0; /* 39 + 1 40 bytes */
/* image dimension fields - describes image sizes */
short int dim[8]; /* 0 + 16 same */
short int unused8; /* 16 + 2 intent_p1... */
short int unused9; /* 18 + 2 ... */
short int unused10; /* 20 + 2 intent_p2... */
short int unused11; /* 22 + 2 ... */
short int unused12; /* 24 + 2 intent_p3... */
short int unused13; /* 26 + 2 ... */
short int unused14; /* 28 + 2 intent_code */
short int datatype; /* 30 + 2 same */
short int bitpix; /* 32 + 2 same */
short int dim_un0; /* 34 + 2 slice_start */
float pixdim[8]; /* 36 + 32 same */
float vox_offset; /* 68 + 4 same */
float funused1; /* 72 + 4 scl_slope */
float funused2; /* 76 + 4 scl_inter */
float funused3; /* 80 + 4 slice_end, */
/* slice_code, */
/* xyzt_units */
float cal_max; /* 84 + 4 same */
float cal_min; /* 88 + 4 same */
float compressed; /* 92 + 4 slice_duration */
float verified; /* 96 + 4 toffset */
int glmax,glmin; /* 100 + 8 108 bytes */
/* data history fields - optional */
char descrip[80]; /* 0 + 80 same */
char aux_file[24]; /* 80 + 24 same */
char orient; /* 104 + 1 NO GOOD OVERLAP */
char originator[10]; /* 105 + 10 FROM HERE DOWN... */
char generated[10]; /* 115 + 10 */
char scannum[10]; /* 125 + 10 */
char patient_id[10]; /* 135 + 10 */
char exp_date[10]; /* 145 + 10 */
char exp_time[10]; /* 155 + 10 */
char hist_un0[3]; /* 165 + 3 */
int views; /* 168 + 4 */
int vols_added; /* 172 + 4 */
int start_field; /* 176 + 4 */
int field_skip; /* 180 + 4 */
int omax, omin; /* 184 + 8 */
int smax, smin; /* 192 + 8 200 bytes */
} nifti_analyze75; /* total: 348 bytes */
/*****************************************************************************/
/*--------------- Prototypes of functions defined in this file --------------*/
char const * nifti_datatype_string ( int dt ) ;
char const *nifti_units_string ( int uu ) ;
char const *nifti_intent_string ( int ii ) ;
char const *nifti_xform_string ( int xx ) ;
char const *nifti_slice_string ( int ss ) ;
char const *nifti_orientation_string( int ii ) ;
int nifti_is_inttype( int dt ) ;
mat44 nifti_mat44_inverse ( mat44 R ) ;
mat44 nifti_mat44_mul ( mat44 A , mat44 B );
nifti_dmat44 nifti_dmat44_inverse( nifti_dmat44 R ) ;
int nifti_mat44_to_dmat44(mat44 * fm, nifti_dmat44 * dm);
int nifti_dmat44_to_mat44(nifti_dmat44 * dm, mat44 * fm);
nifti_dmat44 nifti_dmat44_mul ( nifti_dmat44 A , nifti_dmat44 B );
nifti_dmat33 nifti_dmat33_inverse( nifti_dmat33 R ) ;
nifti_dmat33 nifti_dmat33_polar ( nifti_dmat33 A ) ;
double nifti_dmat33_rownorm( nifti_dmat33 A ) ;
double nifti_dmat33_colnorm( nifti_dmat33 A ) ;
double nifti_dmat33_determ ( nifti_dmat33 R ) ;
nifti_dmat33 nifti_dmat33_mul ( nifti_dmat33 A , nifti_dmat33 B ) ;
mat33 nifti_mat33_inverse( mat33 R ) ;
mat33 nifti_mat33_polar ( mat33 A ) ;
float nifti_mat33_rownorm( mat33 A ) ;
float nifti_mat33_colnorm( mat33 A ) ;
float nifti_mat33_determ ( mat33 R ) ;
mat33 nifti_mat33_mul ( mat33 A , mat33 B ) ;
void nifti_swap_2bytes ( int64_t n , void *ar ) ;
void nifti_swap_4bytes ( int64_t n , void *ar ) ;
void nifti_swap_8bytes ( int64_t n , void *ar ) ;
void nifti_swap_16bytes( int64_t n , void *ar ) ;
void nifti_swap_Nbytes ( int64_t n , int siz , void *ar ) ;
int nifti_datatype_is_valid (int dtype, int for_nifti);
int nifti_datatype_from_string (const char * name);
const char * nifti_datatype_to_string(int dtype);
int nifti_header_version (const char * buf, size_t nbytes);
int64_t nifti_get_filesize( const char *pathname ) ;
void swap_nifti_header ( void * hdr , int ni_ver ) ;
void old_swap_nifti_header( struct nifti_1_header *h , int is_nifti );
void nifti_swap_as_analyze( nifti_analyze75 *h );
void nifti_swap_as_nifti1( nifti_1_header *h );
void nifti_swap_as_nifti2( nifti_2_header *h );
/* main read/write routines */
nifti_image *nifti_image_read_bricks(const char *hname , int64_t nbricks,
const int64_t *blist, nifti_brick_list * NBL);
int nifti_image_load_bricks(nifti_image *nim , int64_t nbricks,
const int64_t *blist, nifti_brick_list * NBL);
void nifti_free_NBL( nifti_brick_list * NBL );
nifti_image *nifti_image_read ( const char *hname , int read_data);
int nifti_image_load ( nifti_image *nim);
void nifti_image_unload ( nifti_image *nim);
void nifti_image_free ( nifti_image *nim);
int64_t nifti_read_collapsed_image( nifti_image * nim,
const int64_t dims[8], void ** data);
int64_t nifti_read_subregion_image(nifti_image *nim, const int64_t *start_index,
const int64_t *region_size, void ** data);
void nifti_image_write ( nifti_image * nim ) ;
void nifti_image_write_bricks(nifti_image * nim,
const nifti_brick_list * NBL);
void nifti_image_infodump( const nifti_image * nim ) ;
void nifti_disp_lib_hist( int ver ) ; /* to display library history */
void nifti_disp_lib_version( void ) ; /* to display library version */
int nifti_disp_matrix_orient( const char * mesg, nifti_dmat44 mat );
int nifti_disp_type_list( int which );
char * nifti_image_to_ascii ( const nifti_image * nim ) ;
nifti_image *nifti_image_from_ascii( const char * str, int * bytes_read ) ;
int64_t nifti_get_volsize(const nifti_image *nim) ;
/* basic file operations */
int nifti_set_filenames(nifti_image * nim, const char * prefix, int check,
int set_byte_order);
char * nifti_makehdrname (const char * prefix, int nifti_type, int check,
int comp);
char * nifti_makeimgname (const char * prefix, int nifti_type, int check,
int comp);
int is_nifti_file (const char *hname);
char * nifti_find_file_extension(const char * name);
int nifti_is_complete_filename(const char* fname);
int nifti_validfilename(const char* fname);
int disp_nifti_1_header(const char * info, const nifti_1_header * hp ) ;
int disp_nifti_2_header( const char * info, const nifti_2_header * hp ) ;
void nifti_set_debug_level( int level ) ;
void nifti_set_skip_blank_ext( int skip ) ;
void nifti_set_allow_upper_fext( int allow ) ;
int nifti_get_alter_cifti( void );
void nifti_set_alter_cifti( int alter_cifti );
int nifti_alter_cifti_dims(nifti_image * nim);
int valid_nifti_brick_list(nifti_image * nim , int64_t nbricks,
const int64_t * blist, int disp_error);
/* znzFile operations */
znzFile nifti_image_open(const char * hname, char * opts, nifti_image ** nim);
znzFile nifti_image_write_hdr_img(nifti_image *nim, int write_data,
const char* opts);
znzFile nifti_image_write_hdr_img2( nifti_image *nim , int write_opts ,
const char* opts, znzFile imgfile, const nifti_brick_list * NBL);
int64_t nifti_read_buffer(znzFile fp, void* dataptr, int64_t ntot,
nifti_image *nim);
int nifti_write_all_data(znzFile fp, nifti_image * nim,
const nifti_brick_list * NBL);
int64_t nifti_write_buffer(znzFile fp, const void * buffer, int64_t numbytes);
nifti_image *nifti_read_ascii_image(znzFile fp, const char *fname, int flen,
int read_data);
znzFile nifti_write_ascii_image(nifti_image *nim, const nifti_brick_list * NBL,
const char * opts, int write_data, int leave_open);
void nifti_datatype_sizes( int datatype , int *nbyper, int *swapsize ) ;
void nifti_dmat44_to_quatern(nifti_dmat44 R ,
double *qb, double *qc, double *qd,
double *qx, double *qy, double *qz,
double *dx, double *dy, double *dz, double *qfac);
nifti_dmat44 nifti_quatern_to_dmat44( double qb, double qc, double qd,
double qx, double qy, double qz,
double dx, double dy, double dz, double qfac );
nifti_dmat44 nifti_make_orthog_dmat44( double r11, double r12, double r13 ,
double r21, double r22, double r23 ,
double r31, double r32, double r33 ) ;
void nifti_mat44_to_quatern( mat44 R ,
float *qb, float *qc, float *qd,
float *qx, float *qy, float *qz,
float *dx, float *dy, float *dz, float *qfac ) ;
mat44 nifti_quatern_to_mat44( float qb, float qc, float qd,
float qx, float qy, float qz,
float dx, float dy, float dz, float qfac );
mat44 nifti_make_orthog_mat44( float r11, float r12, float r13 ,
float r21, float r22, float r23 ,
float r31, float r32, float r33 ) ;
int nifti_short_order(void) ; /* CPU byte order */
/* Orientation codes that might be returned from nifti_mat44_to_orientation().*/
#define NIFTI_L2R 1 /* Left to Right */
#define NIFTI_R2L 2 /* Right to Left */
#define NIFTI_P2A 3 /* Posterior to Anterior */
#define NIFTI_A2P 4 /* Anterior to Posterior */
#define NIFTI_I2S 5 /* Inferior to Superior */
#define NIFTI_S2I 6 /* Superior to Inferior */
void nifti_mat44_to_orientation( mat44 R , int *icod, int *jcod, int *kcod ) ;
void nifti_dmat44_to_orientation( nifti_dmat44 R,
int *icod, int *jcod, int *kcod ) ;
/*--------------------- Low level IO routines ------------------------------*/
char * nifti_findhdrname (const char* fname);
char * nifti_findimgname (const char* fname , int nifti_type);
int nifti_is_gzfile (const char* fname);
char * nifti_makebasename(const char* fname);
/* other routines */
int nifti_convert_nim2n1hdr(const nifti_image* nim, nifti_1_header * hdr);
int nifti_convert_nim2n2hdr(const nifti_image* nim, nifti_2_header * hdr);
nifti_1_header * nifti_make_new_n1_header(const int64_t arg_dims[], int arg_dtype);
nifti_2_header * nifti_make_new_n2_header(const int64_t arg_dims[], int arg_dtype);
void * nifti_read_header(const char *hname, int *nver, int check);
nifti_1_header * nifti_read_n1_hdr(const char *hname, int *swapped, int check);
nifti_2_header * nifti_read_n2_hdr(const char *hname, int *swapped, int check);
nifti_image * nifti_copy_nim_info(const nifti_image * src);
nifti_image * nifti_make_new_nim(const int64_t dims[], int datatype,
int data_fill);
nifti_image * nifti_simple_init_nim(void);
nifti_image * nifti_convert_n1hdr2nim(nifti_1_header nhdr,const char *fname);
nifti_image * nifti_convert_n2hdr2nim(nifti_2_header nhdr,const char *fname);
int nifti_looks_like_cifti(nifti_image * nim);
int nifti_hdr1_looks_good (const nifti_1_header * hdr);
int nifti_hdr2_looks_good (const nifti_2_header * hdr);
int nifti_is_valid_datatype (int dtype);
int nifti_is_valid_ecode (int ecode);
int nifti_nim_is_valid (nifti_image * nim, int complain);
int nifti_nim_has_valid_dims (nifti_image * nim, int complain);
int is_valid_nifti_type (int nifti_type);
int nifti_test_datatype_sizes (int verb);
int nifti_type_and_names_match (nifti_image * nim, int show_warn);
int nifti_update_dims_from_array(nifti_image * nim);
void nifti_set_iname_offset (nifti_image *nim, int nifti_ver);
int nifti_set_type_from_names (nifti_image * nim);
int nifti_add_extension(nifti_image * nim, const char * data, int len,
int ecode );
int nifti_compiled_with_zlib (void);
int nifti_copy_extensions (nifti_image *nim_dest,const nifti_image *nim_src);
int nifti_free_extensions (nifti_image *nim);
int64_t * nifti_get_int64list(int64_t nvals , const char *str);
int * nifti_get_intlist (int nvals , const char *str);
char * nifti_strdup (const char *str);
int valid_nifti_extensions(const nifti_image *nim);
int nifti_valid_header_size(int ni_ver, int whine);
/*-------------------- Some C convenience macros ----------------------------*/
/* NIfTI-1.1 extension codes:
see http://nifti.nimh.nih.gov/nifti-1/documentation/faq#Q21 */
#define NIFTI_ECODE_IGNORE 0 /* changed from UNKNOWN, 29 June 2005 */
#define NIFTI_ECODE_DICOM 2 /* intended for raw DICOM attributes */
#define NIFTI_ECODE_AFNI 4 /* Robert W Cox: rwcox@nih.gov
https://afni.nimh.nih.gov/afni */
#define NIFTI_ECODE_COMMENT 6 /* plain ASCII text only */
#define NIFTI_ECODE_XCEDE 8 /* David B Keator: dbkeator@uci.edu
http://www.nbirn.net/Resources
/Users/Applications/
/xcede/index.htm */
#define NIFTI_ECODE_JIMDIMINFO 10 /* Mark A Horsfield:
mah5@leicester.ac.uk
http://someplace/something */
#define NIFTI_ECODE_WORKFLOW_FWDS 12 /* Kate Fissell: fissell@pitt.edu
http://kraepelin.wpic.pitt.edu
/~fissell/NIFTI_ECODE_WORKFLOW_FWDS
/NIFTI_ECODE_WORKFLOW_FWDS.html */
#define NIFTI_ECODE_FREESURFER 14 /* http://surfer.nmr.mgh.harvard.edu */
#define NIFTI_ECODE_PYPICKLE 16 /* embedded Python objects
http://niftilib.sourceforge.net
/pynifti */
/* LONI MiND codes: http://www.loni.ucla.edu/twiki/bin/view/Main/MiND */
#define NIFTI_ECODE_MIND_IDENT 18 /* Vishal Patel: vishal.patel@ucla.edu*/
#define NIFTI_ECODE_B_VALUE 20
#define NIFTI_ECODE_SPHERICAL_DIRECTION 22
#define NIFTI_ECODE_DT_COMPONENT 24
#define NIFTI_ECODE_SHC_DEGREEORDER 26 /* end LONI MiND codes */
#define NIFTI_ECODE_VOXBO 28 /* Dan Kimberg: www.voxbo.org */
#define NIFTI_ECODE_CARET 30 /* John Harwell: john@brainvis.wustl.edu
http://brainvis.wustl.edu/wiki
/index.php/Caret:Documentation
:CaretNiftiExtension */
#define NIFTI_ECODE_CIFTI 32 /* CIFTI-2_Main_FINAL_1March2014.pdf */
#define NIFTI_ECODE_VARIABLE_FRAME_TIMING 34
/* 36 is currently unassigned, waiting on NIFTI_ECODE_AGILENT_PROCPAR */
#define NIFTI_ECODE_EVAL 38 /* Munster University Hospital */
/* http://www.mathworks.com/matlabcentral/fileexchange/42997-dicom-to-nifti-converter */
#define NIFTI_ECODE_MATLAB 40 /* MATLAB extension */
/* Quantiphyse extension
https://quantiphyse.readthedocs.io/en/latest/advanced/nifti_extension.html*/
#define NIFTI_ECODE_QUANTIPHYSE 42 /* Quantiphyse extension */
#define NIFTI_MAX_ECODE 42 /******* maximum extension code *******/
/* nifti_type file codes */
#define NIFTI_FTYPE_ANALYZE 0 /* old ANALYZE */
#define NIFTI_FTYPE_NIFTI1_1 1 /* NIFTI-1 */
#define NIFTI_FTYPE_NIFTI1_2 2
#define NIFTI_FTYPE_ASCII 3
#define NIFTI_FTYPE_NIFTI2_1 4 /* NIFTI-2 */
#define NIFTI_FTYPE_NIFTI2_2 5
#define NIFTI_MAX_FTYPE 5 /* this should match the maximum code */
/*------------------------------------------------------------------------*/
/*-- the rest of these apply only to nifti2_io.c, check for _NIFTI2_IO_C_ */
#ifdef _NIFTI2_IO_C_
typedef struct {
int debug; /*!< debug level for status reports */
int skip_blank_ext; /*!< skip extender if no extensions */
int allow_upper_fext; /*!< allow uppercase file extensions */
int alter_cifti; /*!< convert CIFTI dimensions */
} nifti_global_options;
typedef struct {
int type; /* should match the NIFTI_TYPE_ #define */
int nbyper; /* bytes per value, matches nifti_image */
int swapsize; /* bytes per swap piece, matches nifti_image */
char const * const name; /* text string to match #define */
} nifti_type_ele;
#undef LNI_FERR /* local nifti file error, to be compact and repetative */
#define LNI_FERR(func,msg,file) \
fprintf(stderr,"** ERROR (%s): %s '%s'\n",func,msg,file)
#undef swap_2
#undef swap_4
#define swap_2(s) nifti_swap_2bytes(1,&(s)) /* s: 2-byte short; swap in place */
#define swap_4(v) nifti_swap_4bytes(1,&(v)) /* v: 4-byte value; swap in place */
/***** isfinite() is a C99 macro, which is
present in many C implementations already *****/
#undef IS_GOOD_FLOAT
#undef FIXED_FLOAT
#ifdef isfinite /* use isfinite() to check floats/doubles for goodness */
# define IS_GOOD_FLOAT(x) isfinite(x) /* check if x is a "good" float */
# define FIXED_FLOAT(x) (isfinite(x) ? (x) : 0) /* fixed if bad */
#else
# define IS_GOOD_FLOAT(x) 1 /* don't check it */
# define FIXED_FLOAT(x) (x) /* don't fix it */
#endif
#undef ASSIF /* assign v to *p, if possible */
#define ASSIF(p,v) if( (p)!=NULL ) *(p) = (v)
#undef MSB_FIRST
#undef LSB_FIRST
#undef REVERSE_ORDER
#define LSB_FIRST 1
#define MSB_FIRST 2
#define REVERSE_ORDER(x) (3-(x)) /* convert MSB_FIRST <--> LSB_FIRST */
#define LNI_MAX_NIA_EXT_LEN 100000 /* consider a longer extension invalid */
#undef NIFTI_IS_16_BIT_INT
#define NIFTI_IS_16_BIT_INT(x) ((x) <= 32767 && (x) >= -32768)
#endif /* _NIFTI2_IO_C_ section */
/*------------------------------------------------------------------------*/
/*=================*/
#ifdef __cplusplus
}
#endif
/*=================*/
#endif /* _NIFTI2_IO_HEADER_ */