forked from tkkarjal/magia
-
Notifications
You must be signed in to change notification settings - Fork 0
/
magia_get_freesurfer_roi_info.m
151 lines (149 loc) · 5.34 KB
/
magia_get_freesurfer_roi_info.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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
function roi_info = magia_get_freesurfer_roi_info(specs)
if(strcmp(specs.magia.roi_set,'tracer_default'))
roi_set = magia_get_tracer_default_roi_set(specs.study.tracer);
else
roi_set = specs.magia.roi_set;
end
switch roi_set
case 'rs1'
roi_info.labels = {'amy' 'cau' 'cer' 'dacc' 'inftemp' 'ins' 'midtemp' 'nacc' 'ofc' 'parsop' 'pcc' 'put' 'racc' 'supfront' 'suptemp' 'tempol' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[8 47] % cer
[1002 2002] % dacc
[1009 2009] % inftemp
[1035 2035] % ins
[1015 2015] % midtemp
[26 58] % nacc
[1012 1014 2012 2014] % ofc
[1018 2018] % parsop
[1023 2023] % pcc
[12 102 51 111] % put
[1026 2026] % racc
[1028 2028] % supfront
[1030 2030] % suptemp
[1033 2033] % tempol
[9 10 48 49] % tha
};
case 'rs2'
roi_info.labels = {'amy' 'cau' 'cer' 'dacc' 'hip' 'inftemp' 'ins' 'medul' 'midbr' 'midtemp' 'nacc' 'ofc' 'parsop' 'pcc' 'pons' 'put' 'racc' 'supfront' 'suptemp' 'tempol' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[8 47] % cer
[1002 2002] % dacc
[17 53] % hip
[1009 2009] % inftemp
[1035 2035] % ins
[175] % medul
[173] % midbr
[1015 2015] % midtemp
[26 58] % nacc
[1012 1014 2012 2014] % ofc
[1018 2018] % parsop
[1023 2023] % pcc
[174] % pons
[12 102 51 111] % put
[1026 2026] % racc
[1028 2028] % supfront
[1030 2030] % suptemp
[1033 2033] % tempol
[9 10 48 49] % tha
};
case 'rs3'
roi_info.labels = {'amy' 'cau' 'nacc' 'pal' 'put' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[26 58] % nacc
[13 52] % pal
[12 102 51 111] % put
[9 10 48 49] % tha
};
case 'rs4'
roi_info.labels = {'amy' 'cau' 'dacc' 'hip' 'ins' 'medul' 'midbr' 'nacc' 'oc' 'ofc' 'pal' 'parhip' 'pcc' 'pons' 'put' 'racc' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[1002 2002] % dacc
[17 53] % hip
[1035 2035] % ins
[175] % medul
[173] % midbr
[26 58] % nacc
[1011 2011] % oc
[1012 1014 2012 2014] % ofc
[13 52] % pal
[1016 2016] % parhip
[1023 2023] % pcc
[174] % pons
[12 102 51 111] % put
[1026 2026] % racc
[9 10 48 49] % tha
};
case 'rs5'
roi_info.labels = {'amy' 'cau' 'dacc' 'hip' 'ins' 'nacc' 'ofc' 'parhip' 'pcc' 'pcun' 'put' 'racc' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[1002 2002] % dacc
[1035 2035] % ins
[17 53] % hip
[26 58] % nacc
[1012 1014 2012 2014] % ofc
[1016 2016] % parhip
[1023 2023] % pcc
[1025 2025] % pcun
[12 102 51 111] % put
[1026 2026] % racc
[9 10 48 49] % tha
};
case 'rs6'
roi_info.labels = {'CGA' 'CGP' 'LOC' 'LTC' 'MTC' 'PARCALL' 'PIB1PFCALL' 'PIBcomp' 'PREC' 'STR'};
roi_info.codes = {
[1002 2002 1026 2026] % anterior cingulate
[1010 2010 1023 2023] % posterior cingulate
[1011 2011] % lateral occipital cortex
[1009 2009 1015 2015 1030 2030 1033 2033 1034 2034] % lateral temporal cortex
[17 53 18 54 1006 2006] % medial temporal cortex
[1008 2008 1029 2029 1031 2031] % parietal cortex
[1003 2003 1012 2012 1014 2014 1018 2018 1019 2019 1020 2020 1027 2027 1028 2028 1032 2032] % prefrontal cortex
[1003 2003 1012 2012 1014 2014 1018 2018 1019 2019 1020 2020 1027 2027 1028 2028 1032 2032 1009 2009 1015 2015 1030 2030 1033 2033 1034 2034 1008 2008 1029 2029 1031 2031 1025 2025 1002 2002 1026 2026 1010 2010 1023 2023] % pibcomp
[1025 2025] % precuneus
[11 50 12 51] % striatum
};
case 'rs7'
roi_info.labels = {'amy' 'cau' 'dacc' 'hip' 'ins' 'nacc' 'oc' 'ofc' 'pal' 'parhip' 'pcc' 'put' 'racc' 'tha'};
roi_info.codes = {
[18 54] % amy
[11 50] % cau
[1002 2002] % dacc
[17 53] % hip
[1035 2035] % ins
[26 58] % nacc
[1012 1014 2012 2014] % ofc
[1011 2011] % oc
[13 52] % pal
[1016 2016] % parhip
[1023 2023] % pcc
[12 102 51 111] % put
[1026 2026] % racc
[9 10 48 49] % tha
};
case 'striat_unil'
roi_info.labels = {'cau_l' 'cau_r' 'nacc_l' 'nacc_r' 'pal_l' 'pal_r' 'put_l' 'put_r'};
roi_info.codes = {
11 % cau_l
50 % cau_r
26 % nacc_l
58 % nacc_r
13 % pal_l
52 % pal_r
[12 102] % put_l
[51 111] % put_r
};
otherwise
error('The roi_set %s has not been defined. Please use another roi_set.',roi_set);
end
end