Skip to content.

AFNI/NIfTI Server

Sections

• Home
• News
• board
• afni
• suma
• nifti
• sscc
• old

Personal tools

• SSCC Staff
• Log in

You are here: Home » AFNI » Documentation » source

            Search  

---

nifti1_io.h File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include "nifti1.h"
#include "znzlib/znzlib.h"

Go to the source code of this file.

Data Structures
struct	mat33
struct	mat44
struct	nifti_brick_list
struct	nifti_image
Defines
#define	DONT_INCLUDE_ANALYZE_STRUCT   /*** not needed herein ***/
#define	NIFTI_L2R   1
#define	NIFTI_R2L   2
#define	NIFTI_P2A   3
#define	NIFTI_A2P   4
#define	NIFTI_I2S   5
#define	NIFTI_S2I   6
#define	LNI_FERR(func, msg, file)   fprintf(stderr,"** ERROR (%s): %s '%s'\n",func,msg,file)
#define	swap_2(s)   swap_2bytes(1,&(s))
#define	swap_4(v)   swap_4bytes(1,&(v))
#define	IS_GOOD_FLOAT(x)   1
#define	FIXED_FLOAT(x)   (x)
#define	ASSIF(p, v)   if( (p)!=NULL ) *(p) = (v)
#define	LSB_FIRST   1
#define	MSB_FIRST   2
#define	REVERSE_ORDER(x)   (3-(x))
#define	LNI_MAX_NIA_EXT_LEN   100000
Functions
char *	nifti_datatype_string (int dt)
char *	nifti_units_string (int uu)
char *	nifti_intent_string (int ii)
char *	nifti_xform_string (int xx)
char *	nifti_slice_string (int ss)
char *	nifti_orientation_string (int ii)
int	nifti_is_inttype (int dt)
mat44	mat44_inverse (mat44 R)
mat33	mat33_inverse (mat33 R)
mat33	mat33_polar (mat33 A)
float	mat33_rownorm (mat33 A)
float	mat33_colnorm (mat33 A)
float	mat33_determ (mat33 R)
mat33	mat33_mul (mat33 A, mat33 B)
void	swap_2bytes (int n, void *ar)
void	swap_4bytes (int n, void *ar)
void	swap_8bytes (int n, void *ar)
void	swap_16bytes (int n, void *ar)
void	swap_Nbytes (int n, int siz, void *ar)
void	swap_nifti_header (struct nifti_1_header *h, int is_nifti)
unsigned int	get_filesize (char *pathname)
nifti_image *	nifti_image_read_bricks (char *hname, int nbricks, int *blist, nifti_brick_list *NBL)
int	nifti_image_load_bricks (nifti_image *nim, int nbricks, int *blist, nifti_brick_list *NBL)
void	free_NBL (nifti_brick_list *NBL)
nifti_image *	nifti_image_read (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)
void	nifti_image_write (nifti_image *nim)
void	nifti_image_infodump (nifti_image *nim)
void	nifti_disp_lib_hist (void)
void	nifti_disp_lib_version (void)
int	nifti_disp_matrix_orient (char *mesg, mat44 mat)
char *	nifti_image_to_ascii (nifti_image *nim)
nifti_image *	nifti_image_from_ascii (char *str, int *bytes_read)
size_t	nifti_get_volsize (nifti_image *nim)
int	is_nifti_file (char *hname)
char *	nifti_find_file_extension (char *name)
int	nifti_validfilename (char *fname)
int	disp_nifti_1_header (char *info, nifti_1_header *hp)
void	nifti_set_debug_level (int level)
int	print_hex_vals (char *data, int nbytes, FILE *fp)
int	valid_nifti_brick_list (nifti_image *nim, int nbricks, int *blist, int disp_error)
void	nifti_datatype_sizes (int datatype, int *nbyper, int *swapsize)
void	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	quatern_to_mat44 (float qb, float qc, float qd, float qx, float qy, float qz, float dx, float dy, float dz, float qfac)
mat44	make_orthog_mat44 (float r11, float r12, float r13, float r21, float r22, float r23, float r31, float r32, float r33)
int	unescape_string (char *str)
char *	escapize_string (char *str)
int	short_order (void)
void	mat44_to_orientation (mat44 R, int *icod, int *jcod, int *kcod)
char *	nifti_findhdrname (char *fname)
char *	nifti_findimgname (char *fname, int nifti_type)
int	nifti_is_gzfile (char *fname)
char *	nifti_makebasename (char *fname)
int	nifti_read_extensions (nifti_image *nim, znzFile fp)
int	nifti_read_next_extension (nifti1_extension *nex, nifti_image *nim, znzFile fp)
int	nifti_add_exten_to_list (nifti1_extension *new_ext, nifti1_extension **list, int new_length)
int	nifti_write_extensions (znzFile fp, nifti_image *nim)
int	nifti_valid_extension (nifti_image *nim, int size, int code)
void	update_nifti_image_for_brick_list (nifti_image *nim, int nbricks)
int	nifti_load_NBL_bricks (nifti_image *nim, int *slist, int *sindex, nifti_brick_list *NBL, znzFile fp)
int	nifti_alloc_NBL_mem (nifti_image *nim, int nbricks, nifti_brick_list *nbl)
int	nifti_copynsort (int nbricks, int *blist, int **slist, int **sindex)
int	int_force_positive (int *list, int nel)
znzFile	nifti_image_write_hdr_img (nifti_image *nim, int write_data, char *opts)
znzFile	nifti_image_write_hdr_img2 (nifti_image *nim, int write_data, char *opts, znzFile *imgfile)
znzFile	nifti_image_open (char *hname, char *opts, nifti_image **nim)
znzFile	nifti_image_load_prep (nifti_image *nim)
size_t	nifti_read_buffer (znzFile fp, void *datatptr, size_t ntot, nifti_image *nim)
void	nifti_write_all_data (znzFile fp, nifti_image *nim)
size_t	nifti_write_buffer (znzFile fp, void *buffer, size_t numbytes)
nifti_image *	nifti_copy_nim_info (nifti_image *src)
nifti_image *	nifti_simple_init_nim ()
nifti_1_header	nifti_convert_nim2nhdr (nifti_image *nim)
nifti_image *	nifti_convert_nhdr2nim (struct nifti_1_header nhdr, char *fname)
void	nifti_set_iname_offset (nifti_image *nim)

---

Define Documentation

#define ASSIF (  p, v    )     if( (p)!=NULL ) *(p) = (v)

Definition at line 307 of file nifti1_io.h. Referenced by mat44_to_quatern(), MCW_widget_geom(), and nifti_datatype_sizes().

#define DONT_INCLUDE_ANALYZE_STRUCT   /*** not needed herein ***/

Definition at line 11 of file nifti1_io.h.

#define FIXED_FLOAT (  x    )     (x)

Definition at line 303 of file nifti1_io.h. Referenced by nifti_convert_nhdr2nim().

#define IS_GOOD_FLOAT (  x    )     1

Definition at line 302 of file nifti1_io.h. Referenced by nifti_convert_nhdr2nim(), nifti_read_buffer(), thd_complexscan(), and thd_floatscan().

#define LNI_FERR (  func, msg, file    )     fprintf(stderr,"** ERROR (%s): %s '%s'\n",func,msg,file)

Definition at line 284 of file nifti1_io.h. Referenced by nifti_image_load_prep(), and nifti_image_read().

#define LNI_MAX_NIA_EXT_LEN   100000

Definition at line 316 of file nifti1_io.h. Referenced by nifti_valid_extension().

#define LSB_FIRST   1

Definition at line 312 of file nifti1_io.h.

#define MSB_FIRST   2

Definition at line 313 of file nifti1_io.h.

#define NIFTI_A2P   4

Definition at line 223 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define NIFTI_I2S   5

Definition at line 224 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define NIFTI_L2R   1

Definition at line 220 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define NIFTI_P2A   3

Definition at line 222 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define NIFTI_R2L   2

Definition at line 221 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define NIFTI_S2I   6

Definition at line 225 of file nifti1_io.h. Referenced by mat44_to_orientation(), and nifti_orientation_string().

#define REVERSE_ORDER (  x    )     (3-(x))

Definition at line 314 of file nifti1_io.h.

#define swap_2 (  s    )     swap_2bytes(1,&(s))

Definition at line 289 of file nifti1_io.h. Referenced by ART_send_volume(), disp_numeric_data(), ge4_swap_all_bytes(), main(), mri_read3D_analyze75(), mri_read_analyze75(), nifti_convert_nhdr2nim(), swap_analyze_hdr(), swap_nifti_header(), THD_open_analyze(), THD_open_ctfmri(), and write_data_to_file().

#define swap_4 (  v    )     swap_4bytes(1,&(v))

Definition at line 290 of file nifti1_io.h. Referenced by disp_numeric_data(), ge4_swap_all_bytes(), ge_header(), is_nifti_file(), main(), mri_imcount_siemens(), mri_read(), mri_read_siemens(), nifti_convert_nhdr2nim(), read_ge_header(), read_ge_image(), read_socket(), swap_analyze_hdr(), swap_nifti_header(), THD_open_ctfmri(), THD_open_ctfsam(), and write_data_to_file().

---

Function Documentation

int disp_nifti_1_header (  char *    info, nifti_1_header *    hp )

Definition at line 1947 of file nifti1_io.c. References nifti_1_header::aux_file, nifti_1_header::bitpix, c, nifti_1_header::cal_max, nifti_1_header::cal_min, nifti_1_header::data_type, nifti_1_header::datatype, nifti_1_header::db_name, nifti_1_header::descrip, nifti_1_header::dim, nifti_1_header::dim_info, nifti_1_header::extents, nifti_1_header::glmax, nifti_1_header::glmin, nifti_1_header::intent_code, nifti_1_header::intent_name, nifti_1_header::intent_p1, nifti_1_header::intent_p2, nifti_1_header::intent_p3, nifti_1_header::magic, nifti_1_header::pixdim, print_hex_vals(), nifti_1_header::qform_code, nifti_1_header::qoffset_x, nifti_1_header::qoffset_y, nifti_1_header::qoffset_z, nifti_1_header::quatern_b, nifti_1_header::quatern_c, nifti_1_header::quatern_d, nifti_1_header::regular, nifti_1_header::scl_inter, nifti_1_header::scl_slope, nifti_1_header::session_error, nifti_1_header::sform_code, nifti_1_header::sizeof_hdr, nifti_1_header::slice_code, nifti_1_header::slice_duration, nifti_1_header::slice_end, nifti_1_header::slice_start, nifti_1_header::srow_x, nifti_1_header::srow_y, nifti_1_header::srow_z, nifti_1_header::toffset, nifti_1_header::vox_offset, and nifti_1_header::xyzt_units. Referenced by nifti_convert_nhdr2nim(). { int c; fputs( "-------------------------------------------------------\n", stderr ); if ( info ) fputs( info, stderr ); if ( !hp ){ fputs(" ** no nifti_1_header to display!\n",stderr); return 1; } fprintf(stderr," nifti_1_header :\n" " sizeof_hdr = %d\n" " data_type[10] = ", hp->sizeof_hdr); print_hex_vals(hp->data_type, 10, stderr); fprintf(stderr, "\n" " db_name[18] = "); print_hex_vals(hp->db_name, 18, stderr); fprintf(stderr, "\n" " extents = %d\n" " session_error = %d\n" " regular = 0x%x\n" " dim_info = 0x%x\n", hp->extents, hp->session_error, hp->regular, hp->dim_info ); fprintf(stderr, " dim[8] ="); for ( c = 0; c < 8; c++ ) fprintf(stderr," %d", hp->dim[c]); fprintf(stderr, "\n" " intent_p1 = %f\n" " intent_p2 = %f\n" " intent_p3 = %f\n" " intent_code = %d\n" " datatype = %d\n" " bitpix = %d\n" " slice_start = %d\n" " pixdim[8] =", hp->intent_p1, hp->intent_p2, hp->intent_p3, hp->intent_code, hp->datatype, hp->bitpix, hp->slice_start); /* break pixdim over 2 lines */ for ( c = 0; c < 4; c++ ) fprintf(stderr," %f", hp->pixdim[c]); fprintf(stderr, "\n "); for ( c = 4; c < 8; c++ ) fprintf(stderr," %f", hp->pixdim[c]); fprintf(stderr, "\n" " vox_offset = %f\n" " scl_slope = %f\n" " scl_inter = %f\n" " slice_end = %d\n" " slice_code = %d\n" " xyzt_units = 0x%x\n" " cal_max = %f\n" " cal_min = %f\n" " slice_duration = %f\n" " toffset = %f\n" " glmax = %d\n" " glmin = %d\n", hp->vox_offset, hp->scl_slope, hp->scl_inter, hp->slice_end, hp->slice_code, hp->xyzt_units, hp->cal_max, hp->cal_min, hp->slice_duration, hp->toffset, hp->glmax, hp->glmin); fprintf(stderr, " descrip = '%s'\n" " aux_file = '%s'\n" " qform_code = %d\n" " sform_code = %d\n" " quatern_b = %f\n" " quatern_c = %f\n" " quatern_d = %f\n" " qoffset_x = %f\n" " qoffset_y = %f\n" " qoffset_z = %f\n" " srow_x[4] = %f, %f, %f, %f\n" " srow_y[4] = %f, %f, %f, %f\n" " srow_z[4] = %f, %f, %f, %f\n" " intent_name = %-16s\n" " magic = %-4s\n", hp->descrip, hp->aux_file, hp->qform_code, hp->sform_code, hp->quatern_b, hp->quatern_c, hp->quatern_d, hp->qoffset_x, hp->qoffset_y, hp->qoffset_z, hp->srow_x[0], hp->srow_x[1], hp->srow_x[2], hp->srow_x[3], hp->srow_y[0], hp->srow_y[1], hp->srow_y[2], hp->srow_y[3], hp->srow_z[0], hp->srow_z[1], hp->srow_z[2], hp->srow_z[3], hp->intent_name, hp->magic); fputs( "-------------------------------------------------------\n", stderr ); return 0; }

char* escapize_string (  char *    str )

Definition at line 3573 of file nifti1_io.c. References calloc, CR, LF, and nifti_strdup(). Referenced by nifti_image_to_ascii(). { int ii,jj , lstr,lout ; char *out ; if( str == NULL || (lstr=strlen(str)) == 0 ){ /* 0 length */ out = nifti_strdup("''") ; return out ; /* string?? */ } lout = 4 ; /* initialize length of output */ for( ii=0 ; ii < lstr ; ii++ ){ /* count characters for output */ switch( str[ii] ){ case '&': lout += 5 ; break ; /* replace '&' with "&amp;" */ case '<': case '>': lout += 4 ; break ; /* replace '<' with "&lt;" */ case '"' : case '\'': lout += 6 ; break ; /* replace '"' with "&quot;" */ case CR: case LF: lout += 6 ; break ; /* replace CR with "&#x0d;" LF with "&#x0a;" */ default: lout++ ; break ; /* copy all other chars */ } } out = (char *)calloc(1,lout) ; /* allocate output string */ out[0] = '\'' ; /* opening quote mark */ for( ii=0,jj=1 ; ii < lstr ; ii++ ){ switch( str[ii] ){ default: out[jj++] = str[ii] ; break ; /* normal characters */ case '&': memcpy(out+jj,"&amp;",5) ; jj+=5 ; break ; case '<': memcpy(out+jj,"&lt;",4) ; jj+=4 ; break ; case '>': memcpy(out+jj,"&gt;",4) ; jj+=4 ; break ; case '"' : memcpy(out+jj,"&quot;",6) ; jj+=6 ; break ; case '\'': memcpy(out+jj,"&apos;",6) ; jj+=6 ; break ; case CR: memcpy(out+jj,"&#x0d;",6) ; jj+=6 ; break ; case LF: memcpy(out+jj,"&#x0a;",6) ; jj+=6 ; break ; } } out[jj++] = '\'' ; /* closing quote mark */ out[jj] = '\0' ; /* terminate the string */ return out ; }

void free_NBL (  nifti_brick_list *    NBL )

Definition at line 317 of file nifti1_io.c. References nifti_brick_list::bricks, nifti_brick_list::bsize, c, free, and nifti_brick_list::nbricks. Referenced by nifti_image_load_bricks(). { int c; if( NBL->bricks ){ for( c = 0; c < NBL->nbricks; c++ ) if( NBL->bricks[c] ) free(NBL->bricks[c]); free(NBL->bricks); NBL->bricks = NULL; } NBL->nbricks = NBL->bsize = 0; }

unsigned int get_filesize (  char *    pathname )

Definition at line 1628 of file nifti1_io.c. Referenced by nifti_image_load_prep(), and nifti_image_read(). { struct stat buf ; int ii ; if( pathname == NULL || *pathname == '\0' ) return 0 ; ii = stat( pathname , &buf ); if( ii != 0 ) return 0 ; return (unsigned int)buf.st_size ; }

int int_force_positive (  int *    list, int    nel )  [static]

Referenced by valid_nifti_brick_list().

int is_nifti_file (  char *    hname )

if using compression libraries, .gz files are looked for first * Definition at line 1886 of file nifti1_io.c. References free, gni_debug, nifti_findhdrname(), NIFTI_ONEFILE, nifti_validfilename(), NIFTI_VERSION, nifti_1_header::sizeof_hdr, and swap_4. { struct nifti_1_header nhdr ; znzFile fp ; int ii ; char *tmpname; /* bad input name? */ if( !nifti_validfilename(hname) ) return -1 ; /* open file */ tmpname = nifti_findhdrname(hname); if( tmpname == NULL ){ if( gni_debug > 0 ) fprintf(stderr,"** no header file found for '%s'\n",hname); return -1; } fp = znzopen( tmpname , "rb" , 1 ) ; free(tmpname); if (znz_isnull(fp)) return -1 ; /* bad open? */ /* read header, close file */ ii = znzread( &nhdr , 1 , sizeof(nhdr) , fp ) ; znzclose( fp ) ; if( ii < sizeof(nhdr) ) return -1 ; /* bad read? */ /* check for NIFTI-ness */ if( NIFTI_VERSION(nhdr) != 0 ){ return ( NIFTI_ONEFILE(nhdr) ) ? 1 : 2 ; } /* check for ANALYZE-ness (sizeof_hdr field == 348) */ ii = nhdr.sizeof_hdr ; if( ii == sizeof(nhdr) ) return 0 ; /* matches */ /* try byte-swapping header */ swap_4(ii) ; if( ii == sizeof(nhdr) ) return 0 ; /* matches */ return -1 ; /* not good */ }

mat44 make_orthog_mat44 (  float    r11, float    r12, float    r13, float    r21, float    r22, float    r23, float    r31, float    r32, float    r33 )

INPUT MATRIX IS: * Definition at line 1104 of file nifti1_io.c. References mat33::m, mat44::m, and mat33_polar(). { mat44 R ; mat33 Q , P ; double val ; R.m[3][0] = R.m[3][1] = R.m[3][2] = 0.0l ; R.m[3][3] = 1.0l ; Q.m[0][0] = r11 ; Q.m[0][1] = r12 ; Q.m[0][2] = r13 ; /* load Q */ Q.m[1][0] = r21 ; Q.m[1][1] = r22 ; Q.m[1][2] = r23 ; Q.m[2][0] = r31 ; Q.m[2][1] = r32 ; Q.m[2][2] = r33 ; /* normalize row 1 */ val = Q.m[0][0]*Q.m[0][0] + Q.m[0][1]*Q.m[0][1] + Q.m[0][2]*Q.m[0][2] ; if( val > 0.0l ){ val = 1.0l / sqrt(val) ; Q.m[0][0] *= val ; Q.m[0][1] *= val ; Q.m[0][2] *= val ; } else { Q.m[0][0] = 1.0l ; Q.m[0][1] = 0.0l ; Q.m[0][2] = 0.0l ; } /* normalize row 2 */ val = Q.m[1][0]*Q.m[1][0] + Q.m[1][1]*Q.m[1][1] + Q.m[1][2]*Q.m[1][2] ; if( val > 0.0l ){ val = 1.0l / sqrt(val) ; Q.m[1][0] *= val ; Q.m[1][1] *= val ; Q.m[1][2] *= val ; } else { Q.m[1][0] = 0.0l ; Q.m[1][1] = 1.0l ; Q.m[1][2] = 0.0l ; } /* normalize row 3 */ val = Q.m[2][0]*Q.m[2][0] + Q.m[2][1]*Q.m[2][1] + Q.m[2][2]*Q.m[2][2] ; if( val > 0.0l ){ val = 1.0l / sqrt(val) ; Q.m[2][0] *= val ; Q.m[2][1] *= val ; Q.m[2][2] *= val ; } else { Q.m[2][0] = Q.m[0][1]*Q.m[1][2] - Q.m[0][2]*Q.m[1][1] ; /* cross */ Q.m[2][1] = Q.m[0][2]*Q.m[1][0] - Q.m[0][0]*Q.m[1][2] ; /* product */ Q.m[2][2] = Q.m[0][0]*Q.m[1][1] - Q.m[0][1]*Q.m[1][0] ; } P = mat33_polar(Q) ; /* P is orthog matrix closest to Q */ R.m[0][0] = P.m[0][0] ; R.m[0][1] = P.m[0][1] ; R.m[0][2] = P.m[0][2] ; R.m[1][0] = P.m[1][0] ; R.m[1][1] = P.m[1][1] ; R.m[1][2] = P.m[1][2] ; R.m[2][0] = P.m[2][0] ; R.m[2][1] = P.m[2][1] ; R.m[2][2] = P.m[2][2] ; R.m[0][3] = R.m[1][3] = R.m[2][3] = 0.0 ; return R ; }

float mat33_colnorm (  mat33    A )

Definition at line 1220 of file nifti1_io.c. References mat33::m. Referenced by mat33_polar(). { float r1,r2,r3 ; r1 = fabs(A.m[0][0])+fabs(A.m[1][0])+fabs(A.m[2][0]) ; r2 = fabs(A.m[0][1])+fabs(A.m[1][1])+fabs(A.m[2][1]) ; r3 = fabs(A.m[0][2])+fabs(A.m[1][2])+fabs(A.m[2][2]) ; if( r1 < r2 ) r1 = r2 ; if( r1 < r3 ) r1 = r3 ; return r1 ; }

float mat33_determ (  mat33    R )

INPUT MATRIX: * Definition at line 1192 of file nifti1_io.c. References mat33::m. Referenced by mat33_polar(), and mat44_to_orientation(). { double r11,r12,r13,r21,r22,r23,r31,r32,r33 ; /** INPUT MATRIX: **/ r11 = R.m[0][0]; r12 = R.m[0][1]; r13 = R.m[0][2]; /* [ r11 r12 r13 ] */ r21 = R.m[1][0]; r22 = R.m[1][1]; r23 = R.m[1][2]; /* [ r21 r22 r23 ] */ r31 = R.m[2][0]; r32 = R.m[2][1]; r33 = R.m[2][2]; /* [ r31 r32 r33 ] */ return r11*r22*r33-r11*r32*r23-r21*r12*r33 +r21*r32*r13+r31*r12*r23-r31*r22*r13 ; }

mat33 mat33_inverse (  mat33    R )

Definition at line 1161 of file nifti1_io.c. References mat33::m. Referenced by mat33_polar(). { double r11,r12,r13,r21,r22,r23,r31,r32,r33 , deti ; mat33 Q ; /** INPUT MATRIX: **/ r11 = R.m[0][0]; r12 = R.m[0][1]; r13 = R.m[0][2]; /* [ r11 r12 r13 ] */ r21 = R.m[1][0]; r22 = R.m[1][1]; r23 = R.m[1][2]; /* [ r21 r22 r23 ] */ r31 = R.m[2][0]; r32 = R.m[2][1]; r33 = R.m[2][2]; /* [ r31 r32 r33 ] */ deti = r11*r22*r33-r11*r32*r23-r21*r12*r33 +r21*r32*r13+r31*r12*r23-r31*r22*r13 ; if( deti != 0.0l ) deti = 1.0l / deti ; Q.m[0][0] = deti*( r22*r33-r32*r23) ; Q.m[0][1] = deti*(-r12*r33+r32*r13) ; Q.m[0][2] = deti*( r12*r23-r22*r13) ; Q.m[1][0] = deti*(-r21*r33+r31*r23) ; Q.m[1][1] = deti*( r11*r33-r31*r13) ; Q.m[1][2] = deti*(-r11*r23+r21*r13) ; Q.m[2][0] = deti*( r21*r32-r31*r22) ; Q.m[2][1] = deti*(-r11*r32+r31*r12) ; Q.m[2][2] = deti*( r11*r22-r21*r12) ; return Q ; }

mat33 mat33_mul (  mat33    A, mat33    B )

Definition at line 1234 of file nifti1_io.c. References i, and mat33::m. Referenced by mat44_to_orientation(). { mat33 C ; int i,j ; for( i=0 ; i < 3 ; i++ ) for( j=0 ; j < 3 ; j++ ) C.m[i][j] = A.m[i][0] * B.m[0][j] + A.m[i][1] * B.m[1][j] + A.m[i][2] * B.m[2][j] ; return C ; }

mat33 mat33_polar (  mat33    A )

Definition at line 1251 of file nifti1_io.c. References mat33::m, mat33_colnorm(), mat33_determ(), mat33_inverse(), and mat33_rownorm(). Referenced by make_orthog_mat44(), and mat44_to_quatern(). { mat33 X , Y , Z ; float alp,bet,gam,gmi , dif=1.0 ; int k=0 ; X = A ; /* force matrix to be nonsingular */ gam = mat33_determ(X) ; while( gam == 0.0 ){ /* perturb matrix */ gam = 0.00001 * ( 0.001 + mat33_rownorm(X) ) ; X.m[0][0] += gam ; X.m[1][1] += gam ; X.m[2][2] += gam ; gam = mat33_determ(X) ; } while(1){ Y = mat33_inverse(X) ; if( dif > 0.3 ){ /* far from convergence */ alp = sqrt( mat33_rownorm(X) * mat33_colnorm(X) ) ; bet = sqrt( mat33_rownorm(Y) * mat33_colnorm(Y) ) ; gam = sqrt( bet / alp ) ; gmi = 1.0 / gam ; } else { gam = gmi = 1.0 ; /* close to convergence */ } Z.m[0][0] = 0.5 * ( gam*X.m[0][0] + gmi*Y.m[0][0] ) ; Z.m[0][1] = 0.5 * ( gam*X.m[0][1] + gmi*Y.m[1][0] ) ; Z.m[0][2] = 0.5 * ( gam*X.m[0][2] + gmi*Y.m[2][0] ) ; Z.m[1][0] = 0.5 * ( gam*X.m[1][0] + gmi*Y.m[0][1] ) ; Z.m[1][1] = 0.5 * ( gam*X.m[1][1] + gmi*Y.m[1][1] ) ; Z.m[1][2] = 0.5 * ( gam*X.m[1][2] + gmi*Y.m[2][1] ) ; Z.m[2][0] = 0.5 * ( gam*X.m[2][0] + gmi*Y.m[0][2] ) ; Z.m[2][1] = 0.5 * ( gam*X.m[2][1] + gmi*Y.m[1][2] ) ; Z.m[2][2] = 0.5 * ( gam*X.m[2][2] + gmi*Y.m[2][2] ) ; dif = fabs(Z.m[0][0]-X.m[0][0])+fabs(Z.m[0][1]-X.m[0][1]) +fabs(Z.m[0][2]-X.m[0][2])+fabs(Z.m[1][0]-X.m[1][0]) +fabs(Z.m[1][1]-X.m[1][1])+fabs(Z.m[1][2]-X.m[1][2]) +fabs(Z.m[2][0]-X.m[2][0])+fabs(Z.m[2][1]-X.m[2][1]) +fabs(Z.m[2][2]-X.m[2][2]) ; k = k+1 ; if( k > 100 || dif < 3.e-6 ) break ; /* convergence or exhaustion */ X = Z ; } return Z ; }

float mat33_rownorm (  mat33    A )

INPUT MATRIX: * Definition at line 1206 of file nifti1_io.c. References mat33::m. Referenced by mat33_polar(). { float r1,r2,r3 ; r1 = fabs(A.m[0][0])+fabs(A.m[0][1])+fabs(A.m[0][2]) ; r2 = fabs(A.m[1][0])+fabs(A.m[1][1])+fabs(A.m[1][2]) ; r3 = fabs(A.m[2][0])+fabs(A.m[2][1])+fabs(A.m[2][2]) ; if( r1 < r2 ) r1 = r2 ; if( r1 < r3 ) r1 = r3 ; return r1 ; }

mat44 mat44_inverse (  mat44    R )

Definition at line 1040 of file nifti1_io.c. References mat44::m, and v1. Referenced by nifti_convert_nhdr2nim(), and nifti_image_from_ascii(). { double r11,r12,r13,r21,r22,r23,r31,r32,r33,v1,v2,v3 , deti ; mat44 Q ; /** INPUT MATRIX IS: **/ r11 = R.m[0][0]; r12 = R.m[0][1]; r13 = R.m[0][2]; /* [ r11 r12 r13 v1 ] */ r21 = R.m[1][0]; r22 = R.m[1][1]; r23 = R.m[1][2]; /* [ r21 r22 r23 v2 ] */ r31 = R.m[2][0]; r32 = R.m[2][1]; r33 = R.m[2][2]; /* [ r31 r32 r33 v3 ] */ v1 = R.m[0][3]; v2 = R.m[1][3]; v3 = R.m[2][3]; /* [ 0 0 0 1 ] */ deti = r11*r22*r33-r11*r32*r23-r21*r12*r33 +r21*r32*r13+r31*r12*r23-r31*r22*r13 ; if( deti != 0.0l ) deti = 1.0l / deti ; Q.m[0][0] = deti*( r22*r33-r32*r23) ; Q.m[0][1] = deti*(-r12*r33+r32*r13) ; Q.m[0][2] = deti*( r12*r23-r22*r13) ; Q.m[0][3] = deti*(-r12*r23*v3+r12*v2*r33+r22*r13*v3 -r22*v1*r33-r32*r13*v2+r32*v1*r23) ; Q.m[1][0] = deti*(-r21*r33+r31*r23) ; Q.m[1][1] = deti*( r11*r33-r31*r13) ; Q.m[1][2] = deti*(-r11*r23+r21*r13) ; Q.m[1][3] = deti*( r11*r23*v3-r11*v2*r33-r21*r13*v3 +r21*v1*r33+r31*r13*v2-r31*v1*r23) ; Q.m[2][0] = deti*( r21*r32-r31*r22) ; Q.m[2][1] = deti*(-r11*r32+r31*r12) ; Q.m[2][2] = deti*( r11*r22-r21*r12) ; Q.m[2][3] = deti*(-r11*r22*v3+r11*r32*v2+r21*r12*v3 -r21*r32*v1-r31*r12*v2+r31*r22*v1) ; Q.m[3][0] = Q.m[3][1] = Q.m[3][2] = 0.0l ; Q.m[3][3] = (deti == 0.0l) ? 0.0l : 1.0l ; /* failure flag if deti == 0 */ return Q ; }

void mat44_to_orientation (  mat44    R, int *    icod, int *    jcod, int *    kcod )

Definition at line 1318 of file nifti1_io.c. References i, mat33::m, mat44::m, mat33_determ(), mat33_mul(), NIFTI_A2P, NIFTI_I2S, NIFTI_L2R, NIFTI_P2A, NIFTI_R2L, NIFTI_S2I, p, and r. Referenced by nifti_disp_matrix_orient(), and nifti_image_to_ascii(). { float xi,xj,xk , yi,yj,yk , zi,zj,zk , val,detQ,detP ; mat33 P , Q , M ; int i,j,k=0,p,q,r , ibest,jbest,kbest,pbest,qbest,rbest ; float vbest ; if( icod == NULL || jcod == NULL || kcod == NULL ) return ; /* bad */ *icod = *jcod = *kcod = 0 ; /* error returns, if sh*t happens */ /* load column vectors for each (i,j,k) direction from matrix */ /*-- i axis --*/ /*-- j axis --*/ /*-- k axis --*/ xi = R.m[0][0] ; xj = R.m[0][1] ; xk = R.m[0][2] ; yi = R.m[1][0] ; yj = R.m[1][1] ; yk = R.m[1][2] ; zi = R.m[2][0] ; zj = R.m[2][1] ; zk = R.m[2][2] ; /* normalize column vectors to get unit vectors along each ijk-axis */ /* normalize i axis */ val = sqrt( xi*xi + yi*yi + zi*zi ) ; if( val == 0.0 ) return ; /* stupid input */ xi /= val ; yi /= val ; zi /= val ; /* normalize j axis */ val = sqrt( xj*xj + yj*yj + zj*zj ) ; if( val == 0.0 ) return ; /* stupid input */ xj /= val ; yj /= val ; zj /= val ; /* orthogonalize j axis to i axis, if needed */ val = xi*xj + yi*yj + zi*zj ; /* dot product between i and j */ if( fabs(val) > 1.e-4 ){ xj -= val*xi ; yj -= val*yi ; zj -= val*zi ; val = sqrt( xj*xj + yj*yj + zj*zj ) ; /* must renormalize */ if( val == 0.0 ) return ; /* j was parallel to i? */ xj /= val ; yj /= val ; zj /= val ; } /* normalize k axis; if it is zero, make it the cross product i x j */ val = sqrt( xk*xk + yk*yk + zk*zk ) ; if( val == 0.0 ){ xk = yi*zj-zi*yj; yk = zi*xj-zj*xi ; zk=xi*yj-yi*xj ; } else { xk /= val ; yk /= val ; zk /= val ; } /* orthogonalize k to i */ val = xi*xk + yi*yk + zi*zk ; /* dot product between i and k */ if( fabs(val) > 1.e-4 ){ xk -= val*xi ; yk -= val*yi ; zk -= val*zi ; val = sqrt( xk*xk + yk*yk + zk*zk ) ; if( val == 0.0 ) return ; /* bad */ xk /= val ; yk /= val ; zk /= val ; } /* orthogonalize k to j */ val = xj*xk + yj*yk + zj*zk ; /* dot product between j and k */ if( fabs(val) > 1.e-4 ){ xk -= val*xj ; yk -= val*yj ; zk -= val*zj ; val = sqrt( xk*xk + yk*yk + zk*zk ) ; if( val == 0.0 ) return ; /* bad */ xk /= val ; yk /= val ; zk /= val ; } Q.m[0][0] = xi ; Q.m[0][1] = xj ; Q.m[0][2] = xk ; Q.m[1][0] = yi ; Q.m[1][1] = yj ; Q.m[1][2] = yk ; Q.m[2][0] = zi ; Q.m[2][1] = zj ; Q.m[2][2] = zk ; /* at this point, Q is the rotation matrix from the (i,j,k) to (x,y,z) axes */ detQ = mat33_determ( Q ) ; if( detQ == 0.0 ) return ; /* shouldn't happen unless user is a DUFIS */ /* Build and test all possible +1/-1 coordinate permutation matrices P; then find the P such that the rotation matrix M=PQ is closest to the identity, in the sense of M having the smallest total rotation angle. */ /* Despite the formidable looking 6 nested loops, there are only 3*3*3*2*2*2 = 216 passes, which will run very quickly. */ vbest = -666.0 ; ibest=pbest=qbest=rbest=1 ; jbest=2 ; kbest=3 ; for( i=1 ; i <= 3 ; i++ ){ /* i = column number to use for row #1 */ for( j=1 ; j <= 3 ; j++ ){ /* j = column number to use for row #2 */ if( i == j ) continue ; for( k=1 ; k <= 3 ; k++ ){ /* k = column number to use for row #3 */ if( i == k || j == k ) continue ; P.m[0][0] = P.m[0][1] = P.m[0][2] = P.m[1][0] = P.m[1][1] = P.m[1][2] = P.m[2][0] = P.m[2][1] = P.m[2][2] = 0.0 ; for( p=-1 ; p <= 1 ; p+=2 ){ /* p,q,r are -1 or +1 */ for( q=-1 ; q <= 1 ; q+=2 ){ /* and go into rows #1,2,3 */ for( r=-1 ; r <= 1 ; r+=2 ){ P.m[0][i-1] = p ; P.m[1][j-1] = q ; P.m[2][k-1] = r ; detP = mat33_determ(P) ; /* sign of permutation */ if( detP * detQ <= 0.0 ) continue ; /* doesn't match sign of Q */ M = mat33_mul(P,Q) ; /* angle of M rotation = 2.0*acos(0.5*sqrt(1.0+trace(M))) */ /* we want largest trace(M) == smallest angle == M nearest to I */ val = M.m[0][0] + M.m[1][1] + M.m[2][2] ; /* trace */ if( val > vbest ){ vbest = val ; ibest = i ; jbest = j ; kbest = k ; pbest = p ; qbest = q ; rbest = r ; } }}}}}} /* At this point ibest is 1 or 2 or 3; pbest is -1 or +1; etc. The matrix P that corresponds is the best permutation approximation to Q-inverse; that is, P (approximately) takes (x,y,z) coordinates to the (i,j,k) axes. For example, the first row of P (which contains pbest in column ibest) determines the way the i axis points relative to the anatomical (x,y,z) axes. If ibest is 2, then the i axis is along the y axis, which is direction P2A (if pbest > 0) or A2P (if pbest < 0). So, using ibest and pbest, we can assign the output code for the i axis. Mutatis mutandis for the j and k axes, of course. */ switch( ibest*pbest ){ case 1: i = NIFTI_L2R ; break ; case -1: i = NIFTI_R2L ; break ; case 2: i = NIFTI_P2A ; break ; case -2: i = NIFTI_A2P ; break ; case 3: i = NIFTI_I2S ; break ; case -3: i = NIFTI_S2I ; break ; } switch( jbest*qbest ){ case 1: j = NIFTI_L2R ; break ; case -1: j = NIFTI_R2L ; break ; case 2: j = NIFTI_P2A ; break ; case -2: j = NIFTI_A2P ; break ; case 3: j = NIFTI_I2S ; break ; case -3: j = NIFTI_S2I ; break ; } switch( kbest*rbest ){ case 1: k = NIFTI_L2R ; break ; case -1: k = NIFTI_R2L ; break ; case 2: k = NIFTI_P2A ; break ; case -2: k = NIFTI_A2P ; break ; case 3: k = NIFTI_I2S ; break ; case -3: k = NIFTI_S2I ; break ; } *icod = i ; *jcod = j ; *kcod = k ; return ; }

void 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 )

Definition at line 917 of file nifti1_io.c. References a, ASSIF, c, mat33::m, mat44::m, and mat33_polar(). { double r11,r12,r13 , r21,r22,r23 , r31,r32,r33 ; double xd,yd,zd , a,b,c,d ; mat33 P,Q ; /* offset outputs are read write out of input matrix */ ASSIF(qx,R.m[0][3]) ; ASSIF(qy,R.m[1][3]) ; ASSIF(qz,R.m[2][3]) ; /* load 3x3 matrix into local variables */ r11 = R.m[0][0] ; r12 = R.m[0][1] ; r13 = R.m[0][2] ; r21 = R.m[1][0] ; r22 = R.m[1][1] ; r23 = R.m[1][2] ; r31 = R.m[2][0] ; r32 = R.m[2][1] ; r33 = R.m[2][2] ; /* compute lengths of each column; these determine grid spacings */ xd = sqrt( r11*r11 + r21*r21 + r31*r31 ) ; yd = sqrt( r12*r12 + r22*r22 + r32*r32 ) ; zd = sqrt( r13*r13 + r23*r23 + r33*r33 ) ; /* if a column length is zero, patch the trouble */ if( xd == 0.0l ){ r11 = 1.0l ; r21 = r31 = 0.0l ; xd = 1.0l ; } if( yd == 0.0l ){ r22 = 1.0l ; r12 = r32 = 0.0l ; yd = 1.0l ; } if( zd == 0.0l ){ r33 = 1.0l ; r13 = r23 = 0.0l ; zd = 1.0l ; } /* assign the output lengths */ ASSIF(dx,xd) ; ASSIF(dy,yd) ; ASSIF(dz,zd) ; /* normalize the columns */ r11 /= xd ; r21 /= xd ; r31 /= xd ; r12 /= yd ; r22 /= yd ; r32 /= yd ; r13 /= zd ; r23 /= zd ; r33 /= zd ; /* At this point, the matrix has normal columns, but we have to allow for the fact that the hideous user may not have given us a matrix with orthogonal columns. So, now find the orthogonal matrix closest to the current matrix. One reason for using the polar decomposition to get this orthogonal matrix, rather than just directly orthogonalizing the columns, is so that inputting the inverse matrix to R will result in the inverse orthogonal matrix at this point. If we just orthogonalized the columns, this wouldn't necessarily hold. */ Q.m[0][0] = r11 ; Q.m[0][1] = r12 ; Q.m[0][2] = r13 ; /* load Q */ Q.m[1][0] = r21 ; Q.m[1][1] = r22 ; Q.m[1][2] = r23 ; Q.m[2][0] = r31 ; Q.m[2][1] = r32 ; Q.m[2][2] = r33 ; P = mat33_polar(Q) ; /* P is orthog matrix closest to Q */ r11 = P.m[0][0] ; r12 = P.m[0][1] ; r13 = P.m[0][2] ; /* unload */ r21 = P.m[1][0] ; r22 = P.m[1][1] ; r23 = P.m[1][2] ; r31 = P.m[2][0] ; r32 = P.m[2][1] ; r33 = P.m[2][2] ; /* [ r11 r12 r13 ] */ /* at this point, the matrix [ r21 r22 r23 ] is orthogonal */ /* [ r31 r32 r33 ] */ /* compute the determinant to determine if it is proper */ zd = r11*r22*r33-r11*r32*r23-r21*r12*r33 +r21*r32*r13+r31*r12*r23-r31*r22*r13 ; /* should be -1 or 1 */ if( zd > 0 ){ /* proper */ ASSIF(qfac,1.0) ; } else { /* improper ==> flip 3rd column */ ASSIF(qfac,-1.0) ; r13 = -r13 ; r23 = -r23 ; r33 = -r33 ; } /* now, compute quaternion parameters */ a = r11 + r22 + r33 + 1.0l ; if( a > 0.5l ){ /* simplest case */ a = 0.5l * sqrt(a) ; b = 0.25l * (r32-r23) / a ; c = 0.25l * (r13-r31) / a ; d = 0.25l * (r21-r12) / a ; } else { /* trickier case */ xd = 1.0 + r11 - (r22+r33) ; /* 4*b*b */ yd = 1.0 + r22 - (r11+r33) ; /* 4*c*c */ zd = 1.0 + r33 - (r11+r22) ; /* 4*d*d */ if( xd > 1.0 ){ b = 0.5l * sqrt(xd) ; c = 0.25l* (r12+r21) / b ; d = 0.25l* (r13+r31) / b ; a = 0.25l* (r32-r23) / b ; } else if( yd > 1.0 ){ c = 0.5l * sqrt(yd) ; b = 0.25l* (r12+r21) / c ; d = 0.25l* (r23+r32) / c ; a = 0.25l* (r13-r31) / c ; } else { d = 0.5l * sqrt(zd) ; b = 0.25l* (r13+r31) / d ; c = 0.25l* (r23+r32) / d ; a = 0.25l* (r21-r12) / d ; } if( a < 0.0l ){ b=-b ; c=-c ; d=-d; a=-a; } } ASSIF(qb,b) ; ASSIF(qc,c) ; ASSIF(qd,d) ; return ; }

int nifti_add_exten_to_list (  nifti1_extension *    new_ext, nifti1_extension **    list, int    new_length )  [static]

Referenced by nifti_read_extensions().

int nifti_alloc_NBL_mem (  nifti_image *    nim, int    nbricks, nifti_brick_list *    nbl )  [static]

Referenced by nifti_image_load_bricks().

nifti_image* nifti_convert_nhdr2nim (  struct nifti_1_header    nhdr, char *    fname )

Definition at line 2035 of file nifti1_io.c. References nifti_1_header::aux_file, nifti_image::aux_file, nifti_image::byteorder, nifti_1_header::cal_max, nifti_image::cal_max, nifti_1_header::cal_min, nifti_image::cal_min, calloc, nifti_image::datatype, nifti_1_header::datatype, nifti_1_header::descrip, nifti_image::descrip, nifti_image::dim, nifti_1_header::dim, nifti_1_header::dim_info, DIM_INFO_TO_FREQ_DIM, DIM_INFO_TO_PHASE_DIM, DIM_INFO_TO_SLICE_DIM, disp_nifti_1_header(), nifti_image::dt, DT_BINARY, DT_UNKNOWN, nifti_image::du, nifti_image::dv, nifti_image::dw, nifti_image::dx, nifti_image::dy, nifti_image::dz, ERREX, nifti_image::ext_list, FIXED_FLOAT, nifti_image::fname, free, nifti_image::freq_dim, gni_debug, nifti_image::iname, nifti_image::iname_offset, nifti_1_header::intent_code, nifti_image::intent_code, nifti_1_header::intent_name, nifti_image::intent_name, nifti_1_header::intent_p1, nifti_image::intent_p1, nifti_1_header::intent_p2, nifti_image::intent_p2, nifti_1_header::intent_p3, nifti_image::intent_p3, IS_GOOD_FLOAT, mat44::m, mat44_inverse(), nifti_image::nbyper, nifti_image::ndim, nifti_datatype_sizes(), nifti_disp_matrix_orient(), nifti_findimgname(), NIFTI_ONEFILE, nifti_strdup(), nifti_image::nifti_type, NIFTI_VERSION, NIFTI_XFORM_UNKNOWN, nifti_image::nt, nifti_image::nu, nifti_image::num_ext, nifti_image::nv, nifti_image::nvox, nifti_image::nw, nifti_image::nx, nifti_image::ny, nifti_image::nz, nifti_image::phase_dim, nifti_image::pixdim, nifti_1_header::pixdim, nifti_image::qfac, nifti_image::qform_code, nifti_1_header::qform_code, nifti_1_header::qoffset_x, nifti_image::qoffset_x, nifti_1_header::qoffset_y, nifti_image::qoffset_y, nifti_1_header::qoffset_z, nifti_image::qoffset_z, nifti_image::qto_ijk, nifti_image::qto_xyz, nifti_1_header::quatern_b, nifti_image::quatern_b, nifti_1_header::quatern_c, nifti_image::quatern_c, nifti_1_header::quatern_d, nifti_image::quatern_d, quatern_to_mat44(), REVERSE_ORDER, nifti_1_header::scl_inter, nifti_image::scl_inter, nifti_1_header::scl_slope, nifti_image::scl_slope, nifti_image::sform_code, nifti_1_header::sform_code, short_order(), nifti_1_header::sizeof_hdr, nifti_1_header::slice_code, nifti_image::slice_code, nifti_image::slice_dim, nifti_1_header::slice_duration, nifti_image::slice_duration, nifti_1_header::slice_end, nifti_image::slice_end, nifti_1_header::slice_start, nifti_image::slice_start, nifti_1_header::srow_x, nifti_1_header::srow_y, nifti_1_header::srow_z, nifti_image::sto_ijk, nifti_image::sto_xyz, swap_2, swap_4, swap_nifti_header(), nifti_image::swapsize, nifti_image::time_units, nifti_1_header::toffset, nifti_image::toffset, nifti_1_header::vox_offset, nifti_image::xyz_units, XYZT_TO_SPACE, XYZT_TO_TIME, and nifti_1_header::xyzt_units. Referenced by nifti_copy_nim_info(), nifti_image_read(), and nifti_simple_init_nim(). { int ii , doswap , ioff, ndim, nvox ; int is_nifti , is_onefile ; short ss ; char *iname=NULL; nifti_image *nim; nim = (nifti_image *) calloc( 1 , sizeof(nifti_image) ) ; if( !nim ) ERREX("failed to allocate nifti image"); nim->fname=NULL; nim->iname=NULL; /** check if have to swap bytes **/ doswap = 0 ; /* swap data flag */ ss = nhdr.dim[0] ; if( ss != 0 ){ /* check dim[0] for good value */ if( ss < 0 || ss > 7 ){ swap_2(ss) ; if( ss < 0 || ss > 7 ) ERREX("bad dim[0]") ; doswap = 1 ; } } else { /* dim[0] == 0 is illegal, but does occur */ ii = nhdr.sizeof_hdr ; /* so check sizeof_hdr field instead */ if( ii != sizeof(nhdr) ){ swap_4(ii) ; if( ii != sizeof(nhdr) ) ERREX("bad sizeof_hdr") ; doswap = 1 ; } } /** determine if this is a NIFTI-1 compliant header **/ is_nifti = NIFTI_VERSION(nhdr) ; if( doswap ) { if ( gni_debug > 2 ) disp_nifti_1_header("-d ni1 pre-swap: ", &nhdr); swap_nifti_header( &nhdr , is_nifti ) ; } if ( gni_debug > 1 ) disp_nifti_1_header("-d nhdr2nim : ", &nhdr); if( nhdr.datatype == DT_BINARY || nhdr.datatype == DT_UNKNOWN ) ERREX("bad datatype") ; if( nhdr.dim[1] <= 0 ) ERREX("bad dim[1]") ; for( ii=2 ; ii <= 7 ; ii++ ) if( nhdr.dim[ii] <= 0 ) nhdr.dim[ii] = 1 ; /* fix bad dim[] values */ /** get number of dimensions (ignoring dim[0] now) **/ for( ii=7 ; ii >= 2 ; ii-- ) /* loop backwards until we */ if( nhdr.dim[ii] > 1 ) break ; /* find a dim bigger than 1 */ ndim = ii ; /* set bad grid spacings to 1.0 */ for( ii=1 ; ii <= 7 ; ii++ ){ if( nhdr.pixdim[ii] == 0.0 || !IS_GOOD_FLOAT(nhdr.pixdim[ii]) ) nhdr.pixdim[ii] = 1.0 ; } is_onefile = is_nifti && NIFTI_ONEFILE(nhdr) ; if( is_nifti ) nim->nifti_type = (is_onefile) ? 1 : 2 ; else nim->nifti_type = 0 ; ii = short_order() ; if( doswap ) nim->byteorder = REVERSE_ORDER(ii) ; else nim->byteorder = ii ; /** dimensions of data array **/ nim->ndim = nim->dim[0] = ndim ; nim->nx = nim->dim[1] = nhdr.dim[1]; nvox = nim->nx; nim->ny = nim->dim[2] = nhdr.dim[2]; nvox *= nim->ny; nim->nz = nim->dim[3] = nhdr.dim[3]; nvox *= nim->nz; nim->nt = nim->dim[4] = nhdr.dim[4]; nvox *= nim->nt; nim->nu = nim->dim[5] = nhdr.dim[5]; nvox *= nim->nu; nim->nv = nim->dim[6] = nhdr.dim[6]; nvox *= nim->nv; nim->nw = nim->dim[7] = nhdr.dim[7]; nvox *= nim->nw; nim->nvox = nvox; /** type of data in voxels and how many bytes per voxel */ nim->datatype = nhdr.datatype ; nifti_datatype_sizes( nim->datatype , &(nim->nbyper) , &(nim->swapsize) ) ; if( nim->nbyper == 0 ){ free(nim); ERREX("bad datatype"); } /** grid spacings **/ nim->dx = nim->pixdim[1] = nhdr.pixdim[1] ; nim->dy = nim->pixdim[2] = nhdr.pixdim[2] ; nim->dz = nim->pixdim[3] = nhdr.pixdim[3] ; nim->dt = nim->pixdim[4] = nhdr.pixdim[4] ; nim->du = nim->pixdim[5] = nhdr.pixdim[5] ; nim->dv = nim->pixdim[6] = nhdr.pixdim[6] ; nim->dw = nim->pixdim[7] = nhdr.pixdim[7] ; /** compute qto_xyz transformation from pixel indexes (i,j,k) to (x,y,z) **/ if( !is_nifti || nhdr.qform_code <= 0 ){ /** default transformation **/ nim->qto_xyz.m[0][0] = nim->dx ; /* grid spacings */ nim->qto_xyz.m[1][1] = nim->dy ; /* along diagonal */ nim->qto_xyz.m[2][2] = nim->dz ; /* off diagonal is zero */ nim->qto_xyz.m[0][1]=nim->qto_xyz.m[0][2]=nim->qto_xyz.m[0][3] = 0.0; nim->qto_xyz.m[1][0]=nim->qto_xyz.m[1][2]=nim->qto_xyz.m[1][3] = 0.0; nim->qto_xyz.m[2][0]=nim->qto_xyz.m[2][1]=nim->qto_xyz.m[2][3] = 0.0; /* last row is always [ 0 0 0 1 ] */ nim->qto_xyz.m[3][0]=nim->qto_xyz.m[3][1]=nim->qto_xyz.m[3][2] = 0.0; nim->qto_xyz.m[3][3]= 1.0 ; nim->qform_code = NIFTI_XFORM_UNKNOWN ; } else { /** NIFTI: quaternion-specified transformation **/ nim->quatern_b = FIXED_FLOAT( nhdr.quatern_b ) ; nim->quatern_c = FIXED_FLOAT( nhdr.quatern_c ) ; nim->quatern_d = FIXED_FLOAT( nhdr.quatern_d ) ; nim->qoffset_x = FIXED_FLOAT(nhdr.qoffset_x) ; nim->qoffset_y = FIXED_FLOAT(nhdr.qoffset_y) ; nim->qoffset_z = FIXED_FLOAT(nhdr.qoffset_z) ; nim->qfac = (nhdr.pixdim[0] < 0.0) ? -1.0 : 1.0 ; /* left-handedness? */ nim->qto_xyz = quatern_to_mat44( nim->quatern_b, nim->quatern_c, nim->quatern_d, nim->qoffset_x, nim->qoffset_y, nim->qoffset_z, nim->dx , nim->dy , nim->dz , nim->qfac ) ; nim->qform_code = nhdr.qform_code ; if( gni_debug > 1 ) nifti_disp_matrix_orient("-d qform orientations:\n", nim->qto_xyz); } /** load inverse transformation (x,y,z) -> (i,j,k) **/ nim->qto_ijk = mat44_inverse( nim->qto_xyz ) ; /** load sto_xyz affine transformation, if present **/ if( !is_nifti || nhdr.sform_code <= 0 ){ /** no sto transformation **/ nim->sform_code = NIFTI_XFORM_UNKNOWN ; } else { /** sto transformation from srow_*[] **/ nim->sto_xyz.m[0][0] = nhdr.srow_x[0] ; nim->sto_xyz.m[0][1] = nhdr.srow_x[1] ; nim->sto_xyz.m[0][2] = nhdr.srow_x[2] ; nim->sto_xyz.m[0][3] = nhdr.srow_x[3] ; nim->sto_xyz.m[1][0] = nhdr.srow_y[0] ; nim->sto_xyz.m[1][1] = nhdr.srow_y[1] ; nim->sto_xyz.m[1][2] = nhdr.srow_y[2] ; nim->sto_xyz.m[1][3] = nhdr.srow_y[3] ; nim->sto_xyz.m[2][0] = nhdr.srow_z[0] ; nim->sto_xyz.m[2][1] = nhdr.srow_z[1] ; nim->sto_xyz.m[2][2] = nhdr.srow_z[2] ; nim->sto_xyz.m[2][3] = nhdr.srow_z[3] ; /* last row is always [ 0 0 0 1 ] */ nim->sto_xyz.m[3][0]=nim->sto_xyz.m[3][1]=nim->sto_xyz.m[3][2] = 0.0; nim->sto_xyz.m[3][3]= 1.0 ; nim->sto_ijk = mat44_inverse( nim->sto_xyz ) ; nim->sform_code = nhdr.sform_code ; if( gni_debug > 1 ) nifti_disp_matrix_orient("-d qform orientations:\n", nim->sto_xyz); } /* miscellaneous NIFTI stuff */ if( is_nifti ){ nim->scl_slope = FIXED_FLOAT( nhdr.scl_slope ) ; nim->scl_inter = FIXED_FLOAT( nhdr.scl_inter ) ; nim->intent_code = nhdr.intent_code ; nim->intent_p1 = FIXED_FLOAT( nhdr.intent_p1 ) ; nim->intent_p2 = FIXED_FLOAT( nhdr.intent_p2 ) ; nim->intent_p3 = FIXED_FLOAT( nhdr.intent_p3 ) ; nim->toffset = FIXED_FLOAT( nhdr.toffset ) ; memcpy(nim->intent_name,nhdr.intent_name,15); nim->intent_name[15] = '\0'; nim->xyz_units = XYZT_TO_SPACE(nhdr.xyzt_units) ; nim->time_units = XYZT_TO_TIME (nhdr.xyzt_units) ; nim->freq_dim = DIM_INFO_TO_FREQ_DIM ( nhdr.dim_info ) ; nim->phase_dim = DIM_INFO_TO_PHASE_DIM( nhdr.dim_info ) ; nim->slice_dim = DIM_INFO_TO_SLICE_DIM( nhdr.dim_info ) ; nim->slice_code = nhdr.slice_code ; nim->slice_start = nhdr.slice_start ; nim->slice_end = nhdr.slice_end ; nim->slice_duration = FIXED_FLOAT(nhdr.slice_duration) ; } /* Miscellaneous ANALYZE stuff */ nim->cal_min = FIXED_FLOAT(nhdr.cal_min) ; nim->cal_max = FIXED_FLOAT(nhdr.cal_max) ; memcpy(nim->descrip ,nhdr.descrip ,79) ; nim->descrip [79] = '\0' ; memcpy(nim->aux_file,nhdr.aux_file,23) ; nim->aux_file[23] = '\0' ; /** will read image data from file 'iname' starting at offset 'ioff' **/ is_onefile = is_nifti && NIFTI_ONEFILE(nhdr) ; if( is_onefile ){ ioff = (int)nhdr.vox_offset ; if( ioff < sizeof(nhdr) ) ioff = sizeof(nhdr) ; } else { ioff = (int)nhdr.vox_offset ; } nim->iname_offset = ioff ; /* deal with file names if set */ if (fname!=NULL) { nim->fname = nifti_strdup(fname); /* determine name of image, if not already set */ if (nim->iname==NULL) { if (is_onefile) { iname = nifti_strdup(nim->fname); } else { iname = nifti_findimgname(nim->fname,nim->nifti_type); if (iname==NULL) ERREX("bad filename"); } /* don't free iname, as now nim->iname is using this storage */ nim->iname = iname ; /* save image filename */ } } else { nim->fname = NULL; nim->iname = NULL; } /* clear extension fields */ nim->num_ext = 0; nim->ext_list = NULL; return nim; }

struct nifti_1_header nifti_convert_nim2nhdr (  nifti_image *    nim )

load the ANALYZE-7.5 generic parts of the header struct * Definition at line 3119 of file nifti1_io.c. References nifti_1_header::aux_file, nifti_1_header::bitpix, nifti_1_header::cal_max, nifti_1_header::cal_min, nifti_1_header::datatype, nifti_1_header::descrip, nifti_1_header::dim, nifti_1_header::dim_info, FPS_INTO_DIM_INFO, nifti_1_header::intent_code, nifti_1_header::intent_name, nifti_1_header::intent_p1, nifti_1_header::intent_p2, nifti_1_header::intent_p3, nifti_1_header::magic, nifti_1_header::pixdim, nifti_1_header::qform_code, nifti_1_header::qoffset_x, nifti_1_header::qoffset_y, nifti_1_header::qoffset_z, nifti_1_header::quatern_b, nifti_1_header::quatern_c, nifti_1_header::quatern_d, nifti_1_header::regular, nifti_1_header::scl_inter, nifti_1_header::scl_slope, nifti_1_header::sform_code, nifti_1_header::sizeof_hdr, nifti_1_header::slice_code, nifti_1_header::slice_duration, nifti_1_header::slice_end, nifti_1_header::slice_start, SPACE_TIME_TO_XYZT, nifti_1_header::srow_x, nifti_1_header::srow_y, nifti_1_header::srow_z, nifti_1_header::toffset, nifti_1_header::vox_offset, and nifti_1_header::xyzt_units. Referenced by nifti_copy_nim_info(), and nifti_image_write_hdr_img2(). { struct nifti_1_header nhdr; memset(&nhdr,0,sizeof(nhdr)) ; /* zero out header, to be safe */ /** load the ANALYZE-7.5 generic parts of the header struct **/ nhdr.sizeof_hdr = sizeof(nhdr) ; nhdr.regular = 'r' ; /* for some stupid reason */ nhdr.dim[0] = nim->ndim ; nhdr.dim[1] = nim->nx ; nhdr.dim[2] = nim->ny ; nhdr.dim[3] = nim->nz ; nhdr.dim[4] = nim->nt ; nhdr.dim[5] = nim->nu ; nhdr.dim[6] = nim->nv ; nhdr.dim[7] = nim->nw ; nhdr.pixdim[0] = 0.0 ; nhdr.pixdim[1] = nim->dx ; nhdr.pixdim[2] = nim->dy ; nhdr.pixdim[3] = nim->dz ; nhdr.pixdim[4] = nim->dt ; nhdr.pixdim[5] = nim->du ; nhdr.pixdim[6] = nim->dv ; nhdr.pixdim[7] = nim->dw ; nhdr.datatype = nim->datatype ; nhdr.bitpix = 8 * nim->nbyper ; if( nim->cal_max > nim->cal_min ){ nhdr.cal_max = nim->cal_max ; nhdr.cal_min = nim->cal_min ; } if( nim->scl_slope != 0.0 ){ nhdr.scl_slope = nim->scl_slope ; nhdr.scl_inter = nim->scl_inter ; } if( nim->descrip[0] != '\0' ){ memcpy(nhdr.descrip ,nim->descrip ,79) ; nhdr.descrip[79] = '\0' ; } if( nim->aux_file[0] != '\0' ){ memcpy(nhdr.aux_file ,nim->aux_file ,23) ; nhdr.aux_file[23] = '\0' ; } /** Load NIFTI specific stuff into the header **/ if( nim->nifti_type > 0 ){ if( nim->nifti_type == 1 ) strcpy(nhdr.magic,"n+1") ; /* 1 file */ else strcpy(nhdr.magic,"ni1") ; /* 2 files */ /* for nifti files, enforce positive pixdims! */ nhdr.pixdim[1] = fabs(nhdr.pixdim[1]) ; nhdr.pixdim[2] = fabs(nhdr.pixdim[2]) ; nhdr.pixdim[3] = fabs(nhdr.pixdim[3]) ; nhdr.pixdim[4] = fabs(nhdr.pixdim[4]) ; nhdr.pixdim[5] = fabs(nhdr.pixdim[5]) ; nhdr.pixdim[6] = fabs(nhdr.pixdim[6]) ; nhdr.pixdim[7] = fabs(nhdr.pixdim[7]) ; nhdr.intent_code = nim->intent_code ; nhdr.intent_p1 = nim->intent_p1 ; nhdr.intent_p2 = nim->intent_p2 ; nhdr.intent_p3 = nim->intent_p3 ; if( nim->intent_name[0] != '\0' ){ memcpy(nhdr.intent_name,nim->intent_name,15) ; nhdr.intent_name[15] = '\0' ; } nhdr.vox_offset = (float) nim->iname_offset ; nhdr.xyzt_units = SPACE_TIME_TO_XYZT( nim->xyz_units, nim->time_units ) ; nhdr.toffset = nim->toffset ; if( nim->qform_code > 0 ){ nhdr.qform_code = nim->qform_code ; nhdr.quatern_b = nim->quatern_b ; nhdr.quatern_c = nim->quatern_c ; nhdr.quatern_d = nim->quatern_d ; nhdr.qoffset_x = nim->qoffset_x ; nhdr.qoffset_y = nim->qoffset_y ; nhdr.qoffset_z = nim->qoffset_z ; nhdr.pixdim[0] = (nim->qfac >= 0.0) ? 1.0 : -1.0 ; } if( nim->sform_code > 0 ){ nhdr.sform_code = nim->sform_code ; nhdr.srow_x[0] = nim->sto_xyz.m[0][0] ; nhdr.srow_x[1] = nim->sto_xyz.m[0][1] ; nhdr.srow_x[2] = nim->sto_xyz.m[0][2] ; nhdr.srow_x[3] = nim->sto_xyz.m[0][3] ; nhdr.srow_y[0] = nim->sto_xyz.m[1][0] ; nhdr.srow_y[1] = nim->sto_xyz.m[1][1] ; nhdr.srow_y[2] = nim->sto_xyz.m[1][2] ; nhdr.srow_y[3] = nim->sto_xyz.m[1][3] ; nhdr.srow_z[0] = nim->sto_xyz.m[2][0] ; nhdr.srow_z[1] = nim->sto_xyz.m[2][1] ; nhdr.srow_z[2] = nim->sto_xyz.m[2][2] ; nhdr.srow_z[3] = nim->sto_xyz.m[2][3] ; } nhdr.dim_info = FPS_INTO_DIM_INFO( nim->freq_dim , nim->phase_dim , nim->slice_dim ) ; nhdr.slice_code = nim->slice_code ; nhdr.slice_start = nim->slice_start ; nhdr.slice_end = nim->slice_end ; nhdr.slice_duration = nim->slice_duration ; } return nhdr; }

nifti_image* nifti_copy_nim_info (  nifti_image *    src )

Definition at line 3446 of file nifti1_io.c. References calloc, nifti_image::fname, nifti_convert_nhdr2nim(), and nifti_convert_nim2nhdr(). { struct nifti_1_header nhdr; nifti_image *dest; dest = (nifti_image *)calloc(1,sizeof(nifti_image)); nhdr = nifti_convert_nim2nhdr(src); dest = nifti_convert_nhdr2nim(nhdr,src->fname); return dest; }

int nifti_copynsort (  int    nbricks, int *    blist, int **    slist, int **    sindex )  [static]

Referenced by nifti_image_load_bricks().

void nifti_datatype_sizes (  int    datatype, int *    nbyper, int *    swapsize )

Definition at line 807 of file nifti1_io.c. References ASSIF, DT_COMPLEX128, DT_COMPLEX256, DT_COMPLEX64, DT_FLOAT128, DT_FLOAT32, DT_FLOAT64, DT_INT16, DT_INT32, DT_INT64, DT_INT8, DT_RGB24, DT_UINT16, DT_UINT32, DT_UINT64, and DT_UINT8. Referenced by nifti_convert_nhdr2nim(), nifti_image_from_ascii(), and nifti_simple_init_nim(). { int nb=0, ss=0 ; switch( datatype ){ case DT_INT8: case DT_UINT8: nb = 1 ; ss = 0 ; break ; case DT_INT16: case DT_UINT16: nb = 2 ; ss = 2 ; break ; case DT_RGB24: nb = 3 ; ss = 0 ; break ; case DT_INT32: case DT_UINT32: case DT_FLOAT32: nb = 4 ; ss = 4 ; break ; case DT_COMPLEX64: nb = 8 ; ss = 4 ; break ; case DT_FLOAT64: case DT_INT64: case DT_UINT64: nb = 8 ; ss = 8 ; break ; case DT_FLOAT128: nb = 16 ; ss = 16 ; break ; case DT_COMPLEX128: nb = 16 ; ss = 8 ; break ; case DT_COMPLEX256: nb = 32 ; ss = 16 ; break ; } ASSIF(nbyper,nb) ; ASSIF(swapsize,ss) ; return ; }

char* nifti_datatype_string (  int    dt )

Definition at line 630 of file nifti1_io.c. References dt, DT_BINARY, DT_COMPLEX128, DT_COMPLEX256, DT_COMPLEX64, DT_FLOAT128, DT_FLOAT32, DT_FLOAT64, DT_INT16, DT_INT32, DT_INT64, DT_INT8, DT_RGB24, DT_UINT16, DT_UINT32, DT_UINT64, DT_UINT8, and DT_UNKNOWN. Referenced by nifti_image_to_ascii(). { switch( dt ){ case DT_UNKNOWN: return "UNKNOWN" ; case DT_BINARY: return "BINARY" ; case DT_INT8: return "INT8" ; case DT_UINT8: return "UINT8" ; case DT_INT16: return "INT16" ; case DT_UINT16: return "UINT16" ; case DT_INT32: return "INT32" ; case DT_UINT32: return "UINT32" ; case DT_INT64: return "INT64" ; case DT_UINT64: return "UINT64" ; case DT_FLOAT32: return "FLOAT32" ; case DT_FLOAT64: return "FLOAT64" ; case DT_FLOAT128: return "FLOAT128" ; case DT_COMPLEX64: return "COMPLEX64" ; case DT_COMPLEX128: return "COMPLEX128" ; case DT_COMPLEX256: return "COMPLEX256" ; case DT_RGB24: return "RGB24" ; } return "**ILLEGAL**" ; }

void nifti_disp_lib_hist (  void    )

Definition at line 120 of file nifti1_io.c. References gni_history. Referenced by main(). { fputs(gni_history, stdout); }

void nifti_disp_lib_version (  void    )

Definition at line 125 of file nifti1_io.c. References gni_version. Referenced by main(). { printf("%s, compiled %s\n", gni_version, __DATE__); }

int nifti_disp_matrix_orient (  char *    mesg, mat44    mat )

Definition at line 600 of file nifti1_io.c. References i, mat, mat44_to_orientation(), and nifti_orientation_string(). Referenced by nifti_convert_nhdr2nim(). { int i, j, k; if ( mesg ) fputs( mesg, stderr ); /* use stdout? */ mat44_to_orientation( mat, &i,&j,&k ); if ( i <= 0 || j <= 0 || k <= 0 ) return -1; /* so we have good codes */ fprintf(stderr, " i orientation = '%s'\n" " j orientation = '%s'\n" " k orientation = '%s'\n", nifti_orientation_string(i), nifti_orientation_string(j), nifti_orientation_string(k) ); return 0; }

char* nifti_find_file_extension (  char *    name )

check input file(s) for sanity * Definition at line 1699 of file nifti1_io.c. References name. Referenced by nifti_findhdrname(), nifti_makebasename(), and nifti_validfilename(). { char * ext; int len; if ( ! name ) return NULL; len = strlen(name); if ( len < 4 ) return NULL; ext = name + len - 4; if ( (strcmp(ext, ".hdr") == 0) || (strcmp(ext, ".img") == 0) || (strcmp(ext, ".nia") == 0) || (strcmp(ext, ".nii") == 0) ) return ext; #ifdef HAVE_ZLIB if ( len < 7 ) return NULL; ext = name - 7; if ( (strcmp(ext, ".hdr.gz") == 0) || (strcmp(ext, ".img.gz") == 0) || (strcmp(ext, ".nia.gz") == 0) || (strcmp(ext, ".nii.gz") == 0) ) return ext; #endif return NULL; }

char* nifti_findhdrname (  char *    fname )  [static]

Referenced by is_nifti_file(), and nifti_image_read().

char* nifti_findimgname (  char *    fname, int    nifti_type )  [static]

Referenced by nifti_convert_nhdr2nim(), and nifti_image_write_hdr_img2().

size_t nifti_get_volsize (  nifti_image *    nim )

Definition at line 1652 of file nifti1_io.c. References nifti_image::nbyper, and nifti_image::nvox. Referenced by nifti_image_load(), and nifti_image_load_prep(). { return (size_t)(nim->nbyper) * (size_t)(nim->nvox) ; /* total bytes */ }

void nifti_image_free (  nifti_image *    nim )

Definition at line 2966 of file nifti1_io.c. References nifti_image::data, nifti_image::fname, free, and nifti_image::iname. Referenced by nifti_image_from_ascii(), nifti_image_read_bricks(), THD_load_nifti(), THD_open_analyze(), and THD_open_nifti(). { if( nim == NULL ) return ; if( nim->fname != NULL ) free(nim->fname) ; if( nim->iname != NULL ) free(nim->iname) ; if( nim->data != NULL ) free(nim->data ) ; free(nim) ; return ; }

nifti_image* nifti_image_from_ascii (  char *    str, int *    bytes_read )

Definition at line 3899 of file nifti1_io.c. References nifti_image::byteorder, calloc, nifti_image::datatype, nifti_image::dim, dt, nifti_image::dt, nifti_image::du, nifti_image::dv, nifti_image::dw, nifti_image::dx, nifti_image::dy, nifti_image::dz, nifti_image::fname, nifti_image::iname, mat44::m, mat44_inverse(), nifti_image::nbyper, nifti_image::ndim, nifti_datatype_sizes(), nifti_image_free(), nifti_strdup(), nifti_image::nifti_type, nifti_image::nt, nifti_image::nu, nifti_image::nv, nifti_image::nvox, nifti_image::nw, nifti_image::nx, nifti_image::ny, nz, nifti_image::nz, nifti_image::pixdim, nifti_image::qfac, nifti_image::qform_code, QNUM, nifti_image::qoffset_x, nifti_image::qoffset_y, nifti_image::qoffset_z, QQNUM, QSTR, nifti_image::qto_ijk, nifti_image::qto_xyz, nifti_image::quatern_b, nifti_image::quatern_c, nifti_image::quatern_d, quatern_to_mat44(), nifti_image::sform_code, short_order(), nifti_image::sto_ijk, nifti_image::sto_xyz, nifti_image::swapsize, and unescape_string(). Referenced by nifti_image_read(). { char lhs[1024] , rhs[1024] ; int ii , spos, nn , slen ; nifti_image *nim ; /* will be output */ if( str == NULL || *str == '\0' ) return NULL ; /* bad input!? */ /* scan for opening string */ spos = 0 ; slen = strlen(str) ; ii = sscanf( str+spos , "%1023s%n" , lhs , &nn ) ; spos += nn ; if( ii == 0 || strcmp(lhs,"<nifti_image") != 0 ) return NULL ; /* create empty image struct */ nim = (nifti_image *) calloc( 1 , sizeof(nifti_image) ) ; nim->nx = nim->ny = nim->nz = nim->nt = nim->nu = nim->nv = nim->nw = 1 ; nim->dx = nim->dy = nim->dz = nim->dt = nim->du = nim->dv = nim->dw = nim->qfac = 1.0 ; nim->byteorder = short_order() ; /* starting at str[spos], scan for "equations" of the form lhs = 'rhs' and assign rhs values into the struct component named by lhs */ while(1){ while( isspace(str[spos]) ) spos++ ; /* skip whitespace */ if( str[spos] == '\0' ) break ; /* end of string? */ /* get lhs string */ ii = sscanf( str+spos , "%1023s%n" , lhs , &nn ) ; spos += nn ; if( ii == 0 || strcmp(lhs,"/>") == 0 ) break ; /* end of input? */ /* skip whitespace and the '=' marker */ while( isspace(str[spos]) || str[spos] == '=' ) spos++ ; if( str[spos] == '\0' ) break ; /* end of string? */ /* if next character is a quote ', copy everything up to next ' otherwise, copy everything up to next nonblank */ if( str[spos] == '\'' ){ ii = spos+1 ; while( str[ii] != '\0' && str[ii] != '\'' ) ii++ ; nn = ii-spos-1 ; if( nn > 1023 ) nn = 1023 ; memcpy(rhs,str+spos+1,nn) ; rhs[nn] = '\0' ; spos = (str[ii] == '\'') ? ii+1 : ii ; } else { ii = sscanf( str+spos , "%1023s%n" , rhs , &nn ) ; spos += nn ; if( ii == 0 ) break ; /* nothing found? */ } unescape_string(rhs) ; /* remove any XML escape sequences */ /* Now can do the assignment, based on lhs string. Start with special cases that don't fit the QNUM/QSTR macros. */ if( strcmp(lhs,"nifti_type") == 0 ){ if( strcmp(rhs,"ANALYZE-7.5") == 0 ) nim->nifti_type = 0 ; else if( strcmp(rhs,"NIFTI-1+") == 0 ) nim->nifti_type = 1 ; else if( strcmp(rhs,"NIFTI-1") == 0 ) nim->nifti_type = 2 ; else if( strcmp(rhs,"NIFTI-1A") == 0 ) nim->nifti_type = 3 ; } else if( strcmp(lhs,"header_filename") == 0 ){ nim->fname = nifti_strdup(rhs) ; } else if( strcmp(lhs,"image_filename") == 0 ){ nim->iname = nifti_strdup(rhs) ; } else if( strcmp(lhs,"sto_xyz_matrix") == 0 ){ sscanf( rhs , "%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f" , &(nim->sto_xyz.m[0][0]) , &(nim->sto_xyz.m[0][1]) , &(nim->sto_xyz.m[0][2]) , &(nim->sto_xyz.m[0][3]) , &(nim->sto_xyz.m[1][0]) , &(nim->sto_xyz.m[1][1]) , &(nim->sto_xyz.m[1][2]) , &(nim->sto_xyz.m[1][3]) , &(nim->sto_xyz.m[2][0]) , &(nim->sto_xyz.m[2][1]) , &(nim->sto_xyz.m[2][2]) , &(nim->sto_xyz.m[2][3]) , &(nim->sto_xyz.m[3][0]) , &(nim->sto_xyz.m[3][1]) , &(nim->sto_xyz.m[3][2]) , &(nim->sto_xyz.m[3][3]) ) ; } else if( strcmp(lhs,"byteorder") == 0 ){ if( strcmp(rhs,"MSB_FIRST") == 0 ) nim->byteorder = MSB_FIRST ; if( strcmp(rhs,"LSB_FIRST") == 0 ) nim->byteorder = LSB_FIRST ; } else QQNUM(image_offset,iname_offset) ; else QNUM(datatype) ; else QNUM(ndim) ; else QNUM(nx) ; else QNUM(ny) ; else QNUM(nz) ; else QNUM(nt) ; else QNUM(nu) ; else QNUM(nv) ; else QNUM(nw) ; else QNUM(dx) ; else QNUM(dy) ; else QNUM(dz) ; else QNUM(dt) ; else QNUM(du) ; else QNUM(dv) ; else QNUM(dw) ; else QNUM(cal_min) ; else QNUM(cal_max) ; else QNUM(scl_slope) ; else QNUM(scl_inter) ; else QNUM(intent_code) ; else QNUM(intent_p1) ; else QNUM(intent_p2) ; else QNUM(intent_p3) ; else QSTR(intent_name,15) ; else QNUM(toffset) ; else QNUM(xyz_units) ; else QNUM(time_units) ; else QSTR(descrip,79) ; else QSTR(aux_file,23) ; else QNUM(qform_code) ; else QNUM(quatern_b) ; else QNUM(quatern_c) ; else QNUM(quatern_d) ; else QNUM(qoffset_x) ; else QNUM(qoffset_y) ; else QNUM(qoffset_z) ; else QNUM(qfac) ; else QNUM(sform_code) ; else QNUM(freq_dim) ; else QNUM(phase_dim) ; else QNUM(slice_dim) ; else QNUM(slice_code) ; else QNUM(slice_start) ; else QNUM(slice_end) ; else QNUM(slice_duration) ; else QNUM(num_ext) ; } /* end of while loop */ if( bytes_read ) *bytes_read = spos+1; /* "process" last '\n' */ /** do miscellaneous checking and cleanup **/ if( nim->ndim <= 0 ){ nifti_image_free(nim); return NULL; } /** bad! **/ nifti_datatype_sizes( nim->datatype, &(nim->nbyper), &(nim->swapsize) ); if( nim->nbyper == 0 ){ nifti_image_free(nim); return NULL; } /** bad! **/ nim->dim[0] = nim->ndim ; nim->dim[1] = nim->nx ; nim->pixdim[1] = nim->dx ; nim->dim[2] = nim->ny ; nim->pixdim[2] = nim->dy ; nim->dim[3] = nim->nz ; nim->pixdim[3] = nim->dz ; nim->dim[4] = nim->nt ; nim->pixdim[4] = nim->dt ; nim->dim[5] = nim->nu ; nim->pixdim[5] = nim->du ; nim->dim[6] = nim->nv ; nim->pixdim[6] = nim->dv ; nim->dim[7] = nim->nw ; nim->pixdim[7] = nim->dw ; nim->nvox = nim->nx * nim->ny * nim->nz * nim->nt * nim->nu * nim->nv * nim->nw ; if( nim->qform_code > 0 ) nim->qto_xyz = quatern_to_mat44( nim->quatern_b, nim->quatern_c, nim->quatern_d, nim->qoffset_x, nim->qoffset_y, nim->qoffset_z, nim->dx , nim->dy , nim->dz , nim->qfac ) ; else nim->qto_xyz = quatern_to_mat44( 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , nim->dx , nim->dy , nim->dz , 0.0 ) ; nim->qto_ijk = mat44_inverse( nim->qto_xyz ) ; if( nim->sform_code > 0 ) nim->sto_ijk = mat44_inverse( nim->sto_xyz ) ; return nim ; }

void nifti_image_infodump (  nifti_image *    nim )

Definition at line 2979 of file nifti1_io.c. References free, and nifti_image_to_ascii(). Referenced by main(), and nifti_image_read(). { char *str = nifti_image_to_ascii( nim ) ; /* stdout -> stderr 2 Dec 2004 [rickr] */ if( str != NULL ){ fputs(str,stderr) ; free(str) ; } return ; }

int nifti_image_load (  nifti_image *    nim )

open image data file * Definition at line 2821 of file nifti1_io.c. References calloc, nifti_image::data, free, gni_debug, nifti_get_volsize(), nifti_image_load_prep(), and nifti_read_buffer(). Referenced by nifti_image_read(). { /* set up data space, open data file and seek, then call nifti_read_buffer */ size_t ntot , ii ; znzFile fp ; fp = nifti_image_load_prep( nim ); if( fp == NULL ){ if( gni_debug > 0 ) fprintf(stderr,"** nifti_image_load, failed load_prep\n"); return -1; } ntot = nifti_get_volsize(nim); /* if the data pointer is already set, use the given space */ if( nim->data == NULL ) { nim->data = (void *)calloc(1,ntot) ; /* create image memory */ if( nim->data == NULL ){ if( gni_debug > 0 ) fprintf(stderr,"** failed to alloc %d bytes for image data\n", (int)ntot); znzclose(fp); return -1; } } /* now that everything is set up - do the reading */ ii = nifti_read_buffer(fp,nim->data,ntot,nim); if( ii < ntot ){ znzclose(fp) ; free(nim->data) ; nim->data = NULL ; return -1 ; /* errors were printed in nifti_read_buffer() */ } znzclose( fp ) ; return 0 ; }

int nifti_image_load_bricks (  nifti_image *    nim, int    nbricks, int *    blist, nifti_brick_list *    NBL )

Definition at line 260 of file nifti1_io.c. References free, free_NBL(), gni_debug, nifti_brick_list::nbricks, nifti_alloc_NBL_mem(), nifti_copynsort(), nifti_image_load_prep(), nifti_load_NBL_bricks(), and valid_nifti_brick_list(). Referenced by nifti_image_read_bricks(). { int * slist = NULL, * sindex = NULL, rv; znzFile fp; /* we can have blist == NULL */ if( !nim || !NBL ){ fprintf(stderr,"** nifti_image_load_bricks, bad params (%p,%p)\n", nim, NBL); return -1; } if( blist && ! valid_nifti_brick_list( nim, nbricks, blist, 1 ) ) return -1; /* for efficiency, let's read the file in order */ if( blist && nifti_copynsort( nbricks, blist, &slist, &sindex ) != 0 ) return -1; /* open the file and position the FILE pointer */ fp = nifti_image_load_prep( nim ); if( !fp ){ if( gni_debug > 0 ) fprintf(stderr,"** nifti_image_load_bricks, failed load_prep\n"); if( blist ){ free(slist); free(sindex); } return -1; } /* this will flag to allocate defaults */ if( !blist ) nbricks = 0; if( nifti_alloc_NBL_mem( nim, nbricks, NBL ) != 0 ){ if( blist ){ free(slist); free(sindex); } znzclose(fp); return -1; } rv = nifti_load_NBL_bricks(nim, slist, sindex, NBL, fp); if( rv != 0 ){ free_NBL( NBL ); /* failure! */ NBL->nbricks = 0; /* repetative, but clear */ } if( slist ){ free(slist); free(sindex); } znzclose(fp); return NBL->nbricks; }

znzFile nifti_image_load_prep (  nifti_image *    nim )  [static]

Referenced by nifti_image_load(), and nifti_image_load_bricks().

znzFile nifti_image_open (  char *    hname, char *    opts, nifti_image **    nim )  [static]

nifti_image* nifti_image_read (  char *    hname, int    read_data )

open image data file * Definition at line 2349 of file nifti1_io.c. References calloc, free, get_filesize(), gni_debug, ilen, LNI_FERR, nifti_convert_nhdr2nim(), nifti_findhdrname(), nifti_image_from_ascii(), nifti_image_infodump(), nifti_image_load(), nifti_is_gzfile(), and nifti_read_extensions(). Referenced by main(), nifti_image_open(), nifti_image_read_bricks(), THD_load_nifti(), THD_open_analyze(), and THD_open_nifti(). { struct nifti_1_header nhdr ; nifti_image *nim ; znzFile fp ; int ii , ilen; char fname[] = { "nifti_image_read" }; char buf[16] , *workingname=NULL; /* determine file name to use for header */ workingname = nifti_findhdrname(hname); if( workingname == NULL ){ if( gni_debug > 0 ) LNI_FERR(fname,"failed to find header file for", hname); return NULL; /* check return */ } else if( gni_debug > 1 ) fprintf(stderr,"-d %s: found header filename '%s'\n",fname,workingname); fp = znzopen(workingname,"rb",1); if( znz_isnull(fp) ){ if( gni_debug > 0 ) LNI_FERR(fname,"failed to open header file",workingname); free(workingname); return NULL; } /** test if header file starts with ASCII string "<nifti_image"; if so, read the dataset that special way and return now (this is NOT part of the NIFTI-1 standard!) **/ ii = znzread( buf , 1 , 12 , fp ) ; if( ii < 12 ){ znzclose( fp ) ; if( gni_debug > 0 ) LNI_FERR(fname,"short header read",workingname); free(workingname); return NULL; } znzrewind(fp) ; buf[12] = '\0' ; if( strcmp(buf,"<nifti_image") == 0 ){ /* have an ASCII header! */ int slen, txt_size; char *sbuf ; if( nifti_is_gzfile(workingname) ){ LNI_FERR(fname, "compressed file with negative offset", workingname); free(workingname); return NULL; } slen = get_filesize( workingname ) ; if( gni_debug > 1 ) fprintf(stderr,"-d %s: have ASCII image of size %d\n",fname,slen); if( slen > 65530 ) slen = 65530 ; sbuf = (char *)calloc(sizeof(char),slen+1) ; znzread( sbuf , 1 , slen , fp ) ; nim = nifti_image_from_ascii( sbuf, &txt_size ) ; free( sbuf ) ; if( nim == NULL ){ LNI_FERR(fname,"failed nifti_image_from_ascii()",workingname); free(workingname); return NULL; } nim->nifti_type = 3 ; /* read extensions (but reposition file pointer, first */ znzseek(fp, txt_size, SEEK_SET); (void) nifti_read_extensions(nim, fp); znzclose( fp ) ; nim->iname_offset = -1 ; /* now check from the end of the file */ if( read_data ) ii = nifti_image_load( nim ) ; else nim->data = NULL ; free(workingname); /* check for nifti_image_load() failure, maybe bail out */ if( read_data && ii != 0 ){ if( gni_debug > 1 ) fprintf(stderr,"-d failed image_load, free nifti image struct\n"); free(nim); return NULL; } return nim ; } /**================ Normal case: read binary header ===================*/ ii = znzread( &nhdr , 1 , sizeof(nhdr) , fp ) ; /* read the thing */ /* rcr - test header struct for consistency here? */ /* keep the file open for a while so we can check for extensions * after nifti_convert_nhdr2nim() */ if( ii < sizeof(nhdr) ){ if( gni_debug > 0 ){ LNI_FERR(fname,"bad binary header read for file", workingname); fprintf(stderr," - read %d of %d bytes\n",ii, (int)sizeof(nhdr)); } free(workingname); return NULL; } /*=== create output image struct and set it up ===*/ /** convert all nhdr fields to nifti_image fields **/ nim = nifti_convert_nhdr2nim(nhdr,workingname); if( nim == NULL ){ znzclose( fp ) ; /* close the file */ if( gni_debug > 0 ) LNI_FERR(fname,"failed to create nifti image from header",workingname); free(workingname); /* had to save this for debug message */ return NULL; } if( gni_debug > 1 ){ fprintf(stderr,"+d nifti_image_read(), have nifti image:\n"); nifti_image_infodump(nim); } /* check for extensions, let us just ignore any errors here */ (void) nifti_read_extensions(nim, fp); znzclose( fp ) ; /* close the file */ free(workingname); /* check that the data file is OK */ if (!nifti_is_gzfile(nim->iname)) { /* can't use get_filesize on a compressed file :( */ ilen = get_filesize(nim->iname) ; if( ilen <= nim->iname_offset ){ fprintf(stderr,"ERROR %s: file size less than offset (%d < %d)\n", fname,ilen,nim->iname_offset); free(nim); return NULL; } } /** read the data if desired, then bug out **/ if( read_data ) ii = nifti_image_load( nim ) ; else nim->data = NULL ; return nim ; }

nifti_image* nifti_image_read_bricks (  char *    hname, int    nbricks, int *    blist, nifti_brick_list *    NBL )

Definition at line 178 of file nifti1_io.c. References nifti_image_free(), nifti_image_load_bricks(), nifti_image_read(), and update_nifti_image_for_brick_list(). { nifti_image * nim; if( !hname || !NBL ){ fprintf(stderr,"** nifti_image_read_bricks: bad params (%p,%p)\n", hname, NBL); return NULL; } if( blist && nbricks <= 0 ){ fprintf(stderr,"** nifti_image_read_bricks: bad nbricks, %d\n", nbricks); return NULL; } nim = nifti_image_read(hname, 0); /* read header, but not data */ if( !nim ) return NULL; /* errors were already printed */ /* if we fail, free image and return */ if( nifti_image_load_bricks(nim, nbricks, blist, NBL) <= 0 ){ nifti_image_free(nim); return NULL; } if( blist ) update_nifti_image_for_brick_list(nim, nbricks); return nim; }

char* nifti_image_to_ascii (  nifti_image *    nim )

Definition at line 3631 of file nifti1_io.c. References nifti_image::aux_file, nifti_image::byteorder, nifti_image::cal_max, nifti_image::cal_min, calloc, nifti_image::datatype, nifti_image::descrip, nifti_image::dt, nifti_image::du, nifti_image::dv, nifti_image::dw, nifti_image::dx, nifti_image::dy, nifti_image::dz, escapize_string(), nifti_image::fname, free, nifti_image::freq_dim, i, nifti_image::iname, nifti_image::iname_offset, nifti_image::intent_code, nifti_image::intent_name, nifti_image::intent_p1, nifti_image::intent_p2, nifti_image::intent_p3, mat44::m, mat44_to_orientation(), nifti_image::nbyper, nifti_image::ndim, nifti_datatype_string(), nifti_intent_string(), nifti_orientation_string(), nifti_slice_string(), nifti_image::nifti_type, nifti_units_string(), nifti_xform_string(), nifti_image::nt, nifti_image::nu, nifti_image::num_ext, nifti_image::nv, nifti_image::nvox, nifti_image::nw, nifti_image::nx, nifti_image::ny, nifti_image::nz, nifti_image::phase_dim, nifti_image::qfac, nifti_image::qform_code, nifti_image::qoffset_x, nifti_image::qoffset_y, nifti_image::qoffset_z, nifti_image::qto_ijk, nifti_image::qto_xyz, nifti_image::quatern_b, nifti_image::quatern_c, nifti_image::quatern_d, realloc, nifti_image::scl_inter, nifti_image::scl_slope, nifti_image::sform_code, nifti_image::slice_code, nifti_image::slice_dim, nifti_image::slice_duration, nifti_image::slice_end, nifti_image::slice_start, nifti_image::sto_ijk, nifti_image::sto_xyz, nifti_image::time_units, nifti_image::toffset, and nifti_image::xyz_units. Referenced by nifti_image_infodump(), and nifti_image_write_hdr_img2(). { char *buf , *ebuf ; int nbuf ; if( nim == NULL ) return NULL ; /* stupid caller */ buf = (char *)calloc(1,65534); nbuf = 0; /* longer than needed, to be safe */ sprintf( buf , "<nifti_image\n" ) ; /* XML-ish opener */ sprintf( buf+strlen(buf) , " nifti_type = '%s'\n" , (nim->nifti_type == 1) ? "NIFTI-1+" :(nim->nifti_type == 2) ? "NIFTI-1" :(nim->nifti_type == 3) ? "NIFTI-1A" : "ANALYZE-7.5" ) ; /** Strings that we don't control (filenames, etc.) that might contain "weird" characters (like quotes) are "escaped": - A few special characters are replaced by XML-style escapes, using the function escapize_string(). - On input, function unescape_string() reverses this process. - The result is that the NIFTI ASCII-format header is XML-compliant. **/ ebuf = escapize_string(nim->fname) ; sprintf( buf+strlen(buf) , " header_filename = %s\n",ebuf); free(ebuf); ebuf = escapize_string(nim->iname) ; sprintf( buf+strlen(buf) , " image_filename = %s\n", ebuf); free(ebuf); sprintf( buf+strlen(buf) , " image_offset = '%d'\n" , nim->iname_offset ); sprintf( buf+strlen(buf), " ndim = '%d'\n", nim->ndim); sprintf( buf+strlen(buf), " nx = '%d'\n", nim->nx ); if( nim->ndim > 1 ) sprintf( buf+strlen(buf), " ny = '%d'\n", nim->ny ); if( nim->ndim > 2 ) sprintf( buf+strlen(buf), " nz = '%d'\n", nim->nz ); if( nim->ndim > 3 ) sprintf( buf+strlen(buf), " nt = '%d'\n", nim->nt ); if( nim->ndim > 4 ) sprintf( buf+strlen(buf), " nu = '%d'\n", nim->nu ); if( nim->ndim > 5 ) sprintf( buf+strlen(buf), " nv = '%d'\n", nim->nv ); if( nim->ndim > 6 ) sprintf( buf+strlen(buf), " nw = '%d'\n", nim->nw ); sprintf( buf+strlen(buf), " dx = '%g'\n", nim->dx ); if( nim->ndim > 1 ) sprintf( buf+strlen(buf), " dy = '%g'\n", nim->dy ); if( nim->ndim > 2 ) sprintf( buf+strlen(buf), " dz = '%g'\n", nim->dz ); if( nim->ndim > 3 ) sprintf( buf+strlen(buf), " dt = '%g'\n", nim->dt ); if( nim->ndim > 4 ) sprintf( buf+strlen(buf), " du = '%g'\n", nim->du ); if( nim->ndim > 5 ) sprintf( buf+strlen(buf), " dv = '%g'\n", nim->dv ); if( nim->ndim > 6 ) sprintf( buf+strlen(buf), " dw = '%g'\n", nim->dw ); sprintf( buf+strlen(buf) , " datatype = '%d'\n" , nim->datatype ) ; sprintf( buf+strlen(buf) , " datatype_name = '%s'\n" , nifti_datatype_string(nim->datatype) ) ; sprintf( buf+strlen(buf) , " nvox = '%d'\n" , nim->nvox ) ; sprintf( buf+strlen(buf) , " nbyper = '%d'\n" , nim->nbyper ) ; sprintf( buf+strlen(buf) , " byteorder = '%s'\n" , (nim->byteorder==MSB_FIRST) ? "MSB_FIRST" : "LSB_FIRST" ) ; if( nim->cal_min < nim->cal_max ){ sprintf( buf+strlen(buf) , " cal_min = '%g'\n", nim->cal_min ) ; sprintf( buf+strlen(buf) , " cal_max = '%g'\n", nim->cal_max ) ; } if( nim->scl_slope != 0.0 ){ sprintf( buf+strlen(buf) , " scl_slope = '%g'\n" , nim->scl_slope ) ; sprintf( buf+strlen(buf) , " scl_inter = '%g'\n" , nim->scl_inter ) ; } if( nim->intent_code > 0 ){ sprintf( buf+strlen(buf) , " intent_code = '%d'\n", nim->intent_code ) ; sprintf( buf+strlen(buf) , " intent_code_name = '%s'\n" , nifti_intent_string(nim->intent_code) ) ; sprintf( buf+strlen(buf) , " intent_p1 = '%g'\n" , nim->intent_p1 ) ; sprintf( buf+strlen(buf) , " intent_p2 = '%g'\n" , nim->intent_p2 ) ; sprintf( buf+strlen(buf) , " intent_p3 = '%g'\n" , nim->intent_p3 ) ; if( nim->intent_name[0] != '\0' ){ ebuf = escapize_string(nim->intent_name) ; sprintf( buf+strlen(buf) , " intent_name = %s\n",ebuf) ; free(ebuf) ; } } if( nim->toffset != 0.0 ) sprintf( buf+strlen(buf) , " toffset = '%g'\n",nim->toffset ) ; if( nim->xyz_units > 0 ) sprintf( buf+strlen(buf) , " xyz_units = '%d'\n" " xyz_units_name = '%s'\n" , nim->xyz_units , nifti_units_string(nim->xyz_units) ) ; if( nim->time_units > 0 ) sprintf( buf+strlen(buf) , " time_units = '%d'\n" " time_units_name = '%s'\n" , nim->time_units , nifti_units_string(nim->time_units) ) ; if( nim->freq_dim > 0 ) sprintf( buf+strlen(buf) , " freq_dim = '%d'\n",nim->freq_dim ) ; if( nim->phase_dim > 0 ) sprintf( buf+strlen(buf) , " phase_dim = '%d'\n",nim->phase_dim ) ; if( nim->slice_dim > 0 ) sprintf( buf+strlen(buf) , " slice_dim = '%d'\n",nim->slice_dim ) ; if( nim->slice_code > 0 ) sprintf( buf+strlen(buf) , " slice_code = '%d'\n" " slice_code_name = '%s'\n" , nim->slice_code , nifti_slice_string(nim->slice_code) ) ; if( nim->slice_start >= 0 && nim->slice_end > nim->slice_start ) sprintf( buf+strlen(buf) , " slice_start = '%d'\n" " slice_end = '%d'\n" , nim->slice_start , nim->slice_end ) ; if( nim->slice_duration != 0.0 ) sprintf( buf+strlen(buf) , " slice_duration = '%g'\n", nim->slice_duration ) ; if( nim->descrip[0] != '\0' ){ ebuf = escapize_string(nim->descrip) ; sprintf( buf+strlen(buf) , " descrip = %s\n",ebuf) ; free(ebuf) ; } if( nim->aux_file[0] != '\0' ){ ebuf = escapize_string(nim->aux_file) ; sprintf( buf+strlen(buf) , " aux_file = %s\n",ebuf) ; free(ebuf) ; } if( nim->qform_code > 0 ){ int i,j,k ; sprintf( buf+strlen(buf) , " qform_code = '%d'\n" " qform_code_name = '%s'\n" " qto_xyz_matrix = '%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g'\n" , nim->qform_code , nifti_xform_string(nim->qform_code) , nim->qto_xyz.m[0][0] , nim->qto_xyz.m[0][1] , nim->qto_xyz.m[0][2] , nim->qto_xyz.m[0][3] , nim->qto_xyz.m[1][0] , nim->qto_xyz.m[1][1] , nim->qto_xyz.m[1][2] , nim->qto_xyz.m[1][3] , nim->qto_xyz.m[2][0] , nim->qto_xyz.m[2][1] , nim->qto_xyz.m[2][2] , nim->qto_xyz.m[2][3] , nim->qto_xyz.m[3][0] , nim->qto_xyz.m[3][1] , nim->qto_xyz.m[3][2] , nim->qto_xyz.m[3][3] ) ; sprintf( buf+strlen(buf) , " qto_ijk_matrix = '%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g'\n" , nim->qto_ijk.m[0][0] , nim->qto_ijk.m[0][1] , nim->qto_ijk.m[0][2] , nim->qto_ijk.m[0][3] , nim->qto_ijk.m[1][0] , nim->qto_ijk.m[1][1] , nim->qto_ijk.m[1][2] , nim->qto_ijk.m[1][3] , nim->qto_ijk.m[2][0] , nim->qto_ijk.m[2][1] , nim->qto_ijk.m[2][2] , nim->qto_ijk.m[2][3] , nim->qto_ijk.m[3][0] , nim->qto_ijk.m[3][1] , nim->qto_ijk.m[3][2] , nim->qto_ijk.m[3][3] ) ; sprintf( buf+strlen(buf) , " quatern_b = '%g'\n" " quatern_c = '%g'\n" " quatern_d = '%g'\n" " qoffset_x = '%g'\n" " qoffset_y = '%g'\n" " qoffset_z = '%g'\n" " qfac = '%g'\n" , nim->quatern_b , nim->quatern_c , nim->quatern_d , nim->qoffset_x , nim->qoffset_y , nim->qoffset_z , nim->qfac ) ; mat44_to_orientation( nim->qto_xyz , &i,&j,&k ) ; if( i > 0 && j > 0 && k > 0 ) sprintf( buf+strlen(buf) , " qform_i_orientation = '%s'\n" " qform_j_orientation = '%s'\n" " qform_k_orientation = '%s'\n" , nifti_orientation_string(i) , nifti_orientation_string(j) , nifti_orientation_string(k) ) ; } if( nim->sform_code > 0 ){ int i,j,k ; sprintf( buf+strlen(buf) , " sform_code = '%d'\n" " sform_code_name = '%s'\n" " sto_xyz_matrix = '%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g'\n" , nim->sform_code , nifti_xform_string(nim->sform_code) , nim->sto_xyz.m[0][0] , nim->sto_xyz.m[0][1] , nim->sto_xyz.m[0][2] , nim->sto_xyz.m[0][3] , nim->sto_xyz.m[1][0] , nim->sto_xyz.m[1][1] , nim->sto_xyz.m[1][2] , nim->sto_xyz.m[1][3] , nim->sto_xyz.m[2][0] , nim->sto_xyz.m[2][1] , nim->sto_xyz.m[2][2] , nim->sto_xyz.m[2][3] , nim->sto_xyz.m[3][0] , nim->sto_xyz.m[3][1] , nim->sto_xyz.m[3][2] , nim->sto_xyz.m[3][3] ) ; sprintf( buf+strlen(buf) , " sto_ijk matrix = '%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g'\n" , nim->sto_ijk.m[0][0] , nim->sto_ijk.m[0][1] , nim->sto_ijk.m[0][2] , nim->sto_ijk.m[0][3] , nim->sto_ijk.m[1][0] , nim->sto_ijk.m[1][1] , nim->sto_ijk.m[1][2] , nim->sto_ijk.m[1][3] , nim->sto_ijk.m[2][0] , nim->sto_ijk.m[2][1] , nim->sto_ijk.m[2][2] , nim->sto_ijk.m[2][3] , nim->sto_ijk.m[3][0] , nim->sto_ijk.m[3][1] , nim->sto_ijk.m[3][2] , nim->sto_ijk.m[3][3] ) ; mat44_to_orientation( nim->sto_xyz , &i,&j,&k ) ; if( i > 0 && j > 0 && k > 0 ) sprintf( buf+strlen(buf) , " sform_i_orientation = '%s'\n" " sform_j_orientation = '%s'\n" " sform_k_orientation = '%s'\n" , nifti_orientation_string(i) , nifti_orientation_string(j) , nifti_orientation_string(k) ) ; } sprintf( buf+strlen(buf) , " num_ext = '%d'\n", nim->num_ext ) ; sprintf( buf+strlen(buf) , "/>\n" ) ; /* XML-ish closer */ nbuf = strlen(buf) ; buf = (char *)realloc((void *)buf, nbuf+1); /* cut back to proper length */ return buf ; }

void nifti_image_unload (  nifti_image *    nim )

Definition at line 2954 of file nifti1_io.c. References nifti_image::data, and free. { if( nim != NULL && nim->data != NULL ){ free(nim->data) ; nim->data = NULL ; } return ; }

void nifti_image_write (  nifti_image *    nim )

Definition at line 3436 of file nifti1_io.c. References nifti_image_write_hdr_img(). Referenced by main(). { nifti_image_write_hdr_img(nim,1,"wb"); }

znzFile nifti_image_write_hdr_img (  nifti_image *    nim, int    write_data, char *    opts )  [static]

Referenced by nifti_image_write().

znzFile nifti_image_write_hdr_img2 (  nifti_image *    nim, int    write_data, char *    opts, znzFile *    imgfile )  [static]

Referenced by nifti_image_write_hdr_img().

char* nifti_intent_string (  int    ii )

Definition at line 724 of file nifti1_io.c. References NIFTI_INTENT_BETA, NIFTI_INTENT_BINOM, NIFTI_INTENT_CHI, NIFTI_INTENT_CHISQ, NIFTI_INTENT_CHISQ_NONC, NIFTI_INTENT_CORREL, NIFTI_INTENT_DIMLESS, NIFTI_INTENT_DISPVECT, NIFTI_INTENT_ESTIMATE, NIFTI_INTENT_EXTVAL, NIFTI_INTENT_FTEST, NIFTI_INTENT_FTEST_NONC, NIFTI_INTENT_GAMMA, NIFTI_INTENT_GENMATRIX, NIFTI_INTENT_INVGAUSS, NIFTI_INTENT_LABEL, NIFTI_INTENT_LAPLACE, NIFTI_INTENT_LOG10PVAL, NIFTI_INTENT_LOGISTIC, NIFTI_INTENT_LOGPVAL, NIFTI_INTENT_NEURONAME, NIFTI_INTENT_NORMAL, NIFTI_INTENT_POINTSET, NIFTI_INTENT_POISSON, NIFTI_INTENT_PVAL, NIFTI_INTENT_QUATERNION, NIFTI_INTENT_SYMMATRIX, NIFTI_INTENT_TRIANGLE, NIFTI_INTENT_TTEST, NIFTI_INTENT_TTEST_NONC, NIFTI_INTENT_UNIFORM, NIFTI_INTENT_VECTOR, NIFTI_INTENT_WEIBULL, and NIFTI_INTENT_ZSCORE. Referenced by nifti_image_to_ascii(). { switch( ii ){ case NIFTI_INTENT_CORREL: return "Correlation statistic" ; case NIFTI_INTENT_TTEST: return "T-statistic" ; case NIFTI_INTENT_FTEST: return "F-statistic" ; case NIFTI_INTENT_ZSCORE: return "Z-score" ; case NIFTI_INTENT_CHISQ: return "Chi-squared distribution" ; case NIFTI_INTENT_BETA: return "Beta distribution" ; case NIFTI_INTENT_BINOM: return "Binomial distribution" ; case NIFTI_INTENT_GAMMA: return "Gamma distribution" ; case NIFTI_INTENT_POISSON: return "Poisson distribution" ; case NIFTI_INTENT_NORMAL: return "Normal distribution" ; case NIFTI_INTENT_FTEST_NONC: return "F-statistic noncentral" ; case NIFTI_INTENT_CHISQ_NONC: return "Chi-squared noncentral" ; case NIFTI_INTENT_LOGISTIC: return "Logistic distribution" ; case NIFTI_INTENT_LAPLACE: return "Laplace distribution" ; case NIFTI_INTENT_UNIFORM: return "Uniform distribition" ; case NIFTI_INTENT_TTEST_NONC: return "T-statistic noncentral" ; case NIFTI_INTENT_WEIBULL: return "Weibull distribution" ; case NIFTI_INTENT_CHI: return "Chi distribution" ; case NIFTI_INTENT_INVGAUSS: return "Inverse Gaussian distribution" ; case NIFTI_INTENT_EXTVAL: return "Extreme Value distribution" ; case NIFTI_INTENT_PVAL: return "P-value" ; case NIFTI_INTENT_LOGPVAL: return "Log P-value" ; case NIFTI_INTENT_LOG10PVAL: return "Log10 P-value" ; case NIFTI_INTENT_ESTIMATE: return "Estimate" ; case NIFTI_INTENT_LABEL: return "Label index" ; case NIFTI_INTENT_NEURONAME: return "NeuroNames index" ; case NIFTI_INTENT_GENMATRIX: return "General matrix" ; case NIFTI_INTENT_SYMMATRIX: return "Symmetric matrix" ; case NIFTI_INTENT_DISPVECT: return "Displacement vector" ; case NIFTI_INTENT_VECTOR: return "Vector" ; case NIFTI_INTENT_POINTSET: return "Pointset" ; case NIFTI_INTENT_TRIANGLE: return "Triangle" ; case NIFTI_INTENT_QUATERNION: return "Quaternion" ; case NIFTI_INTENT_DIMLESS: return "Dimensionless number" ; } return "Unknown" ; }

int nifti_is_gzfile (  char *    fname )  [static]

Referenced by nifti_image_load_prep(), and nifti_image_read().

int nifti_is_inttype (  int    dt )

Definition at line 658 of file nifti1_io.c. References dt, DT_BINARY, DT_COMPLEX128, DT_COMPLEX256, DT_COMPLEX64, DT_FLOAT128, DT_FLOAT32, DT_FLOAT64, DT_INT16, DT_INT32, DT_INT64, DT_INT8, DT_RGB24, DT_UINT16, DT_UINT32, DT_UINT64, DT_UINT8, and DT_UNKNOWN. { switch( dt ){ case DT_UNKNOWN: return 0 ; case DT_BINARY: return 0 ; case DT_INT8: return 1 ; case DT_UINT8: return 1 ; case DT_INT16: return 1 ; case DT_UINT16: return 1 ; case DT_INT32: return 1 ; case DT_UINT32: return 1 ; case DT_INT64: return 1 ; case DT_UINT64: return 1 ; case DT_FLOAT32: return 0 ; case DT_FLOAT64: return 0 ; case DT_FLOAT128: return 0 ; case DT_COMPLEX64: return 0 ; case DT_COMPLEX128: return 0 ; case DT_COMPLEX256: return 0 ; case DT_RGB24: return 1 ; } return 0 ; }

int nifti_load_NBL_bricks (  nifti_image *    nim, int *    slist, int *    sindex, nifti_brick_list *    NBL, znzFile    fp )  [static]

Referenced by nifti_image_load_bricks().

char* nifti_makebasename (  char *    fname )

Definition at line 1738 of file nifti1_io.c. References nifti_find_file_extension(), and nifti_strdup(). Referenced by main(), nifti_findhdrname(), and nifti_findimgname(). { /* NB: it allocates memory for basename which should be freed when no longer required */ char *basename, *ext; basename=nifti_strdup(fname); ext = nifti_find_file_extension(basename); if ( ext ) *ext = '\0'; /* clear out extension */ return basename; /* in either case */ }

char* nifti_orientation_string (  int    ii )

Definition at line 789 of file nifti1_io.c. References NIFTI_A2P, NIFTI_I2S, NIFTI_L2R, NIFTI_P2A, NIFTI_R2L, and NIFTI_S2I. Referenced by nifti_disp_matrix_orient(), and nifti_image_to_ascii(). { switch( ii ){ case NIFTI_L2R: return "Left-to-Right" ; case NIFTI_R2L: return "Right-to-Left" ; case NIFTI_P2A: return "Posterior-to-Anterior" ; case NIFTI_A2P: return "Anterior-to-Posterior" ; case NIFTI_I2S: return "Inferior-to-Superior" ; case NIFTI_S2I: return "Superior-to-Inferior" ; } return "Unknown" ; }

size_t nifti_read_buffer (  znzFile    fp, void *    datatptr, size_t    ntot, nifti_image *    nim )  [static]

Referenced by nifti_image_load(), and nifti_load_NBL_bricks().

int nifti_read_extensions (  nifti_image *    nim, znzFile    fp )  [static]

Referenced by nifti_image_read().

int nifti_read_next_extension (  nifti1_extension *    nex, nifti_image *    nim, znzFile    fp )  [static]

Referenced by nifti_read_extensions().

void nifti_set_debug_level (  int    level )

Definition at line 1754 of file nifti1_io.c. References gni_debug. Referenced by main(). { gni_debug = level; }

void nifti_set_iname_offset (  nifti_image *    nim )

Definition at line 3246 of file nifti1_io.c. References gni_debug, nifti_image::iname_offset, nifti_extension_size(), nifti_image::nifti_type, and offset. Referenced by nifti_image_write_hdr_img2(). { int offset = nifti_extension_size(nim); /* init based on extensions */ switch( nim->nifti_type ){ default: /* writing into 2 files */ /* we only write files with 0 offset in the 2 file format */ nim->iname_offset = 0 ; break ; case 1: /* NIFTI-1 single binary file */ if (nim->iname_offset < sizeof(struct nifti_1_header)) { offset += sizeof(struct nifti_1_header) + 4 ; /* be sure offset is aligned to a 16 byte word boundary */ if ( ( offset % 16 ) != 0 ) offset = ((offset + 0xf) & ~0xf); if( nim->iname_offset != offset && gni_debug > 1 ){ fprintf(stderr,"+d changing offset from %d to %d\n", nim->iname_offset, offset); } nim->iname_offset = offset; } break ; /* non-standard case: */ case 3:{ /* NIFTI-1 ASCII header + binary data (single file) */ nim->iname_offset = -1 ; /* compute offset from filesize */ } } }

nifti_image* nifti_simple_init_nim (    )

Definition at line 3085 of file nifti1_io.c. References nifti_1_header::bitpix, nifti_1_header::datatype, nifti_1_header::dim, nifti_image::fname, nifti_image::iname, nifti_convert_nhdr2nim(), nifti_datatype_sizes(), NIFTI_TYPE_FLOAT32, nifti_1_header::pixdim, nifti_1_header::regular, and nifti_1_header::sizeof_hdr. { nifti_image *nim; struct nifti_1_header nhdr; int nbyper, swapsize; memset(&nhdr,0,sizeof(nhdr)) ; /* zero out header, to be safe */ /** load the ANALYZE-7.5 generic parts of the header struct **/ nhdr.sizeof_hdr = sizeof(nhdr) ; nhdr.regular = 'r' ; /* for some stupid reason */ nhdr.dim[0] = 3 ; nhdr.dim[1] = 1 ; nhdr.dim[2] = 1 ; nhdr.dim[3] = 1 ; nhdr.dim[4] = 0 ; nhdr.pixdim[0] = 0.0 ; nhdr.pixdim[1] = 1.0 ; nhdr.pixdim[2] = 1.0 ; nhdr.pixdim[3] = 1.0 ; nhdr.datatype = NIFTI_TYPE_FLOAT32 ; nifti_datatype_sizes( nhdr.datatype , &nbyper, &swapsize ); nhdr.bitpix = 8 * nbyper ; nim = nifti_convert_nhdr2nim(nhdr,NULL); nim->fname = NULL; nim->iname = NULL; return nim; }

char* nifti_slice_string (  int    ss )

Definition at line 773 of file nifti1_io.c. References NIFTI_SLICE_ALT_DEC, NIFTI_SLICE_ALT_INC, NIFTI_SLICE_SEQ_DEC, and NIFTI_SLICE_SEQ_INC. Referenced by nifti_image_to_ascii(). { switch( ss ){ case NIFTI_SLICE_SEQ_INC: return "sequential_increasing" ; case NIFTI_SLICE_SEQ_DEC: return "sequential_decreasing" ; case NIFTI_SLICE_ALT_INC: return "alternating_increasing" ; case NIFTI_SLICE_ALT_DEC: return "alternating_decreasing" ; } return "Unknown" ; }

char* nifti_units_string (  int    uu )

Definition at line 687 of file nifti1_io.c. References NIFTI_UNITS_HZ, NIFTI_UNITS_METER, NIFTI_UNITS_MICRON, NIFTI_UNITS_MM, NIFTI_UNITS_MSEC, NIFTI_UNITS_PPM, NIFTI_UNITS_RADS, NIFTI_UNITS_SEC, and NIFTI_UNITS_USEC. Referenced by nifti_image_to_ascii(). { switch( uu ){ case NIFTI_UNITS_METER: return "m" ; case NIFTI_UNITS_MM: return "mm" ; case NIFTI_UNITS_MICRON: return "um" ; case NIFTI_UNITS_SEC: return "s" ; case NIFTI_UNITS_MSEC: return "ms" ; case NIFTI_UNITS_USEC: return "us" ; case NIFTI_UNITS_HZ: return "Hz" ; case NIFTI_UNITS_PPM: return "ppm" ; case NIFTI_UNITS_RADS: return "rad/s" ; } return "Unknown" ; }

int nifti_valid_extension (  nifti_image *    nim, int    size, int    code )  [static]

Referenced by nifti_read_next_extension().

int nifti_validfilename (  char *    fname )

Definition at line 1672 of file nifti1_io.c. References gni_debug, and nifti_find_file_extension(). Referenced by is_nifti_file(), nifti_findhdrname(), nifti_findimgname(), and nifti_image_write_hdr_img2(). { char * ext; int len; /** check input file(s) for sanity **/ if( fname == NULL || *fname == '\0' ) return 0; len = strlen(fname); if( len < 1 ) { if ( gni_debug > 0 ) fprintf(stderr,"-- empty filename in nifti_validfilename()\n"); return 0; } ext = nifti_find_file_extension((char *)fname); if ( ext && ext == fname ) { /* then no filename prefix */ if ( gni_debug > 0 ) fprintf(stderr,"-- no prefix for filename '%s'\n", fname); return 0; } return 1; }

void nifti_write_all_data (  znzFile    fp, nifti_image *    nim )  [static]

Referenced by nifti_image_write_hdr_img2().

size_t nifti_write_buffer (  znzFile    fp, void *    buffer, size_t    numbytes )  [static]

Referenced by nifti_write_all_data(), and nifti_write_extensions().

int nifti_write_extensions (  znzFile    fp, nifti_image *    nim )  [static]

Referenced by nifti_image_write_hdr_img2().

char* nifti_xform_string (  int    xx )

Definition at line 708 of file nifti1_io.c. References NIFTI_XFORM_ALIGNED_ANAT, NIFTI_XFORM_MNI_152, NIFTI_XFORM_SCANNER_ANAT, and NIFTI_XFORM_TALAIRACH. Referenced by nifti_image_to_ascii(). { switch( xx ){ case NIFTI_XFORM_SCANNER_ANAT: return "Scanner Anat" ; case NIFTI_XFORM_ALIGNED_ANAT: return "Aligned Anat" ; case NIFTI_XFORM_TALAIRACH: return "Talairach" ; case NIFTI_XFORM_MNI_152: return "MNI_152" ; } return "Unknown" ; }

int print_hex_vals (  char *    data, int    nbytes, FILE *    fp )

Definition at line 1934 of file nifti1_io.c. References c. Referenced by disp_nifti_1_header(). { int c; if ( !data || nbytes < 1 || !fp ) return -1; fputs("0x", fp); for ( c = 0; c < nbytes; c++ ) fprintf(fp, " %x", data[c]); return 0; }

mat44 quatern_to_mat44 (  float    qb, float    qc, float    qd, float    qx, float    qy, float    qz, float    dx, float    dy, float    dz, float    qfac )

Definition at line 851 of file nifti1_io.c. References a, c, l, and mat44::m. Referenced by nifti_convert_nhdr2nim(), and nifti_image_from_ascii(). { mat44 R ; double a,b=qb,c=qc,d=qd , xd,yd,zd ; /* last row is always [ 0 0 0 1 ] */ R.m[3][0]=R.m[3][1]=R.m[3][2] = 0.0 ; R.m[3][3]= 1.0 ; /* compute a parameter from b,c,d */ a = 1.0l - (b*b + c*c + d*d) ; if( a < 1.e-7l ){ /* special case */ a = 1.0l / sqrt(b*b+c*c+d*d) ; b *= a ; c *= a ; d *= a ; /* normalize (b,c,d) vector */ a = 0.0l ; /* a = 0 ==> 180 degree rotation */ } else{ a = sqrt(a) ; /* angle = 2*arccos(a) */ } /* load rotation matrix, including scaling factors for voxel sizes */ xd = (dx > 0.0) ? dx : 1.0l ; /* make sure are positive */ yd = (dy > 0.0) ? dy : 1.0l ; zd = (dz > 0.0) ? dz : 1.0l ; if( qfac < 0.0 ) zd = -zd ; /* left handedness? */ R.m[0][0] = (a*a+b*b-c*c-d*d) * xd ; R.m[0][1] = 2.0l * (b*c-a*d ) * yd ; R.m[0][2] = 2.0l * (b*d+a*c ) * zd ; R.m[1][0] = 2.0l * (b*c+a*d ) * xd ; R.m[1][1] = (a*a+c*c-b*b-d*d) * yd ; R.m[1][2] = 2.0l * (c*d-a*b ) * zd ; R.m[2][0] = 2.0l * (b*d-a*c ) * xd ; R.m[2][1] = 2.0l * (c*d+a*b ) * yd ; R.m[2][2] = (a*a+d*d-c*c-b*b) * zd ; /* load offsets */ R.m[0][3] = qx ; R.m[1][3] = qy ; R.m[2][3] = qz ; return R ; }

int short_order (  void    )

Strings that we don't control (filenames, etc.) that might contain "weird" characters (like quotes) are "escaped": A few special characters are replaced by XML-style escapes, using the function escapize_string(). On input, function unescape_string() reverses this process. The result is that the NIFTI ASCII-format header is XML-compliant. * Definition at line 3859 of file nifti1_io.c. Referenced by nifti_convert_nhdr2nim(), nifti_image_from_ascii(), nifti_image_write_hdr_img2(), nifti_read_buffer(), nifti_read_next_extension(), and nifti_write_all_data(). { union { unsigned char bb[2] ; short ss ; } fred ; fred.bb[0] = 1 ; fred.bb[1] = 0 ; return (fred.ss == 1) ? LSB_FIRST : MSB_FIRST ; }

void swap_16bytes (  int    n, void *    ar )

Definition at line 1538 of file nifti1_io.c. References sixteenbytes::A, sixteenbytes::a, sixteenbytes::B, sixteenbytes::b, sixteenbytes::C, sixteenbytes::c, sixteenbytes::D, sixteenbytes::d, sixteenbytes::E, sixteenbytes::e, sixteenbytes::F, sixteenbytes::f, sixteenbytes::G, sixteenbytes::g, sixteenbytes::H, sixteenbytes::h, and tt. Referenced by swap_Nbytes(). { register int ii ; register sixteenbytes *tb = (sixteenbytes *)ar ; register unsigned char tt ; for( ii=0 ; ii < n ; ii++ ){ tt = tb[ii].a ; tb[ii].a = tb[ii].A ; tb[ii].A = tt ; tt = tb[ii].b ; tb[ii].b = tb[ii].B ; tb[ii].B = tt ; tt = tb[ii].c ; tb[ii].c = tb[ii].C ; tb[ii].C = tt ; tt = tb[ii].d ; tb[ii].d = tb[ii].D ; tb[ii].D = tt ; tt = tb[ii].e ; tb[ii].e = tb[ii].E ; tb[ii].E = tt ; tt = tb[ii].f ; tb[ii].f = tb[ii].F ; tb[ii].F = tt ; tt = tb[ii].g ; tb[ii].g = tb[ii].G ; tb[ii].G = tt ; tt = tb[ii].h ; tb[ii].h = tb[ii].H ; tb[ii].H = tt ; } return ; }

void swap_2bytes (  int    n, void *    ar )

Definition at line 1485 of file nifti1_io.c. References twobytes::a, twobytes::b, and tt. Referenced by swap_Nbytes(), and swap_nifti_header(). { register int ii ; register twobytes *tb = (twobytes *)ar ; register unsigned char tt ; for( ii=0 ; ii < n ; ii++ ){ tt = tb[ii].a ; tb[ii].a = tb[ii].b ; tb[ii].b = tt ; } return ; }

void swap_4bytes (  int    n, void *    ar )

Definition at line 1501 of file nifti1_io.c. References fourbytes::a, fourbytes::b, fourbytes::c, fourbytes::d, and tt. Referenced by nifti_read_next_extension(), swap_Nbytes(), and swap_nifti_header(). { register int ii ; register fourbytes *tb = (fourbytes *)ar ; register unsigned char tt ; for( ii=0 ; ii < n ; ii++ ){ tt = tb[ii].a ; tb[ii].a = tb[ii].d ; tb[ii].d = tt ; tt = tb[ii].b ; tb[ii].b = tb[ii].c ; tb[ii].c = tt ; } return ; }

void swap_8bytes (  int    n, void *    ar )

Definition at line 1518 of file nifti1_io.c. References eightbytes::A, eightbytes::a, eightbytes::B, eightbytes::b, eightbytes::C, eightbytes::c, eightbytes::D, eightbytes::d, and tt. Referenced by swap_Nbytes(). { register int ii ; register eightbytes *tb = (eightbytes *)ar ; register unsigned char tt ; for( ii=0 ; ii < n ; ii++ ){ tt = tb[ii].a ; tb[ii].a = tb[ii].A ; tb[ii].A = tt ; tt = tb[ii].b ; tb[ii].b = tb[ii].B ; tb[ii].B = tt ; tt = tb[ii].c ; tb[ii].c = tb[ii].C ; tb[ii].C = tt ; tt = tb[ii].d ; tb[ii].d = tb[ii].D ; tb[ii].D = tt ; } return ; }

void swap_Nbytes (  int    n, int    siz, void *    ar )

Definition at line 1560 of file nifti1_io.c. References swap_16bytes(), swap_2bytes(), swap_4bytes(), and swap_8bytes(). Referenced by nifti_read_buffer(), and THD_load_nifti(). { switch( siz ){ case 2: swap_2bytes ( n , ar ) ; break ; case 4: swap_4bytes ( n , ar ) ; break ; case 8: swap_8bytes ( n , ar ) ; break ; case 16: swap_16bytes( n , ar ) ; break ; } return ; }

void swap_nifti_header (  struct nifti_1_header *    h, int    is_nifti )

Definition at line 1578 of file nifti1_io.c. References nifti_1_header::bitpix, nifti_1_header::cal_max, nifti_1_header::cal_min, nifti_1_header::datatype, nifti_1_header::dim, nifti_1_header::extents, nifti_1_header::glmax, nifti_1_header::glmin, nifti_1_header::intent_code, nifti_1_header::intent_p1, nifti_1_header::intent_p2, nifti_1_header::intent_p3, nifti_1_header::pixdim, nifti_1_header::qform_code, nifti_1_header::qoffset_x, nifti_1_header::qoffset_y, nifti_1_header::qoffset_z, nifti_1_header::quatern_b, nifti_1_header::quatern_c, nifti_1_header::quatern_d, nifti_1_header::scl_inter, nifti_1_header::scl_slope, nifti_1_header::session_error, nifti_1_header::sform_code, nifti_1_header::sizeof_hdr, nifti_1_header::slice_duration, nifti_1_header::slice_end, nifti_1_header::slice_start, nifti_1_header::srow_x, nifti_1_header::srow_y, nifti_1_header::srow_z, swap_2, swap_2bytes(), swap_4, swap_4bytes(), nifti_1_header::toffset, and nifti_1_header::vox_offset. Referenced by nifti_convert_nhdr2nim(). { #if 0 /* ANALYZE fields not used by this software */ swap_4(h->sizeof_hdr) ; swap_4(h->extents) ; swap_2(h->session_error) ; swap_4(h->compressed) ; swap_4(h->glmax) ; swap_4(h->glmin) ; #endif /* this stuff is always present, for ANALYZE and NIFTI */ swap_2bytes( 8 , h->dim ) ; swap_4bytes( 8 , h->pixdim ) ; swap_2(h->datatype) ; swap_2(h->bitpix) ; swap_4(h->vox_offset); swap_4(h->cal_max); swap_4(h->cal_min); /* this stuff is NIFTI specific */ if( is_nifti ){ swap_2(h->qform_code); swap_2(h->sform_code); swap_4(h->quatern_b); swap_4(h->quatern_c); swap_4(h->quatern_d); swap_4(h->qoffset_x); swap_4(h->qoffset_y); swap_4(h->qoffset_z); swap_4(h->intent_p1); swap_4(h->intent_p2); swap_4(h->intent_p3); swap_4(h->scl_slope); swap_4(h->scl_inter); swap_4bytes(4,h->srow_x); swap_4bytes(4,h->srow_y); swap_4bytes(4,h->srow_z); swap_2(h->intent_code); swap_4(h->toffset); swap_2(h->slice_start); swap_2(h->slice_end); swap_4(h->slice_duration); } return ; }

int unescape_string (  char *    str )

Definition at line 3476 of file nifti1_io.c. References CR, and LF. Referenced by nifti_image_from_ascii(). { int ii,jj , nn,ll ; if( str == NULL ) return 0 ; /* no string? */ ll = strlen(str) ; if( ll == 0 ) return 0 ; /* scan for escapes: &something; */ for( ii=jj=nn=0 ; ii<ll ; ii++,jj++ ){ /* scan at ii; results go in at jj */ if( str[ii] == '&' ){ /* start of escape? */ if( ii+3 < ll && /* &lt; */ str[ii+1] == 'l' && str[ii+2] == 't' && str[ii+3] == ';' ){ str[jj] = '<' ; ii += 3 ; nn++ ; } else if( ii+3 < ll && /* &gt; */ str[ii+1] == 'g' && str[ii+2] == 't' && str[ii+3] == ';' ){ str[jj] = '>' ; ii += 3 ; nn++ ; } else if( ii+5 < ll && /* &quot; */ str[ii+1] == 'q' && str[ii+2] == 'u' && str[ii+3] == 'o' && str[ii+4] == 't' && str[ii+5] == ';' ){ str[jj] = '"' ; ii += 5 ; nn++ ; } else if( ii+5 < ll && /* &apos; */ str[ii+1] == 'a' && str[ii+2] == 'p' && str[ii+3] == 'o' && str[ii+4] == 's' && str[ii+5] == ';' ){ str[jj] = '\'' ; ii += 5 ; nn++ ; } else if( ii+4 < ll && /* &amp; */ str[ii+1] == 'a' && str[ii+2] == 'm' && str[ii+3] == 'p' && str[ii+4] == ';' ){ str[jj] = '&' ; ii += 4 ; nn++ ; } /* although the comments above don't mention it, we also look for XML style numeric escapes of the forms &#32; (decimal) and &#xfd; (hex) */ else if( ii+3 < ll && str[ii+1] == '#' && isdigit(str[ii+2]) ){ /* &#dec; */ unsigned int val='?' ; int kk=ii+3 ; while( kk < ll && kk != ';' ) kk++ ; sscanf( str+ii+2 , "%u" , &val ) ; str[jj] = (char) val ; ii = kk ; nn++ ; } else if( ii+4 < ll && str[ii+1] == '#' && str[ii+2] == 'x' && isxdigit(str[ii+3]) ){ /* &#hex; */ unsigned int val='?' ; int kk=ii+4 ; while( kk < ll && kk != ';' ) kk++ ; sscanf( str+ii+3 , "%x" , &val ) ; str[jj] = (char) val ; ii = kk ; nn++ ; } /* didn't start a recognized escape, so just copy as normal */ else if( jj < ii ){ str[jj] = str[ii] ; } } else if( str[ii] == CR ) { /* is a carriage return */ if( str[ii+1] == LF ){ str[jj] = LF ; ii++ ; nn++ ; } /* CR LF */ else { str[jj] = LF ; ; nn++ ; } /* CR only */ } else { /* is a normal character, just copy to output */ if( jj < ii ){ str[jj] = str[ii] ; } } /* at this point, ii=index of last character used up in scan jj=index of last character written to (jj <= ii) */ } if( jj < ll ) str[jj] = '\0' ; /* end string properly */ return nn ; }

void update_nifti_image_for_brick_list (  nifti_image *    nim, int    nbricks )  [static]

Referenced by nifti_image_read_bricks().

int valid_nifti_brick_list (  nifti_image *    nim, int    nbricks, int *    blist, int    disp_error )

Definition at line 541 of file nifti1_io.c. References c, nifti_image::dim, gni_debug, int_force_positive(), and nsubs(). Referenced by nifti_image_load_bricks(). { int c, nsubs; if( !nim ){ if( disp_error || gni_debug > 0 ) fprintf(stderr,"** valid_nifti_brick_list: missing nifti image\n"); return 0; } if( nbricks <= 0 || !blist ){ if( disp_error || gni_debug > 1 ) fprintf(stderr,"** valid_nifti_brick_list: no brick list to check\n"); return 0; } /* nsubs sub-brick is nt*nu*nv*nw */ nsubs = nim->dim[4] * nim->dim[5] * nim->dim[6] * nim->dim[7]; if( nsubs <= 0 ){ fprintf(stderr,"** VNBL warning: bad dim list (%d,%d,%d,%d)\n", nim->dim[4], nim->dim[5], nim->dim[6], nim->dim[7]); int_force_positive(nim->dim+4, 4); nsubs = nim->dim[4] * nim->dim[5] * nim->dim[6] * nim->dim[7]; } for( c = 0; c < nbricks; c++ ) if( (blist[c] < 0) || (blist[c] >= nsubs) ){ if( disp_error || gni_debug > 1 ) fprintf(stderr, "-d ** bad sub-brick chooser %d (#%d), valid range is [0,%d]\n", blist[c], c, nsubs-1); return 0; } return 1; /* all is well */ }

                               

This site conforms to the following standards:

• 
• 
• 
• 
•