forked from KrickPP/Brain_Network_Hypergraph
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathdata_read_COBRE.m
More file actions
42 lines (32 loc) · 1.07 KB
/
Copy pathdata_read_COBRE.m
File metadata and controls
42 lines (32 loc) · 1.07 KB
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
%% 读取脑区数据
data_path = 'D:\Data\code_Du\COBRE\data_ROI';
num_subjects = 146;
ROI_num = 32;
time_length = 150;
ROI_time_series = zeros(time_length, ROI_num, num_subjects);
% 前74control, 后72patient
for i = 1:num_subjects
if i==112 % 112的长度为67
continue;
end
load(fullfile(data_path, sprintf('ROI_Subject%03d_Session001.mat', i)));
ROI_time_series(:, :, i) = cell2mat(data(1, 4:35));
fprintf('ROI_Subject%03d_Session001.mat has been read.\n', i);
end
%% 读取去噪后的脑区数据
data_path = 'D:\Data\code_Du\COBRE\denoise';
num_subjects = 145;
ROI_num = 32;
time_length = 150;
ROI_time_series = zeros(time_length, ROI_num, num_subjects);
% 前74control, 后72patient
for i = 1:num_subjects
if i==112 % 112的长度为67
continue;
end
load(fullfile(data_path, sprintf('ROI_Subject%03d_Condition000.mat', i)));
ROI_time_series(:, :, i) = cell2mat(data(1, 4:35));
fprintf('ROI_Subject%03d_Condition000.mat has been read.\n', i);
end
%% 保存
save('COBRE_denoise_32', 'ROI_time_series');