forked from ossadtchi/PSIICOS
-
Notifications
You must be signed in to change notification settings - Fork 0
/
LoadHeadModels.m
48 lines (46 loc) · 2.28 KB
/
LoadHeadModels.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
function ConData = LoadHeadModels(Conditions, ProtocolDir, Protocol, bUseHR)
% -------------------------------------------------------------------------
% LoadHeadModels: loads head models for each condition recording in a
% specified protocol from brainstorm database and stores them into
% ConData structure
% -------------------------------------------------------------------------
% FORMAT:
% ConData = LoadHeadModels(Conditions, ProtocolDir, Protocol, bUseHR)
% INPUTS:
% Conditions - {1 x Nconditions} cell array storing names of
% conditions
% ProtocolDir - folder path to the protocol which we want to
% load from
% Protocol - brainstorm-generated structure with info about
% protocol. Can be acquired via running brainstorm and
% executing the command
% Protocol = bst_get('ProtocolStudies', ProtocolName);
% bUseHR - boolean flag; if True high resolution gain matrices
% will be loaded as well. If false loads only low res.
%
% OUTPUTS:
% ConData - {1 x N_conditions_total} cell array with each element
% containing a structure for a single condition
% recording from protocol. Gain matrices are stored in
% ConData{i}.HM_LR and ConData{i}.HM_HR
% __________________________________________________________________________
% Alex Ossadtchi [email protected], Dmitrii Altukhov, [email protected]
Ncond = length(Conditions);
sc = 1;
for c = 1:Ncond
for s = 1:length(Protocol.Study)
if(strcmp(Protocol.Study(s).Name, Conditions{c}))
fprintf('Found study condition %s \n ', Conditions{c});
for hm = 1:length(Protocol.Study(s).HeadModel)
if(strcmp(Protocol.Study(s).HeadModel(hm).Comment,'Overlapping spheres_HR'))
if bUseHR
ConData{sc}.HM_HR = load([ProtocolDir Protocol.Study(s).HeadModel(hm).FileName]);
end
else
ConData{sc}.HM_LR = load([ProtocolDir Protocol.Study(s).HeadModel(hm).FileName]);
end
end;
sc = sc + 1;
end;
end;
end;