-
Notifications
You must be signed in to change notification settings - Fork 16
/
perceive_stitch_single_interruption_together.m
144 lines (113 loc) · 6.5 KB
/
perceive_stitch_single_interruption_together.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
function data=perceive_stitch_single_interruption_together(recording1, recording2, keeprecordings)
arguments (Input)
recording1 {mustBeFile(recording1)}
recording2 {mustBeFile(recording2)}
keeprecordings {mustBeMember(keeprecordings, ['yes','no'])} = 'no'
end
% For questions contact Jojo Vanhoecke
%
% This is a function to concatenate percept recordings by filling the gaps with NaN's, meant for
% a technical interruption. It reads in the matlab structures, and will create fieldnames.
% If not keeprecordings (default 'no'), the output structure is the same as the input structure.
% If keeprecordings ('yes'):
% The first fieldname is the concatenated file, the second fieldname is the first recording, the third
% fieldname is the second recording.
%
% %% Example:
% % name of the first recording
% recording1 = 'sub-001_ses-Fu12mMedOff03_task-TASK4_acq-StimOff_mod-BrainSenseBip_run-1_part-1.mat';
% % Make sure the recording filename ends on "part-1.mat". Apart from the "part" it must have the same
% % naming as recording2. It needs to be in your path.
% % name of the first recording
% recording2 = 'sub-001_ses-Fu12mMedOff03_task-TASK5_acq-StimOff_mod-BrainSenseBip_run-1_part-2.mat';
% % Make sure the recording filename ends on "part-2.mat". Apart from the "part" it must have the same
% % naming as recording1.
% data=perceive_stitch_interruption_together(recording1, recording2)
load(recording1,'data')
recording1=data;
load(recording2,'data')
recording2=data;
intermission=[recording1.sampleinfo(2)+1 recording2.sampleinfo(1)-1];
intermission_length = recording2.sampleinfo(1) - recording1.sampleinfo(2) + 1;
%sampleinfo = [recording1.sampleinfo intermission recording2.sampleinfo];
if strcmp(keeprecordings, 'no')
data=struct();
assert(strcmp(recording1.datatype,recording2.datatype))
data.datatype=recording1.datatype;
assert(isequal(recording1.label,recording2.label))
data.label=recording1.label;
data.trial={[recording1.trial{1}, nan(size(recording1.trial{1},1),intermission_length), recording2.trial{1}]};
data.time={recording1.time{1}(1):1/recording1.fsample:recording2.time{1}(end)};
assert(isequal(recording1.fsample,recording2.fsample))
data.fsample=recording1.fsample;
data.sampleinfo = [recording1.sampleinfo(1) recording2.sampleinfo(2)];
data.sampleinfo_intermission = intermission;
data.sampleinfo_intermission_length = intermission_length;
if isfield(recording1,'BrainSenseDateTime')
data.BrainSenseDateTime=[recording1.BrainSenseDateTime(1) recording2.BrainSenseDateTime(2)];
data.BrainSenseDateTime_intermission=[recording1.BrainSenseDateTime(end) recording2.BrainSenseDateTime(1)];
end
data.trialinfo = [recording1.trialinfo ; recording2.trialinfo];
assert(strcmp(strrep(recording1.fname,'_part-1',''),strrep(recording2.fname,'_part-2','')))
assert(strcmp(recording1.fname(end-10:end), '_part-1.mat'), 'The file name of recording 1 does not end on _part-1 in data.fname and/or .mat file')
assert(strcmp(recording2.fname(end-10:end), '_part-2.mat'), 'The file name of recording 2 does not end on _part-2 in data.fname and/or .mat file')
data.fname = {strrep(recording1.fname,'_part-1','')};
assert(str2double(recording1.fnamedate) < str2double(recording2.fnamedate))
data.fnamedate = {recording1.fnamedate};
if isfield(recording1,'ecg_cleaned') && isfield(recording2,'ecg_cleaned')
try
data.ecg_cleaned={[recording1.ecg_cleaned, nan(size(recording1.ecg_cleaned,1),intermission_length), recording2.ecg_cleaned]};
catch
data.ecg_cleaned={[]};
end
end
else
data=struct();
data.hdr(2)=recording1;
data.hdr(3)=recording2;
assert(strcmp(recording1.datatype,recording2.datatype))
data.datatype=recording1.datatype;
assert(isequal(recording1.label,recording2.label))
data.label=recording1.label;
data.trial(1)={[recording1.trial{1}, nan(size(recording1.trial{1},1),intermission_length), recording2.trial{1}]};
data.trial(2)=recording1.trial;
data.trial(3)=recording2.trial;
data.time(1)={recording1.time{1}(1):1/recording1.fsample:recording2.time{1}(end)};
data.time(2)=recording1.time;
data.time(3)=recording2.time;
assert(isequal(recording1.fsample,recording2.fsample))
data.fsample=recording1.fsample;
data.sampleinfo(1,:) = [recording1.sampleinfo(1) recording2.sampleinfo(2)];
data.sampleinfo(2,:) = recording1.sampleinfo;
data.sampleinfo(3,:) = recording2.sampleinfo;
data.sampleinfo_intermission = intermission;
data.sampleinfo_intermission_length = intermission_length;
if isfield(recording1,'BrainSenseDateTime')
data.BrainSenseDateTime(1,:)=[recording1.BrainSenseDateTime(1) recording2.BrainSenseDateTime(2)];
data.BrainSenseDateTime(2,:)=recording1.BrainSenseDateTime;
data.BrainSenseDateTime(3,:)=recording2.BrainSenseDateTime;
data.BrainSenseDateTime_intermission=[recording1.BrainSenseDateTime(end) recording2.BrainSenseDateTime(1)];
end
data.trialinfo = [recording1.trialinfo ; recording2.trialinfo];
assert(strcmp(strrep(recording1.fname,'_part-1',''),strrep(recording2.fname,'_part-2','')))
assert(strcmp(recording1.fname(end-10:end), '_part-1.mat'), 'The file name of recording 1 does not end on _part-1 in data.fname and/or .mat file')
assert(strcmp(recording2.fname(end-10:end), '_part-2.mat'), 'The file name of recording 2 does not end on _part-2 in data.fname and/or .mat file')
data.fname(1) = {strrep(recording1.fname,'_part-1','')};
data.fname(2) = {recording1.fname};
data.fname(3) = {recording2.fname};
assert(str2double(recording1.fnamedate) < str2double(recording2.fnamedate))
data.fnamedate(1) = {recording1.fnamedate};
data.fnamedate(2) = {recording1.fnamedate};
data.fnamedate(3) = {recording2.fnamedate};
if isfield(recording1,'ecg')
data.ecg(2) = {recording1.ecg}; end
if isfield(recording2,'ecg')
data.ecg(3) = {recording2.ecg}; end
if isfield(recording1,'ecg_cleaned') && isfield(recording2,'ecg_cleaned')
data.ecg_cleaned(1)={[recording1.ecg_cleaned, nan(size(recording1.ecg_cleaned,1),intermission_length), recording2.ecg_cleaned]}; end
if isfield(recording1,'ecg_cleaned')
data.ecg_cleaned(2)={recording1.ecg_cleaned}; end
if isfield(recording2,'ecg_cleaned')
data.ecg_cleaned(3)={recording2.ecg_cleaned}; end
end
end