-
Notifications
You must be signed in to change notification settings - Fork 1
/
loaddata.m
137 lines (133 loc) · 5.67 KB
/
loaddata.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
% Load paths, GA (in weeks) and manually determined landmarks (in voxel
% coordinates).
function [data,ga,bman,lman,rman] = loaddata()
d = 'data';
data{ 1} = fullfile(d, 'ep2d_01.nii.gz');
data{ 2} = fullfile(d, 'ep2d_02.nii.gz');
data{ 3} = fullfile(d, 'ep2d_03.nii.gz');
data{ 4} = fullfile(d, 'ep2d_04.nii.gz');
data{ 5} = fullfile(d, 'ep2d_05.nii.gz');
data{ 6} = fullfile(d, 'ep2d_06.nii.gz');
data{ 7} = fullfile(d, 'ep2d_07.nii.gz');
data{ 8} = fullfile(d, 'ep2d_08.nii.gz');
data{ 9} = fullfile(d, 'ep2d_09.nii.gz');
data{10} = fullfile(d, 'ep2d_10.nii.gz');
data{11} = fullfile(d, 'ep2d_11.nii.gz');
data{12} = fullfile(d, 'ep2d_12.nii.gz');
data{13} = fullfile(d, 'ep2d_13.nii.gz');
data{14} = fullfile(d, 'ep2d_14.nii.gz');
data{15} = fullfile(d, 'ep2d_15.nii.gz');
data{16} = fullfile(d, 'ep2d_16.nii.gz');
data{17} = fullfile(d, 'ep2d_17.nii.gz');
data{18} = fullfile(d, 'ep2d_18.nii.gz');
data{19} = fullfile(d, 'ep2d_19.nii.gz');
data{20} = fullfile(d, 'ep2d_20.nii.gz');
data{21} = fullfile(d, 'ep2d_21.nii.gz');
data{22} = fullfile(d, 'ep2d_22.nii.gz');
data{23} = fullfile(d, 'ep2d_23.nii.gz');
data{24} = fullfile(d, 'ep2d_24.nii.gz');
data{25} = fullfile(d, 'ep2d_25.nii.gz');
data{26} = fullfile(d, 'ep2d_26.nii.gz');
data{27} = fullfile(d, 'ep2d_27.nii.gz');
data{28} = fullfile(d, 'ep2d_28.nii.gz');
data{29} = fullfile(d, 'ep2d_29.nii.gz');
data{30} = fullfile(d, 'ep2d_30.nii.gz');
data{31} = fullfile(d, 'ep2d_31.nii.gz');
data{32} = fullfile(d, 'ep2d_32.nii.gz');
data{33} = fullfile(d, 'ep2d_33.nii.gz');
data{34} = fullfile(d, 'ep2d_34.nii.gz');
data{35} = fullfile(d, 'ep2d_35.nii.gz');
data{36} = fullfile(d, 'ep2d_36.nii.gz');
data{37} = fullfile(d, 'ep2d_37.nii.gz');
data{38} = fullfile(d, 'ep2d_38.nii.gz');
data{39} = fullfile(d, 'ep2d_39.nii.gz');
data{40} = fullfile(d, 'ep2d_40.nii.gz');
data{41} = fullfile(d, 'ep2d_41.nii.gz');
% Gestational age.
ga( 1) = 31;
ga( 2) = 31;
ga( 3) = 31;
ga( 4) = 31;
ga( 5) = 31;
ga( 6) = 31;
ga( 7) = 31;
ga( 8) = 31;
ga( 9) = 31;
ga(10) = 26;
ga(11) = 32.57;
ga(12) = 35.14;
ga(13) = 27.14;
ga(14) = 27.14;
ga(15) = 31.57;
ga(16) = 31.57;
ga(17) = 32.86;
ga(18) = 32.86;
ga(19) = 32.86;
ga(20) = 36.71;
ga(21) = 36.71;
ga(22) = 36.71;
ga(23) = 29;
ga(24) = 29;
ga(25) = 28;
ga(26) = 28;
ga(27) = 32;
ga(28) = 32;
ga(29) = 28;
ga(30) = 28;
ga(31) = 31;
ga(32) = 31;
ga(33) = 32;
ga(34) = 32;
ga(35) = 27;
ga(36) = 27;
ga(37) = 35;
ga(38) = 33;
ga(39) = 33;
ga(40) = 29.43;
ga(41) = 29.43;
% Voxel coordinates determined with FreeView.
bman{ 1} = [ 62 86 28]; lman{ 1} = [ 70 89 37]; rman{ 1} = [ 62 96 36];
bman{ 2} = [ 73 87 21]; lman{ 2} = [ 87 83 27]; rman{ 2} = [ 87 90 20];
bman{ 3} = [ 63 88 28]; lman{ 3} = [ 75 86 35]; rman{ 3} = [ 67 94 38];
bman{ 4} = [ 74 84 22]; lman{ 4} = [ 84 86 30]; rman{ 4} = [ 89 91 24];
bman{ 5} = [ 66 85 28]; lman{ 5} = [ 78 87 35]; rman{ 5} = [ 69 92 38];
bman{ 6} = [ 75 86 20]; lman{ 6} = [ 86 86 28]; rman{ 6} = [ 87 95 22];
bman{ 7} = [ 65 89 27]; lman{ 7} = [ 71 77 34]; rman{ 7} = [ 64 84 37];
bman{ 8} = [ 63 85 28]; lman{ 8} = [ 73 85 36]; rman{ 8} = [ 67 95 36];
bman{ 9} = [ 79 83 22]; lman{ 9} = [ 86 84 31]; rman{ 9} = [ 91 90 26];
bman{10} = [ 53 56 16]; lman{10} = [ 46 56 7]; rman{10} = [ 46 66 7];
bman{11} = [ 48 39 19]; lman{11} = [ 48 25 11]; rman{11} = [ 37 30 13];
bman{12} = [ 61 65 19]; lman{12} = [ 49 54 23]; rman{12} = [ 52 54 12];
bman{13} = [ 54 84 19]; lman{13} = [ 52 77 31]; rman{13} = [ 45 74 25];
bman{14} = [ 54 83 28]; lman{14} = [ 52 75 16]; rman{14} = [ 43 78 20];
bman{15} = [ 44 84 30]; lman{15} = [ 32 73 37]; rman{15} = [ 39 77 45];
bman{16} = [ 49 85 34]; lman{16} = [ 39 74 28]; rman{16} = [ 43 80 20];
bman{17} = [ 66 72 37]; lman{17} = [ 71 59 35]; rman{17} = [ 79 66 38];
bman{18} = [ 73 87 29]; lman{18} = [ 78 89 19]; rman{18} = [ 86 85 23];
bman{19} = [ 67 72 23]; lman{19} = [ 54 74 26]; rman{19} = [ 61 71 31];
bman{20} = [ 70 96 18]; lman{20} = [ 86 94 25]; rman{20} = [ 76 88 28];
bman{21} = [ 71 97 18]; lman{21} = [ 88 93 23]; rman{21} = [ 78 88 28];
bman{22} = [ 70 96 56]; lman{22} = [ 88 86 54]; rman{22} = [ 79 81 49];
bman{23} = [ 49 53 17]; lman{23} = [ 44 54 33]; rman{23} = [ 52 46 31];
bman{24} = [ 44 55 15]; lman{24} = [ 43 57 31]; rman{24} = [ 38 49 29];
bman{25} = [ 57 49 15]; lman{25} = [ 63 46 28]; rman{25} = [ 69 42 22];
bman{26} = [ 49 83 56]; lman{26} = [ 53 74 46]; rman{26} = [ 51 82 42];
bman{27} = [ 48 83 51]; lman{27} = [ 42 68 51]; rman{27} = [ 33 76 48];
bman{28} = [ 39 42 32]; lman{28} = [ 52 45 43]; rman{28} = [ 57 42 32];
bman{29} = [ 44 68 22]; lman{29} = [ 44 56 28]; rman{29} = [ 53 59 21];
bman{30} = [ 42 51 10]; lman{30} = [ 56 52 16]; rman{30} = [ 51 61 20];
bman{31} = [ 31 38 21]; lman{31} = [ 20 48 27]; rman{31} = [ 28 47 34];
bman{32} = [ 29 26 19]; lman{32} = [ 26 36 32]; rman{32} = [ 36 34 31];
bman{33} = [ 51 71 77]; lman{33} = [ 57 62 72]; rman{33} = [ 48 69 65];
bman{34} = [ 53 74 82]; lman{34} = [ 48 63 66]; rman{34} = [ 46 74 64];
bman{35} = [ 53 84 19]; lman{35} = [ 52 76 31]; rman{35} = [ 45 74 25];
bman{36} = [ 53 83 28]; lman{36} = [ 52 75 16]; rman{36} = [ 42 78 20];
bman{37} = [ 58 94 40]; lman{37} = [ 53 89 24]; rman{37} = [ 68 91 27];
bman{38} = [ 69 95 17]; lman{38} = [ 60 82 21]; rman{38} = [ 71 79 19];
bman{39} = [ 95 100 24]; lman{39} = [ 82 97 17]; rman{39} = [ 78 99 25];
bman{40} = [ 74 98 26]; lman{40} = [ 73 84 25]; rman{40} = [ 72 88 19];
bman{41} = [ 72 99 28]; lman{41} = [ 75 86 30]; rman{41} = [ 67 87 26];
% FreeSurfer/FreeView starts counting voxels at 0, we do at 1.
bman = cellfun(@(x)x+1, bman, 'uniformoutput', 0);
lman = cellfun(@(x)x+1, lman, 'uniformoutput', 0);
rman = cellfun(@(x)x+1, rman, 'uniformoutput', 0);