TPCCLIB
Loading...
Searching...
No Matches
analyze.c
Go to the documentation of this file.
1
7/*****************************************************************************/
8int ANALYZE_TEST;
9/*****************************************************************************/
10#include "libtpcimgio.h"
11/*****************************************************************************/
12
13/*****************************************************************************/
17int anaExists(
19 const char *dbname
20) {
21 char temp[FILENAME_MAX];
22
23 if(dbname==NULL || strlen(dbname)==0) return(0);
24 /* Header file? */
25 strlcpy(temp, dbname, FILENAME_MAX); strlcat(temp, ".hdr", FILENAME_MAX);
26 if(access(temp, 0) == -1) return(0);
27 /* Image data? */
28 strlcpy(temp, dbname, FILENAME_MAX); strlcat(temp, ".img", FILENAME_MAX);
29 if(access(temp, 0) == -1) return(0);
30 /* SIF? */
31 strlcat(temp, ".sif", FILENAME_MAX); if(access(temp, 0) != -1) return(2);
32 strlcpy(temp, dbname, FILENAME_MAX); strlcat(temp, ".sif", FILENAME_MAX);
33 if(access(temp, 0) != -1) return(2);
34 return(1);
35}
36/*****************************************************************************/
37
38/*****************************************************************************/
43int anaExistsNew(
46 const char *filename,
50 char *hdrfile,
54 char *imgfile,
58 char *siffile
59) {
60 char *cptr, basefile[FILENAME_MAX], temp[FILENAME_MAX];
61 ANALYZE_DSR h;
62 int ret;
63
64 if(filename==NULL || strlen(filename)==0) return(0);
65 if(ANALYZE_TEST>1) printf("\nanaExistsNew(%s, *str, *str, *str)\n", filename);
66
67 /* Construct the base file name wo extensions */
68 strcpy(basefile, filename);
69 cptr=strrchr(basefile, '.');
70 if(cptr!=NULL) {
71 if(strncasecmp(cptr, ".HDR", 4)==0 || strncasecmp(cptr, ".IMG", 4)==0 )
72 *cptr=(char)0;
73 }
74 cptr=strrchr(basefile, '.');
75 if(cptr!=NULL) {
76 if(strncasecmp(cptr, ".IMG", 4)==0 )
77 *cptr=(char)0;
78 }
79 if(ANALYZE_TEST>2) printf("\n basefile := %s\n", basefile);
80
81 /* Header file exists? */
82 strcpy(temp, basefile); strcat(temp, ".hdr");
83 if(access(temp, 0) == -1) {
84 strcpy(temp, basefile); strcat(temp, ".img.hdr");
85 if(access(temp, 0) == -1) {
86 if(ANALYZE_TEST) printf("\n hdr file not found or accessible.\n");
87 return(0);
88 }
89 }
90 /* Is this Analyze header file? */
91 if((ret=anaReadHeader(temp, &h))!=0) {
92 if(ANALYZE_TEST) printf("\n %s was not identified as Analyze header file (%d).\n", temp, ret);
93 return(0);
94 }
95 /* Preserve header filename */
96 if(hdrfile!=NULL) strcpy(hdrfile, temp);
97
98 /* Image file exists? */
99 strcpy(temp, basefile); strcat(temp, ".img");
100 if(access(temp, 0) == -1) {
101 if(ANALYZE_TEST) printf("\n %s not found or accessible.\n", temp);
102 return(0);
103 }
104 /* Preserve image filename */
105 if(imgfile!=NULL) strcpy(imgfile, temp);
106
107 /* SIF exists? */
108 strcpy(temp, basefile); strcat(temp, ".sif");
109 if(access(temp, 0) == -1) {
110 strcpy(temp, basefile); strcat(temp, ".img.sif");
111 if(access(temp, 0) == -1) {
112 if(ANALYZE_TEST) printf("\n SIF not found or accessible.\n");
113 if(siffile!=NULL) strcpy(siffile, "");
114 return(1); // but otherwise ok
115 }
116 }
117 /* Preserve SIF filename */
118 if(siffile!=NULL) strcpy(siffile, temp);
119 return(2);
120}
121/*****************************************************************************/
122
123/*****************************************************************************/
129int anaReadHeader(
131 char *filename,
133 ANALYZE_DSR *h
134) {
135 unsigned char buf1[ANALYZE_HEADER_KEY_SIZE];
136 unsigned char buf2[ANALYZE_HEADER_IMGDIM_SIZE];
137 unsigned char buf3[ANALYZE_HEADER_HISTORY_SIZE];
138 int little; /* 1 if current platform is little endian (i386), else 0 */
139 FILE *fp;
140 int ret, nr=0, s1, s2, same_order;
141
142 if(ANALYZE_TEST) printf("anaReadHeader(%s, *dsr)\n", filename);
143
144 /* Check arguments */
145 if(strlen(filename)<1 || h==NULL) return(1);
146 little=little_endian();
147 /* Open file */
148 fp=fopen(filename, "rb"); if(fp==NULL) return(2);
149 /* Get file size */
150 nr=0; while((ret=fgetc(fp))!=EOF) nr++; rewind(fp);
151 if(nr<1) {fclose(fp); return(3);}
152 /* Read Analyze header key */
153 if(fread(buf1, ANALYZE_HEADER_KEY_SIZE, 1, fp)<1) return(3);
154 /* Read Analyze header image dimension */
155 if(fread(buf2, ANALYZE_HEADER_IMGDIM_SIZE, 1, fp)<1) return(3);
156 /* Read Analyze header image data history */
157 memset(buf3, 0, sizeof(ANALYZE_HEADER_HISTORY));
158 ret=fread(buf3, ANALYZE_HEADER_HISTORY_SIZE, 1, fp);
159 if(ANALYZE_TEST>1 && ret<1) printf(" complete data_history not found.\n");
160 /* Close file */
161 fclose(fp);
162 /* Compare file size from header contents to the calculated value */
163 /* to determine whether Analyze file is in little or big endian */
164 memcpy(&s1, buf1+0, 4); s2=s1; swawbip(&s2, 4);
165 if(abs(s1-nr)<abs(s2-nr)) same_order=1; else same_order=0;
166 if(ANALYZE_TEST>1) printf("same byte order: %d (s1=%d s2=%d nr=%d)\n",
167 same_order, s1, s2, nr);
168 if(same_order) h->little=little;
169 else {if(little) h->little=0; else h->little=1;}
170
171 /* Set key header structure contents */
172 if(!same_order) swawbip(buf1+0, 4);
173 memcpy(&h->hk.sizeof_hdr, buf1+0, 4);
174 memcpy(h->hk.data_type, buf1+4, 10);
175 memcpy(h->hk.db_name, buf1+14, 18);
176 if(!same_order) swawbip(buf1+32, 4);
177 memcpy(&h->hk.extents, buf1+32, 4);
178 if(!same_order) swabip(buf1+36, 2);
179 memcpy(&h->hk.session_error, buf1+36, 2);
180 memcpy(&h->hk.regular, buf1+38, 1);
181 memcpy(&h->hk.hkey_un0, buf1+39, 1);
182
183 /* Set image dimension header structure contents */
184 if(!same_order) swabip(buf2+0, 16);
185 memcpy(h->dime.dim, buf2+0, 16);
186 if(!same_order) swabip(buf2+16, 2);
187 memcpy(&h->dime.unused8, buf2+16, 2);
188 if(!same_order) swabip(buf2+18, 2);
189 memcpy(&h->dime.unused9, buf2+18, 2);
190 if(!same_order) swabip(buf2+20, 2);
191 memcpy(&h->dime.unused10, buf2+20, 2);
192 if(!same_order) swabip(buf2+22, 2);
193 memcpy(&h->dime.unused11, buf2+22, 2);
194 if(!same_order) swabip(buf2+24, 2);
195 memcpy(&h->dime.unused12, buf2+24, 2);
196 if(!same_order) swabip(buf2+26, 2);
197 memcpy(&h->dime.unused13, buf2+26, 2);
198 if(!same_order) swabip(buf2+28, 2);
199 memcpy(&h->dime.unused14, buf2+28, 2);
200 if(!same_order) swabip(buf2+30, 2);
201 memcpy(&h->dime.datatype, buf2+30, 2);
202 if(!same_order) swabip(buf2+32, 2);
203 memcpy(&h->dime.bitpix, buf2+32, 2);
204 if(!same_order) swabip(buf2+34, 2);
205 memcpy(&h->dime.dim_un0, buf2+34, 2);
206 if(!same_order) swawbip(buf2+36, 32);
207 memcpy(h->dime.pixdim, buf2+36, 32);
208 if(!same_order) swawbip(buf2+68, 4);
209 memcpy(&h->dime.vox_offset, buf2+68, 4);
210 if(!same_order) swawbip(buf2+72, 4);
211 memcpy(&h->dime.funused1, buf2+72, 4);
212 if(!same_order) swawbip(buf2+76, 4);
213 memcpy(&h->dime.funused2, buf2+76, 4);
214 if(!same_order) swawbip(buf2+80, 4);
215 memcpy(&h->dime.funused3, buf2+80, 4);
216 if(!same_order) swawbip(buf2+84, 4);
217 memcpy(&h->dime.cal_max, buf2+84, 4);
218 if(!same_order) swawbip(buf2+88, 4);
219 memcpy(&h->dime.cal_min, buf2+88, 4);
220 if(!same_order) swawbip(buf2+92, 4);
221 memcpy(&h->dime.compressed, buf2+92, 4);
222 if(!same_order) swawbip(buf2+96, 4);
223 memcpy(&h->dime.verified, buf2+96, 4);
224 if(!same_order) swawbip(buf2+100, 4);
225 memcpy(&h->dime.glmax, buf2+100, 4);
226 if(!same_order) swawbip(buf2+104, 4);
227 memcpy(&h->dime.glmin, buf2+104, 4);
228
229 /* Set data history header structure contents */
230 memcpy(h->hist.descrip, buf3+0, 80);
231 memcpy(h->hist.aux_file, buf3+80, 24);
232 memcpy(&h->hist.orient, buf3+104, 1);
233 memcpy(h->hist.originator, buf3+105, 10);
234 memcpy(h->hist.generated, buf3+115, 10);
235 memcpy(h->hist.scannum, buf3+125, 10);
236 memcpy(h->hist.patient_id, buf3+135, 10);
237 memcpy(h->hist.exp_date, buf3+145, 10);
238 memcpy(h->hist.exp_time, buf3+155, 10);
239 memcpy(h->hist.hist_un0, buf3+165, 3);
240 if(!same_order) swawbip(buf3+168, 4);
241 memcpy(&h->hist.views, buf3+168, 4);
242 if(!same_order) swawbip(buf3+172, 4);
243 memcpy(&h->hist.vols_added, buf3+172, 4);
244 if(!same_order) swawbip(buf3+176, 4);
245 memcpy(&h->hist.start_field, buf3+176,4);
246 if(!same_order) swawbip(buf3+180, 4);
247 memcpy(&h->hist.field_skip, buf3+180, 4);
248 if(!same_order) swawbip(buf3+184, 4);
249 memcpy(&h->hist.omax, buf3+184, 4);
250 if(!same_order) swawbip(buf3+188, 4);
251 memcpy(&h->hist.omin, buf3+188, 4);
252 if(!same_order) swawbip(buf3+192, 4);
253 memcpy(&h->hist.smax, buf3+192, 4);
254 if(!same_order) swawbip(buf3+196, 4);
255 memcpy(&h->hist.smin, buf3+196, 4);
256
257 /* Check header contents */
258 if(h->hk.extents!=16384 && h->hk.extents!=0) {
259 if(ANALYZE_TEST>1) printf("hk.extents := %d\n", h->hk.extents);
260 return(11);
261 }
262 if(h->hk.regular!='r') {
263 if(ANALYZE_TEST>1) {
264 if(isalnum(h->hk.regular)) printf("hk.regular := %c\n", h->hk.regular);
265 else printf("hk.regular := (%d)\n", h->hk.regular);
266 }
267 return(12);
268 }
269
270 return(0);
271}
272/*****************************************************************************/
273
274/*****************************************************************************/
283int anaWriteHeader(
285 char *filename,
287 ANALYZE_DSR *h
288) {
289 unsigned char buf1[ANALYZE_HEADER_KEY_SIZE];
290 unsigned char buf2[ANALYZE_HEADER_IMGDIM_SIZE];
291 unsigned char buf3[ANALYZE_HEADER_HISTORY_SIZE];
292 FILE *fp;
293 int same_order, little;
294
295
296 if(ANALYZE_TEST) printf("anaWriteHeader(%s, *dsr)\n", filename);
297
298 /* Check arguments */
299 if(strlen(filename)<1 || h==NULL) return(1);
300 little=little_endian();
301 if(little==h->little) same_order=1; else same_order=0;
302
303 /* Copy header contents into buffers */
304 /* Header key */
305 memset(buf1, 0, ANALYZE_HEADER_KEY_SIZE);
306 memcpy(buf1+0, &h->hk.sizeof_hdr, 4); if(!same_order) swawbip(buf1+0, 4);
307 memcpy(buf1+4, &h->hk.data_type, 10);
308 memcpy(buf1+14, &h->hk.db_name, 18);
309 memcpy(buf1+32, &h->hk.extents, 4); if(!same_order) swawbip(buf1+32, 4);
310 memcpy(buf1+36, &h->hk.session_error, 2); if(!same_order) swabip(buf1+36, 2);
311 memcpy(buf1+38, &h->hk.regular, 1);
312 memcpy(buf1+39, &h->hk.hkey_un0, 1);
313 /* Image dimension */
314 memset(buf2, 0, ANALYZE_HEADER_IMGDIM_SIZE);
315 memcpy(buf2+0, h->dime.dim, 16); if(!same_order) swabip(buf2+0, 16);
316 memcpy(buf2+16, &h->dime.unused8, 2); if(!same_order) swabip(buf2+16, 2);
317 memcpy(buf2+18, &h->dime.unused9, 2); if(!same_order) swabip(buf2+18, 2);
318 memcpy(buf2+20, &h->dime.unused10, 2); if(!same_order) swabip(buf2+20, 2);
319 memcpy(buf2+22, &h->dime.unused11, 2); if(!same_order) swabip(buf2+22, 2);
320 memcpy(buf2+24, &h->dime.unused12, 2); if(!same_order) swabip(buf2+24, 2);
321 memcpy(buf2+26, &h->dime.unused13, 2); if(!same_order) swabip(buf2+26, 2);
322 memcpy(buf2+28, &h->dime.unused14, 2); if(!same_order) swabip(buf2+28, 2);
323 memcpy(buf2+30, &h->dime.datatype, 2); if(!same_order) swabip(buf2+30, 2);
324 memcpy(buf2+32, &h->dime.bitpix, 2); if(!same_order) swabip(buf2+32, 2);
325 memcpy(buf2+34, &h->dime.dim_un0, 2); if(!same_order) swabip(buf2+34, 2);
326 memcpy(buf2+36, h->dime.pixdim, 32); if(!same_order) swawbip(buf2+36, 32);
327 memcpy(buf2+68, &h->dime.vox_offset, 4); if(!same_order) swawbip(buf2+68, 4);
328 memcpy(buf2+72, &h->dime.funused1, 4); if(!same_order) swawbip(buf2+72, 4);
329 memcpy(buf2+76, &h->dime.funused2, 4); if(!same_order) swawbip(buf2+76, 4);
330 memcpy(buf2+80, &h->dime.funused3, 4); if(!same_order) swawbip(buf2+80, 4);
331 memcpy(buf2+84, &h->dime.cal_max, 4); if(!same_order) swawbip(buf2+84, 4);
332 memcpy(buf2+88, &h->dime.cal_min, 4); if(!same_order) swawbip(buf2+88, 4);
333 memcpy(buf2+92, &h->dime.compressed, 4); if(!same_order) swawbip(buf2+92, 4);
334 memcpy(buf2+96, &h->dime.verified, 4); if(!same_order) swawbip(buf2+96, 4);
335 memcpy(buf2+100, &h->dime.glmax, 4); if(!same_order) swawbip(buf2+100, 4);
336 memcpy(buf2+104, &h->dime.glmin, 4); if(!same_order) swawbip(buf2+104, 4);
337 /* Data history */
338 memset(buf3, 0, ANALYZE_HEADER_HISTORY_SIZE);
339 memcpy(buf3+0, &h->hist.descrip, 80);
340 memcpy(buf3+80, &h->hist.aux_file, 24);
341 memcpy(buf3+104, &h->hist.orient, 1);
342 memcpy(buf3+105, &h->hist.originator, 10);
343 memcpy(buf3+115, &h->hist.generated, 10);
344 memcpy(buf3+125, &h->hist.scannum, 10);
345 memcpy(buf3+135, &h->hist.patient_id, 10);
346 memcpy(buf3+145, &h->hist.exp_date, 10);
347 memcpy(buf3+155, &h->hist.exp_time, 10);
348 memcpy(buf3+165, &h->hist.hist_un0, 3);
349 memcpy(buf3+168, &h->hist.views, 4); if(!same_order) swawbip(buf3+168, 4);
350 memcpy(buf3+172, &h->hist.vols_added,4); if(!same_order) swawbip(buf3+172, 4);
351 memcpy(buf3+176, &h->hist.start_field,4); if(!same_order) swawbip(buf3+176,4);
352 memcpy(buf3+180, &h->hist.field_skip,4); if(!same_order) swawbip(buf3+180, 4);
353 memcpy(buf3+184, &h->hist.omax, 4); if(!same_order) swawbip(buf3+184, 4);
354 memcpy(buf3+188, &h->hist.omin, 4); if(!same_order) swawbip(buf3+188, 4);
355 memcpy(buf3+192, &h->hist.smax, 4); if(!same_order) swawbip(buf3+192, 4);
356 memcpy(buf3+196, &h->hist.smin, 4); if(!same_order) swawbip(buf3+196, 4);
357
358 /* Open header file for write */
359 fp=fopen(filename, "wb"); if(fp==NULL) return(2);
360 /* Write header key */
361 if(fwrite(buf1, 1, ANALYZE_HEADER_KEY_SIZE, fp) != ANALYZE_HEADER_KEY_SIZE) {
362 fclose(fp); return(3);}
363 /* Write image dimension */
364 if(fwrite(buf2, 1, ANALYZE_HEADER_IMGDIM_SIZE, fp) != ANALYZE_HEADER_IMGDIM_SIZE) {
365 fclose(fp); return(4);}
366 /* Write data history */
367 if(fwrite(buf3, 1, ANALYZE_HEADER_HISTORY_SIZE, fp) != ANALYZE_HEADER_HISTORY_SIZE) {
368 fclose(fp); return(5);}
369 fclose(fp);
370
371 return(0);
372}
373/*****************************************************************************/
374
375/*****************************************************************************/
380int anaPrintHeader(
382 ANALYZE_DSR *h,
384 FILE *fp
385) {
386 if(fp==NULL || h==NULL) return(1);
387 fprintf(fp, "original_byte_order := %d (1=little, 0=big)\n", h->little);
388 /* Key */
389 fprintf(fp, "header_key.sizeof_hdr := %d\n", h->hk.sizeof_hdr);
390 fprintf(fp, "header_key.data_type := %.10s\n", h->hk.data_type);
391 fprintf(fp, "header_key.db_name := %.18s\n", h->hk.db_name);
392 fprintf(fp, "header_key.extents := %d\n", h->hk.extents);
393 fprintf(fp, "header_key.session_error := %d\n", h->hk.session_error);
394 fprintf(fp, "header_key.regular := %d (%c)\n",
395 (int)h->hk.regular, h->hk.regular);
396 fprintf(fp, "header_key.hkey_un0 := %d\n", (int)h->hk.hkey_un0);
397 /* Image dimension */
398 fprintf(fp, "header_image_dimension.dim :=");
399 for(int i=0; i<8; i++) fprintf(fp, " %d", h->dime.dim[i]);
400 fprintf(fp, "\n");
401 fprintf(fp, "header_image_dimension.unused8 := %d\n", h->dime.unused8);
402 fprintf(fp, "header_image_dimension.unused9 := %d\n", h->dime.unused9);
403 fprintf(fp, "header_image_dimension.unused10 := %d\n", h->dime.unused10);
404 fprintf(fp, "header_image_dimension.unused11 := %d\n", h->dime.unused11);
405 fprintf(fp, "header_image_dimension.unused12 := %d\n", h->dime.unused12);
406 fprintf(fp, "header_image_dimension.unused13 := %d\n", h->dime.unused13);
407 fprintf(fp, "header_image_dimension.unused14 := %d\n", h->dime.unused14);
408 fprintf(fp, "header_image_dimension.datatype := %d\n", h->dime.datatype);
409 fprintf(fp, "header_image_dimension.bitpix := %d\n", h->dime.bitpix);
410 fprintf(fp, "header_image_dimension.dim_un0 := %d\n", h->dime.dim_un0);
411 fprintf(fp, "header_image_dimension.pixdim :=");
412 for(int i=0; i<8; i++) fprintf(fp, " %g", h->dime.pixdim[i]);
413 fprintf(fp, "\n");
414 fprintf(fp, "header_image_dimension.vox_offset := %g\n", h->dime.vox_offset);
415 fprintf(fp, "header_image_dimension.funused1 := %g\n", h->dime.funused1);
416 fprintf(fp, "header_image_dimension.funused2 := %g\n", h->dime.funused2);
417 fprintf(fp, "header_image_dimension.funused3 := %g\n", h->dime.funused3);
418 fprintf(fp, "header_image_dimension.cal_max := %g\n", h->dime.cal_max);
419 fprintf(fp, "header_image_dimension.cal_min := %g\n", h->dime.cal_min);
420 fprintf(fp, "header_image_dimension.compressed := %g\n", h->dime.compressed);
421 fprintf(fp, "header_image_dimension.verified := %g\n", h->dime.verified);
422 fprintf(fp, "header_image_dimension.glmax := %d\n", h->dime.glmax);
423 fprintf(fp, "header_image_dimension.glmin := %d\n", h->dime.glmin);
424 /* Data history */
425 fprintf(fp, "header_data_history.descrip := %s.80\n", h->hist.descrip);
426 fprintf(fp, "header_data_history.aux_file := %.24s\n", h->hist.aux_file);
427 fprintf(fp, "header_data_history.orient := %d\n", (int)h->hist.orient);
428 fprintf(fp, "header_data_history.originator := %.10s\n", h->hist.originator);
429 fprintf(fp, "header_data_history.generated := %.10s\n", h->hist.generated);
430 fprintf(fp, "header_data_history.scannum := %.10s\n", h->hist.scannum);
431 fprintf(fp, "header_data_history.patient_id := %.10s\n", h->hist.patient_id);
432 fprintf(fp, "header_data_history.exp_date := %.10s\n", h->hist.exp_date);
433 fprintf(fp, "header_data_history.exp_time := %.10s\n", h->hist.exp_time);
434 fprintf(fp, "header_data_history.hist_un0 := %.3s\n", h->hist.hist_un0);
435 fprintf(fp, "header_data_history.views := %d\n", h->hist.views);
436 fprintf(fp, "header_data_history.vols_added := %d\n", h->hist.vols_added);
437 fprintf(fp, "header_data_history.start_field := %d\n", h->hist.start_field);
438 fprintf(fp, "header_data_history.field_skip := %d\n", h->hist.field_skip);
439 fprintf(fp, "header_data_history.omax := %d\n", h->hist.omax);
440 fprintf(fp, "header_data_history.omin := %d\n", h->hist.omin);
441 fprintf(fp, "header_data_history.smax := %d\n", h->hist.smax);
442 fprintf(fp, "header_data_history.smin := %d\n", h->hist.smin);
443
444 return(0);
445}
446/*****************************************************************************/
447
448/*****************************************************************************/
454int anaReadImagedata(
456 FILE *fp,
458 ANALYZE_DSR *h,
460 int frame,
462 float *data
463) {
464 int little;
465 long int dimNr, dimx, dimy, dimz=1, dimt=1;
466 char *mdata, *mptr;
467 float f, *fptr;
468 short int *sptr;
469 int *iptr;
470 double d;
471
472
473 if(ANALYZE_TEST) printf("anaReadImagedata(fp, h, %d, data)\n", frame);
474
475 /* Check the arguments */
476 if(frame<=0 || fp==NULL || h==NULL || data==NULL) return(1);
477
478 /* Get the image dimensions from header */
479 dimNr=h->dime.dim[0]; if(dimNr<2) return(2);
480 dimx=h->dime.dim[1];
481 dimy=h->dime.dim[2];
482 if(dimNr>2) dimz=h->dime.dim[3];
483 if(dimNr>3) dimt=h->dime.dim[4];
484 if(frame>dimt) return(3);
485 long long pxlNr=dimx*dimy*dimz; if(pxlNr<1) return(4);
486
487 /* Allocate memory for the binary data */
488 if(h->dime.bitpix==0 && h->dime.datatype>0) { /* Fix bitpix if necessary */
489 if(h->dime.datatype==ANALYZE_DT_BINARY) h->dime.bitpix=1;
490 else if(h->dime.datatype==ANALYZE_DT_UNSIGNED_CHAR) h->dime.bitpix=8;
491 else if(h->dime.datatype==ANALYZE_DT_SIGNED_SHORT) h->dime.bitpix=16;
492 else if(h->dime.datatype==ANALYZE_DT_SIGNED_INT) h->dime.bitpix=32;
493 else if(h->dime.datatype==ANALYZE_DT_FLOAT) h->dime.bitpix=32;
494 else if(h->dime.datatype==ANALYZE_DT_DOUBLE) h->dime.bitpix=64;
495 else if(h->dime.datatype==ANALYZE_DT_RGB) h->dime.bitpix=24;
496 }
497 if(h->dime.bitpix<8) return(5); /* We don't support bit data */
498 long long rawSize=pxlNr*(h->dime.bitpix/8); if(rawSize<1) return(5);
499 if(ANALYZE_TEST>0) printf(" pxlNr=%lld rawSize=%lld\n", pxlNr, rawSize);
500 mdata=(char*)malloc(rawSize); if(mdata==NULL) return(11);
501
502 /* Seek the start of current frame data */
503 long long start_pos=(frame-1)*rawSize;
504 {
505 long int n=(long int)h->dime.vox_offset;
506 if((n>0 && frame==1) || (n<0)) start_pos+=labs(n);
507 }
508 if(ANALYZE_TEST>2) printf("start_pos=%lld\n", start_pos);
509 fseeko(fp, start_pos, SEEK_SET);
510 if(ftello(fp)!=start_pos) {
511 if(ANALYZE_TEST>5) printf("could not move to start_pos\n");
512 free(mdata); return(7);
513 }
514
515 /* Read the data */
516 mptr=mdata;
517 {
518 size_t n=fread(mptr, rawSize, 1, fp);
519 if(n<1) {
520 if(ANALYZE_TEST>5) printf("could read only %zu bytes when request was %lld\n", n, rawSize);
521 free(mdata); return(8);
522 }
523 }
524
525 /* Convert byte order if necessary */
526 little=little_endian(); mptr=mdata;
527 if(little!=h->little) {
528 if(ANALYZE_TEST>0) printf("byte conversion\n");
529 switch(h->dime.bitpix) {
530 case 8: /* no conversion needed */ break;
531 case 16: swabip(mptr, rawSize); break;
532 case 32: swawbip(mptr, rawSize); break;
533 case 64: swawbip(mptr, rawSize); break;
534 default:
535 if(ANALYZE_TEST>5) printf("unsupported anahdr.dime.bitpix := %d\n", h->dime.bitpix);
536 free(mdata); return(5);
537 }
538 }
539
540 /* Get scale factor */
541 f=1.0;
542 if(h->dime.funused1>0.0) f*=h->dime.funused1;
543
544 /* Copy data to float pixel values */
545 mptr=mdata; fptr=data;
546 switch(h->dime.datatype) {
547 case ANALYZE_DT_UNSIGNED_CHAR:
548 if(h->dime.bitpix!=8) {
549 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
550 h->dime.datatype, h->dime.bitpix);
551 free(mdata); return(5);
552 }
553 for(long long i=0; i<pxlNr; i++, mptr++, fptr++)
554 *fptr=f*(float)(unsigned char)(*mptr);
555 break;
556 case ANALYZE_DT_SIGNED_SHORT:
557 if(h->dime.bitpix!=16) {
558 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
559 h->dime.datatype, h->dime.bitpix);
560 free(mdata); return(5);
561 }
562 for(long long i=0; i<pxlNr; i++, mptr+=2, fptr++) {
563 sptr=(short int*)mptr; *fptr=f*(float)(*sptr);
564 }
565 break;
566 case ANALYZE_DT_SIGNED_INT:
567 if(h->dime.bitpix!=16 && h->dime.bitpix!=32) {
568 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
569 h->dime.datatype, h->dime.bitpix);
570 free(mdata); return(5);
571 }
572 if(h->dime.bitpix==16) {
573 for(long long i=0; i<pxlNr; i++, mptr+=2, fptr++) {
574 iptr=(int*)mptr; *fptr=f*(float)(*iptr);
575 }
576 } else if(h->dime.bitpix==32) {
577 for(long long i=0; i<pxlNr; i++, mptr+=4, fptr++) {
578 iptr=(int*)mptr; *fptr=f*(float)(*iptr);
579 }
580 }
581 break;
582 case ANALYZE_DT_FLOAT:
583 if(h->dime.bitpix!=16 && h->dime.bitpix!=32) {
584 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
585 h->dime.datatype, h->dime.bitpix);
586 free(mdata); return(5);
587 }
588 if(h->dime.bitpix==16) {
589 memcpy(fptr, mptr, pxlNr*2);
590 for(long long i=0; i<pxlNr; i++, fptr++) *fptr*=f;
591 } else if(h->dime.bitpix==32) {
592 memcpy(fptr, mptr, pxlNr*4);
593 for(long long i=0; i<pxlNr; i++, fptr++) *fptr*=f;
594 }
595 break;
596 case ANALYZE_DT_COMPLEX:
597 if(h->dime.bitpix!=32) {
598 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
599 h->dime.datatype, h->dime.bitpix);
600 free(mdata); return(5);
601 }
602 if(h->dime.bitpix==32) {
603 memcpy(fptr, mptr, pxlNr*4);
604 for(long long i=0; i<pxlNr; i++, fptr++) *fptr*=f;
605 }
606 break;
607 case ANALYZE_DT_DOUBLE:
608 /* Add support for 64-bit double, if needed */
609 if(h->dime.bitpix!=32) {
610 if(ANALYZE_TEST>5) printf("invalid combination of datatype and bitpix (%d, %d)\n",
611 h->dime.datatype, h->dime.bitpix);
612 free(mdata); return(5);
613 }
614 for(long long i=0; i<pxlNr; i++, mptr+=8, fptr++) {
615 memcpy(&d, mptr, 8); *fptr=f*d;
616 }
617 break;
618 default:
619 if(ANALYZE_TEST>5) printf("unsupported anahdr.dime.datatype := %d\n", h->dime.datatype);
620 free(mdata); return(5);
621 }
622
623 free(mdata);
624 if(ANALYZE_TEST>1) printf("anaReadImagedata() succeeded\n");
625 return(0);
626}
627/*****************************************************************************/
628
629/*****************************************************************************/
635int anaFlipping() {
636 int ret;
637 char *cptr;
638
639 /* Is there an environment variable name for flipping? */
640 cptr=getenv("ANALYZE_FLIP");
641 if(cptr==NULL) cptr=getenv("ANALYZE_FLIPPING");
642 if(cptr==NULL) cptr=getenv("analyze_flip");
643 if(cptr==NULL) cptr=getenv("analyze_flipping");
644 if(cptr==NULL) {
645 if(ANALYZE_TEST>1) printf("ANALYZE_FLIP = not defined\n");
646 ret=ANALYZE_FLIP_DEFAULT; /* if not, then use default value */
647 } else {
648 if(ANALYZE_TEST>1) printf("ANALYZE_FLIP = '%s'\n", cptr);
649 if(*cptr=='y' || *cptr=='Y' || *cptr=='1') ret=1;
650 else if(*cptr=='n' || *cptr=='N' || *cptr=='0') ret=0;
651 else ret=ANALYZE_FLIP_DEFAULT;
652 }
653 if(ANALYZE_TEST) printf("anaFlipping()=%d\n", ret);
654 return(ret);
655}
656/*****************************************************************************/
657
658/*****************************************************************************/
666int anaRemove(
668 const char *dbname
669) {
670 char datfile[FILENAME_MAX], hdrfile[FILENAME_MAX], siffile[FILENAME_MAX];
671
672 if(ANALYZE_TEST) printf("anaRemove(%s)\n", dbname);
673 if(anaDatabaseExists(dbname, hdrfile, datfile, siffile)==0) return 0;
674 if(ANALYZE_TEST>2) printf(" removing %s and %s\n", hdrfile, datfile);
675 if(remove(hdrfile)!=0) return 1;
676 if(remove(datfile)!=0) return 2;
677 return 0;
678}
679/*****************************************************************************/
680
681/*****************************************************************************/
687void anaRemoveFNameExtension(
689 char *fname
690) {
691 char *cptr;
692 cptr=strrchr(fname, '.'); if(cptr==NULL) return;
693 if(strcasecmp(cptr, ".")==0 || strcasecmp(cptr, ".img")==0 ||
694 strcasecmp(cptr, ".hdr")==0 || strcasecmp(cptr, ".sif")==0)
695 *cptr=(char)0;
696}
697/*****************************************************************************/
698
699/*****************************************************************************/
704int anaDatabaseExists(
707 const char *dbname,
710 char *hdrfile,
713 char *imgfile,
716 char *siffile
717) {
718 char temp[FILENAME_MAX], database[FILENAME_MAX];
719 int checked=0;
720
721 if(ANALYZE_TEST)
722 printf("\nanaDatabaseExists(%s, *hdrfile, *imgfile, *siffile)\n", dbname);
723
724 /* Check the input */
725 if(hdrfile!=NULL) strcpy(hdrfile, "");
726 if(imgfile!=NULL) strcpy(imgfile, "");
727 if(siffile!=NULL) strcpy(siffile, "");
728 if(dbname==NULL || strlen(dbname)==0) return(0);
729
730 strcpy(database, dbname);
731 while(1) {
732 /* Header file? */
733 strcpy(temp, database); strcat(temp, ".hdr");
734 if(access(temp, 0) != -1) {
735 /* Also image file? */
736 strcpy(temp, database); strcat(temp, ".img");
737 if(access(temp, 0) != -1) {
738 if(hdrfile!=NULL) sprintf(hdrfile, "%s.hdr", database);
739 if(imgfile!=NULL) sprintf(imgfile, "%s.img", database);
740 /* Even SIF? */
741 if(anaMakeSIFName(database, temp)==0) { /* yes! */
742 if(siffile!=NULL) strcpy(siffile, temp);
743 return(2);
744 }
745 /* Image and header files did exist anyway */
746 return(1);
747 }
748 }
749 if(checked==1) break;
750 /* Try to remove extension */
751 anaRemoveFNameExtension(database);
752 checked=1;
753 } /* try again */
754 return(0);
755}
756/*****************************************************************************/
757
758/*****************************************************************************/
763int anaMakeSIFName(
765 const char *dbname,
767 char *siffile
768) {
769 if(dbname==NULL || siffile==NULL) return(1);
770 sprintf(siffile, "%s.sif", dbname); if(access(siffile, 0) != -1) return(0);
771 sprintf(siffile, "%s.SIF", dbname); if(access(siffile, 0) != -1) return(0);
772 sprintf(siffile, "%s.img.sif", dbname); if(access(siffile, 0) != -1) return(0);
773 sprintf(siffile, "%s.IMG.SIF", dbname); if(access(siffile, 0) != -1) return(0);
774 sprintf(siffile, "%s.sif", dbname); return(2);
775}
776/*****************************************************************************/
777
778/*****************************************************************************/
783int anaEditHeader(
785 ANALYZE_DSR *h,
787 char *field,
789 char *value
790) {
791 int ii;
792 short int si;
793 float f;
794
795 si=atoi(value); ii=atoi(value); f=atof(value);
796 /* Header keys */
797 if(strcmp(field, "header_key.sizeof_hdr")==0 || strcmp(field, "sizeof_hdr")==0) {
798 h->hk.sizeof_hdr=ii;
799 } else if(strcmp(field, "header_key.data_type")==0 || strcmp(field, "data_type")==0) {
800 strlcpy(h->hk.data_type, value, 10);
801 } else if(strcmp(field, "header_key.db_name")==0 || strcmp(field, "db_name")==0) {
802 strlcpy(h->hk.db_name, value, 18);
803 } else if(strcmp(field, "header_key.extents")==0 || strcmp(field, "extents")==0) {
804 h->hk.extents=ii;
805 } else if(strcmp(field, "header_key.session_error")==0 || strcmp(field, "session_error")==0) {
806 h->hk.session_error=si;
807 } else if(strcmp(field, "header_key.regular")==0 || strcmp(field, "regular")==0) {
808 h->hk.regular=value[0];
809 } else if(strcmp(field, "header_key.hkey_un0")==0 || strcmp(field, "hkey_un0")==0) {
810 h->hk.hkey_un0=value[0];
811 /* Header imgdim */
812 } else if(strcmp(field, "header_image_dimension.dim")==0 || strcmp(field, "dim")==0) {
813 sscanf(value, "%hd %hd %hd %hd %hd %hd %hd %hd",
814 h->dime.dim+0, h->dime.dim+1, h->dime.dim+2, h->dime.dim+3,
815 h->dime.dim+4, h->dime.dim+5, h->dime.dim+6, h->dime.dim+7
816 );
817 } else if(strcmp(field, "header_image_dimension.unused8")==0 || strcmp(field, "unused8")==0) {
818 h->dime.unused8=si;
819 } else if(strcmp(field, "header_image_dimension.unused9")==0 || strcmp(field, "unused9")==0) {
820 h->dime.unused9=si;
821 } else if(strcmp(field, "header_image_dimension.unused10")==0 || strcmp(field, "unused10")==0) {
822 h->dime.unused10=si;
823 } else if(strcmp(field, "header_image_dimension.unused11")==0 || strcmp(field, "unused11")==0) {
824 h->dime.unused11=si;
825 } else if(strcmp(field, "header_image_dimension.unused12")==0 || strcmp(field, "unused12")==0) {
826 h->dime.unused12=si;
827 } else if(strcmp(field, "header_image_dimension.unused13")==0 || strcmp(field, "unused13")==0) {
828 h->dime.unused13=si;
829 } else if(strcmp(field, "header_image_dimension.unused14")==0 || strcmp(field, "unused14")==0) {
830 h->dime.unused14=si;
831 } else if(strcmp(field, "header_image_dimension.datatype")==0 || strcmp(field, "datatype")==0) {
832 h->dime.datatype=si;
833 } else if(strcmp(field, "header_image_dimension.bitpix")==0 || strcmp(field, "bitpix")==0) {
834 h->dime.bitpix=si;
835 } else if(strcmp(field, "header_image_dimension.dim_un0")==0 || strcmp(field, "dim_un0")==0) {
836 h->dime.dim_un0=si;
837 } else if(strcmp(field, "header_image_dimension.pixdim")==0 || strcmp(field, "pixdim")==0) {
838 sscanf(value, "%f %f %f %f %f %f %f %f",
839 h->dime.pixdim+0, h->dime.pixdim+1, h->dime.pixdim+2, h->dime.pixdim+3,
840 h->dime.pixdim+4, h->dime.pixdim+5, h->dime.pixdim+6, h->dime.pixdim+7
841 );
842 } else if(strcmp(field, "header_image_dimension.vox_offset")==0 || strcmp(field, "vox_offset")==0) {
843 h->dime.vox_offset=f;
844 } else if(strcmp(field, "header_image_dimension.funused1")==0 || strcmp(field, "funused1")==0) {
845 h->dime.funused1=f;
846 } else if(strcmp(field, "header_image_dimension.funused2")==0 || strcmp(field, "funused2")==0) {
847 h->dime.funused2=f;
848 } else if(strcmp(field, "header_image_dimension.funused3")==0 || strcmp(field, "funused3")==0) {
849 h->dime.funused3=f;
850 } else if(strcmp(field, "header_image_dimension.cal_max")==0 || strcmp(field, "cal_max")==0) {
851 h->dime.cal_max=f;
852 } else if(strcmp(field, "header_image_dimension.cal_min")==0 || strcmp(field, "cal_min")==0) {
853 h->dime.cal_min=f;
854 } else if(strcmp(field, "header_image_dimension.compressed")==0 || strcmp(field, "compressed")==0) {
855 h->dime.compressed=f;
856 } else if(strcmp(field, "header_image_dimension.verified")==0 || strcmp(field, "verified")==0) {
857 h->dime.verified=f;
858 } else if(strcmp(field, "header_image_dimension.glmax")==0 || strcmp(field, "glmax")==0) {
859 h->dime.glmax=ii;
860 } else if(strcmp(field, "header_image_dimension.glmin")==0 || strcmp(field, "glmin")==0) {
861 h->dime.glmin=ii;
862 /* Header history */
863 } else if(strcmp(field, "header_data_history.descrip")==0 || strcmp(field, "descrip")==0) {
864 strlcpy(h->hist.descrip, value, 80);
865 } else if(strcmp(field, "header_data_history.aux_file")==0 || strcmp(field, "aux_file")==0) {
866 strlcpy(h->hist.aux_file, value, 24);
867 } else if(strcmp(field, "header_data_history.orient")==0 || strcmp(field, "orient")==0) {
868 h->hist.orient=value[0];
869 } else if(strcmp(field, "header_data_history.originator")==0 || strcmp(field, "originator")==0) {
870 strlcpy(h->hist.originator, value, 10);
871 } else if(strcmp(field, "header_data_history.generated")==0 || strcmp(field, "generated")==0) {
872 strlcpy(h->hist.generated, value, 10);
873 } else if(strcmp(field, "header_data_history.scannum")==0 || strcmp(field, "scannum")==0) {
874 strlcpy(h->hist.scannum, value, 10);
875 } else if(strcmp(field, "header_data_history.patient_id")==0 || strcmp(field, "patient_id")==0) {
876 strlcpy(h->hist.patient_id, value, 10);
877 } else if(strcmp(field, "header_data_history.exp_date")==0 || strcmp(field, "exp_date")==0) {
878 strlcpy(h->hist.exp_date, value, 10);
879 } else if(strcmp(field, "header_data_history.exp_time")==0 || strcmp(field, "exp_time")==0) {
880 strlcpy(h->hist.exp_time, value, 10);
881 } else if(strcmp(field, "header_data_history.hist_un0")==0 || strcmp(field, "hist_un0")==0) {
882 memcpy(h->hist.hist_un0, value, 3);
883 } else if(strcmp(field, "header_data_history.views")==0 || strcmp(field, "views")==0) {
884 h->hist.views=ii;
885 } else if(strcmp(field, "header_data_history.vols_added")==0 || strcmp(field, "vols_added")==0) {
886 h->hist.vols_added=ii;
887 } else if(strcmp(field, "header_data_history.start_field")==0 || strcmp(field, "start_field")==0) {
888 h->hist.start_field=ii;
889 } else if(strcmp(field, "header_data_history.field_skip")==0 || strcmp(field, "field_skip")==0) {
890 h->hist.field_skip=ii;
891 } else if(strcmp(field, "header_data_history.omax")==0 || strcmp(field, "omax")==0) {
892 h->hist.omax=ii;
893 } else if(strcmp(field, "header_data_history.omin")==0 || strcmp(field, "omin")==0) {
894 h->hist.omin=ii;
895 } else if(strcmp(field, "header_data_history.smax")==0 || strcmp(field, "smax")==0) {
896 h->hist.smax=ii;
897 } else if(strcmp(field, "header_data_history.smin")==0 || strcmp(field, "smin")==0) {
898 h->hist.smin=ii;
899 } else
900 return(1);
901
902 return(0);
903}
904/*****************************************************************************/
905
906/*****************************************************************************/
anaMakeSIFName
int anaMakeSIFName(const char *dbname, char *siffile)
Definition analyze.c:763
anaEditHeader
int anaEditHeader(ANALYZE_DSR *h, char *field, char *value)
Definition analyze.c:783
anaWriteHeader
int anaWriteHeader(char *filename, ANALYZE_DSR *h)
Definition analyze.c:283
anaExistsNew
int anaExistsNew(const char *filename, char *hdrfile, char *imgfile, char *siffile)
Definition analyze.c:43
anaFlipping
int anaFlipping()
Definition analyze.c:635
anaPrintHeader
int anaPrintHeader(ANALYZE_DSR *h, FILE *fp)
Definition analyze.c:380
anaDatabaseExists
int anaDatabaseExists(const char *dbname, char *hdrfile, char *imgfile, char *siffile)
Definition analyze.c:704
anaRemoveFNameExtension
void anaRemoveFNameExtension(char *fname)
Definition analyze.c:687
ANALYZE_TEST
int ANALYZE_TEST
Definition analyze.c:8
anaExists
int anaExists(const char *dbname)
Definition analyze.c:17
anaReadHeader
int anaReadHeader(char *filename, ANALYZE_DSR *h)
Definition analyze.c:129
anaReadImagedata
int anaReadImagedata(FILE *fp, ANALYZE_DSR *h, int frame, float *data)
Definition analyze.c:454
anaRemove
int anaRemove(const char *dbname)
Definition analyze.c:666
libtpcimgio.h
Header file for libtpcimgio.
ANALYZE_DT_SIGNED_INT
#define ANALYZE_DT_SIGNED_INT
Definition libtpcimgio.h:57
ANALYZE_DT_COMPLEX
#define ANALYZE_DT_COMPLEX
Definition libtpcimgio.h:61
ANALYZE_DT_SIGNED_SHORT
#define ANALYZE_DT_SIGNED_SHORT
Definition libtpcimgio.h:55
ANALYZE_DT_BINARY
#define ANALYZE_DT_BINARY
Definition libtpcimgio.h:51
ANALYZE_DT_RGB
#define ANALYZE_DT_RGB
Definition libtpcimgio.h:65
ANALYZE_DT_UNSIGNED_CHAR
#define ANALYZE_DT_UNSIGNED_CHAR
Definition libtpcimgio.h:53
ANALYZE_DT_DOUBLE
#define ANALYZE_DT_DOUBLE
Definition libtpcimgio.h:63
ANALYZE_DT_FLOAT
#define ANALYZE_DT_FLOAT
Definition libtpcimgio.h:59
swabip
void swabip(void *buf, long long int size)
Definition swap.c:72
swawbip
void swawbip(void *buf, long long int size)
Definition swap.c:93
strlcpy
size_t strlcpy(char *dst, const char *src, size_t dstsize)
Definition strext.c:244
little_endian
int little_endian()
Definition swap.c:14
strlcat
size_t strlcat(char *dst, const char *src, size_t dstsize)
Definition strext.c:205
ANALYZE_DSR
Definition libtpcimgio.h:193
ANALYZE_DSR::hist
ANALYZE_HEADER_HISTORY hist
Definition libtpcimgio.h:199
ANALYZE_DSR::little
int little
Definition libtpcimgio.h:201
ANALYZE_DSR::hk
ANALYZE_HEADER_KEY hk
Definition libtpcimgio.h:195
ANALYZE_DSR::dime
ANALYZE_HEADER_IMGDIM dime
Definition libtpcimgio.h:197
ANALYZE_HEADER_HISTORY
Definition libtpcimgio.h:150
ANALYZE_HEADER_HISTORY::scannum
char scannum[10]
Definition libtpcimgio.h:164
ANALYZE_HEADER_HISTORY::exp_time
char exp_time[10]
Definition libtpcimgio.h:170
ANALYZE_HEADER_HISTORY::generated
char generated[10]
Definition libtpcimgio.h:162
ANALYZE_HEADER_HISTORY::aux_file
char aux_file[24]
Definition libtpcimgio.h:154
ANALYZE_HEADER_HISTORY::exp_date
char exp_date[10]
Definition libtpcimgio.h:168
ANALYZE_HEADER_HISTORY::hist_un0
char hist_un0[3]
Definition libtpcimgio.h:172
ANALYZE_HEADER_HISTORY::vols_added
int vols_added
Definition libtpcimgio.h:176
ANALYZE_HEADER_HISTORY::omax
int omax
Definition libtpcimgio.h:182
ANALYZE_HEADER_HISTORY::smax
int smax
Definition libtpcimgio.h:186
ANALYZE_HEADER_HISTORY::smin
int smin
Definition libtpcimgio.h:188
ANALYZE_HEADER_HISTORY::field_skip
int field_skip
Definition libtpcimgio.h:180
ANALYZE_HEADER_HISTORY::orient
char orient
Definition libtpcimgio.h:158
ANALYZE_HEADER_HISTORY::omin
int omin
Definition libtpcimgio.h:184
ANALYZE_HEADER_HISTORY::patient_id
char patient_id[10]
Definition libtpcimgio.h:166
ANALYZE_HEADER_HISTORY::views
int views
Definition libtpcimgio.h:174
ANALYZE_HEADER_HISTORY::start_field
int start_field
Definition libtpcimgio.h:178
ANALYZE_HEADER_HISTORY::descrip
char descrip[80]
Definition libtpcimgio.h:152
ANALYZE_HEADER_HISTORY::originator
char originator[10]
Definition libtpcimgio.h:160
ANALYZE_HEADER_IMGDIM::pixdim
float pixdim[8]
Definition libtpcimgio.h:124
ANALYZE_HEADER_IMGDIM::unused14
short int unused14
Definition libtpcimgio.h:110
ANALYZE_HEADER_IMGDIM::dim_un0
short int dim_un0
Definition libtpcimgio.h:122
ANALYZE_HEADER_IMGDIM::funused1
float funused1
Definition libtpcimgio.h:130
ANALYZE_HEADER_IMGDIM::unused9
short int unused9
Definition libtpcimgio.h:100
ANALYZE_HEADER_IMGDIM::datatype
short int datatype
Definition libtpcimgio.h:118
ANALYZE_HEADER_IMGDIM::funused2
float funused2
Definition libtpcimgio.h:132
ANALYZE_HEADER_IMGDIM::funused3
float funused3
Definition libtpcimgio.h:134
ANALYZE_HEADER_IMGDIM::compressed
float compressed
Definition libtpcimgio.h:140
ANALYZE_HEADER_IMGDIM::cal_max
float cal_max
Definition libtpcimgio.h:136
ANALYZE_HEADER_IMGDIM::glmax
int glmax
Definition libtpcimgio.h:144
ANALYZE_HEADER_IMGDIM::cal_min
float cal_min
Definition libtpcimgio.h:138
ANALYZE_HEADER_IMGDIM::verified
float verified
Definition libtpcimgio.h:142
ANALYZE_HEADER_IMGDIM::glmin
int glmin
Definition libtpcimgio.h:146
ANALYZE_HEADER_IMGDIM::unused11
short int unused11
Definition libtpcimgio.h:104
ANALYZE_HEADER_IMGDIM::dim
short int dim[8]
Definition libtpcimgio.h:96
ANALYZE_HEADER_IMGDIM::unused10
short int unused10
Definition libtpcimgio.h:102
ANALYZE_HEADER_IMGDIM::bitpix
short int bitpix
Definition libtpcimgio.h:120
ANALYZE_HEADER_IMGDIM::unused13
short int unused13
Definition libtpcimgio.h:108
ANALYZE_HEADER_IMGDIM::unused8
short int unused8
Definition libtpcimgio.h:98
ANALYZE_HEADER_IMGDIM::unused12
short int unused12
Definition libtpcimgio.h:106
ANALYZE_HEADER_IMGDIM::vox_offset
float vox_offset
Definition libtpcimgio.h:128
ANALYZE_HEADER_KEY::extents
int extents
Definition libtpcimgio.h:82
ANALYZE_HEADER_KEY::regular
char regular
Definition libtpcimgio.h:87
ANALYZE_HEADER_KEY::sizeof_hdr
int sizeof_hdr
Definition libtpcimgio.h:75
ANALYZE_HEADER_KEY::hkey_un0
char hkey_un0
Definition libtpcimgio.h:89
ANALYZE_HEADER_KEY::db_name
char db_name[18]
Definition libtpcimgio.h:79
ANALYZE_HEADER_KEY::data_type
char data_type[10]
Definition libtpcimgio.h:77
ANALYZE_HEADER_KEY::session_error
short int session_error
Definition libtpcimgio.h:84
Generated on for TPCCLIB by doxygen 1.18.0