Homec4science
Diffusion ISFC clustering (master) Browse

ISFCClustering/Matlab_Scripts/Functions/MakeGraphVideo_NetworkLevel.m
257e25cb59ccmaster

Last Change
Raw File

MakeGraphVideo_NetworkLevel.m
View Options

%% Creates a matlab video object to show how the nodes evolve in activity
% level with time (one per group)
%
% Inputs:
% - Data has size n_ROI x n_ROI x n_TP x n_subjects
% - TC has size n_clusters x n_timepoints
function [X] = MakeGraphVideo_NetworkLevel(Data,CodeBook,Group_Label,Exp_Sphere,View,TC_CTR,TC_ASD)
% Number of time points
n_TP = size(Data,3);
time = 19:2:351;
% I want to average the data across subjects: I get dimensions of n_ROI
% x n_ROI x n_TP
A_CTR = mean(Data(:,:,:,Group_Label==1),4);
A_ASD = mean(Data(:,:,:,Group_Label==2),4);
% Here I have n_ROI x n_DP
D_CTR = squeeze(sum(abs(A_CTR),1));
D_ASD = squeeze(sum(abs(A_ASD),1));
% This is a normalizing factor to have my nodes displayed with possible
% comparison across cases
Fnorm = max(max([D_CTR,D_ASD]));
switch View
case 2
% CTR group
% First frame
PlotBrainGraph_Kmeans(A_CTR(:,:,1),D_CTR(:,1),CodeBook,0,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','jet',1,0.5,[-1 1]);
PlotBrainGraph_Kmeans(zeros(270,270),sum(abs(tmp_plot3),2),...
CodeBook,0,scale_overall/2,scale_overall,1,2,'hot','autumn',1,1,[]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
% Axis for the text
axes('position',[0.1 0.65 0.3 0.1]);
axis off
text(0.1,0.5,[num2str(time(1)),' [s]'],'FontSize',20,'FontWeight','bold');
% Axis for the display of the global number of excursions (or
% for some key connections)
axes('position',[0.15 0.1 0.75 0.15]);
tmp_CM = colorcube(n_clusters+1);
for j = 1:n_clusters
plot(time(1),TC_CTR(j,1),'color',tmp_CM(j,:));
hold on;
end
ylim([-1 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
xlabel('Movie time [s]');
ylabel('Subj. frac.');
drawnow
F = getframe(gcf);
X_CTR(:,:,:,1) = frame2im(F);
clear F
close all
% For each time point, I want to create my frame
for t = 2:size(D_CTR,2)
disp(['Processing CTR frame ',num2str(t),'...']);
PlotBrainGraph_Kmeans(A_CTR(:,:,t),D_CTR(:,t),CodeBook,0,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
% Axis for the text
axes('position',[0.1 0.65 0.3 0.1]);
axis off
text(0.1,0.5,[num2str(time(t)),' [s]'],'FontSize',20,'FontWeight','bold');
% Axis for the display of the clustering outcomes
axes('position',[0.15 0.1 0.75 0.15]);
for j = 1:n_clusters
plot(time(1:t),TC_CTR(j,1:t),'color',tmp_CM(j,:));
hold on;
end
ylim([-1 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
xlabel('Movie time [s]');
ylabel('Subj. frac.');
drawnow
F = getframe(gcf);
X_CTR(:,:,:,t) = frame2im(F);
clear F
close all
end
% ASD group
% First frame
PlotBrainGraph_Kmeans(A_ASD(:,:,1),D_ASD(:,1),CodeBook,0,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
% Axis for the text
axes('position',[0.1 0.65 0.3 0.1]);
axis off
text(0.1,0.5,[num2str(time(1)),' [s]'],'FontSize',20,'FontWeight','bold');
% Axis for the display of the clustering outcomes
axes('position',[0.15 0.1 0.75 0.15]);
for j = 1:n_clusters
plot(time(1),TC_ASD(j,1),'color',tmp_CM(j,:));
hold on;
end
ylim([-1 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
xlabel('Movie time [s]');
ylabel('Subj. frac.');
drawnow
F = getframe(gcf);
X_ASD(:,:,:,1) = frame2im(F);
clear F
close all
% For each time point, I want to create my frame
for t = 2:size(D_CTR,2)
disp(['Processing CTR frame ',num2str(t),'...']);
PlotBrainGraph_Kmeans(A_ASD(:,:,t),D_ASD(:,t),CodeBook,0,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
% Axis for the text
axes('position',[0.1 0.65 0.3 0.1]);
axis off
text(0.1,0.5,[num2str(time(t)),' [s]'],'FontSize',20,'FontWeight','bold');
% Axis for the display of the clustering outcomes
axes('position',[0.15 0.1 0.75 0.15]);
for j = 1:n_clusters
plot(time(1:t),TC_ASD(j,1:t),'color',tmp_CM(j,:));
hold on;
end
ylim([-1 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
xlabel('Movie time [s]');
ylabel('Subj. frac.');
drawnow
F = getframe(gcf);
X_ASD(:,:,:,t) = frame2im(F);
clear F
close all
end
case 1
% CTR group
% First frame
PlotBrainGraph_Kmeans(A_CTR(:,:,1),D_CTR(:,1),CodeBook,0,Fnorm/2,Fnorm/2,2,1,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
drawnow
F = getframe(gcf);
X_CTR(:,:,:,1) = frame2im(F);
clear F
close all
% For each time point, I want to create my frame
for t = 2:size(D_CTR,2)
disp(['Processing CTR frame ',num2str(t),'...']);
PlotBrainGraph_Kmeans(A_CTR(:,:,t),D_CTR(:,t),CodeBook,0,Fnorm/2,Fnorm/2,2,1,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
drawnow
F = getframe(gcf);
X_CTR(:,:,:,t) = frame2im(F);
clear F
close all
end
% ASD group
% First frame
PlotBrainGraph_Kmeans(A_ASD(:,:,1),D_ASD(:,1),CodeBook,0,Fnorm/2,Fnorm/2,2,1,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
drawnow
F = getframe(gcf);
X_ASD(:,:,:,1) = frame2im(F);
clear F
close all
% For each time point, I want to create my frame
for t = 2:size(D_CTR,2)
disp(['Processing CTR frame ',num2str(t),'...']);
PlotBrainGraph_Kmeans(A_ASD(:,:,t),D_ASD(:,t),CodeBook,0,Fnorm/2,Fnorm/2,2,1,'hot','jet',1,0.5,[-1 1]);
set(gcf,'Position',[100 100 700 900]);
set(gcf,'color','w');
axis off
drawnow
F = getframe(gcf);
X_ASD(:,:,:,t) = frame2im(F);
clear F
close all
end
end
end
c4science · Help