forked from fangq/jsonlab
/
niicodemap.m
93 lines (80 loc) · 3.72 KB
/
niicodemap.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
function newval = niicodemap(name, value)
%
% newval=niicodemap(name, value)
%
% Bi-directional conversion from NIFTI codes to human-readable JNIfTI
% header string values
%
% author: Qianqian Fang (q.fang <at> neu.edu)
%
% input:
% name: a header name as a string, currently support the below nii
% headers: 'intent_code', 'slice_code', 'datatype', 'qform',
% 'sform' and 'xyzt_units' and their corresponding JNIfTI
% headers:
% 'Intent','SliceType','DataType','QForm','SForm','Unit'
% value:the current header value, if it is a code, newval will
% output the string version; if it is a string, newval will
% return the code
%
% output:
% newval: the converted header value
%
% For the detailed nii header codes, please see
% https://nifti.nimh.nih.gov/pub/dist/src/niftilib/nifti1.h
%
% example:
% newval=niicodemap('slice_code', '')
% newval=niicodemap('datatype', 'uint64')
% newval=niicodemap('datatype', 2)
%
% this file was initially developed for the MCX project: https://github.com/fangq/mcx/blob/master/utils/mcxloadnii.m
%
% this file is part of JNIfTI specification: https://github.com/NeuroJSON/jnifti
%
% License: Apache 2.0, see https://github.com/NeuroJSON/jnifti for details
%
% code to name look-up-table
if (~exist('containers.Map'))
newval = value;
return
end
lut.intent_code = containers.Map([0, 2:24 1001:1011 2001:2005], ...
{'', 'corr', 'ttest', 'ftest', 'zscore', 'chi2', 'beta', ...
'binomial', 'gamma', 'poisson', 'normal', 'ncftest', ...
'ncchi2', 'logistic', 'laplace', 'uniform', 'ncttest', ...
'weibull', 'chi', 'invgauss', 'extval', 'pvalue', ...
'logpvalue', 'log10pvalue', 'estimate', 'label', 'neuronames', ...
'matrix', 'symmatrix', 'dispvec', 'vector', 'point', 'triangle', ...
'quaternion', 'unitless', 'tseries', 'elem', 'rgb', 'rgba', 'shape'});
lut.slice_code = containers.Map(0:6, {'', 'seq+', 'seq-', 'alt+', 'alt-', 'alt2+', 'alt-'});
lut.datatype = containers.Map([0, 2, 4, 8, 16, 32, 64, 128, 256, 512, 768, 1024, 1280, 1536, 1792, 2048, 2304], ...
{'', 'uint8', 'int16', 'int32', 'single', 'complex64', 'double', 'rgb24', 'int8', ...
'uint16', 'uint32', 'int64', 'uint64', 'double128', 'complex128', ...
'complex256', 'rgba32' });
lut.xyzt_units = containers.Map([0:3 8 16 24 32 40 48], ...
{'', 'm', 'mm', 'um', 's', 'ms', 'us', 'hz', 'ppm', 'rad'});
lut.qform = containers.Map(0:4, {'', 'scanner', 'aligned', 'talairach', 'mni'});
lut.unit = lut.xyzt_units;
lut.sform = lut.qform;
lut.slicetype = lut.slice_code;
lut.intent = lut.intent_code;
% inverse look up table
tul.intent_code = containers.Map(values(lut.intent_code), keys(lut.intent_code));
tul.slice_code = containers.Map(values(lut.slice_code), keys(lut.slice_code));
tul.datatype = containers.Map(values(lut.datatype), keys(lut.datatype));
tul.xyzt_units = containers.Map(values(lut.xyzt_units), keys(lut.xyzt_units));
tul.qform = containers.Map(values(lut.qform), keys(lut.qform));
tul.sform = tul.qform;
tul.slicetype = tul.slice_code;
tul.intent = tul.intent_code;
tul.unit = tul.xyzt_units;
% map from code to name, or frmo name to code
if (~isfield(lut, lower(name)))
error('property can not be found');
end
if (~(ischar(value) || isa(value, 'string')))
newval = lut.(lower(name))(value);
else
newval = tul.(lower(name))(value);
end