Homec4science
Diffusion Research Scripts (Thomas Bolton) (master) Browse

Scripts_Sharing/ISFC/Utils/MakeGraphVideo.m
bdfe7aadccccmaster

Last Change
Raw File

MakeGraphVideo.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_conn x n_conn x n_TP x n_subjects
% - TC has size n_clusters x n_timepoints
function [X_CTR,X_ASD] = MakeGraphVideo(Data,CodeBook,Group_Label,Exp_Sphere,View,TC_CTR,TC_ASD,Threshold_CTR,Threshold_ASD)
n_TP = size(Data,3);
time = 19:2:351;
n_clusters = size(TC_CTR,1);
% I want to average the data across subjects: size n_conn2 x n_conn2 x
% n_time for each adjacency
A_CTR = mean(Data(:,:,:,Group_Label==1),4);
A_ASD = mean(Data(:,:,:,Group_Label==2),4);
D_CTR = squeeze(sum(A_CTR,1));
D_ASD = squeeze(sum(A_ASD,1));
Fnorm = max(max([D_CTR,D_ASD]));
switch View
case 2
% CTR group
% First frame
PlotBrainGraph_Kmeans(A_CTR(:,:,1),D_CTR(:,1),CodeBook,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray',1,0.5);
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]);
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([0 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
plot([time(1) time(end)],repmat(Threshold_CTR,1,2),'k--');
xlabel('Movie time [s]');
ylabel('Cluster expression (CTR)');
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,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray',1,0.5);
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([0 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
plot([time(1) time(end)],repmat(Threshold_CTR,1,2),'k--');
xlabel('Movie time [s]');
ylabel('Cluster expression (CTR)');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray',1,0.5);
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([0 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
plot([time(1) time(end)],repmat(Threshold_ASD,1,2),'k--');
xlabel('Movie time [s]');
ylabel('Cluster expression (ASD)');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray',1,0.5);
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([0 1]);
xlim([time(1) time(end)]);
set(gca,'Box','off');
plot([time(1) time(end)],repmat(Threshold_ASD,1,2),'k--');
xlabel('Movie time [s]');
ylabel('Cluster expression (ASD)');
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,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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
otherwise
% CTR group
% First frame
PlotBrainGraph_Kmeans(A_CTR(:,:,1),D_CTR(:,1),CodeBook,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_CTR,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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,Threshold_ASD,Fnorm/2,Fnorm/2,Exp_Sphere,View,'hot','gray');
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