forked from benjaminthuerer/Sourcereconstruction
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathUiO_process_MIDA.m
148 lines (108 loc) · 3.35 KB
/
UiO_process_MIDA.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
function [] = UiO_process_MIDA( MIDA, MIDApath, MIDAwrite )
% Remove air, downsample MIDA (which should contain filename and/or path of MIDA model, and location to write to disk and write the TPMs for segmentation
% This script loads and preprocesses MIDA and calculates tissue propability maps
% (TPM) of the MIDA which will be used in a different script to normalize
% MIDA to subject space.
% load MIDA define TPMs
MIDA = ft_read_mri([MIDApath MIDA]);
cfg = [];
cfg.spmversion = 'spm12';
cfg.method = 'interactive';
cfg.coordsys = 'acpc';
cfg.unit = 'mm';
[MIDA] = ft_volumerealign(cfg,MIDA);
%% Make TPM segments to segment the subject
Tissues = {'Soft','Background','Bone','Gray','White','CSF'};
Soft = [33:35,37,39,51,85,86,55:59,38,42,60,61,63:84,88:96,98,43,62];
Background = [50];
Bone = [52,36,40,41,44:49,53,54,87];
Gray = [2,3:5,7,8,10,16,17,20,21,99,116];
White = [9, 12,18,22,23,100:115];
CSF = [6,32];
for i = 1:length(Tissues)
Segment = ismember(MIDA.anatomy,eval(Tissues{i}));
eval([Tissues{i} ' = Segment;']);
end
Seg = struct;
Seg.Soft = Soft;
Seg.Bone = Bone;
Seg.Gray = Gray;
Seg.White = White;
Seg.CSF = CSF;
Seg.Background = Background;
Seg.dim = size(MIDA.anatomy);
Seg.coordsys = 'acpc';
Seg.transform = MIDA.transform;
Seg.unit = 'mm';
%% write a downsampled version of MIDA without background
MIDA_withoutair = MIDA;
MIDA_withoutair.anatomy(Seg.Background)=0;
cfg=[];
cfg.downsample=2;
cfg.method='nearest';
cfg.smooth='no';
cfg.spmversion='spm12';
MIDA_withoutair=ft_volumedownsample(cfg,MIDA_withoutair);
ft_write_mri([MIDAwrite 'MIDA_withoutair.nii'],MIDA_withoutair,'dataformat','nifti');
%% Write and smooth TPM segments
myTPM=MIDA;
myTPM.anatomy=zeros(myTPM.dim);
myTPM.anatomy(Seg.Gray)=0.999;
myTPM.anatomy(~Seg.Gray)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_gray.nii'],myTPM,'dataformat','nifti');
myTPM=MIDA;
myTPM.anatomy=zeros(myTPM.dim);
myTPM.anatomy(Seg.White)=0.999;
myTPM.anatomy(~Seg.White)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_white.nii'],myTPM,'dataformat','nifti');
myTPM=MIDA;
myTPM.anatomy=zeros(myTPM.dim);
myTPM.anatomy(Seg.Soft)=0.999;
myTPM.anatomy(~Seg.Soft)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_soft.nii'],myTPM,'dataformat','nifti');
myTPM=MIDA;
myTPM.anatomy=zeros(myTPM.dim);
myTPM.anatomy(Seg.Bone)=0.999;
myTPM.anatomy(~Seg.Bone)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_bone.nii'],myTPM,'dataformat','nifti');
myTPM=MIDA;
myTPM.anatomy=zeros(myTPM.dim);
myTPM.anatomy(Seg.CSF)=0.999;
myTPM.anatomy(~Seg.CSF)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_CSF.nii'],myTPM,'dataformat','nifti');
myTPM=MIDA;
myTPM.anatomy=ones(myTPM.dim)*0.999;
myTPM.anatomy(~Seg.Background)=0.001;
cfg=[];
cfg.downsample=2;
cfg.smooth=500;
cfg.spmversion='spm12';
myTPM=ft_volumedownsample(cfg,myTPM);
ft_write_mri([MIDAwrite 'TPM_background.nii'],myTPM,'dataformat','nifti');
disp('Done. All TPMs and MIDA are written into nifti.');
end