Codebase list octave-iso2mesh / lintian-fixes/main niicodemap.m
lintian-fixes/main

Tree @lintian-fixes/main (Download .tar.gz)

niicodemap.m @lintian-fixes/main — raw · history · blame 

 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/fangq/jnifti
%
%    License: Apache 2.0, see https://github.com/fangq/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