-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathdraw_surface_bluewhitered_gallery_dull.m
100 lines (94 loc) · 3.61 KB
/
draw_surface_bluewhitered_gallery_dull.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
function fig = draw_surface_bluewhitered_gallery_dull(surface_to_plot, data_to_plot, hemisphere, medial_wall, with_medial)
% draw_surface_bluewhitered_gallery_dull.m
%
% Draw multiple data on surface using blue-white-red colormap for
% negative-zero-positive values with dull plot lighting
%
% Inputs: surface_to_plot : surface structure with fields
% vertices - vertex locations [Vx3], V = number of vertices
% faces - which vertices are connected [Fx3], F = number of faces
% data_to_plot : data to plot [VxP]
% P = number of independent data
% hemisphere : which hemisphere (string)
% lh - left hemisphere
% rh - right hemisphere
% medial_wall : indices of the medial wall (vector)
% with_medial : draw medial wall view (boolean)
%
% Output: fig : figure handle
%
% Original: James Pang, Monash University, 2022
%%
if nargin<5
with_medial = 0;
end
if nargin<4
medial_wall = [];
end
num_modes = size(data_to_plot,2);
if with_medial
factor_x_small = 1.02;
factor_x_big = 2.1;
init_x = 0.005;
init_y = 0.01;
length_x = (1-2*init_x)/(factor_x_small + factor_x_big*(num_modes-1) + 1);
length_y = (1-2*init_y);
fig = figure('Position', [200 200 num_modes*400 150]);
for mode=1:num_modes
data_to_plot(medial_wall,mode) = min(data_to_plot(:,mode))*1.1;
clims = [min(data_to_plot(:,mode)), max(data_to_plot(:,mode))];
if clims(2)<=0
clims(2) = 0.01;
end
ax1 = axes('Position', [init_x+factor_x_big*length_x*(mode-1) init_y length_x length_y]);
obj1 = patch(ax1, 'Vertices', surface_to_plot.vertices, 'Faces', surface_to_plot.faces, 'FaceVertexCData', data_to_plot(:,mode), ...
'EdgeColor', 'none', 'FaceColor', 'interp');
if strcmpi(hemisphere, 'lh')
view([-90 0]);
elseif strcmpi(hemisphere, 'rh')
view([90 0]);
end
caxis(clims)
camlight('headlight')
material dull
colormap(ax1,[0.5,0.5,0.5; bluewhitered])
axis off
axis image
ax2 = axes('Position', [init_x+factor_x_small*length_x+factor_x_big*length_x*(mode-1) init_y length_x length_y]);
obj2 = patch(ax2, 'Vertices', surface_to_plot.vertices, 'Faces', surface_to_plot.faces, 'FaceVertexCData', data_to_plot(:,mode), ...
'EdgeColor', 'none', 'FaceColor', 'interp');
if strcmpi(hemisphere, 'lh')
view([90 0]);
elseif strcmpi(hemisphere, 'rh')
view([-90 0]);
end
caxis(clims)
camlight('headlight')
material dull
colormap(ax2,[0.5,0.5,0.5; bluewhitered])
axis off
axis image
end
else
factor_x = 1.05;
init_x = 0.01;
init_y = 0.01;
length_x = (1-2*init_x)/(factor_x*(num_modes-1) + 1);
length_y = (1-2*init_y);
fig = figure('Position', [200 200 num_modes*200 150]);
for mode=1:num_modes
ax = axes('Position', [init_x+factor_x*length_x*(mode-1) init_y length_x length_y]);
obj1 = patch(ax, 'Vertices', surface_to_plot.vertices, 'Faces', surface_to_plot.faces, 'FaceVertexCData', data_to_plot(:,mode), ...
'EdgeColor', 'none', 'FaceColor', 'interp');
if strcmpi(hemisphere, 'lh')
view([-90 0]);
elseif strcmpi(hemisphere, 'rh')
view([90 0]);
end
camlight('headlight')
material dull
colormap(ax,[bluewhitered])
axis off
axis image
end
end